Abstract:Forward modeling is the basis of inversion imaging and quantitative interpretation for DC resistivity exploration. Currently, a numerical model of the DC resistivity method must be finely divided to obtain a highly accurate solution under complex conditions, resulting in a long calculation time and large storage. Therefore, we propose a 3D numerical simulation method in a mixed space-wavenumber domain to overcome this challenge. The partial differential equation about abnormal potential is transformed into many independent ordinary differential equations with different wavenumbers using a 2D Fourier transform along the x axis and y axis direction. In this way, a large-scale 3D numerical simulation problem is decomposed into several 1D numerical simulation problems, which significantly reduces the computational and storage requirements. In addition, these ordinary 1D differential equations with different wavenumbers are independent of each other and high parallelelism of the algorithm. They are solved using a finite-element algorithm combined with a chasing method, and the obtained solution is modified using a contraction operator. In this method, the vertical direction is reserved as the spatial domain, then grid size can be determined flexibly based on the underground current density distribution, which considers the solution accuracy and calculation efficiency. In addition, for the first time, we use the contraction operator in the integral equation method to iterate the algorithm. The algorithm takes advantage of the high efficiency of the standard Fourier transform and chasing method, as well as the fast convergence of the contraction operator. We verified the accuracy of the algorithm and the convergence of the contraction operator. Compared with a volume integral method and goal-oriented adaptive finite-element method, the proposed algorithm has lower memory requirements and high computational efficiency, making it suitable for calculating a model with large-scale nodes. Moreover, different examples are used to verify the high adaptability and parallelism of the proposed algorithm. The findings show that the 3D numerical simulation method of DC resistivity method in a mixed space-wavenumber domain is highly efficient, precise, and parallel.

On the Phase-Field Algorithm for Distinguishing Connected Regions in Digital Model

Phase-Field Simulation of Three-Dimensional Dendritic Growth

Enhanced conservative phase field method for moving contact line problems

Consistently and unconditionally energy-stable linear method for the diffuse-interface model of narrow volume reconstruction

Multiphase image segmentation using a phase-field model

Second-order accurate and unconditionally stable algorithm with unique solvability for a phase-field model of 3D volume reconstruction

A novel segmentation algorithm for complex 3D mesh model in computer vision.

Zone Partitioning Sequential Calculation Method for Three-Dimensional Phase-Field Material Simulation During Directional Solidification

Efficiently Solving the Piecewise Constant Mumford-Shah Model using Graph Cuts

Numerical Simulation of Three-Dimensional Dendritic Growth Using Phase-Field Method

Multi-Phase Segmentation Using Modified Complex Cahn-Hilliard Equations

3D Adaptive Finite Element Method for a Phase Field Model for the Moving Contact Line Problems

On the phase field based model for the crystalline transition and nucleation within the Lagrange multiplier framework

A Modified and Efficient Phase Field Model for the Biological Transport Network

An End-to-End Deep Learning Method for Solving Nonlocal Allen-Cahn and Cahn-Hilliard Phase-Field Models

On the Stability and Accuracy of Partially and Fully Implicit Schemes for Phase Field Modeling

A Multi-Mesh Adaptive Finite Element Approximation to Phase Field Models

Numerical modeling of 3D DC resistivity method in the mixed space-wavenumber domain

Phase Field Modeling and Numerical Algorithm for Two-Phase Dielectric Fluid Flows

A ternary phase-field model for two-phase flows in complex geometries

Discrete Maximum Principle and Energy Stability of the Compact Difference Scheme for Two-Dimensional Allen-Cahn Equation