Journal article
A divergence-free vector finite-element method for efficient 3D magnetotelluric forward modeling
Geophysics, Vol.89(1), pp.E1-E11
01/01/2024
DOI: 10.1190/geo2023-0037.1
Abstract
For large-scale magnetotelluric (MT) forward modeling using vector finite-element methods, iterative solvers are commonly applied. However, as frequency decreases close to zero, the iterative solution process struggles to converge. This is mainly caused by the fact that the weak conductivity term in the curl-curl equation governing the electromagnetic (EM) diffusion in the earth fades during the iterative solution for EM fields. This can lead to violation of the divergence-free condition for current and false jumps in the calculated fields. We develop a regularization technique for vector finite-element MT forward modeling, in which a scaled grad-div operator for electrical fields is included in the curl-curl equation to enforce the divergence-free condition explicitly. Because of its use of basis functions, the direct edge element discretization of the scaled term can lead to zero coefficients. To address this, an equivalent substitute of the scaled grad-div operator based on potential representation of electrical fields is used, discretized with node elements. For one specific element, this is finally reduced to a scaled gradient of the surface integral of current across all the surfaces of the element, resulting in better connectivity of the linear system of equations. The correctness of our algorithm is verified with two synthetic models and an inversion model from real data. The numerical performance is compared to the traditional iterative divergence correction technique (and without the application of the divergence correction for one case) for each of the three models. The results indicate that our algorithm is generally more efficient and stable compared to the traditional technique for all models at all periods considered, with significant improvement at long periods.
Details
- Title: Subtitle
- A divergence-free vector finite-element method for efficient 3D magnetotelluric forward modeling
- Creators
- Yongfei Wang - Central South UniversityRongwen Guo - Central South UniversityJianxin Liu - Central South UniversityJian Li - Central South UniversityRong Liu - Central South UniversityHang Chen - Boise State UniversityXun Cao - Central South UniversityZihui Yin - Central South UniversityChuanghua Cao - Shanghai Institute of Geological Survey
- Resource Type
- Journal article
- Publication Details
- Geophysics, Vol.89(1), pp.E1-E11
- DOI
- 10.1190/geo2023-0037.1
- ISSN
- 0016-8033
- eISSN
- 1942-2156
- Grant note
- 2023JJ60164 / Natural Science Foundation of Hunan Province (http://data.elsevier.com/vocabulary/SciValFunders/501100004735) HNGSTP202201 / Institute of Geological Survey of Hunan Province 42074165; 42130810; 42174171; 42374093 / National Natural Science Foundation of China (http://data.elsevier.com/vocabulary/SciValFunders/501100001809)
- Language
- English
- Date published
- 01/01/2024
- Academic Unit
- Earth and Environmental Sciences
- Record Identifier
- 9984962630202771
Metrics
1 Record Views