1887
Volume 50, Issue 5
  • ISSN: 0812-3985
  • E-ISSN: 1834-7533

Abstract

ABSTRACT

The spectral element method (SEM) based on high-order complete orthogonal polynomials is an accurate and efficient numerical method for electromagnetic modelling due to its spectral accuracy and exponential convergence. The SEM combines the flexibility of the finite-element method and the high accuracy of the spectral method. In this paper, we introduce SEM into three-dimensional frequency-domain airborne electromagnetic forward modelling. Starting from Maxwell's equations, we obtain a vector Helmholtz equation for the electric field. We use the Galerkin method to discretise the Helmholtz equation, in which the curl-conforming Gauss–Lobatto–Chebyshev polynomials are used as basis functions. The GLC polynomials help to derive the analytical expressions of entries in the system matrix and thus guarantee the modelling accuracy. Finally, we use the direct solver MUMPS to solve for the electric field and calculate the magnetic field by interpolation. For numerical experiments, we first compare our results with the semi-analytical solutions of a homogeneous half-space to verify the accuracy of our algorithm. We then analyse the characteristics of SEM by assuming different orders of interpolation polynomials and meshes. We also compare our method with the finite element method and SEM based on Gauss–Lobatto–Legendre polynomials. The results show that SEM is an efficient and effective method for electromagnetic modelling, it can deliver very accurate results and is less sensitive to mesh quality than the finite element method.

© 2019 Australian Society of Exploration Geophysics
Loading

Article metrics loading...

/content/journals/10.1080/08123985.2019.1614162
2019-09-03
2026-01-23

  • From This Site

    /content/journals/10.1080/08123985.2019.1614162
    dcterms_title,dcterms_subject,pub_keyword
    -contentType:Journal -contentType:Contributor -contentType:Concept -contentType:Institution
    10
    5
Loading full text...

Full text loading...

