ALBERT

Notes: The changes in porosity and elastic moduli of YSZ-containing nickel-based anode materials for solid oxide fuel cells were studied as a function of the fraction of reduced NiO. Anode samples were reduced in a gas mixture of 4% hydrogen and 96% argon for different periods of time at 800°C and their Young's and shear moduli were determined afterward at room temperature using resonant ultrasound spectroscopy and impulse excitation. It was found that the magnitude of Young's and shear moduli decreased significantly with increasing fraction of reduced NiO and that the magnitude of the elastic moduli of a fully reduced Ni–YSZ anode was ∼45% lower than that of unreduced NiO–YSZ. Because the elastic moduli of NiO are close to those of Ni, the observed decrease in the magnitude of the elastic moduli was found to be caused mainly by the significant increase in the porosity of the sample as a result of NiO reduction. Expressions are presented for the amount of porosity and the magnitude of the elastic moduli as a function of the fraction of reduced NiO.

Notes: Two- and three-dimensional SiCr/SiC composites have been prepared starting from Tyranno SA(tm) fiber preforms. Preform densification has been performed by a modified preceramic polymer impregnation and pyrolysis (PIP) process consisting of filling the preform large interbundle voids with SiC powder before the PIP process. This step was accomplished by low-pressure infiltration of a SiC powder dilute slurry through the preform thickness. Specimens were further processed with polymer impregnation and pyrolysis to determine the effects on structural, thermal, and mechanical properties of the obtained composites. High-temperature pyrolysis treatment, which promoted polymer derived SiC matrix crystallization, markedly increased thermal diffusivity.

Notes: Carbon is commonly added to sulfate-fined silicate-glass batches to enhance the fining process. Reactions between carbon and Na2SO4 modify the SOx emissions from Na2SO4 decomposition. Fourier-transform infrared (FTIR) spectrometry is used to analyze the emission of air pollutants from the isothermal decomposition of Na2SO4 + C undertaken using thermogravimetric analysis (TGA). The FTIR spectrometer is calibrated using standard gas mixtures containing CO, CO2, SO2, NO, and NO2. The collected spectra are quantified using the classical least-squares (CLS) approximation. The TGA-FTIR system provides SOx, and COx, concentrations versus time data from the isothermal decomposition of Na2SO4, in the presence of a carbon black. Mass spectrometry (MS) complements FTIR by being able to detect SO(g).

Notes: Different ratios of the precursor phases of SrFeO3–x (SFO) and SrMoO4 (SMO) were used to prepare Sr2FeMoO6 (SFMO) by a solid-state reaction. X-ray diffraction was used to identify the phases. A residual SMO was observed to exist in the sample with an SFO/SMO ratio of 0.9:1. The sample with a residual SMO phase had higher resistivity, lower magnetization, but higher low-field magnetoresistance (LFMR). High-resolution transmission electron microscopy was used to identify the compositions and phases. Nanometer-sized amorphous-like clusters of SMO phase were located inside the grains rather than at grain boundaries; however, some boundaries were rich in the strontium ion. The possible mechanisms for the conduction and the increase of LFMR of SFMO are discussed.

Notes: Submicro- and nano-sized liquid-phase-sintered SiC ceramics were mechanically tested by nanoindentation in the peak load range 5–400 mN. The submicro-sized sample showed a marked indentation size effect which the nano-sized samples did not exhibit. The relevance of indentation depth with respect to the microstructural scale has been outlined. In the investigated grain-size range, the hardness dependence on the grain size could be described by a load-dependent inverse Hall–Petch relation. Young's modulus was less microstructure- and load-dependent. Because of the very fine microstructure, the nano-sized SiC materials gave lower elastic values than the submicro-sized SiC ceramic.

Notes: Samples of composition Ba1−xLaxTi1−x/4O3, x= 0, 0.003, 0.03, and 0.10, were prepared by an alkoxide sol–gel route with final firing of ceramics at 1100°C, 2 h in air. All samples showed bulk insulating behavior with no evidence of semiconductivity caused by either direct donor doping or oxygen loss.

Notes: We have studied the rheological property evolution and hydration behavior of white and ordinary portland cement (type I) pastes and concentrated cement–polyelectrolyte suspensions. Cement composition had a marked effect on the elastic property evolution (G′(t)) and hydration behavior of these suspensions in the presence of poly(acrylic acid)/poly(ethylene oxide) copolymer (PAA/PEO), even though their affinity to adsorb such species was nearly identical. Both white and ordinary portland cement pastes exhibited G′0 values of ∼104 Pa and fully reversible G′(t) behavior until the onset of the acceleratory period (t= 2 h), where the pastes stiffened irreversibly. In contrast, cement–PAA/PEO suspensions exhibited G′0 values of ∼1 Pa and G′(t) behavior comprised of both reversible and irreversible features. Interestingly, ordinary portland cement–PAA/PEO suspensions experienced a gel-to-fluid transition on high shear mixing at short hydration times (〈1 h), and the particle network did not rebuild until ∼24 h of hydration. In sharp contrast, white portland cement–PAA/PEO suspensions remained weakly gelled throughout the initial stage of hydration even after high shear mixing. At longer hydration times (〉1 h), both cement–PAA/PEO suspensions exhibited G′i(t) ∼ exp(t/τc) with τc values of 5.6 and 1.3 h for ordinary and white portland cement, respectively. Our observations suggest that hydration phenomena impact interparticle forces during early stage hydration and, ultimately, lead to initial setting through the formation of solid bridges at the contact points between particles within the gelled network.

Notes: A new ceramic freeze-casting technique capable of manufacturing near room temperature with a sublimable vehicle was accomplished. Fluid-concentrated slurries of Al2O3 powder in molten camphene (C10H16) were prepared at 55°C. These slurries were quickly solidified (frozen) at room temperature to yield rigid solid green bodies, followed by frozen camphene removal by sublimation (freeze-drying) with negligible shrinkage. Sintering without any special binder burnout process yielded sintered bodies with over 98% theoretical density. The proposed advantages include (1) elimination of extremely cold temperatures, (2) elimination of troublesome binder burnout process, and (3) fast manufacturing cycle due to quick solidification.

Notes: Dense composites in the Ti-B-N system have been produced by reactive hot pressing of titanium and BN powders. The effect of the addition of a small amount of nickel (1–3 wt%) on the reaction kinetics and densification of TiN–TiB2 (40 vol%) composite has been studied. Composites of ∼99% of theoretical density have been produced at 1600°C under 40 MPa for 30 min with 1% nickel. The hardness and fracture toughness of these composites are 24.5 ± 0.97 GPa and 6.53 ± 0.27 MPa·m1/2, respectively. The microstructural studies on samples produced at lower temperatures indicate the formation of a transient liquid phase, which enhances the kinetics of the reaction and densification of the composite.

Notes: This paper presents new findings on ultrasonic spray pyrolysis of zirconium hydroxyl acetate precursor drops whose sizes were precisely measured using laser light diffraction technique. Precursor concentration plays a predominant role in determination of product particle size. At 0.01 wt% precursor concentration, conventional spray pyrolysis at 750°C using precursor drops 5–8 μm in diameter, generated by an ultrasonic nebulizer at 2.66 MHz, yielded uniform spherical yttria-stabilized zirconia (YSZ) particles 73 nm in diameter measured by scanning electron microscopy. The YSZ particle diameters were much smaller than those predicted by the one-particle-per-drop mechanism. Under similar reaction conditions, the high-throughput ultrasound-modulated two-fluid (UMTF) spray pyrolysis of larger precursor drops (28-μm peak diameter) also yielded spherical dense particles; they were significantly smaller in size than those produced by the low-throughput conventional ultrasonic spray pyrolysis of smaller drops (6.8-μm peak diameter).

Notes: Thermoelectric elements consisting of the layered polycrystalline materials of Al-doped ZnO and NaCo2O4 were prepared using the pulse electric-current sintering (PECS) method at 900°C for 3 min. Direct contact between the polycrystalline Al-doped ZnO and the NaCo2O4 was obtained in a single-step process for the stacked powders. The electrical conductivities of the polycrystalline materials prepared by PECS were higher than those of materials prepared by conventional sintering, despite their porous structure. The thermoelectric voltage of the 1-mol%-Al-doped ZnO and NaCo2O4 polycrystalline element (measuring ∼6 mm × 3 mm × 15 mm) was 83 mV at dT= 500 K, when the junction of the elements was at 800°C.

Notes: Pb(Zr,Ti)O3–Pb(Mn1/3 Nb2/3)O3 (PZT–PMnN) system has been studied for high-power piezoelectric applications. This study investigates this system to find out the composition with high-power density piezoelectric characteristics and low tem-perature coefficient of resonance frequency (TCF). It was found that the composition 0.9PZT–0.1PMnN (Zr/Ti = 0.51/0.49) modified with 6 mol% Sr exhibits a TCF of −8 ppm/°C (−20 to +80°C). Further, the dielectric and piezoelectric properties of this composition are as follows: kp= 0.53; Qm= 800; d33= 274; ε33/ε0= 1290 and tan δ=1.1%, which shows the suitability of this composition for ultrasonic devices used under fluctuating thermal environment.

Notes: Beam bending is an excellent method for measuring low permeabilities (≤10−18 m2) in homogeneous materials, because it is fast, requires no high pressure, and provides a concurrent measurement of the modulus of the material. The method was previously analyzed and substantiated for cylindrical or square beams. Recently, the analysis was extended to include isotropic and transversely isotropic rectangular beams. In this paper, the analysis is applied to measurements performed on cement paste, and it is shown that the solution for isotropic rectangular beams accounts for changes in the hydrodynamic behavior caused by changing the aspect ratio of the sample. The permeability and elastic modulus results are verified through comparison to previous measurements on cylindrical beams.

