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: 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: The widespread use of common depth point techniques has emphasized the need for accurate static corrections. Manual interpretation methods can give excellent results, but a computer technique is desirable because of the great volumn of data recorded in common depth point shooting. The redundancy inherent in common depth point data may be used to compute a statistical estimate of the static corrections. The corrections are assumed to be time-invarient, surface-consistent, and independent of frequency. Surface consistency implies that all traces from a particular shot will receive the same shot static correction and all traces from a particular receiver position will receive the same receiver correction.Time shifts are computed for all input traces using crosscorrelation functions between common depth point traces. The time shift for each trace is composed of a shot static, a receiver static, residual normal moveout if present, and noise. Estimates of the shot and receiver static corrections are obtained by averaging different sets of the measured time shifts. Time shifts which are greatly in error are detected and removed from the computations.The method is useful for data which has a moderate to good signal to noise ratio. Residual normal moveout should be corrected before estimating the statics. The program estimates the statics for correctly stacking common depth point traces but it is not sensitive to constant or very slowly changing static errors.

Notes: The principle of the mise-á-la-masse method is to earth one current electrode of a pair in a conducting mineral show (in a borehole, in an outcrop etc.) and measure the resulting distribution of electric potentials. The distribution will, to some extent, reflect the geometry of the ore mass of which the mineral show forms a part.In a mise-á-la-masse survey of a lead-zinc deposit in Central Sweden electric potentials were measured on the surface of the ground as well as in some 25 drillholes, in either case with earthings (successively) in three different parts of the irregular ore deposit. Besides this, measurements were made in drillholes with earthings in two further drillholes.Geologic correlation between the drillholes is difficult in the present case on account of the irregular geometry of the ore deposit. However, the mise-á-la-masse measurements clearly showed the dip and the pitch of the ore body, established connections between the different ore widths encountered in the various holes, and yielded information about the shape of the ore mass.In particular, the survey showed that the ore lenses must be crescent-shaped rather than tabular, and the dip was indicated to be westerly, instead of easterly as originally presumed.Three-dimensional models of equipotential surfaces were constructed from the observed drillhole and surface potentials (using transparent plastic sheets and thin copper wire) and these helped to elucidate the mass geometry further.The surface and underground potential data collected in the present case should be of nterest to geophysicists working on analytic continuation problems.

Notes: The author comments the techniques and methods used in continuous seismic profiling by the Monaco Oceanographic Museum during the last five years.Pinger probe give informations on the Mud layering and thickness. On the shelf or in abyssal plains, boomer and sparker are employed and allow to rise various geological conclusions. South coast of France, continental shelf in the Gulf of Lion, abyssal plain of the Ligurian Sea and Lake of Geneva are the areas concerned to with the results.

Notes: The author proposes the parameter Δ2/ΔT1 for possible application in aeromagnetic surveying making use of already available gradiometer systems equipped with sensors separated vertically.

Notes: Geophysical field equipment has undergone rapid changes in the past decade; from simple AGC amplifiers and galvo cameras to binary gain switching amplifiers and digital recorders, all in an attempt to keep pace with the new geophysical interpretive methods developed, and the growing acceptance of the terminology, methods, and philosophy of communication theory.The additional tools of the digital recorder and digital computer make it possible to utilize these new techniques in geophysical processing. Accomplishing these new techniques demands severe requirements on the digital field recording process in handling the decreasing energy return from the seismometer, and to fully realize the capabilities of digital techniques in reducing data. Simple automatic gain control may be used. However, in the more sophisticated interpretive methods, such as autoregression and deconvolution, it is necessary to reconstruct the actual energy levels in the computer. Recording the control signal used in master AGC or programmed gain control may prove satisfactory; however the accuracy of control versus gain is limited to I% for such analog instrumentation.To utilize the computer to its fullest extent, and to accurately perform these new techniques, requires an accuracy of I% or better. This accuracy is obtainable by using a step gain control where the gain is increased by fixed steps in which each step represents a gain in amplification by a constant factor. The accuracy in this case can be made dependent only on the tolerance of resistors used as attenuators or feedback elements. Preferably the constant factor of gain steps should be a number easily handled by the computer. By using 6 db steps it becomes a simple matter to shift binary numbers, such as multiplying or dividing by 10 in the decimal system.The requirements or parameters for such an amplifier system, and the elements of the amplifier necessary to achieve these requirements are presented.

