Abstract:Summary Geophysicists today face the challenge of quickly and reliably interpreting extensive controlled-source electromagnetic (CSEM) datasets to map subsurface conductivity structures within realistic geological environments. An ideal 3D CSEM inversion algorithm using tetrahedral grids should be capable of distinguishing different resolution requirements between forward modeling and inversion grids, have an optimal parallel strategy that fully exploits the inherent independence of CSEM datasets while also possessing the capability to handle large-scale geo-electrical models, and incorporate conductivity anisotropy which should be a common characteristic in realistic subsurface environments. However, existing tools in the geo-electromagnetic community often fall short of these three demands. Addressing this gap, our study introduces a scalable and parallel anisotropic inversion technique for CSEM data, capitalizing on the potential of unstructured tetrahedral grids. We first apply the tetrahedral longest-edge bisection method to create a refined dense, heterogeneous forward modeling grid from a coarse inversion grid. This refinement, focused on areas around transmitters and receivers, is seamlessly integrated within the coarser inversion grid’s topology, enabling precise conductivity mapping and preserving electromagnetic response accuracy during model updates. We further innovate with a source-mesh double-level parallel strategy, utilizing the message passing interface technique for parallel handling of independent CSEM datasets and large-scale geo-electrical models. Externally, we dedicate a processor for inversion model updates employing the Limited-memory Broyden-Fletcher-Goldfarb-Shanno optimization algorithm and divide other processors into groups, each associated with specific transmitting sources and frequencies. Internally, in each group, we employ a domain-decomposition based scalable and robust iterative solvers using the Auxiliary-Space Maxwell preconditioner to parallel quickly calculate the electromagnetic responses from its assigned source-frequency set. Additionally, recognizing the potential for electrical conductivity anisotropy in field data, we incorporate the case of vertical transverse isotropy. We validate the effectiveness of our method through examples, including an isotropic land model with undulating topography, an anisotropic marine model, and a real-field data case. Results from both synthetic and field data inversions underscore our method’s significant advancements in efficiency and practicality, particularly in addressing large-scale 3D CSEM datasets inversion challenges in realistic geological environments.

Quasi-3D modeling of airborne TEM data by spatially constrained inversion

3D Inversion of Large-Scale Frequency-Domain Airborne Electromagnetic Data Using Unstructured Local Mesh

Quasi-2D inversion of surface large fixed-loop transient electromagnetic sounding data

Inversion for 3D Conductivity and Chargeability Models Using EM Data Acquired by the New Airborne TargetEM System in Ontario, Canada

Three-Dimensional Separate and Joint Inversions of Multi-Component Frequency-Domain Airborne Electromagnetic Data: Synthetic Model Studies

3D Time-Domain Airborne Electromagnetic Inversion Based on Secondary Field Finite-Volume Method

Three-dimensional Inversion of Semi-Airborne Transient Electromagnetic Data Based on Finite Element Method

Quasi-2-D Bayesian Inversion of Central Loop Transient Electromagnetic Data Using an Adaptive Voronoi Parametrization

Sparse-Promoting 3-D Airborne Electromagnetic Inversion Based on Shearlet Transform

Fast 2.5D and 3D Inversion of Transient Electromagnetic Surveys Using the Octree-Based Finite-Element Method

3-D Inversion of Semi-Airborne Transient Electromagnetic Data Based on Decoupled Mesh

3-D Structurally Constrained Inversion of the Controlled-Source Electromagnetic Data Using Octree Meshes

3D Structurally Constrained Inversion of the Controlled-Source Electromagnetic Data Using Octree Meshes

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

3D Inversion of Frequency-Domain Marine CSEM Data in VTI Media

3-D Parallel Anisotropic Inversion of Controlled-Source Electromagnetic Data Using Nested Tetrahedral Grids

Multiscale 3-D Stochastic Inversion of Frequency-Domain Airborne Electromagnetic Data

3D parallel anisotropic inversion of controlled-source electromagnetic data using nested tetrahedral grids

3D Near-to-surface Conductivity Reconstruction by Inversion of VETEM Data Using the Distorted Born Iterative Method

A fast 3-D inversion for airborne EM data using pre-conditioned stochastic gradient descent

3D near-to-surface conductivity reconstruction by inversion of VETEM data using the distorted Born iterative method : Electromagnetic characterization of buried obstacles