ALBERT

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

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

Notes: In mathematical statistical filtering the deconvolution problem can be solved by two different methods:〈list xml:id="l1" style="custom"〉1by inverse filtering2by calculating the prediction error.Both methods are well known in the theory of Wiener filters.If, however, the generating process of the signal is known and can be described by a set of linear first order differential equations, then the Kalman filter can also be used to solve the deconvolution problem. In the case of the inverse filtering method this was shown by Bayless and Brigham (1970). But, while their method can only be used if the original signal is a colored random process, this paper shows that in the case of a white process the prediction error filtering method is a more appropriate approach. The method is extremely efficient and simple. This can be demonstrated by an example which maybe of special interest for seismic exploration.

Notes: The use of the electrical image method in calculating the potential of a source in the presence of a buried conductive sphere proves to be difficult because of the great number of terms making up the potential formula. The disadvantages contained by this method are avoided here by resorting to the theory of graphs. In view of this aim the “propagation” of the electrical images due to the sphere has been represented by an infinite graph to which an infinite matrix has been attached.Taking this as a basis several FORTRAN IV programs for the computation of resistivity pseudo-sections for various electrode configurations have been coded, and a great deal of such pseudo-sections have already been calculated with an IBM-360 computer. Using these results, recommendations regarding both the methodology of the geoelectrical prospecting of spherical shaped conductors and the most suitable way of building resistivity pseudo-sections are developed.

Notes: An account is given of the use of the square array technique in investigating the surface effects of rotational anisotropy when the axes of anisotropy are inclined to the surface. It is shown that, as with other arrays, two anisotropy parameters 〈inlineGraphic alt="inline image" href="urn:x-wiley:00168025:GPR249:GPR_249_mu1" location="equation/GPR_249_mu1.gif"/〉 and n can be derived by varying the array orientation.On the basis of these considerations, the effects of such anisotropy on the values of the mean apparent resistivity and azimuthal inhomogeneity ratio normally obtained in square array measurements is reviewed.Particular attention is paid to the variation of resistivity with orientation and it is noted that, in areas of moderate anisotropy, this variation is lower for the square than for the Schlumberger array.In addition to this advantage, the azimuthal inhomogeneity ratio obtained from square array measurements may be used to indicate the severity of anisotropy in an area and two field examples of this use are given.Where anisotropy is severe, gross variations of apparent resistivity with orientation are obtained with either square or collinear arrays. In these circumstances, the use of crossed measurements is considered and the particular stability of the crossed square array demonstrated.

Notes: Electrical resistivity sounding was attempted with success in 1960 on the surface of the Nile water in Khor Kundi El-Bahari, the starting place for the construction of the offset channel of the High Dam. The method was applied to determine the depth to the upper surface of the granite below the alluvium and Nile water. The success of the method conducted at two stations in the Khor is attributed to the favourable geological conditions of the section which consists of nearly homogeneous layers besides the good earthing conditions of electrodes dipped in water. The high resistivity contrast between water and granite facilitates interpretation of sounding curves. Factors interfering with the results of interpretation are mentioned.

Notes: This paper presents a dot chart for calculating gravity effects from two-dimensional bodies with the addition of a supplementary set of curves which can be used to determine the “end corrections” to modify the calculated effects when the length of the body is not relatively long with respect to the horizontal dimensions of the cross-section. The use of the chart is demonstrated by different examples. The results of theoretical cases amenable to mathematical treatment and computer results both compare favorably to those computed using the present chart.A fundamental incentive of constructing this chart, in a final form ready to use, is to make it available to gravity interpreters for utilization. Future effort in chart construction is thereby avoided.

Notes: By applying dynamic corrections a seismic trace recorded at a distance x from the energy source should be varied in such a way as to obtain a trace which would be recorded at zero-distance, i.e. at the source itself. Only such a zero-offset-trace contains the correct sequence of reflection coefficients (reflectivity function), whilst all other traces contain a distorted reflectivity function. In the simplest case, the reflectivity function is compressed over a shorter time whereas in more complicated cases a partial inversion of the reflectivity function results. This happens when some of the reflection hyperbolae intersect one another.The reconstruction of the true zero-offset reflectivity function by the application of dynamic corrections can only be an approximative process. In the first case mentioned we must expect a decrease in accuracy of the corrected trace in comparison with a zero-offset-trace. In the second case, where intersections of the hyperbolae occur, accurate reconstruction is clearly impossible.The problems are discussed with the help of theoretical and practical examples.

Notes: Fourier transform techniques have been used to calculate the theoretical filter (amplitude) response function of Nth order vertical derivative continuation operation. The amplitude response functions of the vertical gradient and its continuation follow from the same. These response functions are subsequently used to calculate the weighting coefficients suitable for two dimensional equispaced data. A shortening operator has been incorporated to limit the extent of the operator. For comparative study, some of the developed coefficient sets and the one presented in this paper are analysed in the frequency domain and their merits and demerits are discussed.

