ALBERT

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 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: 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: This is the first of a series of papers giving the solution of the inverse problem in seismic exploration. The acoustic approximation is used together with the assumption that the velocity field has the form 〈displayedItem type="mathematics" xml:id="mu1" numbered="no"〉〈mediaResource alt="image" href="urn:x-wiley:00168025:GPR998:GPR_998_mu1"/〉 . The forward problem is then linearized (thus neglecting multiple reflected waves) and the inverse problem of estimating δ is set up. Its rigorous solution can be obtained using an iterative algorithm, each step consisting of a classical Kirchhoff migration (hyperbola summation) plus a classical forward modeling step (circle summation).

Notes: The response of a seismic model to excitation by a source can be represented in terms of the action of reflection and transmission operators for portions of the structure. This approach provides a flexible framework for both modeling and processing problems.The operator development provides a physical description of the wave propagation process and, via the expansion of reverberation operators, gives a mechanism for assessing the accuracy of approximate developments. The representation suggests new ways of developing modeling algorithms by balancing the computational effort expended on minor and major features of the model.For processing problems, the operator representation shows the relation of processing stages to the seismic wave field and thereby indicates effective sequences of operations. For migration it is possible to specify an ideal pre-stack migration procedure in terms of the inverse of the propagation operators and to examine the problems which need to be overcome by practical algorithms.

Notes: A numerical method is presented for direct interpretation of resistivity sounding measurements. The early part of the resistivity transform curve derived from field observations by standard methods is approximated by a two-layer curve. The resistivity of the first layer is determined from the arithmetic mean of the successive computations which are carried on each of three successive discrete values of the resistivity transform curve. Using this mean value of the resistivity, the thickness of the first layer is computed from the sample values in pairs of the resistivity transform curve. After these determinations, the top layer is removed by Pekeris's reduction equation. The parameters of the second layer are obtained from the discrete values of the reduced transform curve (which corresponds to the second part of the resistivity transform curve) by the same procedure as described for the first layer.The same computational scheme is repeated until the parameters of all intermediate layers are obtained. The resistivity of the substratum is determined from the reduction equation.

Notes: In odd-depth structure the two-way traveltime to each boundary is constrained to be an odd integer. The odd-depth property of a model is exposed to possible refutation under a seismogram test. Test function is a simple transformation of a synthetic seismogram. For an odd-depth model the test function has identically the value 1.The testability of a synthetic seismogram over an odd-depth structure provides a method of deterministic deconvolution. There is no need of specialized assumptions, like the minimum-phase property, about the source wavelet. The deconvolution may be performed in the absence of the early segment of a seismogram.

Notes: Zero-offset-source VSP surveys provide information about the subsurface only within the Fresnel zone centered at the well. Offsetting the source location moves the reflection zones away from the well thus providing lateral cover.Conventional processing of this type of data gives rise to a distorted image of the subsurface. Using a simple ray-tracing scheme, this image may be reconstructed into the more familiar coordinate system of the surface seismic section. This simple data-independent mapping is based on the assumption of horizontal layering and requires a vertical velocity profile.The technique of placing the source away from the borehole was first applied to the single-offset-source VSP survey. However, data from any survey geometry (such as deviated well with rig source, walkaway VSP, etc.) can be mapped to the coordinate system defined by the appropriate seismic section.To obtain the best results from this type of survey the target area must be defined and simple modeling techniques used to optimize the source location(s). These pre-survey modeling methods may also be used to anticipate—and hence avoid a number of problem areas which experience has highlighted.The data from any VSP survey is the result of a realizable experiment and as such obeys the wave equation. This implies that the wave equation may be used to migrate the data to its true subsurface location. Theoretically, such a process is more secure than ray-tracing techniques, although its practice presents many difficulties.

Notes: A mathematical expression for potential of a direct current point source in an inhomoge-neous anisotropic earth is derived. The coefficient of anisotropy is given by f= (σr/σz), where σr and σz are the conductivities parallel and perpendicular to the bedding plane. It is assumed that σz varies with depth, whereas σr varies transversely. This potential may be useful in interpretation of geoelectrical data in specified geological situations. Master curves for Wenner and Schlumberger configurations are presented

