Application of the ${\chi }^2$ principle and unbiased predictive risk estimator for determining the regularization parameter in 3D focusing gravity inversion

Abstract

The ${\chi }^2$ principle and the unbiased predictive risk estimator are used to determine optimal regularization parameters in the context of 3D focusing gravity inversion with the minimum support stabilizer. At each iteration of the focusing inversion the minimum support stabilizer is determined and then the fidelity term is updated using the standard form transformation. Solution of the resulting Tikhonov functional is found efficiently using the singular value decomposition of the transformed model matrix, which also provides for efficient determination of the updated regularization parameter each step. Experimental 3D simulations using synthetic data of a dipping dike and a cube anomaly demonstrate that both parameter estimation techniques outperform the Morozov discrepancy principle for determining the regularization parameter. Smaller relative errors of the reconstructed models are obtained with fewer iterations. Data acquired over the Gotvand dam site in the south-west of Iran are used to validate use of the methods for inversion of practical data and provide good estimates of anomalous structures within the subsurface.

References in corpus (1)

• Regularization Parameter Estimation for Underdetermined problems by the ${\chi }^2$ principle with application to $2D$ focusing gravity inversion

Cited by in corpus (7)

• Hybrid and iteratively reweighted regularization by unbiased predictive risk and weighted GCV for projected systems
• 3-D Projected $L_1$ inversion of gravity data
• A fast algorithm for regularized focused 3-D inversion of gravity data using the randomized SVD
• Total variation regularization of the $3$-D gravity inverse problem using a randomized generalized singular value decomposition
• A fast methodology for large-scale focusing inversion of gravity and magnetic data using the structured model matrix and the $2D$ fast Fourier transform
• IGUG: A MATLAB package for $3$D inversion of gravity data using graph theory
• Convergence of Regularization Parameters for Solutions Using the Filtered Truncated Singular Value Decomposition