水平n层导电介质中时谐垂直电偶极子的电磁场

对于包含n层线性、均匀、各向同性的水平成层导电介质,本文应用汉克尔变换及其逆变换推导求得了位于顶层导电媒质中的时谐垂直电偶极子在该层产生的电磁场的表达通式．通过这些表达式能够很方便地求得包含任意多层导体半空间中的电磁场．文中最后计算了导电半体空间和二层导电介质中时谐垂直电偶极子产生的电磁场,验正了本文所得结果的正确性。     

裂缝性储层的电各向异性响应特征研究

以水平裂缝分布的孔隙介质模型为基础,建立了水平裂缝面的平行裂缝分布和带粗糙表面的点接触裂缝模型,导出了两类裂缝模型的水平和垂直电阻率响应关系,分析了两类裂缝模型在没有围压和存在围压条件下的水平电阻率、垂直电阻率变化规律和电各向异性系数的变化特征.讨论了裂缝开度、裂缝密度和裂缝粗糙度等裂缝特征参数对裂缝性储层电各向异性的影响.为简化讨论,所有电性响应特征的分析都忽略了裂缝和孔隙表面的导电性和极化的影响.利用平行分布水平裂缝面模型和带粗糙面的裂缝模型,考察了围压条件下的裂缝性岩石的电各向异性响应特征,得到了对实际应用有意义的结果.     

ELECTROMAGNETIC FIELD OF SOURCES AT THE SURFACE OF A HOMOGENEOUS CONDUCTING HALFSPACE WITH HORIZONTAL ANISOTROPY: APPLICATION TO FISSURED MEDIA*

The electric and magnetic fields generated by horizontal electric and vertical magnetic dipoles lying on the surface of a conducting medium with horizontal anisotropy are investigated. Full expressions of their Fourier transforms are given, and the fields for a vertical magnetic dipole are calculated numerically. The radial and vertical magnetic components are found to be independent of the receiver-transmitter direction, whereas the other magnetic and electric components strongly vary with this direction. These results give useful criteria for defining the direction and amplitude of anisotropy from ground data; a ground experiment on fissured limestone was found to confirm the expected variations of the various field components. It is believed that this electromagnetic method can be used in order to provide information about the direction and amplitude of electric anisotropy.     

Methods and results of monitoring of the natural telluric field in the Baikal rift zone

Variations in horizontal and vertical components of the telluric field are studied within the range of periods from a few tens of seconds to diurnal and longer variations. Variations caused by an external ionospheric source are used for studying the geoelectric heterogeneity of the medium. For this purpose, an electric tipper is introduced by analogy with the geomagnetic field and its frequency response is examined. The magnetotelluric impedance is studied using data of remote geomagnetic observatories. According to the data of monitoring of the electric tipper and magnetotelluric impedance, anomalous effects are revealed that can be related to earthquakes with K ≥ 12. A method is developed for the extraction of electric field anomalies produced by intraterrestrial electric sources.     

Two-Dimensional Magnetotelluric Modelling of Ore Deposits: Improvements in Model Constraints by Inclusion of Borehole Measurements