Notes: In this paper it is shown how one may obtain a generalized Ohm's law which relates the induced polarization electric field to the steady-state current density through the introduction of a fictitious resistivity defined as the product of the chargeability and the resistivity of a given medium.The potential generated by the induced polarization is calculated at any point in a layered earth by the same procedure as used for calculating the potential due to a point source of direct current.On the basis of the definition of the apparent chargeability ma, the expressions of ma for different stratigraphie situations are obtained, provided the IP measurements are carried out on surface with an appropriate AMNB array. These expressions may be used to plot master curves for IP vertical soundings.Finally some field experiments over sedimentary formations and the quantitative interpretation procedure are reported.

Notes: In the linear filter method of interpreting resistivity sounding data, as developed by Ghosh (1971), it appears that the filter function in the x-domain approaches an oscillating function for both large positive and large negative abscissa values. In the present note the reason for this oscillating behaviour is derived, and a possible practical application is indicated.

Notes: Amplitude spectra of input FM signals used in the vibratory source method of seismic exploration often show undesirable oscillations near the initial and terminal frequencies. These oscillations have an effect on the correlation background and distort the output signal. Considerable improvement in reducing the amplitude of these oscillations is obtained using a proper taper fuction. Attention is given to the relation between the tapering time and bandwidth of the spectrum.Analyses of the spectra of the received data from vibratory sources show considerable attenuation in comparison with the original field sweep. Since the matched filtering process will result in a series of waveforms which have the shape of the autocorrelation of the input signal, consideration is given to the autocorrelation function and its zero-lag coefficient of the FM signal in the presence of attenuation. A method has been developed which compensates for the attenuation and recovers the distortion of waveforms when the received data is correlated.The design of a waveform shaping filter for vibratory source data is given to reduce the influence of phase distortion on the received waveforms as well as to increase S/N ratio resolution. Parameters used for this filter are based on the properties of the FM signal and its autocorrelation function.Several examples from field data are presented to illustrate the methods. The results indicate that the use of the above techniques yields sections with good frequency resolution and improved S/N ratio.

Notes: The problem of a plane wave incident on a non-isotropic dipping layer lying over an isotropic conducting substratum has been studied and some numerical results are presented to show the effects of〈list xml:id="l1" style="custom"〉1) degree of anisotropy m,2) conductivity contrast between the substratum to the upper layer b,3) angle of inclination of the axis of anisotropy α,on the variation of amplitude and phase difference of the wave impedance with the normalized skin-depth (d/δ1).

Notes: Consideration of the spectral representation of gravity and magnetic fields shows the field to be the result of the convolution of factors depending on the parameters of the field. Wiener filters, calculated using model transformations, provide an optimum method for altering these factors and hence effecting field transformations.

Notes: A discussion is given of the requirements, the advantages and the methods to be considered in attempting the quantitative interpretation of gravity or magnetic fields from computed maps of the vertically derived field.The transform which is used here is the first vertical derivative (or vertical gradient) with or without downward continuation, but the computed maps are in fact obtained by a controllable Fourier method in which two kinds of operations can be simultaneously performed in complete independance: on one hand the separation of any part of the data by frequency cut-off, and on the other hand the transformation by vertical derivation or continuation of the part which is retained.Taking as raw data either actual surveys or artificially constructed maps, it is first shown how separated and transformed maps of this type can be efficiently obtained under quite flexible conditions, using a special computer program. It is further seen that for correctly controlled filterings the accuracy of the computed maps actually permits to take them as the basis for quantitative interpretation.To effect this, any one of the conventional methods which make use of equivalent model computations may in the first place be adapted to the interpretation of gradients, with the benefits, however, of an enhanced lateral separation of the anomalies and of a large attenuation of the regional effects.Particularly, the delineation of horizontal contours for even fairly complex models can often be made directly in a sufficiently safe way on the anomalies as they show on the vertical gradient maps. This greatly accelerates the process of determining equivalent model bodies.More special methods of interpretation can also be designed by taking into account first the fact that the vertical derivation of the field amounts to an operation of separation on the field's sources themselves, and in addition the availability of the frequency form of the information as a result of using a Fourier method of transform computation.Trial utilisations of various such interpretation processes of either the conventional or the less conventional type, are presented especially in connection with an interpretation study on transformed maps of parts of an offshore aeromagnetic survey (English Channel).

Notes: Generation of electrical power through the use of natural steam obtained from the Earth has the potential to be a significant source of pollution-free power. Successful future exploitation of geothermal power depends upon the development of cost-effective exploration techniques. A field experiment in the Imperial Valley of California indicates that there is a close empirical relationship between hot water deposits at depth and an anomalously high seismic background level at the surface. If this relationship proves to be a general one, it will serve as the basis of a simple and inexpensive method of geothermal exploration.

