Effect of clayey media parameters on the negative response of a coincident loop1

The aim of this work is to consider the phenomenon of negative response and its detectability in a polarizable half-space formed of Clay-water mixture which commonly exists in water aquifers and which can be used as a simple indicator of the presence of Underground water. We make use of an analytical treatment for the transient voltage induced in a coincident loop lying on a half-space as a basis for our computations. A Cole-Cole model is used to represent the Clay with appropriate parameters. As saturated Clay is characterized by high conductivity and moderate chargeability, it is hardly detectable relative to the practical noise level. However, there are optimum Clay parameters that offer a maximum negative response. It is also observed that for a very short time constant τ in the Cole-Cole model, the amplitude of the negative response decreases with the decrease of the time constant due to the fast decay of the polarization current. Finally, it is found that there is an optimum value for the loop radius which gives the largest value of the negative response and this loop radius depends on the model parameters.     

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.     

用柱状成层各向异性介质的并矢Green函数模拟偏心对多分量感应测井响应的影响

采用偏心状态下柱状成层各向异性(横向各向同性)介质中并矢Green函数的解析表达式高效模拟多分量感应测井仪器在井眼中偏心时的响应.为提高精度,在模拟时考虑到了金属心轴、绝缘保护层的存在以及各分量线圈系的具体形状.数值模拟结果表明,当井眼内钻井液电导率相对较高、地层电导率相对较低时,偏心对仪器响应的影响较大,尤其是对短线圈距线圈系的影响更为明显,必须进行偏心效应校正.当钻井液电导率相对较低时,偏心对线圈系响应的影响可忽略不计.对位于相对低电导率井眼中的线圈系而言仪器方位角的影响可忽略不计,而当线圈系位于相对高电导率井眼中时仪器方位角的影响极为明显.当仪器偏心率较小时线圈系的响应随仪器方位角的变化较小,仪器偏心率越大线圈系的响应随仪器方位角的变化越明显.     

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

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

Time domain electromagnetic‐induced polarisation: extracting more induced polarisation information from grounded source time domain electromagnetic data

Electrical induced polarisation surveys are used to detect chargeable materials in the earth. For interpretation of time domain electrical‐induced polarisation data a common procedure is to first invert the direct current data (electric current on time) to recover conductivity and then invert the induced polarisation data (current off‐time) to recover chargeability. This direct current‐induced polarisation inversion procedure assumes that the off time data are free of secondary electromagnetic induction effects. To comply with this, early time data are often discarded or not recorded. For mid‐time data, an electromagnetic decoupling technique, which removes electromagnetic induction in the observations, needs to be implemented. Usually, responses from a half‐space or a layered earth are subtracted. Recent capability in three‐dimensional time domain electromagnetic forward modelling and inversion allows to revisit these procedures. In a Time domain electromagnetic‐induced polarisation survey, a high sampling rate allows early time channels of the electromagnetic data to be recorded. The recovery of chargeability then follows a three‐step workflow: (i) invert early time channel time domain electromagnetic data to recover the three‐dimensional conductivity; (ii) use that conductivity to compute the time domain electromagnetic response at later time channels and subtract this fundamental response from the observations to extract the induced polarisation responses, and (iii) invert the induced polarisation responses to recover a three‐dimensional chargeability. This workflow effectively removes electromagnetic induction effects in the observations and produces better chargeability and conductivity models compared with conventional approaches. In a synthetic example involving a gradient array, we show that the conductivity structure obtained from the early time channel data, which are usually discarded, is superior to that obtained from the steady state direct current voltages. This adds a further reason to collect these electromagnetic data.     

基于加窗Hilbert变换的复偏振分析方法及其应用

考虑到传统Hilbert变换所隐含的无限长序列假设的实际局限性, 本文将多窗分析法引入到短时序列Hilbert变换中, 通过构建复协方差矩阵和对该矩阵的特征值与特征向量的求解, 获得三分量地震记录的时变偏振参数;基于不同波型的偏振特性与实测偏振参数, 采用线性_余弦权重函数或高斯权重函数自适应空间滤波, 识别与分离具有特定偏振特性的不同地震响应. 针对天然地震三分量记录的处理结果表明, 该方法在识别、分离具有特定偏振特性的不同地震响应方面具有一定的潜力.     

