Magnetic inversion to recover the subsurface block structures based onL1 norm and total variation regularization
- 1 September 2021
- journal article
- research article
- Published by Oxford University Press (OUP) in Geophysical Journal International
- Vol. 228 (1), 510-537
- https://doi.org/10.1093/gji/ggab355
Abstract
This paper presents a new sparse inversion method based on L1 norm regularization for 3-D magnetic data. In isolation, L1 norm regularization yields model elements which are unconstrained by the input data to be exactly zero, leading to a sparse model with compact and focused structure. Here, we complement the L1 norm with a penalty minimizing total variation, the L1 norm of the model gradients; it is expected that the sharp boundaries of the subsurface structure are not compromised by incorporating this penalty. Although this penalty is widely used in the geophysical inversion studies, it is often replaced by an alternative quadratic penalty to ease solution of the penalized inversion problem; in this study, the original definition of the total variation, that is form of the L1 norm of the model gradients, is used. To solve the problem with this combined penalty of L1 norm and total variation, this study introduces alternative direction method of multipliers, which is a primal-dual optimization algorithm that solves convex penalized problems based on the optimization of an augmented Lagrange function. To improve the computational efficiency of the algorithm to make this method applicable to large-scale magnetic inverse problems, this study applies matrix compression using the wavelet transform and the preconditioned conjugate gradient method. The inversion method is applied to both synthetic tests and real data, the synthetic tests demonstrate that, when subsurface structure is blocky, it can be reproduced almost perfectly.Funding Information
- JSPS (JP19K04967)
This publication has 57 references indexed in Scilit:
- The solution path of the generalized lassoThe Annals of Statistics, 2011
- Reduction of Magnetic Anomaly Observations from Helicopter Surveys at Varying ElevationsExploration Geophysics, 2006
- A comparison of smooth and blocky inversion methods in 2D electrical imaging surveysExploration Geophysics, 2003
- Efficient determination of multiple regularization parameters in a generalized L-curve frameworkInverse Problems, 2002
- Constraints in 3D gravity inversionGeophysical Prospecting, 2001
- Early Cretaceous Sr‐rich silicic magmatism by slab melting in the Kitakami Mountains, northeast JapanJournal of Geophysical Research, 1994
- Partialling out the Spatial Component of Ecological VariationEcology, 1992
- A theory for multiresolution signal decomposition: the wavelet representationIEEE Transactions on Pattern Analysis and Machine Intelligence, 1989
- Orogeny and relative plate motions: Example of the Japanese IslandsTectonophysics, 1986
- Petrochemistry of the Sakainokami plutonic body, Kitakami mountains, northeastern Japan.The Journal of the Japanese Association of Mineralogists, Petrologists and Economic Geologists, 1976