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.     

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.     

一维电阻率各向异性对海洋可控源电磁响应的影响研究

地下介质的电阻率常常表现为各向异性,海底裂隙地层和层状沉积序列可能形成宏观电阻率各向异性.在解释海洋电磁资料时,电阻率各向异性的影响不应该被忽略,否则可能会得到错误的海底地电模型.作者编写了电阻率任意各向异性一维层状介质海洋可控源电磁场计算程序,计算了电阻率各向异性层状模型的海洋可控源电磁响应,讨论了覆盖层和高阻储层分别具有电阻率各向异性时的电磁场响应特征.     

强震前地下水活动引起的垂线变化

本文分析了引起强震前天文测时测纬异常现象的原因问题,认为这种异常主要是地下水活动引起垂线变化的结果.垂线变化包括重力场水平分量变化和垂直分量变化,前者可引起测时测纬异常.地下水活动的主要形式是垂直迁移,这种垂直迁移又沿着一定的水平方向传递.从测时测纬的异常形态和地下水的活动规律,可以分析出应力场变化的某些规律.地下水的这种活动,引起的重力场水平分量变化远远大于垂直分量的变化,两者的变化形式也很不同.      

Biaxial anisotropy in geoelectric prospecting

In applied studies on electromagnetic well-logging and electric prospecting, a medium is supposed to be either isotropic or with conventional anisotropy of conductivity. However, in the meantime, there is a clear hint that a medium may exhibit even biaxial anisotropy; i.e., the resistances along all three directions X, and Y, Z are different. Based on the analytical solution by the method of separation of variables, the paper considers an algorithm for the calculation of the electromagnetic field in a layered medium with biaxial conductivity anisotropy involving an arbitrary direction of horizontal conductivities in each layer. In this case, the theoretical solution and the algorithm display substantial peculiarities, and the numerical implementation involves many complexities. These problems were examined and solved mainly due to the constant comparison with the calculations carried by the finite element method and their analysis. Programs for the calculation of electromagnetic fields in the frequency and time domains were developed, and results of geoelectric interest were obtained. For example, the excitement of a horizontally-layered section by a magnetic dipole produces a vertical electric component of the field.     

RESPONSE OF A TWO-LAYER GROUND TO AN INFINITE CABLE CARRYING ALTERNATING CURRENT WITH REFERENCE TO TURAM ANOMALIES *

A formula is developed for the vertical magnetic field due to an alternating current passing through a long horizontal cable placed on a two-layer ground. The variations of the phase and amplitude ratio of the vertical field, along profiles perpendicular to the cable line are investigated. Within the limits encountered in practice, it is found that if the upper layer is conductive, as compared to the lower layer, the phase may vary appreciably whereas the amplitude ratio changes as in the case of vacuum.     

QUANTITATIVE INTERPRETATION OF TRANSIENT ELECTROMAGNETIC FIELDS OVER A LAYERED HALF SPACE*

The transient fields from a finite horizontal loop excited by a half sine wave current pulse have been computed numerically for a particular source receiver configuration at a height of 100 meters above a layered ground. The amplitude of the vertical component of the magnetic field has been chosen for the interpretation. Curves of apparent conductivity vs. time, plotted during the off-time of the signal, show that layering is easily resolved, that resonance effects are present and that polarization effects are detectable for certain types of polarization.     

Variations in the schumann resonance polarization ellipse in the horizontal and vertical planes according to observations at the Barentsburg and Lovozero observatories

The diurnal variations in the amplitude and parameters of polarization ellipses of the first Schumann resonance according to three magnetic field components, observed on December 1–10, 2007, at the Barentsburg and Lovozero observatories, have been studied. Ellipses have been constructed in the (H, D) and (H, Z) planes. The value of the minor axis, inclination of the major axis, ellipticity, and the rotation direction have been estimated. The vertical magnetic component of the Schumann resonance is three to four times as small as the horizontal component. The difference in the diurnal variations in the ellipse parameters between both observatories has been found. The effect can be caused by a difference in the Earth’s conductivity in the vicinity of the observatories and by conductivity anisotropy. The major axis inclination and ellipticity have semidiurnal components. The polarization vector mostly rotates from D to H in the horizontal plane and from Z to H in the vertical plane at both stations.     

SOTEM响应特性分析与最佳观测区域研究

电性源短偏移距瞬变电磁法(SOTEM)是目前研究和应用较为广泛的一种人工源时间域电磁法工作装置,对深部资源地球物理精细探测具有一定的实际意义.为了深入理解方法内涵并更好地进行推广应用,本文基于电性源瞬变电磁一维正演理论,研究了SOTEM地下感应电流扩散、多分量电磁响应平面分布、多偏移距衰减等特性,然后根据上述特性研究了SOTEM的最佳观测区域.研究结果表明:电性源在地下可以产生水平和垂直两个方向的感应电流.其中,水平感应电流又分为上部水平感应电流和下部水平感应电流(又称作返回电流),水平感应电流的极大值主要集中于发射源附近并垂直向下扩散;垂直感应电流极大值沿与地面呈45°角的方向向下、向外扩散,并且具有较低的振幅和较快的扩散速度.电性源激发的六个方向的电磁场分量都具有一定的探测能力,但是考虑到地面观测的方便性和各分量的传播、分布特点,大多数情况仅利用垂直磁场分量Hz(B/t)和水平电场分量Ex.其中,Hz仅对低阻目标体敏感,且敏感区域位于赤道向区域,并集中在发射源附近;Ex既对低阻体敏感也对高阻体敏感,对低阻体的敏感区域位于赤道向区域,而对高阻体的敏感区域位于轴向区域,并且敏感区域距发射源的距离与目标体埋深和围岩电性有关.     