Notes: Total magnetic intensity contour maps for the study region (between 2°E to 10°E and 56°N to 60°N) were digitized and converted to a regular grid of 285 × 285 points. The study area measures approximately 444 km × 444 km and the grid spacing is thus 1. 56 km. The International Geomagnetic Reference Field for 1975 was gridded for the above-used net, and from the two data sets a further grid of the ▵T field was generated. A large number of profiles were constructed which were suitable for depth determinations. The regular grid ▵T data is also convenient for the computation of the second vertical derivative. Using the method of vertical prisms of Vacquier et al. (1963), a large suite of curvature-depth indices was measured to complement the depths obtained from the intensity slopes and from boreholes which reach the crystalline basement. The depth to the magnetic basement has been contoured, and the resulting map is shown to be in good agreement with what is known about the deeper geology of the study area.The work reported here is part of a research project supported by Amoco Norway, BP Petroleum Development Ltd, Elf Aquitaine, Esso Exploration and Production, Norwegian Gulf, Norsk Hydro, Mobil Exploration Norway, Norwegian Petroleum Directorate, Royal Norwegian Council for Scientific and Industrial Research (NTNF), Norske Shell, and Statoil.

Notes: Geophysical inversion involves the estimation of the parameters of a postulated earth model from a set of observations. Since the associated model responses can be nonlinear functions of the model parameters, nonlinear least-squares techniques prove to be useful for performing the inversion. A common type of inversion applies iterative damped linear least squares through use of the Marquardt-Levenberg method. Traditionally, this method has been implemented by solving the associated normal equations in conventional ways. However, Singular Value Decomposition (SVD) produces significant improvements in computational precision when applied to the same system of normal equations. Iterative least-squares modeling finds application in a wide variety of geophysical problems. Two examples illustrate the approach: (1) seismic wavelet deconvolution, and (2) the location of a buried wedge from surface gravity data. More generally, nonlinear least-squares inversion can be used to estimate earth models for any set of geophysical observations for which an appropriate mathematical description is available.

Notes: In this study we derive expressions for particle displacement or particle velocity anywhere inside a stratified earth and at its surface due to horizontal torque source located in the top layer. Equivalently, invoking Green's function reciprocity theorem, the solution applies also to the case of a surface or subsurface source when the resulting displacement or velocity is measured within the top layer.In order to evaluate the closed-form analytical solution economically and accurately it is advisable to introduce inelastic attenuation. Causal inelastic attenuation also lends the necessary realism to the computed seismic trace. To provide proof that the analytical solution is indeed correct and applicable to the multilayer case, a thick uniform overburden was assumed to consist of many thin layers. The correctness of the computed particle velocity response can be very simply verified by inspection. The computed response can also serve as a check on other less accurate methods of producing synthetic seismograms, such as the techniques of finite differences, finite elements, and various sophisticated ray-tracing techniques.It is not difficult to construct horizontal surface torque source. It appears that such source is well suited for seismic exploration in areas with a high-velocity surface layer. A realistic source function is analyzed in detail and normalized displacement response evaluated at different incidence angles in the near and the far fields.In an effort to distinguish the features of an SH torque seismogram from a pressure seismogram two models with identical layerings and layer parameters have been set up. As expected the torque seismogram is very different from the compressional seismogram. One desirable feature of a torque seismogram is the fast decay of multiples.Exact synthetic seismograms have many uses; some of them, such as the study of complex interference phenomena, phase change at wide angle reflection, channeling effects, dispersion (geometrical and material), absolute gain, and inelastic attenuation, can be carried out accurately and effortlessly. They can also be used to improve basic processing techniques such as deconvolution and velocity analysis.The numerical evaluation of the analytical solution of the wave equation as described in this paper has a long history. Most of the work leading to this paper was carried out by one of us (M. J. K.) in the years 1957 to 1968 at the Geophysical Research Corporation. However, the full testing of the various computer codes was carried out only very recently at the Phillips Petroleum Company.

Notes: In South and Central Goa iron ore occurs in two parallel belts with the general NW-SE Dharawar trend. The ore occurrence, however, is not continuous. There are barren zones as well as zones of very high concentration in some of which there are mining activities.Landsat MSS data have been interpreted over a zone covering both mineral belts in order to delineate the ore occurrences. As a guide line a known ore-bearing area has been considered along with the unknown zones.On the basis of two-dimensional plotting of gray level values it has been found that the MSS bands 4 and 7 are most suitable for the studies over iron-rich areas in Goa. Two techniques are described here for the processing of the MSS data; the separation of residual from the regional and MSS band-ratioing. It is observed that (i) the gray level residual maps of MSS bands 4 and 7 are of use in demarcating the iron-ore-bearing zones, and (ii) an existing mine, an abandoned mine, and a proved iron ore zone could be delineated by MSS band-ratioing. On the basis of the latter technique, a few areas with ore occurrence potential have been indicated.

