3D inversion of time domain electromagnetic data using finite elements and a triple mesh formulation

Bo Zhang*, Kim Wann Engebretsen, Gianluca Fiandaca, Hongzhu Cai, Esben Auken

*Corresponding author for this work

Research output: Contribution to journal/Conference contribution in journal/Contribution to newspaperJournal articleResearchpeer-review


Over several decades, much research has been done to develop 3D electromagnetic inversion algorithms. Due to the computational complexity and the memory requirements for 3D time-domain electromagnetic (TEM) inversion algorithms, many real-world surveys are inverted within one dimension. To speed up calculations and manage memory for 3D inversions of TEM data, we have developed an approach using three uncoupled meshes: an inversion mesh, a forward-model mesh, and a mesh for Jacobian calculations. The inversion mesh is a coarse regular and structured mesh, such that constraints are easily enforced between the model parameters. Forward responses are calculated on a dense unstructured mesh to obtain accurate electromagnetic fields, whereas the Jacobian is calculated on a coarse unstructured mesh. We found that using a coarse mesh for the Jacobian is sufficient for the inversion to converge and, equally important, that it provides a significant speed boost in the overall inversion process, compared to calculating it on the forward-modeling mesh. The unstructured meshes are made of tetrahedral elements, and the electromagnetic fields are calculated using the finiteelement method. The inversion optimization uses a standard Gauss-Newton formulation. For further speed up and memory optimizing of the inversion, we use domain decomposition for calculating the responses for each transmitter separately and parallelize the problem over domains using OpenMP. Compared to a 1D solution, the accuracy for the Jacobian is 1%-5% for the dense mesh and 2%-7% for the coarse mesh, but the calculation time is approximately five times faster for the coarse mesh. We also examined the algorithm on a small ground-based TEM data set acquired in an area where a 3D earth distorts the electromagnetic fields to such a degree that a 1D inversion is not feasible.

Original languageEnglish
Pages (from-to)E257-E267
Number of pages11
Publication statusPublished - May 2021




Dive into the research topics of '3D inversion of time domain electromagnetic data using finite elements and a triple mesh formulation'. Together they form a unique fingerprint.

Cite this