SCALE MODEL ELECTROMAGNETIC RESPONSE TO INLINE AND BROADSIDE SYSTEMS AT SKEW TRAVERSES OF A DIPPING HALF PLANE EMBEDDED IN A CONDUCTING HOST ROCK*

Laboratory scale model experiments have been performed to obtain the electromagnetic response of a finitely conducting half plane embedded in resistive/conductive surrounding and excited by an oscillating magnetic dipole. Inphase and quadrature profiles are presented for two horizontal coplanar transmitter-receiver systems (inline and broadside) for normal and skew traverses and for different dips of the conductor. It is observed that the broadside system is more diagnostic in delineating the strike and dip of the conductor and is more sensitive to the conducting host rock. The broadside profile over a vertical or dipping half plane is characterized, when traversing perpendicular to strike, by two positive peaks flanking a zero response when the coils are over the top edge of the conductor. For skew traverses a negative peak replaces the zero response. An increasing asymmetry in the anomalies is caused by changing the dip of the conductor from the vertical in both the systems, but it is more pronounced for the broadside system. The quadrature response in the broadside system changes in a characteristic way when the target is surrounded by a conducting host rock. The comparative results of the two systems may, therefore, be useful in the induction prospecting for ore deposits approximated by a half plane, especially in delineating the strike, dip, and effect of conductive host rock.     

THEORETICAL BROADSIDE RESISTIVITY PROFILES OVER AN OUTCROPPING DYKE

Theoretical broadside resistivity profiles over an outcropping dyke of infinite depth extent with three Wenner configurations (namely Alpha, Beta and Gamma), Schlumberger, and two-electrode configurations of various electrode spacings and for various reflection factors are presented. The broadside profiles qualitatively indicate that the shape of the anomaly is invariant with the electrode configurations. The various electrode configurations can be arranged in decreasing order of magnitude of anomalies as Beta Wenner, Alpha Wenner, Schlumberger, Gamma Wenner, and two-electrode. The broadside Wenner profiles also show larger anomalies compared to those in inline profiles.     

THE VERTICAL MAGNETIC FIELD OF AN ALTERNATING CURRENT DIPOLE FOR HORIZONTALLY STRATIFIED MEDIA *

For the computation of the vertical component H_z of the magnetic field of a horizontal A.C. dipole lying on the earth's surface, a recurrence formula is presented for a horizontally stratified half space, to obtain the (n+ 1)-layer case from the w-layer case. By means of several computed diagrams for the two-layer case, H_z can be determined for different ratios of conductivity of the subsoil and that of the overburden. Thereby the distance from the dipole as well as the layer thickness h are expressed in terms of the wave length A of a plain wave in the overburden. Assuming a sufficiently large conductivity difference, the results show that evidence about the subsurface conditions can be obtained if the distance between the measuring coil and the dipole is of the order of A/3, and if the thickness h of the layer varies within the range A/100 < h < A/6. As an example for the 3-layer case, a nonconducting intermediate layer is assumed.     

MEASURED UNDERWATER NEAR-FIELD E-PATTERNS OF A PULSED,HORIZONTAL DIPOLE ANTENNA IN AIR: COMPARISON WITH THE THEORY OF THE CONTINUOUS WAVE,INFINITESIMAL ELECTRIC DIPOLE1

At Delft Geotechnics the technique of ground-penetrating radar is in use for the detection of buried objects such as pipes. To enable us to give our ‘measurements in the field’ a more quantitative interpretation than can be deduced from these alone, a series of experiments has been started under well-defined conditions. A cylindrical vessel containing water simulates wet soil. Mounted horizontally above the water surface is a pulsed triangular half-wave dipole which is used as a transmitting antenna (TA). It has a carrier-frequency of about 160 MHz and a pulse repetition-frequency of about 50 kHz. A movable receiving dipole (‘probe’) in the water measures the transverse, mutually orthogonal E_φ,- and E_θ-components of the pulses as a function of probe-position (r, θ, φ) and of the height h of the TA above the water surface. When these pulses are Fourier-transformed, the transverse electric fields E_φ and E_θ at 200 MHz are obtained. The resulting field patterns are compared with computational results on the basis of the theory of the continuous wave, infinitesimal electric dipole (‘point dipole’). It can be concluded that:

• 1 Far-field conditions have not fully developed at a depth of about 2.50 m, the largest value of the radius r at which field patterns were measured, although it represents a distance of about 15 wavelengths.
• 2 The attenuation constant of the tapwater used, as deduced from E-field measurements for θ= 0, 2.50 m < r < 2.75 m, is slightly less than the value measured using a network analyser and air line combination, in agreement with (1).
• 3 E _φ field patterns calculated using the value of the conductivity σ corresponding to the former value of the attenuation constant agree reasonably well with the measured patterns for r≤ 2.50 m and for θ < 20° at all antenna heights considered. Calculated Eφ patterns do not agree so well with the measured patterns when h is close to zero. With increasing height the agreement inproves.
• 4 In accordance with the theory of the point-dipole, the angular distribution of the radiation patterns of the TA becomes wider as the frequency decreases.
• 5 The normalized underwater pulse-spectra shift to lower frequencies with increasing r. This can be explained since the attenuation constant of the water rises with rising frequency.

     

A numerical comparison of 2D resistivity imaging with 10 electrode arrays

Numerical simulations are used to compare the resolution and efficiency of 2D resistivity imaging surveys for 10 electrode arrays. The arrays analysed include pole‐pole (PP), pole‐dipole (PD), half‐Wenner (HW), Wenner‐α (WN), Schlumberger (SC), dipole‐dipole (DD), Wenner‐β (WB), γ‐array (GM), multiple or moving gradient array (GD) and midpoint‐potential‐referred measurement (MPR) arrays. Five synthetic geological models, simulating a buried channel, a narrow conductive dike, a narrow resistive dike, dipping blocks and covered waste ponds, were used to examine the surveying efficiency (anomaly effects, signal‐to‐noise ratios) and the imaging capabilities of these arrays. The responses to variations in the data density and noise sensitivities of these electrode configurations were also investigated using robust (L₁‐norm) inversion and smoothness‐constrained least‐squares (L₂‐norm) inversion for the five synthetic models. The results show the following. (i) GM and WN are less contaminated by noise than the other electrode arrays. (ii) The relative anomaly effects for the different arrays vary with the geological models. However, the relatively high anomaly effects of PP, GM and WB surveys do not always give a high‐resolution image. PD, DD and GD can yield better resolution images than GM, PP, WN and WB, although they are more susceptible to noise contamination. SC is also a strong candidate but is expected to give more edge effects. (iii) The imaging quality of these arrays is relatively robust with respect to reductions in the data density of a multi‐electrode layout within the tested ranges. (iv) The robust inversion generally gives better imaging results than the L₂‐norm inversion, especially with noisy data, except for the dipping block structure presented here. (v) GD and MPR are well suited to multichannel surveying and GD may produce images that are comparable to those obtained with DD and PD. Accordingly, the GD, PD, DD and SC arrays are strongly recommended for 2D resistivity imaging, where the final choice will be determined by the expected geology, the purpose of the survey and logistical considerations.     

Approximate depth of penetration in EM dipole prospecting

Summary FollowingDoll, approximate expressions for the relative contribution to the total signal by various regions of a homogeneous earth have been derived for the nine electro-magnetic dipole-dipole systems. For the three systems with parallel dipole axes, there exists a vertical focussing effect centered around a depth equal to 0.34 times the dipole spacing. For dipoles with perpendicular axes and non-zero signal, there is no focussing and the near-surface material will have an undesirably large influence on the fields measured. Considering all aspects, the system with colinear axes seems to be the best for em prospecting. For this system, the portion of the ground between approximate depths 0.1L and 1.0L accounts for the bulk of the signal,L being the dipole spacing.NGRI Contribution No. NGRI-70-160.     

LA SECHERESSE 1976 EN FRANCE: ASPECTS CLIMATOLOGIQUES ET CONSEQUENCES / The 1976 drought in France: climatological aspects and consequences

Abstract

The author recalls the four classical definitions of drought:

• — pluviometrie drought, numerically characterized by a pluviometric ratio

• — climatic drought, depending on the potential water budget

• — agricultural drought, referring to the easily available soil water storage

• — hydrological drought, recognized by an anomaly in the water supply to the water courses.

The first two types are analysed with respect to the 1976 drought in France, and the results are plotted on maps. The meteorological causes of the drought are then examined, using data from Europe and North Africa. Finally the effects of the drought on surface water (groundwater does not appear to have been affected), river transport, pollution, hydroelectric power supply, and above all on agriculture, are given.     

海洋可控源电磁法发射源姿态影响研究

海洋电磁理论是建立在水平发射偶极的基础上.然而,由于海底洋流作用海洋可控源发射偶极可能产生水平摆动、倾斜和水平旋转,造成测量数据与理想状态的水平电偶极子的电磁响应有较大的偏差.为了研究发射源姿态变化对场的影响,本文提出利用欧拉旋转将发射偶极由源坐标系转换到地球坐标系,进而分析发射源姿态变化造成海洋电磁观测数据的误差分布...     

