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

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

水平层状介质中基于DTA的三维电磁波逆散射快速模拟算法

为提高水平层状介质中三维电磁波散射和逆散射数值模拟的效率,在对角张量近似(DTA)的基础上根据不同回代方式得到了求解积分方程的DTA1和DTA2两种近似. 这两种近似可以作为计算积分方程稳定型双共轭梯度快速Fourier变换(BCGS-FFT)算法的初始猜测值和预条件因子,从而形成效率更高的混合DTA-BCGS算法. 散射实例说明了DTA2的高精度和混合DTA-BCGS算法尤其是混合DTA2-BCGS算法的高效率. 由于DTA2近似程度更高,将DTA2与变型Born迭代反演方法(DBIM)相结合形成了一种对三维异常体进行重构的快速电磁波逆散射技术. 文中的逆散射实例说明所开发的逆散射技术对重构水平层状介质中的任意三维异常体是非常有效的.     

Electromagnetic fields generated by finite‐length wire sources: comparison with point dipole solutions

In present‐day land and marine controlled‐source electromagnetic (CSEM) surveys, electromagnetic fields are commonly generated using wires that are hundreds of metres long. Nevertheless, simulations of CSEM data often approximate these sources as point dipoles. Although this is justified for sufficiently large source‐receiver distances, many real surveys include frequencies and distances at which the dipole approximation is inaccurate. For 1D layered media, electromagnetic (EM) fields for point dipole sources can be computed using well‐known quasi‐analytical solutions and fields for sources of finite length can be synthesized by superposing point dipole fields. However, the calculation of numerous point dipole fields is computationally expensive, requiring a large number of numerical integral evaluations. We combine a more efficient representation of finite‐length sources in terms of components related to the wire and its end points with very general expressions for EM fields in 1D layered media. We thus obtain a formulation that requires fewer numerical integrations than the superposition of dipole fields, permits source and receiver placement at any depth within the layer stack and can also easily be integrated into 3D modelling algorithms. Complex source geometries, such as wires bent due to surface obstructions, can be simulated by segmenting the wire and computing the responses for each segment separately. We first describe our finite‐length wire expressions and then present 1D and 3D examples of EM fields due to finite‐length sources for typical land and marine survey geometries and discuss differences to point dipole fields.     

埋地目标体矢量电磁散射的一种快速正演算法

利用积分方程方法以及半空间并矢格林函数的快速算法对埋地目标体矢量电磁散射进行正演计算。首先,利用半空间电磁并矢格林函数建立起埋地目标体的体积分方程。然后通过将空间偏导转移至格林函数谱域积分的积分号之外,并采用离散复镜像方法来近似余下的零阶索末菲积分,进一步得到并矢格林函数各个分量的闭合形式。由于避免了对索末菲积分的繁琐数值计算,使得生成反应矩阵和计算散射场时由半空间并矢格林函数计算带来的瓶颈问题得到较好的克服,因而极大地提高了埋地目标体电磁响应正演计算效率,同时也能保证足够的精度。     

瞬变电磁法正演计算进展

详细介绍了瞬变电磁法正演计算的方法、现状和发展趋势.瞬变电磁法一维正演计算需要将电磁场从频率域转换至时间域,转换方法有三种,分别是Gaver-Stehfest算法、余弦变换和Guptasarma算法.在这三种方法中,使用较多的是Gaver-Stehfest算法和余弦变换,Gaver-Stehfest算法速度较快,但精度不及余弦变换.瞬变电磁法的数值模拟主要集中于2.5维和三维,使用的数值计算方法有积分方程法、有限差分法、有限单元法和SLDM法.积分方程法主要在三维数值模拟中使用,现已很少使用;有限差分法和有限单元法是目前瞬变电磁法2.5维和三维数值模拟的主要方法;SLDM法主要应用于三维数值模拟.我国瞬变电磁法正演计算成果主要集中在回线源激发的瞬变电磁场一维数值计算和利用有限单元法进行2.5维和三维数值模拟.瞬变电磁法正演计算的发展趋势有:数值算法的改进、提高计算效率和研究地形对瞬变电磁场的影响规律.     

