含激电效应的CSAMT一维正演研究

地电体对频率域电磁波激发源的响应为电磁感应和激电效应的综合响应.传统CSAMT法进行数据正反演时认为大地介质电阻率是与频率无关的实数,而实际上因为激电效应,地下可极化体的电阻率是一个与频率相关的复数.为推进二者总体响应研究,并扩展激电法的应用范围,同时提高电磁法勘探的精度,本文基于Dias模型,以复电阻率代替不考虑地电体极化效应的直流电阻率,对CSAMT场源一维层状模型进行了正演模拟,为提取CSAMT信号中所含激电信息提供理论基础.结果表明,考虑激电参数后,视电阻率及相位响应曲线出现明显异常（包括远场、过渡场、近场）;极化前后振幅比值异常峰值、相位差值异常峰值可直观体现激电异常;异常峰值与极化层层厚、埋深以及电阻率变化有连续的对应关系.认为从频率域电磁法信号中提取激电信息有乐观的前景.     

INDUCED POLARIZATION RESPONSE OF A HORIZONTALLY MULTILAYERED EARTH WITH NO RESISTIVITY CONTRAST*

The induced polarization response of a horizontally multilayered earth with no resistivity contrast can rapidly be calculated on a desk calculator or minicomputer for any electrode array. The formulation is a simple series summation of the products of weighting coefficients and the true induced polarization responses for each of the layers. The coefficients are directly derivable from the corresponding resistivity model. This series approach to IP formulation was originally described by Seigel but has not been treated extensively in the present-day geophysical literature. This method can be applied to either time or frequency domain induced polarization measurements. Once the coefficients are known, apparent induced polarization response can readily be obtained by judicious substitution of known, suspected, or assumed values of the true induced polarization of each layer. Basic formulation is presented for the IP potential coefficients (pole-pole or two array) with no resistivity contrast between the layers. From these coefficients, response of any number of layers for any electrode array can be obtained by suitable differentiation. Some examples of Wenner array for a three-layered earth and dipole-dipole array for a four-layered earth are used to illustrate the application. The results of this technique are valid for many natural situations of modest resistivity contrast. However, they definitely cannot be used if there are highly contrasting resistivity layers present. Such an approach is conceptually simple and is useful for survey planning, checking or setting the “depth-of-penetration”of a given array. For field induced polarization data that fits reasonably well to the no-resistivity-contrast model, this simple approach facilitates quantitative interpretation.     

A Method to Remove Electromagnetic Coupling from Induced Polarization Data for an "Exponential" Earth Model

--The electromagnetic (EM) coupling effect in induced polarization (IP) data is an important problem. In many works it has been computed only considering homogeneous or layered earth models with discretely uniform conductivity. In this study, an algorithm has been developed to compute the EM coupling effect in IP data measured on the earth, whose conductivity varies (increases or decreases) exponentially with depth. The EM coupling effects for Percent Frequency Effect (PFE) and phase data are computed for a dipole-dipole array with different separations, however the method can be applied to any electrode array. The results obtained for the cases of increasing and decreasing conductivity as a function of depth indicate that the EM coupling effect strongly depends on the subsurface resistivity and the dipole length. Here an "exponential" earth model is considered to remove EM coupling from the IP data in frequency and phase domain. For this purpose, first, the region of pseudo-section is divided into segments, and within each segment a typical average apparent resistivity („_a) curve is constructed. An exponential conductivity model is fitted to average „_a data. The conductivity model is then used to compute EM responses. Next the data are corrected for the EM coupling contribution. This decoupling process is applied to field data from a galenite-pyrite mineralization area at the Dolluk site, in western Turkey. The results from the decoupling method developed here are compared with other techniques.     

时间域航空电磁法激电效应对电磁扩散的影响