Evaluation of transport and storage properties in the soil and groundwater zone from induced polarization measurements1

Spectral induced polarization as well as complex electrical measurements are used to estimate, on a non-invasive basis, hydraulic permeability in aquifers. Basic laboratory measurements on a variety of shaly sands, silts and clays showed that the main feature of their conductivity spectra in the frequency range from 10^-3 to 10³ Hertz is a nearly constant phase angle. Thus, a constant-phase-angle model of electrical conductivity is applied to interpret quantitatively surface and borehole spectral induced polarization measurements. The model allows for the calculation of two independent electrical parameters from only one frequency scan and a simple separation of electrical volume and interface effects. The proposed interpretation algorithm yields the true formation factor, the cation exchange capacity and the surface-area-to-porosity ratio, which corresponds to the inverse hydraulic radius. Using a Kozeny–Carman-like equation, the estimation of hydraulic permeability is possible.     

Current and voltage source induced polarization transients: a comparative consideration

The article discusses the excitation of transient induced polarization responses using current and voltage sources. The first method has found a wide application in induced polarization surveys and—directly or indirectly—in the theory of the induced polarization method. Typically, rectangular current pulses are injected into the earth via grounding electrodes, and decaying induced polarization voltage is measured during the pauses between pulses. In this case, only the secondary field is recorded in the absence of the primary field, which is an important advantage of this method. On the other hand, since the current injected into the ground is fully controlled by the source, this method does not allow studying induced polarization by measuring the current in the transmitter line or associated magnetic field. When energising the earth with voltage pulses, the measured quantity is the transient induced polarization current. In principle, this method allows induced polarization studies to be done by recording the transmitter line current, the associated magnetic field, or its rate of change. The decay of current in a grounded transmitter line depends not only on the induced polarization of the earth but also on the polarization of the grounding electrodes. This problem does not occur when induced polarization transients in the earth are excited inductively. A grounded transmitter line is a mixed‐type source; hence, for a purely inductive excitation of induced polarization transients, one should use an ungrounded loop, which is coupled to the earth solely by electromagnetic induction.     

VERTICAL FREQUENCY EFFECT SOUNDING IN INDUCED POLARIZATION AND GALVANIC RESISTIVITY METHODS *

It has been shown that for a polarizable layer as a transmission medium there is an indirect proportionality between the frequency effect (fe) in induced polarization (IP) and the wave number of the electrical induced field. Making use of this relationship for a two layered earth a polarization transform function has been obtained. Since the mathematical expression for the polarization transform function is the same as that of the resistivity transform function, it is possible to make direct interpretation for IP frequency effect field curves. Thus, AA or QQ type resistivity sequences can be interpreted from induced polarization response of a horizontally stratified earth without resistivity extrema. A depth factor has been defined in order to obtain the true depth using the apparent depth. In this way, some electromagnetic effects between horizontal layers with different polarizabilities can partly be eliminated.     

Time-Domain Spectral Induced Polarization Based on Pseudo-random Sequence

To reduce noise during electrical prospecting, we hereby propose a new method using correlation identification technology and conventional electrical exploration devices. A correlation operation can be carried out with the transmitted pseudo-random sequence and received time signal to suppress the random noise, and the time-domain impulse response and frequency response of the frequency domain of the underground media can be obtained. At the same time, using a dual Cole–Cole model to fit a complex resistivity spectrum, which is close to the frequency response, we can get a variety of induced polarization parameters and electromagnetic parameters of subsurface, which can provide more useful information for the exploration of mineral resources. This time domain prospecting method can effectively improve the efficiency of the spectral induced polarization method. In this article, we have carried out theoretical calculations and a simulation to prove the feasibility of such a method.     

层状导电介质中地面核磁共振响应特征理论研究

水平层状导电大地中层状含水层地面核磁共振响应的数值计算技术涉及到导电介质中回线源 磁场的计算以及地下含水层中质子磁矩在线圈中产生感应信号的体积积分. 文中采用直接数 值积分技术,对具有振荡核函数的Hankel变换进行积分,以求取回线源磁场的径向与垂向分 量,计算并研究了回线产生的径向与垂向磁场分量随空间位置的变化规律. 基于磁场空间分 布特点,利用不等间距空间剖分技术计算地面核磁共振的体积积分,模拟了不同模型的地面 核磁共振响应并讨论了其影响因素. 结果表明,结合能对任意层状导电介质中磁场进行稳定 快速计算的直接数值积分技术与不等距空间模型剖分技术,可正确模拟地面核磁共振响应. 导电性是产生地面核磁共振信号相位的先决条件,但影响响应振幅强度与相位的因素还有含 水层的埋深、厚度以及装置大小等.     