Notes: The deformation behavior of boron- and carbon-doped β-silicon carbide (B,C-SiC) with an average grain size of 260 ± 18 nm containing 1 wt% boron was investigated by compression testing at elevated temperatures. Extensive grain growth during deformation was observed. The stress–strain curves were compensated for grain growth by assuming power-law type of dependence on grain size and strain rate. The stress exponent n was ∼1.3 and the grain size exponent p was ∼2.7 at temperatures ranging from 1593° to 1758°C. The apparent activation energy of deformation Qd was ∼760 kJ/mol, which was lower than the activation energy for lattice diffusion of silicon and carbon in SiC and higher than that for grain-boundary diffusion of carbon in SiC. These results suggest that the deformation mechanism of the fine-grained B,C-SiC is grain-boundary sliding accommodated by the grain-boundary diffusion.

Notes: Chromium-containing stainless steel (SS) is a prospective material for use as an interconnect in solid oxide fuel cells (SOFCs). However, during operations at high temperatures, the growth of oxide scales causes the performance of the interconnect and SOFC as a whole to deteriorate. The coating of SS 446 with a conducting perovskite is a potential method of slowing the growth of oxide scale and, therefore, improving overall SOFC performance. In the present research, the structural characterization of a pure LaCrO3 thin film on the SS 446 substrates has been performed as a model material that can be used as a barrier coating for the metallic interconnect. The deposition of an amorphous La-Cr-O thin film on SS 446 was performed using radio-frequency (rf) magnetron sputtering. The deposited amorphous film was annealed in air to form the desired perovskite phase. The film underwent an amorphous to LaCrO4 phase transition during annealing at 500°C with further transformation to LaCrO3 orthorhombic phase during annealing at 700°C. A self-organized dendritic structure was reported as a result of the perovskite-phase formation. Although formation of various oxides, such as Fe2O3 and Fe3O4, was observed during the annealing of uncoated SS 446 in air, the coating of SS 446 surface with LaCrO3 film prevented formation of various oxide phases at the interconnect surface. The structural characterization of the films and SS 446 surfaces was accomplished using scanning electron microscopy with energy-dispersive X-ray analysis, X-ray diffractometry, micro-Raman spectroscopy, and nanoindentation.

Notes: The density, surface tension, and viscosity of the melts from the PbO-B2O3-SiO2 system have been measured at temperatures in the range 1073–1473 K. The effect of composition on these properties was also investigated. The density of the melt was found to increase linearly with increasing PbO content. Molar volume was derived from the density data, and its deviation from the additivity of partial molar volumes was calculated. These deviations in molar volume from those obtained from additivity rules have been used along with the ratio of various coordination numbers of boron (as reported by Bray) to discuss the structure of the melts. The surface tension was found to decrease with decreasing SiO2/B2O3 ratio, and to increase in the range of the PbO content between 30 and 60 mol%, showing a maximum at ∼60 mol% PbO, and then decreased with further additions. This result suggested that the surface tension would be affected primarily by the B2O3 content in the range of the PbO content between 30–60 mol%, and mainly by the PbO content in the range of the PbO content 〉60 mol%, respectively. The viscosity of the melt was found to decrease linearly with increasing PbO content. The results obtained indicate that the increase in viscosity with B2O3 was half that of SiO2 (on a molar basis), and an empirical equation has been proposed for the viscosity as a function of mole fraction.

Notes: Core/shell structures have been prepared via a mechanofusion system by employing several kinds of spherical polymers as a core material and Al2O3 powder or a mixture of Al2O3 and SiO2 powders as a shell material. The effect of the kind of core polymers on the quality of the resulting hollow alumina microspheres has been discussed on the basis of the thermal decomposition behavior of spherical polymers used as a core material. A large fraction of hollow alumina microspheres reflecting the shape and the particle size distribution of the core polymer could be fabricated after sintering at 1600°3C for 3 h, when highly cross-linked poly(methyl methacrylate) (PMMA) microspheres with a gel fraction of 99.03% were used as a core polymer, and abrupt firing at temperatures higher than 500°3C was adopted to remove the PMMA microspheres. The addition of 5 mass% SiO2 to the Al2O3 shell layer was found to be useful for maintaining the spherical shell structure during the firing process and for fabricating a large fraction of hollow alumina microspheres after the sintering.

Notes: Well-defined and stoichiometric spherical particles of BaTiO3 of narrow size distribution were produced at 82° and 92°3C by precipitation from chloride solutions in a strong alkaline environment. The size of the particles can be tailored in the range from ≅103 to 70–80 nm by increasing the barium concentration from ≅0.07 to 0.7 mol/L. The particles are composed of tight aggregates resulting from the assembly of several nanocrystals. The size of the nanocrystals decreases from 200–300 to 30–40 nm by increasing reactant concentration. At low barium concentration (≤0.07 mol/L at 82°3C, ≤0.06 mol/L at 92°3C), formation of BaTiO3 is strongly slowed down and nonstoichiometric, Ti-rich powders are produced. Under these conditions, the particles have the tendency to develop a dendritic-like morphology.

Notes: Crystals of δ-Y2Si2O7 (space group P121/c1) were examined using high-temperature powder X-ray diffractometry to determine their unit-cell dimensions from 296 to 1473 K. The lattice deformation induced by thermal expansion was investigated using matrix algebra analysis to determine the directions and magnitudes of the principal distortions (Λi, i= 1,2, and 3). The directions of Λ1 and Λ3 were defined by the acute angle Λ1c, which linearly decreased from 5(2)° to —5.5(3)° with increased temperature from 504 to 1473 K. The Λ2-axis invariably coincided with the crystallographic b-axis. The magnitudes of Λ1 and Λ2 steadily increased to, respectively, 1.0061(1) and 1.0068(1) during heating to 1473 K, while Λ3 remained almost constant for the entire temperature range. The mean principal distortion, Λm (= (Λ1+Λ2+Λ3)/3), steadily increased to 1.0044(1) with increased temperature to 1473 K. The coefficient of mean linear thermal expansion (α) was derived from the mean principal strain (Λm - 1) as α= (Λm - 1)/ΔT. The temperature dependence was determined to be α= 2.03 times 103+ 1.36(T - 296) (10-9 K-1). Provided that the rule-of-mixtures holds for the Y2Si2O7/Y2SiO5 composites as protective coating on SiC substrates, the volume fractions of 0.72-0.77 (70–75 mass%) would be necessary for the Y2Si2O7 component to match the α-values of both materials.

Notes: Microcellular silicon oxycarbide open cell ceramic foams were fabricated from a silicone resin. Microcellular foams, with a cell size ranging from ∼1–80 μm, were fabricated using poly(methyl methacrylate) microbeads as sacrificial templates. The compression strength of the foams decreased with increasing cell size.

Notes: Electroconductive zirconia-toughened mullite (TiN/ZTM) intragranular nanocomposite was fabricated by hot-pressing a powder mixture of nano-sized TiN, ZrO2(2Y), and mullite gel. The material showed a good sinterability and could be highly densified at a low temperature of 1300°3C. Sintering temperature strongly influenced the microstructure and electrical resistivity of the material. The electrical resistivity increased monotonously from 20 Ω-cm to 1.5 times 106°3Cm, as the sintering temperature was increased from 1300° to 1500°3C. TEM results indicated that such a phenomenon could be ascribed to the changes in the microstructure of the material, which led to a decrease in the connectivity of the TiN network in the sample as the sintering temperature was increased.

Notes: Effects of fluids on material removal rate, chipping damage, and surface roughness in the simulated clinical-dental machining of a dental-type glass ceramic were investigated. Significant differences in removal rate were obtained among the fluids investigated, but only a 4 wt% boric acid solution gave a higher removal rate than conventionally used water. Chipping damage was substantially lower for the boric acid and an oil-emulsion coolant compared with other fluids tested. Surface roughness was independent of the fluids used. The results indicate that improvement can be achieved in both material removal rate and machining damage by the appropriate selection of coolant chemistry.

Notes: This paper reviews the structures and properties of 10 binary, ternary, and quaternary crystals within the equilibrium phase diagram of the SiO2–Y2O3–Si3N4 system. They are binary compounds SiO2, Y2O3, Si3N4; ternary compounds Si2N2O, Y2Si2O7, and YSi2O5; and quaternary crystals Y2Si3N4O3 (M-melilite), Y4Si2O7N2, (N-YAM), YSiO2N (wallastonite), and Y10(SiO4)6N2 (N-apatite, N-APT). Although the binary compounds are well-known and extensively studied, the ternary and the quaternary crystals are not. Most of the ternary and the quaternary crystals simply have been referenced as secondary phases in the processing of nitrogen ceramics. Their crystal structures are complex and not precisely determined. In the quaternary crystals, there exists O/N disorder in that the exact atomic positions of the anions cannot be uniquely determined. It is envisioned that a variety of cation–anion bonding configurations exist in these complex crystals. The electronic structure and bonding in these crystals are, therefore, of great interest and are indispensable for a fundamental understanding of structural ceramics. We have used ab initio methods to study the structure and bonding properties of these 10 crystals. For crystals with unknown or incomplete structural information, we use an accurate total energy relaxation scheme to obtain the most likely atomic positions. Based on the theoretically modeled structures, the electronic structure and bonding in these crystals are investigated and related to various local cation–anion bonding configurations. These results are presented in the form of atom-resolved partial density of states, Mulliken effective charges, and bond order values. It is shown that Y–O and Y–N bonding are not negligible and should be a part of the discussion of the overall bonding schemes in these crystals. Spectroscopic properties in the form of complex, frequency-dependent dielectric functions, X-ray absorption near-edge structure (XANES), and the electron energy-loss near-edge structure (ELNES) spectra in these crystals also are calculated and compared. These results are discussed in the context of specific bonding configurations between cations (silicon and yttrium) and anions (oxygen and nitrogen) and their implications on intergranular thin films in polycrystalline Si3N4 containing rare-earth elements.

