位于导电介质中三维导体的电磁响应

本文给出一种用积分方程法配合有限差分法计算位于导电介质中三维导体的电磁响应的新途径,在二层大地条件下给出了具体计算方法和结果。详细讨论了围岩介质的导电性对导体异常的影响。这对于开展低阻覆盖层地区的电磁法工作具有一定的指导意义。     

位于导电介质中三维导体的电磁响应

本文给出一种用积分方程法配合有限差分法计算位于导电介质中三维导体的电磁响应的新途径,在二层大地条件下给出了具体计算方法和结果。详细讨论了围岩介质的导电性对导体异常的影响。这对于开展低阻覆盖层地区的电磁法工作具有一定的指导意义。     

Two field applications of the rotating current EM method1

The rotating current EM method has been applied to the delineation of two conductive orebodies, Elura near Cobar, NSW, and Thalanga near Charter's Towers, Queensland. The field data were collected in the form of observations of the vertical magnetic field strength ratio and phase difference using a Turam-style receiver with twin vertical coils. By reconstituting this data back to the ring source field and phase, i.e. the observed H_z, phasor, it is possible to present contoured maps of the EM field. Anomaly phasors are obtained by subtracting theoretical phasors from the observed phasors in the complex plane of the H_z phasor. The theoretical phasors for the finite source are based on horizontally layered, half-space earth models, computed at each point of the survey grids, then normalized to selected points of the observed fields. Use is made of the intrinsic circular symmetry of the method in X–Y plots of field versus source-receiver distance to ascertain geoelectric parameters for the earth models. A steel picket fence at Thalanga is modelled by a line source grounded at each end and its H_z, phasor is removed by the same process. A considerable improvement in anomaly delineation is gained over previous Turam-style anomalies and the two survey examples illustrate the limitations of the method in the presence of a conductive overburden (Elura) and its abilities in the absence of a conductive overburden (Thalanga).     

Theoretical central frequency sounding response curves over a generalized three-layer earth

The digital linear filter method is used to compute the normalized vertical magnetic field for a circular loop in CFS system. Three-layer earth models with resistive and conductive basement are considered. The corresponding field expressions are suitably written, and the multifrequency response is computed and presented in convenient forms. Analysis of theoretical data indicates that for highly resistive basement, the variation in layer conductivity and intermediate layer thickness is well reflected on three-layer amplitude response curves at low frequencies and at high conductivity contrasts between first and second layers. This, however, is not true in the case of conductive basement, where the resolution of the intermediate layer is observed to be comparatively poor. The resolution of an intermediate conductive layer in a three-layer sequence is found to be satisfactory.     

Comparison of Theoretical and Physical Model Studies of the Responses of Moving Source and Fixed Loop Electromagnetic Exploration Systems

Theoretical and physical scale model responses for moving source and fixed loop electromagnetic exploration systems are presented for the case of a vertical, thin, plate conductor located first in free-space and then in a conductive host. The results are presented in the form of anomaly index diagrams plotted in complex space. The two modeling approaches show general agreement on the behavior of the current gathering effect that operates when the conductive host is present. Both approaches show an onset of a strong positive effect at higher frequencies in the responses provided by both horizontal coplanar coil systems and by fixed loop systems. Agreement is also shown on an unanticipated negative effect that a conductive host causes in the responses provided by the vertical coplanar coincident coil system. Both modeling approaches demonstrate that the responses provided by fixed loop systems are predominantly caused by the current gathering effect. The extent of this effect for fixed loop systems depends on the size of the transmitter loop with respect to the target and on the depth of the target. The theoretical modeling demonstrates that both moving source and fixed loop systems operated over a target located in a conductive host, provide responses that become almost identical and independent of the conductance of the target at high frequency.     

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.     

A model for the coast-effect