Notes: Spurious reflections showing a residual move-out are generally present on both corrected seismic cross-sections and common depth point composites.We propose to determine a space-time filter satisfying the following conditions: to attenuate as much as possible reflections presenting a given move-out, to retain intact reflections whose time gradient is zero (or has a predetermined value), to be applicable efficiently even to a small number of traces, not to amplify random noise unduly.After briefly indicating the design principle of these filters, we shall give the results of their application to the theoretical examples, in order to bring out: the influence of sampling interval and filter length (possibility of achieving efficient filtering using an apparatus with a small number of terms), the effect of filtering on dipping reflections (reduction in amplitude and distortion increasing with the time gradient), the range of reflection move-out that can effectively be filtered with the same apparatus, the disturbing influence of random noise.〈section xml:id="abs1-2"〉〈title type="main"〉RESUMEOn constate généralement, tant sur les films sismiques corrigés que sur les résultats de couverture multiple, la présence de réflexions parasites présentant une courbure résiduelle d'indicatrice.Nous nous sommes proposé de déterminer un opérateur de filtrage spatio-temporel satisfaisant aux conditions suivantes: filtrer le plus possible les réflexions présentant une courbure donnée, conserver intégralement les réflexions dont le gradient temps est nul (ou a une valeur déterminée), pourvoir être appliqué efficacement à un petit nombre de traces, ne pas amplifier outre mesure les bruits inorganicés.Après avoir indiqué brièvement le principe du calcul de ces filtres, nous présenterons les résultats de leur application à des exemples théoriques, afin de mettre en évidence:ľ'influence du pas d'échantillonnage et de la longueur du filtre (possibilityé d'obtenir un filtrage efficace avec un opérateur comportant un petit nombre de termes),ľ'effet du filtrage sur les réflexions pentées (réduction d'amplitude et distorsion croissant avec le gradient temps),ľ'ouverture de ľ'éventail des réflexions courbes que ľ'on peut filtrer efficacement avec un même opérateur,ľ'influence perturbatrice de bruits inorganicés.

Notes: An iterative process is proposed that computes a distribution of masses giving rise to a certain gravitational field. The possibility of applying a similar procedure to magnetic and electric fields is also discussed.

Notes: The method presented takes the influence of refraction rigidly into account. All travel-times available are used in an application of the method of least squares. The solution is illustrated by a computer flow diagram.

Notes: Charts are presented for the computation of the vertical component of the attraction of horizontal laminas and of horizontal cylinders bounded by irregular w-sided polygons. Contrary to older methods no restricting conditions are made for the directions of the sides of the polygons.

