Abstract:The Controlled-source electromagnetic (CSEM) method in the frequency domain has evolved into an established technique for mineral resources prospecting and hydrocarbon exploration. An increasing trend is to carry out three-dimensional (3D) CSEM surveys as EM explorations now are increasingly conducted in complex geological environments in order to improve the spatial resolution of subsurface conductivity structure.Quantitative interpretation of large-scale CSEM data in the frequency domain requires efficient and stable 3D forward modeling and inversion codes. Considerable efforts have been contributed to developing numerical algorithms concerning 3D inversion of CSEM data that can accurately and efficiently recover subsurface electrical structure over last decade. Most existing 3D inversion algorithms employ Krylov subspace iterative methods as their forward solvers. Iterative techniques usually need little memory due to only matrix-vector products storage required, and they are fast for computation of single field solution. However, there are two main issues when using iterative solvers for 3D CSEM problems: (1) The ill-conditioning of linear systems arising from discretization of Helmholtz equation can lead to poor behavior of iterative solvers and even divergence in some cases. (2) Iterative solvers are very time-consuming for multi-source problems. These difficulties become major impediments when solving large-scale multi-source CSEM problems using iterative solvers.Given the availability of more powerful workstations or computer clusters, direct methods have been increasingly used for solving 3D CSEM problems. Compared to iterative methods, direct solvers have several distinct advantages: (1) They provide more accurate solutions. (2) They are less prone to ill-conditioning of matrix, making them more robust for highly heterogeneous models or non-uniform grids. And (3) they separate the solving of linear system into an expensive matrix factorization and comparably inexpensive forward-backward substitution steps, which make them more suitable for multi-source CSEM surveying.In this paper, we present an efficient inversion algorithm for 3D inversion of multi-frequency and multi-source CSEM data, which is based on a direct solver for solving the forward problem. The Gauss-Newton (GN) optimization algorithm is applied considering its high convergence rate, thus limiting the number of expensive matrix factorization required when using a direct solver. A preconditioned conjugate gradient solver (PCG) is used to solve the normal equation resulted from linearization at each GN iterate, in order to avoid computing and storing sensitivity matrix explicitly. This scheme only requires matrix-vector products of Jocabian and its transpose with vectors, which are equivalent to one forward and one adjoint problem. In addition, the matrix factorization obtained when solving forward problem can be used in subsequent PCG iterates, which dramatically speeds up PCG iteration and reduces computational costs. Numerical experiments on synthetic data from land and marine CSEM surveys show that our inversion algorithm has an excellent convergence rate and only tens to several tens of iterates are needed to reach desired misfit, demonstrating its efficiency and stability.

An Alternating Direction Multiplier Method for the Inversion of FDEM Data

An Alternating Direction Method of Multipliers Algorithm for 1-D Magnetotelluric Anisotropic Inversion Using Fourier Series Expansion

Artificial Bee Colony Algorithm with Adaptive Parameter Space Dimension: A Promising Tool for Geophysical Electromagnetic Induction Inversion

Rapid Surrogate Modeling of Electromagnetic Data in Frequency Domain Using Neural Operator.

FDEMtools: a MATLAB package for FDEM data inversion

Solving Electromagnetic Inverse Problem Using Adaptive Gradient Descent Algorithm.

Robust elastic full waveform inversion using alternating direction method of multipliers with reconstructed wavefields

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

Effective 3D-Transient Electromagnetic Inversion Using Finite-Element Method with a Parallel Direct Solver

3d Inversion Of Frequency-Domain Csem Data Based On Gauss-Newton Optimization

Parallelized 3-D Csem Inversion With Secondary Field Formulation And Hexahedral Mesh

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

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

Probabilistic One-Dimensional Inversion of Frequency-Domain Electromagnetic Data Using a Kalman Ensemble Generator

Robust 3D Joint Inversion of Gravity and Magnetic Data: A High-Performance Computing Approach

Iterative ADMM for Inverse FE–BI Problem: A Potential Solution to Radio Tomography of Asteroids

Parallel Inversion of 3D Airborne Transient Electromagnetic Data Using an Approximate Jacobi Matrix

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

Efficient Solution Scheme for Large-Scale Anisotropic Forward Modeling of 3-D Magnetotelluric Data

Application of Supervised Descent Method to Nonlinear Inversion for Transient EM Data

2D Magnetotelluric Inversion Using Linear Finite Element Methods and a Discretize‐Last Strategy With First and Second–Order Anisotropic Regularization