EM AND IP RESPONSE OF A STEEL WELL CASING FOR A FOUR-ELECTRODE SURFACE ARRAY. PART I: THEORY*

The basic formulation for the mutual coupling between grounded circuits in the presence of a vertical conductor is used to examine the full electromagnetic and induced polarization response for a homogeneous half-space model. Both interfacial polarization at the steel/electrolyte boundary and induced eddy currents in the casing are accounted for in the analysis. It is shown that for mild steel casings the eddy currents, hitherto neglected, play an important role because of their strong frequency dependence. Calculated curves of the appar- ent complex resistivity as a function of frequency show that the overall response can be influenced by the state of corrosion at the surface of the casing.     

The influence of induced polarization on long-offset transient electromagnetic data

The diffusion of electromagnetic fields is dependent not only on conductivity, but also on magnetic permeability, dielectric permittivity and polarizability, i.e. dispersive conductivity. The long‐offset transient electromagnetic (LOTEM) method is mainly used to determine the spatial distribution of conductivity in the subsurface. However, earlier work on loop‐loop TEM suggests that transient EM methods can also be affected by induced polarization (IP). Numerous 1D forward calculations were carried out to study the IP effect on LOTEM data, using the Cole‐Cole relaxation model to simulate the polarizability of the ground. Besides the polarizability of each layer, the IP effect depends on the LOTEM field set‐up and the spatial distribution of conductivity in the ground. In particular, near‐surface layers with high chargeabilities can significantly distort the late time transients of the electric field components in the vicinity of the transmitter. The influence of polarizable layers on the magnetic field components can be neglected under normal circumstances. In 1997 and 1999, LOTEM measurements were carried out at Mt. Vesuvius in Italy to explore the geological structure of the volcano. Sensitivity studies on the effect of polarizable layers suggest that high chargeabilities in connection with conductive layers at greater depths would result in a detectable distortion of the electric field transients. Although the simultaneous IP measurements revealed high chargeabilities in a near‐surface layer, no evidence of IP effects could be found in the measured LOTEM data. We conclude that the observed chargeabilities are local and that 3D effects are probably present in the data. Another aspect is the measurement of the system response, which is usually measured by placing a receiver very close to the transmitter. Therefore, large distortions can be expected if near‐surface polarizable layers exist. This was verified in practice by field measurements in an area with high chargeabilities in Longerich, Cologne.     

Geomagnetic sounding of conductivity anomalies in the lower crust and uppermost mantle

Knowledge of the structure of the lower crust and uppermost mantle is of special importance in understanding plate tectonics. Electrical conductivity of this region has been measured recently in various parts of the world. Transfer functions are still the most widely used quantity in data analysis and model fitting. Anomalies in the horizontal magnetic components in combination with anomalies in the vertical component have been found very useful in locating conductivity contrasts. With these, when the cause of the anomaly is a concentrated line current, both the position of the line current and its depth can be directly located. The method of hypothetical event analysis is another new technique and this is highly suited for areas having complex subsurface geology or areas under non-uniform source fields or both. The anomaly is more suitable for modelling geophysical structures when it is separated into regional and local components. Model calculations still are not very satisfactory and the importance of one dimensional calculations must be emphasized for they give direct information on the variation of conductivity with depth, which is the purpose of GDS. We need more results, especially from tectonically active areas, before the underlying physical processes can be completely understood.     

Electric field variations measured continuously in free air over a conductive thin zone in the tilted Lias-epsilon black shales near Osnabrück,Northwest Germany

