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.

Magnetotelluric Responses of a Horizontally Nonhomogeneous Model of the Earth with Anisotropic Conductivity Structures

Two-dimensional Magnetotelluric Modelling by the Boundary Element Method.

Analytic Solution to the MT Responses of a Two-Segment Model with Axially Anisotropic Conductivity Structures Overlying a Perfect Insulator

Magnetotelluric Responses of an Anisotropic 1-D Earth with a Layer of Exponentially Varying Conductivity

Modeling the 3D Terrain Effect on MT by the Boundary Element Method

Magnetotelluric Responses of a Vertical Inhomogeneous and Anisotropic Resistivity Structure with a Transitional Layer

Quasi-analytic Solution of 2-D Magnetotelluric Fields on an Axially Anisotropic Infinite Fault

Analytic Solution to the Magnetotelluric Response over Anisotropic Medium and Its Discussion

Theoretical magnetotelluric response of stratiform earth consisting of alternative homogeneous and transitional layers

Magnetotelluric Soundings on a Stratified Earth with Two Transitional Layers

Unstructured Grid Finite Element Modeling of the Three-Dimensional Magnetotelluric Responses in a Model With Arbitrary Conductivity and Magnetic Susceptibility Anisotropies

Transient Responses of a Magneto-Electro-elastic Hollow Sphere for Fully Coupled Spherically Symmetric Problem

The Transient Responses of a Special Non-Homogeneous Magneto-Electro-elastic Hollow Cylinder for Axisymmetric Plane Strain Problem

A novel Fourier domain scheme for three‐dimensional magnetotelluric modelling in anisotropic media

Three-Dimensional Unstructured Finite Element Modeling of Magnetotelluric Problems Allowing for Continuous Variation of Conductivity in Each Block

A Review of Subsurface Electrical Conductivity Anomalies in Magnetotelluric Imaging

A damped forward EMI model for a horizontally stratified earth

Investigating the Magnetotelluric Responses in Electrical Anisotropic Media

An Approach to Model Earth Conductivity Structures with Lateral Changes for Calculating Induced Currents and Geoelectric Fields during Geomagnetic Disturbances

Three-Dimensional Marine Magnetotelluric Parallel Forward Modeling in Conductive and Magnetic Anisotropic Medium Using Finite-Element Method Based on Secondary Field

Sensitivity study on two-dimensional and three-dimensional magnetotelluric inversions