Notes: In massive rocks ground waters mainly flow in fracture zones whose permeability greatly changes depending on their filling. When studying ground-water flow in fissures, the results of observations of the electric fields of filtrational origin—which, in this case, considerably differ from those in porous media—can be used. Therefore the authors have made experiments on fissured media models. The measurements have been made in a special filtration tube with the fissured media simulated by a system of quartz glass plates. The spacings between plates were regarded as fissure widths. The observations have been made in fully open “fissures” and in those partially filled with sand or sandy-clay material.These experiments have resulted in establishing a dependence between the values of streaming potentials and pressure drops. The SP values have particularly been found to decrease with the opening of “fissures”. The most intensive electro-filtrational fields were observed at 40 per cent filling with medium grained sand. Additions of argillaceous material to the sand filler brought about sharp reductions in the intensity of electro-filtrational potentials.

Notes: During the last few years the airpulser, or air gun, has become very common as an energy source for marine seismic surveys. This paper describes the physical processes which take place during the operation of the pulser and develops theoretical results concerning the energy and frequency of the radiated signal and the amplitude decay of the secondary bubble pulses. The theory takes into account the presence of the airpulser itself which is assumed to be a rigid sphere within the bubble of released air. The theoretical results are combined and compared with measurements made of the pressure within the airpulser, the acceleration of the body of the pulser, and the amplitude and frequency of the signal radiated into the surrounding water. A formula for calculating the bubble frequency is given and a diagram made of the energy partition between mechanical losses, radiated energy, etc. Finally, a comparison is made of the energy release from the airpulser with that from TNT.

Notes: An equivalent stratum model is used to provide an explanation for the power spectrum characteristics of potential fields.The power spectrum of observed magnetic fields is found to consist of effects due to sources which can be represented by apparent monopoles and dipoles. Similarly, the power spectrum of observed gravity fields is found to consist of two groups of effects. A matched filter is proposed whereby the effects due to the two components in the potential field can be separated. As a consequence of the established theoretical expressions a scheme is suggested whereby the aliasing power of sampled data can be estimated. Also, the concepts of downward continuation, reduction to the pole, and reduction to pseudo-gravity of magnetic fields are re-examined in light of the theoretical expressions due to the equivalent stratum model.

Notes: The oscillatory character of the solution obtained by Qureshi and Mula's method for the direct determination of two-dimensional mass distributions from gravity anomalies is examined.

Notes: A number of resistivity profiles over dolerite and metadolerite dykes have been taken using Wenner and Schlumberger configurations around Dhanbad in order to find the potential of these for water accumulation. These dykes give rise to both resistivity highs and lows as compared to the background. The results are discussed in the light of known theories. Some suitable sites for water accumulation have been delineated.

Notes: Differential and difference soundings of resistivity for real, imaginary and complex values of the newly introduced parameters ν and μ are considered. The transverse and longitudinal differential soundings, known up to the present, correspond with the values ν= 1 and ν=— 1.The point of departure is the generalization of the concept of Dar Zarrouk parameters and the dimensional analogy between D.Z. resistivity and transverse and longitudinal differential resistivities. Properties of generalized D.Z. curves as well as the method of their construction are given. The examples of differential and difference curves are presented and the properties of these curves are discussed; particulary the effective spacings are determined. The problems of the realization of differential and difference soundings are considered; the bisymmetrical and asymmetrical arrays are proposed as a practical ones.

Notes: The effect of a conducting overburden on the electromagnetic response of sulphide ore-bodies has been studied with the help of quantitative electromagnetic model experiments. These experiments were conducted at a fixed, crystal controlled, frequency of 100 kHz using a number of transmitter-receiver configurations, though the results discussed here mostly pertain to a horizontal coplanar system. An analysis of the anomaly profiles—after accounting for the regional anomaly—indicates a general enhancement of the response. This is attributed (a) to the concentration of the current lines in the host-rock effected by the embedded target and (b) to the change in the phase and space orientation of the field vectors, brought about by the conducting overburden. The enhancement is most pronounced in the case of non-symmetrical bodies such as veins and ribbon-like structures represented by sheet models, but is also observed to a lesser degree for isometric structures. These results are expected to pave the way for a more accurate interpretation of the induction prospecting data.

