3D Wavelet Finite-Element Modeling of Frequency-Domain Airborne EM Data Based on B-Spline Wavelet on the Interval Using Potentials

We present a wavelet finite-element method (WFEM) based on B-spline wavelets on the interval (BSWI) for three-dimensional (3D) frequency-domain airborne EM modeling using a secondary coupled-potential formulation. The BSWI, which is constructed on the interval (0, 1) by joining piecewise B-spline polynomials between nodes together, has proved to have excellent numerical properties of multiresolution and sparsity and thus is utilized as the basis function in our WFEM. Compared to conventional basis functions, the BSWI is able to provide higher interpolating accuracy and boundary stability. Furthermore, due to the sparsity of the wavelet, the coefficient matrix obtained by BSWI-based WFEM is sparser than that formed by general FEM, which can lead to shorter solution time for the linear equations system. To verify the accuracy and efficiency of our proposed method, we ran numerical experiments on a half-space model and a layered model and compared the results with one-dimensional (1D) semi-analytic solutions and those obtained from conventional FEM. We then studied a synthetic 3D model using different meshes and BSWI basis at different scales. The results show that our method depends less on the mesh, and the accuracy can be improved by both mesh refinement and scale enhancement.

[1]  W. Sweldens.  The Lifting Scheme: A Custom - Design Construction of Biorthogonal Wavelets "Industrial Mathematics , 1996 .

[2]  Toshihiro Uchida,et al.  3D magnetotelluric modeling using the T‐Ω finite‐element method , 2004 .

[3]  C. Chui,et al.  On solving first-kind integral equations using wavelets on a bounded interval , 1995 .

[4]  T. Belytschko,et al.  Finite element derivative recovery by moving least square interpolants , 1994 .

[5]  I. J. Schoenberg.  Contributions to the problem of approximation of equidistant data by analytic functions. Part B. On the problem of osculatory interpolation. A second class of analytic approximation formulae , 1946 .

[6]  José María Cela,et al.  Parallel 3-D marine controlled-source electromagnetic modelling using high-order tetrahedral Nédélec elements , 2019, Geophysical Journal International.

[7]  A. Barbosa,et al.  Implementation of an hybrid-mixed stress model based on the use of wavelets , 2006 .

[8]  Zhengyong Ren,et al.  An investigation into the electromagnetic response of porous media with GPR using stochastic processes and FEM of B-spline wavelet on the interval , 2019, Journal of Applied Geophysics.

[9]  L. Castro.  Polynomial wavelets in hybrid‐mixed stress finite element models , 2009 .

[10]  K. Yee.  Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media , 1966 .

[11]  Hanbo Chen,et al.  3-D marine controlled-source electromagnetic modelling in an anisotropic medium using a Wavelet–Galerkin method with a secondary potential formulation , 2019, Geophysical Journal International.

[12]  Bo Zhang,et al.  A goal-oriented adaptive finite-element method for 3D scattered airborne electromagnetic method modeling , 2016 .

[13]  Y. Meyer.  Principe d'incertitude, bases hilbertiennes et algèbres d'opérateurs , 1986 .

[14]  I. Daubechies.  Orthonormal bases of compactly supported wavelets , 1988 .

[15]  J. Morlet,et al.  Wave propagation and sampling theory—Part II: Sampling theory and complex waves , 1982 .

[16]  Bernhard Siemon,et al.  Airborne Electromagnetic, Magnetic, and Radiometric Surveys at the German North Sea Coast Applied to Groundwater and Soil Investigations , 2020, Remote. Sens..

[17]  N. Christensen.  Interpretation attributes derived from airborne electromagnetic inversion models using the continuous wavelet transform , 2018, Near Surface Geophysics.

[18]  C. Chui,et al.  Wavelets on a Bounded Interval , 1992 .

[19]  Jean-Yves L'Excellent,et al.  On Computing Inverse Entries of a Sparse Matrix in an Out-of-Core Environment , 2012, SIAM J. Sci. Comput..

[20]  M. L. Barton,et al.  New vector finite elements for three‐dimensional magnetic field computation , 1987 .

[21]  Xinming Zhang,et al.  A Wavelet Finite Element Method for 2‐D Wave Equation in the Fluid‐Saturated Porous Media , 2005 .

[22]  B. Zhang,et al.  3D frequency-domain airborne EM forward modelling using spectral element method with Gauss–Lobatto–Chebyshev polynomials , 2019, Exploration Geophysics.

[23]  Chester J. Weiss,et al.  Project APhiD: A Lorenz-gauged A-Φ decomposition for parallelized computation of ultra-broadband electromagnetic induction in a fully heterogeneous Earth , 2013, Comput. Geosci..

[24]  Guillaume Houzeaux,et al.  A parallel finite-element method for three-dimensional controlled-source electromagnetic forward modelling , 2013 .

[25]  D. C. Fraser,et al.  The effect of the electrical anisotropy on the response of helicopter‐borne frequency‐domain electromagnetic systems , 2004 .

[26]  S. Mallat.  Multiresolution approximations and wavelet orthonormal bases of L^2(R) , 1989 .

[27]  Zhengjia He,et al.  The construction of plane elastomechanics and Mindlin plate elements of B-spline wavelet on the interval , 2006 .

[28]  G. Newman,et al.  Frequency‐domain modelling of airborne electromagnetic responses using staggered finite differences , 1995 .

[29]  Kevin Amaratunga,et al.  WAVELET BASED GREEN'S FUNCTION APPROACH TO 2D PDEs , 1993 .

[30]  Changchun Yin,et al.  Wavelet-based 3-D inversion for frequency-domain airborne EM data , 2018 .

[31]  K. Preis,et al.  On the use of the magnetic vector potential in the finite-element analysis of three-dimensional eddy currents , 1989 .

[32]  José Márcio Machado,et al.  Wavelet-Galerkin method for computations of electromagnetic fields-computation of connection coefficients , 2000 .

[33]  Changchun Yin,et al.  3D anisotropic modeling for airborne EM systems using finite-difference method , 2014 .

[34]  Oszkar Biro,et al.  Finite-element analysis of controlled-source electromagnetic induction using Coulomb-gauged potentials , 2001 .

[35]  A. Grossmann,et al.  DECOMPOSITION OF HARDY FUNCTIONS INTO SQUARE INTEGRABLE WAVELETS OF CONSTANT SHAPE , 1984 .

[36]  Zhengjia He,et al.  Adaptive multiresolution finite element method based on second generation wavelets , 2007 .

[37]  I. Johnstone,et al.  Ideal spatial adaptation by wavelet shrinkage , 1994 .

[38]  D. Okaya,et al.  Frequency‐time decomposition of seismic data using wavelet‐based methods , 1995 .