Notes: It was found in Part I of this paper that approximating the sharp cut-off frequency characteristic best in a mean square sense by an impulse response of finite length M produced a characteristic whose slope on a linear frequency scale was proportional to the length of impulse response, but whose maximum overshoot of ±9% was independent of this length (Gibbs' phenomenon). Weighting functions, based on frequency tapering or arbitrarily chosen, were used in Part II to modify the truncated impulse response of the sharp cut-off frequency characteristic, and thereby obtain a trade-off between the value of maximum overshoot and the sharpness of the resulting characteristic. These weighting functions, known as apodising functions, were dependent on the time-bandwidth product Mξ, where 2ξ, corresponded to the tapering range of frequencies.Part III now deals with digital filters where the number 2N–1 of coefficients is directly related to the finite length M of the continuous impulse response. The values of the filter coefficients are taken from the continuous impulse response at the sampling instants, and the resulting characteristic is approximately the same as that derived in Part II for the continuous finite length impulse response. Corresponding to known types of frequency tapering, we now specify a filter characteristic which is undefined in the tapering range, and determine the filter coefficients according to a mean square criterion over the rest of the frequency spectrum. The resulting characteristic is dependent on the time bandwidth product Mξ= (N–1/2)ξ up to a maximum value of 2, beyond which undesirable effects occur. This optimum partially specified characteristic is an improvement on the previous digital filters in terms of the trade-off ratio for values of maximum overshoot less than 1%. Similar to the previous optimum characteristic is the optimum partially specified weighted digital filter, where greater “emphasis is placed on reducing the value of maximum overshoot than of maximum undershoot”. Such characteristics are capable of providing better trade-off ratios than the other filters for maximum overshoots greater than 1/2%. However these filters have critical maximum numbers 2.NC–1 of coefficients, beyond which the resulting characteristics have unsuitable shapes. This type of characteristic differs from the others in not being a biassed odd function about its cut-off frequency.

Notes: In Part I of this paper, we examined the properties of the best mean square approximation to the sharp cut-off frequency characteristic by an impulse response of finite length. It was found that the sharpness of cut-off for the resulting frequency characteristic depended on the length of the impulse response–but because of the discontinuous nature of the specified frequency characteristic, this best mean square approximation always had a maximum overshoot of ± 9%, independent of the length of the impulse response (Gibbs phenomenon).In Part II, we investigate ways of reducing this ± 9% overshoot at the expense of a reduced sharpness of cut-off. The discontinuous frequency characteristic is first approximated by a continuous characteristic with linear or cosine frequency tapering. The impulse response for such tapered characteristics consists of the impulse response of the discontinuous frequency characteristic weighted by a certain function corresponding to the type of tapering employed. The best mean square approximation to the tapered characteristic by an impulse response of finite length M will produce a frequency characteristic whose properties are now dependent on the time-band width product Mζ, where 2ζ is the tapering range.A trade-off exists between the maximum overshoot and the sharpness of cut-off for the resulting characteristic for both forms of frequency tapering. Instead of considering other forms of tapering in the frequency domain, we now investigate arbitrarily chosen weighting functions in the time domain to determine the minimum length of impulse response for a minimum value of maximum overshoot and a maximum value of sharpness of cut-off.Part III will discuss the digital realization of the above finite length impulse responses together with the optimum partially specified digital filter approximation to the desired frequency characteristic.

Notes: Book reviewed in this article:Le Filtrage en Sismique, Tome IJ. AUBOUIN “Geosynclines” (Developments in Geotectonics 1). Elsevier Publishing CompanyEarth and Planetary Science Letters Vol. 1. Nr. 2H. Ramberg Gravity, Deformation and the Earth's Crust Academic Press“Potassium Argon Dating” Compiled by O. A. Schaeffer and J. Zähringer

Notes: Deep seismic sounding was performed along two profiles which cross at the Dinarides area right angles. One of the profiles goes far into the Adriatic Sea.Besides considerations on the lithophysical conditions, characteristics of the registered waves are analysed. The amplitude curves and curves of amplitude ratios are shown. Special attention was paid to the frequency of the registered waves.In order to obtain a better knowledge of the registered wave pattern three-component recordings of waves were carried out. The analysis of the records obtained is given, with particular regard to the possibility of creating converted waves.The Earth's crust structure along the profiles II and III is given.

Notes: Nomograms for solving equations in multilayer and dipping layer cases are presented. The nomograms constructed are used to solve the following equations: I. Intercept-time formula. 2. Critical distance formula. 3. Critical angle formula. 4. Critical angle and dip angle formula. 5. Vertical depth formula.

Notes: Consideration is given to the use of a configuration of four electrodes set in a square array for resistivity measurements.It is found that, by passing current successively between different pairs of electrodes, an apparent resistivity can be determined which is both more sensitive to the position of the array centre and less dependent on orientation than the measures usually obtained with colinear arrays of electrodes. At the same time the observations made enable the degree of the departure of local conditions from conditions of lateral homogeneity to be assessed.Theoretical and practical examples of the use of this electrode system are given and the use of the system both as a tool in mapping and in depth investigations is considered.It is shown that provided electrode spacings are suitably arranged the results of a probe carried out using the square array can be interpreted by conventional methods.The system is shown to have particular advantages in the investigation of lateral resistivity variations and the reduced dependence on orientation makes possible the recasting of interpretation data in an orientationally invariant form with a consequent drastic reduction in the number of type curves required for a particular problem.

Notes: A large self-potential anomaly was outlined in 1963–64 at Tantahuatay near Hualgayoc in the Andes of northern Peru. Peak value recorded was–1842 millivolts–thus making it one of the strongest, or perhaps the strongest, SP anomaly ever measured. A lack of detailed geological data precludes the formulation of an adequate explanation for the Tantahuatay anomaly, but geological and mineralogical similarities with the well-documented Venencocha anomaly near Cerro de Pasco, Peru (Kruger and Lacy, 1949), suggest that the anomaly arises from sulfuric acid associated with the mineral alunite. The anomaly obviously cannot be explained by the half-cell mechanism of Sato and Mooney (1960), who place a limit of 700 millivolts on self-potential anomalies over sulfide bodies. Further study of the Tantahuatay anomaly would be of interest in understanding self-potential mechanisms in general.

Notes: Book review in this articleF. S. Grant and G. F. West, Interpretation Theory in Applied Geophysics, McGraw-Hill Comp., New York, 1965, 583 pp., price $ 17.50.Archaeo-Physica. Technische und naturwissenschaftliche Beiträge zur Feldarchäologie. Mit Beitragen von E. Le Borgne, I. Scollar, J. D. Mudie, J. Görier, H.-J. Appel, G. Strunk-Lichtenberg, Böhlau Verlag, Köln, Graz, 1965, 202 pp., 124 fig.H. Takeuchi, Theory of the Earth's Interior, Blaisdell Publishing Company, Waltham (Mass.), Toronto, London, 1966, 131 pp.G. D. Garland, The Earth's Shape and Gravity, The Commonwealth and International Library, Geophysics Division, Pergamon Press, Oxford, London, Edinburgh, New York, Paris, Frankfurt, 1965, 183 pp., price: 17s 6d.

Notes: According to a study of seismic velocities in the Alpine Foreland of Eastern Switzerland, the initial velocity is rather high in comparison with other Tertiary basins and shows an exceptionally high increase rate. When analysing the average and the interval velocities, especially of Tertiary strata, and when comparing them with velocities of strata of the same age and a very similar facies of the Northern Rheintalgraben, it has been found that the increase of velocity is closely related to the distance to the Subalpine Molasse. The conclusion is that the velocity of the Tertiary strata is strongly influenced by the folding pressure of the Subalpine Molasse or of the Alps. The same method has been applied to a relatively large number of wells in the area of the “German Molasse”. Not only the results in Eastern Switzerland have been confirmed, but also it has been proved that the diagenesis of the Tertiary strata and, hence, their velocities are influenced only partially by the specific depth of the basin. Velocities increase towards the Folded Molasse or the Alps. Consequently they depend on lateral folding pressure, which decreases from west to east with the increasing width of the basin. The tertiary strata of the basin have been affected by lateral folding pressurefrom south to north. However, structures with lateral compression have not been discovered yet in the German Alpine Foreland.

Notes: In a paper by Koefoed and Kegge (1968), which was based on previous work of Wesley (1958), the electrical current pattern has been derived that is induced by an oscillating magnetic dipole in a semi-infinite thin plate of infinitesimal resistivity.In the present paper, the range of validity of the assumptions, on which the work of Wesley is based, is subjected to a theoretical analysis. It is found that the decisive factor for the validity of Wesley's derivation is the quotient of the square of the penetration depth of the electrical current over the product of the thickness of the plate and a distance that is indicative of the size of the current loops in the plate. Wesley's derivation is shown to be valid only when this factor is negligible. It is also shown that in this condition the imaginary component of the anomaly must be negligible.Model experiments are described in which the electrical current pattern is studied also in the range in which the derivation of Wesley is not valid. The procedure used in these model experiments was to measure the tangential component of the magnetic field strength very close to the metal plate that simulated the conductive dyke.In order to express the results of the measurements in terms of the imaginary to real ratio, these results are compared with an interpretation graph for field measurements that was published by Hedström and Parasnis (1958). It is found that the current pattern in the plate is essentially the same as that which follows from Wesley's derivation, provided that the imaginary to real ratio is less than one third. The measurements do not permit to draw conclusions regarding the current pattern in the plate in conditions that correspond to larger values of the imaginary to real ratio.

Notes: Electrical resistivity anomalies of a symmetrical four-electrode co-linear system across a single vertical discontinuity are treated in relation to: a) the ratio of potential electrode separation to the current electrode separation that are employed in the system and b) the angle which the electrode alignment makes with the discontinuity. Several conclusions are extracted from this treatment and methods for obtaining an optimum sensitivity of the system, with respect to these parameters, are shown. Disadvantages of special arrangements, such as the Wenner configuration, are indicated. Methods are outlined to utilise variations in the apparent resistivity plot for determining the angle between the electrode alignment and the discontinuity, quantitatively or qualitatively. These variations include certain deviations from the standard curves obtained in longitudinal traverses made at right angles to the discontinuity. Also, a comparison is made between longitudinal and cross traverses, in relation to the discontinuity.

