Three-Dimensional Unstructured Finite Element Modeling of Magnetotelluric Problems Allowing for Continuous Variation of Conductivity in Each Block
Nian Yu,Hongye Zhang,Huang Chen,Wenxin Kong,Xiao Feng,Yunyi Qian
DOI: https://doi.org/10.1109/tgrs.2024.3398601
IF: 8.2
2024-05-25
IEEE Transactions on Geoscience and Remote Sensing
Abstract:At present, the 3-D magnetotelluric (MT) unstructured finite element modeling approaches often adopt the block-wise uniform conductivity parameterization method, in which the Earth conductivity model is discretized into plenty of tetrahedrons with uniform conductivity. However, that block-wise uniform parameterization method is not only difficult to accurately simulate actual continuous subsurface conductivity distribution, but also would substantially increase the unknowns in the subsequent inversion, which largely aggravates the non-uniqueness of inversion. In this study, we propose a 3-D unstructured finite element modeling strategy of MT problems allowing for continuous variation of conductivity in each tetrahedron. First, the curl–curl governing differential equation of electrical field for the 3-D MT problems and the corresponding integral weak form are derived successively. Next, the open-source tetrahedral mesh generator Tetgen is utilized to discretize the computational domain, both the conductivity and the electric field within each tetrahedral element are approximated by using linear interpolation with node and edge vector basis functions, respectively. Then, the closed-form expression of the resulting elemental integral is derived and the Newman-type boundary condition is loaded. The resulting system of linear equations is solved by using the direct LU decomposition method. Further, to enhance the accuracy of the solution, a goal-oriented adaptive technology based on the continuity of the normal component of the current density is adopted. Finally, the correctness of the proposed algorithm is verified by using a layered model with continuous conductivity distribution and a 3-D topographic model. While the necessity of considering the continuous variation of conductivity in blocks is proved by using two 3-D continuum models. The results show that the proposed algorithm can obviously improve the accuracy of the responses at the sites above the region with continuous conductivity distribution and will greatly reduce the unknowns of the subsequent inversion, compared with the traditional algorithm with block-wise uniform parameterization.
imaging science & photographic technology,remote sensing,engineering, electrical & electronic,geochemistry & geophysics