CA2100248C - Method and apparatus for surveying terrain resistivity - Google Patents
Method and apparatus for surveying terrain resistivityInfo
- Publication number
- CA2100248C CA2100248C CA002100248A CA2100248A CA2100248C CA 2100248 C CA2100248 C CA 2100248C CA 002100248 A CA002100248 A CA 002100248A CA 2100248 A CA2100248 A CA 2100248A CA 2100248 C CA2100248 C CA 2100248C
- Authority
- CA
- Canada
- Prior art keywords
- vlf
- vertical
- magnetic field
- survey line
- terrain
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Electromagnetism (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
Terrain resistivity is mapped by collecting data as to the vertical magnetic field ratio along a survey line by VLF survey techniques, and precision integrating rather than differentiating the data obtained. The received VLF signal is demodulated and low-pass filtered in a manner such as to preserve the DC component of the demodulated signal, and is then digitized and spatially integrated to obtain pseudo-resistivity readings along the survey line.
Description
2 1 0 0 2 ~ ~ PCT/CA91/00401 METHOD AND APPARATUS FOR ~uKv~Yl~G TERRAIN RESI~llvl~r 8ACKGROUND OF THE lNV :N l'lON
FIF~n OF T~ lNV~:NllON
This invention relates to geophysical ~ ying, and more particularly to te~h~iques for surveying terrain resistivity.
REVIEW OF THE ART
Various techniques have been utilized for the dete~tion of subsurface geological anomalies, characterised by different resistivity from ~uL~o~nding terrain. One such techni~ue relies upon the detection of s~on~ry signals produced as a result of ex~u~le of terrain being mapped to the coherent electromagnetic field produced by a distant VLF radio station. Some what analogous techniques utilize non-coherent sources such as naturally occurring magneto-telluric and audio frequency magneto-telluric electro-magnetic fields.
Available VLF electro-magnetic fields have been used for over twenty years to map subsurface geology. In the one mode of implementation the horizontal (primary) magnetic field component is used as a phase and amplitude ~;UBSTITUT_ SHEET
2~ 0 ~ 248 _ - 2 -reference against which to measure the normalized in-phase and quadrature phase components of the secondary vertical magnetic field. In some cases the tilt angle (or tangent of the tilt angle, usually expressed as a percentage) and elipticity of the polarization ellipse is measured instead, but in the usual case where the anomalous or secondary VLF
magnetic field components are relatively small compared with the primary field this measurement technique also yields quantities which are simply related to the in-phase and quadrature phase components of the vertical magnetic field.
In the early days of the VLF technique, survey interpretation was based on the assumption that, even in an electrically conductive earth, the time-varying primary magnetic field induced vortex or eddy currents in vertical conductors. These vortex currents generated the secondary magnetic field which was used to locate and characterise the conductor, using well established but flawed principles governing the ratio and phase of the secondary field components as functions of target size, depth and conductivity. Since only targets with appreciable vertical extent would be excited by the horizontal primary magnetic field the expected (and indeed measured) anomalies often took the form of a lateral S-curve with a zero-crossing.
The location of the target was under the zero-crossing or, in the case where the response from nearby targets was also strong enough to contribute to the total response, the location was under the point of maximum inflection of the central limb of the curve.
Such multiple-target responses occur frequently and the resultant interference between the anomalies sometimes makes survey data interpretation quite difficult.
In order to rectify this situation Fraser ("Contouring of VLF-EM data", Geophysics 34, page 958-967) proposed WO93~10473 210 ~ 2 ~ ~ PCT/CA91/00401 ~_ - 3 -applying a simple numerical filter to the survey data which would simultaneously (1) numerically differentiate the data horizontally so as to transform the cross-overs into peaks, which could then be contoured.
(2) high-pass filter the data so as to attenuate the response from slowly varying anomalies and relatively enhance the response from localized anomalies.
FIF~n OF T~ lNV~:NllON
This invention relates to geophysical ~ ying, and more particularly to te~h~iques for surveying terrain resistivity.
REVIEW OF THE ART
Various techniques have been utilized for the dete~tion of subsurface geological anomalies, characterised by different resistivity from ~uL~o~nding terrain. One such techni~ue relies upon the detection of s~on~ry signals produced as a result of ex~u~le of terrain being mapped to the coherent electromagnetic field produced by a distant VLF radio station. Some what analogous techniques utilize non-coherent sources such as naturally occurring magneto-telluric and audio frequency magneto-telluric electro-magnetic fields.
Available VLF electro-magnetic fields have been used for over twenty years to map subsurface geology. In the one mode of implementation the horizontal (primary) magnetic field component is used as a phase and amplitude ~;UBSTITUT_ SHEET
2~ 0 ~ 248 _ - 2 -reference against which to measure the normalized in-phase and quadrature phase components of the secondary vertical magnetic field. In some cases the tilt angle (or tangent of the tilt angle, usually expressed as a percentage) and elipticity of the polarization ellipse is measured instead, but in the usual case where the anomalous or secondary VLF
magnetic field components are relatively small compared with the primary field this measurement technique also yields quantities which are simply related to the in-phase and quadrature phase components of the vertical magnetic field.
In the early days of the VLF technique, survey interpretation was based on the assumption that, even in an electrically conductive earth, the time-varying primary magnetic field induced vortex or eddy currents in vertical conductors. These vortex currents generated the secondary magnetic field which was used to locate and characterise the conductor, using well established but flawed principles governing the ratio and phase of the secondary field components as functions of target size, depth and conductivity. Since only targets with appreciable vertical extent would be excited by the horizontal primary magnetic field the expected (and indeed measured) anomalies often took the form of a lateral S-curve with a zero-crossing.
The location of the target was under the zero-crossing or, in the case where the response from nearby targets was also strong enough to contribute to the total response, the location was under the point of maximum inflection of the central limb of the curve.
Such multiple-target responses occur frequently and the resultant interference between the anomalies sometimes makes survey data interpretation quite difficult.
In order to rectify this situation Fraser ("Contouring of VLF-EM data", Geophysics 34, page 958-967) proposed WO93~10473 210 ~ 2 ~ ~ PCT/CA91/00401 ~_ - 3 -applying a simple numerical filter to the survey data which would simultaneously (1) numerically differentiate the data horizontally so as to transform the cross-overs into peaks, which could then be contoured.
(2) high-pass filter the data so as to attenuate the response from slowly varying anomalies and relatively enhance the response from localized anomalies.
(3) not emphasize random geological noise.
(4) be simple to apply for untrained field personnel.
The resultant "Fraser filter" has been widely used in VLF data processing. Examples of its application will be considered hereinafter.
More recently Karous and Hjelt ("Linear filtering of VLF dip-angle measurements", Geophysical Prospecting 31 (1983), pages 782-794) proposed a similar numerical filter, the output of which is interpreted in terms of the equivalent current density at a specified depth that would cause the measured VLF magnetic field strength. This filter also gives contourable anomaly peaks rather than cross-overs.
Finally Lankston and Hecker ("Enchancing VLF-EM
data through application of frequency domain operators", 2S Proc. 2nd National Outdoor Action Conference on Aquifer Restoration, Groundwater Monitoring and Geophysical Methods, Las Vegas (1988), pages 655-673) suggest use of the Hilbert transform to improve the geological signal-to-noise ratio over that achieved by the Fraser filter, but at the expense of greater mathematical complexity.
From the above it will be noted that use of numerical data processing (i.e. filtering) to enhance VLF
data and simplify interpretation is a common practice.
~ 4 ~ 2~Q~48 SUMMARY OF THE INVENTION
I have found as a result of studies on the role of vortex currents in VLF responses that, contrary to general belief, it is the galvanic current component (also known as current gathering) that dominates the response from typical anomalies generated by VLF radio waves, at all but the highest background resistivities (>5,000 nm).
These currents, which are generated by the horizontal electric field in the earth, flow in a horizontal direction radial to the transmitter. The distribution of this current flow, as it would occur in a uniform earth, is distorted by the presence of conductive (and resistive) anomalies, in turn causing an anomalous secondary magnetic field. Interestingly, the secondary magnetic field due to galvanic current response from a vertical conductive plate aligned with its long axis pointing to the source is very much like that from the notional vortex current response already discussed above, since both vortex and galvanic responses arise principally from the enhanced current flow in the top of the conductor. However the response from galvanic current flow behaves very differently (compared with vortex current flow) as a function of target size, depth, conductivity contrast and dip, among other factors.
Indeed a thin horizontal plate conductor will generate virtually no vortex current flow anomaly (being in null coupling with the horizontal primary magnetic field) whereas it can and does generate a large galvanic current flow anomaly.
Furthermore other geological features, such as a vertical contact, also generate a distinctive and measurable anomaly as a result of galvanic current flow.
The reason for this is that the variation of galvanic current flow with depth is determined entirely by skin depth considerations. At high values of terrain WO93/10473 2 1 Q 0 2 ~ 8 PCT/CA91/0~01 _ - 5 -resistivity, the skin depth (defined as ~ = (2/~aw)'h where is a constant, o is the ground conductivity and w is the angular frequency of the source) is large, and current flows to relatively great depth, whereas at small terrain resistivity the skin depth is small and current flows only to a shallow depth. The galvanic current behaviour in the vicinity of a vertical contact between zones of different resistivity and the resultant magnetic field anomaly produce a response curve which, although there is no cross-over, is somewhat similar to part of the anomaly producedby a vertical conductive plate as discussed above.
Furthermore the sharp anomaly decrease on the conductive side of the anomaly means that the differentiating action of the various filters described above creates an lS asymmetrical but large anomaly at the location of the contact which can easily be mistaken for the anomaly from a vertical plate. Despite the fact that the filtered anomaly from the contact-has a different shape from-that for the vertical plate (or other anomaly) it turns out in actual practice that the data from VLF surveys contain so many anomalies, of so many slightly different shapes, that it is often impossible, in either raw or filtered data sets, to further identify the source of the anomaly.
This has been a serious shortcoming of the VLF
technique. What is required is a measured or derived quantity that is indicative of those variations of the resistivity of the ground which are causing the observed anomalies. Variations in the ground resistivity can presently be measured using VLF radio waves by inserting two electrodes in the ground, aligned so that the array axis is parallel to the direction to the transmitter, and measuring the VLF potential difference between these two electrodes as described by Collett and Becker in Canadian Patent No. 795,919. However this procedure is very time consuming and detracts from the most attractive feature of the measurement of the magnetic field components, which is _ - 6 -that no electrical contact with the ground is required and thus these measurements can be carried out very quickly and inexpensively.
I have adopted an alternative approach based upon my determination of the nature of the galvanic current flow which is in actuality responsible for virtually all VLF
anomalies. In this case the effective depth of the current flow i.e. skin depth, will more or less continually reflect local resistivity, and a lateral change in resistivity will result in a change in depth of the current flow to accommodate the resulting change in skin depth. Any change in depth of the current flow in a direction perpendicular to that flow will result in a vertical component of the magnetic field, the amplitude of which will vary according to the rate of change of terrain resistivity with distance.
According to the invention, a method of surveying the resistivity characteristics of terrain comprises measuring a series of values indicative of local relative magnitudes of a vertical magnetic field component as compared to a horizontal magnetic field component due to electromagnetic radiation from a distant stationary VLF
source whilst moving along a survey line generally perpendicular to the direction of the source, so as to provide data capable of providing a record of changes in said values along said survey line, and spatially integrating the data so obtained with respect to displacement along the survey line to provide a record indicative of the resistivity characteristics of the terrain along said survey line.
Values which may be measured include the ratio of the amplitude of the in-phase component of the vertical magnetic field to the amplitude of the horizontal magnetic field occasioned by radiation from the VLF source, the ratio of the amplitudes of the total vertical and ~',' WO93/10473 2 ~ ~ 0 2 ~ g PCT/CA91/00401 _ - 7 -horizontal fields occasioned by radiation from the VLF
source, and the tilt angle or tangent of the tilt angle of the ellipse of polarization of the magnetic field occasioned by radiation from the VLF source.
The invention extends to apparatus for carrying out the above method.
Further features of the invention will be apparent from the following description with reference to the accompanying drawings.
SHORT DESCRIPTION OF THE DRAWINGS
Figures la, lb are diagrams illustrating traditional theory as to the generation of electric and magnetic fields by a distant VLF radio transmitter, along with flow of vortex currents within a radially oriented vertical plate conductor;
Figure 2a is fragmentary section through an area of terrain further illustrating traditional theory as to currents and fields in the vicinity of the conductor shown in Figure l, Figure 2b illustrating the resulting vertical magnetic field anomaly, and Figure 2c illustrating how this anomaly may be distorted by an adjacent anomaly;
Figure 3 is a diagram corresponding to Figure lb, but illustrating the galvanic current flow which is in practice dominant in the vicinity of the buried conductive plate;
Figure 4a is a diagram corresponding to Figure 3, in which the conductive plate is replaced by a contact area between zones of resistive and conductive material, and Figure 4b illustrates the vertical magnetic field anomaly in the vicinity of the contact;
- 8 - ~ 4 8 -Figure 5a is a vertical section through terrain containing various common subterranean features, Figure 5b illustrates the profile of the resulting vertical magnetic field anomaly, Figure 5c illustrates the profile obtained by difference filtering the profile of Figure 5b, and Figure 5d is a profile illustrating the effect of continued spatial integration of the profile of Figure 5b from left to right;
Figure 6 is a block diagram illustrating apparatus for generating data forming a profile such as that shown in Figure 5d;
Figures 7a and 7b are families of curves illustrating the response obtained after, spatially integrating the vertical magnetic field ratio along survey lines over an idealised vertical dike and an idealised vertical contact respectively;
Figures 8a - 8e illustrates the manipulation of actual VLF field survey data, Figure 8a plotting collected VLF vertical magnetic field data, Figure 8b showing the data horizontally differenced, Figure 8c plotting a corresponding resistivity profile prepared from field collected resistivity data utilizing spaced electrodes, Figure 8d illustrating the actual geology along the survey line, and Figure 8e plotting the data of Figure 8a spatially integrated along the survey line to simulate the operation of the method and apparatus of the invention;
Figures 9 and lOa - lOd illustrate a further example based on actual field survey data, Figure 9 showing a geological section illustrative of the type of terrain surveyed, Figure lOa being a plot of collected VLF vertical magnetic field data, Figure lOb a plot of the data after Fraser filtering, Figure lOc a plot of field surveyed conductivity data, and Figure lOd a plot of spatially Gi WO93/10473 210 0 2 i~ 8 PCT/CA91/0~01 _ g integrated VLF vertical magnetic field data along the survey line; and Figure 11 is a flow diagram illustrating one possible implementation of the invention in microprocessor controlled VLF survey equipment.
DESCRIPTION OF TH~ PREFERRED EMBODIMENTS
Referring to Figure la, the primary VLF field components generated at a point P on a conductive half-space by a distant VLF transmitter T are vertical and radial electric field components Ep and E" and a horizontal magnetic field component ~. Assuming the presence at P of - a subsurface radially and vertically extending conductive plate C (see Figures lb and 2a) the component ~ will induce a vortex current flow V in the plate which in turn will generate a secondary magnetic field ~ which will have vertical components on either-side of the plate, producing a vertical magnetic field anomaly as shown in Figure 2b, which plots the vertical magnetic field component ~ along the same section line as Figure 2a. Figure 2c illustrates how the curve of Figure 2b may be distorted if the vertical field components generated by plate C overlap with those (shown in broken lines) produced as a result of another nearby conductive body.
In practice I have found that, except under conditions where large steeply dipping and highly conductive bodies are found in terrain of generally very high resistivity, the above discussed generation of vertical field components due to vortex current flow is much less significant than the vertical field components produced by distortions of the galvanic current flow through the terrain resulting from the radial electric field component E" as illustrated in Figure 3, in which the presence of the conductor C causes a local reduction in skin depth and a 'bunching' effect as shown. As it ~ ~ ~ Q ~ ~ 8 ~
_, -- 10 --happens, this bunching in the example shown produces a vertical magnetic field anomaly very similar to that of Figure 2b.
Considering now the feature shown in Figure 4a, in which the plate C is replaced by a contact zone D
between relatively resistive material (on the left) and relatively conductive material (on the right), the galvanic current flow to the right of the contact will be concentrated nearer to the surface than that to the left of the contact. The distortion of the current flow near the contact zone produces a vertical field anomaly as shown in Figure 4b, which in real life conditions could readily be confused with a distorted anomaly as shown in Figure 2c.
Figure 5a shows schematically a section through an imaginary geological structure of somewhat greater complexity, in which, facing towards a distant VLF
transmitter, a highly conductive vertical plate like target T is flanked by zones of moderately conductive material S
which in turn contact zones of moderately resistive material R. The symbols ~ represent galvanic currents passing into the page, whilst the arrows show the direction of the resulting magnetic field. Figure 5b plots the intensity of the vertical magnetic field component ~, and Figure 5c shows the effect of differencing this field component, the plot being inverted so that the response over the wanted conductor is upward. Neither plot provides any clear picture of the underlying structure.
The effect of continually spatially integrating the vertical magnetic field profile from left to right is shown schematically in Figure 5d. The process of integration makes the presence of the two contacts, the presence of the enclosed block of conductive material, and the presence of the conductive target much more obvious;
indeed the profile of Figure Sd can be thought of as a WO93/10473 21 ~ 0 2 ~18 PCT/CA91/00401 profile of "pseudo-conductivity" and as such is a useful supplement to the raw and filtered VLF data. The profile is not, of course, an actual conductivity profile since it is generated by the magnetic field which is generated, in turn, by currents whose varying depth is defined by the skin depth, which is a function of the square root of resistivity. However, it is much more useful than the difference filtered profile, shown in Figure 5c, which yields no information about the background conductivity against which the anomalies are being detected and also makes correct identification of the source of any specific anomaly quite difficult.
Figure 6 illustrates apparatus for generating such pseudo-conductivity profiles. A generally conventional VLF receiver 2 is configured and tuned to receive at least the vertical component of anomalous secondary magnetic fields locally generated as a result of transmissions from a distant VLF transmitter located approximately perpendicular to a survey line over which the receiver is conveyed. In a preferred embodiment the relative magnitude of the in-phase vertical and horizontal components of these fields is required, and the quadrature components can be ignored, thus simplifying the receiver, which may be constructed in accordance with established practice in the VLF survey art. Reference may if needed be had to Canadian Patent No. 937,286 (Ronka) for further description of such a receiver, although in the present case simplification may be, as stated, possible where only the real component of the vertical field is to be detected.
In practice, as discussed further below, the signal which is wanted is a value indicative of the relative magnitude of the vertical magnetic field with respect to the horizontal magnetic field, in order to compensate for changes in the horizontal field.
~- - 12 - 2 ~ 9 As in the Ronka patent, the vertical and horizontal signals are received by separate channels and separately detected by a detector means 3, which must be constructed so as to preserve DC components of the demodulated signals, which are low pass filtered by filter means 4, conveniently effected using so-called integrator circuits which have a sufficiently long time constant to smooth out short term fluctuations in the output (such as those caused by modulation of the source transmitter) without obscuring terrain induced features of the signal profile which it is desired to examine. The filtered signals are combined by determining the magnitude of the vertical component as a fraction of the horizontal component, and the combined signal is passed to an analog-to-digital converter 6, having a sampling rate controlled by a timing generator producing signals related to increments of distance covered by a carrier 1, such as an aircraft, or portable instrument case, in which the receiver is carried along a survey line. If the carrier moves at a constant speed, the distance increments will be constant if the period between samples is constant. The sampling rate should be sufficient to permit the signal profile to be adequately resolved.
Rather than using an analog low pass filter, the sampling rate may be increased, and digital low pass filtering applied to produce a desired spectral content by summing groups of successive samples to provide samples with the required band-pass characteristics.
The resulting digital data is stored in a suitable memory device 8 for further processing. The data so stored may be subjected to further digital processing using similar algorithms to those discussed above in relation to prior art VLF techniques, to obtain digitized signal profiles which can be plotted or otherwise displayed by a display unit 10 in known manner to permit .. ..
WO93/10473 2 i ~ ~ 2 ~ 8 PCT/CA91/0~01 ~- - 13 -interpretation. Additionally or instead, successive samples are added in an accumulator 12 to obtain the left-to-right integration discussed above, the accumulator total after each addition being again stored in memory device 8 to build up a pseudo-conductivity profile which can be suitable scaled and plotted or displayed by display unit 10 .
The data required for operation of the invention are values indicative of the relative magnitude of a vertical field component ~ as compared to the horizontal magnetic field component ~ produced by the VIF excitation.
Examples of such values are the ratio of the in-phase vertical field component ~(I) to ~, the ratio of the total vertical field component ~ to ~, or the tilt angle or tangent of the tilt angle of the polarization ellipse of the field produced. Those values involving only in-phase components are preferred, particularly the ratio ~(I)/~.
Assuming that data representing any of the above values is measured at a series of stations spaced along the survey line at an interval of ~x metres, and that the value measured at station i in the series where i>0 is Vi, all ratios being expressed as decimal fractions, then the result of integration is X(i)=X(i-l)~V(i).~x.R
where X(-l)=0, X(i) is the spatially integrated value after station i, and R is an arbitrary constant or scale factor, and successive values for X(i) are plotted at plot points x = (i+0.5). x. If the units of ~x are metres, then a convenient value for R is unity, in which case the data can be expressed approximately in terms of "pseudoconductivity"
expressed in milliSiemens/metre (mS/m). The choice of R=l (where x is in metres) will usually result in the plotted values being of a similar order of magnitude to the actual conductivity. A notional example is shown in the following table:
~_ - 14 - ~ 4 8 -i xi=i.~x Vl Vi.~x Xi Plot Point(i+~ x 2 200 0.1 10 10 250 3 300 0.2 20 30 350 4 400 0.3 30 60 450 500 o.o 0 60 550 6 600 -0.3 -30 30 650 7 700 -0.2 -20 10 750 8 800 -0.1 -10 0 850 9 900 0.O O o 950 It should be appreciated that the apparatus used may be largely conventional, provided that it has the capability of receiving and separately demodulating, preferably synchronously, the vertical and horizontal field components. Synchronous demodulation is preferred in order to recover the vertical component which is in-phase with the horizontal component. The apparatus must of course be modified to include the integration function provided by accumulator 12. Thus the system must be able to preserve accurately the DC component of the amplitude profile of the vertical magnetic field component, relative to the amplitude of the horizontal component, and the integration of the signal must also be achieved accurately and without drift. For this reason, it is preferred that the detected signal representing the relative vertical magnetic field ratio be digitized prior to integration, which is then performed numerically. It may be advantageous to utilize digital techniques to demodulate the received VLF signal components and to perform any necessary low pass filtering required to exclude unwanted short term fluctuations in the signal. The demodulation and temporal integration technique disclosed in the Figure 5 embodiment of the Ronka ~093/10473 PCT/CA9l/0040l - 15 2 1 ~ 3 2 4 8 patent is not suitable, not only because the integration is temporal rather than spatial, but also because it is difficult to maintain sufficient DC stability for the purposes of the present invention in an analog integrator S such as that disclosed by Ronka. Furthermore, the Ronka arrangement deliberately utilizes negative feedback to modify the DC level at the output of the integrator and thus in effect greatly reduces the direct current gain of the integrating amplifier. The integrating amplifier is thus caused to act as a bandpass filter which, whilst appropriate for recording data as to the amplitude of the vertical component, renders accurate integration impossible and essentially negates the function of the amplifier as anything other than a filter.
Digital processing is advantageous since it facilitates accurate maintenance of DC levels.
Furthermore, if signal processing is performed under microprocessor control, the required accumulation or precision integration step may be performed by a simple addition to the control program which program may include routines for performing optimal filtering techniques as well as the method of my invention. A flow diagram of an exemplary routine forms Figure 11, and assumes that the samples from the digital to analog computer have been placed in a data storage file. The routine is self-explanatory, and results in the generation of a second file containing samples which represent a pseudo-conductivity plot along the same survey, in the manner illustrated in Table 1.
Before proceeding to some theoretical and field examples further illustrating the invention, two points should be made. Since integration is achieved by adding the values of data taken at discreet measurement points which are separated by a fixed interval x, the value of the integral at any distance x will be a function of 1/ x, - 16 ~ 21~024~
which is why the correction factor ~x is four~ in the equation presented above. Further, the true zer~ ~evel of the integrated data is not known. Whilst - ~s s a disadvantage, the most significant information is ~e ~ted by relative changes in the integrated signal alor 1 rvey line.
Figures 7a and 7b present results o lined from integrating the theoretical response from t- important geological targets, viz a vertical sheet-like conductor (dike) (Figure 7A) and a vertical contact (Figure 7b). The curves were calculated from theoretical data giving the normal vertical in-phase VLF magnetic field resulting from two dimensional inhomogeneities oriented with their long axis parallel to the direction to the VLF transmitter.
Consider first a vertical dike, of thickness t, and conductance S = ot, buried in a half space of resistivity P2. The integrated response becomes a simple bell-shaped curve centered over the target, with maximum value a function of S and P2, and with width proportional to 20 P2. For the extremely resistive half-space shown here (P2 = 8000 nm) the anomaly is quite wide (full width to half-maximum of about 300m), but in the more typical case of P2-1000 nm the full width to half maximum value (FWHM) becomes about lO0 m. Indeed it can be shown that, in general, the FWHM is approximately equal to the skin depth in the half-space. Furthermore the high value of P2 also yields a high peak value of the integrated anomaly; reduction of the half space resistivity results in a substantial reduction of anomaly amplitude as well as anomaly width.
Similar results but for a contact between zones of resistivity p~ and P2 respectively are shown in Figure 7b. The apparent conductivity increases as one traverses from the resistive to the conductive side of the contact by an amount which is a function of ratio Pl/P2 and which WO93/10473 2 ~ O ~ 2 4 ~ PCT/CA91/00401 ~_ - 17 -saturates for large values of this ratio. Similar curves will apply for other values of Pl/P2- The usual principles of electrical similitude apply. For example, the ratio Pl/P2 for curve B is lO. If now p~ = lO00 Qm and P2 = lO0 nm 5 (Pl/P2 = lO again), the new curve would look like curve (B) except that both the vertical and horizontal scales would be divided by (lO)~ = 3.2. The anomaly will thus be smaller and more closely confined to the contact. Such an anomaly is shown by curve (C).
The responses from these two anomalies, viz a vertical conductive dike and a vertical contact, are useful in interpreting field survey results.
(Prakla Seismos - location not known) The first case history is selected because it shows very well the relationship between VLF vertical magnetic field measurements, the integrated "pseudoconductivity" data calculated as described above, and terrain resistivity. Figure 8a shows the vertical magnetic field component measured at lO meter intervals along an 850 m survey line over the geology illustrated in Figure 8d. Figure 8b is the simple horizontal difference of the in-phase data extracted from the data shown in Figure 8a. The reversal of polarity usually applied to the 2S differenced data has not been carried out here since we wish to compare this panel with the resistivity data of Figure 8c which is the result of Schlumberger resistivity profiling with half-current electrode spacing of 205 m and half-voltage electrode spacing of 5 m. The generally good correlation between lows in the differenced in-phase vertical magnetic field data and the Schlumberger apparent resistivity lows is immediately obvious.
~~ - 18 -The results of spatially integrating the profile of Figure 8a are shown in Figure 8e. It is immediately apparent, by comparison with Figure 7b, that there is a contact at the llOm station, with resistive material to the left, conductive material to the right and about a 10/1 conductivity contact. The presence of this contact is completely missed by the raw and filtered data. In essence the data of Figure 8b is that of Figure 8e, doubly horizontally differentiated, which is exactly what it looks like. All the peaks of Figure 8e are accentuated in Figure 8b, but the all-important slowly varying component, which is related to the terrain conductivity, has been completely suppressed.
Comparing the resistivity curve shown in Figure 8c with Figure 8e, the low resistivity zones at stations 130, 220, 500 - 600 and 700 correspond extremely well to the pseudoconductivity highs at these same locations, and also the resistivity highs at stations 330, 460, 620 and 740 correspond well with the pseudoconductivity lows.
Furthermore the pseudoconductivity trace of Figure 8e depicts very clearly the fact that there is a contact with an adjacent conductor at the llOm station, a fact which no interpreter would have picked from the raw or differenced VLF data.
FIELD EXAMPLE #2 (Niagara, N.Y.) A typical geological section for the area surveyed in the second case history is shown in Figure 9.
Highly resistive dolomite bedrock 20 has at some time in the past been fractured, and the fractured zones 28 have subsequently been differentially weathered (the fracturing allows water to enter the rock and expand on freezing) to form bedrock topographic lows or troughs 22. More recently (post glacially) the bedrock surface has been covered by a lake with the result that electrically conductive ~, W093/10473 21 C ~ 2 4 ~ PCT/CA91/0~1 _ -- 19 --lacustrine (lake) clays 24 have been deposited on the bedrock to an average thickness of 3 - 8 meters. Thus the geoelectric section consists of an electrically conductive layer of variable thickness overlying an electrically resistive bedrock; fracture zones will often exist in the bedrock at the base of the troughs and these fractures may be electrically conductive as a result of infilling with water-, conductive clays or mineralization. A similar type of bedrock structure is particularly important in many countries for locating adequate sources of groundwater, in which case the interpreter wants to pick out those fracture zones which lie in the deepest troughs so that there will be maximum water available for pumping (see the line 26 indicating the level of the water table.) lS A VLF survey line was run over a geological section of this general type; survey results are shown in Figures lOa - lOd. The ~v~y line was 430 m long, and the raw VLF data shown in Figure lOa indicates three conductors located at A, B and C respectively. These conductors were picked (as is customary when anomalies are so close together as to cause each other to be vertically off-set) by choosing the points of maximum inflection. Figure lOb, which shows the Fraser filtered data, confirms this choice of conductors.
Figure lOc shows results of an additional conductivity survey over the area using a Geonics (trade-mark) EM34-3 Ground Conductivity Meter, constructed in accordance with my U.S. Patent No. 4,070,612. This survey profile shows a conductivity high located between stations 100 and 200. Since the clays are more conductive, this conductivity high corresponds to a thickening of the clay layer, i.e. to a depression or trough in the bedrock since the ground surface is virtually flat. Knowing that conductors A and B are located in the trough would make them the first priority drill targets for groundwater ~~ - 20 - ~ 4 8 -~
exploration. Target C, on the other hand, located on the flank of the trough, would be lower priority.
Unfortunately, neither the raw or Fraser filtered VLF data gives us the essential conductivity information to allow rational selection of any of the three conductors as having more groundwater potential than the others. On the other hand, the horizontally integrated VLF data, or pseudoconductivity, shown in the fourth panel, clearly shows that the conductors at stations 125 and 200 are located in a conductivity high or bedrock trough and that the conductor at station 260 is on the flank, completely confirming the results of the conductivity data, but without requiring the additional survey with a ground conductivity meter.
Whilst the invention has been described primarily with reference to the use of distant VLF radio transmitters as a VLF source, other naturally occurring or artificial sources of electromagnetic fields with frequency components in a similar or lower range may be exploited as is known in the art. Since these fields are typically non-coherent, the absolute magnitude of the vertical field component must be determined in such a case. For the purpose of interpreting the appended claims, such a source should be considered to be within the term 'VLF source'. Moreover, instead of a receiver that resolves horizontal and vertical components of the magnetic field, a receiver that directly determines a function of the tilt angle such as the angle itself or its tangent, of the polarization ellipse of the field may be utilized; the tangent of this angle is proportional to the magnitude of the vertical component for small values of the latter.
V
The resultant "Fraser filter" has been widely used in VLF data processing. Examples of its application will be considered hereinafter.
More recently Karous and Hjelt ("Linear filtering of VLF dip-angle measurements", Geophysical Prospecting 31 (1983), pages 782-794) proposed a similar numerical filter, the output of which is interpreted in terms of the equivalent current density at a specified depth that would cause the measured VLF magnetic field strength. This filter also gives contourable anomaly peaks rather than cross-overs.
Finally Lankston and Hecker ("Enchancing VLF-EM
data through application of frequency domain operators", 2S Proc. 2nd National Outdoor Action Conference on Aquifer Restoration, Groundwater Monitoring and Geophysical Methods, Las Vegas (1988), pages 655-673) suggest use of the Hilbert transform to improve the geological signal-to-noise ratio over that achieved by the Fraser filter, but at the expense of greater mathematical complexity.
From the above it will be noted that use of numerical data processing (i.e. filtering) to enhance VLF
data and simplify interpretation is a common practice.
~ 4 ~ 2~Q~48 SUMMARY OF THE INVENTION
I have found as a result of studies on the role of vortex currents in VLF responses that, contrary to general belief, it is the galvanic current component (also known as current gathering) that dominates the response from typical anomalies generated by VLF radio waves, at all but the highest background resistivities (>5,000 nm).
These currents, which are generated by the horizontal electric field in the earth, flow in a horizontal direction radial to the transmitter. The distribution of this current flow, as it would occur in a uniform earth, is distorted by the presence of conductive (and resistive) anomalies, in turn causing an anomalous secondary magnetic field. Interestingly, the secondary magnetic field due to galvanic current response from a vertical conductive plate aligned with its long axis pointing to the source is very much like that from the notional vortex current response already discussed above, since both vortex and galvanic responses arise principally from the enhanced current flow in the top of the conductor. However the response from galvanic current flow behaves very differently (compared with vortex current flow) as a function of target size, depth, conductivity contrast and dip, among other factors.
Indeed a thin horizontal plate conductor will generate virtually no vortex current flow anomaly (being in null coupling with the horizontal primary magnetic field) whereas it can and does generate a large galvanic current flow anomaly.
Furthermore other geological features, such as a vertical contact, also generate a distinctive and measurable anomaly as a result of galvanic current flow.
The reason for this is that the variation of galvanic current flow with depth is determined entirely by skin depth considerations. At high values of terrain WO93/10473 2 1 Q 0 2 ~ 8 PCT/CA91/0~01 _ - 5 -resistivity, the skin depth (defined as ~ = (2/~aw)'h where is a constant, o is the ground conductivity and w is the angular frequency of the source) is large, and current flows to relatively great depth, whereas at small terrain resistivity the skin depth is small and current flows only to a shallow depth. The galvanic current behaviour in the vicinity of a vertical contact between zones of different resistivity and the resultant magnetic field anomaly produce a response curve which, although there is no cross-over, is somewhat similar to part of the anomaly producedby a vertical conductive plate as discussed above.
Furthermore the sharp anomaly decrease on the conductive side of the anomaly means that the differentiating action of the various filters described above creates an lS asymmetrical but large anomaly at the location of the contact which can easily be mistaken for the anomaly from a vertical plate. Despite the fact that the filtered anomaly from the contact-has a different shape from-that for the vertical plate (or other anomaly) it turns out in actual practice that the data from VLF surveys contain so many anomalies, of so many slightly different shapes, that it is often impossible, in either raw or filtered data sets, to further identify the source of the anomaly.
This has been a serious shortcoming of the VLF
technique. What is required is a measured or derived quantity that is indicative of those variations of the resistivity of the ground which are causing the observed anomalies. Variations in the ground resistivity can presently be measured using VLF radio waves by inserting two electrodes in the ground, aligned so that the array axis is parallel to the direction to the transmitter, and measuring the VLF potential difference between these two electrodes as described by Collett and Becker in Canadian Patent No. 795,919. However this procedure is very time consuming and detracts from the most attractive feature of the measurement of the magnetic field components, which is _ - 6 -that no electrical contact with the ground is required and thus these measurements can be carried out very quickly and inexpensively.
I have adopted an alternative approach based upon my determination of the nature of the galvanic current flow which is in actuality responsible for virtually all VLF
anomalies. In this case the effective depth of the current flow i.e. skin depth, will more or less continually reflect local resistivity, and a lateral change in resistivity will result in a change in depth of the current flow to accommodate the resulting change in skin depth. Any change in depth of the current flow in a direction perpendicular to that flow will result in a vertical component of the magnetic field, the amplitude of which will vary according to the rate of change of terrain resistivity with distance.
According to the invention, a method of surveying the resistivity characteristics of terrain comprises measuring a series of values indicative of local relative magnitudes of a vertical magnetic field component as compared to a horizontal magnetic field component due to electromagnetic radiation from a distant stationary VLF
source whilst moving along a survey line generally perpendicular to the direction of the source, so as to provide data capable of providing a record of changes in said values along said survey line, and spatially integrating the data so obtained with respect to displacement along the survey line to provide a record indicative of the resistivity characteristics of the terrain along said survey line.
Values which may be measured include the ratio of the amplitude of the in-phase component of the vertical magnetic field to the amplitude of the horizontal magnetic field occasioned by radiation from the VLF source, the ratio of the amplitudes of the total vertical and ~',' WO93/10473 2 ~ ~ 0 2 ~ g PCT/CA91/00401 _ - 7 -horizontal fields occasioned by radiation from the VLF
source, and the tilt angle or tangent of the tilt angle of the ellipse of polarization of the magnetic field occasioned by radiation from the VLF source.
The invention extends to apparatus for carrying out the above method.
Further features of the invention will be apparent from the following description with reference to the accompanying drawings.
SHORT DESCRIPTION OF THE DRAWINGS
Figures la, lb are diagrams illustrating traditional theory as to the generation of electric and magnetic fields by a distant VLF radio transmitter, along with flow of vortex currents within a radially oriented vertical plate conductor;
Figure 2a is fragmentary section through an area of terrain further illustrating traditional theory as to currents and fields in the vicinity of the conductor shown in Figure l, Figure 2b illustrating the resulting vertical magnetic field anomaly, and Figure 2c illustrating how this anomaly may be distorted by an adjacent anomaly;
Figure 3 is a diagram corresponding to Figure lb, but illustrating the galvanic current flow which is in practice dominant in the vicinity of the buried conductive plate;
Figure 4a is a diagram corresponding to Figure 3, in which the conductive plate is replaced by a contact area between zones of resistive and conductive material, and Figure 4b illustrates the vertical magnetic field anomaly in the vicinity of the contact;
- 8 - ~ 4 8 -Figure 5a is a vertical section through terrain containing various common subterranean features, Figure 5b illustrates the profile of the resulting vertical magnetic field anomaly, Figure 5c illustrates the profile obtained by difference filtering the profile of Figure 5b, and Figure 5d is a profile illustrating the effect of continued spatial integration of the profile of Figure 5b from left to right;
Figure 6 is a block diagram illustrating apparatus for generating data forming a profile such as that shown in Figure 5d;
Figures 7a and 7b are families of curves illustrating the response obtained after, spatially integrating the vertical magnetic field ratio along survey lines over an idealised vertical dike and an idealised vertical contact respectively;
Figures 8a - 8e illustrates the manipulation of actual VLF field survey data, Figure 8a plotting collected VLF vertical magnetic field data, Figure 8b showing the data horizontally differenced, Figure 8c plotting a corresponding resistivity profile prepared from field collected resistivity data utilizing spaced electrodes, Figure 8d illustrating the actual geology along the survey line, and Figure 8e plotting the data of Figure 8a spatially integrated along the survey line to simulate the operation of the method and apparatus of the invention;
Figures 9 and lOa - lOd illustrate a further example based on actual field survey data, Figure 9 showing a geological section illustrative of the type of terrain surveyed, Figure lOa being a plot of collected VLF vertical magnetic field data, Figure lOb a plot of the data after Fraser filtering, Figure lOc a plot of field surveyed conductivity data, and Figure lOd a plot of spatially Gi WO93/10473 210 0 2 i~ 8 PCT/CA91/0~01 _ g integrated VLF vertical magnetic field data along the survey line; and Figure 11 is a flow diagram illustrating one possible implementation of the invention in microprocessor controlled VLF survey equipment.
DESCRIPTION OF TH~ PREFERRED EMBODIMENTS
Referring to Figure la, the primary VLF field components generated at a point P on a conductive half-space by a distant VLF transmitter T are vertical and radial electric field components Ep and E" and a horizontal magnetic field component ~. Assuming the presence at P of - a subsurface radially and vertically extending conductive plate C (see Figures lb and 2a) the component ~ will induce a vortex current flow V in the plate which in turn will generate a secondary magnetic field ~ which will have vertical components on either-side of the plate, producing a vertical magnetic field anomaly as shown in Figure 2b, which plots the vertical magnetic field component ~ along the same section line as Figure 2a. Figure 2c illustrates how the curve of Figure 2b may be distorted if the vertical field components generated by plate C overlap with those (shown in broken lines) produced as a result of another nearby conductive body.
In practice I have found that, except under conditions where large steeply dipping and highly conductive bodies are found in terrain of generally very high resistivity, the above discussed generation of vertical field components due to vortex current flow is much less significant than the vertical field components produced by distortions of the galvanic current flow through the terrain resulting from the radial electric field component E" as illustrated in Figure 3, in which the presence of the conductor C causes a local reduction in skin depth and a 'bunching' effect as shown. As it ~ ~ ~ Q ~ ~ 8 ~
_, -- 10 --happens, this bunching in the example shown produces a vertical magnetic field anomaly very similar to that of Figure 2b.
Considering now the feature shown in Figure 4a, in which the plate C is replaced by a contact zone D
between relatively resistive material (on the left) and relatively conductive material (on the right), the galvanic current flow to the right of the contact will be concentrated nearer to the surface than that to the left of the contact. The distortion of the current flow near the contact zone produces a vertical field anomaly as shown in Figure 4b, which in real life conditions could readily be confused with a distorted anomaly as shown in Figure 2c.
Figure 5a shows schematically a section through an imaginary geological structure of somewhat greater complexity, in which, facing towards a distant VLF
transmitter, a highly conductive vertical plate like target T is flanked by zones of moderately conductive material S
which in turn contact zones of moderately resistive material R. The symbols ~ represent galvanic currents passing into the page, whilst the arrows show the direction of the resulting magnetic field. Figure 5b plots the intensity of the vertical magnetic field component ~, and Figure 5c shows the effect of differencing this field component, the plot being inverted so that the response over the wanted conductor is upward. Neither plot provides any clear picture of the underlying structure.
The effect of continually spatially integrating the vertical magnetic field profile from left to right is shown schematically in Figure 5d. The process of integration makes the presence of the two contacts, the presence of the enclosed block of conductive material, and the presence of the conductive target much more obvious;
indeed the profile of Figure Sd can be thought of as a WO93/10473 21 ~ 0 2 ~18 PCT/CA91/00401 profile of "pseudo-conductivity" and as such is a useful supplement to the raw and filtered VLF data. The profile is not, of course, an actual conductivity profile since it is generated by the magnetic field which is generated, in turn, by currents whose varying depth is defined by the skin depth, which is a function of the square root of resistivity. However, it is much more useful than the difference filtered profile, shown in Figure 5c, which yields no information about the background conductivity against which the anomalies are being detected and also makes correct identification of the source of any specific anomaly quite difficult.
Figure 6 illustrates apparatus for generating such pseudo-conductivity profiles. A generally conventional VLF receiver 2 is configured and tuned to receive at least the vertical component of anomalous secondary magnetic fields locally generated as a result of transmissions from a distant VLF transmitter located approximately perpendicular to a survey line over which the receiver is conveyed. In a preferred embodiment the relative magnitude of the in-phase vertical and horizontal components of these fields is required, and the quadrature components can be ignored, thus simplifying the receiver, which may be constructed in accordance with established practice in the VLF survey art. Reference may if needed be had to Canadian Patent No. 937,286 (Ronka) for further description of such a receiver, although in the present case simplification may be, as stated, possible where only the real component of the vertical field is to be detected.
In practice, as discussed further below, the signal which is wanted is a value indicative of the relative magnitude of the vertical magnetic field with respect to the horizontal magnetic field, in order to compensate for changes in the horizontal field.
~- - 12 - 2 ~ 9 As in the Ronka patent, the vertical and horizontal signals are received by separate channels and separately detected by a detector means 3, which must be constructed so as to preserve DC components of the demodulated signals, which are low pass filtered by filter means 4, conveniently effected using so-called integrator circuits which have a sufficiently long time constant to smooth out short term fluctuations in the output (such as those caused by modulation of the source transmitter) without obscuring terrain induced features of the signal profile which it is desired to examine. The filtered signals are combined by determining the magnitude of the vertical component as a fraction of the horizontal component, and the combined signal is passed to an analog-to-digital converter 6, having a sampling rate controlled by a timing generator producing signals related to increments of distance covered by a carrier 1, such as an aircraft, or portable instrument case, in which the receiver is carried along a survey line. If the carrier moves at a constant speed, the distance increments will be constant if the period between samples is constant. The sampling rate should be sufficient to permit the signal profile to be adequately resolved.
Rather than using an analog low pass filter, the sampling rate may be increased, and digital low pass filtering applied to produce a desired spectral content by summing groups of successive samples to provide samples with the required band-pass characteristics.
The resulting digital data is stored in a suitable memory device 8 for further processing. The data so stored may be subjected to further digital processing using similar algorithms to those discussed above in relation to prior art VLF techniques, to obtain digitized signal profiles which can be plotted or otherwise displayed by a display unit 10 in known manner to permit .. ..
WO93/10473 2 i ~ ~ 2 ~ 8 PCT/CA91/0~01 ~- - 13 -interpretation. Additionally or instead, successive samples are added in an accumulator 12 to obtain the left-to-right integration discussed above, the accumulator total after each addition being again stored in memory device 8 to build up a pseudo-conductivity profile which can be suitable scaled and plotted or displayed by display unit 10 .
The data required for operation of the invention are values indicative of the relative magnitude of a vertical field component ~ as compared to the horizontal magnetic field component ~ produced by the VIF excitation.
Examples of such values are the ratio of the in-phase vertical field component ~(I) to ~, the ratio of the total vertical field component ~ to ~, or the tilt angle or tangent of the tilt angle of the polarization ellipse of the field produced. Those values involving only in-phase components are preferred, particularly the ratio ~(I)/~.
Assuming that data representing any of the above values is measured at a series of stations spaced along the survey line at an interval of ~x metres, and that the value measured at station i in the series where i>0 is Vi, all ratios being expressed as decimal fractions, then the result of integration is X(i)=X(i-l)~V(i).~x.R
where X(-l)=0, X(i) is the spatially integrated value after station i, and R is an arbitrary constant or scale factor, and successive values for X(i) are plotted at plot points x = (i+0.5). x. If the units of ~x are metres, then a convenient value for R is unity, in which case the data can be expressed approximately in terms of "pseudoconductivity"
expressed in milliSiemens/metre (mS/m). The choice of R=l (where x is in metres) will usually result in the plotted values being of a similar order of magnitude to the actual conductivity. A notional example is shown in the following table:
~_ - 14 - ~ 4 8 -i xi=i.~x Vl Vi.~x Xi Plot Point(i+~ x 2 200 0.1 10 10 250 3 300 0.2 20 30 350 4 400 0.3 30 60 450 500 o.o 0 60 550 6 600 -0.3 -30 30 650 7 700 -0.2 -20 10 750 8 800 -0.1 -10 0 850 9 900 0.O O o 950 It should be appreciated that the apparatus used may be largely conventional, provided that it has the capability of receiving and separately demodulating, preferably synchronously, the vertical and horizontal field components. Synchronous demodulation is preferred in order to recover the vertical component which is in-phase with the horizontal component. The apparatus must of course be modified to include the integration function provided by accumulator 12. Thus the system must be able to preserve accurately the DC component of the amplitude profile of the vertical magnetic field component, relative to the amplitude of the horizontal component, and the integration of the signal must also be achieved accurately and without drift. For this reason, it is preferred that the detected signal representing the relative vertical magnetic field ratio be digitized prior to integration, which is then performed numerically. It may be advantageous to utilize digital techniques to demodulate the received VLF signal components and to perform any necessary low pass filtering required to exclude unwanted short term fluctuations in the signal. The demodulation and temporal integration technique disclosed in the Figure 5 embodiment of the Ronka ~093/10473 PCT/CA9l/0040l - 15 2 1 ~ 3 2 4 8 patent is not suitable, not only because the integration is temporal rather than spatial, but also because it is difficult to maintain sufficient DC stability for the purposes of the present invention in an analog integrator S such as that disclosed by Ronka. Furthermore, the Ronka arrangement deliberately utilizes negative feedback to modify the DC level at the output of the integrator and thus in effect greatly reduces the direct current gain of the integrating amplifier. The integrating amplifier is thus caused to act as a bandpass filter which, whilst appropriate for recording data as to the amplitude of the vertical component, renders accurate integration impossible and essentially negates the function of the amplifier as anything other than a filter.
Digital processing is advantageous since it facilitates accurate maintenance of DC levels.
Furthermore, if signal processing is performed under microprocessor control, the required accumulation or precision integration step may be performed by a simple addition to the control program which program may include routines for performing optimal filtering techniques as well as the method of my invention. A flow diagram of an exemplary routine forms Figure 11, and assumes that the samples from the digital to analog computer have been placed in a data storage file. The routine is self-explanatory, and results in the generation of a second file containing samples which represent a pseudo-conductivity plot along the same survey, in the manner illustrated in Table 1.
Before proceeding to some theoretical and field examples further illustrating the invention, two points should be made. Since integration is achieved by adding the values of data taken at discreet measurement points which are separated by a fixed interval x, the value of the integral at any distance x will be a function of 1/ x, - 16 ~ 21~024~
which is why the correction factor ~x is four~ in the equation presented above. Further, the true zer~ ~evel of the integrated data is not known. Whilst - ~s s a disadvantage, the most significant information is ~e ~ted by relative changes in the integrated signal alor 1 rvey line.
Figures 7a and 7b present results o lined from integrating the theoretical response from t- important geological targets, viz a vertical sheet-like conductor (dike) (Figure 7A) and a vertical contact (Figure 7b). The curves were calculated from theoretical data giving the normal vertical in-phase VLF magnetic field resulting from two dimensional inhomogeneities oriented with their long axis parallel to the direction to the VLF transmitter.
Consider first a vertical dike, of thickness t, and conductance S = ot, buried in a half space of resistivity P2. The integrated response becomes a simple bell-shaped curve centered over the target, with maximum value a function of S and P2, and with width proportional to 20 P2. For the extremely resistive half-space shown here (P2 = 8000 nm) the anomaly is quite wide (full width to half-maximum of about 300m), but in the more typical case of P2-1000 nm the full width to half maximum value (FWHM) becomes about lO0 m. Indeed it can be shown that, in general, the FWHM is approximately equal to the skin depth in the half-space. Furthermore the high value of P2 also yields a high peak value of the integrated anomaly; reduction of the half space resistivity results in a substantial reduction of anomaly amplitude as well as anomaly width.
Similar results but for a contact between zones of resistivity p~ and P2 respectively are shown in Figure 7b. The apparent conductivity increases as one traverses from the resistive to the conductive side of the contact by an amount which is a function of ratio Pl/P2 and which WO93/10473 2 ~ O ~ 2 4 ~ PCT/CA91/00401 ~_ - 17 -saturates for large values of this ratio. Similar curves will apply for other values of Pl/P2- The usual principles of electrical similitude apply. For example, the ratio Pl/P2 for curve B is lO. If now p~ = lO00 Qm and P2 = lO0 nm 5 (Pl/P2 = lO again), the new curve would look like curve (B) except that both the vertical and horizontal scales would be divided by (lO)~ = 3.2. The anomaly will thus be smaller and more closely confined to the contact. Such an anomaly is shown by curve (C).
The responses from these two anomalies, viz a vertical conductive dike and a vertical contact, are useful in interpreting field survey results.
(Prakla Seismos - location not known) The first case history is selected because it shows very well the relationship between VLF vertical magnetic field measurements, the integrated "pseudoconductivity" data calculated as described above, and terrain resistivity. Figure 8a shows the vertical magnetic field component measured at lO meter intervals along an 850 m survey line over the geology illustrated in Figure 8d. Figure 8b is the simple horizontal difference of the in-phase data extracted from the data shown in Figure 8a. The reversal of polarity usually applied to the 2S differenced data has not been carried out here since we wish to compare this panel with the resistivity data of Figure 8c which is the result of Schlumberger resistivity profiling with half-current electrode spacing of 205 m and half-voltage electrode spacing of 5 m. The generally good correlation between lows in the differenced in-phase vertical magnetic field data and the Schlumberger apparent resistivity lows is immediately obvious.
~~ - 18 -The results of spatially integrating the profile of Figure 8a are shown in Figure 8e. It is immediately apparent, by comparison with Figure 7b, that there is a contact at the llOm station, with resistive material to the left, conductive material to the right and about a 10/1 conductivity contact. The presence of this contact is completely missed by the raw and filtered data. In essence the data of Figure 8b is that of Figure 8e, doubly horizontally differentiated, which is exactly what it looks like. All the peaks of Figure 8e are accentuated in Figure 8b, but the all-important slowly varying component, which is related to the terrain conductivity, has been completely suppressed.
Comparing the resistivity curve shown in Figure 8c with Figure 8e, the low resistivity zones at stations 130, 220, 500 - 600 and 700 correspond extremely well to the pseudoconductivity highs at these same locations, and also the resistivity highs at stations 330, 460, 620 and 740 correspond well with the pseudoconductivity lows.
Furthermore the pseudoconductivity trace of Figure 8e depicts very clearly the fact that there is a contact with an adjacent conductor at the llOm station, a fact which no interpreter would have picked from the raw or differenced VLF data.
FIELD EXAMPLE #2 (Niagara, N.Y.) A typical geological section for the area surveyed in the second case history is shown in Figure 9.
Highly resistive dolomite bedrock 20 has at some time in the past been fractured, and the fractured zones 28 have subsequently been differentially weathered (the fracturing allows water to enter the rock and expand on freezing) to form bedrock topographic lows or troughs 22. More recently (post glacially) the bedrock surface has been covered by a lake with the result that electrically conductive ~, W093/10473 21 C ~ 2 4 ~ PCT/CA91/0~1 _ -- 19 --lacustrine (lake) clays 24 have been deposited on the bedrock to an average thickness of 3 - 8 meters. Thus the geoelectric section consists of an electrically conductive layer of variable thickness overlying an electrically resistive bedrock; fracture zones will often exist in the bedrock at the base of the troughs and these fractures may be electrically conductive as a result of infilling with water-, conductive clays or mineralization. A similar type of bedrock structure is particularly important in many countries for locating adequate sources of groundwater, in which case the interpreter wants to pick out those fracture zones which lie in the deepest troughs so that there will be maximum water available for pumping (see the line 26 indicating the level of the water table.) lS A VLF survey line was run over a geological section of this general type; survey results are shown in Figures lOa - lOd. The ~v~y line was 430 m long, and the raw VLF data shown in Figure lOa indicates three conductors located at A, B and C respectively. These conductors were picked (as is customary when anomalies are so close together as to cause each other to be vertically off-set) by choosing the points of maximum inflection. Figure lOb, which shows the Fraser filtered data, confirms this choice of conductors.
Figure lOc shows results of an additional conductivity survey over the area using a Geonics (trade-mark) EM34-3 Ground Conductivity Meter, constructed in accordance with my U.S. Patent No. 4,070,612. This survey profile shows a conductivity high located between stations 100 and 200. Since the clays are more conductive, this conductivity high corresponds to a thickening of the clay layer, i.e. to a depression or trough in the bedrock since the ground surface is virtually flat. Knowing that conductors A and B are located in the trough would make them the first priority drill targets for groundwater ~~ - 20 - ~ 4 8 -~
exploration. Target C, on the other hand, located on the flank of the trough, would be lower priority.
Unfortunately, neither the raw or Fraser filtered VLF data gives us the essential conductivity information to allow rational selection of any of the three conductors as having more groundwater potential than the others. On the other hand, the horizontally integrated VLF data, or pseudoconductivity, shown in the fourth panel, clearly shows that the conductors at stations 125 and 200 are located in a conductivity high or bedrock trough and that the conductor at station 260 is on the flank, completely confirming the results of the conductivity data, but without requiring the additional survey with a ground conductivity meter.
Whilst the invention has been described primarily with reference to the use of distant VLF radio transmitters as a VLF source, other naturally occurring or artificial sources of electromagnetic fields with frequency components in a similar or lower range may be exploited as is known in the art. Since these fields are typically non-coherent, the absolute magnitude of the vertical field component must be determined in such a case. For the purpose of interpreting the appended claims, such a source should be considered to be within the term 'VLF source'. Moreover, instead of a receiver that resolves horizontal and vertical components of the magnetic field, a receiver that directly determines a function of the tilt angle such as the angle itself or its tangent, of the polarization ellipse of the field may be utilized; the tangent of this angle is proportional to the magnitude of the vertical component for small values of the latter.
V
Claims (8)
PROPERTY OR PRIVILEGE IS CLAIMED ARE DEFINED AS FOLLOWS:
1. A method of surveying the resistivity characteristics of terrain, comprising measuring a series of values indicative of local relative magnitudes of a vertical VLF magnetic field component as compared to a horizontal VLF magnetic field component due to electromagnetic radiation at a particular frequency from a preexisting distant stationary VLF source whilst moving along a survey line generally perpendicular to the source, including the step of spatially integrating the series of values so obtained with respect to displacement along the survey line to provide a series of integrated values forming a record indicative of the resistivity characteristics of the terrain along said survey line.
2. A method according to claim 1, wherein the values detected are the ratio of the in-phase vertical to horizonal components of the magnetic field.
3. A method according to claim 1, wherein the values detected are a function of the tilt angle of the ellipse of polarization of the in-phase magnetic field.
4. A method according to claim 1, wherein samples of the relative magnitude of the vertical magnetic field data are digitized at substantially constant horizontal spatial intervals along the survey line, and the digitized samples are successively summed to provide said spatially integrated record.
5. A method of surveying the resistivity characteristics of terrain comprising a static pattern of electric ground currents at a VLF frequency, comprising detecting at multiple spaced points along a survey line values indicative of the relative values of the vertical and horizontal components of a magnetic field associated with the ground currents, and spatially integrating the values obtained along said survey line to provide integrated values indicative of the resistivity characteristics of the terrain at multiple spaced points along the survey line.
6. A method according to claim 5 wherein the points are equally spaced, the data is spatially integrated by summing, and successive summations following each sample are output as integrated values providing a measure of terrain resistivity relative to the terrain resistivity at the beginning of the survey line.
7. Apparatus for surveying the resistivity characteristics of terrain, comprising a VLF receiver constructed to quantify the magnitude of a signal induced in the receiver by a vertical VLF magnetic field component due to electromagnetic radiation from a distant preexisting stationery source of a VLF frequency relative to a horizontal VLF component due to the same source, means to move the receiver along a survey line generally perpendicular to the direction of the source, means to detect said signal whilst preserving the DC
component of the detected signal, and a low-pass filter receiving the detected signal and removing unwanted quickly changing source modulation, the apparatus further comprising an analog-to-digital converter sampling said signal at substantially constant spatial intervals along said survey line, and producing digitized samples, and a digital spatial integrator receiving and summing said samples and outputting corresponding spatially integrated samples as a function of displacement along the survey line to provide a measure of terrain resistivity at each interval relative to the terrain resistivity at the beginning of the survey line.
component of the detected signal, and a low-pass filter receiving the detected signal and removing unwanted quickly changing source modulation, the apparatus further comprising an analog-to-digital converter sampling said signal at substantially constant spatial intervals along said survey line, and producing digitized samples, and a digital spatial integrator receiving and summing said samples and outputting corresponding spatially integrated samples as a function of displacement along the survey line to provide a measure of terrain resistivity at each interval relative to the terrain resistivity at the beginning of the survey line.
8. Apparatus according to claim 7, wherein said receiver includes separate channels and means for synchronously demodulating signals induced by said horizontal and vertical field components, whereby to recover a vertical field component in-phase with the horizontal component.
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US07/561,126 US5175500A (en) | 1990-08-01 | 1990-08-01 | Method and apparatus for surveying terrain resistivity utilizing available VFL electromagnetic fields |
| PCT/CA1991/000401 WO1993010473A1 (en) | 1990-08-01 | 1991-11-11 | Method and apparatus for surveying terrain resistivity |
| EP91919746A EP0566578B1 (en) | 1990-08-01 | 1991-11-11 | Method and apparatus for surveying terrain resistivity |
| CA002100248A CA2100248C (en) | 1990-08-01 | 1991-11-11 | Method and apparatus for surveying terrain resistivity |
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US07/561,126 US5175500A (en) | 1990-08-01 | 1990-08-01 | Method and apparatus for surveying terrain resistivity utilizing available VFL electromagnetic fields |
| PCT/CA1991/000401 WO1993010473A1 (en) | 1990-08-01 | 1991-11-11 | Method and apparatus for surveying terrain resistivity |
| CA002100248A CA2100248C (en) | 1990-08-01 | 1991-11-11 | Method and apparatus for surveying terrain resistivity |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CA2100248A1 CA2100248A1 (en) | 1993-05-12 |
| CA2100248C true CA2100248C (en) | 1999-07-13 |
Family
ID=25676365
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CA002100248A Expired - Lifetime CA2100248C (en) | 1990-08-01 | 1991-11-11 | Method and apparatus for surveying terrain resistivity |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US5175500A (en) |
| EP (1) | EP0566578B1 (en) |
| CA (1) | CA2100248C (en) |
| WO (1) | WO1993010473A1 (en) |
Families Citing this family (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| ES2088005T3 (en) * | 1990-06-06 | 1996-08-01 | Univ Partnerships Pty Limited | INFRASONIC MAGNETIC MEASUREMENT INSTRUMENT. |
| US7047134B2 (en) * | 2004-03-23 | 2006-05-16 | Chemical Lime Company | Method of selecting a soil stabilization protocol for a construction site |
| US8538699B2 (en) | 2006-09-13 | 2013-09-17 | Exxonmobil Upstream Research Company | Rapid inversion of electromagnetic reconnaissance survey data |
| US9110183B2 (en) | 2006-12-06 | 2015-08-18 | Technoimaging, Llc | Systems and methods for remote electromagnetic exploration for mineral and energy resources using stationary long-range transmitters |
| US7969152B2 (en) * | 2006-12-06 | 2011-06-28 | Technoimaging, Llc | Systems and methods for measuring sea-bed resistivity |
| US8564296B2 (en) * | 2006-12-06 | 2013-10-22 | Technoimaging, Llc. | Systems and methods for remote electromagnetic exploration for mineral and energy resources |
| GB2468446B (en) | 2007-12-12 | 2011-09-21 | Exxonmobil Upstream Res Co | Method and apparatus for evaluating submarine formations |
| US20110144472A1 (en) * | 2009-12-11 | 2011-06-16 | Technoimaging, Llc | Methods of gravity and/or magnetic holographic imaging using vector and/or tensor data |
| WO2012015542A1 (en) | 2010-07-27 | 2012-02-02 | Exxonmobil Upstream Research Company | Inverting geophysical data for geological parameters or lithology |
| US8624969B2 (en) | 2010-08-02 | 2014-01-07 | Technoimaging, Llc. | Methods of electromagnetic migration imaging of geologic formation |
| WO2012024025A1 (en) | 2010-08-16 | 2012-02-23 | Exxonmobil Upstream Research Company | Reducing the dimensionality of the joint inversion problem |
| US9335432B2 (en) | 2010-08-30 | 2016-05-10 | King Abdulaziz City For Science And Technology | Semi-permeable terrain geophysical data acquisition |
| WO2012166228A1 (en) | 2011-06-02 | 2012-12-06 | Exxonmobil Upstream Research Company | Joint inversion with unknown lithology |
| EP2721478A4 (en) | 2011-06-17 | 2015-12-02 | Exxonmobil Upstream Res Co | Domain freezing in joint inversion |
| US9322910B2 (en) | 2011-07-15 | 2016-04-26 | Technoimaging, Llc | Method of real time subsurface imaging using electromagnetic data acquired from moving platforms |
| EP2734866B1 (en) | 2011-07-21 | 2020-04-08 | Exxonmobil Upstream Research Company | Adaptive weighting of geophysical data types in joint inversion |
| RU2472131C1 (en) * | 2011-07-21 | 2013-01-10 | Учреждение Российской академии наук Институт горного дела Сибирского отделения РАН | Method of generation of electromagnetic radiation of strained solid of structural material |
| US10591638B2 (en) | 2013-03-06 | 2020-03-17 | Exxonmobil Upstream Research Company | Inversion of geophysical data on computer system having parallel processors |
| US9846255B2 (en) | 2013-04-22 | 2017-12-19 | Exxonmobil Upstream Research Company | Reverse semi-airborne electromagnetic prospecting |
| CN103616725B (en) * | 2013-12-10 | 2016-05-04 | 江苏大学 | A kind of method of utilizing error log raising electromagnetic surveying quality in measured data |
| RU2542075C1 (en) * | 2014-01-09 | 2015-02-20 | Федеральное государственное бюджетное учреждение науки Институт горного дела им. Н.А. Чинакала Сибирского отделения Российской академии наук | Forecasting method of roof breaks and caving formations within section of preliminary mine working, which is not fixed with support (versions) |
Family Cites Families (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA678765A (en) * | 1964-01-28 | O. Seigel Harold | Apparatus and method of electromagnetic geophysical prospecting | |
| US3500175A (en) * | 1966-06-09 | 1970-03-10 | Geonics Ltd | Electromagnetic prospecting apparatus utilizing vlf radio fields |
| US3521153A (en) * | 1968-07-11 | 1970-07-21 | American Smelting Refining | Geophysical prospecting with electromagnetic waves of extremely low frequency |
| US3763419A (en) * | 1969-03-06 | 1973-10-02 | Barringer Research Ltd | Geophysical exploration method using the vertical electric component of a vlf field as a reference |
| US3614600A (en) * | 1969-07-31 | 1971-10-19 | Geonics Ltd | Electromagnetic prospecting apparatus for detecting electrically or magnetically responsive ore bodies |
| US3801896A (en) * | 1971-09-21 | 1974-04-02 | Co D Applic Et De Diffusion De | Method for determining the electrical resistivity of the subsoil by electromagnetic surface waves and apparatus for subsoil prospecting by measuring electromagnetic fields |
| GB1370187A (en) * | 1972-04-11 | 1974-10-16 | Barringer Research Ltd | Geophysical prospecting apparatus |
| US3950695A (en) * | 1973-04-09 | 1976-04-13 | Barringer Research Limited | Geophysical prospecting method utilizing correlation of received waveforms with stored reference waveforms |
| US3866111A (en) * | 1973-07-18 | 1975-02-11 | Exxon Production Research Co | Method of mineral exploration by detecting electromagnetic energy at power line frequency |
| US3936728A (en) * | 1973-11-29 | 1976-02-03 | Mcphar Geophysics Limited | Method and means for investigating the distribution of electrical conductivity in the ground |
| FR2272442B1 (en) * | 1974-05-21 | 1977-03-11 | Aquitaine Petrole | |
| CA1036221A (en) * | 1976-06-02 | 1978-08-08 | Geonics Limited | Method and apparatus for measuring terrain resistivity |
| US4825165A (en) * | 1978-02-08 | 1989-04-25 | Helms Ronald L | Method and apparatus for detecting a transient phenomenon by monitoring variations of an alternating component of a vertical current emanating from the earth's surface |
| US4199720A (en) * | 1978-07-27 | 1980-04-22 | Mcneill J Duncan | Method for investigating the conductivity of layered terrain |
| US4523148A (en) * | 1980-10-20 | 1985-06-11 | Gearhart Industries, Inc. | Method and apparatus for determining the spontaneous earth potential log from downhole gradient measurements |
| US4686476A (en) * | 1981-11-23 | 1987-08-11 | Mobil Oil Corporation | Plotting of magnetotelluric impedance to determine structure and resistivity variation of anomalies |
| US4686475A (en) * | 1984-07-12 | 1987-08-11 | Denver Mineral Exploration Corp. | Passive geophysical prospection system based upon the detection of the vertical electric field component of telluric currents and method therefor |
| US4591791A (en) * | 1984-07-31 | 1986-05-27 | Board Of Regents, University Of Texas System | Electromagnetic array profiling survey method |
| US4675606A (en) * | 1985-07-29 | 1987-06-23 | Suman Ganguly | Magnetometers for detecting metallic objects in earth's magnetic field |
| US4841250A (en) * | 1986-01-24 | 1989-06-20 | J. R. Jackson | Passive geophysical survey apparatus based upon the detection of the DC component of the vertical electrical potential of natural earth currents |
| US4939460A (en) * | 1987-12-14 | 1990-07-03 | Patla Louis J | Indigenous electromagnetic radiation based method for detecting geologic deposits |
-
1990
- 1990-08-01 US US07/561,126 patent/US5175500A/en not_active Expired - Fee Related
-
1991
- 1991-11-11 EP EP91919746A patent/EP0566578B1/en not_active Expired - Lifetime
- 1991-11-11 CA CA002100248A patent/CA2100248C/en not_active Expired - Lifetime
- 1991-11-11 WO PCT/CA1991/000401 patent/WO1993010473A1/en active IP Right Grant
Also Published As
| Publication number | Publication date |
|---|---|
| US5175500A (en) | 1992-12-29 |
| EP0566578B1 (en) | 1995-10-11 |
| WO1993010473A1 (en) | 1993-05-27 |
| CA2100248A1 (en) | 1993-05-12 |
| EP0566578A1 (en) | 1993-10-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA2100248C (en) | Method and apparatus for surveying terrain resistivity | |
| Wilt et al. | Crosswell electromagnetic tomography: System design considerations and field results | |
| US5563513A (en) | Electromagnetic imaging device and method for delineating anomalous resistivity patterns associated with oil and gas traps | |
| Hoover et al. | Audiomagnetotelluric sounding as a reconnaissance exploration technique in Long Valley, California | |
| Sundararajan et al. | VES and VLF—an application to groundwater exploration, Khammam, India | |
| CN109782360A (en) | A kind of low-resistance area of coverage deep sand body detection method | |
| Ibrahim et al. | Frequency-Domain Electromagnetic and VLF Data with an application to modeling of basement structures: A case study within Ibadan Area, Southwestern Nigeria | |
| Nimeck et al. | A progressive geophysical exploration strategy at the Shea Creek uranium deposit | |
| US20040000919A1 (en) | Moving source dipole electromagnetic exploration device for deeper and poorer conductors and a method of detecting such conductors | |
| Abdullahi et al. | Geophysical imaging of municipal solid waste contaminant pathways | |
| Seigel et al. | Mapping earth conductivities using a multifrequency airborne electromagnetic system | |
| Wilt et al. | Electromagnetic sounding in the Columbia Basin, Yakima, Washington | |
| Fraser¹ | Layered-earth resistivity mapping | |
| Desvignes et al. | Simultaneous interpretation of magnetic and electromagnetic prospecting for characterization of magnetic features | |
| Knödel et al. | Electromagnetic methods | |
| Dias et al. | Contribution to oil exploration and development-a successful inductive multi-frequency em survey on-shore brazil | |
| Huang et al. | Limitations of small EM sensors in resistive terrain | |
| Johnson et al. | A novel ground electromagnetic system | |
| RU2279106C1 (en) | Method for geo-electro-surveying with focusing of electric current (variants) | |
| Leite et al. | Resistivity surveys employed to study coastal aquifers in the state of Bahia, Brazil | |
| Paine et al. | Identifying and assessing ground water in the Lower Rio Grande Valley, Texas, using airborne electromagnetic induction | |
| VAN ZIJL | On the uses and abuses of the electrical resistivity method | |
| Seaborne et al. | Examples of horizontal loop electromagnetic anomalies controlled by geological faulting | |
| Chen | Application of CSAMT Method for Gold-Copper Deposits, Chinkuashih Area, Northern Taiwan | |
| Kumar et al. | Detection of fractured zones in part of hard rock area of Mirzapur District, Uttar Pradesh, using integrated geophysical method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| EEER | Examination request | ||
| MKEX | Expiry |