考虑关断时间的回线源激发TEM三维时域有限差分正演

从麦克斯韦旋度方程出发可以直接导出瞬变电磁场扩散方程,然而扩散方程不含电场对时间的一阶导数,不能构成显式的时域有限差分方程,借鉴du Fort-Frankel有限差分离散方法引入虚拟位移电流项构建显式时域有限差分方程.对Wang和Hohmann的经典时域算法进行了两点改进:第一,通过将矩形回线源电流密度加入麦克斯韦方程组的安培环路定理方程,实现回线源瞬变电磁激发源加入;第二,在计算中考虑关断时间.第一点改进使时域有限差分方程考虑了一次场的计算,并且源的计算不再依赖均匀半空间模型响应作为初始条件,使算法能够适应表层电阻率不均匀时的三维复杂模型.由于实际观测中不可能出现阶跃电流的关断形式,第二点改进可以方便设置发射电流下降沿.采用改进的三维时域有限差分正演算法对均匀半空间模型、四类三层模型、均匀半空间中含有低阻块体模型进行了计算并分别与解析解、线性数字滤波解、积分方程解和Wang的三维时域有限差分解进行了对比验证.以H模型为例,采用建立的三维时域有限差分正演算法计算了不同关断时间的斜阶跃脉冲回线源瞬变电磁中心点感应电动势衰减曲线.以实际地质资料为基础,构建包含两层采空区的三维复杂模型,以1 μs的极短关断时间进行了复杂模型定回线源瞬变电磁响应计算,并计算了该复杂模型的视电阻率曲线.     

INTERPRETATION AND DEPTH OF INVESTIGATION OF GRADIENT MEASUREMENTS IN DIRECT CURRENT GEOELECTRICS*

Gradient measurements in a homogeneous electrical primary field can easily be interpreted for simple models. The simplified solution (conducting or resistant body in a homogeneous space in a homogeneous electrical field) is often sufficiently accurate, as comparisons with the exact solution (body of finite resistivity in a homogeneous half-space in a quasihomogeneous electrical field) show. The exact geometry of the body cannot be determined by gradient measurements; the same anomaly of apparent resistivity can be caused by different bodies. In particular, the similarity between a sphere and a cube of the same volume is very high. There is a distinct influence of the resistivity of the overburden: the higher this resistivity is, the stronger is the effect caused by a buried body. If a deviation of 10% of the apparent resistivity is assumed as the lower boundary at which a buried body can be detected by gradient measurements, the depth of investigation for a three-dimensional body is approximately equal to its width; in the two-dimensional case the thickness of the overburden can be twice the width. If the overburden has a resistivity which is higher than the resistivity of the substratum, these depths are greater. The greatest possible depth is approximately three times the width of the body.     

地-井瞬变电磁响应特征数值模拟分析

井中瞬变电磁法(Bore-hole transient electromagnetic method-BHTEM)是指接收线圈在钻井中观测瞬变场响应用以勘查深部矿产资源的勘探方法,其中以地-井(地面激发井中接收)组合方式研究最多、应用最广.本文应用时域有限差分法(FDTD),建立包含薄板导体的均匀半空间二维数学模型,采用线源为激发源,选用Mur吸收边界条件,对矩形回线源在半空间中产生的瞬变电磁场进行数值模拟,计算了低阻板状导体在均质半空间和有低阻覆盖层影响情况下的地-井瞬变电磁异常响应,并对响应的特征及规律进行研究分析,为研究地-井TEM提供参考.     

Surface NMR applied to an electroconductive medium1

An integrated electromagnetic and nuclear magnetic resonance (NMR) method is proposed for investigating highly conductive soil areas. Using a simple model of a homogeneous half-space and the same antenna for both methods, a significant improvement in the NMR data interpretation is obtained. A case study has shown fair agreement between the results from computer modelling, field tests, and data from a nearby observation well. The electromagnetic method and the half-space model were selected for easy integration into an existing instrument used for the NMR method. A more accurate knowledge of the conductivity distribution with depth will further improve the final result.     

Dynamic soil-tunnel interaction in layered half-space for incident plane SH waves