Notes: When, in a two-layer earth, the substratum is insulating or infinitely conducting, the thickness of the top layer can be determined from surface potential measurements along a radius from a single point power electrode. The observed potential needs to be numerically integrated in a direction perpendicular to the said radius, and Figure 4 can then be used to find the thickness. A field example is included.

Notes: In order to increase the seismic efficiency of the Sparker-system developments were made in 1966/67. Results illustrating the main steps of this development are discussed. A new type of electrode called Transploder electrode was field-tested and proved to be promising.

Notes: The theory of electrical dipole soundings proved that this method can produce resistivity measurements, which are comparable with those obtained by electrical soundings of the Wenner or Schlumberger type. Their main advantage is the use of short cable lengths, which is important if the depth of penetration should be large.A considerable disadvantage of the dipole method is the great sensitivity to lateral discontinuities. Though these have an influence on the Schlumberger arrangement as well, they can disturb a dipole sounding to such an extent than an interpretation based on a horizontal layer case is no more possible.There are six different dipole arrays, which differ from each other with respect to the angle enclosed by the two dipole orientations-the current dipole AB and the measuring dipole MN. The theoretical comparison of the dipole arrays with the Schlumberger array concerning their sensitivity to lateral discontinuities is a useful basis for the choice of the most suitable configuration.Considering geological subsurface conditions the right choice of a dipole array can give an optimal result, i.e. a dipole sounding for which the sensitivity to lateral discontinuities is as small as possible under the given circumstances.

Notes: Deconvolution and deringing are well known subjects and it is not necessary to state again their objectives nor the basical methods used to reach them. Let us just remember that, generally, among many others, the two following assumptions are made for simplification purposes:〈list xml:id="l1" style="plain"〉—for deconvolution, it is assumed that the recorded seismic signal is constant, meaning that its shape is the same all along the time interval during which the trace is to be deconvolved;—for de-ringing, it is assumed that the ringing period is constant and that the intensity of the ringing phenomenon is independant of the time.With these two assumptions, a single constant operator can be applied for deconvolving, deringing or both. In most cases, the time variations of the signal or of the ringing are small enough and the error resulting of the application of a constant operator is acceptable. It results into a slight increase of the noise level or into a small residual ringing in the processed trace.When this noise or the residual ringing are too important, the assumption of a constant signal and ringing period must be rejected. This is the case that is examined here according to the following steps:〈list xml:id="l2" style="plain"〉—short definition of the problem;—fast evaluation of some possible solutions;—the selected solution: resulting approximations and how to obviate them, computing method and a remark about the operators;—theoretical example: the efficiency of the process used is evaluated on data in which the results aimed at are known; the influence of the selection of numerical values to be assigned to the parameters is examined;—real cases: comparison of results obtained with the Protee process and with more conventional processes assuming a time invariance or including a weighted composition of several conventional processes each with a different operator.

