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.

A New Method for Spectral Analysis of the Potential Field and Conversion of Derivatives of Gravity Anomalies: Cosine Transform

A New Fast Method for Continuation of Potential Field on an Undulate Terrain

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

Transformations Between Aeromagnetic Gradients in Frequency Domain

Two-dimensional Magnetotelluric Modelling by the Boundary Element Method.

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

High-precision Fourier forward modeling of potential fields

A Fast Interpretation Method of Gravity Gradiometry Data Based on Magnetic Dipole Localization

On Application of Fractal Magnetization in Curie Depth Estimation from Magnetic Anomalies

Geophysical Potential Field Anomaly Separation Method with Optimal Mother Wavelet and Spatial Locating Multiresolution Analysis (MRA)

Correction Of Topographic Distortions In Potential-Field Data - A Fast And Accurate Approach

Determination of vertical derivative of gravity anomalous by upward continuation and Taylor series transform methods: application to the Southwest sub-basin of the East Vietnam Sea

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

Inversion of Gravity Anomaly Based on AlexNet Deep Neural Network

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

Wavelet Transform and Its Application to Decomposition of Gravity Anomalies

Gravity Gradient Parker′s Forward Method and Application Using Cosine Transform

3-Dimensional Inversion for Gravity Anomaly Calculation in Complex Geologic Region

Gravity Anomaly and Its Accuracy Assessment from Multi-source Altimetry Data in the Western Pacific Area

Apparent Depth Filter For Gravity Anomalies And A Case Study

A Signal-Strength Model for Gravity-Field Spectrum Components in a Local Space