地磁感应电流（GIC）的作用与评估

地磁感应电流（GIC）可能对各种人工长距离导电体造成影响与危害.地磁扰动产生的感应电场的强度与地磁场强度、地下电阻率结构相关，在导电系统内生成的GIC的强度则同时与导电系统的内在结构有关.计算了加拿大Manitoba省三个典型地区在2000年7月15日的一个强烈磁暴期间产生的感应电场.通过对地磁活动性的统计分析，估计加拿大魁北克电网可能经受的最大GIC达每相78A（一年一次）和234A（每十年一次）.     

DETECTION OF AN AIR-FILLED DRAINAGE GALLERY BY THE VLF RESISTIVITY METHOD1

VLF wave impedance measurements have been successfully used to detect air-filled drainage galleries near the town of Alcala, Spain. The galleries are detectable by H-polarization electric field measurements due to the electric field anomalies associated with the galleries and overlying gravel deposits. The forced deviation of the primary current flowlines around the 2D void results in a higher-than-normal apparent resistivity and a relative phase low above the gallery. The findings support earlier theoretical predictions that at very low frequencies (VLFs), galvanic current effects may dominate over vortex currents in moderately conductive terrains. Theoretical modelling confirmed that for a resistive target no detectable E-polarization response can be expected from either magnetic or electric field measurements since current line deviations and vortex effects are negligible under such circumstances. The results demonstrate the importance of using at least two orthogonal VLF transmitters in order that anomalies arising from both galvanic and inductive effects may be identified, irrespective of orientation.     

INTERPRETATION OF VLF-EM IN-PHASE DATA USING CURRENT DENSITY PSEUDOSECTIONS1

The interpretation of VLF-EM surveys in terms of buried conductors can be assisted by the application of a linear filter to the observed in-phase component of the vertical magnetic field. One such scheme is examined critically by using the calculated response from a variety of synthetic models to compute theoretical current density pseudosections. The results confirm that this filter technique provides a useful complementary tool for studying the third (i.e. depth) dimension. For single, steeply dipping plates (> 45°) diagnostic information may be derived concerning the depth, size, lateral location, and direction of dip, since the current density maxima seems always to occur within the conductor or at least within one data interval. However, there are some limitations which do not appear to be widely recognized. For single plates the angle of dip cannot be resolved by the current density transformation. Pattern distortions can occur where targets are in close proximity, where the cross-sectional form of the conductor is complex or where the dip is shallow (< 45°). In these latter cases the current density maxima may not occur within the conductive structure and therefore cannot be used to infer depth of burial and/or conductor shape with the same degree of confidence.     

中低纬度电离层对太阳风状态的响应

本文讨论了行星际磁场B₂分量变化时内磁层和中低纬度电离层的响应.指出B₂变化引起的磁层大尺度对流电场的变化在一定条件下有可能透入内磁层,并沿磁力线映射到中低纬度电离层,在那里产生电场和电流体系,从而使S_q电流体系发生畸变,并在地面磁场中反映出来.数值计算表明,当△B₂<0时,S_q电流体系的焦点向东和向高纬移动,地面磁场会观测到数伽马的变化.这就为中低纬地磁观测诊断磁层和太阳风状态提供了一种可能性.此外,本文还用上述物理过程解释了赤道地区一些高空物理现象,如B₂倒转时电离层漂移速度的变化,赤道磁场异常以及赤道q型偶现E层的消失等等.     

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.     

Processing of magnetotelluric data

The processing of magnetotelluric data involves concepts from electromagnetic theory, time series analysis and linear systems theory for reducing natural electric and magnetic field variations recorded at the earth's surface to forms suitable for studying the electrical properties of the earth's interior.The electromagnetic field relations lead to either a scalar transfer impedance which couples an electric component to an orthogonal magnetic component at the surface of a plane-layered earth, or a tensor transfer impedance which couples each electric component to both magnetic components in the vicinity of a lateral inhomogeneity.A number of time series spectral analysis methods can be used for estimating the complex spectral coefficients of the various field quantities. These in turn are used for estimating the nature of the transfer function or tensor impedance. For two dimensional situations, the tensor impedance can be rotated to determine the principal directions of the electrical structure.In general for real data, estimates of the apparent resistivity are more stable when calculated from the tensor elements rather than from simple orthogonal field ratios (Cagniard estimates), even when the fields are measured in the principal coordinates.     