由于激发极化效应的影响,时间域航空电磁晚期道信号经常会出现变号现象.基于电阻率的传统反演方法无法对变号数据进行正确反演,因此通常在数据处理中予以剔除.为深入了解极化介质的电磁扩散特征,认识航空瞬变电磁负响应的产生机理,本文研究时间域航空电磁系统的电磁扩散特征.我们以均匀极化、非极化半空间及层状介质模型为例,通过直接积分的方法求解频率域电场响应,并由欧姆定律得到电流响应,再经过汉克尔变换得到时间域电流响应.通过研究电流随时间在地下极化介质中的传播特征研究电磁扩散过程;通过对比不同激电参数对电磁扩散的影响,研究极化介质中感应电流与极化电流的扩散规律,从而合理地解释极化介质中负响应的产生机理.基于本文研究和分析结果,可加深对时间域航空电磁法中激电效应的认识.     

编码电磁测深

利用逆重复M序列伪随机信号良好的自相关特性,提出相关检测地电传输特性的编码电磁勘探法.分析了编码电磁测深原理及参数提取方法.同时记录发送电流信号和多收发距电磁场响应,对源信号和场信号进行相关运算,解卷积分离接收系统响应后可得到大地的频率特性或时间特性,以此实现地电断面的精细探测.在时间域,大地冲激响应和阶跃响应含有丰富的地电信息.通过冲激响应的峰值时间或阶跃响应的晚期渐近值可估计地电阻率分布.基于层状模型的大地冲激响应和阶跃响应正演计算结果表明,编码电磁测深法对大埋深薄层目标体有精细的分辨能力.可应用于能源、矿产资源、水资源、环境地质及工程地质勘察,有良好的应用前景.     

AN INTERPRETATION THEORY FOR INDUCED POLARIZATION VERTICAL SOUNDINGS (TIME-DOMAIN)*

In this paper it is shown how one may obtain a generalized Ohm's law which relates the induced polarization electric field to the steady-state current density through the introduction of a fictitious resistivity defined as the product of the chargeability and the resistivity of a given medium. The potential generated by the induced polarization is calculated at any point in a layered earth by the same procedure as used for calculating the potential due to a point source of direct current. On the basis of the definition of the apparent chargeability m_a, the expressions of m_a for different stratigraphie situations are obtained, provided the IP measurements are carried out on surface with an appropriate AMNB array. These expressions may be used to plot master curves for IP vertical soundings. Finally some field experiments over sedimentary formations and the quantitative interpretation procedure are reported.     

激发极化法的“饱和效应”问题

西格尔（H.O.Seigel）从理论计算得到导电颗粒含量为25％的极化球体有最大异常,因此认为:激发极化法只宜于寻找中等浸染程度的矿体,对于含量很高的矿体或良导致密型矿体,激发极化法是不行的。据称这种“饱和效应”的理论得到野外结果的证实。 本文根据我国激发极化法在良导致密型矿体上大都能获得明显视极化率异常的事实,说明西格尔的“饱和效应”理论与事实不符。本文认为良导致密型极化体不存在“饱和效应”问题,激发极化法对它们的勘探也是有效的。可能由于西格尔把低含量理论公式,不计条件地推导到高含量、甚至良导矿体,因此他的结论是错误的。本文最后讨论了电阻率因素对浸染型球体视极化率强度的影响。     

瞬变电磁法的探测深度问题

用解析分析、时域有限差分、时-频分析的方法,以地面中心回线装置和阶跃脉冲激励源为例,分析讨论了瞬变电磁测深法的勘探深度问题,以便为野外勘探工作设计提供依据,达到预期的探测目的.解析计算证实了瞬变场在地下以有限速度传播,数值模拟表示出了准静态条件下瞬变场的反射.研究结果表明,由于时间域电磁场遵循因果律,瞬变电磁法的探测深度主要由观测时间决定. 瞬变电磁场的初始传播速度与大地电阻率无关,继后在大地色散作用下,阶跃脉冲前沿逐渐变得平缓,各频率分量的传播速度与电阻率有关,在低阻地层中探测同样的深度需要较长的观测时间. 最大探测深度是在给定时间内电磁波往返地下某一深度的单程距离,最小探测深度受仪器性能的限制,但是埋藏较浅的异常体也有可能在晚时段被观测到.从时-频密度谱中可得到瞬变电磁场信号时间与频率的关系.     

