An edge-based finite element method for 3D marine controlled-source electromagnetic forward modeling with a new type of second-order tetrahedral edge element

Abstract

This paper presents a second order edge-based finite element method based on a new type of second order edge element for numerical modeling of 3D controlled source electromagnetic data in an anisotropic conductive medium. We adopt the anomalous potential field formulation to avoid the effect of source’s singularity. The computational domain was subdivided into unstructured tetrahedral elements. The sparse linear system was solved using the Induced dimensional reduction (IDR(s)) method with an incomplete LU preconditioner. We have applied the developed algorithm to compute a typical MCSEM response over several 3D models, we compare the numerical results from our method and FEM solution based on conventional first and second order element. The results of comparison shown that our algorithm can significantly reduce the computation cost without loss to much accuracy.

Keywords

Marine controlled sources electromagnetic method second order element edge-based finite element method anisotropic modeling

Get full access to this article

View all access options for this article.

References

Hohmann

G.W.

, Three-dimensional induced polarization and electromagnetic modeling, Geophysics 40 (1975), 309–324.

Newman

G.A.

and Hohmann

G.W.

, Transient electromagnetic response of high-contrast prisms in a layered earth, Geophysics 53 (1988), 691–706.

Avdeev

and Knizhnik

, 3D integral equation modeling with a linear dependence on dimensions, Geophysics 74 (2009), F89–F94.

Angiulli

Cacciola

Calcagno

De Carlo

Morabito

C.F.

Sgró

and Versaci

, A numerical study on the performances of the flexible BiCGStab to solve the discretized E-field integral equation, International Journal of Applied Electromagnetics and Mechanics 46(3) (2014), 547–553, doi: 10.3233/JAE-141939.

Mackie

R.M.

Smith

J.T.

and Madden

T.R.

, Three-dimensional electromagnetic modeling using finite difference equations: the magnetotelluric example, Radio Sci. 29 (1994), 923–935.

Newman

G.A.

and Alumbaugh

D.L.

, Frequency-domain modeling of airborne electromagnetic responses using staggered finite differences, Geophys Prospect 43 (1995), 1021–1042.

Haber

Ascher

U.M.

Aruliah

D.M.

and Oldenberg

, Fast simulation of 3D electromagnetic problems using potentials, J. Comput. Phys. 163 (2000), 150–171.

Jahandari

and Farquharson

C.G.

, A finite-volume solution to the geophysical electromagnetic forward problem using unstructured grids, Geophysics 79 (2014), E287–E302.

Jahandari

and Farquharson

C.G.

, Finite-volume modelling of geophysical electromagnetic data on unstructured grids using potentials, Geophys. J. Int. 202 (2015), 1859–1876.

10.

Badea

Everett

Newman

and Biro

, Finite element analysis of controlled-source electromagnetic induction using coulomb-gauged potentials, Geophysics 66 (2001), 786–799.

11.

Schwarzbach

B̈orner

R.U.

and Spitzer

, Three-dimensional adaptive higher order finite element simulation for geo electromagnetics – A marine CSEM example, Geophys. J. Int. 178 (2011), 63–74.

12.

Ansari

and Farquharson

C.G.

, 3D finite-element forward modeling of electromagnetic data using vector and scalar potentials and unstructured grids, Geophysics 79 (2014), E149–E165.

13.

Hussain

Karsiti

M.N.

Jeoti

and Yahya

, Use of wavelets in marine controlled source electromagnetic method for geophysical modeling, International Journal of Applied Electromagnetics and Mechanics (2016), 431–443, doi: 10.3233/JAE-160002.

14.

Cai

Xiong

Han

and Zhdanov

, 3D controlled-source electromagnetic modeling in anisotropic medium using edge-based finite element method, Computers & Geosciences 73 (2014), 164–176.

15.

Cai

Xiong

and Zhdanov

, Three-dimensional marine controlled-source electromagnetic modelling in anisotropic medium using finite element method, Chinese Journal of Geophysics 58 (2015), 2839–2850.

16.

Chen

Wang

and Xiong

, Characteristics analysis of 3D marine CSEM response using the vector finite element method with secondary formulation, International Journal of Applied Electromagnetics and Mechanics (2017), 565–581.

17.

Zhdanov

M.S.

, Geophysical Electromagnetic Theory and Methods, Elsevier, 2009.

18.

Zhdanov

M.S.

Lee

S.K.

and Yoshioka

, Integral equation method for 3D modeling of electromagnetic fields in complex structures with inhomogeneous background conductivity, Geophysics 71 (2006), 333–345.

19.

Streich

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

20.

Jin

, Finite Element Method in Electromagnetics, Wiley-IEEE Press, New York, 2002.

21.

Puzyrev

Koldan

Puente

Houzeaux

Vázquez

and Cela

J.M.

, A parallel finite-element method for three-dimensional controlled-source electromagnetic forward modelling, Geophys. J. Int. 193 (2013), 678–693.

22.

Mukherjee

and Everett

M.E.

, 3D controlled-source electromagnetic edge-based finite element modeling of conductive and permeable heterogeneities, Geophysics 76 (2011), F215–F226.

23.

Ansari

and Farquharson

C.G.

, 3D finite-element forward modeling of electromagnetic data using vector and scalar potentials and unstructured grids, Geophysics 79 (2014), E149–E165.

24.

Weiss

and Constable

, Mapping thin resistors and hydrocarbons with marine EM methods: Modeling and analysis in 3D, Geophysics 71 (2006), G321–G332.

25.

Börner

R.U.

Ernst

G.O.

and Spitzer

, Fast 3-D simulation of transient electromagnetic fields by model reduction in the frequency domain using Krylov subspace projection, Geophysical Journal International 173 (2008), 766–780.

26.

Fujiwara

Nakata

and Ohashi

, Improvement of convergence characteristic of ICCG method for the A-ϕ method using edge elements, IEEE Transactions on Magnetics 32 (1996), 804–807.

27.

Dyczij-Edlinger

Peng

and Lee

, A fast vector-potential method using tangentially continuous vector finite elements, IEEE Transactions on Microwave Theory and Techniques 46 (1998), 863–868.

28.

Biró

and Preis

, Partial tree-gauging of second order edge element vector potential formulations, in: Proc. Compumag Conf., Montreal, QC, Canada, Jun./Jul., 2015.

29.

Sonneveld

and van Gijzen

M.B.

, IDR(s): a family of simple and fast algorithms for solving large nonsymmetric systems of linear equations, SIAM J. Sci. Comput 31(2) (2008/09), 1035–1062, doi: 10.1137/070685804.

30.

Saad

, Iterative methods for sparse linear systems: SIAM, 2003.

31.

Ward

S.H.

and Hohmann

G.W.

, Electromagnetic Theory for Geophysical Applications, 1988, SEG, Oklahoma.

32.

Anderson

W.L.

, A hybrid fast Hankel transform algorithm for electromagnetic modeling, Geophysics 54 (1989), 263–266.