Notes: The most common source of seismic energy is an explosion at some depth in a borehole. The radiated waves are reflected not only at the subsurface layers but also at the free surface. The earth's surface acts as a generator of both P- and S-waves.If the source depth is much less than the dominant wavelength the reflected waves resemble closely the waves generated by a single force. Theoretical seismograms were computed with different methods to look for the relevance of the surface-reflected waves. The numerical experiments show reflected shear waves even for small shotpoint—receiver distances. Due to their polarization these waves can be detected most easily on in-line horizontal geophones. The existence of these waves was examined during a conventional survey in Northern Germany. Conventional data analysis shows a large variability in the νp/νs ratio. The method used here produced a shear-wave section with a rather good signal-to-noise ratio down to 4 s S-wave reflection time.

Notes: We describe the concept of physico-geological models (PGMs) in geophysical exploration. They represent a “general model”, a spatial combination of a set of particular models (disturbing bodies). The modeling is called complete, incomplete or approximate, depending on the degree of characterization of the PGM by parameters such as dimension, shape and petrophysical property. Each of the three modeling types can be realized as a conceptual, and analytical, or a material PGM. Both deterministic and stochastic PGMs exist; deterministic models are mainly used to investigate the possibilities of a geophysical method, while stochastic models serve to substantiate complex geophysical interpretations.Depending on the geological problem, PGMs are subdivided into multi-alternative models (geological mapping, prediction, general prospecting) and double alternative models (specialized prospecting).An exploration-oriented classification of the PGMs of mineral deposits is discussed. According to this classification the variety of known genetic deposit types is reduced to a limited number of generalized PGM types. The development of typical PGMs is illustrated with examples of magnetitic deposits of Siberia.

Notes: An inverse problem is one in which the parameters of a model are determined from measured seismic data. Important to the solution of inverse problems is the issue of whether or not a solution exists. In this paper we show, in a constructive manner, that a solution does exist to the specific inverse problem of determining the parameters of a horizontally stratified, lossless, isotropic and homogeneous layered system that is excited by a non-normal incidence (NNI) plane wave. Mode conversion between P- and S-waves is included.We develop a seven-step layer-recursive procedure for determining all of the parameters for layer j. These parameters are P-wave and S-wave velocities and angles of incidence, density, thickness, traveltimes, and reflection- and transmission-coefficient matrices. Downward continuation of data from the top of one layer to the top of the next lower layer is an important step in our procedure, just as it is in normal incidence (NI) inversion. We show that, in order to compute all parameters of layer j, we need to (and can) compute some parameters for layer j+ 1. This is a non-causal phenomenon that seems to be necessary in NNI inversion but is not present in NI inversion.

Notes: The suitability of VLF ground surveys in the investigation of shallow two-dimensional structures is analyzed. For such structures the polar formalism is derived, necessary in practice since the transmitters are generally not in the structural strike or profile. A simple vertical dike is considered to demonstrate the striking anisotropy which can be expected over such a structure, in particular the high apparent resistivity along the direction of a well-conducting dike and the low resistivity across it. The theory is then confronted with the practical example of an asymmetrical vertical dike resulting from a strike-slip fault. Modelling of the survey results is very successful and yields good confirmation of the polar behaviour. VLF ground surveys thus provide a quick and powerful tool for the study of geological accidents within about 100 m of the surface.

Notes: A new mode of operation for the Turam electromagnetic exploration system is proposed in which the transmitter loop is placed across the expected trend of a conductor and the receiver is operated along lines parallel to one side of the transmitter. The concept appears to offer several benefits which include greatly extended traverse length, the use of large coil spacing, rejection of the effects of conductive environments, and consistency in the indication of target dip.

Notes: Determination of thickness of sediments (usually of high conductivity) overlying a high-resistivity basement is one of the basic problems of electrical exploration methods. This paper proposes to determine horizontal electrical conductance on the basis of impedance calculated from electrical and magnetic fields of distant quasi-static (low-frequency) point sources. Using the proposed method, horizontal conductance of the sediments can be determined also from artificial quasi-static noise-impulses coming from sources of unknown position and intensity. The results of analogue modeling and field examples prove the potential of the proposed technique.