GPR modelling by the Fourier method: improvement of the algorithm

We improve two aspects of the modelling scheme for the simulation of electromagnetic radio waves, based on the Fourier pseudospectral method.
When there are large contrasts in the material properties, use of the standard algorithm (regular grid) causes a series of artefacts, as, for instance, ringing and acausal events. These problems, due to the non-locality of the differential operator, are solved by using the staggered Fourier method (staggered grid).
Realistic radiation patterns can be obtained from simple combinations of magnetic and electric sources. If the directivity pattern of the antenna is known, from either a finite-difference simulation or an analytic evaluation or an experimental characterization, it can then be simulated by a composite-source concept. This effective source is implemented in the modelling algorithm by means of a perturbation technique, which first computes the intensity and directional spectra of the single electromagnetic sources. Their location is optimized to obtain the best fit with a minimum number of sources. The approach is, in principle, valid for the far-field radiation pattern of the antenna.     

复杂地下异常体的可控源电磁法积分方程正演

可控源电磁法具有分辨率高及抗干扰能力强等特点,是一种重要的地电磁勘探方法.目前,可控源电磁法的高精度正演计算一直是其核心研究问题之一.传统积分方程法一般采用近似积分公式、简单矩形网格和近似的奇异性体积分计算技术,制约了体积分方程法处理复杂地下异常体的能力,降低了计算精度.针对上述问题,本文基于完全积分公式、四面体非结构化网格和奇异体积分的精确解析解来高精度求解复杂可控源电磁模型的正演响应.首先,从电场积分公式出发,推导了可控源电磁问题满足的积分方程;其次,借助于非结构化四面体网格离散技术,实现了地下复杂异常体的有效模拟.最后,利用散度定理把强奇异值体积分转换为一系列弱奇异性的面积分公式,并通过推导获得了这些弱奇异性的面积分公式的解析解,从而最终实现三维可控源电磁问题的高精度积分求解.以块状低阻体地电模型为测试模型,采用本文提出的积分方程方法获得的数值解与其他公开数值算法解进行对比分析,其对比结果具有高度的吻合性,验证了算法的正确性;同时,设计了球状及复杂地电模型进行算法收敛性测试,进一步验证算法的正确性以及能够处理地下复杂模型的能力.     

基于改进的接收点插值算法的频率域海洋可控源电磁法2.5维正演

在海洋可控源电磁法勘探中,接收站常置于海底.在进行海洋电磁场模拟时,由于海水和海底介质存在显著电性差异,这给海底接收点处场值的求取带来困难.本文提出一种新的接收点插值算法,该算法考虑到海底电场法向分量不连续性问题,用法向电流分量进行插值以准确求取海底任意接收点处电磁场值.本文利用交错网格有限差分法实现了二维介质中频率域海洋可控源法(CSEM)正演.对构造走向做傅里叶变换,将三维电磁模拟问题转换为波数域2.5维问题,即三维场源激励下针对二维地电模型的电磁模拟问题.使用交错网格有限差分法,基于一次场/二次场分离方法导出波数域二次电场离散形式,并进一步求得波数域电磁场.采用本文提出的改进的插值算法可求得海底任意接收点处波数域电磁场,采用傅里叶逆变换对波数域电磁场进行积分可得到接收点处空间域电磁场.模型算例表明,与常规的线性插值和严格插值算法相比,本文提出的改进的插值算法具有更高的精度.     

层状介质高频电磁场计算方法及结果分析

本文在研究了高密度采样线性滤波算法计算层状介质高频电磁场的电磁响应的基础上,系统比较了Anderson等数值计算方法的速度和精度,给出了不同的地电参数和接收距对该计算方法精度的影响.通过比较表明,高密度采样算法不仅速度快,而且精度高.本文还对水平回线源层状介质当发射线圈位于地面时的表达式的第二项积分进行了公式推导,为电磁场的计算提供了清晰的思路.     

Three‐dimensional modelling of controlled‐source electromagnetic surveys using an edge finite‐element method with a direct solver

