Three-dimension holographic numerical simulation of gravity anomalies

Abstract

In gravity-anomaly-based prospecting, it is difficult to achieve large-scale and high-precision inversion imaging of complex geological models. To address this issue, this paper proposes a three-dimensional holographic numerical simulation method for gravity anomalies. This method transforms the 3D partial differential equation of gravity potential into many independent 1D differential equations with different wavenumbers by performing a 2D Fourier transform along the horizontal direction. It decomposes a large-scale 3D numerical modeling problem into many 1D numerical modeling problems, which greatly reduces the computation and memory requirements of modeling, and each 1D differential equation is independent, so it has high parallelism. The vertical direction is reserved as the spatial domain, which has strict upper and lower boundary conditions; The horizontal 2D Fourier transform uses a holographic Fourier transform (Holo-FT) with a full interval of integration, consistent with the physical boundary in the horizontal direction, thus the physical information of the gravity potential described by the three-dimensional partial differential equation is fully and accurately simulated. Using the high parallelism of the algorithm, the CPU is used for parallel solving ordinary differential equations, while the GPU is used for parallel computing of Holo-FT, which implements the CPU-GPU parallel acceleration scheme and further improves the efficiency of this algorithm.