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.

Cube2sph : A toolkit enabling flexible and accurate continental-scale seismic wave simulations using the SPECFEM3D_Cartesian package

Oscillatory Responses of a Spatially Random Seabed in the Vicinity of Structures: Validation and Application of the Integrated CFD-SFEM

Cubic Lattice-Based Spherical Uniaxial Perfectly Matched Layer For The Fdtd Method

Numerical modeling of the 3-D seismic wavefield with the spectral element method in spherical coordinates

Complete Synthetic Seismograms Based on A Spherical Self-Gravitating Earth Model with an Atmosphere-Ocean-Mantle-Core Structure

Three-dimensional Anisotropic Seismic Wave Modelling in Spherical Coordinates by a Collocated-Grid Finite-Difference Method

Multidomain Finite-Difference and Chebyshev Pseudospectral Hybrid Method for Elastic Wave Simulation: Two-Dimensional Case

ExoCubed: A Riemann-Solver based Cubed-Sphere Dynamic Core for Planetary Atmospheres

Spherical magnetotelluric modeling based on non-uniform source

ExoCubed: A Riemann-solver-based Cubed-sphere Dynamic Core for Planetary Atmospheres

A stable discontinuous Galerkin method for linear elastodynamics in 3D geometrically complex media using physics based numerical fluxes

Amplitude-preserving P-/S-wavefield Separation with the Discontinuous Galerkin Method on Unstructured Meshes

Modelling Seismic Wave Propagation in a Two-Dimensional Cylindrical Whole-Earth Model Using the Pseudospectral Method

Error Propagation and Control in 2D and 3D Hybrid Seismic Wave Simulations for Box Tomography

Nonsplit Complex-Frequency-shifted Perfectly Matched Layer Combined with Symplectic Methods for Solving Second-Order Seismic Wave Equations - Part 2: Wavefield Simulations

A new cube movement test for verification of simulations of contact processes of blocks of different size in geological hazards

Simulation of Seismoelectric Waves Using Time-Domain Finite-Element Method in 2D PSVTM Mode

3D Seismic-Wave Modeling with a Topographic Fluid–Solid Interface at the Sea Bottom by the Curvilinear-Grid Finite-Difference Method

Two Tactics for Accurate Frequency-Domain CSEM Modelling

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

Finite Element Simulation of 3-D Marine Controlled Source Electromagnetic Fields in Anisotropic Media with Unstructured Tetrahedral Grids