The dynamic soil-tunnel interaction is studied by indirect boundary element method (IBEM), using the model of a rigid tunnel in layered half-space, which is simplified to a single soil layer on elastic bedrock, subjected to incident plane SH waves. The accuracy of the results is verified through comparison with the analytical solution. It is shown that soil-tunnel interaction in layered half-space is larger than that in homogeneous half-space and this interaction mechanism is essentially different from that of soil-foundation-superstructure interaction.     

基于电场Helmholtz方程的回线源瞬变电磁法三维正演

正演是电磁法勘探野外工作参数选取、室内资料处理与解释的基础,精确、稳定、高效的三维正演算法尤为重要.本文采取先求解拉普拉斯域电场、再由Gaver-Stehfest算法获得时间域磁场的思路,基于电场异常场Helmholtz方程实现了交错网格有限差分法和有限体积法对回线源瞬变电磁法的三维正演.通过对比低阻块状体的积分方程法、时域有限差分法、矢量有限单元法和SLDM法的数值解,验证了交错网格有限差分法和有限体积法的正确性.由于交错网格有限差分法、有限体积法和基于矩形块单元的矢量有限单元法将待求电场均定义在矩形块单元棱边上,因此三种数值算法可采用相同方法进行电场待求量编码、计算背景场和后处理.然而,与矢量有限单元法相比,交错网格有限差分法和有限体积法的系数矩阵更加稀疏,求解效率更高.通过对水平低阻板状体三维模型的数值模拟,我们发现本研究中交错网格有限差分法比有限体积法精度更高;再利用一维解析法求解相应三层层状地电模型的感应电动势,我们还发现两种数值算法和一维解析法计算的感应电动势等值线形状吻合程度高,只是数值范围略有差异.     

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

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

A finite-element time-domain forward-modelling algorithm for transient electromagnetics excited by grounded-wire sources

A three-dimensional finite-element time-domain forward-modelling algorithm is developed to simulate transient electromagnetics excited by grounded-wire sources. The main advantage of this finite-element time-domain algorithm is that full transmitting-current waveforms and complex-shaped sources resulting from topography can be directly dealt with in this algorithm. The models used to test this algorithm include a homogeneous half-space model, a stratified-medium model, the model of a complex conductor at a vertical contact and the Ovoid Zone massive sulfide deposit at Voisey's Bay, Canada. The homogeneous half-space model is used to determine the truncation boundary for a computational domain, and to compare with the electromagnetic responses excited by step-off, step-on and direct current waveforms. For the stratified-medium model, results demonstrate that full transmitting waveforms have strong effects on the observed electromagnetic responses. The model of a complex conductor at a vertical contact is designed for the grounded electrical source airborne transient electromagnetic method and is also used to examine the effectiveness of the broadside and inline configurations for such a vertical, thin plate embedded in the subsurface. The area of the Ovoid Zone massive sulfide deposit possesses non-negligible topography, the effects of which on the shapes of the grounded-wire sources must be taken into account when implementing the finite-element time-domain solver. The results show that both the broadside and inline electromagnetic responses are strongly affected by the massive conductive ore body.     

深海热液硫化物矿体3D瞬变电磁正演

深海热液硫化物矿体瞬变电磁的正演是考虑深海环境的全空间条件下三维体的涡流电磁响应.采用全空间矢量有限元法模拟计算深海热液硫化物矿的三维瞬变电磁响应,对硫化物矿体采用矩形单元模型剖分,应用Galerkin法推导有限元方程,先计算频率域响应,再通过Fourier反变换将其转换至时间域,得出深海热液硫化物矿矿体的瞬变电磁响应.并用双半空间模型的解析解检验了全空间矢量有限元法模拟计算算法和程序的正确性,最后按照等比例缩小电磁物理实验原则,比对数值计算和物理实验结果论证了全空间3D模型数值的正确性.结果表明:对于海水、矿体以及围岩复杂电磁边界,应用全空间矢量有限元法模拟计算深海热液硫化物矿瞬变电磁响应异常与物理模拟结果一致,而且计算方法简单精确,异常幅值明显,边界清晰.     

多源频率域地空系统三维电磁响应分析

