Three-dimensional forward modeling for magnetotelluric sounding by finite element method

Abstract

A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell’s equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.

Keywords

magnetotelluric sounding / three-dimensional forward modeling / finite element method / general variation principle / divergence condition

References

[1] HohmannG. W.. Three-dimensional induced polarization and electromagnetic modeling [J]. Geophysics, 1975, 40(2): 309-324

[2] WannmakerP. E.. Advances in three-dimensional magnetotelluric modeling using integral equations [J]. Geophysics, 1991, 56(11): 1716-1728

[3] XiongZ.. EM modeling of three-dimensional structures by the method of system iteration using integral equations [J]. Geophysics, 1992, 57(12): 1556-1561

[4] BaoG.-s., ZhangB.-x., JingR.-z., MaoX.-jin.. Numerical modeling for 3-D EM responses using integral equation solution [J]. Journal of Central South University of Technology: Natural Science, 1999, 30(5): 472-474

[5] XuK.-j., LiT.-lin.. Three-dimensional magnetotelluric forward modeling using integral equation [J]. Northwestern Seismological Journal, 2006, 28(2): 104-107

[6] MackieR. L., MaddenT. R., WannamakerP. E.. Three-dimensional magnetotelluric modeling using difference equations-Theory and comparisons to integral equation solutions [J]. Geophysics, 1993, 58(2): 215-226

[7] XiongZ.. Domain decomposition for 3D electromagnetic modeling [J]. Earth Planets Space, 1999, 51(4): 1013-1018

[8] FomenkoE. Y., MogiT.. A new computation method for a staggered grid of 3D EM field conservative modeling [J]. Earth Planets Space, 2002, 54(7): 499-509

[9] TanH. D., YuQ. F., WeiW. B., BookerJ.. Magnetotelluric three-dimensional modeling using the staggered-grid finite difference method [J]. Chinese Journal of Geophysics, 2003, 46(5): 705-711

[10] ZhdanovM. S., SpichakV. V.Mathematical modeling of electromagnetic fields in three-dimensional inhomogeneous media [M], 1992, Moscow, Nauka Publishing House

[11] MehaneeS. A.Multidimensional finite difference electromagnetic modeling and inversion based on the balance method [D], 2003, Salt lake City, University of Utah

[12] MogiT.. Three-dimensional modeling of magnetotelluric data using finite-element method [J]. Journal of Applied Geophysics, 1996, 35(2): 185-189

[13] ZysermanF. I., SantosJ. E.. Parallel finite element algorithm with domain decomposition for three-dimensional magnetotelluric modeling [J]. Journal of Applied Geophysics, 2000, 44(4): 337-351

[14] HuangL.-p., DaiS.-kun.. Finite element calculation method of 3D electromagnetic field under complex condition [J]. Earth Science-Journal of China University of Geosciences, 2002, 27(6): 775-779

[15] MitsuhataY., UchidaT.. 3D magnetotelluric modeling using the T-Ω finite-element method [J]. Geophysics, 2004, 69(7): 108-119

[16] XuS.-zhe.Finite element method for geophysics [M], 1994, Beijing, Science Press

[17] SmithJ. T.. Conservative modeling of 3D electromagnetic fields; Part: II Biconjugate gradient solution and an accelerator [J]. Geophysics, 1996, 61(10): 1319-1324

[18] HaberE., AscherU. M., Aruliah, OlderburgD. W.. Fast simulation of 3D electromagnetic problems using potentials [J]. Journal of Computational Physics, 2000, 163(2): 150-171

[19] ChenL.-s., LiuR., WangT.-sheng.The magnetotelluric sounding data processing and interpretation methods [M], 1989, Beijing, Petroleum Industry Press

[20] ZhdanovM. S., VarentovI. M.. Methods for modeling electromagnetic fields: Results from COMMEMI-The international project on the comparison of modeling methods for electromagnetic induction [J]. Journal of Applied Geophysics, 1997, 37(3): 133-271