Abstract:To solve the non-uniqueness problem of gravity gradient inversion, we proposed a folding calculation method based on preconditioned conjugate gradient inversion. Compared with the original algorithm in which the entire study area is taken as the research subject and all grids are used simultaneously in the inversion, the proposed folding method divides the research area into several sub-areas. A prism unit from any of the four corner grids is selected for the first iteration, whose density anomaly result is taken as the initial density anomaly for the next iteration of the same sub-area. The folding in the left–right and up–down directions takes turns during the calculation until the inversion calculation has covered the entire research area. This folding algorithm demonstrates strong regularity. The inversion results of multiple synthetic models show that the folding calculation method performs multiple parameter corrections in the initial model. Meanwhile, on the basis that the relative fitting of standard errors of the observed values satisfies the convergence condition, the model errors are constrained, and thus, the model errors of the inversion results are consequently reduced. The 3-D density anomaly pattern over the Craton area in North China was obtained via a joint inversion of the four components (\(\varvec{T}_{\text{xx}} ,\) \(\varvec{T}_{\text{xz}} ,\) \(\varvec{T}_{\text{yy}} ,\) and \(\varvec{T}_{\text{zz}}\)) of the Gravity field and steady-state Ocean Circulation Explorer L2 gravity gradient after preprocessing. We compared the inversion results from the folding calculation method and the original method, and performed detailed analysis and discussion on the inversion results with existing geological and geophysical data. Our analysis shows that the improved calculation method is effectively applicable to the inversion of measured gravity gradient data, and the inversion results provide more detailed and reliable pattern information for the density anomaly.

Adaptive Mesh-Free Approach for Gravity Inversion Using Modified Radial Basis Function

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

An Improved Meshfree Scheme Based on Radial Basis Functions for Solving Incompressible Navier–Stokes Equations

Stable solution of normalized total gravity gradient

THE DENSITY INTERFACE INVERSION OF HIGH-PRECISION GRAVITY DATA

Gravity Inversion in Spherical Coordinates With Dynamic Re-Weighting Matrix

Inversion of Gravity Anomaly Based on AlexNet Deep Neural Network

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

Lévy Gradient Descent: Augmented Random Search for Geophysical Inverse Problems

A local level-set method for 3D inversion of gravity-gradient data

Global Gravity Inversion of Bodies with Arbitrary Shape

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

Highly efficient density inversion of gravity data using nonlinear density polynomial fitting

Joint Gravity and Gravity Gradient Inversion for Subsurface Object Detection

A Folding Calculation Method Based on the Preconditioned Conjugate Gradient Inversion Algorithm of Gravity Gradient Tensor

A Joint Inversion Algorithm for the Establishment of High-Accuracy 3-D Marine Gravity Field

3D inversion of magnetic gradient data based on equivalent source weighting method

Joint Inversion of Gravity and Gravity Gradient Data Using Smoothed L0 Norm Regularization Algorithm with Sensitivity Matrix Compression

Joint Inversion of Gravity and Gravity Gradient Data Based on Cross-Gradient Function

3D gravity inversion based on deep learning