Notes: Mit der vorliegenden Arbeit wird das Echo-Log als ein geophysikalisches Meßverfahren zur Bestimmung von Form und Volumensverteilung von Aussolungshohlräumen im Salzgebirge vorgestellt. Hierzu werden neue terminologische Begriffe wie Echo, Echostrahl, Echostrahler, Echozeit, Echogramm usw. erläutert. Weiter wird über die Meßausrüstung, Durchführung der Vermessung sowie die Ergebnisdarstellung berichtet.In einem letzten Abschnitt wird zwischen Echos, die auf einem direkten Weg von der Hohlraumfläche zurückgeworfen werden und anderen bei der Registrierung erfaßten Empfangssignalen unterschieden.In this paper the Echo-Log is presented as geophysical measuring method for determining form and volume spread of leached cavities in salt formations. Therefore, new terminological definitions like Echo, Echo Beam, Echo Transducer, Echo Time, Echogram etc. will be explained. Furthermore, also the measuring equipment, execution of a survey and the presentation of the results will be described.In the last paragraph the difference between Echo directly reflected from cavity walls and other in the cavities registered signals will be discussed.

Notes: The theory of potential field continuation is studied with a view of extending it to include continuation between non-linear surfaces. This theoretical extension is regarded as useful in reducing observed potential fields along one general surface to another general surface.It is demonstrated that the continuation operators considered in past geophysical literature are special cases of the generalized operators. In view of the tremendous growth in the art of high-speed computing it is possible to consider applications of the generalized operators. In worked examples upward continuation between a general surface and a datum is considered for the gravity field due to geometrically simple sources; this in order to test the accuracy of the digital applications. It is indicated that for observations on an undulating surface it is possible to account for considerable errors when, during interpretation, the observation locations are taken to be along a datum. For aids to profile interpretation in areas of considerable topographic relief some simple convolution operators for continuation between linear, non-parallel surfaces are proposed.

Notes: The ground follow-up of a magnetic and radiometric air survey had to cope with the usual dilemma to check a great number of anomalies within a short time. A limited aggregate of magnetic anomalies, expected to correspond to magnetite quartzites was statistically selected for ground identification by this method: the ratios length/width times amplitude were listed for all coherent contours and the calculation of the standard deviation per unit area resulted in different key-numbers for a new lithological “Salem Unit” and for the charnockitic or gneissic environment. The ground work thus directed and reduced by 85% yielded a substantial potential of iron ore. This was supported by by abundant determinations of the magnetic susceptibilities, confirming the sources of anomalies and revealing the amenabilities of iron ores for the magnetic separation process.The lines of truncation of anomalies were found to represent a system of local and regional faults and shear zones, which segregated the area into different tectonic blocks. These sutures have also provided the ways of intrusion for alkaline and basic magmae in the style of a “Rift” structure, housing several carbonatites and impregnations of metal sulphides.The airborne radiometrics obtained many uranium indications by gamma ray spectrometry. However, they led only to disseminated uranium-silicates, associated with syenites, granites and pegmatites, each emanating a characteristic photon energy spectrum. But in general the radiation of thorium prevails, marking northern Madras as a “Thorium Province”.

Notes: In the interpretation of magnetic anomalies and in paleomagnetism, the anisotropy of magnetic susceptibility is commonly neglected. Nevertheless, this property has basic significance, because, owing to susceptibility anisotropy, the directions of the vectors of induced and remanent magnetization are deflected from the direction of the Earth's magnetic field. Almost all rock types investigated possess higher or lower degree of the susceptibility anisotropy. Effusive and sedimentary rocks have the lowest degree of anisotropy. For the latter, the “masking effect” of the paramagnetic mineral components has some influence on the anisotropy degree due to the low mean susceptibility of sedimentary rocks. Metamorphic and plutonic rocks usually exhibit a considerable degree of anisotropy. The highest degree of anisotropy has been found in the rocks containing ferromagnetic minerals with mimetic fabric.The dependence of the degree of the susceptibility anisotropy on the degree of metamorphism proved to be very complicated; of the rock sequence from slates to gneisses, the transient rocks (roofing slates and mica-schist-gneisses) showed the highest degree of anisotropy. This result can be used in geology for reliable determination of these rock types.

Notes: A technique utilizing the convolution theorem is developed to obtain the time-domain electromagnetic response for a number of input pulses from the given step-function response. For illustration an idealized model of a homogeneously conducting non-permeable sphere placed in a uniform field is considered. The nature of the responses due to different types of pulses and their usefulness in estimating the constitution parameter of the sphere are discussed. To show the applicability of the suggested technique to generalized systems, calculations for a conducting permeable sphere are presented in the appendix.

Notes: De préférence, les sismographes seront disposés dans un lieu sans bruit, oú il n'y aurait pas besoin de patternes. Par contre, si les bruits sont tels qu'une patterne est nécessaire, l'objectif de la patterne est d'améliorer le rapport signal/bruit et donc de rendre le signal le plus utile possible. La patterne sera dessinée en fonction des charactéristiques du signal, des charactéristiques du bruit, et de la géologie du site.Il a été démontré que, dans le sens pratique, l'opération optimum sur les données d'une patterne est représentée par le proces “delai et somme”. L'augmentation du nombre N de senseurs sur une surface donnée diminue les espacements entre les senseurs et peut augmenter la coherence entre les bruits enrigistrés aux senseurs voisins, donnant donc moins d'amélioration au rapport signal/bruit que la √N anticipée. L'augmentation du nombre de senseurs par l'augmentation proportionnelle de la surface de la patterne peut resulter en la détérioration du signal, et peut également donner moins d'amélioration que la √N Ces deux effets, ainsi que l'élément économique, combinent à limiter le nombre de senseurs qui peuvent être employés.Bien que les données sur lesquelles l'on a basé ces conclusions soient tirées de la sismologie conventionnelle, ces principes sont également valables pour l'exploration sismique et pour d'autres mesures géophysiques qui emploient des dispositifs dans lesquelles plusieurs senseurs sont requis.