Depth of detection of buried resistive targets with some electrode arrays in electrical prospecting

Assuming the minimum detectable anomaly to be 10%, depths of detection of a 2D vertical resistive sheet of thickness t are found to be 4.0t, 3.0t, 4.0t and 4.0t with Wenner, two-electrode, three-electrode and dipole-dipole (β-Wenner) arrays, respectively, when the array spread is in-line. On the other hand, the depths of detection obtained with a broadside spread of the arrays right over the sheet are much less and are correspondingly 2.5t, 2.0t, 2.5t and 2.5t. An increase in the depth extent W of the sheet from 10t to 20t does not increase its depth of detection with the arrays. The depths of detection of an infinitely resistive horizontal cylinder of radius R are respectively 1.5R, 1.8R, 2.0R and 2.0R with the above-listed arrays when the array spread is in-line. With broadside spread of any of the arrays, the depth of detection is seen to be 2.5R. In the case of a spherical target of radius R, the detection depths of any of the arrays are found to be small and to vary between 0.8R and 1.1R. Comparatively, the detection depths of resistive targets are much lower than those of conductive targets of the same size and shape, with any electrode array. Among all the arrays studied, the two-electrode array performs worst in the detection of resistive targets while it performs best in detecting conducting targets of limited lateral extent. In the case of a spherical target, either resistive or conductive, there is no distinct change in its detection depth with array.     

DETERMINATION OF ELECTRIC ROCK PARAMETERS FROM THE FIELD STRUCTURE OF A MAGNETIC DIPOLE*

In the theoretical part of the present paper, formulas have been analyzed for a magnetic dipole in a homogeneous and unbounded medium. The magnetic field is elliptically polarized in the region between the quasistatic zone and the far field. Since the position and the shape of the polarization ellipses depend on the complex wave number, k, it is possible to determine k by measuring the polarization ellipses. From k, the conductivity and the dielectric constant of the medium are easily calculated. The functions required for the measuring method have been computed and plotted in graphs. In the experimental part it was examined how far the theory may be applied to measurements of propagation through rock at frequencies ranging from 100-1000 kHz. These measurements showed that reasonably defined mean values of rock parameters can be given only if the deviations of the field from the theoretically expected field are not too high. These deviations have been named field distortions and have been examined by means of statistical methods (variance ratio tests). Gallery cavity and inhomogeneity or anisotropy of the medium account for these distortions.     

FREQUENCY ELECTROMAGNETIC SOUNDING USING A VERTICAL MAGNETIC DIPOLE*

A horizontal transmitter loop (vertical magnetic dipole) is used for frequency electromagnetic (FEM) soundings. The frequency ranges from approximately 6 Hz to about 4000 Hz. The vertical and radial magnetic field components are measured for 20 frequencies per decade several hundred meters from the transmitter loop. A very small bandwidth is selected for amplification using a reference signal. An Apple computer is used for data acquisition. A computer program for Marquardt inversion optimizes the parameters for the n-layer case: the resistivities and thicknesses of individual beds and a correction factor for the primary magnetic field. Interpretation of each component individually yields practically the same parameters. Examples from the field are given with interpretation; comparison with dc resistivity measurements is provided. The ratio of vertical/radial magnetic field components vs. frequency can be transformed simply into apparent resistivity vs. apparent depth. This can be done in the field to obtain an estimation of the depth of the layer boundaries. FEM results are compared with Schlumberger d.c. sounding obtained at the same site.     

MODEL TANK EXPERIMENTS ON RESOLUTION OF RESISTIVITY ANOMALIES OBTAINED OVER BURIED CONDUCTING DYKES — INLINE AND BROADSIDE PROFILING*

Model tank inline resistivity profiles with Wenner and two-electrode systems obtained over two vertical, parallel, infinitely conducting dykes submerged in water show that the two-electrode resistivity anomalies resolve better than the four-electrodes Wenner. anomalies. In contrast, the broadside resistivity anomalies obtained with a Wenner configuration resolve much better than those obtained with two-electrode configuration over an identical ground structure.     

VERTICAL MAGNETIC DIPOLE OVER TRANSITIONAL EARTH*

Harmonic steady-state solutions for magnetic and electric vector components are obtained for a layered transitional earth in a vertical magnetic dipole field. The conduction currents dominate the displacement currents. The time factor e^iωt is implied. MKS units are used.     