Notes: Simultaneous synthesis and sintering of hexagonal α-Ti1−x-Alx(N) (0 ≤x≤ 0.08) solid solutions, which contain a small amount of nitrogen, have been performed by a self-propagating high-temperature combustion method under a nitrogen pressure of 4 MPa. Dense materials (∼99% of theoretical) prepared directly from a mixture of elemental (Ti and Al) powders reveal homogeneous microstructure composed of fine grains (12–16 μm). α-Ti1−xAlx(N) (x= 0.02; Ti0.98Al0.02N0.26) exhibits a three-point bending strength σb of 390 MPa, a Vickers hardness Hv of 9.24 GPa, and a fracture toughness KIC of 4.89 MPa·m1/2; their mechanical properties are much improved by doping Al into α-Ti(N), in comparison with those (σb= 245 MPa, Hv= 9.02 GPa, and KIC= 3.77 MPa·m1/2) of α-Ti(N) fabricated under the same conditions.

Notes: Electrophoretic deposition has been used to synthesize nickel–alumina, functionally graded materials from NiO and alumina suspensions in ethanol. Suspension stability and the kinetics of deposition were studied to determine optimum conditions. Deposition starts with an alumina suspension into which a stream of NiO suspension is injected at various flow rates to obtain the desired profiles. The latter were controlled by varying the deposition current density and component flow rate. The green bodies obtained were sintered in Ar/H2 atmosphere to reduce the NiO to nickel. Various gradation profiles were obtained illustrating the facility of EPD to synthesize continuously graded materials. NiO was used as the precursor for nickel to alleviate settling and rough columnar deposit problems.

Notes: X-ray diffraction (XRD) patterns from nominally β-SiC specimens often differ from those expected for the cubic crystal structure. These differences include the presence of additional peaks, enhanced background intensities, peak broadening, changes in relative peak heights, and shifts in peak positions. It has long been recognized that they are due to the presence of stacking faults, and models relating the experimental observations to stacking fault population have continued to evolve. The presence and relative magnitude of these features vary among different β-SiC specimens. In this work, computer simulations were used to show that the variations are closely related to differences in the type and spatial distribution of stacking faults in each specimen. In these simulations, stacking sequences were generated using a selectively activated 1-D Ising model with a Boltzmann-type probability function for specifying errors, which allows a wide variety of fault configurations to be generated. Direct correlations between different features in the XRD data to the underlying fault population are demonstrated, which are discussed in this paper. It is also shown that this computer model is general, in the sense that many of the models presented in prior work can be interpreted as limiting cases of it.

Notes: This second part deals with a subject analogous to the one developed in the first part, but related to surface irregularities with a non-cylindrical structure. In this part we have taken into consideration both the effects due to potential electrodes and the effects due to current electrodes.The study of these cases has been made with the help of mathematical theory, but it has been necessary to extend this theory from the case of cylindrical bodies to the case of non-cylindrical ones.Finally, as conclusion to both parts of the paper, some pages have been dedicated to the details of the calculation procedure, in order to demonstrate its practical htness for the study of this problem.

Notes: The usual computation methods in refraction seismic in the Sahara often involve considerable shot distances. Therefore, in addition to a high consumption of explosives, difficulties of operation and interpretation arise.In order to cut down the shot distances, it was necessary to study with great care the refracted waves having the shortest possible traveling path along the marker. Such paths correspond to the zone of great relative energy resulting, on one hand, from the sudden increase of reflected energy beyond the critical angle (total reflection) and, on the other hand, from the sudden appearance of refracted energy. The proposed method is founded on correlations of traces located in this “critical zone” on the different records shot on a profile; in order to get around certain difficulties due to insufficient knowledge of the marker velocity, the interpreter is led to select a constant correlation distance, at least for sections of the profile.The paper describes in detail the setting up and utilization of the proposed method, analyses a few examples and gives an estimate of the advantages and the restrictions of this method.

Notes: It was felt that Faust's first law was the best guide in classifying data concerning velocities and their correlation with lithology, even if, on a local basis, other laws fit the observed values better.Data have been re-examined which were provided chiefly by borehole measurements, both conventional and continuous, but also by some particularly reliable velocity profiles. “Faust's coefficient” has been computed for a series of geological formations, i.e. the ratio of true velocity to the value resulting from Faust's law, according to the age and present depth of the formation.Some graphs are shown, in the first place concerning the Sahara, and especially the Palaeozoic. It can be seen that the only formation showing values notably less than I disappears.On the other hand, values greater than I can be noticed for the quartzitic Ordovician of limited thickness.This comparison has been continued on a broader basis, using, in particular, data published by Faust and Wyrobek. The data concerning the very homogeneous Westphalian result in an average value slightly less than I.The Aquitanian basin has been studied in greater detail, owing to the closer spacing of the measurements.Generally speaking, the Middle and Upper Jurassic, the Cretaceous and the Eocene-Oligocene show values notably higher than the Palaeozoic in the Sahara. This seems to be connected with the higher proportion of limestones and dolomites. Values as high as 1.7 are observed in compact dolomites.Four maps are given, concerning the Oligocene, the Eocene-Paleocene (one of them excluding the Danian) and the Upper Cretaceous (exclusive of Danian), and a cross-section showing the link with lithology.Some discrepancies with regard to data from borehole samples are discussed.The question may be raised as to whether such studies can be carried out at the very beginning of exploration of a basin, on the basis of only a few logged wells. In the Aquitaine basin, data given by velocity profiles are in close agreement with those from wells as regards the Tertiary; but as regards the Mesozoic, information is meagre. In basins of the Sahara type, synthetic records have shown that interference of true and multiple reflections with each other is a major drawback. But on the other hand, the possibilities of a combined use of reflection and refraction in order to define lateral changes in velocities may not have been fully appreciated. Some computed examples are shown, since the results of a survey at present being carried out on this basis cannot yet be published.Finally, mention is made of the part played by the determination of anisotropy, and an appeal is made for more coöperation between geophysicists in this sphere.

Notes: The vertical component of the attraction exerted by a horizontal homogeneous polygon is obtained as the sum of terms proportional to the difference, for each of the vertices, between the external angle and its spherical perspective. A set of computation charts is presented, each for a fixed value of the angle of the polygon (22° 30′, 30°, 45°, 67° 30’ and 90°). These charts give the difference between the external angle and its spherical perspective as a function of the quotient of the distance from the sides to the vertical through the station over the depth. The chart for the right angle is an alignment chart, permitting a more accurate reading; the others are graphs in rectangular coördinates.For a remote station the result would be obtained as a small difference between relatively large quantities, and consequently with a not quite satisfactory precision. For these stations it is preferable to assume the entire mass to be concentrated in the centre of gravity.

Notes: In the laboratories of the agip mineraria a new apparatus has been developed which, utilizing reflected energy collected by several groups of seismometers, provides on a single trace an indication of the total energy level. This indication is controled in time according to an exponential characteristic, thus compensating for the attenuation which is effected in the subsurface. The apparatus enables indications of successive energy levels to be obtained for different shots, thus furnishing a time section even during the routine field work.

Notes: The factors affecting the response of a velocity-sensitive well geophone are discussed and it is shown that the geophone movement in response to a seismic pulse is primarily due to the pressure difference between its ends. Several practical examples are calculated to show that the first arrival waveform recorded by a geophone depends on the distance from an interface and that these changes in recorded waveform may give errors of up to 3 msec in well-velocity surveys. These errors may be minimised by using pressure geophones or by avoiding measurements at geological interfaces.

Notes: Surface charges in seismic refraction work in the sand dunes of the Sahara have been necessitated by the many problems involved in shot hole driving. These include lost circulation, unconsolidated surface sands and gravels, caving, poor shot efficiency etc.Tests were conducted in September 1958 in an attempt to evaluate some of the problems involved in surface shooting. This method is complicated by the large number of variables, 8 major ones, influencing the shot pattern.It is difficult to evaluate each of such a large number of variables by an exhaustive and systematic test procedure due to the prohibitive number of shots that this would involve. However, two series of tests contributed greatly to an understanding of the problem, although they left some points to be cleared up.Firstly shots were fired using small charges (less than 25 kg). These gave interesting results, however, they are not applicable to the heavy charges currently in use (500-1500 kg). These results refer primarily to the use of suspended charges which gave superior results to those obtained from surface charges layed on the ground. They also tend to indicate the influence of the weight of the charge. This is manifested by a saturation in seismic energy, as the size of the charge increases.The results obtained from the shooting (200 kg) of several individual charges are not as pronounced. The following conclusions are drawn:Nitrate explosives, in current use in the Sahara for several years, give good results in refraction.– The use of suspended shots does not result in sufficient gain in energy over the ground shots to justify their more complicated utilization: preparation, partial destruction of supports, safety etc. Hence further shooting of suspended charges has been abandoned.– The detonation of individual charges by use of primacord layed parallel to the line, in preference to simultaneous electrical detonation, gave superior results from the energy point of view (especially for high frequency arrivals) and simplified lay-out procedures. Hence this method was adopted.– As for the small charges the relation between seismic energy and weight of the individual charges shows a saturation effect which is less pronounced.–The distance between the individual charges appears to be the most important variable. The relation, seismic energy-distance between charges, indicates an increase in a linear manner for distances of 5 to 20 meters. It would appear as though the beneficial effect for greater distances diminishes beyond this point.Finally, it would appear as though it would be preferable for a given overall charge to reduce the weight of the individual charges and increase the distance between the individual charges. Obviously this rule has certain practical and economic limits. If the cost of the unit of seismic energy is evaluated it is preferable to utilize large distances between charges (20 m) but also relatively large individual charges (about 25 kg).The testing program was followed by nine months of surveying in different types of terrain. The parameters were etablished from the results of the testing program.Economically the following conclusions may be drawn: The average explosives consumption increased from 1000 kg for the shot–hole shots to 1500 kg for the surface shots. This increase, costwise, represents almost exactly the cost of shot hole drilling at the present price of explosives. Efforts are being made to have the price of explosives reduced and should this realize the net result will be a definite cost advantage in favor of surface shooting.In addition, the introduction of magnetic recording in refraction work to the Sahara in 1959 appears to indicate a substantial reduction in the size of the charges needed, hence, a further advantage to surface shooting.