References

  1. Amestoy, P. R., I. S. Duff, J. Y. L’Excellent, and J. Koster 2001 A fully asynchronous multifrontal solver using distributed dynamic scheduling. SIAM Journal on Matrix Analysis and Applications23: 15–41. doi:10.1137/S0895479899358194.
     https://doi.org/10.1137/S0895479899358194 [Google Scholar]
  2. Amestoy, P. R., I. S. Duff, J. Y. L’Excellent, Y. Robert, F. H. Rouet, and B. Uçar 2012 On computing inverse entries of a sparse matrix in an out-of-core environment. SIAM Journal on Scientific Computing34: A1975–A1999. doi:10.1137/100799411.
     https://doi.org/10.1137/100799411 [Google Scholar]
  3. Avdeev, D. B. 2005 Three-dimensional electromagnetic modelling and inversion from theory to application. Surveys in Geophysics26: 767–799. doi:10.1007/s10712-005-1836-2. doi: 10.1007/s10712‑005‑1836‑x
     https://doi.org/10.1007/s10712-005-1836-x [Google Scholar]
  4. Börner, R. U. 2010 Numerical modelling in geo-electromagnetics: Advances and challenges. Surveys in Geophysics31: 225–245. doi:10.1007/s10712‑009‑9087‑x.
     https://doi.org/10.1007/s10712-009-9087-x [Google Scholar]
  5. Grayver, A. V., and T. V. Kolev 2015 Large-scale 3D geoelectromagnetic modeling using parallel adaptive high-order finite element method. Geophysics80: E277–E291. doi:10.1190/GEO2015‑0013.1.
     https://doi.org/10.1190/GEO2015-0013.1 [Google Scholar]
  6. Henderson, R. D., and G. E. Karniadakis 1995 Unstructured spectral element methods for simulation of turbulent flows. Journal of Computational Physics122: 191–217. doi:10.1006/jcph.1995.1208.
     https://doi.org/10.1006/jcph.1995.1208 [Google Scholar]
  7. Hesthaven, J. S. 1998 From electrostatics to almost optimal nodal sets for polynomial interpolation in a simplex. SIAM Journal on Numerical Analysis35: 655–676. doi:10.1137/S003614299630587X.
     https://doi.org/10.1137/S003614299630587X [Google Scholar]
  8. Ilic, M. M., and B. M. Notaros 2003 Higher order hierarchical curved hexahedral vector finite elements for electromagnetic modeling. IEEE Transactions on Microwave Theory and Techniques51: 1026–1033. doi:10.1109/TMTT.2003.808680.
     https://doi.org/10.1109/TMTT.2003.808680 [Google Scholar]
  9. Komatitsch, D., and J. Tromp 1999 Introduction to the spectral element method for three-dimensional seismic wave propagation. Geophysical Journal International139: 806–822. doi:10.1046/j.1365‑246x.1999.00967.x.
     https://doi.org/10.1046/j.1365-246x.1999.00967.x [Google Scholar]
  10. Kudela, P., A. Żak, M. Krawczuk, and W. Ostachowicz 2007 Modelling of wave propagation in composite plates using the time domain spectral element method. Journal of Sound and Vibration302: 728–745. doi:10.1016/j.jsv.2006.12.016.
     https://doi.org/10.1016/j.jsv.2006.12.016 [Google Scholar]
  11. Lee, J. H., and Q. H. Liu 2007 A 3-D spectral-element time-domain method for electromagnetic simulation. IEEE Transactions on Microwave Theory and Techniques55: 983–991. doi:10.1109/TMTT.2007.895398.
     https://doi.org/10.1109/TMTT.2007.895398 [Google Scholar]
  12. Lee, J. H., T. Xiao, and Q. H. Liu 2006 A 3-D spectral-element method using mixed-order curl conforming vector basis functions for electromagnetic fields. IEEE Transactions on Microwave Theory and Techniques54: 437–444. doi:10.1109/TMTT.2005.860502. doi: 10.1109/TMTT.2006.886157
     https://doi.org/10.1109/TMTT.2006.886157 [Google Scholar]
  13. Li, Y., and X. K. Li 2016 The Chebyshev spectral element approximation with exact quadratures. Journal of Computational and Applied Mathematics296: 320–333. doi:10.1016/j.cam.2015.09.021.
     https://doi.org/10.1016/j.cam.2015.09.021 [Google Scholar]
  14. Liu, N., G. Cai, C. Zhu, Y. Tang, and Q. H. Liu 2015 The mixed spectral-element method for anisotropic, lossy, and open waveguides. IEEE Transactions on Microwave Theory and Techniques63: 3094–3102. doi:10.1109/TMTT.2015.2472416.
     https://doi.org/10.1109/TMTT.2015.2472416 [Google Scholar]
  15. Liu, Y., and C. Yin 2014 3D anisotropic modeling for airborne EM systems using finite-difference method. Journal of Applied Geophysics109: 186–194. doi:10.1016/j.jappgeo.2014.07.003.
     https://doi.org/10.1016/j.jappgeo.2014.07.003 [Google Scholar]
  16. Maday, Y., A. T. Patera, and E. M. Ronquist 1987 A well-posed optimal spectral element approximation for the Stokes problem. National Aeronautics & Space Administration Report. [Web document]. https://ntrs.nasa.gov/archive/nasa/casi.ntrs.nasa.gov/19870018523.pdf (accessed October 31, 2017).
  17. Magnoni, F., E. Casarotti, A. Michelini, A. Piersanti, D. Komatitsch, D. Peter, and J. Tromp 2014 Spectral-element simulations of seismic waves generated by the 2009 L’Aquila earthquake. Bulletin of the Seismological Society of America104: 73–94. doi:10.1785/0120130106.
     https://doi.org/10.1785/0120130106 [Google Scholar]
  18. Nabighian, M. N. 1988 Electromagnetic methods in applied geophysics. Vol. 2. Tulsa: SEG Books.
  19. Nédélec, J. C. 1980 Mixed finite elements in ℝ3. Numerische Mathematik35: 315–341. doi:10.1007/BF01396415.
     https://doi.org/10.1007/BF01396415 [Google Scholar]
  20. Newman, G. A., and D. L. Alumbaugh 1995 Frequency-domain modelling of airborne electromagnetic responses using staggered finite differences. Geophysical Prospecting43: 1021–1042. doi:10.1111/j.1365‑2478.1995.tb00294.x.
     https://doi.org/10.1111/j.1365-2478.1995.tb00294.x [Google Scholar]
  21. Patera, A. T. 1984 A spectral element method for fluid dynamics: Laminar flow in a channel expansion. Journal of Computational Physics54: 468–488. doi:10.1016/0021‑9991(84)90128‑1.
     https://doi.org/10.1016/0021-9991(84)90128-1 [Google Scholar]
  22. Peter, D., D. Komatitsch, Y. Luo, R. Martin, N. Le Goff, E. Casarotti, P. Le Loher, et al. 2011 Forward and adjoint simulations of seismic wave propagation on fully unstructured hexahedral meshes. Geophysical Journal International186: 721–739. doi:10.1111/j.1365‑246X.2011.05044.x.
     https://doi.org/10.1111/j.1365-246X.2011.05044.x [Google Scholar]
  23. Phillips, G. M. 2003 Interpolation and approximation by polynomials. New York: Springer Science & Business Media.
  24. Raiche, A., F. Sugeng, and G. Wilson 2007 Practical 3D EM inversion? The P223F software suite. 19th Geophysical Conference, ASEG, Extended Abstracts, 1–5.
  25. Ren, Q., L. E. Tobón, Q. Sun, and Q. H. Liu 2015 A new 3-D nonspurious discontinuous galerkin spectral element time-domain (DG-SETD) method for Maxwell’s equations. IEEE Transactions on Antennas and Propagation63: 2585–2594. doi:10.1109/TAP.2015.2417891.
     https://doi.org/10.1109/TAP.2015.2417891 [Google Scholar]
  26. Runge, C. 1901 Über empirische Funktionen und die Interpolation zwischen äquidistanten Ordinaten. Zeitschrift für Mathematik und Physik46: 224–243.
     [Google Scholar]
  27. Schwarzbach, C., R. U. Börner, and K. Spitzer 2011 Three-dimensional adaptive higher order finite element simulation for geo-electromagnetics—a marine CSEM example. Geophysical Journal International187: 63–74. doi:10.1111/j.1365‑246X.2011.05127.x.
     https://doi.org/10.1111/j.1365-246X.2011.05127.x [Google Scholar]
  28. Seriani, G., and E. Priolo 1991 High-order spectral element method for acoustic wave modeling. 61th Annual International Meeting, SEG, Expanded Abstracts, 1561-1564. doi:10.1190/1.1888989.
     https://doi.org/10.1190/1.1888989
  29. Sugeng, F. 1998 Modeling the 3D TDEM response using the 3D full-domain finite-element method based on the hexahedral edge-element technique. Exploration Geophysics29: 615–619. doi:10.1071/EG998615.
     https://doi.org/10.1071/EG998615 [Google Scholar]
  30. Tape, C., Q. Liu, A. Maggi, and J. Tromp 2010 Seismic tomography of the southern California crust based on spectral-element and adjoint methods. Geophysical Journal International180: 433–462. doi:10.1111/j.1365‑246X.2009.04429.x.
     https://doi.org/10.1111/j.1365-246X.2009.04429.x [Google Scholar]
  31. Yin, C. C., W. Huang, and F. Ben 2013 The full-time electromagnetic modeling for time-domain airborne electromagnetic systems. Chinese Journal of Geophysics56: 3153–3162. doi:10.6038/cjg20130928.
     https://doi.org/10.6038/cjg20130928 [Google Scholar]
  32. Yin, C., X. Huang, Y. Liu, and J. Cai 2017 3-D modeling for airborne EM using the spectral-element method. Journal of Environmental and Engineering Geophysics22 no. 1: 13–23. doi:10.2113/JEEG22.1.13.
     https://doi.org/10.2113/JEEG22.1.13 [Google Scholar]
  33. Yin, C., Y. Qi, Y. Liu, and J. Cai 2016a 3D time-domain airborne EM forward modeling with topography. Journal of Applied Geophysics134: 11–22. doi:10.1016/j.jappgeo.2016.08.002.
     https://doi.org/10.1016/j.jappgeo.2016.08.002 [Google Scholar]
  34. Yin, C., B. Zhang, Y. Liu, and J. Cai 2016b A goal-oriented adaptive finite-element method for 3D scattered airborne electromagnetic method modeling. Geophysics81: E337–E346. doi:10.1190/GEO2015‑0580.1.
     https://doi.org/10.1190/GEO2015-0580.1 [Google Scholar]
  35. Zhdanov, M. S., S. K. Lee, and K. Yoshioka 2006 Integral equation method for 3D modeling of electromagnetic fields in complex structures with inhomogeneous background conductivity. Geophysics71: G333–G345. doi:10.1190/1.2358403.
     https://doi.org/10.1190/1.2358403 [Google Scholar]
/content/journals/10.1080/08123985.2019.1614162
Loading
3D frequency-domain airborne EM forward modelling using spectral element method with Gauss–Lobatto–Chebyshev polynomials
Exploration Geophysics 50, 461 (2019); https://doi.org/10.1080/08123985.2019.1614162
/content/journals/10.1080/08123985.2019.1614162
/content/journals/10.1080/08123985.2019.1614162
Loading

Data & Media loading...

  • Article Type: Research Article
Keyword(s): 3D modelling; airborne EM; frequency-domain; Spectral element method

Most Cited This Month Most Cited RSS feed

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error