A combination of magnetotelluric (MT) measurements on the surface and in boreholes (without metal casing) can be expected to enhance resolution and reduce the ambiguity in models of electrical resistivity derived from MT surface measurements alone. In order to quantify potential improvement in inversion models and to aid design of electromagnetic (EM) borehole sensors, we considered two synthetic 2D models containing ore bodies down to 3000 m depth (the first with two dipping conductors in resistive crystalline host rock and the second with three mineralisation zones in a sedimentary succession exhibiting only moderate resistivity contrasts). We computed 2D inversion models from the forward responses based on combinations of surface impedance measurements and borehole measurements such as (1) skin-effect transfer functions relating horizontal magnetic fields at depth to those on the surface, (2) vertical magnetic transfer functions relating vertical magnetic fields at depth to horizontal magnetic fields on the surface and (3) vertical electric transfer functions relating vertical electric fields at depth to horizontal magnetic fields on the surface. Whereas skin-effect transfer functions are sensitive to the resistivity of the background medium and 2D anomalies, the vertical magnetic and electric field transfer functions have the disadvantage that they are comparatively insensitive to the resistivity of the layered background medium. This insensitivity introduces convergence problems in the inversion of data from structures with strong 2D resistivity contrasts. Hence, we adjusted the inversion approach to a three-step procedure, where (1) an initial inversion model is computed from surface impedance measurements, (2) this inversion model from surface impedances is used as the initial model for a joint inversion of surface impedances and skin-effect transfer functions and (3) the joint inversion model derived from the surface impedances and skin-effect transfer functions is used as the initial model for the inversion of the surface impedances, skin-effect transfer functions and vertical magnetic and electric transfer functions. For both synthetic examples, the inversion models resulting from surface and borehole measurements have higher similarity to the true models than models computed exclusively from surface measurements. However, the most prominent improvements were obtained for the first example, in which a deep small-sized ore body is more easily distinguished from a shallow main ore body penetrated by a borehole and the extent of the shadow zone (a conductive artefact) underneath the main conductor is strongly reduced. Formal model error and resolution analysis demonstrated that predominantly the skin-effect transfer functions improve model resolution at depth below the sensors and at distance of \(\sim \) 300–1000 m laterally off a borehole, whereas the vertical electric and magnetic transfer functions improve resolution along the borehole and in its immediate vicinity. Furthermore, we studied the signal levels at depth and provided specifications of borehole magnetic and electric field sensors to be developed in a future project. Our results suggest that three-component SQUID and fluxgate magnetometers should be developed to facilitate borehole MT measurements at signal frequencies above and below 1 Hz, respectively.     

Near surface geophysical surveys with a high frequency mutual impedance measuring system

Relationships for the mutual impedance of horizontal coplanar, vertical coplanar, perpendicular and vertical coaxial loop antenna arrays operating at a high frequency, located over a layered conducting medium, have been derived. The relationships take into account the occurrence of displacement currents in both the free space and the conducting medium. The relations can be used for calculating the earth parameters in geophysical prospecting. The use of high frequency antenna array makes it possible to determine distributions of apparent resistivity as well as apparent electrical permittivity. The potential application of the measurements of mutual impedance of loop antennas operating at a high frequency to the geophysical examinations is presented. The measuring system developed for this purpose is described. A simple method of interpreting the mutual impedance measurements to obtain the traces of apparent resistivity and apparent permittivity of the investigated earth is given.     

TUTORIAL DISTRIBUTION OF ALTERNATING ELECTRICAL CHARGES IN A CONDUCTING MEDIUM*

The main features of the distribution of volume and surface charges in a conducting medium can be described separately for direct and alternating electromagnetic fields. The density of charges depends on the conductivity of a medium and on the electrical field. The relation is particularly simple for the quasi-stationary field, i.e., when the influence of displacement currents is negligible. Conditions are formulated under which electrical charges arise in a conducting medium: electrical charges are shown to exist for direct and quasi-stationary fields when there is a component of electric field parallel to the gradient of conductivity. The density of these charges is proportional to the applied electric field.     

时变海面与其上方导体平板的复合电磁散射研究

应用当前的数值方法求解动态海面与目标的复合电磁散射,由于海面的变化,在不同时刻需要对阻抗矩阵各元素及海面表面电流重新求解,因而要耗费大量内存和运算量.为了克服这一问题,本文应用物理光学(PO)近似求解了导体海面表面电流及导体平板的一阶散射场,应用基尔霍夫近似给出了海面的后向散射场,同时借助互易性定理降低了求解平板和海面之间二次耦合散射场的难度,讨论了平板尺寸、风速等对后向复合散射场的影响.另外,本文还推导出了耦合散射场多普勒谱频移的理论公式,详细分析了复合后向散射场的Doppler 频谱特性.     

Correlation of the anomalies in the electric field and electric conductivity of the lithosphere to earthquakes in Kamchatka

The data of long-term electromagnetic monitoring are used for studying the dynamics of electric conductivity of the medium and the electric field of the terrestrial sources. The electric conductivity of the medium is estimated from the magnetotelluric transfer functions (impedance tensor and telluric tensor). The electric field of terrestrial sources is identified by filtering the variations of the observed electric field of the Earth. The magnetotelluric parameters and the electric field of terrestrial sources feature anomalous changes of supposedly earthquake-related origin. The anomalies associated with the same earthquake are not simultaneous. It is shown that these anomalies are generated by processes occurring at different depths. The strong earthquake is preceded by the appearance of surface anomalies several months before the event and accompanied by a deep coseismic anomaly. The probable nature of the recognized anomalies is discussed.     