Three-dimensional inversion of induced polarization data from simulated waste

The Idaho National Laboratory (INEL) Cold Test Pit (CTP) has been carefully constructed to simulate buried hazardous waste sites. An induced polarization (IP) survey of the CTP shows a very strong polarization and a modest resistivity response associated with the simulated waste. A three-dimensional (3-D) inversion algorithm based on the simultaneous iterative reconstruction technique (SIRT) and finite difference forward modelling has been applied to generate a subsurface model of complex resistivity. The lateral extents of the waste zone are well resolved. Limited depth extent is recognized, but the bottom of the waste appears too deep. With a modelling experiment, the intrinsic polarizability of the waste material is determined. Since IP is a technique for detection of diffuse occurrences of metallic material, this method holds promise as a method to distinguish buried waste from conductive soil material.     

Combined straightforward inversion of resistivity and induced polarization (time-domain) sounding data

It is proposed that the Straightforward Inversion Scheme (SIS) developed by the authors for 1D inversion of resistivity sounding and magneto-telluric sounding data can also be used in similar fashion for time-domain induced polarization sounding data. The necessary formulations based on dynamic dipole theory are presented. It is shown that by using induced polarization potential, measured at the instant when steady state current is switched off, an equation can be developed for apparent ‘chargeability–resistivity’ which is similar to the one for apparent resistivity. The two data sets of apparent resistivity and apparent chargeability–resistivity can be inverted in a combined manner, using SIS for a common uniform thickness layer earth model to estimate the respective subsurface distributions of resistivity and chargeability–resistivity. The quotient of the two profiles will give the sought after chargeability profile. A brief outline of SIS is provided for completeness. Three theoretical models are included to confirm the efficacy of SIS software by inverting only the synthetic resistivity sounding data. Then one synthetic data set based on a geological model and three field data sets (combination of resistivity and IP soundings) from diverse geological and geographical regions are included as validation of the proposal. It is hoped that the proposed scheme would complement the resistivity interpretation with special reference to shaly sand formations.     

1D and 2D Cole-Cole-inversion of time-domain induced-polarization data

A new method for the 2D inversion of induced polarization (IP) data in the time domain has been developed. The entire IP transients were observed and inverted into 2D Cole-Cole earth models, including resistivity, chargeability, relaxation time and the frequency constant. Firstly, a modified 1D time-domain electromagnetic algorithm was used to calculate the response of a layered polarizable ground. The transient signals were then inverted using the Marquardt method to derive the Cole-Cole parameters of each layer. However, model calculations showed that the EM effects could be neglected for the time range (>1 ms) and for the transmitter–receiver distances (<50 m) used in this study. Therefore, the induction effects were not considered for the solution of the 2D inverse problem and a DC solution was applied. An approximative forward algorithm was introduced in order to calculate the IP transients directly in the time domain and in order to speed up the inverse procedure. The approximation is highly accurate, and this is demonstrated by comparing the approximations with their exact solutions up to 3D. The inverse algorithm presented consists of two steps. The transient voltages of an array data set were inverted separately into a two-dimensional resistivity model for each time channel. The time-dependent resistivity of each cell was then interpreted as the response of a homogeneous half-space. In the 2D inversion algorithm, a 3D DC algorithm was used as a forward operator. The method only requires a standard 2D DC inversion and a homogenous half-space Cole-Cole inversion. The developed algorithm has been successfully applied to synthetic data sets and to a field data set obtained from a waste site situated close to Düren in Germany.     

从瞬变电磁响应中提取IP信息的研究(英文)