VLF ANOMALIES FROM A PERFECTLY CONDUCTING HALF PLANE BELOW AN OVERBURDEN*

A theoretical solution to the electromagnetic problem of a perfectly conducting half plane below a conducting overburden has been obtained. The VLF anomalies have been computed for different overburden conductivity and thickness and also for different dip angles of the half plane. In the computations the contribution to the secondary magnetic field from the electric Hertz potential has been neglected. The anomaly curves which are displayed as EM 16 readings, show a fairly complicated behaviour. This is mainly due to the phase shift and attenuation of the field caused by the conductivity of the overburden and the host rock. From the anomaly curves it is possible to define the apparent depth to the top of the conductor as the distance between the peak value and the cross-over of the real component. The apparent depth is usually larger than the actual depth, but it is possible to determine the actual depth to the conductor from the relation between the peak-to-peak anomaly and the apparent depth. When the peak-to-peak anomaly is fairly large, it is also possible to make estimates of the dip angle. However, a complete set of master curves will be a necessary tool for interpretation of VLF data when there is need to obtain more accurate estimates of the half plane parameters. In a specific case the theoretical calculations are shown to be in good agreement with measured data.     

A METHOD FOR CALCULATING IP ANOMALIES FOR MODELS WITH SURFACE POLARIZATION*

A numerical method is given for calculating resistivity and induced polarization anomalies produced by a surface polarization model. Surface polarization is generated when a purely electronic conductor is located in an electrolyte environment. The system that develops on the boundary between the conductor and the electrolyte is described macroscopically by a net surface charge distribution and an electric double layer. An integral equation is derived for the potential by assuming that the electronic conductor forms an equipotential system and that the polarization impedance across the boundary is linear. The integral equation is solved by means of the method of subsections. As an application some numerical modeling results are presented. The surface impedance values used in calculations are based on laboratory measurements that are briefly described. Implications of the results for scale modeling are discussed.     

覆岩变形破坏电阻率响应特征的模拟实验研究

地下采矿活动引起上覆岩层移动变形和破坏,对覆岩变形破坏进行监测,可以为采矿设计和安全生产提供可靠的技术保障. 本文从对应于实际地层介质的相似材料物理模型实验出发,探讨了模拟实验的相似条件和可行性,进行了采矿活动不同时期覆岩电阻率特征变化的动态测试,结合实验数据进行反演计算,最后获得覆岩变形破坏与电阻率变化的响应关系. 实验结果表明,覆岩变形和破坏必然造成地层介质电性特征的变化,不同破坏带的破坏程度不同,电阻率变化响应特征不一样,两者是一一对应的. 利用电阻率方法监测覆岩变形和破坏是完全可行的.     

Effect of anisotropy of the earth on the impedance measurements

Summary Electromagnetic field response of a homogeneous uniaxially anisotropic conducting medium has been investigated in the presence of a vertical electric dipole placed at a heighth above the conducting surface. Expressions for the horizontal components of electric and magnetic fields and hence that of vertical impedance of the field, as a function of the distance of the place of observation from the source, have been derived. Numerical computations have been done and curves are given to illustrate the influence of anisotropy on the amplitude values of the impedance on the surface of the conducting medium.N.G.R.I. Contribution number 71–279.     

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.     

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.     

Grounded-source TEM modelling of some deep-seated 3D resistivity structures

