ALBERT

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: 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.

Notes: The paper discusses the variations of Elevation Correction Factors (E.C.F.) across various stratigraphical and structural occurrences with particular reference to the dipping bed, fault and horst. Graphs are prepared giving E.C.F. variations for various angles of dip and a method is suggested for application in the field. Errors in the linear variation that is commonly used for small angles of dip are discussed.The second part of the paper consists of the application to interpretation of (I) Change of gravity dip and (2) The gravity inversion. It is shown how from a gravity map the dip, position of outcrop and thickness of a dipping bed can be determined. The last part shows how anticlines can give rise to negative anomalies and that in the direction of dip similar structures can at one point have no surface expression and afterwards invert to positive anomalies. The danger of attempting to predict the size of structures in such areas is emphasised.

Notes: When gravimetric prospecting is concerned with a region featuring broken relief the calculation of topographical corrections becomes vitally important, for these rectify apparent anomalies completely. It is essential that they should be calculated as exactly as possible. By the classic method this takes an extremely long time and may in the long run involve an outlay actually higher than that expended on operations on the terrain.It has thus been felt necessary to seek methods which would speed up calculation, and, bearing that in mind, the use of a polar integrator has been envisaged. With this the contour lines of a map can be followed between two fixed radii, the pole corresponding to the location of the position. Two integrator roulettes will give the values of the first two terms of the expansion in terms of the quotient of the height divided by the distance, the second of the terms resulting from the curvature of the earth. Unless there are escarpments very close to the position, the use of these two terms will be sufficient.Polar integration will deal successively with different contour lines and a subsidiary graph will be plotted, giving the value of each of the integrals in terms of the height. By means of this subsidiary graph the values of the corrections will be obtained with a planimeter and a moment integrator. It is envisaged that the same operation could be repeated in the case of two maps with different scales, for instance, 1–20,000 and 1–200,000 the effect of the areas in the immediate neighbourhood of the position being taken into account by the observer. It is hoped that in this way the calculation of topographical corrections will be made quicker, more exact and less laborious.

Notes: The problems connected with reflexion surveys in areas of rough topography consist of choosing proper equipment and of making careful near-surface corrections. This paper illustrates some solutions to these problems that have been adopted in a survey which has recently started in a hily coastal area, in Southern Italy.As the area is characterized by a youthful and very actively eroded topography of Post-Miocene clays and litoral clastic sediments, both rotary and percussion drilling were necessary whilst all the seismic equipment had to be light to allow straight lines to be shot as far as possible. Details are given of the truckmounted, light equipment which has been designed for this special purpose.In connection with the youthful topography of the area being surveyed, the criteria for the proper computation of the near-surface corrections is discussed. The great differences in elevation, combined with variations in the surface and near-surface materials, necessitated a careful analysis before deciding the thickness of the correction zone. The results are shown of the preliminary tests carried out in order to determine the average thickness of the low-velocity zone. The solution adopted for drawing the reference surface determining the base of the correction zone for the whole of the area is subseqently explained.Reference datum is drawn so as to follow gently, and with segments of constant slopes, the general features of the rough topographic relief at a depth ranging from 30 to 60 metres from the surface. Criteria for the routine checking of velocity variations in the correction zone are illustrated. Such checks are based on both uphole time analyses and refraction “first breaks” plots, with deeper shot holes at constant intervals.Unsatisfactory results necessitate longer spreads for “low velocity tests” and uphole shooting from deeper holes and the reference datum is deepened for the particular area where this becomes necessary.The final cross-sections are shown and although, as yet, no deep holes, are available to check the seismic data, the general appearance of the cross-sections seems to indicate that subsurface data are free from the influence of topography and of the velocities in the near-surface formations.

Notes: The paper describes the determination of the chlorine content of deep ground water by means of the direct resistivity method using the Wenner electrode configuration. The purpose was to investigate the distribution of salinity in the brackish waters within a permeable sand of thickness 250 metres and which forms the deeper subsoil of the Ysselmeer areas.Such a knowledge is of great value, firstly because the seepage of salt water into the proposed deep “polders” may reduce the crop yield, and secondly because it allows areas, where potable waters may occur, to be delineated.To test the method, measurements were made near existing boreholes which extended to a depth of 350 metres. From these borings water samples from various depths were collected and it was found that the resistivity pl of the porewater was closely related to their chlorine content. The resistivity p of the saturated sand was determined from the geo-electrical measurements. It was found that the resistivity p was given by p= 4, 5 pl. From this relation, which is independent of the salinity, the chlorine content of the porewater could be estimated. These results are in close agreement with those of other observers. In the area p varied from 1 to 200 ohmmetres corresponding to a chlorine content changing from 12.000 to 15 mg per litre.The method was adapted for working over water. Each of fourteen insulated conductors in a special cable terminated in an electrode, and these electrodes were spaced along the cable so that nine different Wenner spacings could be selected. The electrode separations varied from 44 metres up to 306 metres.

