Three-dimensional land FD-CSEM forward modeling using edge finite-element method

Abstract

A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys, based on a finite-element discretization of the Helmholtz equation for the electric fields, has been developed. The main difference between our modeling method and those previous works is edge finite-element approach applied to solving the three-dimensional land frequency-domain electromagnetic responses generated by horizontal electric dipole source. Firstly, the edge finite-element equation is formulated through the Galerkin method based on Helmholtz equation of the electric fields. Secondly, in order to check the validity of the modeling code, the numerical results are compared with the analytical solutions for a homogeneous half-space model. Finally, other three models are simulated with three-dimensional electromagnetic responses. The results indicate that the method can be applied for solving three-dimensional electromagnetic responses. The algorithm has been demonstrated, which can be effective to modeling the complex geo-electrical structures. This efficient algorithm will help to study the distribution laws of 3-D land frequency-domain controlled-source electromagnetic responses and to setup basis for research of three-dimensional inversion.

Keywords

three-dimensional model / frequency-domain electromagnetic method / horizontal electric dipole / forward modeling / edge finite-element

References

[1] WhiteD, BoernerD, WuJ, LucasS, BerrerE, HannilaJ, SomervilleR. Mineral exploration in the Thompson nickel belt, Manitoba, Canada, using seismic and controlled-source EM methods [J]. Geophysics, 2000, 65(6): 1871-1881

[2] ConstableS, SrnkaL J. Special section-marine controlled-source electromagnetic methods an introduction to marine controlled-source electromagnetic methods for hydrocarbon exploration [J]. Geophysics, 2007, 72(2): WA3-WA12

[3] MacgregorL, SinhaM, ConstableS. Electrical resistivity structure of the Valu Fa Ridge, Lau Basin, from marine controlled-source electromagnetic sounding [J]. Geophysical Journal International, 2001, 146(1): 247-236

[4] AvdeevD B, KuvshinovA V, PankratovO V, NewmanG A. Three-dimensional induction logging problems, Part 1: An integral equation solution and model comparisons [J]. Geophysics, 2002, 67(2): 413-426

[5] HursanG, ZhdanovM S. Contraction integral equation method in 3-D electromagnetic modeling [J]. Radio Science, 2002, 37(6): 1089

[6] LiuJ-x, JiangP-f, TongX-z, XuL-h, XieW, WangHao. Application of BICGSTAB algorithm with incomplete LU decomposition preconditioning to two-dimensional magnetotelluric forward modeling [J]. Journal of Central South University: Science and Technology, 2009, 404485-490

[7] ZhdanovM S, LeeS K, YoshiokaK. Integral equation method for 3Dmodeling of electromagnetic fields in complex structures with inhomogeneous background conductivity [J]. Geophysics, 2006, 71(6): G333-G345

[8] HaberE, AscherU, AruliahD, OldenburgD. Fast simulation of 3D electromagnetic problems using potentials [J]. Journal of Computational Physics, 2000, 163(1): 150-171

[9] HanN, NamM J, KimH J, SongY, SuhJ H. A comparison of accuracy and computation time of threedimensional magnetotelluric modelling algorithms [J]. Journal of Geophysics and Engineering, 2009, 6(2): 136-145

[10] MehaneeS, ZhdanovM. A quasi-analytical boundary condition for three-dimensional finite difference electromagnetic modeling [J]. Radio Science, 2004, 39(6): RS6014

[11] NewmanG A, AlumbaughD L. Frequency-domain modeling of airborne electromagnetic responses using staggered finite differences [J]. Geophysical Prospecting, 1995, 4381021-1042

[12] StreichR. 3D finite-difference frequency-domain modeling of controlled-source electromagnetic data: Direct solution and optimization for high accuracy [J]. Geophysics, 2009, 74(5): F95-F105

[13] WeissC J, ConstableS. Mapping thin resistors and hydrocarbons with marine EM methods, part II-modeling and analysis in 3D [J]. Geophysics, 2006, 71(6): G321-G332

[14] BörnerR U, ErnstO G, SpitzerK. Fast 3-D simulation of transient electromagnetic fields by model reduction in the frequency domain using Krylov subspace projection [J]. Geophysical Journal International, 2008, 173(3): 766-780

[15] da SilvaN V, MorganJ V, MacgregorL, WarnerM. A finite element mutifrontal method for 3D CSEM modeling in the frequency domain [J]. Geophysics, 2012, 77(2): E101-E115

[16] MukherjeeS, EverettM E. 3Dcontrolled-source electromagnetic edge-based finite element modeling of conductive and permeable heterogeneities [J]. Geophysics, 2011, 76(4): F215-F226

[17] NamM J, KimH J, SongY, LeeT J, SonJ S, SuhJ H. 3D magnetotelluric modelling including surface topography [J]. Geophysical Prospecting, 2007, 55(2): 277-287

[18] UmE, HarrisJ M, AlumbaughD L. 3D time-domain simulation of electromagnetic diffusion phenomena: A finiteelement electric-field approach [J]. Geophysics, 2010, 75(4): F115-F126

[19] FarquharsonC G, MiensopustM P. Threedimensional finite-element modelling of magnetotelluric data with a divergence correction [J]. Journal of Applied Geophysics, 2011, 75(4): 699-710

[20] JinJ MThe finite element method in electromagnetics [M], 2002, New York, John Wiley & Sons

[21] SaadYIterative methods for sparse linear systems [M], 2002, London, Cambridge University Press

[22] KaufmanA A, KellerG VFrequency and Transient Soundings (Methods in Geochemistry and Geophysics) [M], 1983

[23] CaiH-z, XiongB, MichaelZ. Threedimensional marine controlled-source electromagnetic modeling in anisotropic medium using finite element method [J]. Chinese Journal of Geophysics, 2015, 58(8): 2839-2850