Electromagnetic analogue model measurements and finite-difference numerical calculations of the response of a conducting slab to three different source fields

This work compares experimental analogue model measurements and finite-difference numerical calculations of the electric and magnetic fields for a highly conducting slab embedded in a poorly conducting host earth for three different source field configurations. Measurements and calculations were carried out for a uniform source, a sheet current source with a y exp(?ay) current intensity distribution, and a horizontal magnetic dipole source. The results indicate reasonable agreement with some exceptions between the analogue and numerical methods. The source field is found to have an important effect on the field anomalies at the interface of the highly conducting slab and the poorly conducting host medium.     

Vertical and horizontal components of the lake Baikal electrotelluric field and their relation to the electric conductivity

Synchronous observations of the electrotelluric field were performed on the lake ice and the coast. The ice-based vertical/horizontal field ratio at periods of a few seconds to a few minutes ranges from a few hundredths to a few thousandths, which is comparable to the noise level. The results of these observations favor a plane wave model. At the coastal station, the vertical and horizontal fields are similar, which is related to geoelectric heterogeneity of the medium. An electric tipper is used to control the behavior of heterogeneity in the spectral-time domain. Anomalous changes associated with earthquakes of K ≥ 12 are revealed from data of the tipper monitoring at epicentral distances of up to 300 km. A longitudinal MTS curve free from the effect of local geoelectric heterogeneities is obtained from ice-based data. Deep conducting layers are recognized as minimums in this curve.     

Magnetic and Electric Fields Produced in the Sea During Geomagnetic Disturbances

— To understand geomagnetic effects on systems with long conductors it is necessary to know the electric field those systems experience. For surface conductors such as power systems and pipelines this can easily be calculated from the magnetic field variations at the surface using the surface impedance of the earth. However, for calculating the electric fields in pipelines and submarine cables at the seafloor it is necessary to take account of the attenuating effect of the conducting seawater. Assuming that the fields are vertically propagating plane waves, we derive the transfer functions between the electric and magnetic fields at the seafloor and the magnetic field variations at the sea surface. These transfer functions are then used, with surface magnetic field data, to determine the power spectra of the seafloor magnetic and electric fields in a shallow sea (depth 100 m) and in the deep ocean (depth 5 km) for different values of the Kp magnetic activity index. For the period range considered (2 min to 3 hrs) the spectral characteristics of the seafloor magnetic and electric fields for a 100 m deep sea are very similar to those of the surface fields. For the deep ocean the seafloor spectra show a faster decrease in spectral density with increasing frequency compared to the surface fields. The results obtained are shown to be consistent with seafloor observations. Assessment of the seafloor electric fields produced by different levels of geomagnetic activity can be useful in the design of the power feed equipment for submarine cables and cathodic protection for undersea pipelines.     

SOME REMARKS ON THE BEHAVIOR OF THE MAGNETOTELLURIC PHASE*

The magnetotelluric phase in the medium of incidence over a good conductor varies from its value at the conductor surface, usually close to 45°, to the characteristic standing-wave phaseshift of 90° over a distance of a few conductor skin-depths, much shorter than the vacuum wavelength at the frequencies considered. Over the top of a graben or at the surface of a two-layer structure the basement can be looked upon as the matrix. If the graben or the overburden is a better conductor it will carry an extra current density in phase with the electric field, thus giving rise to an additional magnetic field also in phase with the electric field. This drives the argument of the surface impedance toward smaller values. For a more resistive graben or overburden the current is depressed, corresponding to a reduction of the surface magnetic field by an amount in phase with the electric field. This effect drives the phase of the impedance to values above the characteristic 45° of a uniform conductor.     

Reduction of reflections from above surface objects in GPR data

