Three‐Dimensional Magnetotelluric Modeling using the Staggered‐Grid Finite Difference Method
Three‐Dimensional Magnetotelluric Modeling using the Staggered‐Grid Finite Difference Method
复制标题
DOI：
10.1002/cjg2.420
复制
发表时间：
2003-09
期刊：
影响因子：
--
通讯作者：
H. Tan;Qinfan Yu;Booker John;Wenbo Wei
中科院分区：
文献类型：
--
作者：
H. Tan;Qinfan Yu;Booker John;Wenbo Wei
本文详细介绍了交错网格有限差分法三维正演建模的关键问题。它们包括交错网格、麦克斯韦方程积分形式的离散化、边界条件、求解线性代数方程和计算三维张量阻抗。由于给出了更明确的边界条件，并采用双共轭梯度稳定法求解大系数矩阵的线性代数方程，得到了快速有效地计算整个空间内高精度电场和磁场的算法。通过比较3 - D正演模型解与邻近四分之一空间的解析解，以及2 - D有限元法对2 - D棱镜的2 - D正演模型，证明了这一点。该算法为三维反演的研究奠定了基础。
The crucial problems of 3‐D forward modeling using the staggered‐grid finite difference method are described in detail in this paper. They include staggered‐grid, discretization of integral form of Maxwell equation, boundary condition, solving linear algebra equations and calculating 3‐D tensor impedance. Because of giving the more explicit boundary conditions and using Bi‐conjugate gradients stabilized method to solve the linear algebra equation with large coefficient matrix, we get the fast algorithm to calculate effectively the high‐precision electrical and magnetic fields in the whole space. This has been proved by comparing the 3‐D forward modeling solutions with analytic solution to abutting quarter‐spaces, and 2‐D forward modeling to 2‐D prism using 2‐D finite element method. This efficient algorithm has setup basis for doing research on 3‐D inversion.