Analysis and interpretation of anomalous conductivity and magnetic permeability effects in time domain electromagnetic data: Part II: Sμ-inversion

In the paper by Pavlov and Zhdanov [J. Appl. Geophys. (2001)], we demonstrated that anomalous magnetic permeability of an ore body could result in measurable anomalous effects in TDEM data, which can be used in geophysical exploration. In the current paper, we develop a new method of TDEM data interpretation, which allows simultaneous reconstruction of both electrical and magnetic properties of the rocks. The method is based on a generalization of the S-inversion technique [66th Ann. Int. Mtg., Soc. Expl. Geophys., Expanded Abstr. (1996) 1306; 68th Ann. Int. Mtg., Soc. Expl. Geophys., Expanded Abstr. (1998) 473; J. Appl. Geophys. (1999)] for models containing thin sheets with anomalous conductivity and anomalous magnetic permeability. We call this method Sμ-inversion. Numerical 3-D modeling results demonstrate that Sμ-inversion provides useful information about both subsurface conductivity and magnetic permeability distributions. It can be used as a new tool for imaging TDEM data in mineral exploration.     

Influence of dispersive soil electromagnetic properties on hand‐held time domain electromagnetic sensors

Although metal detectors remain the workhorses of humanitarian demining, it is well established that the performance of both continuous wave (frequency domain) and pulsed induction (time domain) detectors can be severely compromised by so‐called ‘soil‐effects’. Generally, problem soils reduce the signal‐to‐noise ratio and increase the false‐detection rate. In certain locations, the soil‐effect is so severe as to render the detector practically inoperable. The current study is part of an ongoing international effort to establish and quantify the influence of soil electromagnetic properties on the operation of metal detectors and related sensor technologies. In particular, we examine the relative influence of soil electrical conductivity, magnetic susceptibility and associated frequency dependence on the time domain electromagnetic (TDEM) response of pulsed induction metal detectors and related small‐scale TDEM sensors.     

Magnetotelluric profiling of vertical conductors

The presence of lateral contrasts of electrical conductivity modifies the original pattern of electromagnetic fields radiated from remote sources. A magnetic transverse plane wave field, interacting with a vertical conductive and outcropping dike placed between two quarter-spaces of unequal electrical conductivities, creates an anomalous vertical component of the magnetic field. This anomalous field has been analysed by computation, and drafting of master curves. Two case histories are presented to illustrate the application and the effectiveness of the solution. It is concluded that: (i) the response is higher for intermediate values of the conducting body induction number; (ii) the curves can be used for the inter-pretation of magnetotelluric, AFMAG, and VLF exploration data; (iii) it is necessary to develop solutions taking into account the vertical as well as the lateral variation of conductivity.     

Synthetic radargrams from electrical conductivity and magnetic permeability variations

The objective of this work is to assess the importance of electrical conductivity and magnetic permeability variations in ground penetrating radar (GPR) reflections commonly interpreted only in terms of permittivity variations. We use the matrix propagator approach to obtain the surface electric field associated with a horizontally layered model whose three electromagnetic properties vary from layer to layer. The solution is based on the plane wave boundary value problem using inverse Fourier transformation to accommodate particular GPR pulses. Our results indicate that while magnetic permeability is unimportant, reflections from electrical conductivity variations can be of the same order as those associated with electrical permittivity boundaries. In particular, we show that a realistic ground model composed of thin conductive layers can produce radargrams similar to those caused by a lossless permittivity contrast.     

Geophysical surveys for identifying saline groundwater in the semi‐arid region of the central Altiplano,Bolivia