A fully three‐dimensional finite‐element algorithm has been developed for simulating controlled‐source electromagnetic surveys. To exploit the advantages of geometric flexibility, frequency‐domain Maxwell's equations of the secondary electric field were discretised using edge‐based finite elements while the primary field was calculated analytically for a horizontally layered‐earth model. The resulting system of equations for the secondary field was solved using a parallel version of direct solvers. The accuracy of the algorithm was successfully verified by comparisons with integral‐equations and iterative solutions, and the applicability to models containing large conductivity contrasts was verified against published data. The advantages of geometry‐conforming meshes have been demonstrated by comparing different mesh systems to simulate an inclined sheet model. A comparison of the performance between direct and iterative solvers demonstrated the superior efficiency of direct solvers, particularly for multisource problems.     

Electromagnetic response of a buried cylindrical conductor for sheet and line current sources

A major problem in electromagnetic induction studies in regions of localized source fields, such as the auroral and equatorial electrojet regions, is the source effect. Using an analytical model, the electromagnetic response of a buried conducting cylinder to sheet current and line current excitations has been studied for the period rangeT=5 s to 24 h. The validity of the numerical results obtained from the analytical model are compared with the numerical results obtained from a finite difference model. The results show that for periods less than 30 min, there is no significant difference in the response of the cylinder to both source fields. However, significant differences are observed at longer periods. It was also observed that the equivalent height at which a uniform sheet current at 100 km above the earth's surface can be approximated by a line current varies as a function of the source period.     

基于电场矢量波动方程的三维可控源电磁法有限单元法数值模拟

从可控源电磁法的基本原理出发,推导了基于电场矢量波动方程的三维边值问题,利用广义变分原理,把边值问题转换为变分问题,并引入散度条件,避免了伪解的出现,使有限元计算在理论上更加完备.在准静态近似条件下,把水平电偶极子在空中和大地的远区电场闭合表达式作为有限元计算中的区域外边界条件,解决了边界条件加载的困难;把应用于地震模拟中的伪delta函数引入到可控源电磁法中的三维有限元模拟中,消除了源点的奇异性,提高了方程组的稳定性.通过对均匀大地和层状介质模型的模拟,检验了程序的正确性,并对典型的地质体模型进行了数值模拟,分析了其变化规律.     

基于电场总场矢量有限元法的接地长导线源三维正演

电磁场数值模拟的背景场/异常场算法是三维正演的有效策略之一,优点为采用解析法计算电磁场背景场代替场源项、克服了场源奇异性,缺点为不适用于发射源布置于起伏地表或背景模型复杂的情形.总场算法是直接对电磁场总场开展数值模拟,其难点是有效加载场源、保证近区与过渡区数值解精度.本文以水平电偶源形式分段加载接地长导线源,并以电场总场Helmholtz方程为矢量有限元法控制方程,实现了基于非结构化四面体网格剖分的接地长导线源频率域电磁法三维正演.通过与均匀全空间中水平电偶源产生的电场解析解对比,验证了本文算法的正确性,并分析了四面体外接圆半径与其最短棱边的最大比值和四面体二面角最小值对数值解精度的影响规律.通过与块状高导体地电模型的积分方程法、有限体积法和基于磁矢量势Helmholtz方程的有限元法数值解对比,进一步验证了本文算法正确性,同时说明了非结构化四面体网格能够更加精细地剖分电性异常体,利于获得精确数值解.     

Numerical and analogue model results for electromagnetic induction for an island situated near a coastline

The behaviour of time-varying electromagnetic fields near an island situated in a shallow ocean is investigated using both a three-dimensional finite-difference numerical method and a scaled analogue model method. The effect of a coastline located at some distance from the island is included in the study. The numerical model results and the scale model results are compared for various traverses across the island. The results indicate a high degree of compatibility between the two methods for studying problems involving three-dimensional conductivity structures.     

A comparison of numerical,analogue model,and field-station vertical magnetic-fields for the Vancouver Island region