Long-offset transient electromagnetics (LOTEM) is now regarded as a suitable electrical method for deep exploration along with magnetotellurics (MT). In this method, the vertical magnetic-field impulse response and, occasionally, the horizontal electric-field step response of a grounded-wire source on the surface of the earth are measured. Here, these two responses are computed for 3D models of three deep resistivity structures of interest in hydrocarbon exploration: (i) a faulted graben in a resistive basement rock at a depth of 4 km beneath a conductive overburden; (ii) a facies change in a resistive layer buried at a depth of 2 km in the conductive overburden above a resistive basement; and (iii) an anticlinal uplift of a resistive layer at a depth of 1 km in the conductive overburden above a resistive basement. The results show that the sensitivity of the electric-field response to model perturbation is generally greater than that of the magnetic-voltage response. Further, the electric-field sensitivity is confined to early and intermediate times while that of the magnetic-voltage response is confined to intermediate and late times. Hence it is recommended that both electric and magnetic recordings are made in a LOTEM survey so that the final results can be presented as apparent-resistivity curves derived from the two responses jointly as well as separately.     

NUMERICAL MODELLING OF A NEW MULTI-FREQUENCY ELECTROMAGNETIC SYSTEM*

The amplitude of the horizontal magnetic field in the ground between two parallel wires, both carrying an alternating current in the same direction, is likely to have a saddle point if the separation between the wires is small and the frequency is low. The amplitude has a maximum in the vertical direction and a minimum in the horizontal. Rectangular geological structures in the ground which are centered between the wires have a varying effect on the magnetic fields at the surface. In general, the vertical magnetic field “crosses over” at the center of the structure. A shallow and flat lying conductor displays a broad flat type of profile when the horizontal magnetic field between the wires is measured. Changing the structure to a narrower but more conducting one at depth will provide a more pointed but still broad profile. The phase of the horizontal field is also increased. When the structure is a thin vertical dyke, the amplitude of the horizontal magnetic field anomaly due to the dyke rapidly decreases as the depth of the dyke is increased. The phase of the horizontal field is less sensitive to changes in depth of the dyke but is more sensitive to the conductivity ratio of the dyke and the half-space. The amplitude of the vertical magnetic field anomaly due to the dyke is only slightly influenced by conductivity contrast or the depth of the dyke. The phase of the vertical magnetic field, however, is strongly influenced by the conductivity contrast, particularly if the conductivity frequency product is greater than hundred. In essence, the field behaves like that of the conventional vertical loop source, but the fields are uniform over much larger areas. This suggests the possibility of using dip angle measurements for rapid reconnaissance.     

THE DISPERSION RELATION, OF EARTH''S ELECTROMAGNETIC RESPONSE FUNCTIONS AND THE JOINT INVERSION OF DATA

On the basis of the dispersion relations of MT field, the necessity and applied prospects of the joint inversions using a pair of MT response functions which are correlative with the dispersion relations, are infered. A filter coefficient algorithm is made, with which the corresponding impedance phase data can be estimated using a set of apparent resistivities. The tests for the observed MT data show that when comparing the impedance phase estimated using the dispersion relation with the ob served phase, it can be checked whether the dispersion relation between observed apparent resistivity and phase data is satisfied or not, and that the use of the phase data corrected using the dispersion relation in the joint inversion is advantageous to obtain more confident results. It is shown that joint inversions are more advantageous than single parameter inversions, and that in the most case the joint inversion using the apparent resistivities of impedance real and imaginary parts is more advantageous than the jointinversion using the normal apparent resistivity and impedance phase. The existence of the dipersion relations between the ratio apparent resistivity and corresponding impedance phase of the orthogonal electric and magnetic field horizontal Components in the frequency EM sounding with horizontal electric dipole(FEMS) are discussed, the better effect of the joint inversion using the pair of EM response functions is obtained. The problems on the one-dimensional joint inversion for the MT and FEMS apparent resistivities, for which the observed frequency bands partly overlape each other, are studied. It is shown that this joint inversion is applicable and effective:the joint inversions of the practical data for two kinds of EM methods at two sites give the results well corresponding to the drilling data. The simulated MT inversions for the data of two kinds of EM methods are made, and more confident results also are obtained.     

Specific features of magnetic barrier formation near the magnetopause