A model for the coast-effect of geomagnetism is presented, in which the horizontal magnetic field induces currents in a circuit including a thin finite ocean. The currents flow horizontally across the ocean, vertically down into the earth, back through the deep interior of the earth, and vertically up to the ocean to complete the current loop. The upper layers of the earth are given non-zero conductivity, allowing the possibility of such current loops.A two-dimensional model involving such currents has been worked out analytically, and it is found that a significant induced magnetic field at the seafloor can be obtained with a reasonable conductivity in the earth's upper layers. A three-dimensional model has also been worked out numerically. It is found that the induced vertical component of magnetic field is of comparable magnitude to the horizontal component induced normal to the coast, whereas the horizontal component parallel to the coast is small. These relations are required to explain the observation of Parkinson arrows.     

基于印模法的地面磁电阻率数据三维非线性共轭梯度反演

本文提出了能提高异常体分辨能力,同时得到绝对电导率的地面磁电阻率数据三维反演方法.磁电阻率响应用准直流的低频磁场代替;数值模拟由频率域电场满足的Helmholtz方程出发,采用三维交错网格有限差分法;长直导线源作为发射源,其中源的计算包含在背景场中;结合地面磁电阻率数据各分量的特点,选择y分量进行反演研究;反演采用三维非线性共轭梯度反演技术,为了提高异常体的深度分辨能力,进行迭代重构反演;用印模法对初始模型进行重构,采用的是辅模型在浅部,元模型在深部的组合方式.从合成数据和实际数据的反演结果可以得到以下的认识:(1)由频率域麦克斯韦方程组出发,低频磁场数据反演可以直接得到电导率,而不是相对电导率之比;(2)采用印模法组合初始模型,进行迭代重构反演,可以提高地面磁电阻率数据反演对异常体的分辨能力,确定埋深位置,同时不会丧失对于浅部异常体的分辨能力;(3)在结合印模法的地面磁电阻率数据三维反演中,深部异常体的分辨能力受地表不均匀导电体影响较小;(4)确定印模深度可以采用上一次重构反演结束时的模型变化量,通过相邻两次重构反演结束时的模型变化量之差来确定迭代重构是否终止.因为静磁场与重力场在数学上的相似性,本文的反演方法可以被运用到重力场等位场的地面数据的反演中.     

三维地形频率域人工源电磁场的边界元模拟方法

提出了一种用边界元法计算频率域人工源三维地形电磁场的数值模拟方法.首先用矢量积分理论和电磁场边界条件，将上半空间（空气）和下半空间（地下介质）两个区域电磁场边值问题变为仅对地形界面的两个矢量面积分方程.然后将对地形界面的积分剖分为一系列的三角单元积分.在三角单元积分中，假设单元中电磁场为无限大气空间电磁场与地形影响的叠加,并假设地形影响为常项,这样既保证了计算精度又使得计算方法简便.通过分解和计算，每一个矢量面积分方程分解为对应三个坐标方向的三个常量线性方程，这些线性方程组成了对角占优的线性方程组，可用SSOR方法求解.文中给出了垂直磁偶源的垂直磁场地形影响的例子.     

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.     

Relative performances of central and dipole frequency soundings over a layered earth

Variable frequency soundings in the audio-range replaces shallow conventional direct current methods for determination of layer parameters when surface layer resistivity is high. Central frequency soundings (CFS) is one such method that involves measurement of the existing vertical magnetic field component induced at the centre of a horizontal circular or square loop. Dipole method of frequency sounding using small horizontal coplanar loops (abbreviated DFS) measuring the same field component is also considered. Theoretical studies on CFS and DFS over two- and three-layer horizontally stratified earth are carried out and the response characteristics computed and analysed.Theoretical response curves for CFS and DFS over two- and three-layer earth models are presented in convenient forms. Response curves under similar geological-physical conditions are compared. The study indicates that the relative superiority of a method is controlled largely by the nature of the conductivity contrast and the ratio of the first layer thickness to loop radius or dipole separation. While CFS shows a better resolution of conductivity contrast, DFS works better at high frequencies. For resistive substratum, however, both lack proper resolution.     

