时间谱电阻率法的二维正演算法

时间谱电阻率(TSR)法的二维正演, 是计算可极化二维地电构造上三维电流偶极源的电场瞬变响应, 属所谓2.5维时间域电磁场数值模拟问题, 是目前国际上未妥善解决的计算地球物理疑难问题.针对现有算法的局限性, 建立了新的算法, 其特点是: (1)采用二维有限单元算法, 在矩形网格中增加两对角线形成三角网格剖分, 同时用高斯消元法消除矩形网格中心结点的待求未知量.这样, 既可较准确地模拟任意二维复杂地电断面, 又可节省计算量.(2)采用直接计算二次场的新算法, 只需计算电场和磁场沿地电构造走向两个分量的一次场, 因而, 不但计算精度较高, 而且不显著增加计算量.(3)采用G -S变换法作逆拉氏变换, 并利用拉氏变换延迟定理在倍增的时间间隔中插值, 从而实现对密集采样时间瞬变过程的快速计算.(4)能对可极化和导电大地(即同时包括IP和EM效应). 

2-DIMENSIONAL FORWARD ALGORITHM FOR TIME SPECTRAL RESISTIVITY

The 2 dimensional forward calculation for time spectral resistivity (TSR) is to compute the electric field transient response caused by a 3 D electric current dipole on the 2 D polarizable earth surface. As part of the so called 2.5 D numerical simulation for time domain electromagnetic field, this 2 dimensional forward calculation is still a difficult problem in geophysics. In this paper, a new algorithm is established for this problem with the following four features: (1) A 2 D finite element method is adopted to subdivide a triangular lattice by adding two diagonal lines in the rectangular lattice, and at the same time to eliminate unknown variables to be calculated at the central node of the rectangular lattice by using the Gauss elimination. Therefore, this method can be used not only to simulate accurately any complex 2 D geoelectric section, but also to save the calculation load. (2) A new algorithm for the computation of the secondary field is directly employed to compute only a primary field of the two components of electric and magnetic fields along the geoelectric structure. Such a method is advantageous not only in the relatively high computation precision, but also in the absence of much increase in calculation load. (3) The new G S transformation method is used to make the inverse Laplace transformation, and, at the same time, the delay theorem of the Laplace transformation is used to calculate the interpolation at the intervals of the multiplication time, resulting in the rapid calculation of the transient electromagnetic responses to the dense samples. (4) This algorithm can be used to compute the transient electromagnetic response to the polarizable and conductive earth (including the IP and EM effects).