The numerical three-dimensional MHD model is used to study the formation of the magnetic barrier in the inner part of the magnetosheath near the magnetopause. The set of the quasistationary solutions for several characteristic directions of the interplanetary magnetic field (IMF) has been obtained: for northward and southward IMF, for the direction along the Parker helix (at an angle of 45° with respect to the Sun-Earth line), and for the predominantly radial direction (at an angle of 20° with respect to the Sun-Earth line). The mechanism used to take into account the effect of magnetic reconnection at the magnetopause on a flow in the magnetosheath is introduced in the case of southward IMF. The results of the calculations indicate that the magnetic field absolute value in the magnetic barrier reaches its maximal value when IMF is northward. The introduction of magnetic reconnection at southward IMF can result in an insignificant decrease in the field value. However, the model predicts that a decrease in the magnetic field is much more substantial when the IMF direction is close to radial.     

THE AMPLITUDE AND PHASE RESPONSE OF A SEISMIC VIBRATOR*

The amplitude and phase response of a simple model is compared with the performance of a real vibrator working in the field. The field results show a characteristic phase response which confirms that the real drive force applied to the baseplate and its load impedance is faithfully represented by the acceleration of the reaction mass. It follows that all the parameters necessary to calculate the load impedance and the true power dissipated in the earth can be measured at the output of the vibrator. It also follows that the current method of baseplate phase compensation should be reconsidered.     

Conductivity anomalies: lower crust or asthenosphere?

The indication from surface measurements of a zone of relatively high conductivity (resistivity<200 ohm-m) at depths between 20 and 50 km has become so general over the Earth that regions without this zone can be considered anomalous. However, the depths indicated often span the lower crust and the uppermost mantle, so that before any effect can be definitely attributed to one or the other, the depth resolution in the electromagnetic measurements must be carefully considered. This paper applies the eigenvector decomposition of generalized linear inverse theory to soundings by Schlumberger resistivity, by magnetotellurics, by man-made electromagnetic fields formed by controlled current flow in grounded electric transmission lines, and by natural magnetic field variation studies to improve the bounds on depth, thickness and conductivity of a conductive layer. It is shown that many of the methods are capable of giving the depth to the top of a conductor with remarkable accuracy. Joint inversion of two or more of them offers an advantage in the separation of thickness and conductivity of both conductive and resistive layers. Natural geomagnetic field transfer functions, while accurately mapping the position of the edge of a conductor, do not provide the resolution of the other techniques, largely because the frequencies that can be practically measured at present are much too low.     

Horizontal electrojet associated with the sun-aligned arcs

Analysis of two strong sun-aligned arcs over the Canadian Eureka Observatory (89° CGM) near the north magnetic pole and accompanying related ground magnetometer measurements has identified an electrojet current of ∼1×10⁴ A flowing within the arcs in a sunward direction. The electrojet current was carried by low energy electrons created by impacting precipitation drifting at E×B/B² velocity within the arcs, where E is the dawn-to-dusk electric field. One of the arcs moved rapidly in a dawn to dusk direction. The measured arc velocity was 365 m/s. This agrees well with the velocity of 380 m/s inferred from the magnetic field signature of the electrojet current. This study suggests that such an electrojet is present whenever a polar arc is set up. However, a few conditions are required to observe clear ground magnetic signatures of the electrojet: (1) quiet local magnetic conditions; (2) a single sun aligned arc near or moving across zenith; and (3) arc excitation by ≥1 keV electrons.     

Local time and longitude dependence of the equatorial electrojet magnetic effects

An empirical model of the equatorial electrojet (EEJ), including local time and longitude dependence, has been constructed based on the surface magnetic data recorded at 26 stations located in six different longitude sectors that were set up or augmented during the international equatorial electrojet year (IEEY). The model reproduces the characteristic signatures of the EEJ-associated horizontal and vertical magnetic components at ground level. The model-predicted variations at the orbit of the POGS satellite are generally in good agreement with the onboard magnetic signatures, although strong discrepancies are also often seen. The nature of the differences suggests that the global scale magnetospheric or field-aligned current systems may sometimes dominate the satellite data. The nature of the longitudinal inequalities in the EEJ strength indicates that the equatorial electrojet is strongest in South America (80°–100°W) and weakest in the Indian sector (75°E) with a secondary minimum and a maximum centered, respectively, in the Atlantic Ocean (30°W) and in western Africa (10°E). The EEJ strength is shown to be inversely correlated with the main field intensity along the dip-equator.