In the central part of the Bolivian Altiplano, the shallow groundwater presents electrical conductivities ranging from 0·1 to 20 mS/cm. In order to study the origin of this salinity pattern, a good knowledge is required of the geometry of the aquifer at depth. In this study, geophysics has been used to complement the sparse data available from drill holes. One hundred time‐domain electromagnetic (TDEM) soundings were carried out over an area of 1750 km². About 20 geological logs were available close to some of the TDEM soundings. Three intermediate results were obtained from the combined data: (i) the relationship between the electrical conductivity of the groundwater and the formation resistivity, (ii) geoelectrical cross‐sections and (iii) geoelectrical maps at various depths. The limited data set shows a relationship between resistivity and the nature of the rock. From the cross‐sections, a conductive substratum with a resistivity of less than 1 Ω·m was identified at most of the sites at depths ranging from 50 to 350 m. This substratum could be a clay‐rich formation containing brines. Using derived relationships, maps of the nature of the formation (sandy, intermediate and clayey sediments) were established at depths of 10 and 50 m. Discrimination between sand and clays was impossible where groundwater conductivity is high (>3 mS/cm). In the central part of the area, where the groundwater conductivity is low, sandy sediments are likely to be present from the surface to a depth of more than 200 m. Clayey sediments are more likely to be present in the south‐east and probably constitute a hydraulic barrier to groundwater flow. In conclusion, the study demonstrates the efficiency of the TDEM sounding method to map conductive zones. Copyright © 2001 John Wiley & Sons, Ltd.     

全张量磁梯度数据的斜导数特征值边界识别方法研究

本文提出了一种全张量磁梯度数据的斜导数边界识别方法,该方法利用全张量磁梯度数据定义了北向斜导数、东向斜导数和垂向斜导数.模型试验表明,北向和东向斜导数无法有效识别磁性异常体边界,异常形态复杂虚假异常较多,而垂向斜导数在高纬度斜磁化或者垂直磁化条件下可以清晰而准确得识别地质体的边界,具有一定的深度识别能力.与现有磁张量梯度的边界识别方法相比,识别效果好,可以有效均衡不同深度地质体的响应.将上述方法应用于大兴安岭地区实测航磁三分量数据转换得到的全张量数据,获得了研究区浅部磁性体的分布,并与三分量数据进行联合解释,取得良好的地质效果.     

井中磁源瞬变电磁响应特征研究

井中瞬变电磁波勘探是一个全空间地球物理场问题.采用Gaver-Stehfest逆拉氏变换方法,正演计算了瞬变信号激励下接收线圈上的电磁场响应.分析了包含井眼泥浆、套管、水泥环和地层的轴对称多层介质模型的电磁场响应特征,考察了各层介质参数对井中瞬变电磁响应的影响.不同电导率井眼泥浆的电磁场响应衰减曲线表明,井眼泥浆电导率...     

Global electromagnetic induction in the Moon and planets

A summary of experiments and analyses concerning electromagnetic induction in the Moon and other extraterrestrial bodies is presented. Magnetic step-transient measurements made on the lunar dark side show the eddy current response to be the dominant induction mode of the Moon. Analysis of the poloidal field decay of the eddy currents has yielded a range of monotonic conductivity profiles for the lunar interior: the conductivity rises from 3·10^?4 mho/m at a depth of 170 km to 10^?2 mho/m at 1000 km depth. The static magnetization field induction has been measured and the whole-Moon relative magnetic permeability has been calculated to be $\text{μ}\text{μ}{}_0\text{= 1.01 ± 0.06}$. The remanent magnetic fields, measured at Apollo landing sites, range from 3 to 327 γ. Simultaneous magnetometer and solar wind spectrometer measurements show that the 38-γ remanent field at the Apollo 12 site is compressed to 54 γ by a solar wind pressure increase of 7·10^?8 dyn/cm². The solar wind confines the induced lunar poloidal field; the field is compressed to the surface on the lunar subsolar side and extends out into a cylindrical cavity on the lunar antisolar side. This solar wind confinement is modeled in the laboratory by a magnetic dipole enclosed in a superconducting lead cylinder; results show that the induced poloidal field geometry is modified in a manner similar to that measured on the Moon. Induction concepts developed for the Moon are extended to estimate the electromagnetic response of other bodies in the solar system.     