Anomalies in the electric field and electric conductivity of the Earth’s crust in relation to the Kultuk Earthquake on Lake Baikal

We analyze the variations in the horizontal and vertical components of the Earth??s electromagnetic field measured at nine points onshore along the southern part of Lake Baikal. It is found that variations with a period of a few tens of seconds have decayed 20 minutes before the earthquake. This is a regional effect. It is observed in the variations of the horizontal magnetic field. The intensity of the effect in the electric field is site-specific, which is due to the nonuniform distribution of electric conductivity in a medium, as revealed by the deep magnetotelluric sounding. The electric field also experiences a coseismic variation associated with the seismic wave. This effect is stronger pronounced in the vertical electric field. The bay-like anomalies that precede and accompany the Kultuk earthquake are revealed in the electric field of intraterrestrial origin and in the behavior of the electric tipper. These anomalies are concurrent with the rise of the ground water level in the borehole. The probable nature of the revealed anomalies is discussed.     

垂直向地震作用对节理岩体失稳破坏的影响

基于线弹性断裂力学理论分析了垂直向地震作用对节理岩体地震动力破坏的影响。在仅考虑峰值时,最不利的单向地震动加速度方向是水平倾向坡外,双向则依据破裂机制是拉剪或压剪,加速度分别是水平倾向坡外与向下或向上的组合。地震动的幅值、作用方向及双向地震动的组合都可使岩体的破坏机制发生转化,并且是突变的、不可逆的。较低峰值的双向地震动产生的应力强度因子可能大于较高峰值的单向地震动所产生的应力强度应子。在岩体节理分布特征和静态应力场一定的初始条件下,第一个导致岩体中产生破裂的地震动加速度幅值及其方向的组合唯一地决定了岩体不可逆破坏发展的方向、机制及最终的破坏特征,其复杂性远大于静力作用时的情况。对岩体地震动力破坏问题的认识应充分考虑垂直向地震动的重要影响。     

江苏地区地电暴特征及差异性

利用FFT方法对江苏区域地电场、地磁以及有关台站的钻孔应变数据进行处理,提取频谱特征,对该区各地电台站的地电暴特征进行分析,研究结果表明：①江苏区域地电场各测项地电暴数据与地磁H分量变化趋势接近,但各台地电场不同测项的突跳方向存在各向异性,4个地电台站EW向数据具有一致同向现象,NS向数据有差异性;②由各分量优势周期对应性可知,江苏地区地电场NS向含有的地磁场H分量高于体应变成分,地电场EW向倾向于地磁场H分量与体应变叠加抵消的结果;③当磁暴发生时,地磁H分量变化幅度相同的台站,地电场变化幅度不同,地电台站的地电暴变化特征与其地下结构有关,地下深层阻抗越大,感应电场越大,地电场日变化幅度越大。     

电性源地-井瞬变电磁法全分量响应特性与探测能力分析

结合地面电性源瞬变电磁法（TEM）探测深度大和地-井TEM探测精度高的优点,本文提出对电性源地-井TEM进行研究.基于一维正演理论,本文对电性源在地下激发的六个电磁场分量的扩散、分布特性和探测能力进行分析研究.研究结果表明地下瞬变场各分量的分布特征与地面情况类似,但是对地层的探测能力与地面情况存在较大差别.各分量的探测能力与它们的扩散特性有关,其中受积累电荷作用,垂直电场E_z和水平磁场Ḣ_x分量对目标层的反映最为明显,E_x和Ḣ_y分量在地下扩散时受"返回电流"的影响会出现变号现象,使得异常体带来的影响被掩盖,不易被分辨.本文研究结果为发展电性源地-井TEM的施工技术、数据处理与解释,建立了相应的理论基础,获得了对该方法的初步认识.     

CONSTRUCTION OF COMPONENT MAPS FROM AEROMAGNETIC TOTAL FIELD ANOMALY MAPS1

Total field anomalies as defined from normal aeromagnetic surveys give good approximations of the anomalous components along the direction of the main geomagnetic field, which is generally much larger than the anomalous field. Using the relations between vertical and horizontal components of the field, the total field anomaly is related to any vertical or horizontal component and the corresponding horizontal and vertical derivatives. An example from the Siljan impact structure indicates that such directional filters may be applied to extract useful information from magnetic maps.     

Spatial fine structure of model electric pulses in the mesosphere above a Γ-shaped stroke of lightning