Notes: Sharp cut-off frequency filtering is carried out in the discrete time domain on digital computers. A convolution of the digital filter impulse response with the sampled input yields the output. For practical reasons, the length of the filter inpulse response, corresponding to the number of filter coefficients, is limited, and consequently the resulting frequency characteristic will no longer be identical to that originally specified. This is analogous to synthesising some specified frequency characteristic with a finite number of resistive, capacitative and inductive components.In Part I of this paper, we examine the effect of approximating the sharp cut-off frequency characteristic best in a mean square sense by an impulse response of finite length. The resulting frequency characteristic corresponds to the truncated impulse response of the specified frequency characteristic. It has a cut-off slope proportional to, and a mean square error inversely proportional to, the length of the impulse response, and is a biassed odd function about the cut-off frequency point. Because of the Gibbs phenomenon for discontinuous functions, the resulting frequency characteristic will always have a maximum overshoot with respect to the specified characteristic of ± 9%, regardless of the length of the corresponding impulse response. Equal length truncated impulse responses of specified filters with different cut-off frequencies yield frequency characteristics which are almost identical about their respective cut-off points. Now on a log frequency scale (as against a linear frequency scale implied previously) such characteristics may be made almost identical about the respective cut-off points by having the truncated impulse responses composed of an equal number of zero crossings. Results for the low-pass filter are applicable to the high-pass and band-pass characteristics.In the latter case, the mean square error is double that for a single slope characteristic (low-pass or high-pass) and the slopes at both edges of the passband are approximately equal in magnitude to the length of the impulse response (linear frequency scale).Part II of this paper is concerned with reducing the ± 9% overshoot that results from the discontinuous nature of the sharp cut-off frequency characteristic and which is not dependent on the length of the truncated impulse response. The reduction is achieved, at the expense of the steepness of cut-off for the resulting frequency characteristic, by the use of functions which weight the truncated impulse response of the specified frequency characteristic. These functions are called apodising functions. Among other variables, the length of the truncated weighted impulse response will determine the amount of maximum overshoot since the effective frequency characteristic being approximated is no longer a discontinuous function. The digital realization of the finite length impulse responses of Parts I and II is discussed in Part III, together with the optimum partially specified digital filter approximation to the desired frequency characteristic.

Notes: The purpose of this report is to show a method of determining the top of a refractor departing from the times and slopes of the direct and inverse dromocrones. The method does not need topographical correction and can be applied without knowledge of the distance between the geophone and the shot point.These results having been obtained, the commonly accepted point of view is upset: instead of looking for two points on the surface corresponding to one point of the refractor, we try to etablish, starting with only one point from the surface, the two corresponding points from the top of the refractor.This method can be applied to isolated points and does not demand interpretative hypotheses of any kind, excluding the velocity evaluation of the overburden and of the refractor.The necessary calculations can be easily executed by means of a digital computer to which the dromocrone times and the distances between the geophones must be given. These calculations can also be executed by a person having no knowledge of refraction seismology.This report also examines the validity of the approximations involved in the method proposed.

Notes: Determining the shape of a pulse generated by an explosion solely from the data provided by the recorded seismic trace is a difficult and even ambitious task.Knowledge of parameters such as length and number of “arches” of the pulse under study is, in fact, indispensable in solving this problem.These parameters cannot be found directly in the seismic trace, which nevertheless contains a great amount of information. Autocorrelation, with its mathematical and statistical properties, is an efficient way of making the best of this information.We compute all the autocorrelations of reflections having a given number of arches which fulfil certain conditions determined in advance. Then, after statistical testing of some parameters pertaining to the autocorrelations (abcissae of zeros, of extrema …), we select only those with a maximum likelihood. It is sufficient to consider only the reflections whose autocorrelations have been selected and to arrange them in groups according to their shape and arch number in order to obtain average pulses.In so doing several solutions are arrived at, but when considering a given number of traces, a single record for instance, it is possible by comparing these results with each other to considerably reduce their number.In the last part of the paper the nature of the impulse obtained with our method is examined in order to find out whether it is “minimum phase” for carrying out deconvolutions.