Notes: Recent improvements in instrumentation and field techniques have made time domain electromagnetic methods more acceptable. This acceptance has prompted further theoretical work for use in the interpretation of field data.The asymptotic solutions for the transient electromagnetic field components in the vicinity of a fault zone separating two media of high resistivity contrast are obtained for low frequency or late time. Excitation is by normally incident plane waves at the earth's surface. Because of the slow convergence of the asymptotic time series expansions, a numerical polygonal inversion is performed on the real part, or in-phase term, of the time-harmonic surface expressions for the electric and magnetic field components.For both impulsive and step excitation the transient electric field normal to the contact is more sensitive to changes in the structural attitude of the fault plane than the transient electric field parallel to the contact. The transient anomalous vertical magnetic field for either impulsive or step excitation appears to be most diagnostic of dip angle, although waveform shape does not seem to be significantly dependent on the slope of the fault. For dip angles greater than 90 degrees, as measured on the poorly conducting side of the contact, all field components become more insensitive as indicators of the structural attitude. The results presented here should be useful in obtaining several geologic parameters descriptive of a fault zone or lateral resistivity inhomogeneity from transient electromagnetic soundings; they should also provide an aid to differentiating with available geological information between layering effects on transient electromagnetic responses and effects largely due to lateral changes in resistivity.

Notes: Velekin and Bulgakov (1967) in an interesting model experiment while studying the transient electromagnetic response of a conductive sphere placed below a thin conductive sheet found that at the earlier stages of the transience, the composite system response corresponded to the response due to the overlying sheet alone and at the later stages, it corresponded to that of the sphere alone. To examine whether such a separation of responses due to individual components can be analytically studied and applied to other source configurations, we have analyzed an idealized model consisting of two spherical shells. We find that in corroboration with the above results, the general nature of the curve consists of two humps representing the responses dominated by the outer shell and the inner shell respectively. In addition, however, we find that the two humps gradually disappear to yield a smooth decay curve for increasing values of the ratio σd1b/σ2d2a (where σ1, σ2 are the conductivities, d1, d2 are the thicknesses of the outer and inner shells respectively, and b and a are their respective distances from the centre) and the effect of inner shell on the composite system response is considerably reduced.

Notes: It has been necessary to resort to the use of “long-line” refraction marine operations in certain areas where it proved impossible to eliminate singing from reflection records despite the number and variety of programs at our disposal for this purpose.Experience has shown that manual processing of offshore refraction records takes a disproportionate length of time in comparison with the surveys themselves, and this is incompatible with the requirements for choosing the site of an exploration well. It thus became necessary to find an “industrial approach” to the solution of this processing problem.It was apparent that automatic picking could also facilitate the interpretation of land refraction data, and that in the case of both marine and land work the interpretations would be more accurate when factors were taken into account which could not be considered when working without the aid of a digital computer.For these reasons a set of programs was developed for automatic picking and interpretation of refracted arrivals.The picking itself consists in searching for the maximum values of the normalized cross-correlation functions of the traces with a “model” trace. The first results thus supplied are: “picked” times, intercept times, maximum values of the correlations, and the values of the tie constants between overlapping spreads.Next, the construction of the relative intercept time curves is performed; a statistical analysis of these curves then allows the determination of the offset distance.From these elements,⊙ either the delay time curve is produced, after ensuring correct reciprocal times by means of additional minor corrections.Such work is carried out in order to enable the geophysicist to gain a sound idea of the quality of the interpretation. To assist in this aim, part of the trace on both sides of the pick is plotted on the final documents. Valid groupings of several traces involving the same amount of refraction data are thus possible.⊙ or the refractor depth is constructed with the wavefront method, making use of the relative intercept times.Such a procedure, which is normally applied to first breaks, can also be used for later arrivals exhibiting slight interference and should represent an important step towards processing secondary arrivals with high interference.The development of this package, in response to a need which is shared by both SNPA and CGG, is the result of the joint efforts of the Geophysical Group of SNPA's Pau Research Center and CGG's Technical and Scientific Departments.

