3-D Projected $L_1$ inversion of gravity data

Abstract

Sparse inversion of gravity data based on $L_1$-norm regularization is discussed. An iteratively reweighted least squares algorithm is used to solve the problem. At each iteration the solution of a linear system of equations and the determination of a suitable regularization parameter are considered. The LSQR iteration is used to project the system of equations onto a smaller subspace that inherits the ill-conditioning of the full space problem. We show that the gravity kernel is only mildly to moderately ill-conditioned. Thus, while the dominant spectrum of the projected problem accurately approximates the dominant spectrum of the full space problem, the entire spectrum of the projected problem inherits the ill-conditioning of the full problem. Consequently, determining the regularization parameter based on the entire spectrum of the projected problem necessarily over compensates for the non-dominant portion of the spectrum and leads to inaccurate approximations for the full-space solution. In contrast, finding the regularization parameter using a truncated singular space of the projected operator is efficient and effective. Simulations for synthetic examples with noise demonstrate the approach using the method of unbiased predictive risk estimation for the truncated projected spectrum. The method is used on gravity data from the Mobrun ore body, northeast of Noranda, Quebec, Canada. The $3$-D reconstructed model is in agreement with known drill-hole information.

References in corpus (2)

• Hybrid and iteratively reweighted regularization by unbiased predictive risk and weighted GCV for projected systems
• Application of the ${\chi }^2$ principle and unbiased predictive risk estimator for determining the regularization parameter in 3D focusing gravity inversion

Cited by in corpus (9)

• Generalized L${}_p$-norm joint inversion of gravity and magnetic data using cross-gradient constraint
• 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
• Inversion of multiconfiguration complex EMI data with minimum gradient support regularization: A case study
• Large-scale focusing joint inversion of gravity and magnetic data with Gramian constraint
• 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
• Unique reconstruction of simple magnetizations from their magnetic potential