Ground penetrating radar inversion in 1-D: an approach for the estimation of electrical conductivity, dielectric permittivity and magnetic permeability

This paper presents a method for inverting ground penetrating radargrams in terms of one-dimensional profiles. We resort to a special type of linearization of the damped E-field wave equation to solve the inverse problem. The numerical algorithm for the inversion is iterative and requires the solution of several forward problems, which we evaluate using the matrix propagation approach. Analytical expressions for the derivatives with respect to physical properties are obtained using the self-adjoint Green's function method. We consider three physical properties of materials; namely dielectrical permittivity, magnetic permeability and electrical conductivity. The inverse problem is solved minimizing the quadratic norm of the residuals using quadratic programming optimization. In the iterative process to speed up convergence we use the Levenberg–Mardquardt method. The special type of linearization is based on an integral equation that involves derivatives of the electric field with respect to magnetic permeability, electrical conductivity and dielectric permittivity; this equation is the result of analyzing the implication of the scaling properties of the electromagnetic field. The ground is modeled using thin horizontal layers to approximate general variations of the physical properties. We show that standard synthetic radargrams due to dielectric permittivity contrasts can be matched using electrical conductivity or magnetic permeability variations. The results indicate that it is impossible to differentiate one property from the other using GPR data.     

POTENTIAL FIELD OF A STATIONARY ELECTRIC CURRENT USING FREDHOLM'S INTEGRAL EQUATIONS OF THE SECOND KIND*

An integral equation method is described for solving the potential problem of a stationary electric current in a medium that is linear, isotropic and piecewise homogeneous in terms of electrical conductivity. The integral equations are Fredholm's equations of the ‘second kind’ developed for the potential of the electric field. In this method the discontinuity-surfaces of electrical conductivity are divided into ‘sub-areas’ that are so small that the value of their potential can be regarded as constant. The equations are applied to 3-D galvanic modeling. In the numerical examples the convergence is examined. The results are also compared with solutions derived with other integral equations. Examples are given of anomalies of apparent resistivity and mise-a-la-masse methods, assuming finite conductivity contrast. We show that the numerical solutions converge more rapidly than compared to solutions published earlier for the electric field. This results from the fact that the potential (as a function of the location coordinate) behaves more regularly than the electric field. The equations are applicable to all cases where conductivity contrast is finite.     

Fast and accurate three‐dimensional controlled source electromagnetic modelling†

We present a fast approximate method for three‐dimensional low frequency controlled source electro‐magnetic modeling. We apply the method to a synthetic model in a typical marine controlled source electromagnetic scenario, where conductivity and permittivity are different from the known background medium. For 3D configurations, fast computational methods are relevant for both forward and inverse modelling studies. Since this problem involves a large number of unknowns, it has to be solved efficiently to obtain results in a timely manner, without compromising accuracy. For this reason, the Born approximation, extended Born approximation and iterative extended Born approximation are implemented and compared with the full solution of the conjugate gradient fast Fourier transformation method. These methods are based on an electric field domain integral equation formulation. It is shown here how well the iterative extended Born approximation method performs in terms of both accuracy and speed with different configurations and different source positions. The improved accuracy comes at virtually no additional computational cost. With the help of this method, it is now possible to perform sensitivity analysis using 3D modelling in a timely manner, which is vital for controlled source electromagnetic applications. For forward modeling the solution at the sea‐bottom is of interest, because that is where the receivers are usually located. For inverse modeling, the accuracy of the solution in the target zone is important to obtain reasonably accurate conductivity values from the inversion using this approximate solution method. Our modelling studies show that the iterative extended Born approximation method is fast and accurate for both forward and inverse modelling. Sensitivity analysis as a function of the source position and different reservoir sizes validate the accuracy of the iterative extended Born approximation.     

二维介质柱的电磁成像研究

从体等效原理出发得到用介质柱介电参数表示的体等效电流,并得到一个积分方程组求解散射问题．利用变分原理,对此方程组求变分,并化为矩阵方程组,消去含介质柱内总场变化的项,得到散射场变化与介质柱介电参数变化之间关系的矩阵方程．即反演方程．用算例检验此法的反演精度、收敛性能、多目标反演能力、抗噪声性能及利用不完备测量散射信息反演目标的能力．结果表明,此法是一种实用性很强的二维介质柱电磁成像方法．     

井间电磁测量的2．5维层析成像方法