Notes: This paper is concerned with an application of the Kalman filter to integrated navigation systems. It is shown that this filter is well adapted to resetting the dead-reckoning position by means of a radio-hyperbolic system or of a Doppler sonar. The algorithm can also be used for resetting the dead-reckoning position by means of a satellite navigation system.The first part is briefly concerned with stochastic processes and system theory. The canonical equations for navigation problems are written. The second part is devoted to the Kalman filter. The third part discusses the application of the algorithm to navigation problems and presents some simulation results.

Notes: The expression for the potential due to a point source of current, placed on the surface of a horizontal three-layer earth, has been derived when conductivity in the second layer varies linearly with depth and changes abruptly at the boundaries. Master curves for Wenner and Schlumberger configurations have been presented for an insulating basement and for one value of conductivity gradient in the second layer.

Notes: Using a method and definition given earlier (Roy and Apparao, 1971), this paper computes the depths of investigation in homogeneous ground for (a) the Wenner α, β and γ configurations, (b) the three electrode system and (c) the dipole-dipole arrangements when the dipole lengths are not infinitesimally small. The results for (a) and (b) have been summarised in a table, while those for (c) are shown as contour diagrams. In all the dipolar arrangements examined in this paper, except the equatorial, the depth of investigation decreases (and the vertical resolution increases) with increase in any or both of the dipole lengths. For the equatorial set up, this decrease (or increase) is very small.

Notes: The unsuitability of a particular method of computation of second derivative to varying types of field problems evinces the need of a method by which it may be possible to obtain a wide range of filters. This requirement is achieved by incorporating a mathematical smoothing operator, e-λρ2, λ and ρ being the smoothing parameter and radial frequency respectively, in the theoretical frequency response of second derivative operation. A few sets of weights for various smoothing parameters, are presented. The analysis of Bouguer gravity anomaly map of Los Angeles Basin, California, U.S.A., is carried out by Elkins's (1951) and Rosenbach's (1953) methods and the sets of weights given in this paper which clearly establishes their filtering equivalence.

Notes: The computation method described in this paper is based on the existence of a linear relationship between the mutual coupling ratio and the kernel function in the integral expression for it. Accordingly, the mutual coupling ratio can be determined by first computing sample values of the kernel function and then subjecting these to a digital linear filter. In the present paper the appropriate sampling distance is determined and the values of the digital filter coefficients are computed, both for electromagnetic sounding with horizontal coils and for electromagnetic sounding with perpendicular coils.

Notes: Design procedures and characteristics of three stacking filters are discussed which may find application in various three-dimensional velocity filtering problems. These filters are derived in the time-domain as optimum multichannel Wiener filters. Random stationary functions are assumed as stochastic models for the seismic traces. All power and crosspower spectra which are the basic elements of the multichannel normal equations are statistically averaged according to specific three-dimensional considerations. Various properties of the input traces may be incorporated in the design of the optimum filters. With fairly general assumptions about the input these filters are deterministic in the sense that they are applicable to a broad class of input traces with similar statistics in amplitudes and arrival times of signals.

Notes: An airborne magnetic survey of an area of 580 km2 centred on Gosses Bluff was flown in 1968. This survey formed part of a joint project by the Australian Bureau of Mineral Resources and the United States Geological Survey to evaluate the hypothesis that Gosses Bluff is of impact origin.Analyses of the magnetic data show that shallow sources of magnetic disturbance flank the Bluff on its southern side. Interpretation of the magnetic anomalies indicates their source to be shock-melted breccias of possible Jurassic age.

Notes: The interpretation of the gravity anomaly on a horizontal plane over a causative body having circular symmetry about a vertical axis is considered from a rather unconventional approach. As the analytical expression for the gravity effect of a circular body assumes a closed form only on the axis of symmetry, the interpretation in this approach is carried out with the anomaly profile along the axis—which leads to simpler and faster computation.A numerical method is developed for computation of the anomaly profile along the vertical axis from the horizontal radial profile of the symmetric anomaly by upward continuation. Provision is also made for an end correction when the radial profile has only a limited extension.Some simple geometrical shapes are assumed for the causative body. Its parameters are then determined from least squares fitting of its gravity effects to the observed (upward continued) vertical profile (i) by the steepest descent method and (ii) by the Newton-Raphson method.Some applications of these methods are demonstrated.

Notes: In comparing the similarity between time, frequency, and phase measurements as used in the induced polarization method of geophysical prospecting, parameters must first be defined and the basis for comparison established. The conditions for mathematical equivalence then follow in a natural sequence. Laboratory measurements are made on a variety of rock samples to indicate the type of comparison and correlation expected from field measurements.Results indicate that frequency domain and phase measurements are equivalent, but two frequency domain measurements are needed to produce the same amount of polarization information as a single phase measurement. Frequency and time domain measurements will be similar but never equivalent due to the basic differences in their respective parameter definitions.

