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.

INVERSION OF GRAVITY ANOMALY FOR MULTI-LAYERED DENSITIES AND THEIR INTERFACES BY RIDGE REGRESSION PROCEDURE

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

THE DENSITY INTERFACE INVERSION OF HIGH-PRECISION GRAVITY DATA

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

Application of Apparent Density Inversion Method in the East China Sea and Its Adjacent Area

Inversion of Gravity Anomaly Based on AlexNet Deep Neural Network

Stable solution of normalized total gravity gradient

Stable Inversion of Apparent Density in Undulating Layers in Frequency Domain

The Inversion of Two-Dimensional Mono-Density Interface with Several Known Control Points

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

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

Joint Inversion of Gravity Gradient Tensor Data Based on L1 and L2 Norms

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

Structure-guided gravity inversion for layered density modeling with an application in the Chezhen Depression, Bohai Bay Basin

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

Gravity Inversion in Spherical Coordinates With Dynamic Re-Weighting Matrix

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

Study on multi-scale density inversion of gravity anomalies

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

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

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