Notes: Since the early years of seismic surveying, field engineers and observers have been faced with the appearance of interference at power-line frequency (and harmonics) on seismic records: 60 Hz in the Americas, 50 Hz in Eurasia; and 16 2/3 Hz (the train supply frequency in Germany, Austria, and Switzerland) in part of Europe. Such interference is picked up from the ground surface through electric leakage between geophone leads and ground.To minimize power-line interference, the autobalancer automatically performs a balancing procedure by adjusting two potentiometers on each channel—essentially the same operation as manually performed by seismic operators for decades. The first 24-channel system became operational in 1977.The first-generation design comprised two 256-position electronic switching potentiometers. Over 4 years of experience with this system in all parts of the world has shown its value under a variety of field conditions. Improvements up to 30 dB were often found, sometimes even up to 40 dB. Balancing time was about 30 s.A second-generation design was field-tested in 1981. Potentiometer stepping is now performed at three levels: coarse, medium, and fine. This leads to faster operation and better resolution: balancing time is about 10 s; improvements up to 60 dB are occasionally found.A later modification gives a further reduction of the balancing time to 5 s in 60 Hz areas or 6 s in 50 Hz areas.

Notes: The technique of digital linear filtering is used for transformation of apparent resistivity data from one electrode configuration into another. Usually filter spectra are determined via the discrete Fourier transforms of input and output functions: the filter characteristic is the quotient of the spectra of the output function and input function.In this paper, the transformation of the apparent resistivities is presented for four electrode configurations (Wenner, the two-electrode, Schlumberger, and dipole configurations). In our method, there is no need to use the discrete Fourier transform of the input and output functions in order to determine the filter spectrum for converting apparent resistivity in one electrode configuration to any other configuration. Sine responses for determination of the derivative of apparent resistivities are given in analytical form.If the filter spectrum for converting the apparent resistivity to the resistivity transform for one electrode configuration is known, the filter spectra for transforming the apparent resistivity to the resistivity transform for any electrode configurations can be calculated by using newly derived expressions.

Notes: In the western coal-mining area of Ruhrkohle AG, reflection seismic prospecting for the Carboniferous coal measures is severely impaired by structures with halokinetic features. These structures make the interface between Mesozoic and Paleozoic layers, i.e., the top of Zechstein in general, very rugged. Unfortunately the velocity contrast at this interface is very high in that area, the ratio of velocities being 1.5 to 2.0. Therefore, migration and stacking become a problem.Three types of migration are presented:〈list xml:id="l1" style="custom"〉1(f, x)-time-migration with vertical time-to-depth conversion as a second step.2Kirchhoff migration down to a level determined approximately by the highest points of the top of Paleozoics, i.e., 0.35 s, and Kirchhoff-downward continuation for all times exceeding 0.35 s. Intermediate static corrections for these latter times with subsequent (f, k)-time-migration and final vertical time-to-depth conversion.3Direct depth migration in the (f, x)-domain using three interval velocities.In all cases an intermediate picking of the velocity interfaces is necessary. In case 2 this occurs at an earlier stage of the process than in case 1, and in case 3 at a still earlier stage.The results of the second and third migration procedures are superior to those of the first. Possibilities for misinterpretation of faults are reduced considerably when the second or third migration procedure is applied.

Notes: Airborne electromagnetic methods are most commonly used in mineral exploration. However, new developments, such as multifrequency capability and digital on-board field recording, as well as improvements in instrumentation resulting in high signal-to-noise ratios in recorded data, have made their application in geological mapping possible.A three-frequency airborne EM survey carried out over an area northwest of Timmins, Ontario, was interpreted in terms of thickness and resistivity of the layers of a two-layer earth section. Since both in-phase and quadrature components are measured, this provides six independent parameters at each point in space. Based on prior geological information and a preliminary interpretation of the field records, two two-layer models of the subsurface seemed to be appropriate for most of the survey area. An automatic computerized interpretation procedure was devised to interpret the field data at each point in terms of thickness and resistivity parameters of those two models. When the geology is more complex, the data do not fit the models and no interpretations are made. Two maps illustrating the variation of the resistivity and the thicknesses of the layers were constructed from the interpreted data. These maps agree with the known geological information about the distribution of glacial clay in the area. Areas where the layered models do not fit are known to be areas where the geology is complex with a large number of dykes and other lateral inhomogeneities. The study shows that multifrequency airborne EM surveys can be very useful in geological mapping over inaccessible terrain and can significantly help the mapping geologist where outcrops are scarce.

Notes: Shear (= S) wave studies in addition to compressional (= P) wave surveys have revealed that S-wave velocities are much more variable than P-wave velocities. This strongly affects Poisson's ratio σ, especially in young sediments. It is shown that σ has a great influence on the directivity pattern, i.e., on the radiation and receiving characteristic of horizontal sources and receivers. For their calculation, well-known surface boundary conditions and White's statement of the reciprocity relation (White 1965) are used. They also form the basis for the combination of source and receiver patterns essential for practical field work. The various combinations of horizontally and vertically arranged sources and receivers are investigated with regard to their directivity pattern, their amplitude ratio and, in general, to their usefulness in seismic prospecting.