在对浙江省某铜矿进行瞬变电磁勘查时,由于感应激发极化效应的影响,造成瞬变电磁晚期测道的数据发生倒转。采用常规的瞬变电磁数据处理方法难以处理,造成晚期测道数据不可用,影响TEM的探测深度及精度。针对这个情况,采用Cole.Cole模型分析了均匀半空问模型中直流电阻率、充电率、时间常数以及频率相关系数对磁性源瞬变电磁响应的影响规律。利用奇异值分解法(The Singular Value Decomposition,简称SVD)对实测TEM数据进行分析及反演,并且从瞬变电磁响应数据中分离出勘探区瞬变电磁测量数据Cole—Cole模型的各参数,将分离出来的参数用于探测结果的辅助解释。经过分析得出充电率和频率相关系数对瞬变电磁响应的影响较大,而直流电阻率和时间常数对其影响较小的结论。通过常规方法和奇异值分解法对实测数据分别进行处理,发现后者视电阻率断面图的异常更加突出,从而较准确的圈定了铜矿体地范围,此结论得到_『钻孔的验证,与实际地质情况吻合较好。     

水下激发极化测深异常特征研究

利用有限单元法对水下激发极化测深进行数值模拟以研究其异常特征和适用条件.水下的电阻率测深和激电测深均可对水底地质体进行近距离探测,提高探测分辨率.而水底地形对视电阻率测深的影响很大.水底极化体有限元正演模拟的结果表明,水底地形对激发极化测深结果没有影响.采用水下激发极化法进行地质体探测是可行、有效的.通过计算不同水体和水底岩石电阻率对视极化率幅值的影响发现,水体电阻率是影响水下极化率测深的主要因素.由于海水的电阻率过低,观测视极化率异常微弱而不适合开展水下激电工作.     

时间域航空电磁扩散特征和成像深度研究

为了研究时间域航空电磁系统扩散特征,基于连续性边界条件将电磁场向下延拓得到地下介质中各点处的频率域响应,通过傅里叶变换将其变换到时间域.计算阶跃波形发射时垂直磁偶极子和水平磁偶极子在均匀半空间介质和两层介质中产生的电磁场,通过电流密度矢量图和随时间变化的电流密度等值线展示电磁场在地下介质中的扩散过程.电磁场在地下的扩散受电阻率结构的影响.在良导地区电磁场扩散慢、衰减快,而在高阻地区电磁场扩散快、衰减慢.垂直磁偶极子在地下产生的感应电流形成一个电流环随时间向下、向外扩散.水平磁偶极子在地下产生的感应电流形成两个互相叠加的电流环.时间域电磁场扩散代表着真正物理意义上的电磁扩散,因为它展示了地下电磁场空间分布随时间的变化和场强的衰减.基于对感应电流环的研究,我们探究了扩散深度和时间域航空电磁数据成像深度的关系.时间域电磁场扩散规律的研究不仅加深了对电磁扩散特征的理解,更有助于提高航空电磁数据解释水平.     

Potential Field,Geoelectrical and Reflection Seismic Investigations for Massive Sulphide Exploration in the Skellefte Mining District,Northern Sweden

Multi-scale geophysical studies were conducted in the central Skellefte district (CSD) in order to delineate the geometry of the upper crust (down to maximum ～ 4.5 km depth) for prospecting volcanic massive sulphide (VMS) mineralization. These geophysical investigations include potential field, resistivity/induced polarization (IP), reflection seismic and magnetotelluric (MT) data which were collected between 2009 and 2010. The interpretations were divided in two scales: (i) shallow (～ 1.5 km) and (ii) deep (～4.5 km). Physical properties of the rocks, including density, magnetic susceptibility, resistivity and chargeability, were also used to improve interpretations. The study result delineates the geometry of the upper crust in the CSD and new models were suggested based on new and joint geophysical interpretation which can benefit VMS prospecting in the area. The result also indicates that a strongly conductive zone detected by resistivity/IP data may have been missed using other geophysical data.     

