Abstract:When using traditional Euler deconvolution optimization strategies, it is difficult to distinguish between anomalies and their corresponding Euler tails (those solutions are often distributed outside the anomaly source, forming "tail"-shaped spurious solutions, i.e., misplaced Euler solutions, which must be removed or marked) with only the structural index. The nonparametric estimation method based on the normalized B-spline probability density (BSS) is used to separate the Euler solution clusters and mark different anomaly sources according to the similarity and density characteristics of the Euler solutions. For display purposes, the BSS needs to map the samples onto the estimation grid at the points where density will be estimated in order to obtain the probability density distribution. However, if the size of the samples or the estimation grid is too large, this process can lead to high levels of memory consumption and excessive computation times. To address this issue, a fast linear binning approximation algorithm is introduced in the BSS to speed up the computation process and save time. Subsequently, the sample data are quickly projected onto the estimation grid to facilitate the discrete convolution between the grid and the density function using a fast Fourier transform. A method involving multivariate B-spline probability density estimation based on the FFT (BSSFFT), in conjunction with fast linear binning appropriation, is proposed in this paper. The results of two random normal distributions show the correctness of the BSS and BSSFFT algorithms, which is verified via a comparison with the true probability density function (pdf) and Gaussian kernel smoothing estimation algorithms. Then, the Euler solutions of the two synthetic models are analyzed using the BSS and BSSFFT algorithms. The results are consistent with their theoretical values, which verify their correctness regarding Euler solutions. Finally, the BSSFFT is applied to Bishop 5X data, and the numerical results show that the comprehensive analysis of the 3D probability density distributions using the BSSFFT algorithm, derived from the Euler solution subset of x0,y0,z0, can effectively separate and locate adjacent anomaly sources, demonstrating strong adaptability.

3-D Kernel Density Imaging Based on the Euler Deconvolution of Tensor Gravity Data

3-D probability density imaging of Euler solutions using gravity data: a case study of Mount Milligan, Canada

Kernel Density Derivative Estimation of Euler Solutions

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

The inversion of apparent density based on the separation and continuation of potential field

Ratio‐Euler Deconvolution and Its Applications

3D density imaging using gravity and gravity gradient in the wavenumber domain and its application in the Decorah

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

3-D data-space joint inversion of gravity and magnetic data using a correlation-analysis constraint

Frechet Kernels for Imaging Regional Earth Structure Based on Three-Dimensional Reference Models

Improved Approaches for 3D Gravity and Gradient Imaging Based on Potential Field Separation: Application to the Magma Chamber in Wudalianchi Volcanic Field, Northeastern China

FFT-Based Probability Density Imaging of Euler Solutions

3-D Time-Space Joint Deconvolution for Enhancing Recognition Accuracy of Seismic Microstructure

3-D Gravity Inversion Based on Cross Dimensional Deep Neural Networks

A Decomposed Karhunen–Loève Expansion Scheme for the Discretization of Multidimensional Random Fields in Geotechnical Variability Analysis

Subsurface geology detection from application of the gravity-related dimensionality constraint

Gravity Inversion in Spherical Coordinates With Dynamic Re-Weighting Matrix

Gravity and Magnetic Data Inversion With Random Projection

An Adaptive Conjugate Gradient Least-Squares Regularization (ACGLSR) Method for 3D Gravity Density Inversion

Gravity and magnetic joint imaging based on Gramian constraints

3D gravity inversion based on deep learning