Abstract:SUMMARY For forward modelling of realistic 3-D land-based controlled-source electromagnetic (EM) problems, we develop a parallel spectral element approach, blending the flexibility and versatility of the finite element method in using unstructured grids with the accuracy of the spectral method. Complex-shaped structures and topography are accommodated by using unstructured hexahedral meshes, in which the elements can have curved edges and non-planar faces. Our code is the first spectral element algorithm in EM geophysics that uses the total field formulation (here that of the electric field). Combining unstructured grids and a total field formulation provides advantages in dealing with topography, in particular, when the transmitter is located on rough surface topography. As a further improvement over existing spectral element methods, our approach does not only allow for arbitrary distributions of conductivity, but also of magnetic permeability and dielectric permittivity. The total electric field on the elements is expanded in terms of high-order Lagrangian interpolants, and element-wise integration in the weak form of the boundary value problem is accomplished by Gauss–Legendre–Lobatto quadrature. The resulting complex-valued linear system of equations is solved using the direct solver MUMPS, and, subsequently, the magnetic field is computed at the points of interest by Faraday’s law. Five numerical examples comprehensively study the benefits of this algorithm. Comparisons to semi-analytical and finite element results confirm accurate representation of the EM responses and indicate low dependency on mesh discretization for the spectral element method. A convergence study illuminates the relation between high order polynomial approximation and mesh size and their effects on accuracy and computational cost revealing that high-order approximation yields accurate modelling results for very coarse meshes but is accompanied by high computational cost. The presented numerical experiments give evidence that 2nd and 3rd degree polynomials in combination with moderately discretized meshes provide better trade-offs in terms of computational resources and accuracy than lowest and higher order spectral element methods. To our knowledge, our final example that includes pronounced surface topography and two geometrically complicated conductive anomalies represents the first successful attempt at using 2nd order hexahedral elements supporting curved edges and non-planar faces in controlled-source EM geophysics.

Forward calculation of 3D CSEM responses based on joint application of secondary field and coupled potential formulations

Forward Calculation of 3D Controlled-Source Electromagnetic Responses Based on Joint Application of Secondary Field and Coupled Potential Formulations

Joint Application of Secondary Field and Coupled Potential Formulations to Unstructured Meshes for 3-D CSEM Forward Modeling.

Forward modeling method of 2D transient electromagnetic field induced by line source

3D frequency-domain modeling of marine controlled source electromagnetic responses with topography using finite volume method

A general forward solver for 3D CSEMs with multitype sources and operating environments

3-D Marine CSEM Forward Modeling in Anisotropic Media Using a Space-Wavenumber Domain Integral Equation Method

Three-Dimensional Inversion of CSEM Data Using Finite Element Method in Data Space

Solution of large-scale 3D controlled-source electromagnetic modeling problem using efficient iterative solvers

Fast 3-D Controlled-Source Electromagnetic Modeling Combining UPML and Rational Krylov Method

3D Controlled‐source electromagnetic modelling in anisotropic media using secondary potentials and a cascadic multigrid solver

Finite-Element Scheme With Total-Field for 3-D Marine Controlled-Source Electromagnetic Modeling With Bathymetry and Multisources

High-order Finite Element Forward Modeling Algorithm for Three-dimensional Borehole-to-surface Electromagnetic Method using Octree Meshes

Spectral element method for 3-D controlled-source electromagnetic forward modelling using unstructured hexahedral meshes

A Coupled FEM-BEM Methodology to Calculate 3-D Earth Surface Potentials Caused by HVDC Return Currents

An interpolating scaling functions method with low-storage five-stage fourth-order explicit Runge-Kutta schemes for 3D ground penetrating radar simulation

3D finite-difference transient electromagnetic modeling with a whole-space initial field

3D Airborne EM Forward Modeling Based on Time-Domain Spectral Element Method

Goal-Oriented 3-D Time-Domain Marine CSEM Modeling With Anisotropy and Topography

3D Finite-Element Modeling of Earth Induced Electromagnetic Field and Its Potential Applications for Geomagnetic Satellites

A Comparison of Different Divergence-free Solutions for 3D Anisotropic CSEM Modeling Using Staggered Finite Difference Method