导电大地上航空电磁系统校准

航空电磁系统校准是开展实际测量工作的基础,校准情况直接影响数据处理和解释.传统校准方法通常假设在自由空间中进行,忽略导电大地耦合影响.然而,实际工作中很难找到绝对高阻的校准场地,导电大地对系统校准和观测数据的影响无法忽视.本文以频率域航空电磁系统为例,对导电大地上航电系统校准技术和校准误差改正方法进行研究.我们首先推导了层状导电大地上水平共面和直立共轴线圈系统的校准公式,结果表明导电大地对航电系统校准尤其是水平共面装置的高频信号影响很大.针对校准过程中大地电导率已知的情况,本文采用非线性方程求解技术一次性确定校准线圈位置和Q值;在没有任何辅助信息情况下,也可直接利用实测数据计算校正因子进行迭代求解.测试结果表明该方法快速、准确、有效.考虑到系统相位和增益调整直接影响观测数据,本文提出了航空电磁数据校准误差的改正算法.实测数据误差改正结果表明,导电大地对高频信号影响严重,校准误差改正后的航空电磁数据与实际地质资料更好吻合.     

基于水文-地球物理模型的地下水污染磁电阻率异常动态监测仿真

We present a forward-modeling investigation of time-dependent ground magnetometric resistivity （MMR） anomalies associated with transient leachate transport in groundwater systems. Numerical geo-electrical models are constructed based on the hydrological simulation results of leachate plumes from a highly conceptualized landfill system and the resultant MMR responses are computed using a modified finite difference software MMR2DFD. Three transmitter configurations （i.e., single source, MMR-TE, and MMR-TM modes） and two hydrological models （i.e., uniform and faulted porous media） are considered. Our forward modeling results for the uniform porous medium indicates that the magnetic field components perpendicular to the dominant current flow contain the most information of the underground targets and the MMR-TE mode is an appropriate configuration for detecting contaminant plumes. The modeling experiments for the faulted porous medium also confirm that the MMR method is capable of mapping and monitoring the extent of contaminant plumes in aroundwater systems.     

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.     

Direct current electric potential computation in an inhomogeneous and arbitrarily anisotropic layered earthl

We consider the calculation of the electrical field quantities, electric potential and the vertical component of the total volume density of electric current, in a horizontally layered, piecewise homogeneous and arbitrarily anisotropic earth due to a system of direct current point sources. By applying Fourier transformation with respect to the horizontal space coordinates to the static field equations, the field quantities are obtained as the solutions of the system of transform-domain differential equations in the vertical (depth) coordinates. A recurrence scheme has been given to compute the tranform-domain field quantities at any depth. The corresponding space-domain quantities are then obtained by inverse Fast Fourier Transformation (FFT). A complete computer program has been developed for computing the electric potentials at any depth of the layered earth, which is composed of an arbitrary number of anisotropic layers with arbitrary conductivity tensors. By considering the point sources at different depths from the surface, equipotential contours on the surface of arbitrarily anisotropic layered earth models are given.     

Non-local effects in the mean-field disc dynamo I. An asymptotic expansion

Abstract

We reconsider thin-disc global asymptotics for kinematic, axisymmetric mean-field dynamos with vacuum boundary conditions. Non-local terms arising from a small but finite radial field component at the disc surface are consistently taken into account for quadrupole modes. As in earlier approaches, the solution splits into a local part describing the field distribution along the vertical direction and a radial part describing the radial (global) variation of the eigenfunction. However, the radial part of the eigenfunction is now governed by an integro-differential equation whose kernel has a weak (logarithmic) singularity. The integral term arises from non-local interactions of magnetic fields at different radii through vacuum outside the disc. The non-local interaction can have a stronger effect on the solution than the local radial diffusion in a thin disc, however the effect of the integral term is still qualitatively similar to magnetic diffusion.     