利用正则化最小二乘反演方法实现了井间电磁测量数据的层析成像,对井间地层电阻率进行了重建。在成像算法中,我们假设了井间电磁的激发与接收采用电磁偶极子源,井间介质仅在二维（xoz）平面内变化。在数值模拟中,通过对构造走向（y方向）的Fourier变换,将三维电磁场问题转化为一系列二维问题,用等参有限元方法在波数域求解,使实际地层模型的处理得以实现。对于波数域中每个波数对应的电磁场方程采用等参有限元求解,并用高斯积分将波数域解变换为空间域电磁场。利用源与接收器电磁场的互易原理,实现了电磁场响应对电导率分布灵敏度的快速计算。针对正演模拟中源点的奇异性,我们采用具有一定面积的伪艿函数表达源电流分布,使数值解精度得到提高。用层状介质的解析解与数值计算结果的对比,验证了模拟算法的精度。用介质扰动产生的电磁场变化检验互易性定理计算灵敏度的有效性。对简单块状模型、斜向裂缝带模型及“大”字模型的模拟数据成像结果表明,本文介绍的层析成像方法是正确有效的。     

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

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

位场向下延拓的波数域广义逆算法

位场向下延拓是位场数据处理和反演中的重要运算,但是它的不稳定性影响了它在许多处理和反演方法技术中的应用.本文通过把位场向下延拓视为向上延拓的反问题,得到向下延拓的褶积型线性积分方程,再利用Fourier变换矩阵的正交对称特性,并结合矩阵的奇异值分解和广义逆原理,提出了一种稳定的不需要进行求逆运算的位场向下延拓广义逆方法——波数域广义逆算法,解决了位场大深度向下延拓的不稳定性问题.把这种方法用于三维理论模型数据和实际磁场数据的向下延拓获得了理想的结果.     

频率域2维电磁测深资料的RRI反演

本文将大地电磁反演中的RRI方法应用于线源频率测深模型资料的反演中.当背景电导率变化很小时,可用变化前的电场来近似代替变化后的电场,从而在反演方程的推导过程中将有源电磁场中的源项消掉,得到和大地电磁场相同的反演方程,使有源电磁波的反演也可以应用RRI方法.反演过程中所需要的模型资料通过有限元方法得到,该资料不需做近场校正,直接用适合于有源电磁场的RRI方法反演,避免了近场校正带来的误差.数值模型结果证明该方法是可行的.最后,用RRI方法讨论了当源和目标区间存在低阻异常体时只对目标区反演的可行性,对实际工作的解释有一定的指导意义.     

基于二次场电导率分块连续变化的三维可控源电磁有限元数值模拟

从电偶源三维地电断面可控源电磁法的二次电场边值问题及其变分问题出发,采用任意六面体单元对研究区域进行剖分,并且在单元分析中同时对电导率及二次电场进行三线性插值,实现电导率分块连续变化情况下,基于二次场的可控源电磁三维有限元数值模拟.这个新的可控源电磁三维正演方法可以模拟实际勘探中地下任意形状及电性参数连续变化的复杂模型.理论模型的计算结果表明,均匀大地计算的视电阻率误差和相位误差分别为0.002%和0.0005°.分层连续变化模型的有限元计算结果表明,其与对应的分层均匀模型解析结果有明显差异.三维异常体组合模型以及倾斜异常体等复杂模型的有限元计算结果也有效地反映了异常形态.     

一种新的实时电磁逆散射方法

为解决介质圆柱体逆散射问题,提出一种新的在线逆散射方法,通过支持向量机将逆散射问题转化成一个回归估计问题. 该方法可应用于各种逆散射方面, 尤其是目标的几何与电磁参数重构和埋地目标探测. 文中首次将支持向量机方法应用到该领域,设置多个散射场的观测点,通过提取散射场的不同信息作为样本信息训练支持向量机, 建立了介质圆柱体的逆散射模型, 利用该模型重构了介质圆柱体的电磁参数,同时探测了埋地位置. 数值结果显示了该方法的有效性和准确性,为目标的实时逆散射研究提供了一种有效方法.     

3D modelling and sensitivity in DC resistivity using charge density

A three‐dimensional (3D) electrical resistivity modelling code is developed to interpret surface and subsurface data. Based on the integral equation, it calculates the charge density caused by conductivity gradients at each interface of the mesh, allowing the estimation of the potential everywhere without the need to interpolate between nodes. Modelling generates a huge matrix, made up of Green's functions, which is stored by using the method of pyramidal compression. The potential is compared with the analytical and the numerical solutions obtained by finite‐difference codes for two models: the two‐layer case and the vertical contact case. The integral method is more accurate around the source point and at the limits of the domain for the potential calculation using a pole‐pole array. A technique is proposed to calculate the sensitivity (Jacobian) and Hessian matrices in 3D. The sensitivity is based on the derivative with respect to the block conductivity of the potential computed using the integral equation; it is only necessary to compute the electrical field at the source location. A direct extension of this technique allows the determination of the second derivatives. The technique is compared with the analytical solutions and with the calculation of the sensitivity according to the method using the inner product of the current densities calculated at the source and receiver points. Results are very accurate when the Green's function that includes the source image is used. The calculation of the three components of the electric field on the interfaces of the mesh is carried out simultaneously and quickly, using matrix compression.     