Notes: In recent years considerable work has been done to devise a satisfactory non-dynamite seismic system that would replace dynamite in offshore areas. Prior to the advent of digital recording and processing, the non-dynamite sources have generally not provided the depth of penetration or the resolution required for satisfactory seismic interpretation.More recent developments in non-dynamite offshore marine sources include adaptation of the Vibroseis from a land unit to a marine unit, and adaptation of the Dinoseis unit from a land to a marine unit. The SUE (Seismic Underwater Explorer) system is a thermodynamic non-dynamite source utilizing a mixture of propane and oxygen detonated in a special chamber approximately 15 feet below the water surface. This source gives penetration to more than 4 sec in areas typified by Gulf of Mexico type geology and shows deeper penetration than had previously been obtained by dynamite along the western United States in areas with 20 lb charge limitations. A pneumatic source, the airgun, has been in production use in the United States since June 1966. This non-dynamite source provides an intriguing amount of versatility and can be expanded to provide additional energy as necessary to obtain the penetration desired. Tests using systems comprised of from eight to twenty-three airguns show penetration in excess of 5 seconds in many areas. Power spectra comparisons both in amplitude and frequency content demonstrate that this is a controlled source generating a controlled seismic wavelet and a controlled frequency spectrum that can be tailored to fit requirements of particular areas. Sample sections obtained in the Gulf of Mexico and the Pacific Ocean offshore California show adequate penetration to 5.0 seconds reflection time.Quantitative measurements with the airguns demonstrate the effect of:〈list xml:id="l1" style="custom"〉1Variation of the number of guns in the system;2Shaping the frequency spectrum by using different sizes of airguns in the system;3Effects on signal-to-noise ratios as a result of stacking several small energy sources together;4Reproducibility of the initial pulse wavelet from shot to shot.The improvement in record quality as a result of advanced digital processing with non-dynamite sources is comparable to that obtained with dynamite sources. Non-dynamite sources make additional improvements possible where high source multiplicity is advantageous. Excellent dynamic correlations yield accurate velocity control as well as definitions of apparent velocities attributable to multiples and primary-to-multiple amplitude relationships.Non-dynamite sources are being used more and more extensively in offshore exploration. The advent of digital recording and processing provides a means for improving depth of penetration and resolution of many non-dynamite sources.

Notes: The electromagnetic field radiated from a magnetic dipole lying on the ground is considered, in the extremely low frequency range (DC to 20 000 Hertz).Theoretical and experimental data are given on the characteristics of the surface wave (vanishing wave) generated at the air-ground interface, in the case of an homogeneous subsurface.The case of a subsurface with electrical resistivity varying with depth is considered.It is shown how the above-mentioned characteristics may be applied in the quantitative investigation of the electrical resistivity of the subsurface as a function of depth, in a method using the measurements on the ground of all the components of the radiated field: horizontal electrical components, vertical and horizontal magnetic components.

Notes: Book review in this articleHandbuch der Physik/Encyclopedia of Physics, Ed.: S. FlüggeThe Aurora. By S.-I. Akasofu, S. Chapman and A. B. Meinel.Time-Variations of Cosmic Rays. By S. E. Forbush.Magnetic Properties of Rocks and Minerals. By T. Nagata.

Notes: A theory is developed to predict the amplitude and velocity of the first arrival produced by a thin solid refractor embedded in a solid environment. Laboratory experiments on a three dimensional model substantially confirm the theory. The essential conclusion may be found in Figure 8 where it is shown that for a constant wavelength, and for thicknesses up to about one half a wavelength, the attenuation per unit distance decreases monotonically with increase of refractor thickness.The assumption is made that the amplitude, A, attenuates as (x)-n exp(-βx), where x is horizontal distance and n and β are constants. It is emphasised that it is difficult to obtain reliable values for n and β from the initial amplitude measurements.

Notes: Book reviewed in this article: State of Stress in the Earth's Crust, Proceedings of the International ConferenceF. A. Vening Meinesz, The Earth's Crust and Mantle (Development in Solid Earth Geophysics)

Notes: A formula is developed for the vertical magnetic field due to an alternating current passing through a long horizontal cable placed on a two-layer ground. The variations of the phase and amplitude ratio of the vertical field, along profiles perpendicular to the cable line are investigated. Within the limits encountered in practice, it is found that if the upper layer is conductive, as compared to the lower layer, the phase may vary appreciably whereas the amplitude ratio changes as in the case of vacuum.

Notes: Book review in this articleN. A. Ogil'vi, D. I. Fedorovich, Groundwater Seepage Rates. English translation by N. H. Choksy, Consultants Bureau, New York 1966, 33 pp., 12 fig.

Notes: Mapping of fault patterns is an important part of geophysical exploration. A computerized digital template analysis method is described which tests gravity maps for the effects of faults by comparing measured gravity data with calculated master curves. The interpreted gravity data are incorporated in a tectonic map using geological symbols and units, for ready use by the geologists.Tectonics can be investigated by gravity detailing if the smallest undulations on the Bouguer map are taken into consideration. Up to now, residual and derivative gravity maps have explained gravity effects by assuming spherical bodies which are almost unknown in geology. The method discussed here uses tectonic elements, such as fault-blocks and dikes as a basis of interpretation instead of the spherical bodies of the conventional interpretation methods.Gravity data can be easily and relatively cheaply obtained in the early phases of exploration by area wide spot coverage following lines of easy access such as roads, etc. Seismic studies, by contrast, have the disadvantage of being executed along predetermined profile lines. Only after sufficient detailing do seismic profiles permit-if at all -an areal mapping of faults. Thus a tectonic map is obtained only at the end of a geophysical survey instead of being available prior to the planning of costly seismic profiles. The use of gravimetric data and their interpretation by the suggested method provides tectonic detail maps in the early phases of geophysical studies.In addition, this article discusses a general geophysical interpretation method, using the investigation of faults by gravity as an example. Applications of this method for different gravity and magnetics problems as well as for combined interpretations are outlined. Detailed case histories will be published in later articles.

