Abstract:Practical application of 3D magnetotelluric inversion requires efficient forward modeling of electromagnetic (EM) fields in the earth. To resolve realistic 3D structures, large computational domains and extremely large linear systems of equations are required. The iterative solvers, which are almost exclusively used to solve these systems, can be inefficient due to the abundant null space of the curl-curl operator. Multigrid (MG) solvers are considered a potentially efficient technique for solving such problems. However, due to the abundant null solution space and existence of the air layer, MG solvers can still converge slowly or even diverge. We have developed an efficient MG solver for finite-difference frequency-domain EM solution. In this algorithm, the excellent smoothing property of an efficient fourcolor cell-block Gauss-Seidel (GS) is exploited to remove the short-range errors effectively, and the interpolation and prolongation operators are used to handle the long-range errors. They work as a whole to speed the convergence of our algorithm remarkably. Because all of the nodes for the four-color cell block GS are grouped into four colors and the edge components attached to different nodes in each color are completely de coupled, this can be used to develop a highly vectorized or parallelized algorithm. Another important property is that our algorithm is locally current divergence free, effectively eliminating spurious solutions in the null space of the curl-curl operator. The accuracy and efficiency of the algorithm are verified by comparing the numerical solutions obtained with our MG solver to those from the biconjugate gradient stabilized solver with different preconditioners based on synthetic models and a model from 3D inversion. Comparisons, in terms of iteration number and computational time, indicate that our algorithm is extremely stable and efficient relative to the other solvers. Our MG algorithm will be suitable for massively parallel computing as well.

A Generalized Cross Validation Method For The Inverse Problem Of 3-D Maxwell'S Equation

Choosing Near-Optimal Regularization Parameter for the Inverse Problem of Electrocardiography

An adaptive algorithm for solving inverse problems of Maxwell's Equations

A Wavelet Multiscale Method for Inversion of Maxwell Equations

A Multiresolution Method for Distributed Conductivity Estimation of Maxwell Equation

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

Extension of the Regularization Technique to Controlled-Source Electromagnetic Modeling in General Anisotropic Conductivity Media

A compressed implicit Jacobian scheme for 3D electromagnetic data inversion

An Adaptive Conjugate Gradient Least-Squares Regularization (ACGLSR) Method for 3D Gravity Density Inversion

Method of fundamental solutions for a three-dimensional inverse problem in linear elasticity

An Iterative Method for the Inverse Eddy Current Problem with Total Variation Regularization

A compressed implicit Jacobian scheme for 3D regularized Gauss-Newton inversion

An Efficient Multigrid Solver Based on a Four-Color Cell-Block Gauss-Seidel Smoother for 3D Magnetotelluric Forward Modeling

Inverse Problem for the Nonlinear Convection–Diffusion Equation by Using the Multigrid Method and Constraint Data

A conservative hybrid deep learning method for Maxwell–Ampère–Nernst–Planck equations

Gauss–Newton With Preconditioned Conjugate Gradient Magnetotelluric Inversion for 3-D Axial Anisotropic Conductivities

A general numerical method for solving the three-dimensional hyperbolic heat conduction equation on unstructured grids

Iterative Methods for an Inverse Eddy Current Problem with Total Variation Regularization

An efficient cascadic multigrid solver for 3D magnetotelluric forward modeling problems using potentials

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

A Regularization Method in the Electromagnetic Inverse Scattering Problem