Notes: Previously ignored characteristics of the seismic recording instrument are presently experienced as limitations as more sophisticated interpretive methods using wider frequency ranges are developed to extract stratigraphic information from seismic land data for hydrocarbon and mineral exploration. Most of these limitations arise from inadequate characteristics of the first element of the seismic instrument: the geophone. A geophone does not faithfully follow the motion of the earth for higher frequencies due to poor geophone-earth coupling. This filtering effect brings about time shifts that are dependent on the frequency and the soil type. A geophone can also produce spurious outputs, brought about by the motion of the suspended part of the geophone, with a magnitude comparable to that of the desired output. The suspension is made very compliant to obtain the required sensitivity. A compliant suspension, however, gives a large sag. The geophone can therefore only be used in one position, tolerating little tilt. A compliant suspension also widens the traveling range of the movable part. Minor sensitivity changes with travel are then noticeable as nonlinearity, since the surface wave is large with respect to the reflected wave. A compliant suspension is usually realized in the form of thin, spirally shaped spring-spiders. Such suspensions exhibit transverse or rotational resonances that are in or close to the seismic frequency band. Excited by ground roll, they can produce considerable undesirable output.The novel geophone we describe is a light-weight (17 g) acceleration-sensitive transducer which gives good ground coupling and partial correction for the increasing damping in the earth with increasing frequencies. It employs internal hybrid electronics for a magnetodynamic velocity-nulling feedback system. Velocity nulling makes the movable part of the geophone virtually rigid with respect to the housing. This makes the geophone characteristics independent of the suspension. The springs used are stiff in a transverse and rotational direction so that the suspension resonances are well outside the useful frequency band. This suspension also allows the geophone to be used in any orientation while being only sensitive to the vibration component along the main axis. The feedback system makes the sensitivity flat within 1 dB from 2 Hz to 500 Hz, with a phase tolerance smaller than 5°. The geophone is robust, has no moving internal wires, employs a current output [sensitivity 1 mA/(m s−2)] and internal gain so that the signal-to-cable-noise ratio is improved. This type of output allows parallel connection without any interaction between the geophones.

Notes: There are two types of masked layers in seismic refraction work: the velocity reversal (low-velocity layer) and the hidden layer (insufficient velocity contrast or layer thickness). On the basis of an analytical formulation of the general case of a masked layer under an overburden of plane and parallel multiple refractors the two limiting cases are discussed: the solution resulting from an uncritical interpretation of the measured time-distance curve and the blind zone solution. Between these two limiting cases there is a variety of possible masked layer solutions. These no-blind zone solutions—as well as the blind zone solution itself—are formulated separately for the velocity inversion and the hidden layer case.For the evaluation of some no-blind zone solution a diagram is presented which can be used for any case of multiple refractors in the overburden of the masked layer. However, it is only for the three- and the four-layer case that a blind zone interpretation by use of diagrams is advisable. Such diagrams are presented together with the basic sets of formulae which contain as parameters only ratios of velocities and layer thicknesses. As the velocity of the masked layer is usually unknown the diagrams are principally constructed to show the dependence on the masked layer velocity. This is useful for estimation of the largest possible error.

Notes: Wave equation migration techniques have shown the limits of traditional stacking methods with data from tectonically complicated areas. An improved stack can be obtained utilizing the dip-moveout correction technique based on offset continuation. The properties and the limits of the algorithms used are summarized briefly.Several synthetic and real data examples are shown and compared with the results obtained using conventional processing in order to show the focusing effects and the strong improvement in signal-to-noise ratios, both at the stacked and migrated section level. The possibility of exploiting this technique to transform multiple coverage into increased spatial resolution is illustrated with examples.

Notes: The far-field signatures from a comprehensive and systematic airgun pulse test have been analyzed. Empirical relations between the characteristic signature parameters and depth (5–12 m), pressure (100–137 bar = 10–13.7 MPa) and total chamber volume (0.65–9.5 l) have been derived. Also, the influence of using waveshape kits in different positions within the chamber has been tested.The results indicate that:〈list xml:id="l1" style="custom"〉1 The amplitude is proportional to chamber pressure to the power 3/4.2 The bubble period is nearly independent of the position of the waveshape plate.3 The increase in primary/bubble amplitude ratio is inversely proportional to the chamber volume above the waveshape plate.4 The amplitude is independent of airgun depth.Suggestions and comments regarding this work from Dr B. Ursin and Dr A. Ziolkowski are appreciated. The field work was supported by the Norwegian Petroleum Directorate through the Continental Shelf Project at the Seismological Observatory, University of Bergen.An airgun allowing for continuous variation of the chamber volumes was supplied by GECO (Geophysical Company of Norway). The purchase of two airguns was financed by Norske Getty Exploration A/S.

