Mixed Spectral-Element Method for Overcoming the Low-Frequency Breakdown Problem in Subsurface EM Exploration

One fundamental difficulty in low-frequency subsurface electromagnetic exploration is the low-frequency breakdown phenomenon in numerical computation. It makes the discretized linear system very poorly conditioned and thus difficult to solve. This issue is present in both integral equation and partial differential equation solution methods, and thus has attracted many researchers who have proposed various methods to overcome this difficulty. In this paper, we propose a new mixed spectral element method (mixed SEM) to eliminate this low-frequency breakdown problem and apply this method to solve the subsurface electromagnetic exploration problem. Since Gauss’ law is now explicitly enforced in the mixed SEM to make the system matrix well-conditioned even at extremely low frequency, we can solve the linear system from dc to high frequencies. With the proposed method, we study the surface-to-borehole electromagnetic system for hydrocarbon exploration. Numerical examples show that the mixed SEM is accurate and efficient, and has significant advantages over conventional methods.

[1]  Jin-Fa Lee,et al.  A Domain Decomposition Method for Electromagnetic Radiation and Scattering Analysis of Multi-Target Problems , 2008, IEEE Transactions on Antennas and Propagation.

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

[3]  Qing Huo Liu,et al.  Efficient Near-Field Imaging for Single-Borehole Radar With Widely Separated Transceivers , 2015, IEEE Transactions on Geoscience and Remote Sensing.

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

[5]  D. Schötzau,et al.  Preconditioners for saddle point linear systems with highly singular blocks. , 2006 .

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

[7]  Qing Huo Liu,et al.  Three-dimensional reconstruction of objects buried in layered media using Born and distorted Born iterative methods , 2004, IEEE Geosci. Remote. Sens. Lett..

[8]  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.

[9]  Qing Huo Liu,et al.  Mixed Finite Element Method for 2D Vector Maxwell's Eigenvalue Problem in Anisotropic Media , 2014 .

[10]  Siyuan Chen,et al.  Analysis of low frequency scattering from penetrable scatterers , 2001, IEEE Trans. Geosci. Remote. Sens..

[11]  W. Chew,et al.  Integral equation solution of Maxwell's equations from zero frequency to microwave frequencies , 2000 .

[12]  Zhang Yerong,et al.  Hybrid Born iterative method in low-frequency inverse scattering problem , 1998 .

[13]  Mark E. Everett,et al.  3D controlled-source electromagnetic edge-based finite element modeling of conductive and permeable heterogeneities , 2011 .

[14]  Zhanxiang He,et al.  First Borehole to Surface Electromagnetic Survey in KSA: Reservoir Mapping and Monitoring at a New Scale , 2011 .

[15]  Gene H. Golub,et al.  Numerical solution of saddle point problems , 2005, Acta Numerica.

[16]  Weng Cho Chew,et al.  Three-dimensional imaging of buried objects in very lossy earth by inversion of VETEM data , 2003, IEEE Trans. Geosci. Remote. Sens..

[17]  Michael A. Saunders,et al.  LSQR: An Algorithm for Sparse Linear Equations and Sparse Least Squares , 1982, TOMS.

[18]  Q. Liu,et al.  Mixed Spectral-Element Method for 3-D Maxwell's Eigenvalue Problem , 2015, IEEE Transactions on Microwave Theory and Techniques.

[19]  L. Demkowicz,et al.  Modeling of electromagnetic absorption/scattering problems using hp-adaptive finite elements , 1998 .

[20]  V. Simoncini,et al.  Block--diagonal and indefinite symmetric preconditioners for mixed finite element formulations , 1999 .

[21]  J. Nédélec Mixed finite elements in ℝ3 , 1980 .

[22]  Yuefeng Sun,et al.  Rock-physics-based estimation of critical-clay-volume fraction and its effect on seismic velocity and petrophysical properties , 2014 .

[23]  Qing Huo Liu,et al.  A fast solver for vertical electromagnetic profiles of surface to borehole electromagnetic method (SBEM) , 2014 .

[24]  Colin Farquharson,et al.  Three-dimensional finite-element modelling of magnetotelluric data with a divergence correction , 2011 .

[25]  Qing Huo Liu,et al.  Applications of the BCGS-FFT method to 3-D induction well logging problems , 2003, IEEE Trans. Geosci. Remote. Sens..

[26]  Jin-Fa Lee,et al.  Removal of spurious DC modes in edge element solutions for modeling three-dimensional resonators , 2006, IEEE Transactions on Microwave Theory and Techniques.

[27]  Qing Huo Liu,et al.  An efficient solution for the response of electrical well logging tools in a complex environment , 1991, IEEE Trans. Geosci. Remote. Sens..

[28]  Dan Jiao,et al.  A unified finite-element solution from zero frequency to microwave frequencies for full-wave modeling of large-scale three-dimensional on-chip interconnect structures , 2008, 2008 IEEE Antennas and Propagation Society International Symposium.

[29]  Q. Liu,et al.  Mixed Spectral-Element Method for 3-D Maxwell's Eigenvalue Problem , 2015 .

[30]  G. Palacky 3. Resistivity Characteristics of Geologic Targets , 1988 .

[31]  Qing Huo Liu,et al.  A Semianalytical Spectral Element Method for the Analysis of 3-D Layered Structures , 2011, IEEE Transactions on Microwave Theory and Techniques.

[32]  F. Kikuchi,et al.  Mixed and penalty formulations for finite element analysis of an eigenvalue problem in electromagnetism , 1987 .

[33]  J. T. Smith Conservative modeling of 3-D electromagnetic fields, Part II: Biconjugate gradient solution and an accelerator , 1996 .