Coast effect distortion of marine magnetotelluric data: Insights from a pilot study offshore northeastern Japan

We report on strong coast effect distortions observed for broadband marine magnetotelluric (MT) data collected on the forearc offshore northeastern Japan. Eight days of horizontal electric and magnetic fields recorded at eight seafloor stations and the horizontal magnetic fields from a land remote station were processed with a robust multiple-station algorithm, yielding good MT responses and inter-station transfer functions at periods of 7–10,000 s. Transverse electric (TE) mode responses have cusps in apparent resistivity and negative phases at periods around 1000 s, while the transverse magnetic (TM) mode responses are galvanically depressed below the TE responses. An analysis of inter-station transfer functions confirms that the apparent resistivity cusps are a magnetic field, rather than electric field, phenomenon, consisting of an amplitude minimum and rapid phase change around a characteristic frequency. Poynting vectors for a TE coast effect model study illustrate that the anomalous phases are associated with energy diffusing back up to the seafloor from below, after being turned around from its usual downward propagating trajectory by inductive coupling between the conductive ocean and the resistive seafloor along the continental margin. We show that the characteristic frequency and position of the TE mode apparent resistivity cusps are determined by a relatively simple combination of the electrical resistivity of the seafloor, the depth of the ocean, and the distance from the coastline. By including coastlines and bathymetry in 2D inversion, we recover the seafloor conductivity structure along the forearc, demonstrating that broadband data can constrain the thickness of conductive forearc sediments and the underlying high resistivity associated with the mantle wedge and subducting oceanic lithosphere.     

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

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

Fast approximate EM induction modeling of metallic and UXO targets using a permeable prism

The time-domain EM induction response of non-magnetic and magnetic targets can be approximated using a conductive permeable prism composed of six faces of conductive plates, each face being composed of a set of conductive ribbons. The effect of magnetic permeability is included by the use of two “apparent flux gathering” coefficients, and two “effective magnetic permeability” coefficients, in the axial and transverse directions. These four magnetic property coefficients are a function of physical properties and geometry of the target, but are independent of prism orientation relative to a transmitter. The approximation algorithm is computationally fast, allowing inversions for target parameters to be achieved in seconds. The model is tested on profiles acquired with a Geonics EM63 time-domain EM metal detector over a non-magnetic copper pipe target, and a steel artillery shell in horizontal and vertical orientations. Results show that this approximation to a permeable prism has a capability of fitting geometric, conductivity and magnetic parameters at both early and late sample times. The magnetic parameters show strong change from early to late times on the EMI decay curve, indicating that the magnetic properties of the target have non-linear characteristics. It is proposed that these magnetic parameters and the nature of their non-linearity may carry additional discrimination information for distinguishing between intact munitions and scrap in UXO studies.     

核磁共振与瞬变电磁三维联合解释方法

传统核磁共振地下含水量解释多采用基于均匀半空间或层状导电模型的一维反演,分层给出地下含水信息.然而,这些方法忽略了地下复杂电阻率分布信息对结果的影响,也不能很好地反映局部三维含水构造.本文从三维电介质中核磁共振响应的正演理论出发,提出首先利用瞬变电磁数据进行基于等效导电平面法的快速电阻率成像,然后将成像结果作为核磁共振三维反演的电性模型,进行联合解释.激发磁场的分布采用有限元法直接求解,通过引入伪δ源实现电流源的加载,并强加散度条件排除了三维磁场模拟中"弱解"的影响.针对核磁共振灵敏度矩阵的病态性和数据中存在的干扰信号,提出考虑罚项的非线性拟合目标函数,利用线性化方法进行核磁共振反演.模型数据表明该方法能较准确反映地下三维含水构造,实测算例进一步证明了方法的有效性.本研究将促使核磁共振方法在岩溶、裂隙水、孤立水体等复杂水文地质条件及隧道、矿井灾害水源探测等方面得到有效应用.     

