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.

The geophysical model enhancement based on the convolution model

Reflection Multi‐scale Envelope Inversion

Two-dimensional magnetotelluric inversion method based on deep learning

Seismic data deconvolution using Kalman filter based on a new system model

Characterizing the impact of model error in hydrogeologic time series recovery inverse problems

Enhancing electrical structure in magnetotelluric inversion by the constraint of minimum cross-gradient support coupling

Data and Model Dual-Driven Seismic Deconvolution via Error-Constrained Joint Sparse Representation

Error effects on magnetotelluric inversion

Geological Structure-Guided Initial Model Building for Prestack AVO/AVA Inversion

Envelope-Based Sparse-Constrained Deconvolution for Velocity Model Building

High-resolution Seismic Impedance Inversion Integrating the Closed-Loop Convolutional Neural Network and Geostatistics: an Application to the Thin Interbedded Reservoir

Inversion of CSEM Data with Model-based Inversion Algorithm

A Trade-off Solution of Regularized Geophysical Inversion Using Model Resolution and Covariance Matrices

CSEM Inversions Using the Gauss Newton Method and a Model Compression Scheme

Deep-learning-based upscaling method for geologic models via theory-guided convolutional neural network

Joint Gravity and Magnetic Inversion Using CNNs' Deep Learning

Enhancing the seismic response of faults by using a deep learning‐based method

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

High-resolution seismic acoustic impedance inversion with the sparsity-based statistical model

3-D Time-Space Joint Deconvolution for Enhancing Recognition Accuracy of Seismic Microstructure

Inverting Incomplete Fourier Transforms by a Sparse Regularization Model and Applications in Seismic Wavefield Modeling