Notes: The accuracy of short length digital linear filter operators can be substantially increased if the sampling interval as well as the abscissa shift are properly adjusted. This may be done by a trial and error process of adjustment of these parameters until the error made by the filter operator, applied to a suitably chosen test function, is smallest.As an illustration of the application of this method, 7-, 11- and 19-point filters for the calculation of Schlumberger apparent resistivity from a known resistivity transform are designed. Errors with the new 7-point filter are seen to be less than those with a 19-point filter of conventional design. The errors with the new 19-point filter are two to three orders of magnitude smaller than those made by the conventional 19-point filter.The new method should provide digital linear operators that allow significant improvements in accuracy for comparable computation efforts, or substantial reduction in computation for comparable accuracy of results, or something of both.

Notes: A well-known technique for the migration of normal-incidence two-way travel-time maps is extended to common-source-point travel-time data. The travel time and the travel-time gradient are used to compute the parameters defining the tangent plane of the reflecting interface. It is also shown how the curvature matrix of the received wavefront can be used to compute the curvature of the reflecting interface. The method is initially derived for common-source-point data and then extended to common-midpoint data.In a three-dimensional medium the wavefront curvature matrix is computed by solving a 2 × 2 symmetric matrix Riccati equation. In a two-dimensional medium and in a medium with constant velocity gradient, the wavefront curvature matrix is computed by solving a scalar Riccati equation and two linear equations. The migration procedures are also simplified.When the velocity function is unknown, the migration procedures cannot be used. An inverse modeling algorithm which simultaneously performs the migration and estimates the velocity function must then be applied. Two different inversion schemes are discussed briefly.

Notes: Time-domain-induced polarization (IP) laboratory measurements were performed on about 200 fine sediment samples with varying water content. The results permitted an analysis of IP properties of clays, loams, silts, and sands.Particular emphasis has been given to the analysis of the chargeability m as a function of lithotype and the water content.By analyzing decay curves, a new parameter was identified. It is a statistically specific characteristic of the lithotype and is independent of the water content. Therefore, it provides a diagnostic parameter for lithotype identification. In association with the values of chargeability and electrical resistivity, this parameter permits a reliable evaluation of water content and yields useful information about the porosity and permeability of the lithotype.

Notes: Approaches to the reduction of bias in the computation of the elements of the magnetotelluric impedance tensor have been proposed in the past by several authors. In this paper a clear distinction is made between random errors and bias errors. No effort is made to reduce either, but the emphasis is on their estimation. Both types of errors depend critically upon the polarization of the magnetic field. The random error increases with increasing noise-to-signal ratio in the electrical field, and it is rather insensitive to noise in the magnetic field. The bias error increases with increasing noise-to-signal ratio in the magnetic field. Expressions for random errors and maximum bias errors are developed and discussed using a single station set-up. Random errors with a reference station set-up are also calculated.

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: Dipole-dipole induced polarization (IP) data are displayed typically as multi-level profiles, or as contours on vertical sectional plots referred to as pseudo-sections. The dipole-dipole array tends to yield IP anomalies in which the most anomalous values are displaced laterally from the source body. The data patterns are fairly interpretable on pseudo-sections or on multi-level profiles but are sufficiently complex to discourage the contouring of the data in plan.A method was developed for the presentation of dipole-dipole IP data on a contour map. The method consists of a simple averaging of data which can be performed manually if desired. It yields a single output value per station which reflects all levels of the pseudo-section, and is suitable for contouring in plan. The advantage of the technique is that it provides a quantitative picture of IP anomalies in their background or regional setting.

Notes: Constant offset sections can be mapped to a fixed offset and compared in order to provide a method of velocity analysis. The direct mapping to zero offset prior to stack might provide an alternative processing procedure to NMO and stack. The main advantage of such a procedure would lie in the correct treatment of cross-dips, but interpretational advantages might also follow from the performance of partial stacks biased to either high or low offset information.

Notes: Induced polarization and resistivity model studies over thin dykes for varying resistivity contrasts, depth of burial, and dip angles show striking parallelism of the surface apparent resistivity contours with the boundary of the body. This effect may be utilized for the estimation of strike length of the body. Results show that intermediate values of the electrode spacing is satisfactory for detection under widely varying conditions of resistivity contrast and depth. The percentage frequency effect (P.F.E.) and the metal factor (M.F.) responses are found to be more sensitive to the variation in the depth of burial than the resistivity responses. Pseudosections for P.F.E. and M.F. are concentrated much closer to the body than the resistivity pseudo-sections.