Notes: The decay rate of gamma radiation emitted by a Jamaican mineral spring was measured using portable rate meters. The results of these experiments, supported by auxiliary tests, suggested that the main radioelement causing the radioactivity was radon 222. Fluctuations in the intensity of the radiation were studied in situ. An inverse relation between the intensity and the discharge rate of water was established and the effect of rainfall and of earthquakes was investigated. Some suggestions are made concerning the possible source of radon.

Notes: The iron ores of the Siegerland-district are found in lodes occurring in schists of Lower Devonian age. In the course of the geophysical exploration for these ores, an attempt was also made to clarify the structure and the tectonics of the substratum of the sedimentary rock formations by seismic reflection measurements from the surface. Reflections were recorded in great numbers from layers of Silurian and Devonian age. The form of the folding of the strata derived therefrom corresponds well with that deduced from geological investigations. Beyond that, correlatable reflections were recorded from horizons down to a depth of nearly 25 000 m. It turned out that the iron ore lodes known from the mine workings occur almost vertically above those places where greater dislocations appear in the deepest substratum, and which geologically could be interpreted as fracture zones. The results so far obtained seem to be indicative of the fact that the tectonics of the deepest substratum and the occurrence of ore deposits are related to each other.

Notes: Precision calibration tests carried out with a Worden gravimeter in environment conditions as varied as possible have led to the following conclusions. The secular variation curve of the calibration shows that the secular relative variation has been very appreciable in the first year of life of the instrument, precisely of 1.3 × 10−3 in 13 months, while it has been practically zero in the second year. This must be ascribed to the fact that at the beginning of the operations the gravity-meter had just arrived from the factory with its measuring body renewed.On the other side the effect of 24°C of temperature difference between two series of determinations reveals itself in a variation of 0.0022 of the value of the calibration, so that it can be established that the variation of the calibration is I‰ for each 10°C of temperature variation.

Notes: Due weight must be given to three factors (1) its purpose, (2) the limitations imposed by nature, industry and environment, (3) the effect of errors in auxiliary measurements. Together the last two factors give an estimate of the accuracy which can be achieved. If this is sufficient for the instrument to fulfil its purpose design and construction can proceed. If it is not, a detailed survey of the last two factors is necessary in order to estimate the research effort that will be required when directed into the most effective channels. If the project is considered worth the effort it can proceed. These factors are discussed in relation to a borehole gravimeter. It is shown that the purpose will be served only if the accuracy is comparable with the best surface gravimeters. Borehole versions of some surface instruments are considered and it is concluded that none is suitable without considerable modifications. A successful borehole gravimeter is likely to be a completely new instrument.

Notes: Recently, quite good determinations of interfaces have been made by means of seismic well surveys. These were carried out for various oil companies in Northwest Germany, especially on salt-domes, and in Southern Germany, and in the area of the Upper Rhine Valley. It was the purpose to determine the configuration and position of salt-flanks and important planes of stratification in the neighbourhood of deep wells, and thus to reduce the risk of expensive wells.For the delineation of the boundary surfaces the method described by Gardner in “Geophysics 1949” was applied.The problems arising and the results obtained are discussed in connection with some interesting examples. From them we see, above all, that a knowledge as precise as possible of the seismic velocities of all media involved in the measurement is very important. A good idea of the general stratification in the area of deep wells seismically surveyed will make possible a rational and adequate planning of the survey program and will facilitate the interpretation.

Notes: A method is described in which the major part of gravimetric terrain corrections may be estimated by digital computing methods.

Notes: The electric and magnetic field of a direct-current element, embedded in a space consisting of two halfspaces differing in conductivity and permeability, is given in a form suitable for integration along a cable. From these formulae is obtained the field when one end of a D.C. cable is grounded in a highly conducting, pencilformed ore body, considered as a line electrode with uniform current leakage. Numerical examples are given for a horizontal line electrode. The formulae should be useful as a first approximation to the A.C. case. The validity of this is shortly discussed.

Notes: An analogy is established between the propagation of an electromagnetic field in a horizontally stratified conductive medium and that of currents in a linear electric network. This analogy will allow us to know in which case measurements executed at the surface of the ground can provide information about the properties of deeper layers. The situation is the same as when one tries, by means of impedance measurements at the input terminals, to know the properties of electric lines or lumped networks.A symmetrical linear network establishes between:Voltages V1 and V2 (respectively at input and output terminals)Currents I1 and I2 (respectively at input and output terminals), the relationship 〈displayedItem type="mathematics" xml:id="mu1" numbered="no"〉〈mediaResource alt="image" href="urn:x-wiley:00168025:GPR77:GPR_77_mu1"/〉 is called the “phase constant” and Z the characteristic impedance The table 〈displayedItem type="mathematics" xml:id="mu2" numbered="no"〉〈mediaResource alt="image" href="urn:x-wiley:00168025:GPR77:GPR_77_mu2"/〉 is called “characteristic matrix”If one then takes the horizontal components of an electromagnetic field, one can compare the voltage V with the horizontal electric field E (in volts m-1) and the current I with the horizontal magnetic field H (in Amp m−1)For plane waves in vacuum, one gets (Schelkunoff) Z=C (velocity of light)A layer whose thickness is D has a phase constant⌈= 2 πD/λ (λ= wavelength)For plane waves in a conductive, non magnetic medium whose conductivity is δ, one gets 〈displayedItem type="mathematics" xml:id="mu3" numbered="no"〉〈mediaResource alt="image" href="urn:x-wiley:00168025:GPR77:GPR_77_mu3"/〉 For the system of waves created by a point-source one expresses the solution, since it is classical, as a sum or integral of different modes. For each of these modes, one can still write a bilinear relationship similar to (1), and calculate Z and ⌈.A generalization of this result is given.Boundary conditions that exist at the surfaces of separation are met by writing the continuity of E and H. This condition exactly corresponds to the fact that there one V and one I at the junctions of various sections of lines, or networks.The characteristic matrix of a cascade of networks-or here that of a horizontally stratified ground– is the product of the matrixes of each network. If, in addition, one knows the end impedance–here the conductivity λn of the last layer, supposed to be infinite–one can compute the input impedance Re of the whole system. In electromagnetic prospecting, it is Re that conditions the observable field.Inversely, if measurements were infinitely accurate, the. knowledge of the surface field would give all the thicknesses D and conductivities δ of intermediate layers, in the same way that impedance measurements at the input terminals of a cable would allow to locate any failure. The imperfection of measurements causes the elements located too far away (from the point of view of wave attenuation) to escape detection.

Notes: The attention of everyone has already been drawn to the part which accident reports must play (Geophysical Prospecting, March, 1957, PP- 1 to 8). We intend to show here, with supporting examples, that the systematic analysis of the causes, does allow one to obtain practical conclusions.To begin with, it is advisable to widen the notion of geophysical accident as much as possible. This being done, about 200 reports distributed over a long period and a sufficient number of parties, are dealt with. The collected reports are sufficiently numerous to draw conclusions in a general way bearing on the headings: drilling, transport, outbreaks of fire, explosives, falls, stings, shocks, miscellaneous causes.All the quoted figures are referred to the number of accidents reported. Action must be taken to encourage european geophysicists to study accident reports and to analyse the causes. This may be done in such a way that secret information, such as whereabouts of their personnel, is not disclosed.Efforts to promote “Security” must in the first place bear on the quest for the causes, the definitions of ‘accident’ and of ‘geophysicist’ being taken at their widest possible meaning.

Notes: Seismic model experiments are described in which long strips of plexiglass were used as models. One end of the strip was sawn off at an oblique angle and, at the opposite end, the strip was excited by means of a barium titanate transducer. The experiments showed that, if the width of the strip was sufficiently small, an anomalous reflection against the oblique end occurred which travelled in the longitudinal direction of the strip. This anomalous reflection did not occur when the width of the strip was large. These results are explained on the basis of Fresnel's theory. It is inferred that, in the subsurface, refracted waves may be reflected against fault planes without the law of reflection being satisfied, provided that the refracted wave is propagated in a sufficiently thin high velocity layer.

Notes: Changes of shape of seismic waves provide information on the properties of the material in which the waves propagate. Ricker (1953) has attempted to explain the changes of shape on the basis of a simple viscoelastic theory. His conclusions are at variance with those of others who find a dependence of the attenuation on frequency which could be explained only by a much more complicated linear theory or by nonlinear theories.To provide a basis for discussion, the essentials of the theory of viscoelasticity are briefly reviewed. If a relaxation spectrum, rather than one or very few relaxation times, is admitted, a great variety of experimental results can be described by the linear theory of viscoelasticity. A linear theory is indicated when no obvious violations of the principle of superposition occur.Ricker's theory is presented with some modifications which allow for a finite duration of the initial pulse and for the approximate character of his basic assumptions. There do not appear to be serious discrepancies between his theory and his experimental results. Some of the objections to his theory can be met by assuming a finite duration of the initial pulse. However, more direct measurements made under similar circumstances by McDonal et al. (1958) at the same location lead to a conclusion on the nature of the material not in accordance with Ricker's. This casts doubt on the sensitivity of his method.Laboratory measurements usually yield results which are not explainable in terms of simple viscoelastic models. Whether a linear theory with a relaxation spectrum or a nonlinear theory should apply depends much on the experimental conditions. We must also consider the possibility of nonlinear mechanisms which are active at small amplitudes. No stand is taken in this controversy, but it is pointed out that the question linear or nonlinear could be decided experimentally without considering the details of the theories.

