A fast methodology for large-scale focusing inversion of gravity and magnetic data using the structured model matrix and the $2D$ fast Fourier transform

Focusing inversion of potential field data for the recovery of sparse subsurface structures from surface measurement data on a uniform grid is discussed. For the uniform grid the model sensitivity matrices exhibit block Toeplitz Toeplitz block structure, by blocks for each depth layer of the subsurface. Then, through embedding in circulant matrices, all forward operations with the sensitivity matrix, or its transpose, are realized using the fast two dimensional Fourier transform. Simulations demonstrate that this fast inversion algorithm can be implemented on standard desktop computers with sufficient memory for storage of volumes up to size $n \approx 1M$. The linear systems of equations arising in the focusing inversion algorithm are solved using either Golub Kahan bidiagonalization or randomized singular value decomposition algorithms in which all matrix operations with the sensitivity matrix are implemented using the fast Fourier transform. These two algorithms are contrasted for efficiency for large-scale problems with respect to the sizes of the projected subspaces adopted for the solutions of the linear systems. The presented results confirm earlier studies that the randomized algorithms are to be preferred for the inversion of gravity data, and that it is sufficient to use projected spaces of size approximately $m/8$, for data sets of size $m$. In contrast, the Golub Kahan bidiagonalization leads to more efficient implementations for the inversion of magnetic data sets, and it is again sufficient to use projected spaces of size approximately $m/8$. Moreover, it is sufficient to use projected spaces of size $m/20$ when $m$ is large, $m \approx 50000$, to reconstruct volumes with $n \approx 1M$. Simulations support the presented conclusions and are verified on the inversion of a practical magnetic data set that is obtained over the Wuskwatim Lake region in Manitoba, Canada.

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

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

[3]  R. Blakely.  Potential theory in gravity and magnetic applications , 1996 .

[4]  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.

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

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

[7]  Bo Chen,et al.  High-accuracy 3D Fourier forward modeling of gravity field based on the Gauss-FFT technique , 2018 .

[8]  Valéria C. F. Barbosa,et al.  Interactive gravity inversion , 2006 .

[9]  Rosemary A Renaut,et al.  A Tutorial and Open Source Software for the Efficient Evaluation of Gravity and Magnetic Kernels , 2020, Comput. Geosci..

[10]  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.

[11]  Valéria C. F. Barbosa,et al.  Robust 3D Gravity Gradient Inversion By Planting Anomalous Densities , 2011 .

[12]  Lanbo Liu,et al.  Fast and accurate forward modelling of gravity field using prismatic grids , 2018, Geophysical Journal International.

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

[14]  Gene H. Golub,et al.  Matrix computations , 1983 .

[15]  Yau Shu Wong,et al.  BTTB-based numerical schemes for three-dimensional gravity field inversion , 2015 .

[16]  Mark Pilkington,et al.  3D magnetic data-space inversion with sparseness constraints , 2009 .

[17]  Xiong Li,et al.  Three‐dimensional gravity modeling in all space , 1997 .

[18]  D. Oldenburg,et al.  A comparison of automatic techniques for estimating the regularization parameter in non-linear inverse problems , 2004 .

[19]  Douglas W. Oldenburg,et al.  Magnetic forward modelling and inversion for high susceptibility , 2006 .

[20]  J. Nocedal.  Algorithms and software for large scale optimization. Progress report, August 1, 1989--May 5, 1990 , 1990 .

[21]  Hua Xiang,et al.  Regularization with randomized SVD for large-scale discrete inverse problems , 2013 .

[22]  Michael S. Zhdanov,et al.  3D inversion of airborne electromagnetic data using a moving footprint , 2010 .

[23]  M. Zhdanov,et al.  3‐D magnetic inversion with data compression and image focusing , 2002 .

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

[25]  Nathan Halko,et al.  Finding Structure with Randomness: Probabilistic Algorithms for Constructing Approximate Matrix Decompositions , 2009, SIAM Rev..

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

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

[28]  Michael A. Saunders,et al.  LSQR: An Algorithm for Sparse Linear Equations and Sparse Least Squares , 1982, TOMS.

[29]  Qian-Jiang Zhang,et al.  Fast 3D forward modeling of the magnetic field and gradient tensor on an undulated surface , 2018, Applied Geophysics.

[30]  Rosemary A. Renaut,et al.  Research Note: A unifying framework for the widely used stabilization of potential field inverse problems , 2020, Geophysical Prospecting.

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

[32]  M. Nabighian,et al.  Historical development of the gravity method in exploration , 2005 .

[33]  Rosemary A. Renaut,et al.  A fast algorithm for regularized focused 3-D inversion of gravity data using the randomized SVD , 2017 .

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

[35]  D. Rao,et al.  A rapid method for three‐dimensional modeling of magnetic anomalies , 1991 .

[36]  Rosemary A. Renaut,et al.  Improving the use of the randomized singular value decomposition for the inversion of gravity and magnetic data , 2019, GEOPHYSICS.

[37]  Daniel B. Szyld,et al.  An introduction to iterative Toeplitz solvers , 2009, Math. Comput..

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

[39]  Rosemary A. Renaut,et al.  Hybrid and iteratively reweighted regularization by unbiased predictive risk and weighted GCV for projected systems 5 December 2016 , 2016 .

[40]  Colin Farquharson,et al.  Constructing piecewise-constant models in multidimensional minimum-structure inversions , 2008 .

[41]  Saeed Vatankhah,et al.  IGUG: A MATLAB package for 3D inversion of gravity data using graph theory , 2018, Comput. Geosci..

[42]  Tristan van Leeuwen,et al.  Comparing RSVD and Krylov methods for linear inverse problems , 2020, Comput. Geosci..

[43]  Sergey Voronin,et al.  Compression approaches for the regularized solutions of linear systems from large-scale inverse problems , 2014, 1404.5684.