Sensitivity and inversion of marine electromagnetic data in a vertically anisotropic stratified earth

The controlled‐source electromagnetic (CSEM) and magnetotelluric method (MT) are two techniques that can be jointly used to explore the resistivity structure of the earth. Such methods have, in recent years, been applied in marine environments to the exploration and appraisal of hydrocarbons. In many situations the electric properties of the earth are anisotropic, with differences between resistivity in the vertical direction typically much higher than those in the horizontal direction. In cases such as this, the two modes of the time‐harmonic electromagnetic field are altered in different ways, implying that the sensitivity to the earth resistivity may vary significantly from one particular resistivity component (scalar, horizontal or vertical) to another, depending on the measurement configuration (range, azimuth, frequency or water depth). In this paper, we examine the sensitivity of the electromagnetic field to a vertically anisotropic earth for a typical set of configurations, compare inversion results of synthetic data characterizing a vertically anisotropic earth obtained using the isotropic and anisotropic assumptions and show that correctly accounting for anisotropy can prevent artefacts in inversion results.     

Interpretation of VLF Resistivity Data for Ground Water Contamination Surveys

Very low frequency (VLF) military communications systems provide a primary field that can be used for shallow geophysical surveys to locate ground water contamination and vertical geologic contacts. Useful properties that can be easily obtained from the interaction of the earth and the primary field are the magnitude of the vertical secondary magnetic field, the surface impedence, and the phase angle between the electrical and magnetic horizontal components. The variations in the secondary magnetic field can be related to vertical geologic contacts, such as the edges of landfill trenches. The surface impedence yields an apparent terrain conductivity, which can be used to locate low-resistivity anomalies often associated with contaminated ground water. The phase angle gives information on vertical variations in resistivity, phase angles less than 45° indicating increasing resistivity with depth. The depth of penetration of the VLF field is about one skin depth. For a frequency of 20 kHz, the skin depth in meters is approximately equal to 3.67 where p is terrain resistivity in ohmmeters.     

Depth of detection of highly conducting and volume polarizable targets using induced polarization

We define the apparent frequency effect in induced polarization (IP) as the relative difference between apparent resistivities measured using DC excitation on the one hand and high‐frequency excitation (when the IP effect vanishes) on the other. Assuming a given threshold for the minimum detectable anomaly in the apparent frequency effect, the depth of detection of a target by IP can be defined as that depth below which the target response is lower than the threshold for a given electrode array. Physical modelling shows that for the various arrays, the depth of detection of a highly conducting and volume polarizable target agrees closely with the depth of detection of an infinitely conducting and non‐polarized body of the same shape and size. The greatest depth of detection is obtained with a two‐electrode array, followed by a three‐electrode array, while the smallest depth of detection is obtained with a Wenner array when the array spread is in‐line (i.e. perpendicular to the strike direction). The depth of detection with a Wenner array improves considerably and is almost equal to that of a two‐electrode array when the array spread is broadside (i.e. along the strike direction).     

IP RESPONSE FOR A 2D HORIZONTAL CYLINDER1

The induced polarization response for a 2D horizontal cylinder embedded in a half-space is calculated for a uniform electric source. Response curves, in the form of apparent charge-ability taking into account the effect of the air-earth interface, exhibit a sharp decrease in amplitude with an increase in depth of burial of the target. The resistivity contrast between the cylinder and the host plays a dominant role in determining the IP response, i.e. the amplitude decreases considerably with the increase in resistivity contrast. The decrease is due to the defocusing effect caused by the resistive cylinder. The current lines tend to deviate away from the cylindrical target. In the case of a highly conducting cylinder, apparent defocusing takes place as current lines are confined to the surface of the conducting cylinder. An increase in chargeability contrast is reflected as a steady rise in the response. The peak response at the centre is reduced by about half the magnitude when the air–earth interface is not considered. The variation of response along the profile, though noticeable, is not as high as that obtained at the centre.