Notes: A method of refraction computing is described which establishes the two points of incidence at the surface of the ray-paths from any point on the refractor when observed from opposing directions. The distance between these points is used to determine the radius and the centre of a circle to which the refractor must be tangential. The matching of forward and reverse data is intrinsic in the method so that the resulting profile is automatically positioned in depth. It is shown that if the effects of dip are neglected, errors are introduced in the derived radius and also in the position of the centre but these errors are shown to cancel each other almost exactly.

Notes: The wells Parentis 10 and 12 which are located on the eastern edge of the Parentis anticline were found to be considerably lower at the top of the reservoir series than had been expected.Unfortunately the village of Parentis is located in this important zone and no more seismic results could be obtained which might help locate faults and give a better picture of the area. To help in our understanding of the area, we attempted to investigate the subsurface conditions surrounding the wells by shooting shaped charges or primacord in the well and recording the returning seismic energy along profiles located on the surface. Thus far, we have used this method in three wells: Parentis 10, Parentis 12 and Lugos 2.With the data given by a velocity survey in each well, time-distance curves have been calculated for several typical possibilities: horizontal layers, dipping layers and the effect of a fault. A comparison of the calculated and experimental curves shows anomalies which correspond to a change of the average velocity along the travel path, which can correspond to geological changes such as: a sudden change of dip, a fault or a thickening of the section.The anomalies have been studied on the different profiles from each well and a possible interpretation of the subsurface conditions has been given. The directions of the major accidents determined by this method are in general agreement with those found by the regular seismic shooting.

Notes: The directional sensitivities of an array of geophones depend as much on the size of the array as on the directional sensitivities of the single geophones forming the array. In a former paper (Celminš 1956) the effect of the size of the array was investigated. In addition to this the influence of the directional sensitivities of the individual geophones on the sensitivities of an array of geophones is investigated here.In δ 2 some general formulae for the directional sensitivities are derived. In addition the directional sensitivities of a single geophone according to P- and S-waves are investigated (δ3 and δ4). In δ5 the'directional sensitivities for the most important cases of arrays of parallel geophones are computed. The results of the computation (shown in fig. 4 and 5) give an idea of the types of directional sensitivities which occur in the use of geophone arrays.

Notes: Non-explosive seismic sources, notably the falling weight and the eccentric machine, have occasionally proved advantageous in experimental and practical seismology in spite of their relative feebleness.Another alternative is the pulsed vibrator, which offers the advantage of a completely controllable waveform. Recent theory shows that such a source is quite well adapted to the radiation of high frequency waves from the surface of the ground; it is estimated, for example, that in typical circumstances reflections should be obtainable from depths of the order of 200 metres. The theory also shows how resonances may be avoided so as to radiate the impulsive waveform unimpaired.

Notes: Book review in this articleGedenkboek F. A. Veiling Meinesz, Verhandelingen van het Koninklijk Nederlandsch Geologisch-Mijnbouwkundig Genootschap, Geologische Serie, Vol. XVIII, May, 1957.

Notes: A brief description of the procedure usually adopted for ascertaining the accuracy o aeromagnetic data is given as a background to the discussion of factors which affect this accuracy. These factors fall mainly into two groups: spatial positioning, i.e. relation to the ground surface, and effects contributing to the observed relative magnetic values. The evaluation of the observed anomalies in terms of the most probable causes and the anomalies themselves are directly influenced by the amount of the available control of factors mentioned above. The effects of various factors, such as plan positioning, height keeping, diurnal drift, are discussed, the arguments being based on theoretical and practical premises.

Notes: The effect of the composition of the tank walls on apparent resistivity determinations made in model experiments has been investigated for the two extreme cases of an insulating and a perfectly conducting tank wall: The resulting errors have been determined both by calculation and by experiment.Experiments are described which demonstrate that the magnitude of the effect can be considerably reduced by specially constructed tank walls.

Notes: Amplitude measurements have been made of the height of the first peak of an arrival refracted from a shallow refractor. If the amplitude is assumed to decay as the inverse mth power of the distance, the least squares value for m is found to be 2.16 ± .04. Because of this value and because of the character of the recorded event it is concluded that the arrival is a simple critical refraction. After applying the theoretical ‘spread’ factor for critical refraction there remains a residual attenuation of 1.96 ± 0.28 decibels per 1000 feet. The predominant frequency in the pulse is about 20 c.p.s. and this attenuation agrees with the losses found for such a frequency by extrapolation of the published results of other workers. Although no evidence could be seen on the records for a change of pulse frequency with distance, the quoted result would be consistent with a dependence of residual attenuation on the first power of the frequency, and would be inconsistent with a dependence on the second power of the frequency.It is concluded that studies of the amplitudes of refracted events will give useful estimates of the attenuation factors of rocks.

Notes: A method of constructing depth contour maps of arbitrarily curved horizons obtained from seismic reflection surveys is discussed. This method takes into account three dimensional refraction, avoiding the construction of seismic cross sections of any kind. It requires little work even if refraction at several horizons is taken into account. The multiple layer problem is traced back to the single layer case. Discontinuities in velocity are also taken into account.

Notes: 1. A Worden gravity meter has been tested for sinusoidal variations in the external temperature, by determining simultaneously the corresponding variations of temperature inside the Dewar vacuum bottle, and the variations in the readings.A second Worden gravimeter has been tested for rapid variations in external pressure.

Notes: The results of field and laboratory methods of density determination on a series of Coal Measure, Permian and Triassic rooks are presented and the different methods compared. It is concluded that the most satisfactory method is that of measuring the vertical change of gravity in a mine shaft. Nettleton's method is unsatisfactory to us, due to weathering of the rocks (particularly Magnesian Limestone) and possible effects from drift. Laboratory measurements are of variable value depending on the lithology and source of the samples.A method adopted to solve the problem of finding the true densities for use in a local gravity survey in N.E. England is given.

Notes: The electromagnetic array profiling (EMAP) exploration method can be combined with a direct one-dimensional inversion process for conversion to depth to produce a subsurface resistivity cross-section. This cross-section may then be interpreted in parallel with a seismic cross-section to enhance the prediction of rock type and structure. In complex thrust environments and areas of shallow carbonate rocks, the EMAP method is often used to provide additional data either to help the seismic processor and/or to aid the seismic interpretation. In particular, the electromagnetic (EM) data can be used to build an independent seismic velocity file for depth migration.Three EMAP test areas in the western United States are used to demonstrate such a use of EMAP as an expioration tool. The first shows how a velocity file is estimated from resistivity data for seismic depth migration processing in a complex thrust environment. In the second example, the method is applied in layer-cake geology with high seismic velocity rocks at the earth's surface. The third example is another complex thrust environment, but in this case the velocity file derived from the resistivity data is used for stacking the seismic data.

Notes: The iron ores of the Siegerland district are found in lodes occurring in schists of Lower Devonian age. These lodes, normally dipping almost perpendicularly, could be traced only by seismic reflection measurements underground, up to distances of approximately 500 m. The geophones were attached to the walls of the galleries. The frequencies of the recorded longitudinal and transverse waves ranged from 200 to 500 c.p.s. By this method a new siderite lode was discovered.

Notes: Book Reviewed in this article:Lehrbuch der angewandten Geophysik, Teil II, H. Haalck.

Notes: A series of repeat observations during a period of 5 years reveal a variation with time of the Worden gravimeter No. 142 small dial spring system. A possible change in the gravity difference concerned is ruled out. It is shown that the small dial scale factor has diminished by about 0.25 % within the last two years against about 0.5 % four to five years ago. The ratio of L.D. to S.D. has been analysed and the results obtained show confirmation of the change in the small dial system and prove a consistency in the large dial system.

Notes: A study is made of the possibility of calculating electric potentials, and therefore of interpreting apparent resistivity measurements, when the ground presents structural conditions more complicated than plane, parallel stratification.For this work a schematic type of structure has been chosen, in which the discontinuity planes for electric resistivity are both parallel and perpendicular to the level of the ground.The problem has been examinated in its general aspect, and then some criteria for the numerical calculation have been indicated.Some calculated examples are presented, relevant to vertical electric soundings accomplished by conventional technique, for particular types of structure.

Notes: The elastic impulse-from an explosion is propagated through the earth as an expanding shell bounded by the wavefront. This wavefront and its motion are a physical reality while rays or trajectories are much less sharply defined and cannot actually be observed because the wavelengths involved in seismic processes are too long. This means that the visualization and the interpretation of seismic processes can be carried out in a more satisfactory, simple and comprehensible way with wavefronts than with rays.H. R. Thornburgh (1930) introduced the very lucid and simple method of interpreting seismic refraction sections, using Huygen's principle to construct wavefronts from two symmetrical shotpoints by working back from their observed arrival times at the surface. The depths are found first with only a knowledge of the velocities down to the refractor and after that the velocities in the refractor can be determined.The pattern of the two wavefronts at equal time intervals, obtained when carrying out Thornburgh's construction, leads to the very simple, approximative “Plus-Minus” method. The additions of the travel times from two symmetrical shotpoints to-each geophone give the relative depths and their subtractions give the velocities of the refractor. The method extrapolates exact quantities from the refractor boundary up to the surface by the “Plus”- and “Minus” lines. This extrapolation introduces errors that can, more or less, be corrected for. A number of hypothetical cases are presented to give an insight into the applicability and reliability of the “Plus-Minus” method. Its main applications are for fairly shallow refraction investigations and for determining weathering corrections for refraction and reflection work.

Notes: Major savings for operators of seismic parties have been achieved by the use of prilled ammonium nitrate as a substitute for gelatin dynamite in areas where shot holes are dry. Pound for pound, the prills appear to have the same energy yield as 60% gelatin dynamite, and some improvement of records has been noted where prills are used. Safety problems are less with ammonium nitrate than with dynamite, and the development of efficient field handling techniques prevents the loss of production.

