Research Note: A unifying framework for the widely used stabilization of potential field inverse problems

We present a brief review of the widely-used and well-known stabilizers in the inversion of potential field data. These include stabilizers that use L2, L1, and L0 norms of the model parameters, and the gradients of the model parameters. These stabilizers may all be realized in a common setting using two general forms with different weighting functions. Moreover, we show that this unifying framework encompasses the use of additional stabilizations which are not common for potential field inversion. INTRODUCTION Many forms of stabilizers are used for the inversion of the ill-posed potential field data problem. The aim is the reconstruction of subsurface models that provide relevant physical interpretation. In order to reduce the possibility of over-interpretation of the data, it may be desirable to reconstruct a simple model with as little structure as possible, and to eliminate arbitrary discontinuities in the solution (Constable et al., 1987). Such models can be expected to present only the important and large-scale features of the subsurface under the survey area. They can be obtained using a minimum roughness, or equivalently maximum smoothness, stabilizer that employs a L2-norm of the gradient of the model parameters in the inversion algorithm (Constable et al., 1987; Li and Oldenburg, 1996; Pilkington, 1997; Li and Oldenburg, 1998). A

[1]  D. Oldenburg,et al.  A comprehensive study of including structural orientation information in geophysical inversions , 2009 .

[2]  Douglas W. Oldenburg,et al.  3-D inversion of magnetic data , 1996 .

[3]  Per Christian Hansen,et al.  Analysis of Discrete Ill-Posed Problems by Means of the L-Curve , 1992, SIAM Rev..

[4]  Gene H. Golub,et al.  Generalized cross-validation as a method for choosing a good ridge parameter , 1979, Milestones in Matrix Computation.

[5]  R. Parker,et al.  Occam's inversion; a practical algorithm for generating smooth models from electromagnetic sounding data , 1987 .

[6]  D. Oldenburg,et al.  3-D inversion of gravity data , 1998 .

[7]  Rosemary A Renaut,et al.  Regularization parameter estimation for underdetermined problems by the χ 2 principle with application to 2D focusing gravity inversion , 2014, 1402.3365.

[8]  V. Morozov On the solution of functional equations by the method of regularization , 1966 .

[9]  Jiajia Sun,et al.  Adaptive Lp inversion for simultaneous recovery of both blocky and smooth features in a geophysical model , 2014 .

[10]  Saeed Vatankhah,et al.  A fast algorithm for regularized focused 3D inversion of gravity data using randomized singular-value decomposition , 2018, GEOPHYSICS.

[11]  Michael S. Zhdanov,et al.  Minimum support nonlinear parametrization in the solution of a 3D magnetotelluric inverse problem , 2004 .

[12]  Esben Auken,et al.  Generalized focusing of time-lapse changes with applications to direct current and time-domain induced polarization inversions , 2015 .

[13]  Rosemary A. Renaut,et al.  Automatic estimation of the regularization parameter in 2-D focusing gravity inversion: an application to the Safo manganese mine in northwest of Iran , 2013, ArXiv.

[14]  D. Oldenburg,et al.  NON-LINEAR INVERSION USING GENERAL MEASURES OF DATA MISFIT AND MODEL STRUCTURE , 1998 .

[15]  Rick Chartrand,et al.  Exact Reconstruction of Sparse Signals via Nonconvex Minimization , 2007, IEEE Signal Processing Letters.

[16]  Rosemary A. Renaut,et al.  Application of the χ2 principle and unbiased predictive risk estimator for determining the regularization parameter in 3-D focusing gravity inversion , 2014, 1408.0712.

[17]  Rosemary A. Renaut,et al.  A Newton root-finding algorithm for estimating the regularization parameter for solving ill-conditioned least squares problems , 2009 .

[18]  Per Christian Hansen,et al.  Regularization Tools version 4.0 for Matlab 7.3 , 2007, Numerical Algorithms.

[19]  Dominique Fournier,et al.  A cooperative magnetic inversion method with Lp-norm regularization , 2015 .

[20]  Rosemary A. Renaut,et al.  Total variation regularization of the 3-D gravity inverse problem using a randomized generalized singular value decomposition , 2017, 1709.08125.

[21]  Maurizio Fedi,et al.  Inversion of potential field data using the structural index as weighting function rate decay , 2012 .

[22]  D. Oldenburg,et al.  3-D inversion of gravity data , 1998 .

[23]  Brendt Wohlberg,et al.  An Iteratively Reweighted Norm Algorithm for Minimization of Total Variation Functionals , 2007, IEEE Signal Processing Letters.

[24]  Valéria C. F. Barbosa,et al.  Generalized compact gravity inversion , 1994 .

[25]  M. Chouteau,et al.  Constraints in 3D gravity inversion , 2001 .

[26]  C. Vogel Computational Methods for Inverse Problems , 1987 .

[27]  Per Christian Hansen,et al.  REGULARIZATION TOOLS: A Matlab package for analysis and solution of discrete ill-posed problems , 1994, Numerical Algorithms.

[28]  Peter G. Lelièvre,et al.  Forward modelling and inversion of geophysical magnetic data , 2003 .

[29]  Rosemary A. Renaut,et al.  3-D Projected L1 inversion of gravity data using truncated unbiased predictive risk estimator for regularization parameter estimation , 2017 .

[30]  D. Oldenburg,et al.  Incorporating geological dip information into geophysical inversions , 2000 .

[31]  Douglas W. Oldenburg,et al.  3-D inversion of magnetic data , 1996 .

[32]  Jonathan B. Ajo-Franklin,et al.  Applying Compactness Constraints to Differential Traveltime Tomography , 2007 .

[33]  K. Kubik,et al.  Compact gravity inversion , 1983 .

[34]  Michael S. Zhdanov,et al.  Focusing geophysical inversion images , 1999 .

[35]  Douglas W Oldenburg,et al.  Inversion using spatially variable mixed ℓp norms , 2019, Geophysical Journal International.

[36]  D. Oldenburg,et al.  Fast inversion of large-scale magnetic data using wavelet transforms and a logarithmic barrier method , 2003 .

[37]  Mark Pilkington,et al.  3-D magnetic imaging using conjugate gradients , 1997 .

[38]  Donald W. Marquaridt Generalized Inverses, Ridge Regression, Biased Linear Estimation, and Nonlinear Estimation , 1970 .

[39]  Elena Cherkaev,et al.  Non-smooth gravity problem with total variation penalization functional , 2002 .

[40]  A. Guillén,et al.  Gravity and magnetic inversion with minimization of a specific functional , 1984 .

[41]  Peter G. Lelièvre,et al.  Gradient and smoothness regularization operators for geophysical inversion on unstructured meshes , 2013 .