The fine structure is discussed arising in the spatial distribution of the pulsed electric field above a Γ-shaped stroke of lightning. The channel of discharge contains a vertical and a horizontal section. The structured spatial distribution of field appears due to a superposition of three pulses arriving from vertical and horizontal sections of the causative discharge and from reflections from the ground. The details of electric field distribution depend on time and on the stroke orientation relative to an elevated observer. The characteristic size of ‘filaments’ in the transient electric field is about 1 km along the horizontal direction, while it reaches a few tens of kilometers along the vertical direction.     

Balloon measurement of vertical and horizontal atmospheric electric fields

Summary Twenty-eight balloons have been flown in Canada to measure vertical and horizontal electric fields at balloon altitudes. A horizontal electric field of magnitude typically between 10 and 50 millivolts/meter has been measured at night at high latitudes in association with auroras and magnetic disturbances. Its origin is in the magnetosphere and it maps to balloon alitudes with small attenuation according to Maxwell's equations. The vertical electric field at the balloon displays variations as the horizontal ionospheric field changes in order to maintain . Thus, magnetospheric processes affect both vertical and horizontal atmospheric electric fields and the potential differences induced by these processes may be comparable to weather induced potential differences. Weather processes have also been observed to produce large horizontal electric fields at balloon altitudes. Methods of distinguishing horizontal fields of ionospheric origin from those of weather origin are investigated with the conclusion that a determination of the source of a given event can often if not generally be made.This paper was presented byU. Fahleson     

ELECTROMAGNETIC RESPONSE OF TWO-DIMENSIONAL INHOMOGENEITIES IN A DISSIPATIVE HALF-SPACE FOR TURAM INTERPRETATION *

An integral equation technique based on the integral formulation of the electric field of a harmonically excited line-source in the presence of a two-dimensional scattering in-homogeneity is outlined. In this treatment the inhomogeneity is assumed to have an arbitrary cross-sectional shape and to be located in a layered, dissipative half-space. The scattering responses in terms of the horizontal and the vertical magnetic field components observed on the ground surface are studied for a variety of geologic models involving overburden layers, vertical or inclined veins and multiple inhomogeneities. The horizontal component of the total magnetic field appears to be the most diagnostic response parameter. A thin vein shaped target with a cross sectional area of 0.05 δ₂× 0.50 δ₂ (δ₂ being the skin-depth in the lower half-space) located 0.50 δ₂ away from the line-source and under a highly conductive overburden layer (σ₁/σ₂= 300) of thickness 0.02 δ₂ is easily detected with a moderate resolution. The scattering responses are also sensitive to changes in depths of burial, inclination and conductivity contrast with the surrounding half-space of the target inhomogeneity. An interpretation scheme as well as a number of characteristic detectability parameters are developed for the conventional Turam method used in mineral exploration.     

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.     

基于站间天然电磁场单位脉冲响应的大地电磁时间序列去噪方法

针对天然大地电磁场信号在人文活动密集地区易受噪声干扰的问题,本文提出利用两个同步测点天然电磁场时间序列之间的单位脉冲响应,合成本地点受干扰时段的数据,从而去除大地电磁噪声.首先,选择高信噪比时段的数据,采用最小二乘法,估算本地点与参考点之间的单位脉冲响应,再根据卷积定律,结合参考磁场合成本地点的磁场和电场.最后用合成数据替换含噪声时段数据,实现时间域去噪.实测高信噪比数据和含噪数据的处理结果表明,该方法可以高精度合成本地点磁场与电场信号,有效去除本地点电场和磁场噪声,包括相关噪声,提高大地电磁数据质量.     

An analogue model study of electromagnetic induction in the eastern coastal region of North America

The behavior of electric and magnetic field variations over the eastern coastal region of North America is studied using a scaled laboratory electromagnetic analogue model. The model source frequency used simulates a period of 1 h in the geophysical scale. The results indicate that deflection and conductive channelling of induced electric current is important for both the E-polarization (northeast-southwest direction of the electric field of the source) and the H-polarization (northwest-southeast) of the source field. In the model, conductive channelling occurs through the Strait of Belle Isle, Cabot Strait, and in the St. Lawrence River. Current deflection is particularly prevalent around the southeast coast of Newfoundland for both E- and H-polarization, and around the northeast coastline of Nova Scotia for E-polarization. The model results also show current deflection by cape and bay coastal features, as well as by ocean depth contours.A comparison of model measurements for the cases of a uniform source field and a line current source indicate that the nature of the source field has a measurable but surprisingly small effect on the vertical to horizontal magnetic field ratio for both E- and H-polarizations, and negligible effect on the magnetotelluric ratio for coastal regions.The model fields in coastal regions were found to be strongly influenced by induced currents, deflected and channelled by the coastline and ocean bathymetry, and were dependent on the nature and particularly the polarization of the source field. Thus, along the complex coastline of eastern North America, a wide range of electric and magnetic field values should be expected. In some regions the coast effect, measured by the vertical to horizontal magnetic field ratio at the coast, could be expected to be extremely small or absent, while in other regions the ratio could approach a value as large as unity for variations of 1 h period.