Vertical magnetic fields for a three-dimensional numerical model, for a laboratory-analogue model, and from field stations for the Vancouver Island region of British Columbia, Canada, are compared. The numerical results are obtained using a three-dimensional finite-difference numerical technique employing a 25 × 25 × 25 mesh of grid points for a simplified mathematical model of the Vancouver Island region. The calculations are carried out for a source frequency of 0.004 Hz. The analogue model results for four traverses over the Island model and the field station values (obtained from transfer function analysis) for ten locations are those discussed previously by Nienaber et al. (1979a, b). General agreement exists between the numerical, analogue, and field station data, and comparison of results between these methods is important in three-dimensional electromagnetic induction studies of complex geomagnetic induction problems.     

长偶极大功率可控源电磁波响应特征研究

地球物理学中关于电磁波勘探研究通常采用的是地球半空间模型.然而,对于几十公里的有限长电缆源（长偶极源）,远距离电磁波场探测必须要考虑电离层的影响,它是一个全空间问题.关于包含电离层、空气层和地球介质（我们称“地-电离层”模式）的电磁波场特征的研究在国外较少,国内几乎是空白.本文采用全空间积分方程法首先对小尺度的可控源电磁波场特征进行了研究,由于此时电离层的影响可忽略,它应该和半空间成熟的CSAMT模拟结果一致,对比结果表明,二者是一致的,验证了全空间模拟方法的可靠性和有效性.随后进行了50 km长电缆电离层和空气层高度都为100 km的“地-电离层”模式大尺度电磁波场模拟,以探讨大尺度可控源电磁波场的特征.给定频率的“地-电离层”模式电磁场的衰减曲线表明长电缆远距离电磁波场由于受电离层的作用存在衰减逐渐变小的过渡场和衰减变小的波导场.为了探讨复杂介质“地-电离层”模式电磁波特征,对“地-电离层”模式的典型地盾和地台多层介质模型进行了数值模拟,得到了偶极源长度50 km、电流200A、收发距离远达1600 km和2500 km的合理的电磁场结果.最后,对一简单含油储层结构模型进行了长偶极、大功率、远距离电磁波场响应计算.储层横向不均匀复杂结构模拟的结果表明,考虑电离层和大气层的“地-电离层”模式大尺度深层复杂介质模拟时,电磁场对深部目标体仍有很好的异常响应.     

3D Modelling of the electromagnetic response of geophysical targets using the FDTD method1

A publicly available and maintained electromagnetic finite-difference time domain (FDTD) code has been applied to the forward modelling of the response of 1D, 2D and 3D geophysical targets to a vertical magnetic dipole excitation. The FDTD method is used to analyse target responses in the 1 MHz to 100MHz range, where either conduction or displacement currents may have the controlling role. The response of the geophysical target to the excitation is presented as changes in the magnetic field ellipticity. The results of the FDTD code compare favourably with previously published integral equation solutions of the response of 1D targets, and FDTD models calculated with different finite-difference cell sizes are compared to find the effect of model discretization on the solution. The discretization errors, calculated as absolute error in ellipticity, are presented for the different ground geometry models considered, and are, for the most part, below 10% of the integral equation solutions. Finally, the FDTD code is used to calculate the magnetic ellipticity response of a 2D survey and a 3D sounding of complicated geophysical targets. The response of these 2D and 3D targets are too complicated to be verified with integral equation solutions, but show the proper low- and high-frequency responses.     

A numerical study of the topographic effect on electromagnetic fields in a three-dimensional conductivity model

Summary The topographic effect on electromagnetic fields in a three-dimensional conductivity model is investigated using a finite-difference numerical method. The model is that of a conducting block buried in the host Earth beneath a conducting hill. The altitude effect as well as the effect of conductivity of the hill has been studied by comparing the results for the field components at the surface of the hill with those at the surface of a flat half-space. The results indicate that both altitude and the hill conductivity influence the behaviour of the electromagnetic fields. For certain traverses beyond the edge of the block, it is seen that the variations in some field components are mainly due to the topographic effect and that the contribution from the block itself is negligible.