Abstract:Nowadays, there is a growing trend that direct current (DC) field surveys are shifting towards challenging areas characterized by mountainous topography and electrical anisotropy. Given these complex geological settings, there is an urgent need for 3-D DC forward modeling software capable of effectively addressing large-scale problems and delivering accurate modeling results to interpret field data. However, most open -source software packages face certain limitations, such as the high numerical cost to handle complex surface topography, the lack of consideration for anisotropic conductivity, the absence of mesh refinement techniques to guarantee accuracy in forward modeling, and the lack of parallel computing techniques to solve large-scale problems. In this study, we develop an efficient and highly accurate 3-D DC anisotropic forward modeling software, namely DC3DPAFEM, using the adaptive finite element algorithm based on the unstructured tetrahedral mesh. Firstly, we construct a strong compatible boundary value problem (BVP) for 3-D anisotropic DC problems by adopting a specialized secondary potential approach to handle the surface topography efficiently. Then, we develop a goal -oriented adaptive mesh refinement (AMR) technique to ensure accurate forward modeling results, even with a coarse initial mesh. To ensure time and memory efficiency, we employ a robust conjugate gradient (CG) algorithm preconditioned by the algebraic multigrid (AMG) solver to solve the large-scale linear system of equations resulting from complex geological structures. We aim to investigate the performance of the AMG scheme in anisotropic DC cases. Furthermore, we incorporate the domain decomposition technique into the iterative solution scheme for further efficiency gains. This technique significantly improves computing efficiency for large-scale problems in parallel clusters. Finally, we conduct comprehensive performance tests for DC3DPAFEM using a two -layer anisotropic model and a 3-D complex model with undulating terrain. The results of both examples validate the accuracy of DC3DPAFEM, as they closely align with the analytical solutions and the solutions obtained from the existing 3-D DC forward modeling code. Compared to traditional direct solver MUMPS and ILU-preconditioned iterative solvers, DC3DPAFEM exhibits highly scalable performance for large-scale problems, offering significant advantages in terms of memory and time consumption. Overall, DC3DPAFEM demonstrates substantial advances in efficiency, accuracy, and practicality through a series of numerical examples. This open -source code provides an efficient and available tool for developing a 3-D DC inversion method that can deal with large-scale problems involving intricate topography and anisotropic media.

3d Dcip Inversion Using The Finite-Element Approach With Variable Elements

Efficient Inversion of 3D IP Data for Continuous Model with Complex Geometry

Three-Dimensional Induced Polarization Parallel Inversion Using Nonlinear Conjugate Gradients Method

2.5D Direct‐current Resistivity Forward Modelling and Inversion by Finite‐element–infinite‐element Coupled Method

Three-Dimensional Inversion of Borehole-Surface Resistivity Method Based on the Unstructured Finite Element

Three-Dimensional Anisotropic Inversion for Direct Current Resistivity Method Using Unstructured Tetrahedral Discretization

Effective 3D-Transient Electromagnetic Inversion Using Finite-Element Method with a Parallel Direct Solver

Parallelized 3-D Csem Inversion With Secondary Field Formulation And Hexahedral Mesh

3D transient electromagnetic inversion based on explicit finite-difference forward modeling

2.5D Induced Polarization Forward Modeling Using the Adaptive Finite-Element Method

3-D Parallel Anisotropic Inversion of Controlled-Source Electromagnetic Data Using Nested Tetrahedral Grids

3-D CSEM inversion for the complex model with topography using an efficient dual parallel approach

3D parallel anisotropic inversion of controlled-source electromagnetic data using nested tetrahedral grids

Two-dimensional Inversion of Spectral Induced Polarization Data Using MPI Parallel Algorithm in Data Space

Parallelized 3-D Inversion of Controlled-Source Electromagnetic Data Based on Spectral Element Method with Infinite Element Boundary

Three-Dimensional Inversion of CSEM Data Using Finite Element Method in Data Space

Two-dimensional synchronous inversion of TDIP with cross-gradient constraint

3D FEM Numerical Forward Modeling of Direct Current Electrical Resistivity Based on PCG Algorithm

DC3DPAFEM: an Efficient and Accurate 3-D Direct Current Resistivity Anisotropic Forward Modeling Software for Complex Geological Settings

Inversion Imaging of 3D Resistivity Detection Using Adaptive-Weighted Smooth Constraint and PCG Algorithm

Three-dimensional Land FD-CSEM Regularized Inversion Based on Edge Finite-element Method