Notes: A direct interpretation scheme is developed which is capable of determining most of the geological features of a ground which can be assumed to be two dimensional in structure. This scheme extends the earlier work of Pekeris (1940) and Koefoed (1968) to the case where the basal layer of a ground is undulating. It also has a limited use for finding the parameters of a dipping dyke in the lower medium. Though the top and dip of the dyke can be determined, this is not true for the thickness.

Notes: Laboratory studies of absorption-frequency behavior in rocks often use spectral ratios of digitally recorded ultrasonic signals which have been transmitted through a rock sample and a reference sample of very low absorption, respectively. It is proposed to treat the digitally recorded signals as an autoregressive-moving average (ARMA) process which, using recursive filter concepts, can be represented as a ratio of two polynomials in the z-transform variable z. The numerator polynomial contains only that part of the signal that is modified by anelastic effects, whereas the denominator contains the elastic effects of the physical apparatus such as reverberations. Examples are given which show that this separation of the recorded signal greatly facilitates the laboratory investigation of loss mechanisms and absorption-frequency behavior based on spectral ratios.

Notes: The generalized integral equation for the electric potential governed by a quasi-harmonic equation can be derived via a variational formulation. For surface current distributions it is not always a Fredholm integral equation of the second kind. Numerical solutions of the general heterogeneous problem can be obtained with the “reciprocal averaging technique”, where the solution is obtained a second time after exchange of source and field points.

Notes: A method to calculate the resistivity transform of Schlumberger VES curves has been developed. It consists in approximating the field apparent resistivity data by utilizing a linear combination of simple functions, which must satisfy the following requirements: (i) they must be suitable for fitting the resistivity data; (ii) once the fitting function has been obtained they allow the kernel to be determined in an analytic way.The fitting operation is carried out by the least mean squares method, which also accomplishes a useful smoothing of the field curve (and therefore a partial noise filtering). It gives the possibility of assigning different weights to the apparent resistivity values to be approximated according to their different reliability.For several examples (theoretical resistivity curves in order to estimate the precision of the method and with field data to verify the practicality) yield good results with short execution time independent of shape the apparent resistivity curve.

Notes: We present a new method for the extraction and removal of the source wavelet from the reflection seismogram. In contrast to all other methods currently in use, this one does not demand that there be any mathematically convenient relationship between the phase spectrum of the source wavelet and the phase spectrum of the earth impulse response. Instead, it requires a fundamental change in the field technique such that two different seismograms are now generated from each source-receiver pair: the source and receiver locations stay the same, but the source used to generate one seismogram is a scaled version of the source used to generate the other. A scaling law provides the relationship between the two source signatures and permits the earth impulse response to be extracted from the seismograms without any of the usual assumptions about phase.We derive the scaling law for point sources in an homogeneous isotropic medium. Next, we describe a method for the solution of the set of three simultaneous equations and test it rigorously using a variety of synthetic data and two types of synthetic source waveform: damped sine waves and non-minimum-phase air gun waveforms. Finally we demonstrate that this method is stable in the presence of noise.

Notes: In this paper a theorem is demonstrated which allows—after the introduction of a suitable dipole kernel function or dipole resistivity transform function—to write the apparent resistivity function as an Hankel transformable integral expression.As a practical application of the theorem a procedure of quantitative interpretation of dipole soundings is suggested in which the dipole resistivity transform function obtained after inversion of the original dipole apparent resistivity data is used to control the goodness of the set of layering parameters which have been derived with our previous method of transformation of dipole sounding curves into equivalent Schlumberger diagrams.

Notes: The Hankel transform theorem can be applied to the inversion of gravity data for the buried sphere, the horizontal cylinder, and the vertical rod. This new approach leads to exact solutions of the transforms for the assumed bodies. A comparison with the classical procedure by Fourier transform reveals that for the sphere and the vertical rod, the Hankel transform is preferable.