JOINT INTERPRETATION OF GRAVITY AND MAGNETIC DATA OVER AXIAL SYMMETRIC BODIES WITH APPLICATION TO THE DARNLEY BAY ANOMALY,NWT CANADA*

A method is presented for determining bounds of the properties of axial symmetric bodies from a finite number of gravity and magnetic observations based on Parker's theory of ideal bodies. Bounds on the density contrast and the intensity of magnetization are calculated as a function of depth to the top of the anomalous source, restricting the range of smallest possible solutions to fit the data. The model studied is approximated by an array of vertical annuli cylinders, each of uniform density and magnetization. Linear programming algorithms based on the ideal body theory were used to calculate the distribution of these parameters within the body. Simultaneous inversion of gravity and magnetic data is performed assuming a constant ratio between the density contrast and the intensity of magnetization and that a common body is responsible for both observed fields. The parameter k(|J|/δp) provides information about the rock type of the structure. Interpretation of gravity and aeromagnetic data from Darnley Bay, NWT, Canada, indicated the presence of a shallow ultrabasic intrusion.     

Fast Imaging of TDEM data based on S-inversion

Fast S-inversion is a method of interpretation of time-domain electromagnetic (TDEM) sounding data using the thin sheet model approach. Within the framework of this method, the electromagnetic (EM) response measured at the surface of the earth at every time delay is matched with that of a thin sheet model. The conductivity change with depth is obtained using the conductance, S, and depth, d, of the equivalent thin sheet. We analyze two different numerical techniques, the differential S-transformation and the regularized S-inversion, to determine the parameters of the thin sheet. The first technique is a direct differential transformation of the observed data into conductance and depth values. It is fast and requires no iterations or starting model. The second technique uses a regularized inversion scheme to fit the measured response with that of a thin sheet. In both techniques, the retrieved conductance values are differentiated with respect to depth to obtain the conductivity change with depth. We apply S-inversion to three-dimensional synthetic data and we successfully locate the local conductors. We also demonstrate a case history by interpreting TDEM data obtained at the Nojima fault zone in Japan. The results clearly indicate the location of the fault zone.     

Changes in geophysical properties caused by fluid injection into porous rocks: analytical models

Analytical models are provided that describe how the elastic compliance, electrical conductivity, and fluid‐flow permeability of rocks depend on stress and fluid pressure. In order to explain published laboratory data on how seismic velocities and electrical conductivity vary in sandstones and granites, the models require a population of cracks to be present in a possibly porous host phase. The central objective is to obtain a consistent mean‐field analytical model that shows how each modeled rock property depends on the nature of the crack population. The crack populations are described by a crack density, a probability distribution for the crack apertures and radii, and the averaged orientation of the cracks. The possibly anisotropic nature of the elasticity, conductivity, and permeability tensors is allowed for; however, only the isotropic limit is used when comparing to laboratory data. For the transport properties of conductivity and permeability, the percolation effect of the crack population linking up to form a connected path across a sample is modeled. However, this effect is important only in crystalline rock where the host phase has very small conductivity and permeability. In general, the importance of the crack population to the transport properties increases as the host phase becomes less conductive and less permeable.     

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.     

大地电磁相位超象限现象的各向异性模型研究——以上下结构为例

大地电磁观测数据中的相位超象限现象可以由不同的电性结构产生.本文在已实现大地电磁三维任意各向异性有限差分正演的基础上,以具有上下结构关系的三维各向异性模型为例,分析各向异性体的规模以及参数变化对阻抗相位的影响.在下部各向异性体规模明显大于上部各向异性体或表现为层状特征的情况下,上部各向异性体在两个水平方向上的尺度差异较大,可以看作准二维体时容易发生相位超象限;当上部各向异性体在两个水平方向上尺度相近表现为规则三维体时,要产生相位超象限的现象则需要各向异性体具有更高的各向异性比.在同等条件下,增加各向异性体的各向异性比更容易发生相位超象限;而各向异性走向方位角的变化将直接影响到发生相位超象限的范围.对于准二维模型引起相位超象限的条件,沿用二维模型的近似解析分析方法,进一步构建了基于各向异性体的电导率、背景电导率以及各向异性走向角的相位超象限指标函数,从而更加直观地解释在二维或准二维条件下发生相位超象限现象的模型参数特征.     