THE INFLUENCE OF LATERAL INHOMOGENEITIES ON THE DIPOLE METHODS*

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

单一层面高频电偶极子场的近似计算

本文对一个倾斜平面界面的高频电磁波场近似计算问题做了解答。这是钻井高频电磁波勘探方法中的基本问题。文中写出了问题的严格解,并用鞍点法对φ=0的条件,得出其近似计算式和数值结果.这个结果具有便于计算和明确的物理意义。根据计算结果,指出了特别在高导介质中,界面附近的干涉异常,必须予以充分重视。     

电偶源频率电磁测深中磁场比值视电阻率与其相位间的频散关系及资料的联合反演

探讨了在电偶源频率电磁测深中利用磁场比值定义的视电阻率ρｚｙ的物理意义及与其相位间的频散关系,给出了在频散关系的基础上电视电阻率资料估算相位的方法,并证明所估算的相位与理论模型计算值是一致的研究了ρｚｙ与ψｚｙ联合反演问题,并研制出具有较好应用效果的实用化软件。     

"Double low-points" anomaly in daily variation of vertical component of geomagnetic field before the MS8.0 Wenchuan earthquake

The "double low-points" anomaly in daily variation of vertical geomagnetic component was observed on May 9, 2008 at 13 geomagnetic observatories belonging to the geomagnetic observatory network center of China Earthquake Administration. These observatories distribute roughly on three belts with the intersection in western Sichuan. On May 12, three days after the anomaly appearance, the great M_S8.0 Wenchuan earthquake occurred. The "double low-points" anomaly in daily variation of vertical geomagnetic component is an anomalous phenomenon of regional geomagnetism, which does exist objectively. The possible cause is the change of extrinsic eddy current system resulting in geomagnetic daily quiet variation (S_q), or the delay of several hours between the intrinsic and the extrinsic eddy current systems. The relationship between the "double low-points" anomaly of daily geomagnetic variation and the earthquake reveals that the former possibly reflects the accelerative alteration of earthquake gestation in the deep Earth.     

Global force balance of region 1 current system

At times of strong solar wind forcing such as those that produce major magnetic storms, the region 1 current system dominates over the Chapman–Ferraro current system in mediating the transfer of force between the solar wind and the terrestrial system. The global force balance can be broken into two components, one involving the high-altitude part of the region 1 current system that is in contact with the solar wind (labeled here the HRS) and the other involving the low-altitude part of the region 1 current system that lies in the ionosphere (the LRS). Both communicate their J×B force to the geomagnetic dipole via a gradient in the magnetic field that they generate. In the HRS case the force acts to push the dipole away from the sun. This is the region 1 analog of the Chapman–Ferraro mechanism for transmitting the solar wind's force to the Earth. However, in the LRS case, the force (which is much stronger than in the HRS case) acts to push the dipole toward the sun, seemingly paradoxically. The LRS balances the ‘paradoxical’ sunward force on the dipole with an opposite force on the atmosphere. This paper uses MHD simulations to demonstrate the presence of both the normal force-transmitting gradient generated by the Chapman–Ferraro and the counter-Chapman–Ferraro gradient in the magnetic field generated by the region 1 current system.     

Direct interpretation of self-potential anomalies produced by a vertical dipole

In the present paper a new method is proposed for the quantitative interpretation of self-potential anomalies which are produced by a vertical dipole. First the mathematical expression of the wavenumber spectrum of the self-potential anomaly is deduced. It is pointed out that at relatively high wavenumbers the behavior of the amplitude spectrum is controlled by the closer to the surface pole at depth h. On the other hand, the “width” of the amplitude spectrum depends on the depth h and the dipole length L.Making a proper mathematical transformation of the amplitude spectrum, and applying the least squares method, it is possible to calculate the depth to the upper pole. The dipole length may then be calculated, by solving numerically a characteristic algebraic equation, as long as the “width” of the amplitude spectrum has been previously defined.The proposed method is applied on a well known self-potential profile from Greece. The calculated parameters of the polarized body are in good agreement with real data. Experimentation with synthetic models in which random noise was introduced, showed that this method gives reliable results if the noise amplitude is not more than 20% of the signal amplitude. It is clearly more efficient than the methods which are based on the model of the point pole or the dipole with a small length. It can also give good results if the horizontal extensions of the polarized body are not more than a few tenths of the depth of the upper pole. If the polarized body is tilted, the depth of the upper pole can be calculated with satisfactory accuracy.The direct interpretation method which is proposed in the present paper, may be useful in mineral exploration, and particularly if the target of interest is the detection of massive sulfide mineralization.