Abstract:In order to improve the efficiency and accuracy of numerical simulation of gravity anomalies in large-scale complex models, this paper proposes a method for three-dimensional numerical simulation of gravity anomalies under arbitrary terrain, and implements its CPU-GPU parallel acceleration scheme. By performing a two-dimensional Fourier transform along the horizontal direction, three-dimensional partial differential satisfied by gravitational potential are transformed into one-dimensional ordinary differential equations in different wavenumbers, reducing computational and storage requirements. Moreover, the ordinary differential equations in different wavenumbers are independent of each other and have good parallelism. By retaining the vertical direction as the spatial domain and introducing traveling wave decomposition, the effects caused by upper and lower boundaries are eliminated. The vertical grid is arbitrary. A two-dimensional quasi-complete information Fourier transform is utilized to enhance both the accuracy and efficiency of the Fourier transform process, allowing both uniform and non-uniform flexible sampling. Based on the high parallelism, one-dimensional ordinary differential equations are solved in parallel using CPU, and quasi-complete information Fourier transform are computed in parallel using GPU. An abnormal sphere is designed to analyze the wavenumber spectral distribution characteristics of the anomalous potential and summarize the wavenumber sampling rules. The logarithmic uniform sampling in the wavenumber domain has high accuracy. Compared with Gauss-FFT, the quasi-complete information Fourier transform selects fewer wavenumbers and exhibits higher accuracy. The computational efficiency is greatly enhanced through the use of single-precision CPU-GPU parallel scheme, and a model with tens of millions of nodes only takes a few seconds. Two terrain models are designed to verify the algorithm's adaptability to arbitrary complex terrains. This is significant for achieving fine inversion imaging and interpretation of gravity anomalies under large-scale complex conditions.

Efficient 2D Modeling of Gravity Anomaly in Space-wavenumber Mixed Domain

3D Inversion of Gravity Anomaly – an Application of Technology of Apparent Density Imagery

Fast Numerical Simulation of 2D Gravity Anomaly Based on Nonuniform Fast Fourier Transform in Mixed Space-Wavenumber Domain

Generalized Gauss-FFT 3D forward gravity modeling for irregular topographic mass having any 3D variable density contrast

3D Gravity Inversion in Frequency Domain Based on Mixed L1+L2 Norm

Efficient 3D Gravity and Magnetic Modeling

A Fast Forward and Inversion Strategy for Three-Dimensional Gravity Field

High-accuracy 3D Fourier Forward Modeling of Gravity Field Based on the Gauss-FFT Technique

High-accuracy Gravity Field and Gravity Gradient Forward Modelling Based on 3D Vertex-Centered Finite-Element Algorithm

Three-dimensional modeling of magnetic anomaly integral solution in a mixed space-wavenumber domain

High-Efficiency Forward Modeling of Gravitational Fields in Spherical Harmonic Domain with Application to Lunar Topography Correction

Fast 3D gravity and magnetic modelling using midpoint quadrature and 2D FFT

Fast 3-D Magnetic Anomaly Forward Modeling Based on Integral Equation

Gravity compression forward modeling and multiscale inversion based on wavelet transform

Three-dimensional quasi-complete information numerical simulation of gravity anomalies and Parallel Computing on CPU-GPU

Gravity Inversion in Spherical Coordinates With Dynamic Re-Weighting Matrix

Integrated Gravity and Gravity Gradient 3D Inversion Using the Non-Linear Conjugate Gradient

Efficient 3-D Large-Scale Forward Modeling and Inversion of Gravitational Fields in Spherical Coordinates with Application to Lunar Mascons

Spectral-Element Simulations of Gravity Anomalies for 3-D Topography Model

Integrated gravity and gravity gradient data focusing inversion

An Efficient and Adaptive Approach for Modeling Gravity Effects in Spherical Coordinates