Stochastic joint inversion of hydrogeophysical data for salt tracer test monitoring and hydraulic conductivity imaging

The assessment of hydraulic conductivity of heterogeneous aquifers is a difficult task using traditional hydrogeological methods (e.g., steady state or transient pumping tests) due to their low spatial resolution. Geophysical measurements performed at the ground surface and in boreholes provide additional information for increasing the resolution and accuracy of the inverted hydraulic conductivity field. We used a stochastic joint inversion of Direct Current (DC) resistivity and self-potential (SP) data plus in situ measurement of the salinity in a downstream well during a synthetic salt tracer experiment to reconstruct the hydraulic conductivity field between two wells. The pilot point parameterization was used to avoid over-parameterization of the inverse problem. Bounds on the model parameters were used to promote a consistent Markov chain Monte Carlo sampling of the model parameters. To evaluate the effectiveness of the joint inversion process, we compared eight cases in which the geophysical data are coupled or not to the in situ sampling of the salinity to map the hydraulic conductivity. We first tested the effectiveness of the inversion of each type of data alone (concentration sampling, self-potential, and DC resistivity), and then we combined the data two by two. We finally combined all the data together to show the value of each type of geophysical data in the joint inversion process because of their different sensitivity map. We also investigated a case in which the data were contaminated with noise and the variogram unknown and inverted stochastically. The results of the inversion revealed that incorporating the self-potential data improves the estimate of hydraulic conductivity field especially when the self-potential data were combined to the salt concentration measurement in the second well or to the time-lapse cross-well electrical resistivity data. Various tests were also performed to quantify the uncertainty in the inverted hydraulic conductivity field.     

电导率各向异性的海洋电磁三维有限单元法正演

本文提出了一种基于非结构化网格的海洋电磁有限单元正演算法.为了回避场源奇异性,文中选用二次场算法,将背景电阻率设置为水平层状且各向异性,场源在水平层状各向异性介质中所激发的一次场通过汉克尔积分得到.基于Coulomb规范得到二次矢量位和标量位所满足的Maxwell方程组,通过Galerkin加权余量法形成大型稀疏有限元方程,采用不完全LU分解(ILU)预条件因子的quasi-minimum residual(QMR)迭代解法对有限元方程进行求解得到二次矢量位和标量位;进而,利用滑动平均方法得到二次矢量位和标量位在空间的导数,由此得到二次电磁场;通过一维模型对算法的可靠性进行验证,与此同时,针对实际复杂海洋电磁模型,比较有限元模拟结果与积分方程模拟结果,进一步验证算法精度.若干计算结果均表明,文中算法具有良好的通用性,适用于井中电磁、航空电磁,环境地球物理等非均匀且各向异性介质中的电磁感应基础研究.     

一种新的三维大地电磁积分方程正演方法

采用规则六面体单元和并矢Green函数奇异积分等效积分技术,已有的大地电磁积分正演方法具有不能有效模拟地下复杂地质体和计算精度偏低的缺点.本文提出了一种新的三维大地电磁积分方程正演技术,即采用四面体单元、解析的并矢Green函数奇异积分表达式,达到既能模拟地下复杂异常体,又能有效提高已有积分方程法计算精度的目的.首先,采用四面体网格技术离散地下复杂异常体,获得四面体单元上的大地电磁积分方程.然后,利用针对四面体单元开发的新的奇异值积分的解析表达式,准确计算线性方程中的并矢Green函数的奇异积分,从而获得精确的线性方程.借助于PARDISO高性能并行直接求解器,实现了三维大地电磁问题的高精度求解.最后,基于国际标准3D-1模型和六棱柱模型,通过与其他方法结果的对比分析,验证了本文方法的正确性、处理高电导率对比度的能力(1000:1)和处理复杂模型的能力.     

地面可控源频率测深三维非线性共轭梯度反演