Notes: The authors develop an optimum filtering theory for seismic reflection recording. A criterion of record improvement is applied to a mathematical model consisting of a set of wavelets superimposed on a random noise backgruond. This leads to the design of electrical filters and geophone and shot-hole patterns with optimum filtering characteristics.

Notes: For horizontal layering and plane horizontal impulsive wavefronts it is theoretically possible to get rid of multiple reflections by a feedback procedure which can easily be derived using raypath philosophy. To reduce the increase of noise inherent in the method the precise theoretical formulae are altered in such a manner that a practical application becomes possible. For this purpose the autocorrelation will be used.Application of the new process to a CRP section where the shot geophone distances in the field had not been long enough to attenuate multiple reflections effectively gave favourable results.

Notes: Two in situ measurements of stress wave velocity in a granodiorite massif were made. Firstly, the velocities close to the source were investigated (for elimination of the influence of cracks), and secondly, the velocities in different horizontal directions over a large area were studied from the viewpoint of their dependence on the crack system of the medium. A correlation between the mutual orientation of the cracks and the maximum or minimum values of the stress wave velocities was obtained.

Notes: The technique of linear digital filtering as developed for the direct interpretation of Schlumberger and Wenner soundings (Ghosh 1971) has been applied here for the derivation of the resistivity transform function from the field dipole measurements as the first step in directly interpreting dipole data. Filter coefficients for this transformation have been worked out for the radial-polar, perpendicular and parallel (30°) arrays of dipole sounding. The procedure combines speed with accuracy.

Notes: Dipole soundings are more sensitive to noise caused by lateral and superficial inhomogeneities than Schlumberger soundings. However, the former are preferable for deep explorations in view of the relatively short cables required. The simple solution of carrying out the field work by means of dipole spreads, and to transform the dipole resistivity diagrams into Schlumberger ones by means of proper formulae would be valid only for smooth and regular curves; but often, owing to the presence of lateral noises, the dipole data show a considerable scatter. For such cases a “continuous dipole sounding” method is proposed for which all successive dipoles are contiguous, so that all parts of the profiles are covered and interpolation is not necessary. Obviously the moving dipoles have lengths proportional to their distances, so that they appear equal in the usual bilogarithmic scale. It follows that only polar-dipole arrays may be used. The transition from a dipole to the corresponding Schlumberger apparent resistivity diagram requires an integration constant which is not unequivocally determined. Therefore, the solution is not unique, but all possible derived Schlumberger diagrams have a common part. Similarly, they have some common interpretative results, which may be referred to the original dipole diagram obtained in the field.A special measurement technique is required since the dipole-dipole voltages to be determined are noticeably smaller than the Schlumberger ones. This is true also because dipole soundings are used for great depths and for long distances between the two dipoles.

Notes: Basic requirements for a seismic reflection system using a moderate-power controlled source are discussed. General considerations and computer modelling show promise for geophysical prospecting, especially at low penetration and very high noise level. For many applications, ranging from mineral and water prospecting, to civil and military engineering, and even archeology, a relatively cheap, small and portable device of this kind might be more convenient than conventional explosive sources or sophisticated controlled sources of the vibroseis type. The use of a simple, low frequency periodic signal facilitates rapid data processing procedures and suggests that the mechanical generator will be very simple and cheap.

Notes: A method is described to transform a dipole sounding curve, obtained with any one of the common dipole arrays over a horizontally layered earth, to the form of a Schlumberger sounding curve. Starting from the general expression which relates the dipole apparent resistivity to the Schlumberger apparent resistivity and its derivative with respect to the spacing, it is possible with some approximations to derive an easy numerical computation procedure in order to perform the transformation. The applicability of the method is discussed briefly.