During a ground-penetrating radar (GPR) survey, special attention must be paid to objects located above the earth's surface. Due to the low-loss character of electromagnetic propagation in air and high velocity, above-surface reflections or diffractions can overwhelm subsurface events, making the interpretation a difficult task. The relative sensitivity of reflections and diffractions originating from above-surface objects is a function of the antenna radiation characteristics, the lateral and vertical dimensions of the objects and their position with respect to the antennas. The largest amplitude reflections and diffractions are expected when the polarization of the electric field is parallel to the long-axis of the object. Near the surface in the E-plane, the electric field is vertically polarized and has a larger amplitude than the horizontally polarized electric field in the H-plane. Numerical modeling of reflections from three above surface objects (a vertical plane and elongated horizontal and vertical objects) demonstrate that the largest amplitude difference occurs when an elongated vertical object is present in the E- or H-plane. The calculated reflection from the elongated vertical object present in the E-plane was 21 times larger than when it was present in the H-plane. In 60-m long field data sets, reflections from interfering trees present in the E-plane were at several positions >15 times larger and on average 6 times larger than when the trees were present in the H-plane. These large amplitude differences indicate that appropriate orientation of the antennas can be used to minimize the effects of above-surface reflections and diffractions.     

Magnetotelluric response of a two-dimensional sloping contact by the finite element method

Summary The finite element method, with triangular elements, is used to study the effect of a two-dimensional sloping contact on the surface electromagnetic fields. It is found in the case ofH-polarization and small slopes that the electric field and the apparent resistivity near the contact, on the conductive side, are higher than their asymptotic values. In the case ofE-polarization the apparent resistivity and phase values on the conductive side fall off less rapidly to their asymptotic values with decreasing slope resulting in higher apparent resistivity and phase values on the conductive side, than those expected for a vertical contact. The peak in the amplitude and phase of the normalized vertical magnetic field shifts from the resistive side for a vertical contact to the conductive side for a sloping contact. Far from the sloping contact, on the conductive side, higher values are observed for the normalized vertical magnetic field than in the case of a vertical contact.     

SOME EFFECTS OF MULTIPLE LATERAL IN HOMOGENEITIES IN MAGNETOTELLURICS1

The analytical solution, for the H-polarization magnetotelluric impedance, of a series of multiple, vertical, conducting slabs (dikes) embedded in a host medium is extended to an infinite array in order to model an anisotropic layer. The solution is used to study the effects of such strongly anisotropic media on the surface impedance. At low frequencies such vertically laminated structures behave as a bulk material. It is shown that the effective bulk parameters are those expected from d.c. theory. However, practical electrode separations may not be long enough or adequately positioned to correctly obtain these bulk parameters from the measured impedance. Thus, such structures can masquerade as quite different one-dimensional structures. A haphazard use of long electrode spacings will not necessarily produce correct results.     

Anisotropy Versus Heterogeneity in Continental Solid Earth Electromagnetic Studies: Fundamental Response Characteristics and Implications for Physicochemical State

Electrical anisotropy, the effect of current density in a medium being a function of the orientation of the electric field, is being recognized increasingly as an important effect in explaining Earth electromagnetic observations. A consideration of anisotropy, however, in most cases is an admission of spatial aliasing in earth structure, wherein the averaging volume of diffusive EM fields may be greater than the characteristic dimensions of a family of oriented structures, thus leading to a response which is equivalent to a bulk anisotropic medium. Even for two-dimensional geometries, there can be strong non-parallelism of principal axes of vertical magnetic field relative to the impedance over broad areas, as well as impedance phase variations which leave normal quadrants, if there are multiple directions of anisotropy or anisotropy strike distinct from bulk geometric (2D) strike. This paper concentrates on experience with regional field studies in continental settings where bulk anisotropy is apparent. Upper crustal anisotropy may result from preferred orientations of fracture porosity, or lithologic layering, or oriented heterogeneity. Lower crustal anisotropy may result from preferred orientations of fluidized/melt-bearing or graphitized shear zones, but does not necessarily reflect current state of stress per se. In the upper mantle, the prior causes all may act in pertinent domains, but added to these is the possibility of strong electrical anisotropy due to hydrous defects within shear-aligned olivine crystals (solid-state conduction). Several field examples from continental MT investigations will be discussed, which roughly fall into active transpressional, active transtensional, and fossil transpressional regimes. A general challenge in interpreting data with apparent anisotropic effects is to establish the tradeoff between heterogeneity and anisotropy in the inversion of EM responses.     