地空电磁法已经成为深部资源勘探的重要地球物理方法,但对频率域地空系统的三维多源电磁响应特征研究较少.本文设计了多种激励源组合方式,采用非结构化有限元数值模拟方法,对三维地电模型的空中垂直磁场的响应特征进行了研究.首先推导了基于电场的双旋度公式及其变分形式,加入罚项以减少伪解的影响.接着把有限元稀疏矩阵方程转换为频率的函数,采用Krylov子空间投影方法,通过模型降阶算法降低稀疏矩阵的阶数,实现多频点的快速计算.建立了三维低阻体模型、高阻体模型以及两个相邻低阻体模型,分别采用单源、双源、三源和四源激励模式,从垂直磁场的总场、二次场响应和全域视电阻率等方面进行分析比较.结果表明:多源地空电磁法不仅可以增加总场的强度,而且可以改变异常体的二次电磁响应分布规律.各电偶源延长线呈正三角形分布的三源和矩形分布的四源激励模式在增强信号强度以及削弱异常体的边界效应方面具有一定的优势,是一种优化的多源激励方式.     

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.     

Impedances of embedded rigid strip foundations

This paper is concerned with the dynamic response of rigid strip foundations of arbitrary geometry embedded in a homogeneous elastic half-space. The embedded rigid foundation is modelled by an equivalent domain in a uniform half-space which is subjected to an appropriate body force field. The components of the impedance matrix are determined through the solution of a linear simultaneous equation system which is established by invoking rigid body displacements of discrete locations within the equivalent domain and appropriate equilibrium consideration. It is found that high numerical efficiency and flexibility can be achieved using the body force model when compared to boundary integral formulations through the selection of appropriate displacement influence functions and a ‘parent domain’ in the analysis. Numerical results are presented to illustrate the influence of the embedment ratio, frequency of excitation, foundation geometry and Poisson's ratio on the vertical, horizontal, rocking and coupled impedances of a single embedded foundation. The effect on the impedance due to the presence of an adjacent embedment is investigated for various distances between foundations and embedment ratios.     

矿井瞬变电磁法全空间视电阻率解释方法研究

从全空间中心回线源响应公式出发,推导了全空间晚期和全区视电阻率公式;通过与半空间公式比较和对实测资料的分析,研究了全空间视电阻率公式的应用效果.结果显示,全空间视电阻率为半空间视电阻率的(5/2)^2/3倍,其晚期视电阻率与真值的相对误差小于0.56%,而相应的半空间视电阻率与真值误差均大于42%.与常规解释方法相比,采用全空间全区视电阻率换算方法能更好地圈定低阻异常范围,且电阻率差异大.因此,该公式更适用于矿井瞬变电磁法资料处理与解释,并具有更好的效果.     

Scattering of SV waves by a canyon in a fluid-saturated, poroelastic layered half-space, modeled using the indirect boundary element method

The scattering of SV waves by a canyon in a fluid-saturated, poroelastic layered half-space is modeled using the indirect boundary element method in the frequency domain. The free-field responses are calculated to determine the displacements and stresses at the surface of the canyon, and fictitious distributed loads are then applied at the surface of the canyon in the free field to calculate the Green's functions for displacements and stresses. The amplitudes of the fictitious distributed loads are determined from the boundary conditions, and the displacements arising from the waves in the free field and from the fictitious distributed loads are summed to obtain the solution. The effects of fluid saturation, boundary conditions, porosity, and soil layers on the surface displacement amplitudes and phase shifts are discussed, and some useful conclusions are obtained. It is shown that the surface displacement amplitudes due to saturation and boundary conditions, different porosities, or the presence of a soil layer can be very dissimilar, and large phase shifts can be observed. The resulting wavelengths for an undrained saturated poroelastic medium are slightly longer than those for a drained saturated poroelastic medium; and are longer for a drained saturated poroelastic medium than those for a dry poroelastic medium. As porosity increases, the wavelengths become longer; and a layered half-space produces longer wavelengths than a homogeneous half-space.     

任意方向震源二维兰姆问题的一个精确解