Notes: Some mathematical inequalities are proved and their relevance to the interpretation of local magnetic and gravity anomalies is discussed. They can be applied to the readings of either a vertical magnetometer or a total field magnetometer or a gravimeter and yield a maximum possible value for the depth below the earth's surface of the top surface of the body producing the anomaly.

Notes: A fast inversion technique for the interpretation of data from resistivity tomography surveys has been developed for operation on a microcomputer. This technique is based on the smoothness-constrained least-squares method and it produces a two-dimensional subsurface model from the apparent resistivity pseudosection. In the first iteration, a homogeneous earth model is used as the starting model for which the apparent resistivity partial derivative values can be calculated analytically. For subsequent iterations, a quasi-Newton method is used to estimate the partial derivatives which reduces the computer time and memory space required by about eight and twelve times, respectively, compared to the conventional least-squares method. Tests with a variety of computer models and data from field surveys show that this technique is insensitive to random noise and converges rapidly. This technique takes about one minute to invert a single data set on an 80486DX microcomputer.

Notes: High-speed photography has been used visually to study the shape, surface, turbulence and behaviour of an underwater oscillating bubble generated by an airgun. The source was a BOLT airgun with a chamber volume of 1.6cu.in., placed in a 0.85m3 tank at 0.5m depth. Near-field signatures were also recorded in order to compare the instant photographs of the oscillating bubble with the pressure field recorded about 25 cm from the gun. Estimations of the bubble-wall velocity and bubble radius estimated from high-speed film sequences are also presented, and are compared with modelled results. The deviation between the modelled and measured bubble radii was at most 9%. In order to check the capacity for transmission of light through the bubble, a concentrated laser beam was used as illumination. We found that the air bubble is a strong scattering medium of laser light, hence the bubble is opaque.

Notes: Full-wavefield inversion for distributions of acoustic velocity, density and Q on a vertical slice through a25D model is implemented for common-source gathers in a cross-hole geometry. The wavefield extrapolation used is 3D, so all geometrical spreading, scattering, reflection, and transmission effects are correctly and automatically compensated for. In order to keep the number of unknowns tractable, application was limited to 2.5D models of known geometry; the latter assurnes a prior step, such as tomography, to fix the layer geometries. With the model geometry fixed, reliable solutions are obtained using synthetic data from only two independent source locations. Solutions from data with noisy and missing traces are comparable to those from noise-free data, but with higher residuals. When the source locations are spatially widely separated, conunon-source gathers may be summed and treated as a single wavefield to yield the same model estimates as when the individual source wavefields are treated separately, at substantially reduced cost. Inversions for full 3D parameter distributions can be handled with the same software, requiring only solution for more unknowns.

Notes: It is seen that the apparent chargeability (Ma)L anomaly over a 2D graphite body splits into two distinct (Vs)L anomalies which closely follow the apparent resistivity profile. This suggests that the electric field amplitude is distorted due to a superficial inhomogeneity creating a (Vs)L anomaly, which bears no relation to the polarized body. The target depth obtained by continuation of such a profile is therefore, not acceptable.

Notes: The improvement in accuracy and efficiency of wave-equation migration techniques is an ongoing topic of research. The main problem is the correct imaging of steeply dipping reflectors in media with strong lateral velocity variations. We propose an improved migration method which is based on cascading phase-shift and finite-difference operators for downward continuation. Due to these cascaded operators we call this method‘Fourier finite-difference migration’(FFD migration).In our approach we try to generalize and improve the split-step Fourier migration method for strong lateral velocity variations using an additional finite-difference correction term. Like most of the current migration methods in use today, our method is based on the one-way wave equation. It is solved by first applying the square-root operator but using a constant velocity at each depth step which has to be the minimum velocity. In a second step, the approximate difference between the correct square-root operator and this constant-velocity squareroot operator (the error made in the first step) is implemented as an implicit FD migration scheme, part of which is the split-step Fourier correction term.Some practical aspects of the new FFD method are discussed. Its performance is compared with that of split-step and standard FD migration schemes. First applications to synthetic and real data sets are presented. They show that the superiority of FFD migration becomes evident by migrating steeply dipping reflectors with complex overburden having strong lateral velocity variations. If velocity is laterally constant, FFD migration has the accuracy of the phase-shift method. The maximum migration angle is velocity adaptive, in contrast to conventional FD migration schemes. It varies laterally depending on the local level of velocity variation. FFD migration is more efficient than higher-order implicit FD schemes. These schemes use two cascaded downward-continuation steps in order to attain comparable migration performance.

Notes: In order to eliminate the effect of smoothing due to the use of a finite number of grid points, the second derivative is computed by integrating the product of g with a convenient continuous function, which yields the second derivative to the extent to which the first terms of the Taylor expansion of g represent its value correctly. By applying this method to the anomaly caused by an isolated mass, and to that caused by a homogeneous half plane, it is shown that, if the result obtained is interpreted as if it really were a second derivative, erroneous values for the depth and the mass are obtained. If the real depth of the mass is small, a too large apparent depth is obtained. In the case of a half plane the use of a system of grid points gives the same result. These considerations permit the rational choice of the method of computing the second derivative, such that the effects of too shallow mass irregularities are attenuated.

Notes: From gamma ray absorption the density of soil to a depth of i m. can be determined in situ, without taking samples, by a method which has been described in a previous paper (Wendt 1954) and which has been checked with satisfactory results (mean error about 1%). A nearly vertical steel rod with a gamma ray source at its lower end is pushed into the ground, and the gamma radiation is measured at the surface with two counters placed symmetrically with respect to the rod.As a first step in interpreting the observed effects, the assumption had been made that counters are infinitesimally small. This paper shows how a theory can be developed without this assumption, the finite dimensions of the counters being taken into account by an integral expression. For this purpose a determination of the dependence of the counting rate on the angle of incidence of the gamma rays is required. A comparison of the experimental results with the theoretical ones demonstrates that the integral expression is correct.The influence of deviations from the normal geometry of the arrangement (e.g. caused by bending of the probe) on the results is investigated, as well as the influence of variations of the density with depth.The treatment of the basic problems has been carried through in such a manner that the results obtained can be applied to other absorption methods as well.

Notes: In this paper an attempt has been made to evaluate the possibilities of determining subsurface features from temperature measurements made at shallow depths. The temperatures at the observation depth used are subject to seasonable variations. The amount of this variation differs, and these differences have been attributed to structural features. Local differences in temperature are also caused by surface factors such as the thermal conductivity of the soil, the vegetation and the micro-climate. The differences in temperature due to surface factors are so great that the variation in heat flow caused by differing thermal conductivities of the rocks at depth cannot be detected at the surface. It is not possible therefore to use this method to determine the position of subsurface structures such as domes, anticlines or horsts.It is also shown that concentrations of radioactive elements in the rocks do not provide any measurable heating.The method has been successfully applied to the problem of finding fissures, cracks and similar features provided that convective heat transport from depth to the surface has taken place along these features. In order to compute the area through which a given amount of water at a given depth must pass to give a specific temperature increase a model has been considered. From the shape of the temperature anomaly above a fissure carrying such heated water, the places at which the rising water enters the groundwater stream can be seen immediately. It is possible to estimate the amount of rising water by computing the amount of heat energy transferred to the surface.The method is particularly suitable for determining the position of steam deposits in regions of recent vulcanicity. In such a region a location, hitherto unknown, was found where steam rising from depth condenses beneath the surface.The area and form of the anomaly indicate the extent in depth of this steam bearing zone. By means of observations made at the individual points, it was also possible to compute approximately the energy being released at the surface. This value can give an indication of the extent to which the steam deposits regenerate themselves.Finally it has been shown that percolating canal water can be detected thermally near the canal banks. Such anomalies are dependent on the season.

Notes: The following assumptions are made in the mathematical treatment of the problem. Below a plane earth's surface there is a three-layered elastic medium the interfaces of which are parallel to the earth's surface. The uppermost layer represents the weathered layer in which the velocity of propagation of seismic waves increases linearly with depth. The two lower layers, the so-called intermediate layer and the substratum each have a constant velocity. The surface of the earth is acted on simultaneously by a normal pressure N in the form of a Heaviside pulse. The seismic wave thus generated is propagated through the elastic media.The aim of the investigation is to study the shape of the wave1) in the intermediate layer, after the wave has entered it the first time2) at the earth's surface, after the wave has been reflected once at the interface between the intermediate layer and the substratum.The mathematical solutions can in both cases be expressed as series of Bessel functions. Some numerical examples illustrate the quasi-periodic nature of the solutions. The pseudo-frequency is determined by the gradient of velocity in the uppermost layer; it assumes a value of approximately 50 c.p.s. for a gradient of appr.〈inlineGraphic alt="inline image" href="urn:x-wiley:00168025:GPR328:GPR_328_mu1" location="equation/GPR_328_mu1.gif"/〉

Notes: The President of the E.A.E.G. Safety Committee has stated that European geophysicists fail to attach sufficient importance to accident statistics. Nevertheless these statistics are the best means at our disposal for the study of the factors influencing the number and severity of accidents.It appears essential that two points should be stressed, firstly that the accident reports must be completed correctly, and secondly that the interpretation of each report and of the statistical data must be done systematically.If these conditions are fulfilled it should be possible to create, at all levels, a technical and psychological atmosphere which would do much to protect geophysicists against accidents.

Notes: Simulated annealing is a stochastic combinatorial optimization technique, based on ideas from statistical mechanics, thermodynamics and multivariable probability theory. This paper presents the use of simulated annealing as a means of inversion for both linear magnetics and non-linear resistivity problems. The subsurface is viewed as being constructed of smaller elemental blocks which possess either uniform internal magnetization or conductivity, enabling larger structures to be modelled. Simulated annealing is employed to calculate the distribution of the particular physical property which causes a measured anomalous field curve.A general description of simulated annealing and its application is given, followed by specific descriptions of its application to the magnetics and resistivity cases.For the magnetics case the subsurface consists of 2D prismatic elements as the basis for the forward model. Synthetic model data is used to test the algorithm and an example of actual field data; a survey across an igneous dike is used to demonstrate the use of the method with real data. In the resistivity case, the finite-element method is used to generate the forward models. Synthetic vertical profiling data is used to test the application of the simulated annealing method to the resistivity case. Actual data from an archaeological survey is used to show again the use of the method with real data.Simulated annealing is shown to be capable of inverting both the linear and non-linear methods of magnetic surveying and resistivity surveying respectively.