讨论了地面可控源电磁勘探三维非线性共轭梯度反演的可行性以及反演过程中考虑场源的必要性.反演采用非线性共轭梯度反演方法.反演过程中,模型响应利用交错网格有限差分技术计算.反演数据采用与发射源平行的电场x分量E_x.利用层状导电模型作为背景,设计了两个理论模型进行数值试验:第一个模型中包含两个电阻率异常,以检验反演的有效性;第二个模型中,在测区外设置了一个低阻异常,以考察源的信息在反演中的作用.两个模型的反演分别从层状背景模型开始,迭代120次后终止.数值试验结果表明,(1)非线性共轭梯度反演所获得的电阻率分布和理论模型吻合较好;(2)非线性共轭梯度算法收敛速度较慢,需要较多的迭代次数完成反演;(3)对于可控源频率电磁勘探,必须考虑源位置信息.因此,本文采用考虑场源信息的地面可控源非线性共轭梯度反演方法能完成真正意义上的可控源频率电磁测深数据的反演.     

Electric field data in inductive source electromagnetic surveys

Measurement of the electric field data due to an inductive loop source in a controlled source electromagnetic survey is not common, because electric field data, usually involving grounded electrodes, are expensive to acquire and difficult to interpret. With the recently developed capability of versatile three‐dimensional inversion, we revisit the idea of measuring electric field in a large ground loop survey for mineral exploration. The three‐dimensional modelling and inversion approach helps us quantitatively understand the detectability and recoverability of the proposed survey configuration. Our detectability study using forward modelling shows that the relative anomaly (percentage difference) in electric field does not decay with a lower induction number, but the conventional magnetic field data (dB/dt) does. Our recoverability study examines how much and what kind of information can be extracted from electric field data for the reconstruction of a three‐dimensional model. Synthetic inversions show the following observations. (i) Electric field data are good at locating lateral discontinuity, whereas dB/dt has better depth resolution. (ii) Electric field is less sensitive to the background conductivity and, thus, is prone to misinterpretation because of a bad initial model in inversion. We recommend warm‐starting the electric field inversion with an initial model from a separate dB/dt inversion. (iii) Electric field data may be severely contaminated by near‐surface heterogeneity, but an inversion can recover the deep target concealed by the geologic noise. (iv) Even one line of single‐component electric field data can greatly improve the horizontal resolution in a dB/dt inversion. Finally, we investigate a field dataset of both electric field and dB/dt measurements at a uranium deposit. The field example confirms that the electric field and magnetic field data contain independent information that is crucial in the accurate recovery of subsurface conductivity. Our synthetic and field examples demonstrate the benefit of acquiring electric field data along with magnetic field data in an inductive source survey.     

Appreciation of spherically symmetric models of electrical conductivity

A procedure is suggested of a more effective and faster computation of the impedance, the transfer function and amplitudes of the induced field in a spherically symmetric model of the electrical conductivity. The existing induction data have been supplemented by about 80 new values derived from the analysis of daily means. The fit of the existing 1-D models of the electrical conductivity of the mantle to the set of induction data is investigated. The characteristic equation for the free electromagnetic oscillations of a radially inhomogeneous Earth is derived and its possible importance in solving the inverse problem of electric conductivity is pointed out.     

基于精确震源函数的解调包络多尺度全波形反演

本文提出解调包络方法来重构地震记录中缺失的低频信号,同时该方法能够降低全波形反演的非线性程度;提出伴随状态震源函数反演方法来得到精确的震源函数,并推导了梯度计算公式;解调包络方法结合低通滤波技术,实现了从低频到高频的多尺度反演策略,有效缓解了全波形反演的周波跳跃问题.数值算例证明了解调包络、伴随状态震源函数反演方法和低通滤波多尺度反演策略的可行性及优越性.震源函数反演精度测试结果表明:即使观测记录在缺失低频信息的情况下,也能反演得到精确的震源函数.缺失低频测试和抗噪能力测试结果表明:即使地震数据中缺失9Hz以下的低频信号或者信噪比极低的情况下,利用反演得到的精确震源函数进行解调包络多尺度全波形反演,同样可以得到高精度的全波形反演结果.与Hilbert包络全波形反演对比结果表明:解调包络在重构低频和降低伴随震源主频方面具有一定优势.     

基于概率成像技术的低纬度磁异常化极方法

化极转换是磁异常解释的重要基础,为了克服在地磁纬度较低的地区尤其是磁赤道附近化极不稳定的问题,出现了多种化极方法．本文基于概率成像技术提出了一种等效物性的反演方法,实现对地下等效场源的反演成像,取得了对低纬度磁异常稳定化极的效果．化极反演中逐次对剩余异常进行反演成像,实现由概率模型到物性模型的复杂映射,避免了类似反演中需要大型方程组求解等问题,并将概率模型的构制、物性参数的反演和反演评价有机地集成到一起,加速了反演成像的进程,使反演成像速度与目前概率成像的计算速度达到可相比拟的程度．对理论模型和实际资料的计算表明,该方法对低纬度磁异常化极处理是稳定有效的,而且可以较好地压制噪声干扰,能够在噪声干扰条件下进行反演化极．     

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

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