均匀弹性半空间表面或内部震源产生的地震波场的解析解属于Lamb问题,采用Cagniard-deHoop方法,对与水平面呈任意夹角的表面线源,求解了其作用于弹性半空间时的拉普拉斯-傅里叶双积分变换解.以δ-脉冲函数为例,给出了任意方向作用力下波场的构成,并定量解出了P波、S波、首波和Rayleigh波等各波的位移表达式,分析了不同作用方向下各波位移的相对大小.建立数值模型并进行数值模拟,模拟结果验证了理论研究的正确性.研究成果为近地表地震波场的研究提供了理论依据.     

A study on the discrete image method for calculation of transient electromagnetic fields in geological media

We conducted a study on the numerical calculation and response analysis of a transient electromagnetic field generated by a ground source in geological media. One solution method, the traditional discrete image method, involves complex operation, and its digital filtering algorithm requires a large number of calculations. To solve these problems, we proposed an improved discrete image method, where the following are realized: the real number of the electromagnetic field solution based on the Gaver–Stehfest algorithm for approximate inversion, the exponential approximation of the objective kernel function using the Prony method, the transient electromagnetic field according to discrete image theory, and closed-form solution of the approximate coefficients. To verify the method, we tentatively calculated the transient electromagnetic field in a homogeneous model and compared it with the results obtained from the Hankel transform digital filtering method. The results show that the method has considerable accuracy and good applicability. We then used this method to calculate the transient electromagnetic field generated by a ground magnetic dipole source in a typical geoelectric model and analyzed the horizontal component response of the induced magnetic field obtained from the “ground excitation–stratum measurement” method. We reached the conclusion that the horizontal component response of a transient field is related to the geoelectric structure, observation time, spatial location, and others. The horizontal component response of the induced magnetic field reflects the eddy current field distribution and its vertical gradient variation. During the detection of abnormal objects, positions with a zero or comparatively large offset were selected for the drillhole measurements or a comparatively long observation delay was adopted to reduce the influence of the ambient field on the survey results. The discrete image method and forward calculation results in this paper can be used as references for relevant research.     

Dynamic stiffness of rigid rectangular foundations on the half-space

The dynamic soil–structure interaction of a rigid rectangular foundation with the subsoil represents a mixed-boundary value problem. This problem is formulated in terms of a system of coupled Fredholm integral equations of the first kind. The subsoil is modelled by a homogeneous, linear-elastic and isotropic half-space which is perfectly bonded to the rigid, rectangular foundation. An approximate solution for the resultant loads between the foundation and the half-space due to a unit forced displacement or rotation is obtained using the Bubnov–Galerkin method. Using this method the displacement boundary value conditions are exactly satisfied and the contact stress distributions between the foundation and the half-space are approximated by series expansions of Chebyshev polynomials. This method provides a simple means of studying the soil-structure interaction of rectangular foundations with different inertia properties.     

基于神经网络的视电阻率快速算法

本文从瞬变电磁均匀半空间二次磁场响应公式出发,提出了一种基于神经网络的视电阻率快速计算方法.以中心回线为例,根据瞬变响应公式的特点,简化网络结构,选用三层BP神经网络和误差训练算法,用均匀半空间样本数据进行训练,确定了收敛快、误差小的一步正割法和隐含单元数,得到基于不同采样时窗的一组网络参数.用本文方法与二分法、牛顿迭代法做模型计算比较,及最后的实验计算,说明算法的快速,准确.本文方法不依赖初始模型,避开了复杂的电磁场数值计算,实现了视电阻率的快速计算,对瞬变电磁法资料的快速解释有一定的参考价值.     

瞬变电磁法深度探测正演模拟研究

针对当前瞬变电磁法探测能力有限的问题,利用一次磁场理论公式,分析相同磁矩下、不同发射线圈一次磁场分布及衰减规律,从而确定其信号深度影响范围;对水平层状介质相同磁矩下、不同发射线圈的二次场进行数值模拟,采用Guptasarma和Singh滤波算法及余弦变换法,利用均匀半空间解析公式进行验证,并对不同正演模型结果进行分析,总结发射磁矩决定探测深度大小。对磁矩相同条件下产生的结果数据进行拟合,得出有效探测深度与线圈边长2倍相当的结论,可为野外矿井生产与工程勘探提供施工布置与定量解释的指导依据,从而提高野外工作效率。