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.

A Model-Based Inversion Method for Reconstruction of 3d Anisotropic Permittivities

RECONSTRUCTION OF 3D ANISOTROPIC OBJECTS BY VIE AND MODEL-BASED INVERSION METHODS

Model-based Inversion of 3D Anisotropic Objects and Sensitivity Analysis

Retrieving Material Permittivity Under Oblique Incidence

Far-field Characterization of Material Permittivity Based on Nonlinear Inversions

Three-Dimensional Anisotropic Inversion for Direct Current Resistivity Method Using Unstructured Tetrahedral Discretization

Inversion of CSEM Data with Model-based Inversion Algorithm

Multiparametric Electromagnetic Inversion of 3-D Biaxial Anisotropic Objects Embedded in Layered Uniaxial Media Using VBIM Enhanced by Structural Consistency Constraint

FORWARD MODELLING AND INVERSION OF 3-D CROSS-HOLE ELECTROMAGNETIC FIELDS

3D transient electromagnetic inversion based on explicit finite-difference forward modeling

Inversion of Electromagnetic Scattering for 3D Dielectric Objects Through Integral Equation Method with Nyström Discretization

A Three-Dimensional Model-Based Inversion Algorithm Using Radial Basis Functions for Microwave Data

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

3-D 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

3-D Direct Current Resistivity Anisotropic Modelling by Goal-Oriented Adaptive Finite Element Methods

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

Reconstruction of the 2-Dimensional Inhomogeneous Conductivity Profile Using Iterative and Bi-Conjugate Gradient Method

Mathematical and Numerical Study of a Three-Dimensional Inverse Eddy Current Problem

Robust Iterative Permittivity Inversion Method Based on Reduced Order Models

Two-dimensional Inversion of Spectral Induced Polarization Data Using MPI Parallel Algorithm in Data Space