Notes: Approximate deconvolution by means of Wiener filters has become standard practice in seismic data-processing. It is well-known that addition of a certain percentage of noise energy to the autocorrelation of the signal wavelet leads to a filter that does not increase, or even reduces, the noise level on the seismogram. This noise addition will, in general, cause a minimum phase signal to become mixed phase. A technique is presented for the calculation of the optimum-lag shaping filter for a contaminated signal wavelet. The advantages of this method over the more conventional approach are that it needs less arithmetic operations and that it automatically gives the filter with the optimum combination of shaping performance and noise reduction.

Notes: A set of parabolae is defined to approximate the observed gravity profile for specifying the lateral density variation in a two dimensional causative body. This variation separates the “biased” as well as “unbiased” residual anomaly structure. The method may be directly incorporated into existing numerical and graphical techniques of interpretation for subsequent delineation of structural configuration. The simple relations derived by inspection of the observed gravity profile permit universal application of the technique and in particular for the determination of the size and shape of sub-shelf structures. The suitability of the method in various problems of two dimensional interpretation has been indicated through an illustration of the Godavari Basin (India).

Notes: Decomposition of IP discharge curve into a sum of exponentials yields more information than present practice does. Such a work was performed, as a beginning, by hand, plotting the ordinates on semi logarithimic paper. It gave, according to the cases, two or three exponentials, each with an amplitude Ai and the time constant τi.Three examples of ground prospection are given, the third of which shows an anomaly, unnoticed with present technique but obvious on a profile curve with A1/A2 in ordinates.A computer processing technique is described which yields directly the time constants τi by seeking a differential equation which the IP voltage satisfies; the order of the equation is set arbitrarily. On a practical prospection record with only two exponentials, the assumption that the differential equation does not depend on charge duration (everything else being equal) has been checked (the shape of the curves would depend on boundary conditions). The assumption roughly holds for charge durations from 5 seconds to 40 seconds. We hope that a larger number of exponentials will yield better results.〈section xml:id="abs1-2"〉〈title type="main"〉SommaireLa décomposition de la courbe de décharge en une somme d'exponentielles permet d'en tirer plus d'information que la pratique actuelle. Un tel travail a été fait, pour commencer, à la main en portant les ordonnées sur papier logarithmique; il a fourni, suivant les cas, deux ou trois exponentielles avec, pour chacune d'elles, le coefficient Ai et la constante de temps τi.On donne trois exemples de prospection sur le terrain, dont le dernier comporte une anomalie, inaperçue avec la technique actuelle, mais bien mise en évidence en traçant la courbe de A1/A2 le long du profil.On donne ensuite une technique de calcul sur ordinateur des τi, en recherchant l'équation différentielle à laquelle obéit la tension ΔVPP enregistrée, l'ordre de cette équation étant fixé arbitrairement. On a cherchéà vérifier sur un enregistrement de terrain, avec. deux exponentielles seulement, l'hypothèse suivant laquelle l'équation différentielle ne dépend pas de la durée de charge (toutes choses égales d'ailleurs). L'hypothèse se trouve grossièrement vérifiée pour les durées de charge allant de 5 secondes à 40 secondes, et nous espérons qu'un plus grand nombre d'exponentielles donnera de meilleurs résultats.

Notes: The conventional seismic technique is subject to a recording time following each transmission of energy, in which it is forbidden to release any new pulse. The recording time depends on the deepest reflection of interest, and is often 10 seconds or more in actual practice. To each transmission corresponds one record, i.e. a fixed amount of data which cannot be increased in a given time.Pulse coding allows us to go beyond this limit, by transmitting several times during the normal recording time. The procedure gives as many records as there are pulses, but they overlap, each event being repeated every time there is a pulse. It is possible to process the composite record back to its usual appearance with all events in their proper place if the time breaks are accurately known and make up a code such that the unavoidable noise generated by the process be kept, on the final section, below the ambient noise. The processing is quite similar to that of records made from vibrating sources, though faster in practice.The additional information can be devoted to a saving of time and money as the same profile may be recorded in a shorter time; or to an improvement of quality of the section due to a higher order of coverage, a multiplication of the ray paths and a closer sampling of the reflectors. It is also possible to record information in several planes at the same time, and to work out a 3-dimensional restitution, without loss of production.The process applies to all kinds of sources provided they can be triggered according to the code with sufficient accuracy. Depending on the source and conditions of implementation, the method benefits from other advantages such as better resolution, increased flexibility, and better coupling.Two different names have been given to the process, Sosie and Seiscode, which apply to slightly different parameters for the sequence of pulses. Sosie is more useful at sea, while the normal scope for Seiscode is onshore. Both names are trademarks for SNPA.

Notes: Master diagrams for electromagnetic responses of nine dipole dipole systems are presented for various depths of a vertical infinitely conducting vein. The YY system gives the minimum anomaly for both the inline and broadside arrays. Among other inline systems, it is difficult to decide the clearcut superiority of one system over the other, whereas the XX system in broadside array gives maximum anomaly.

Notes: The interpretation of Turam data is treated by means of a development of the conventional approach to this type of data. The development consists of a curve matching technique based upon a line current approximation to the true induced current pattern.The technique requires computation facilities which are readily available in modern programmable calculators.Practical examples of the use of the method illustrate that it can reveal information concerning the geometry of a conducting body which would not be gained by the conventional approach.

Notes: Two-layer type curves of apparent frequency effect for the Wenner configuration are presented. The formulation is based on the normal definition of frequency effect in terms of resistivities measured at different frequencies plus the definition of apparent resistivity over two horizontal layers as a function of first and second layer resistivities. The use of these type curves in the interpretation of multilayer apparent frequency-effect curves is described and some field examples are given.

Notes: In this paper properties of the discrete zero-phase time function are derived and compared with related properties of the discrete minimum-phase time function.The two-sided minimum-length signal is introduced and it is derived that, for any given amplitude spectrum, the two-sided minimum-length signal and the signal with zero-phase spectrum are identical signals. A comparison is made between the one-sided minimum-length signal (minimum-phase signal) and the two-sided minimum-length signal (zero-phase signal).A computational scheme is discussed which determines the zero-phase correspondent of a given signal.A method is proposed to compute zero-phase least-square inverse filters. The efficiency of minimum-phase and zero-phase least-square inverse filters is shown on signals with different phase properties.A criterion is derived which determines whether a symmetric time function has the zero-phase property. The close relationship with the minimum-phase criterion is discussed.Finally the relationship between signal length and resolving power is illustrated on numerical examples.

Notes: Empirical equalities derived from time domain induced polarization scale modelling with the gradient array over simple geometries, and from the potential field functions for equivalent simple charge configurations are similar. The function for the dyke like body is analogous to the magnetic case allowing both total and vertical magnetic field interpretation techniques to be applied to gradient array chargeability anomalies.

Notes: “Remote sensing” techniques have been used for many years as an exploration tool; broadly speaking, any method using sensors not directly connected to the earth's surface falls into this group.It is preferable, however, to use the term solely with reference to a number of new sensors which have recently become available to earth scientists for the detection of various radiations of the electromagnetic spectrum in the wave length range from 0.1 micron (UV) to some centimeters including the visible and infrared regions.Radiations artificially induced in the earth (active sensing) or spontaneous radiations (passive sensing) can be detected. Some of the bands investigated seem promising for application to geological surveys, especially using sensors on aircraft and satellites, which allow synoptic and large scale investigations.A brief account is given of the results already achieved and of future developments such as the ERTS and manned satellite (post Apollo) “earth observation” programmes.The potential role of Europe in the field is discussed.Details are given of some applications of thermal I.R. imagery to geological and vulcanological problems, with special emphasis on the research carried out in Italy.

Notes: A simple calculation procedure that gives the equilibrium temperature in geothermal measurements reduces the reading time in field work from reported 60 minutes to about 10 minutes.

Notes: Boundary-value problems in steady-state current flow were solved numerically in bispherical coordinates for a sphere of arbitrary conductivity in a half-space. Solutions for the potential on the surface of the half-space were examined for the cases where both current sources were on the surface, one source on the surface and the second between the surface and the sphere, and one source on the surface and the other in the sphere. Results show a great similarity with the layered case when the buried electrode is placed between the surface and the conducting region. Such a buried electrode configuration makes it possible to obtain an accurate measurement of the depth to the conductor in both cases. A model with the current electrode placed in a conductive sphere is compared with a three-layered model with the source in a conductive intermediate layer, and results indicate that the lateral extent of a finite anomalous zone can be estimated using these limiting curves. The validity of these theoretical calculations for buried spheres was confirmed experimentally by tests conducted on an analog model.

Notes: A checking method of digital multiple elimination and of deconvolution processing using computers and based on optical autocorrelation is first presented. Comparison between autocorrelograms before and after a single or several processing steps allows to estimate, on one hand, the strength of the deconvolution obtained, known by the study of the central parts which is in fact the signal autocorrelation, on the other hand, the multiple elimination given by the study of side parts of the autocorrelogram.Further, an optical deconvolution procedure, is presented. For this, it is supposed that the signal is known and optically reproduced in the same way as the one of a trace. This is achieved by sphero-cylindrical optics allowing trace to trace processing. Deconvolution is carried out in the spectral domain by inserting a filter in the Fourier plane of the optical unit, the transmission law of which expresses the Fourier transform of the antisignal. This filter device introduces a holographic technique called Fourier holography, in such a way phases as well as amplitudes are preserved.Several results are presented from a synthetic section and also from field sections.

Notes: The separation of regional from local gravity anomalies by means of the application of two-dimensional linear filters is analyzed. It was found that optimization of the filter in the least squares sense leads to filters that produce strong localized concentrations of the error, which may erroneously be interpreted as anomalies. For this reason the maximum absolute value of the error is a more important criterion for the quality of the filter than the root mean square error. This maximum absolute error is minimized by the minimax filter. Intermediate filters are derived which give a transition zone which comes appreciably closer to that of the optimal filter at only a small price in terms of increase of the maximum absolute error.

Notes: The Laboratory of Geophysics of the University of Arizona was presented with an exploration problem by Falconbridge, Ltd. of Canada. Massive mineralized hemispherical “pods” are in the vicinity of a tuff layer of high conductivity and induced polarization response, covered by large thicknesses of resistive volcanics. The initial approach was to utilize electrolytic tank modeling. The extreme resistivity and IP contrasts proved to be difficult to recreate. Two dimensional modeling was attempted next with conductive paper, using copper and silver paint for anomalous masses. This method also proved inadequate. Finally, mathematical equations were solved which could model any arbitrary anomalous body in any steady state electrical field. Plane waves as well as point current sources producing non-plane waves are possible. Finite difference equations were derived for the non-linear partial differential equations under consideration. The equations were solved using a digital computer. Initially, the boundary conditions had to be satisfied at the boundaries of resistivity changes, severely restricting possible geometric shapes for anomalous bodies. The final and successful solution was to apply numerical techniques to obtain solutions of equations which require only that the relative resistivities through the area be specified. The Falconbridge problem and its solution are analyzed.

Notes: For the two and three layer cases geo-electrical sounding graphs can be rapidly and accurately evaluated by comparing them with an adequate set of standard model graphs. The variety of model graphs required is reasonably limited and the use of a computer is unnecessary for this type of interpretation.For more than three layers a compilation of model graphs is not possible, because the variety of curves required in practice increases immensely. To evaluate a measured graph under these conditions, a model graph is calculated by computer for an approximately calculated resistivity profile which is determined, for example, by means of the auxiliary point methods. This model graph is then compared with the measured curve, and from the deviation between the curves a new resistivity profile is derived, the model graph of which is calculated for another comparison procedure, etc. This type of interpretation, although exact, is very inconvenient and time-consuming, because there is no simple method by which an improved resistivity profile can be derived from the deviations between a model graph and a measured graph.The aim of this paper is, on the one hand, to give a simple interpretation method, suitable for use during field work, for multi-layer geo-electrical sounding graphs, and, on the other hand, to indicate an automatic evaluation procedure based on these principles, suitable for use by digital computer.This interpretation system is based on the resolution of the kernel function of Stefanescu's integral into partial fractions. The system consists of a calculation method for an arbitrary multi-layer case and a highly accurate approximation method for determining those partial fractions which are important for interpretation. The partial fractions are found by fitting three-layer graphs to a measured curve. Using the roots and coefficients of these partial fractions and simple equations derived from the kernel function of Stefanescu's integral, the thicknesses and resistivities of layers may be directly calculated for successively increasing depths.The system also provides a simple method for the approximative construction of model graphs.

Notes: The theoretical horizontal resistivity profiles over an outcropping vertical dyke for various parameters-electrode spacing, vein-width and resistivity contrast—with inline alpha, beta and gamma-Wenner electrode systems are described. The resistivity profiles present a most bewildering variety of shapes as compared to those in resistivity soundings. The analysis of type curves suggests suitable electrode configuration for detection of wide, moderately wide, and thin veins. The negative apparent resistivity values on the gamma anomalies for resistive vein over certain vein-widths of higher positive values of resistivity reflection factor contradict the concept of apparent resistivity hitherto held.