On removing the primary field from fixed-wing time-domain airborne electromagnetic data: some consequences for quantitative modelling, estimating bird position and detecting perfect conductors

In the process of removing the primary field from fixed‐wing time‐domain airborne EM data, the response is decomposed into two parts, which are referred to here as the time‐domain ‘in‐phase’ and ‘quadrature’ components. The time‐domain in‐phase component is dominated by the primary field, which varies significantly as the transmitter–receiver separation changes. The time‐domain quadrature component comes solely from the secondary response associated with currents induced in the ground and this is the component that has traditionally been used in the interpretation of data from fixed‐wing towed‐bird time‐domain EM systems. In the off‐time, the quadrature response is very similar to the total secondary response. However, there are large differences in the on‐time and even some small differences in the off‐time.One consequence of these differences is that when airborne EM data are to be interpreted using a synthetic mathematical model, the synthetic data calculated should also be the quadrature component. A second consequence relates to the time‐domain in‐phase component which is sometimes used to estimate the receiver‐sensor (bird) position. The bird‐position estimation process assumes there is no secondary field in the in‐phase component. If the ground is resistive, the secondary contained in the in‐phase component is small, so the bird‐position estimate is accurate to about 30 cm, but in highly conductive areas the secondary contribution can be large and the position estimate can be out by as much as 5 m. A third consequence arises for highly conductive bodies, the response of which is predominantly in‐phase. This means that any response from these types of body is lost in the component that has been removed in the primary‐field extraction procedure. However, if the bird position is measured very accurately, the actual free‐space primary field can be estimated. If this is then subtracted from the estimated primary (actually free‐space primary plus secondary in‐phase response), then the residual is the secondary in‐phase response of the ground. Using this methodology, very conductive ore bodies could be detected. However, a sensitivity analysis shows that detection of a large vertically dipping very conductive body at 150 m depth would require that the bird position be measured to an accuracy of about 1.4 cm and the aircraft attitude to within about 0.01°. Such tolerances are very stringent and not easily attainable with current technology.     

磁激发极化法中水平圆柱体的一些研究结果

本文简要介绍导电柱体在全空间和半空间的均匀电流场中形成的异常电流磁场理论,给出了正、反演理论公式;计算出了磁电阻率法和磁激发极化法理论剖面曲线;得到了相应的模型实验结果以及在浸染型铜矿上的野外试验结果。这些资料证明了所建立的理论。     

Satellite and surface geophysical expression of anomalous crustal structure in Kentucky and Tennessee

An equivalent layer magnetization model obtained from inversion of long-wavelength satellite magnetic anomaly data indicates a very magnetic source region centered in south central Kentucky. The magnetization maximum nearly coincides with a gravity high elongated north-south and extending into Tennessee. Previous refraction profiles suggest that the source of the gravity anomaly is a large mass of rock occupying much of the crustal thickness. The outline of the source delineated by gravity contours is also discernible in aeromagnetic anomaly patterns. Taken together, the geophysical data suggest a large, localized mass of intracrustal rock which is both dense and very magnetic. A simple magnetization/density model is given which accounts for the gravity and long-wavelength aeromagnetic anomalies due to the body. We interpret it as a mafic plutonic complex, and several lines of evidence are consistent with a rift association. The body is, however, clearly related to the inferred position of the Grenville Front. It is bounded on the north by the fault zones of the 38th Parallel Lineament. The inferred mean magnetization (4 A/m) of the body is large, but not inconsistent with values reported by others for deep crustal bodies associated with long-wavelength magnetic anomalies. Such magnetization levels can be achieved with magnetic mineralogies produced by normal oxidation and metamorphic processes and enhanced by viscous build-up, especially in mafic rocks of alkaline character.