Aarhus University Seal / Aarhus Universitets segl

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

Publikation: Bidrag til tidsskrift/Konferencebidrag i tidsskrift /Bidrag til avisTidsskriftartikelForskningpeer review

Standard

3D inversion of time domain electromagnetic data using finite elements and a triple mesh formulation. / Zhang, Bo; Engebretsen, Kim Wann; Fiandaca, Gianluca; Cai, Hongzhu; Auken, Esben.

I: Geophysics, Bind 86, Nr. 3, 05.2021, s. E257-E267.

Publikation: Bidrag til tidsskrift/Konferencebidrag i tidsskrift /Bidrag til avisTidsskriftartikelForskningpeer review

Harvard

APA

CBE

MLA

Vancouver

Author

Bibtex

@article{227a1f51057d4654a9b567a12e93e959,
title = "3D inversion of time domain electromagnetic data using finite elements and a triple mesh formulation",
abstract = "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.",
keywords = "3-D INVERSION, 3-DIMENSIONAL MAGNETOTELLURIC INVERSION, SOLVER, TOPOGRAPHY",
author = "Bo Zhang and Engebretsen, {Kim Wann} and Gianluca Fiandaca and Hongzhu Cai and Esben Auken",
year = "2021",
month = may,
doi = "10.1190/geo2020-0079.1",
language = "English",
volume = "86",
pages = "E257--E267",
journal = "Geophysics",
issn = "0016-8033",
publisher = "Society of Exploration Geophysicists",
number = "3",

}

RIS

TY - JOUR

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

AU - Zhang, Bo

AU - Engebretsen, Kim Wann

AU - Fiandaca, Gianluca

AU - Cai, Hongzhu

AU - Auken, Esben

PY - 2021/5

Y1 - 2021/5

N2 - 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.

AB - 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.

KW - 3-D INVERSION

KW - 3-DIMENSIONAL MAGNETOTELLURIC INVERSION

KW - SOLVER

KW - TOPOGRAPHY

U2 - 10.1190/geo2020-0079.1

DO - 10.1190/geo2020-0079.1

M3 - Journal article

VL - 86

SP - E257-E267

JO - Geophysics

JF - Geophysics

SN - 0016-8033

IS - 3

ER -