Possible configurations of the magnetic field in the outer magnetosphere during geomagnetic polarity reversals

Possible configurations of the magnetic field in the outer magnetosphere during geomagnetic polarity reversals are investigated by considering the idealized problem of a magnetic multipole of order m and degree n located at the centre of a spherical cavity surrounded by a boundless perfect diamagnetic medium. In this illustrative idealization, the fixed spherical (magnetopause) boundary layer behaves as a perfectly conducting surface that shields the external diamagnetic medium from the compressed multipole magnetic field, which is therefore confined within the spherical cavity. For a general magnetic multipole of degree n, the non-radial components of magnetic induction just inside the magnetopause are increased by the factor 1 + [(n + 1)/n] relative to their corresponding values in the absence of the perfectly conducting spherical magnetopause. An exact equation is derived for the magnetic field lines of an individual zonal (m = 0), or axisymmetric, magnetic multipole of arbitrary degree n located at the centre of the magnetospheric cavity. For such a zonal magnetic multipole, there are always two neutral points and n – 1 neutral rings on the spherical magnetopause surface. The two neutral points are located at the poles of the spherical magnetopause. If n is even, one of the neutral rings is coincident with the equator; otherwise, the neutral rings are located symmetrically with respect to the equator. The actual existence of idealized higher-degree (n > 1) axisymmetric magnetospheres would necessarily imply multiple (n + 1) magnetospheric cusps and multiple (n) ring currents. Exact equations are also derived for the magnetic field lines of an individual non-axisymmetric magnetic multipole, confined by a perfectly conducting spherical magnetopause, in two special cases; namely, a symmetric sectorial multipole (m = n) and an antisymmetric sectorial multipole (m = n – 1). For both these non-axisymmetric magnetic multipoles, there exists on the spherical magnetopause surface a set of neutral points linked by a network of magnetic field lines. Novel magnetospheric processes are likely to arise from the existence of magnetic neutral lines that extend from the magnetopause to the surface of the Earth. Finally, magnetic field lines that are confined to, or perpendicular to, either special meridional planes or the equatorial plane, when the multipole is in free space, continue to be confined to, or perpendicular to, these same planes when the perfectly conducting magnetopause is present.Also Honorary Research Associate, Rutherford Appleton Laboratory, Chilton, Didcot, Oxon OX11 0QX, UK and Visiting Reader in Physics. University of Sussex, Falmer, Brighton BN1 9QH, UK     

First in situ measurement of electric field fluctuations during strong spread F in the Indian zone

An RH-560 rocket flight was conducted from Sriharikota rocket range (SHAR) (14°N, 80°E, dip 14°N) along with other experiments, as a part of equatorial spread F (ESF) campaign, to study the nature of irregularities in electric field and electron density. The rocket was launched at 2130 local time (LT) and it attained an apogee of 348 km. Results of vertical and horizontal electric field fluctuations are presented here. Scale sizes of electric field fluctuations were measured in the vertical direction only. Strong ESF irregularities were observed in three regions, viz., 160/190 km, 210/257 km and 290/330 km. Some of the valley region vertical electric field irregularities (at 165 km and 168 km), in the intermediate-scale size range, observed during this flight, show spectral peak at kilometer scales and can be interpreted in terms of the image striation theory suggested by Vickrey et al. The irregularities at 176 km do not exhibit any peak at kilometer scales and appear to be of a new type. Scale sizes of vertical electric field fluctuations showed a decrease with increasing altitude. The most prominent scales were of the order of a few kilometers around 170 km and a few hundred meters around 310 km. Spectra of intermediate-scale vertical electric field fluctuations below the base of the F region (210/257 km) showed a tendency to become slightly flatter (spectral index n = –2.1 ± 0.7) as compared to the valley region (n = –3.6 ± 0.8) and the region below the F peak (n = –2.8 ± 0.5). Correlation analysis of the electron density and vertical electric field fluctuations suggests the presence of a sheared flow of current in 160/330 km region.