A physical approach to computing magnetic fields1

The superposition integral expressing the field due to a magnetic source body is relatively simple to evaluate in the case of a homogeneous magnetization. In practice this generally requires that any remnant component is uniform and the susceptibility of the body is sufficiently low to permit the assumption of a uniform induced magnetization. Under these conditions the anomalous magnetic field due to a polyhedral body can be represented in an intuitive and physically appealing manner. It is demonstrated that the components of the magnetic field H can be expressed as a simple combination of the potentials due to two elementary source distributions. These are, firstly, a uniform double layer (normally directed dipole moment density) located on the planar polygonal faces of the body and, secondly, a uniform line source located along its edges. In practice both of these potentials (and thus the required magnetic field components) are easily computed. The technique is applicable to polyhedra with arbitrarily shaped faces and the relevant expressions for the magnetic field components are suitable for numerical evaluation everywhere except along the edges of the body where they display a logarithmic singularity.     

Fast imaging of subsurface conductors using very low‐frequency electromagnetic data

The study presents a fast imaging technique for the very low‐frequency data interpretation. First, an analytical expression was derived to compute the vertical component of the magnetic field at any point on the Earth's surface for a given current density distribution in a rectangular block on the subsurface. Current density is considered as exponentially decreasing with depth, according to the skin depth rule in a particular block. Subsequently, the vertical component of the magnetic field due to the entire subsurface was computed as the sum of the vertical component of the magnetic field due to an individual block. Since the vertical component of the magnetic field is proportional to the real part of very low‐frequency anomaly, an inversion program was developed for imaging of the subsurface conductors using the real very low‐frequency anomaly in terms of apparent current density distribution in the subsurface. Imaging results from the presented formulation were compared with other imaging techniques in terms of apparent current density and resistivity distribution using a standard numerical forward modelling and inversion technique. Efficacy of the developed approach was demonstrated for the interpretation of synthetic and field very low‐frequency data. The presented imaging technique shows improvement with respect to the filtering approaches in depicting subsurface conductors. Further, results obtained using the presented approach are closer to the results of rigorous resistivity inversion. Since the presented approach uses only the real anomaly, which is not sensitive to very small isolated near‐surface conducting features, it depicts prominent conducting features in the subsurface.     

瞬变电磁场的直接时域数值分析

为了深入了解瞬变电磁场的勘探原理,直接在时间域对负阶跃脉冲激发的二维瞬态场进行了数值分析.采用的方法是从反映电磁场基本规律的麦克斯韦方程组出发,导出时域电场的齐次扩散方程,对所研究的空间区域作差分离散,源作为初始条件加入,利用准静态近似处理空中边界,然后进行时间的逐步递推,由此展现瞬变电磁场在地下扩散随时间发展的全过程.通过模拟计算不同时刻瞬态电场在地下的分布形态及地面上感生电动势相应的变化,揭示了低阻异常体对感应涡流的聚集作用,低阻覆盖层对瞬变场扩散的减速作用,及瞬变场的延时效应.因此,瞬变电磁法对低阻体是敏感的,有上覆低阻层时探测同样的深度需要较长的时间,而延时效应瞬变场的晚期时段可反映埋藏较浅的异常体.     

Stationary approximation of the magnetotelluric field in a halfspace with a vertical boundary and three-dimensional perturbing body

Summary Based on the generalized potential of a double layer, integral formulae have been derived for calculating the stationary approximation of the magnetotelluric field in a halfspace divided by a vertical boundary into two quarterspaces, one of which contains a three-dimensional perturbing body. The appropriate boundary integral equation and other surface integrals have been computed for a perturbing body in the shape of the three-dimensional prism located at the vertical boundary, or in contact with the said boundary. The exciting electrical field is assumed to be homogeneous and perpendicular to the vertical boundary. Isoline graphs of the electrical and magnetic fields on the surface of the halfspace have been plotted and their anomalies are discussed.