Notes: To be able to make use of all existing probe interpretation techniques, it is desirable that field resistivity observations be conducted in such a way that it is possible to construct an equivalent curve for the surface variation of potential about a point source. Further, the usual parallel layer interpretation will only be justified in as much as the observed potential curve is compatible with such a subsurface geometry. Thus, whilst a potential curve may be constructed from suitable finite potential differences obtained using a Wenner configuration of electrodes, it can be shown, using the tri-potential technique of measurements, that these potential differences may arise partly from lateral resistivity variations.In this paper, a ladder network technique is employed to display these lateral effects and an adjustment method proposed to reduce them. The adjusted potential differences obtained form a consecutive series suitable for summation to give the potential function. These values are subject to further minimum adjustments required for them to comply with slope and curvature conditions for a layered medium. After forming the potential curve by summation, a final numerical smoothing process is carried out. From this smoothed potential curve, corresponding Wenner and Schlumberger curves can readily be derived for interpretation.The method proposed thus attempts to extract from the observed data the maximum part which can be reconciled with a purely depth variation of resistivity. A method is also proposed for compounding the ‘errors of closure’ of the network to provide a Lateral Inhomogeneity Index which gives a measure of the departure of the observed data from the basic interpretation requirements.The method involves no great labour, but can readily be programmed for a computer if desired. Examples are given of the application of the method to field observations. In the interest of objectivity, the final smoothing has been confined to a single stage but in certain very extreme examples a further stage may prove desirable.

Notes: Book reviewed in this article:M. R. J. Wyllie, The Fundamentals of Well Log Interpretation.

Notes: Proton precession magnetometers usually operate discontinuously since thy rely one intermittent impulses. The aerial versions of these magnetometers suffer from relatively high background noise. The type of magnetometer discussed in this paper has solved the two problems of discontinuity of measurements and reduced sensitivity due to background noise.Comparative texts with an established flux-gate apparatus have given significant and convincing results.

Notes: Book reviewed in this article:S. I. Kuznetsov, M. V. Ivanov and N. N. Lyalikova, Introduction to Geological Microbiology Organic Geochemistry, Editor I. A. Breger

Notes: The formula for the Bouguer anomaly contains the density of the masses above the datum plane. If we use a wrong value for the density we obtain on the Bouguer map additional anomalies which coincide with the topographic features.If we do not know the Bouguer density we can compute several Bouguer maps based on different arbitrary densities. Then we select that map which shows the least correlation with the topography and obtain in this way a reasonable value for the density. This method has been found by Nettleton.Several authors have formulated this process mathematically and have elaborated numerical methods. This paper will present such a method. It differs from the former methods by the following:〈list xml:id="l1" style="custom"〉1A general formula for the Bouguer density is presented, using a non-specified linear operator. By choosing special definition for this operator, we can obtain the formulae of several of the known methods.2It is shown that the determination of the density by “Nettleton Profiles” gives less reliable results than the statistical investigation of the whole area covered by the gravity survey.3The known methods do not consider the variations of the vertical gradient. It is shown that in areas with a regional gravity anomaly or an extended topographic feature the known methods yield poor results.4It is shown that the vertical gradient of gravity used as operator gives a most exact and simple formula.5It is shown how the electronic computer—after reading the Bouguer values based on an arbitrary initial density—can compute the residuals, derivatives etc., which are based on the corrected Bouguer density.

Notes: An electronic analogue was used to obtain geophone distortion of typical seismic pulses. Some of the results were checked by computation and in all cases agreement was extremely close. For late arrival events estimates of attenuation and phase shift may be made from the steady-state response curves. It is explained why these curves are not applicable to the first peak, for which one must use the results in Table III.Ground-geophone coupling introduces the same distortion as a very lightly damped galvanometer with a natural frequency usually in the range 100c/s—200c/s. In swampy terrain this frequency may reduce to 30c/s—40c/s and coupling effects may be severe.An electronic unit whose characteristics are the inverse of those of a geophone has been used to improve the first peak amplitude of a seismic refraction record. Since 1c/s geophones are often desirable in refraction survey it may well be more economic to use higher frequency geophones followed by their inverse filter.

Notes: The paper describes a new method for the calculation of standard-graphs i.e. resistivity-curves for assumed resistivity profiles. It starts from the integral of Stefanesco and arrives at a system of image poles, similar to, but simpler than earlier developed systems. The advantages of the new method are described, the accuracy is discussed and it is indicated how the accuracy can be improved by a simple artifice.

Notes: In the present study an attempt is made to relate the general pattern of reflection seismograms obtained over a coal basin to the acoustical properties peculiar to the Coal Measures. The Carboniferous is characterized by a very rapid succession of distinct lithological units, resulting in a surprising variability in their physical properties. Coal seams in particular are conspicuous by their low density and velocity and they produce acoustical contrasts of 35–50 % with respect to country rock. The resulting heavy shielding action prevents the seismic energy from penetrating deeply into this type of formation.With the aid of synthetic seismograms, arranged into synthetic profiles, it is shown that only a minor part of the Coal Measures reflections is of primary origin. The majority has a secondary character brought about by reverberations inside the weathering zone. Thus the conclusion is reached that reflection shooting allows of mapping only the structure of the top section of the Carboniferous, which is sufficient for the purposes of the coalmining industry.