Notes: We consider the problem of computing the most probable location of a target based on radar measurements of the subsurface. Our algorithm makes use of the maximum likelihood estimator (MLE), which represents a correlation between the measured data and synthetic data generated for the object of interest at different locations. Previous studies assume a plane-wave acquisition geometry and target object(s) embedded in a uniform background. In this paper, a generalization of the MLE method is presented which is valid for discrete point sources (and receivers) and a 2D model (i.e. a 2.5D acquisition geometry). Within this formulation the treatment of a non-uniform background model is also possible. We concentrate on geotechnical ground investigations and assume that the characteristic dimensions of the target object are in the range 1–2λ, (λ being the wavelength). The potential of the method is demonstrated employing cross-hole radar data acquired in a controlled field experiment. The MLE result is also compared with the image obtained employing a full reconstruction method such as diffraction tomography.

Notes: We present a simple method for estimating an effective source wavelet from the first arrival in marine vertical seismic profiling (VSP) data. The method, which utilizes the free-space Green's function of the Helmholz equation, is simple and very computer efficient. We show examples from synthetic and real offset and walkaway VSP data.In the synthetic examples, we show that data modelled with the estimated wavelet give small residuals when subtracted from the reference data. In the real data examples, we show that when modelling with the wavelet estimated from the real data, in a smooth macromodel, we obtain a good fit between the first arrivals in the real and modelled data.

Notes: An expression which is optimum with respect to the simplicity of the numerical computations is obtained for the magnetic field of a polyhedron with constant magnetization. The high accuracy of the results is illustrated using a realistic numerical model.The existence of the magnetic field at points inside the source and on its boundary is discussed and related to real magnetic data modelling.

Notes: The posterior probability density function (PPD), σ(m|dobs), of earth model m, where dobs are the measured data, describes the solution of a geophysical inverse problem, when a Bayesian inference model is used to describe the problem. In many applications, the PPD is neither analytically tractable nor easily approximated and simple analytic expressions for the mean and variance of the PPD are not available. Since the complete description of the PPD is impossible in the highly multi-dimensional model space of many geophysical applications, several measures such as the highest posterior density regions, marginal PPD and several orders of moments are often used to describe the solutions. Calculation of such quantities requires evaluation of multidimensional integrals. A faster alternative to enumeration and blind Monte-Carlo integration is importance sampling which may be useful in several applications. Thus how to draw samples of m from the PPD becomes an important aspect of geophysical inversion such that importance sampling can be used in the evaluation of these multi-dimensional integrals. Importance sampling can be carried out most efficiently by a Gibbs' sampler (GS). We also introduce a method which we called parallel Gibbs' sampler (PGS) based on genetic algorithms (GA) and show numerically that the results from the two samplers are nearly identical.We first investigate the performance of enumeration and several sampling based techniques such as a GS, PGS and several multiple maximum a posteriori (MAP) algorithms for a simple geophysical problem of inversion of resistivity sounding data. Several non-linear optimization methods based on simulated annealing (SA), GA and some of their variants can be devised which can be made to reach very close to the maximum of the PPD. Such MAP estimation algorithms also sample different points in the model space. By repeating these MAP inversions several times, it is possible to sample adequately the most significant portion(s) of the PPD and all these models can be used to construct the marginal PPD, mean) covariance, etc. We observe that the GS and PGS results are identical and indistinguishable from the enumeration scheme. Multiple MAP algorithms slightly underestimate the posterior variances although the correlation values obtained by all the methods agree very well. Multiple MAP estimation required 0.3% of the computational effort of enumeration and 40% of the effort of a GS or PGS for this problem. Next, we apply GS to the inversion of a marine seismic data set to quantify uncertainties in the derived model, given the prior distribution determined from several common midpoint gathers.

Notes: The long-wavelength propagation and attenuation characteristics of three geological structures that frequently occur in reservoir environments are investigated using a theoretical model that consists of a stack of fine and viscoelastic plane layers, with the layers being either solid or fluid. Backus theory properly describes fine layering and a set of fluid-filled microfractures, under the assumption that interfaces between different materials are bonded. The effects of saturation on wave attenuation are modelled by the relative values of the bulk and shear quality factors.The anisotropic quality factor in a fine-layered system shows a variety of behaviours depending on the saturation and velocities of the single constituents. The wave is less attenuated along the layering direction when the quality factors are proportional to velocity, and vice versa when inversely proportional to velocity. Fractured rocks have very anisotropic wavefronts and quality factors, in particular for the shear modes which are strongly dependent on the characteristics of the fluid filling the microfractures.When the size of the boundary layer is much smaller than the thickness of the fluid layer, the stack of solid-fluid layers becomes a layered porous media of the Biot type. This behaviour is caused by the slip-wall condition at the interface between the solid and the fluid. As in Biot theory, there are two compressional waves, but here the medium is anisotropic and the slow wave does not propagate perpendicular to the layers. Moreover, this wave shows pronounced cusps along the layering direction, like shear waves in a very anisotropic single-phase medium.

Notes: Responses of a multifrequency, multicoil airborne electromagnetic (AEM) system were modelled numerically for 3D electrical conductors embedded in a resistive bedrock and overlain by an overburden of low to moderate conductivity. The results cover a horizontal coplanar coil configuration and two frequencies, 7837 Hz and 51 250 Hz. The models studied are single or multiple, poor conductors (conductance lower than 0.1 S) embedded in a host rock of high but finite resistivity (5000 Ωm) and overlain by a layer of overburden with finite thickness and low to moderate conductivity (conductance up to 2 S).On the basis of the modelling results, limits of detectability for poor conductors have been studied for the various model structures. The results indicate that the anomaly from a steeply dipping, plate-like conductor will decrease significantly when the conductor is embedded in a weakly conductive host rock and is overlain by a conductive overburden. However, an anomaly is obtained, and its magnitude can even increase with increasing overburden conductivity or frequency. The plate anomaly remains practically constant when only the overburden thickness is varied. Changes in overburden conductivity will cause the plate-anomaly values to change markedly. If the plate conductance is less than that of the overburden, a local anomaly opposite in sign to the normal type of anomaly will be recorded. Another major consequence is that conductors interpreted with free-space models will be heavily overestimated in depth or underestimated in conductance, if in reality induction and current channelling in the host rock and overburden make even a slight contribution to the anomalous EM field.The lateral resolution for the horizontal coplanar coil system was found to be about 1.7 times the sensor altitude. Similarly, the lateral extension of a horizontal conductive ribbon, required to reach the semi-infinite (half-space) behaviour, is more than three times the sensor altitude. Finally, screening of a steeply dipping plate, caused by a small, conductive horizontal ribbon, is much more severe than screening of the same plate by an extensive horizontal layer.

Notes: We consider the problem of determining and predicting how the wave speeds in particular directions for a transversely isotropic (TI) medium depend on particular combinations of the density-normalized moduli Aij. The expressions for the qP and qSV velocities are known to depend on four moduli. Normally, we can only determine three independent parameters from qP data, or two from qSZ data, as the others have much lower sensitivity. The resolvable parameters are conveniently described by axial and off-axis parameters: for qP rays, P0°= A11, P90°= A33 and P45°=(A11+ A33)/4 + (A13+2A55)/2; and for qSV rays, S0°= S90°=A55 and S 45°= (A11+ A 33)/4- A13/2. These parameters control the magnitude of the squared-velocities on the axes and at approximately 45°. For an arbitrary TI medium, if the medium is perturbed in a way that preserves a particular parameter, then slowness points in the associated direction and mode witl be approximately preserved in the new medium. we refer to these parameters as ‘push-pins’, i.e. if a parameter is fixed, the associated part of the slowness surface is pinned in place.Because, these five push-pins only contain four independent moduli, we can only fix at most three push-pins. Perturbing one of the other parameters inevitably perturbs the other. Numerical results illustrating the linkage between two push-pins, when three are fixed, are presented.So-called anomalous TI media occur when the roles of the qP and qSV waves are reversed: in some directions the faster ray has transverse polarization. That, in turn, requires anomalous velocities at the push-pins, i.e. S0° 〉 P0°, S45° 〉 P45° and/or S90° 〉 P90° (equivalent to the usual anomalous conditions A11 〈 A55, 〈 0 and/or A33 〈 A55). In the Appendix, we confirm that anomalous sensitivities of the velocities at the five push-pins only occur in such media, although the push-pins still apply if interpreted appropriately. Truly anomalous sensitivities, in which push-pins play no role, only occur in media near the boundary between normal and anomalous.

Notes: Electrical conductivity mapping is a prerequisite tool for hydrogeological or environmental studies. Its interpretation still remains qualitative but advantages can be expected from a quantitative approach. However a full 3D interpretation is too laborious a task in comparison with the limited cost and time which are involved in the majority of such field studies. It is then of value to define the situations where lateral variations are sufficiently smooth for a 1D model to describe correctly the underlying features.For slingram conductivity measurements, criteria allowing an approximate 1D inversion are defined: these mainly consist of a limited rate of variation over three times the intercoil spacing.In geological contexts where the weathering has generated a conductive intermediate layer between the underlying sound rock and the soil, this processing can be applied to determine the thickness of the conductive layer from the apparent resistivity map when the other geoelectrical parameters are known. The examples presented illustrate this application.

