Three-Dimensional Anisotropic Inversion for Controlled-Source Audio-Frequency Magnetotellurics Using Unstructured Tetrahedral Discretization

We propose a new 3-D arbitrary anisotropic inversion algorithm for controlled-source audio-frequency magnetotellurics (CSAMTs) with consideration of angle parameters. In this method, the study area is discretized using unstructured tetrahedral meshes, and vector basis functions are employed to inter...

Full description

Saved in:
Bibliographic Details
Published in:IEEE transactions on geoscience and remote sensing Vol. 63; pp. 1 - 17
Main Authors: Ma, Fengqun, Tan, Handong, Zhu, Depeng, Riahi, Mohamed Kamel, Kong, Wenxin, Wang, Shuya
Format: Journal Article
Language:English
Published: New York IEEE 2025
The Institute of Electrical and Electronics Engineers, Inc. (IEEE)
Subjects:
3-D inversion
Algorithms
Anisotropic
Anisotropic magnetoresistance
Anisotropic media
Anisotropy
Basis functions
Conductivity
controlled-source audio-frequency magnetotellurics (CSAMTs)
Discretization
Electromagnetics
Euler angles
Geology
Image reconstruction
Impedance
Magnetic field
Magnetic fields
Mathematical models
Parameters
Perpendicular magnetic anisotropy
Rocks
Tensors
unstructured grids
ISSN:0196-2892, 1558-0644
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:We propose a new 3-D arbitrary anisotropic inversion algorithm for controlled-source audio-frequency magnetotellurics (CSAMTs) with consideration of angle parameters. In this method, the study area is discretized using unstructured tetrahedral meshes, and vector basis functions are employed to interpolate electric and magnetic field values at the receiver locations. The anisotropic resistivity in the governing electromagnetic equations is described using a tensor, and parameterized by three principal resistivities and three Euler angles. The inversion problem is solved by the Limited-memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) optimization scheme to efficiently find a reasonable model. We designed three representative types of synthetic models to investigate and test the algorithm: 1) anisotropic targets embedded in isotropic surrounding rocks; 2) isotropic targets surrounded by anisotropic media; and 3) complex models in which both targets and surrounding rocks exhibit anisotropy. The results demonstrate that the proposed inversion method can successfully recover these two horizontal principal resistivities, especially addressing the gap left by previous studies, as well as their corresponding rotation angles. For the second type of model, we further examined the influence of various initial models, particularly on the recovery of anisotropy angle parameters. Finally, by the separate inversion using the diagonal and off-diagonal components of the impedance tensor, we confirm the crucial role of off-diagonal components in the accurate reconstruction of the anomaly's shape and, especially, the angle parameters.
Bibliography:ObjectType-Article-1
SourceType-Scholarly Journals-1
ObjectType-Feature-2
content type line 14
ISSN:0196-2892
1558-0644
DOI:10.1109/TGRS.2025.3612987