Notes: A new logging method, which makes use of specific chemical reactions to locate petroleum or other minerals has been developed. The method has been named “Differential Electric Log” (CED).In this paper the application to the CED method of redox (reduction/oxidation) reactions is discussed. A series of organic oxidants was found, capable of undergoing reaction with crude oil, sulphur, mineral sulphides, lignite. A small amount of an oxidant of this type is added to the drilling fluid and pumped into the borehole where reaction takes place in correspondence with layers containing oxidizable minerals. Two SP logs, recorded before and after the oxidation, are compared and the presence of such minerals is thus detected.A theoretical hypothesis explaining the mechanism of the potential change due to the chemical reaction is proposed both in the case of electronically conducting and non conducting minerals. The influence of redox potential on SP is discussed.The results of laboratory experiments on several different oxidizable ores and those of field tests performed on pyrite, lignite, and oil wells are reported on and discussed on the basis of the theoretical hypothesis proposed.

Notes: As an introduction the various ray paths of a refraction wave, which is reflected at a fault, are discussed for the case of an arbitrary angle between the refracting horizon and the fault. Simple geometric considerations lead to the conclusion that the best chances for recording these pulses are encountered, if the angle between the refracting horizon and the fault is either 90° or the critical angle of refraction. In both cases identical travel times of the pulses are to be expected.The case of a fault perpendicular to the refracting horizon is considered in detail for dipping beds. Formulas for the shot point travel time curve and the time contour map are derived. Computed time contour maps show considerable differences between the direction of strike of the contour lines and the strike of the fault, as well as between the recorded apparent velocity and the true velocity of the refracting horizon. Finally, alignment charts and computing procedures are given by which the position of the fault and the velocity of the refracting horizon can be obtained from the recorded shot point travel times or the time contour map.

Notes: It is not easy to record a seismic pulse at distances of interest to oil prospectors without there being reflected or refracted pulses superimposed on the direct arrival. Accordingly the record illustrated is considered worth publishing, although it was taken fortuitously during a normal survey and cannot claim to be a controlled experiment. A comparison with the filtered pulse to be expected from a theoretical Ricker-type wavelet is presented.

Notes: In this paper a method of calculating synthetic seismograms from continuous velocity logs (CVL) is presented and described in detail.The seismic velocities on the CVL are frequently and rapidly changing. Therefore a seismic field-record consists of a great number of reflections with different sign and amplitude. The sum of all reflections forms the seismogram.The first part deals with some relations between the log, the reflection-layers, and the normal field-record, with special emphasis on the influence of the seismic amplifier.Then several assumptions are made for the construction of the synthetic seismograms. Finally some examples and applications are given.

Notes: The Karroo of Madagascar is a thick continental deposit of sandstone with some shale, but no characteristic beds. Where it outcrops, velocities have been measured through bore hole surveys. It seems that they depend primarily on the age, but long refraction soundings show a progressive increase of the velocity with depth and a weak differentiation of the marker beds. These results may possibly be brought in accordance by discriminating between the present position and the deepest position of a bed through its geologic history. Thus we must admit very important erosions which, according to regional geologic data, seem possible. If we make the appropriate correction due to this distinction, all the velocity measurements fall near the curve proposed by Faust for an average shale and sand section.

Notes: Some new formulae for the rapid estimation from the gravity field of the maximum possible depth to the top surface of gravitating bodies of positive mass are stated, discussed and proved.

Notes: The shothole drilling industry has had new developments and its evolution is tied to several causes. First, more and more difficult and varied formations have to be drilled. Second, new procedures have to be used in desert regions recently open to prospecting. Finally, modifications of the seismic methods themselves (multiplication of the charges, shallow holes, etc.) lead the driller to use more economical and lighter tools. In this paper the following types of drilling equipment are described: the Mouton Delmag for spudding through pebbles and loose gravels; a mud mixing water truck and an extensible mast used for drilling multiple medium depth holes through caving sands; a compressed air-drill for regions where the water supply is practically nil; rock drills, such as a Wagon-Jack and the Warsop hammer, for multiple shallow holes in various formations.

Notes: A strong late arrival, with several cycles, has been observed in line refraction shooting, for distances ranging from 5000 to 30,000 feet. Using equations given by Officer (1953) it has been possible to establish that this event is a multiply reflected refraction. The conclusive factor was its very large amplitude which was about 60 times that expected for the simple head wave, at the distances involved. The wave-guide was bounded by the surface of the earth and by the base of the Low Velocity layer, at a depth of about 80 feet.An earlier paper (O'Brien, 1957a) gave a study of the head wave pulse refracted from this interface and the conclusions in the two papers are altogether compatible.