Notes: A method to determine the position and magnetization vector of buried objects producing a magnetic anomaly is described. The data used were collected in boreholes. Since the anomaly is due to a number of objects, a ‘stripping’ procedure is employed for finding them, and therefore the process of inversion for finding all objects causing the anomaly consists of a few inversion steps.In each inversion step, two dipoles are considered as a model which approximates an object. The position and magnetic moments of the dipoles are the unknown parameters. The initial parameters are optimized by minimization of an objective function. The optimization procedure consists of a combination of linear and non-linear inversion. The solution of the linear inversion is obtained by singular value decomposition and that of the non-linear inversion by a six-dimensional simplex method (polytope algorithm). After finding one object, its effect is subtracted (‘stripped’) from the data and a new inversion step is started with new initial models and with a reduced data set. The inversion steps for finding different objects are continued until the absolute norm of the data becomes less than some adjustable value.The data will also be inverted assuming a three-dipole model in order to find the effect of using a more complex model in the inversion.The efficiency of the method is demonstrated using synthetic and real borehole data.

Notes: A simple numerical procedure is described for measuring the distance XY in the generalized reciprocal method when there are significant measurement errors in the refracted arrival times. It is applicable when the XY value is of similar magnitude to or less than the receiver spacing. Such conditions frequently occur in using the reciprocal or generalized reciprocal methods to estimate static corrections from first-break times measured in multifold seismic reflection profiling. The use of the method is illustrated with data from both deep and high-resolution seismic reflection profiles.

Notes: An analytical transient solution is obtained for propagation of compressional waves in a homogeneous porous dissipative medium. The solution, based on a generalization of Biot's poroelastic equations, holds for the low- and high-frequency ranges, and includes viscoelastic phenomena of a very general nature, besides the Biot relaxation mechanism. The viscodynamic operator is used to model the dynamic behaviour associated with the relative motion of the fluid in the pores at all frequency ranges. Viscoelasticity is introduced through the standard linear solid which allows the modelling of a general relaxation spectrum. The solution is used to study the influence of the material properties, such as bulk moduli, porosity, viscosity, permeability and intrinsic attenuation, on the kinematic and dynamic characteristics of the two compressional waves supported by the medium. We also obtain snapshots of the static mode arising from the diffusive behaviour of the slow wave at low frequencies.

Notes: The application of McClellan transformations considerably reduces the computational cost of 3D wavefield depth extrapolation by explicit convolutional methods. The accuracy of migration methods based on McClellan transformation depends on how well the transformation filter (cos !;κ!;) is approximated; errors in this approximation cause anisotropy in the extrapolation operator and frequency dispersion in the migrated results. The anisotropy can be greatly reduced by rotating the approximate filter by 45° and averaging the rotated filter with the original filter. The application of the rotated filter yields a migration method that correctly images very steep dips, with little or no additional computational cost. McClellan migration with the improved circular response enhances the imaging of synthetic and real data.

Notes: It is often very useful to be able to smooth velocity fields estimated from exploration seismic data. For example seismic migration is most successful when accurate but also smooth migration velocity fields are used. Smoothing in one, two and three dimensions is examined using North Sea velocity data.A number of ways for carrying out this smoothing are examined, and the technique of locally weighted regression (LOESS) emerges as most satisfactory. In this method each smoothed value is formed using a local regression on a neighbourhood of points downweighted according to their distance from the point of interest. In addition the method incorporates ‘blending’ which saves computations by using function and derivative information, and ‘weighting and robustness’ which allows the smooth to be biased towards reliable points, or away from unreliable ones.A number of other important factors are also considered: namely, the effect of changing the scales of axes, or of thinning the velocity field, prior to smoothing, as well as the problem of smoothing on to irregular subsurfaces.

Notes: A fast imaging technique is developed to deduce the spatial conductivity distribution in the earth from low-frequency (〉 1 MHz) cross-well electromagnetic measurements. A sinusoidally oscillating, vertically orientated, magnetic dipole employed as a source, and it is assumed that the scattering bodies are azimuthally symmetric about the source dipole axis. The use of this model geometry reduces the 3D vector problem to a more manageable 2D scalar form. Additional efficiency is obtained by using the Born series approximation which is derived from nonlinear integral equations that account for the scattered magnetic fields generated by inhomogeneities embedded in a layered earth. Stabilization of the inversion problem is accomplished through the use of bounding constraints and a regularization method which results in a smooth model that fits the data to the desired noise level.The applicability of cross-well electromagnetics for imaging and monitoring changes caused by subsurface processes has been tested by simulating plumes of conductive fluid with 2D models. The images that result from inverting these synthetic data indicate that the vertical resolution of the method is better than the horizontal, increasing the noise decreases the image resolution, and incorporating a priori knowledge in the form of positivity constraints improves the results. Although higher operating frequencies are usually associated with better resolution, frequencies as low as 100 Hz can produce acceptable images in simulated oilfield environments.The imaging scheme has been applied to data collected during a salt-water injection experiment at the Richmond Field Station test site in Richmond, California. Both the data and the resulting images clearly reveal the presence of the plume and indicate that it is migrating towards the north-northwest rather than spreading symmetrically about the injection well. Applying the imaging code to synthetic data generated by a 3D sheet model verifies the interpretation of these results.

Notes: An accurate analytical expression for shot-gather dip-moveout (DMO) in the timespace log-stretch domain has until now not been published. We present a simpler, alternative derivation of the exact DMO relationships of Black et al. which correctly take account of the repositioning of the midpoint. A new computationally efficient frequency-wavenumber (F-K) DMO operator for shot profiles is then derived, based on these DMO relationships in the time-space log-stretch domain. The newly derived DMO operator is, unlike most other log-stretch DMO operators) accurate for the full range of reflector dips. Along with other schemes which are performed in the log-stretch domain, it offers considerable time savings over conventional DMO processing. We have compared numerically the impulse response of the new operator with those of a number of other shot-gather DMO operators, and found it to be superior and well match to the theoretical elliptical DMO response.

Notes: In order to investigate the possibility of using low-frequency electromagnetic waves to detect and monitor oil contamination of soils, a series of laboratory measurements were performed. A new measurement system to monitor the resistivities of soil and sand samples while samples are being contaminated by diesel oil is presented. The frequency used in measurements is 100 kHz. Since the measurement system is composed of coil-type transmitters and receivers, there is no need for electrodes to be in contact with samples. The contamination process is simulated using diesel oil dripping on top of soil and sand samples. The conductivity distributions in samples along the sample length are recorded as a function of time. Water-wet sand and soil samples were measured during diesel oil contamination. The measured data show that the conductivities of soil and sand samples change during the contamination process. The change in resistivity for measured samples before and after diesel oil contamination is in the range of 20% to 50%, giving a reflection coefficient change in the low-frequency limit of 4.7% to 7%. This amount of change in the reflection coefficient makes it very challenging to detect and monitor oil contamination based on EM reflection from the contaminants. The results suggest that EM methods based on propagation and induction, such as tomography and borehole induction, could be used for this purpose.

Notes: A simple expression ties the midpoint of a surface spread to reflection points on a dipping plane. If we use two coordinate systems, an unprimed one with a z-axis perpendicular to the surface and a primed one with a z-axis perpendicular to the reflector, we have 〈displayedItem type="mathematics" xml:id="mu1" numbered="no"〉〈mediaResource alt="image" href="urn:x-wiley:00168025:GPR831:GPR_831_mu1"/〉where θ is the dip angle, φ is the profile angle, X is the source-to-receiver separation, and D is the depth of the reflector. The reflection point is (x, yp, D) and the surface midpoint is (xc, yc, 0).Using the expression, I show that if complete azimuthal coverage is required at a CMP position, the reflection points lie on an ellipse. Similarly, a fixed reflection point generates a circle of surface midpoints. A circle of CMP positions for fixed θ and φ becomes an ellipse of reflection points and a circle of reflection points becomes an ellipse of midpoints. A user can easily find the shape and location of the reflection area generated by a surface aperture.

Notes: An integrated electromagnetic and nuclear magnetic resonance (NMR) method is proposed for investigating highly conductive soil areas. Using a simple model of a homogeneous half-space and the same antenna for both methods, a significant improvement in the NMR data interpretation is obtained. A case study has shown fair agreement between the results from computer modelling, field tests, and data from a nearby observation well. The electromagnetic method and the half-space model were selected for easy integration into an existing instrument used for the NMR method. A more accurate knowledge of the conductivity distribution with depth will further improve the final result.

Notes: A new method is proposed for the quantitative interpretation of SP field data produced by a polarized ore body and simulated by an inclined sheet. The theoretical concept is based on the study of the amplitude Spectrum. It is shown that the SP amplitude Spectrum is not continuous at zero frequency; this leads to the dip angle determination of the inclined sheet. It is also shown that the SP amplitude spectrum is practically nullified at a characteristic amortization frequency that depends on the depth of the polarized body. The maximum amplitude Spectrum value of the SP gradient is used to estimate the depth to the bottom of the polarized body. Thus, the geometrical parameters h and H, the depths to the top and bottom, respectively, as well as the dip angle of the inclined sheet, can be satisfactorily determined. Some problems may arise in the determination of these parameters, affecting their accuracy, whenever unwanted frequency noise is present.

Notes: Twelve ground penetrating radar (GPR) experiments were conducted on the modern, wave-influenced William River delta, on the Southern shore of Lake Athabasca in northern Saskatchewan, Canada. The delta is a well-sorted, quartzoserich, clean, sand-dominated, water-saturated geomorphic feature which provided an ideal site to test GPR. Penetration depths, resolution and continuity of reflections were compared for different antennae frequencies (25, 50, 100, 200 MHz) and transmitter powers (pulser voltage: 400 V, 1000 V). The data show significant variations in vertical resolution from 0.15 m to 0.76 m (200-25 MHz), depth of penetration from 14 m-28 m (200-25 MHz), and continuity of reflections. Increasing the transmitter power from 400 V to 1000 V increases the depth of penetration by 5 to 14% and improves the continuity of reflections with little effect on the resolution.