Abstract:Magnetotelluric (MT) exploration has been widely used in geological survey, mineral and oil and gas exploration, and the research on deep structures of the earth. With the models in the MT method becoming more and more complicated, conventional modeling methods based on structured grids cannot satisfy the accuracy requirement. To better simulate the complex earth, we adopt unstructured grids in the MT modeling. Since an effective mesh can greatly improve the modeling accuracy and save the cost, we present a goal-oriented unstructured adaptive finite element method for MT modeling. Starting from the frequency-domain Maxwell's equations, a vector Helmholtz equation is obtained for a stable solution with existing air. We use the Galerkin method to discretize the Helmholtz equation. The continuity condition for normal current density is used to evaluate the posterior error, while the weighted coefficient is obtained by solving a dual forward modeling problem. Besides, we multiply the weighted coefficients proposed by previous researches by the conductivity that proves to be more effective for grid refinement. Finally, we define a convergence rule to determine receivers used to calculate the weighted coefficients in the next refinement. The response of a homogeneous half-space model is simulated by both the goal-oriented and global adaptive finite element methods to verify the accuracy and to demonstrate advantages of our goal-oriented adaptive method over the global one. We further calculate MT response for a classical model to check the effectiveness of our method on topography models. Then, we calculate the response for a topographic model with an abnormal body embedded to demonstrate that our method can not only refine meshes near the source on the topographic earth surface but also refine meshes on physical interfaces between the abnormal body and surrounding medium. From the adaptive meshes and modeling results we can draw the following conclusions: 1) Our goal-oriented adaptive finite-element algorithm for MT modeling can obtain higher accuracy with less grids in comparison to the global method. 2) Our algorithm is effective for MT topographic modeling. 3) The algorithm presented in this paper can not only effectively refine the meshes close to the topographic surface with receivers, but also refine the meshes near the abnormal body in the earth. 4) The improved weighted coefficient is more effective than the traditional one in terms of grids refinement.

Highly Efficient Numerical Modeling for Static-Steady Fields in Underground Inhomogeneous Media

Numerical simulation of cutterhead and soil interaction in slurry shield tunneling

Numerical solution of Green's function of potential field of eccentric point source in axisymmetric inhomogeneous media

A General Method for Numerical Green'S Function in Arbitrarily Layered Soils

Numerical Green's Function of a Point Current Source in Horizontal Multilayer Soils by Utilizing the Vector Matrix Pencil Technique

Numerical Analysis Of Thin Skin Depths Of 3-D Eddy-Current Problems Using A Combination Of Finite Element And Meshless Methods

A Well-Balanced Numerical Model for the Simulation of Long Waves over Complex Domains

Numerical simulation of acoustic waves propagation by finite element method based on optimized matrices

A robust numerical method for modeling ventilation through long tunnels in high temperature regions based on 1D pipe model

Finite Analytic Numerical Method for Two-Dimensional Fluid Flow in Heterogeneous Porous Media

3D Vector Finite-Element Electromagnetic Forward Modeling for Large Loop Sources Using a Total-Field Algorithm and Unstructured Tetrahedral Grids

Efficient Triple-Grid Multiscale Finite Element Method for 3d Groundwater Flow Simulation in Heterogeneous Porous Media

Scalable Semismooth Newton Methods with Multilevel Domain Decomposition for Subsurface Flow and Reactive Transport in Porous Media

3D finite-difference transient electromagnetic modeling with a whole-space initial field

A Simplified Numerical Model as a Design Tool for Vertical Single U-Tube Ground Heat Exchangers

Precision of Meshfree Methods and Application to Forward Modeling of Two-Dimensional Electromagnetic Sources

Numerical simulation of cross-hole seismic wave field

A goal-oriented adaptive algorithm for 3D magnetotelluric forward modeling

A general numerical method for solving the three-dimensional hyperbolic heat conduction equation on unstructured grids

Numerical simulation on staggered grids of three-dimensional brinkman-forchheimer flow and heat transfer in porous media

Improved Integral Equation Method for Rapid 3-D Forward Modeling of Magnetotelluric