Notes: Air guns have been used in various applications for a number of years. They were first used in coal-mining operations and were operated at up to 16000 psi charge pressures. Later, single air guns, operated at 2000 psi, found application as an oceanographic survey tool. Air gun arrays were first used in offshore seismic exploration in the mid-1960's. These early arrays were several hundred cubic inches in total volume and were operated at 2000 psi; they were either tuned arrays or several large guns of the same size with wave-shape kits. Today's arrays have total volumes greater than 5000 cu in. and are typically operated at 2000 psi. Recently, higher-pressure, lower-volume arrays operated at 4000–5000 psi have been introduced; guns used in these arrays are descendants of the coal-mining gun.On first thought one would equate increased gun pressure linearly with the amplitude of the initial pulse. This is approximately true for the signature radiated by a “free-bubble” (no confining vessel) and recorded broadband. The exact relation depends on the depth at which the gun is operated; from solution of the free-bubble oscillation equation, the relation is 〈displayedItem type="mathematics" xml:id="mu1" numbered="no"〉〈mediaResource alt="image" href="urn:x-wiley:00168025:GPR700:GPR_700_mu1"/〉 If Pc,1= 6014.7 psia, Pc,2= 2014.7 psia and PO, 1=PO, 2= 25.8 psia (corresponding to absolute pressure at 25 ft water depth), then 〈displayedItem type="mathematics" xml:id="mu2" numbered="no"〉〈mediaResource alt="image" href="urn:x-wiley:00168025:GPR700:GPR_700_mu2"/〉 Experiments were conducted offshore California in deep water to determine the performance of several models of air guns at pressures ranging from 2000 to 6000 psi and gun volumes ranging from 5 to 300 cu in. At a given gun pressure, the initial acoustic pulse Pa correlated with gun volume Vc according to the classical relation 〈displayedItem type="mathematics" xml:id="mu3" numbered="no"〉〈mediaResource alt="image" href="urn:x-wiley:00168025:GPR700:GPR_700_mu3"/〉 For 1 ms sampled data the ratio 〈displayedItem type="mathematics" xml:id="mu4" numbered="no"〉〈mediaResource alt="image" href="urn:x-wiley:00168025:GPR700:GPR_700_mu4"/〉 varied between 4.5 and 5.5 dB depending on gun model. Pulse width of the 2000 psi signatures indicated they are compatible with 2 ms sample-rate recording while pulse width of the 6000 psi signatures was greater, indicating they are less compatible with 2 ms sample-rate recording.Conclusions reached were that 2000 psi air guns are more efficient than higher pressure guns and are more compatible with 2 ms sample-rate requirements.

Notes: Numerous electrodes, already used in geophysics or just perfected by us, have been compared by measuring the three main characteristics which interest the user: noise spectrum, temperature coefficient and polarization with its stability versus time.Among the most used unpolarizable electrodes, silver-silver chloride (Ag-AgCl) are the best ones. But a systematic research of all different possible metal-salt couples, have led us to use lead-lead chloride (Pb-PbCl2) for the following reasons: noise as low as the one of Ag-AgCl at 1 Hz and even lower for the low frequencies (0.4 μV at 1 Hz and 1.2 μV at 0.01 Hz for peak to peak value and ΔF=F), temperature coefficient about ten times weaker (−40 μV/°C instead of −450 μV/°C) and also better long time stability of the polarization (1 mV/month instead of 2 at 10 mV/month).We have been using these electrodes since 1977 as “tube” electrodes which are very easy to use. They allow us to record correctly the fast variations thanks to their low noise, the very slow variations, their low temperature coefficient and their stability, and this with telluric lines only about 100 m long.

Notes: The compression of seismic signals which have propagated through a dispersive medium can be achieved by a frequency-domain transformation. This transformation is formally related to the dispersion characteristic of the medium and is independent of phase and group velocities. By defining a suitable cost-function to measure the degree of dispersion of a time- signal, an iterative technique can be employed to find that transformation which gives minimum dispersion.In this paper, we assume that the inverse of the dispersion characteristic can be adequately approximated by a finite polynomial in the region of maximum signal energy density. The coefficients of this polynomial are the parameters of dispersion of the medium. These parameters can be estimated both in the presence of noise and in the case of signals made up of multiple arrivals.The techniques developed in this paper are applied to seismic signals which have been recorded underground by a set of geophones designed to detect and map discontinuities in coal-seams. Results of dispersion estimation and signal compression are discussed for data collected in the Blackshale seam, Pye Hill Colliery, near Derby, in the United Kingdom.

Notes: Sign reversals have been observed with the one-loop version of the TEM method in an area near Cloncurry, Queensland. This is the only area in which such a response has been recorded in seven years of field surveys by the Australian Bureau of Mineral Resources.The geology of the area consists mainly of pyrrhotitic graphitic shale, which has resistivities as low as 0.1 Ωm and frequency effects up to 30%. The sign reversals could be due to magnetic phenomena, reflections from layers, or complex conductivity effects.