DC3DPAFEM: An efficient and accurate 3-D direct current resistivity anisotropic forward modeling software for complex geological settings

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.