Joint inversion of gravity and multiple components of tensor gravity data

Abstract

Geological structures often exhibit smooth characteristics away from sharp discontinuities. One aim of geophysical inversion is to recover information about the smooth structures as well as about the sharp discontinuities. Because no specific operator can provide a perfect sparse representation of complicated geological models, hyper-parameter regularization inversion based on the iterative split Bregman method was used to recover the features of both smooth and sharp geological structures. A novel preconditioned matrix was proposed, which counteracted the natural decay of the sensitivity matrix and its inverse matrix was calculated easily. Application of the algorithm to synthetic data produces density models that are good representations of the designed models. The results show that the algorithm proposed is feasible and effective.

Keywords

hyper-parameter regularization / full gravity gradient tensor / preconditioned matrix / Occam’s inversion / focusing inversion

References

[1] WangJ-b, XiongS-q, ZhouX-h, GuoZ-hong. The advances in the study of the airborne gravimetry system [J]. Geophysical and Geochemical Exploration, 2009, 33(4): 368-373

[2] ÁlvarezO, GimenezM, BraitenbergC, AndresF. Goce satellite derived gravity and gravity gradient corrected for topographic effect in the south central andes region [J]. Geophysical Journal International, 2012, 190(2): 941-959

[3] LiY-g, OldenburgD W. 3-D inversion of gravity data [J]. Geophysics, 1998, 63(1): 109-119

[4] ZhdanovM, EllisR, MukherjeeS. Three-dimensional regularized focusing inversion of gravity gradient tensor component data [J]. Geophysics, 2004, 69(4): 925-937

[5] FarquharsonC G. Constructing piecewise-constant models in multi-dimensional minimum-structure inversions [J]. Geophysics, 2008, 73(1): K1-K9

[6] VascoD, TaylorC. Inversion of airborne gravity gradient data, southwestern Oklahoma [J]. Geophysics, 1991, 56(1): 90-101

[7] LuW-t, QianJ-liang. A local level-set method for 3d inversion of gravity-gradient data [J]. Geophysics, 2015, 80(1): G35-G51

[8] O’BrienJ, RodriguezA, SixtaD, DaviesM, HoughtonP. Resolving the k-2 salt structure in the gulf of mexico an integrated approach using prestack depth imaging and full tensor gravity gradiometry [J]. The Leading Edge, 2005, 24(4): 404-409

[9] ZhangL-l, YuP, WangJ-l, WuJ-sheng. Smoothest model and sharp boundary based two-dimensional magnetotelluric inversion [J]. Chinese Journal of Geophysics, 2009, 52(6): 1625-1632

[10] WrightS JPrimal-dual interior-point methods[M], 1997PhiladelphiaSociety for Industrial and Applied Mathematics

[11] GoldsteinT, OsherS. The split bregman method for L1-regularized problems [J]. SIAM Journal on Imaging Sciences, 2009, 2(2): 323-343

[12] JiaR-q, ZhaoH-q, ZhaoWei. Convergence analysis of the bregman method for the variational model of image denoising [J]. Applied and Computational Harmonic Analysis, 2009, 27(3): 367-379

[13] GholamiA, SiahkoohiH. Regularization of linear and non-linear geophysical ill-posed problems with joint sparsity constraints [J]. Geophysical Journal International, 2010, 180(2): 871-882

[14] KirkendallB, LiY G, OldenburgD. Imaging cargo containers using gravity gradiometry [J]. Geoscience and Remote Sensing, IEEE Transactions on, 2007, 45(6): 1786-1797

[15] LiX, ChouteauM. Three-dimensional gravity modeling in all space [J]. Surveys in Geophysics, 1998, 19(4): 339-368

[16] HerrmannJ, MoghaddamP, StolkC C. Sparsity-and continuity-promoting seismic image recovery with curvelet frames [J]. Applied and Computational Harmonic Analysis, 2008, 24(2): 150-173

[17] MaJ, AntoniadisA, le DimetF X. Curvelet-based snake for multiscale detection and tracking of geophysical fluids [J]. Geoscience and Remote Sensing, IEEE Transactions on, 2006, 44(12): 3626-3638

[18] SimaD M, van HuffelS, GolubG H. Regularized total least squares based on quadratic eigenvalue problem solvers [J]. BIT Numerical Mathematics, 2004, 44(4): 793-812

[19] PlonkaG, MaJ-wei. Curvelet-wavelet regularized split bregman iteration for compressed sensing [J]. International Journal of Wavelets, Multiresolution and Information Processing, 2011, 9(1): 79-110

[20] LorisI, VerhoevenC. Iterative algorithms for total variation-like reconstructions in seismic tomography [J]. GEM-International Journal on Geomathematics, 2012, 3(2): 179-208

[21] RudinL I, OsherS, FatemiE. Nonlinear total variation based noise removal algorithms [J]. Physica D: Nonlinear Phenomena, 1992, 60(1): 259-268

[22] DonohoD L. Nonlinear solution of linear inverse problems by wavelet-vaguelette decomposition [J]. Applied and Computational Harmonic Analysis, 1995, 2(2): 101-126

[23] MaJ-wei. A single-pixel imaging system for remote sensing by two-step iterative curvelet thresholding [J]. Geoscience and Remote Sensing Letters, 2009, 6(4): 676-680

[24] PilkingtonM. 3-D magnetic imaging using conjugate gradients [J]. Geophysics, 1997, 62(4): 1132-1142

[25] KohK, KimS J, BoydS P. An interior-point method for large-scale L1-regularized logistic regression [J]. Journal of Machine Learning Research, 2007, 8(8): 1519-1555

[26] PortniaguineO, ZhdanovM S. 3-D magnetic inversion with data compression and image focusing [J]. Geophysics, 2002, 67(5): 1532-1541

[27] PilkingtonM. 3D magnetic data-space inversion with sparseness constraints [J]. Geophysics, 2008, 74(1): L7-L15

[28] ZHDANOV M, MUKHERJEE S. Modeling and inversion of 3-d gravity tensor field[EB/OL].[2013-11-19]. http://www.cemi.utah. edu/appls/grav/tensgrav/tensor.html.

[29] PortniaguineO, ZhdanovM. Focusing geophysical inversion images [J]. Geophysics, 1999, 64(3): 874-887

[30] FernandoD, ValériaB, JoãoS. Adaptive learning 3d gravity inversion for salt-body imaging [J]. Geophysics, 2011, 76(3): I49-I57

[31] KrahenbuhlABinary inversion of gravity data for salt imaging [D], 2005Colorado, USAColorado School of Mines