Notes: From seismic surveys zero offset reflection times and root-mean-square velocities are obtained. By use of Dix-Krey's formula, the interval velocities can be calculated.If no well velocity survey exists, the interval velocities and T(o) times are the only available information. The suggested way to get a regionally valid velocity distribution is to select N“leading horizons”, where a major change in the velocity parameters occurs and to compute the parameters of the selected velocity depth function (in most cases linear increase with depth) by a special approximation for the interval between two adjacent “leading horizons”.Herewith all reflection horizons within the interval are taken into account.

Notes: A frequency-domain analysis is outlined for a conducting sphere in a uniform Input field: inequispaced alternating half-sine wave pulses. The Barringer Input air-borne electromagnetic exploration system uses such source fields. Theoretical profiles of Hρ(t), the horizontal magnetic component over the sphere from different elevations and for various conductivity and geometrical factors are presented. Based on these results some useful features such as penetration and detectibility are discussed.

Notes: In this paper an idealized valley of a semi-elliptic cross-section is considered. For a Schlumberger configuration on the axis, sets of master curves are calculated for the ratio of semi-axis a/b= 1, 2, 3, 4, 5 in corresponding to various resistivity ratios of surrounding rocks and valley sediments. For small resistivity ratios, these model curves have the shape of three-layer curves for horizontal bedding and are often equivalent to them within the accuracy of measurements. The axial ratio a/b considerably affects the depth determination of valley sediments. In the special case of a circular cross-section (a/b= 1) the influence of the position of the electrodes on the sounding curve is studied in more detail.The application of the master curves in practice shows that the influence of the specific shape of the valley on soundings should not be neglected. In general, the valleys have a greater “true” depth than can be seen from the interpretation of a sounding by master curves corresponding to the horizontal bedding.

Notes: A method is presented for determining the lower surface of a two-dimensional body producing a magnetic anomaly when its upper surface and the intensity of magnetization are given. The magnetization vector is assumed to lie along a specified direction but the sense of magnetization may be different in different vertical sections of the body and may be regarded as unknown. The method is illustrated using the computed anomalies of some theoretical models.

Notes: Gravity surveys in western Texas with station spacing of about 400 m were complemented by carefully observed vertical gradients of gravity making use of specially designed and automated instrumentation. The presented areal surveys of 51 and 52 stations taken in relatively flat terrain are parts of a large survey of close to 1000 stations complemented by an equal number of vertical gradient measurements. Quite irregular anomalous vertical gravity gradients surpassing 10 microgal/m were often encountered. Assuming the causative density contrasts to be located not more than 10 m below ground surface then free air correction errors of ± 0.1 mgal and more must be accepted. From a practical point of view there seems to be no other way to control such unpredictable errors than to carefully observe the local variation of vertical gravity gradients with adequate equipment. Making use of very closely spaced gravity measurements to derive these anomalous features seems more costly and cumbersome.

Notes: The electrical behaviour of a polarizable sample formerly studied (Time-Domain) through both the classical and the condenser methods is simulated using the theoretical responses given by a three dimensional physical-and-mathematical model, which has been proposed in a previous paper by one of us on the basis of certain hypotheses.The apparent capacitance responses of the sample are qualitatively simulated in order to tackle the problem under this point of view; moreover, an attempt to obtain a quantitative simulation of experimental data within the limits of the present reliability of the model, using theoretically obtained Model Master Curves characterized by dimensionless parameters, is done.This research is carried out in the frame of a broader study referring to a possible scheme of simulation of the so-called induced ionic polarization, concerning a number of idealized pseudo-reversible, potentially irreversible, non-linear phenomena.

Notes: Velocity analysis of reflection seismic data has increased in popularity with petroleum exploration personnel to such an extent that it now constitutes a significant tool in the routine digital processing program. One drawback to such an analysis is its computational cost, resulting from the numerous mathematical calculations required by a digital computer to extract comprehensive velocity information.The subject algorithm affords a computational cost saving of approximately one order of magnitude over a straightforward summational type velocity analysis with little sacrifice in accuracy. The utility of the algorithm is demonstrated on some Offshore Louisiana seismic data and a comparison of the results reveals that they are almost identical with a conventional velocity analysis.The algorithm is based on the concept that a normal-moveout pattern of time shifts may be accurately applied over a considerable span of zero-offset times to a CDP trace gather, rather than only at a single zero-offset time. With each application however, the rms-velocity associated with the moveout pattern is different. A specific relationship which associates an rms-velocity with a zero-offset time for the same moveout time delay pattern is approximated by numerical simulation studies on a digital computer.

Notes: The finite element method, here viewed as a special case of the Galerkin projective method, is applied to the modelling of magnetotelluric problems, and its adaptation to geological profiles is outlined. A novel method for obtaining surface field values, involving matrix representation of the normal derivative operator, is presented in detail.Results obtained by this method are compared with well-known infinite series solutions for the vertical fault and the outcropping dyke. Two profiles containing sulphide zones are also modelled, the results being compared with field data; satisfactory agreement is obtained.

