Spectral-Element Method With Divergence-Free Constraint for 2.5-D Marine CSEM Hydrocarbon Exploration

Rapid simulations of large-scale low-frequency subsurface electromagnetic measurements are still a challenge because of the low-frequency breakdown phenomenon that makes the system matrix extremely poor-conditioned. Hence, significant attention has been paid to accelerate the numerical algorithms for Maxwell’s equations in both integral and partial differential forms. In this letter, we develop a novel 2.5-D method to overcome the low-frequency breakdown problem by using the mixed spectral element method with the divergence-free constraint and apply it to solve the marine-controlled-source electromagnetic systems. By imposing the divergence-free constraint, the proposed method considers the law of conservation of charges, unlike the conventional governing equation for these problems. Therefore, at low frequencies, the Gauss law guarantees the stability of the solution, and we can obtain a well-conditioned system matrix even as the frequency approaches zero. Several numerical experiments show that the proposed method is well suited for solving low-frequency electromagnetic problems.

[1]  Aria Abubakar,et al.  2.5D forward and inverse modeling for interpreting low-frequency electromagnetic measurements , 2008 .

[2]  Qing Huo Liu,et al.  The 3-D multidomain pseudospectral time-domain algorithm for inhomogeneous conductive media , 2003, IEEE Antennas and Propagation Society International Symposium. Digest. Held in conjunction with: USNC/CNC/URSI North American Radio Sci. Meeting (Cat. No.03CH37450).

[3]  Ki Ha Lee,et al.  A numerical solution for the electromagnetic scattering by a two-dimensional inhomogeneity , 1985 .

[4]  Qing Huo Liu,et al.  An efficient 3-D spectral-element method for Schrödinger equation in nanodevice simulation , 2005, IEEE Trans. Comput. Aided Des. Integr. Circuits Syst..

[5]  Dan Jiao,et al.  A Rigorous Solution to the Low-Frequency Breakdown in Full-Wave Finite-Element-Based Analysis of General Problems Involving Inhomogeneous Lossless/Lossy Dielectrics and Nonideal Conductors , 2011, IEEE Transactions on Microwave Theory and Techniques.

[6]  Qing Huo Liu,et al.  Spectral Element Method and Domain Decomposition for Low-Frequency Subsurface EM Simulation , 2016, IEEE Geoscience and Remote Sensing Letters.

[7]  R. N. Edwards,et al.  Transient marine electromagnetics: the 2.5-D forward problem , 1993 .

[8]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[9]  Rita Streich,et al.  2 . 5 D controlled-source EM modeling with general 3 D source geometries , 2011 .

[10]  Fernando L. Teixeira,et al.  General closed-form PML constitutive tensors to match arbitrary bianisotropic and dispersive linear media , 1998 .

[11]  Qing Huo Liu,et al.  The 2.5-D multidomain pseudospectral time-domain algorithm , 2003 .

[12]  Michael S. Zhdanov,et al.  Optimal Synthetic Aperture Method for Marine Controlled-Source EM Surveys , 2015, IEEE Geoscience and Remote Sensing Letters.

[13]  Trond Mannseth,et al.  An Approximate Hybrid Method for Electromagnetic Scattering From an Underground Target , 2013, IEEE Transactions on Geoscience and Remote Sensing.

[14]  Alan D. Chave,et al.  Electromagnetic induction by a finite electric dipole source over a 2-D earth , 1993 .

[15]  Allan McKay,et al.  Anisotropic 2.5D Inversion of Towed Streamer EM Data from Three North Sea Fields Using Parallel Adaptive Finite Elements , 2014 .

[16]  R. Greenfield,et al.  Numerical solutions of the response of a two-dimensional earth to an oscillating magnetic dipole source , 1976 .

[17]  Rita Streich,et al.  2.5D controlled-source EM modeling with general 3D source geometries , 2011 .

[18]  General PML constitutive tensors to match arbitrary bianisotropic and dispersive linear media , 1998, IEEE Antennas and Propagation Society International Symposium. 1998 Digest. Antennas: Gateways to the Global Network. Held in conjunction with: USNC/URSI National Radio Science Meeting (Cat. No.98CH36.

[19]  J. Coggon Electromagnetic and electrical modeling by the finite element method , 1971 .

[20]  Tage Røsten,et al.  A 2.5D finite-element-modeling difference method for marine CSEM modeling in stratified anisotropic media , 2008 .

[21]  Q.H. Liu,et al.  A 3-D spectral-element method using mixed-order curl conforming vector basis functions for electromagnetic fields , 2006, IEEE Transactions on Microwave Theory and Techniques.

[22]  Jeffrey S. Ovall,et al.  A parallel goal-oriented adaptive finite element method for 2.5-D electromagnetic modelling , 2011 .

[23]  Siyuan Chen,et al.  Analyzing low-frequency electromagnetic scattering from a composite object , 2002, IEEE Trans. Geosci. Remote. Sens..