Common studies on the static electric field distribution over a conductivity anomaly use the self-potential method. However, this method is time consuming and requires nonpolarizable electrodes to be placed in the ground. Moreover, the information gained by this method is restricted to the horizontal variations of the electric field. To overcome the limitation in the self-potential technique, we conducted a field experiment using a non conventional technique to assess the static electric field over a conductivity anomaly. We use two metallic potential probes arranged on an insulated boom with a separation of 126 cm. When placed into the electric field of the free air, a surface charge will be induced on each probe trying to equalize with the potential of the surrounding atmosphere. The use of a plasma source at both probes facilitated continuous and quicker measurement of the electric field in the air. The present study shows first experimental measurements with a modified potential probe technique (MPP) along a 600-meter-long transect to demonstrate the general feasibility of this method for studying the static electric field distribution over shallow conductivity anomalies.Field measurements were carried out on a test site on top of the Bramsche Massif near Osnabrück (Northwest Germany) to benefit from a variety of available near surface data over an almost vertical conductivity anomaly. High resolution self-potential data served in a numerical analysis to estimate the expected individual components of the electric field vector.During the experiment we found more anomalies in the vertical and horizontal components of the electric field than self-potential anomalies. These contrasting findings are successfully cross-validated with conventional near surface geophysical methods. Among these methods, we used self-potential, radiomagnetotelluric, electric resistivity tomography and induced polarization data to derive 2D conductivity models of the subsurface in order to infer the geometrical properties and the origin of the conductivity anomaly in the survey area. The presented study demonstrates the feasibility of electric field measurements in free air to detect and study near surface conductivity anomalies. Variations in Ez correlate well with the conductivity distribution obtained from resistivity methods. Compared to the self-potential technique, continuously free air measurements of the electric field are more rapid and of better lateral resolution combined with the unique ability to analyze vertical components of the electric field which are of particular importance to detect lateral conductivity contrasts. Mapping Ez in free air is a good tool to precisely map lateral changes of the electric field distribution in areas where SP generation fails. MPP offers interesting application in other geophysical techniques e.g. in time domain electromagnetics, DC and IP.With this method we were able to reveal a ca. 150 m broad zone of enhanced electric field strength.     

The polarization of Pc5 ULF waves around dawn: A possible ionospheric conductivity gradient effect

In this study over 100 Pc5 events observed on the SAMNET magnetometer array have been catalogued. The ground horizontal polarization of these waves has been investigated and it has been shown that there are large differences between the polarization azimuth distributions on either side of noon. In the morning the azimuth distributions are spread evenly across all possible orientations whereas in the afternoon they are predominantly in a North-South orientation. The difference appears to be the result of a specific polarization pattern that occurred across the SAMNET array for 31 of the Pc5 events. These 31 Pc5 events all occurred in the morning sector, close to the dawn terminator. It is suggested that this special polarization pattern is a consequence of ionospheric conductivity gradients that occur at around dawn. This theory is reinforced by the comparison of the local times of these polarization variations with the local time of sunrise. By employing ionospheric conductivity data from the EISCAT radar for two of these Pc5 events, this connection has been studied in more detail.     

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.     

EM AND IP RESPONSE OF A STEEL WELL CASING FOR A FOUR-ELECTRODE SURFACE ARRAY. PART II: NUMERICAL RESULTS*

Some numerical results for the apparent resistivity of a homogeneous half space containing a vertical steel well casing as measured by a dipole-dipole array are presented. The basic model includes the electromagnetic and induced polarization effects by allowing for frequency dependent eddy currents and an electrochemical boundary of the well casing. It is shown that the well casing produces an apparent complex resistivity response that could be mistaken for an induced polarization anomaly existing above a hydrocarbon deposit. This response of the well casing is strongly dependent on frequency.     

面向目标自适应有限元法的带地形三维大地电磁各向异性正演模拟

传统三维大地电磁各向异性模拟均是基于规则六面体网格,计算精度有限且较难拟合复杂地质条件.本文采用面向目标自适应非结构矢量有限元法,对三维大地电磁各向异性介质进行模拟.首先从电场双旋度方程出发,利用伽辽金方法建立变分方程;然后利用电流密度连续性条件构建适合大地电磁各向异性问题的加权后验误差估计方法,实现面向目标的网格自适应正演;最后通过典型算例分析各向异性对网格自适应和大地电磁响应的影响特征以及各向异性的识别方法.本文算法能够高精度地拟合起伏地表和任意各向异性介质,适用于分析复杂地电条件大地电磁响应特征,为提高大地电磁资料解释水平提供了理论基础.