Notes: Various authors have recently studied the Walsh and Hadamard transforms for functions with one or two variables. Like the Fourier transform they may be obtained by using the Cooley-Tukey algorithm. They can be calculated all the faster because Walsh functions have only values equal to + 1 or − 1. The complex number multiplications in the Fourier transform are replaced by additions and subtractions of real numbers.These transforms are more general than that of Fourier because the concept of sequency, i.e. the mean number of zeros per time unit, is substituted for that of frequency.These transforms already have numerous applications to signal transmission. We describe some in the field of seismic prospecting, in particular the compression of the information contained in seismic traces.

Notes: A detailed analysis of the electromagnetic anomaly due to an elongated mineralized zone of low resistivity in presence of a fixed-transmitter (a long cable carrying an, alternating current IeIωt) is presented where (i) the host rock is not highly resistive and might contain some disseminated mineralization, (ii) the ore-body has an inhomogeneous conductivity, (iii) there is a contrast in magnetic properties of the ore-body and the host rocks, and (iv) source is close to the conducting system i.e., anomalous zone is lying shallow.The numerical computation of the generalized analytical expressions have been made. The selective screening behaviour of the cover is delineated and suitable frequency ranges for maximum detectability of a covered conductor have been obtained. A paradoxical decrease of the in-phase component of the response function with the increase of the core conductivity has been found for large conductivities of the cover. For moderate values of response parameter the inhomogeneity in the conductivity is found to significantly affect the induced field. The present analysis of the various features of the secondary field will aid to the existing interpretation sensitivity of the induction prospecting data for porphyry conducting ore deposits with zonal wall-rock alteration and sulphide distributions.

Notes: This report describes the construction and properties of a Ge(Li) bore hole probe for in situ gamma ray spectrometry. The probe has the form of a cylinder 8.9 cm in diameter and 160 cm long. It contains a 22 cm3 coaxial Ge(Li) detector, cooled by liquid nitrogen. Holding time of the assembled probe is in excess of ten hours, energy resolution is 2.3 keV for 122 keV photons and 2.6 keV for 1333 keV photons.The probe has been used in field experiments in the Ranstad uranium mine and for the study of neutron activation spectra obtained in bore hole models. Results from these experiments are presented and discussed. It is stressed that the Ge(Li) probe gives promise to become a powerful tool in bore hole activation analysis.

Notes: The indirect method of interpreting resistivity data is capable only of limited success because of the difficulty of calculating type curves for complex structures. Consequently a need arises for a direct method of interpretation for complex generalized structures.Such a technique for the direct interpretation of apparent resistivity data, obtained by electrical soundings carried out over two-dimensional structures, has been developed from an examination of the Hankel transform of such data. The method is based on the observation that the asymptotic expansion of the Hankel transform is critically dependent upon the minimum distance between the measuring device and the surface of discontinuity in resistivity.The variation in the across strike direction may be mapped by a sequence of depth soundings made parallel to each other and separated in the across strike direction.

Notes: The process of data smoothing by least squares operators with real and symmetric coefficients is shown to be equivalent to linear digital filtering. Study of the frequency response functions of the least squares operators shows that these operators are essentially low-pass filters. The frequency response functions are characterized by the presence of a main lobe centered about the zero frequency together with a number of minor lobes of smaller amplitudes at higher frequencies. The predominant band of frequencies which is passed by a least squares operator occurs within the main lobe of the frequency response function. The filtering characteristics of a least squares operator are dependent on the degree of the polynomial and the number of data points. Generally speaking, the higher the degree of the polynomial the broader the pass band of the filter will be. But the effect of the number of data points on the pass band is exactly opposite.Several figures showing the plots of the half-power frequency, the quarter-power frequency, and the half-width of the main lobe for a large number of least squares operators with different combinations of the degree of the polynomial and the number of data points are presented in this paper. These figures will be helpful to those people who wish to make a selection of least squares operators suitable for their purposes.

Notes: The interpretation of a gravity anomaly in terms of the shape of the anomalous body is a non-linear problem and may, therefore, be carried out using non-linear optimisation techniques. The formulation is extended to include cases where the density contrast and the regional background are also unknown. For a given model the objective function is provided by the discrepancy between the observed anomaly and the calculated anomaly due to the model. Given an initial model, the optimisation procedure searches for a minimum of the objective function by an iterative adjustment of the parameters. A number of suitable objective functions is given. The behaviour of these functions in the parameter hyperspace is quite complex. Accordingly, direct search methods should be employed at the early stages of the search, changing to gradient methods at later stages. The use of constraints is also necessary to ensure the geological feasibility of the model. The required computer time may be largely reduced by careful programming. Two examples of interpretation by optimisation methods are given.