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.

Stable solution of normalized total gravity gradient

THE DENSITY INTERFACE INVERSION OF HIGH-PRECISION GRAVITY DATA

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

Numerical Modeling of Lock-Exchange Gravity/turbidity Currents by a High-Order Upwinding Combined Compact Difference Scheme

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

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

Gravity Inversion in Spherical Coordinates With Dynamic Re-Weighting Matrix

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

Inversion of Gravity Anomaly Based on AlexNet Deep Neural Network

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

Solving Electromagnetic Inverse Problem Using Adaptive Gradient Descent Algorithm.

Euler Deconvolution of the Analytic Signals of the Gravity Gradient Tensor for the Horizontal Pipeline of Finite Length by Horizontal Cylinder Calculation

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

Normalising Flows for Bayesian Gravity Inversion

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

Joint Gravity and Gravity Gradient Inversion for Subsurface Object Detection

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

Global marine gravity gradient tensor inverted from altimetry-derived deflections of the vertical: CUGB2023GRAD

Gravity compression forward modeling and multiscale inversion based on wavelet transform

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