Fast 3D joint inversion of gravity and magnetic data based on cross gradient constraint
), Abstract
The gravity and magnetic data can be adopted to interpret the internal structure of the Earth. To improve the calculation efficiency during the inversion process and the accuracy and reliability of the reconstructed physical property models, the triple strategy is adopted in this paper to develop a fast cross-gradient joint inversion for gravity and magnetic data. The cross-gradient constraint contains solving the gradients of the physical property models and performing the cross-product calculation of their gradients. The sparse matrices are first obtained by calculating the gradients of the physical property models derived from the first-order finite difference. Then, the triple method is applied to optimize the storages and the calculations related to the gradients of the physical property models. Therefore, the storage compression amount of the calculations related to the gradients of the physical property models and the cross-gradient constraint are reduced to one-fold of the number of grid cells at least, and the compression ratio increases with the increase of the number of grid cells. The test results from the synthetic data and field data prove that the structural coupling is achieved by using the fast cross-gradient joint inversion method to effectively reduce the multiplicity of solutions and improve the computing efficiency.
References
N. Linde, A. Tryggvason, J.E. Peterson, et al., Joint inversion of crosshole radar and seismic traveltimes acquired at the south oyster bacterial transport site, Geophysics 73 (4) (2008) G29–G37.
E. Fregoso, L.A. Gallardo, Cross-gradients joint 3D inversion with applications to gravity and magnetic data, Geophysics 74 (4) (2009) L31–L42.
M. Moorkamp, B. Heincke, M. Jegen, A framework for 3D joint inversion of MT, gravity and seismic refraction data, Geophys J Int 184 (1) (2011a) 477–493.
M. Moorkamp, A.W. Roberts, M. Jegen, et al., Verification of velocity-resistivity relationships derived from structural joint inversion with borehole data, Geophys Res Lett 40 (14) (2013) 3596–3601.
G.H.F. Gardner, L.W. Gardner, A.R. Gregory, Formation velocity and density-the diagnostic basics for stratigraphic traps, Geophysics 39 (6) (1974) 770–780.
M. Sun, S.G. Jin, Multiparameter elastic full waveform inversion of ocean bottom seismic four-component data based on a modified acoustic-elastic coupled equation, Rem Sens 12 (17) (2020) 2816.
J. Sun, Y. Li, Multidomain petrophysically constrained inversion and geology differentiation using guided fuzzy c-means clustering, Geophysics 80 (4) (2015) ID1–ID18.
J. Sun, Y. Li, Joint inversion of multiple geophysical data using guided fuzzy c-means clustering, Geophysics 81 (3) (2016) ID37–ID57.
R. Tenzer, W. Chen, S.G. Jin, Effect of the upper mantle density structure on the Moho geometry, Pure Appl Geophys 172 (6) (2015) 1563–1583.
B. Zhang, S.G. Jin, C. Liu, Z. Guo, X. Liu, Prediction of shear wave velocity based on a statistical rock-physics model and Bayesian theory, J Petrol Sci Eng 195 (2020), 107710.
B. Heincke, M. Jegen, R. Hobbs, Joint inversion of MT, gravity and seismic data applied to sub-basalt imaging, SEG Technical Program Expanded Abstracts, 2006.
D. Colombo, M.D. Stefano, Geophysical modeling via simultaneous joint inversion of seismic, gravity, and electromagnetic data: application to prestack depth imaging, Lead Edge 26 (3) (2007) 326–331.
L.A. Gallardo, M.A. Meju, Characterization of heterogeneous near-surface materials by joint 2D inversion of dc resistivity and seismic data, Geophys Res Lett 30 (13) (2003).
L.A. Gallardo, M.A. Meju, Joint two-dimensional DC resistivity and seismic traveltime inversion with cross-gradients constraints, J Geophys Res 109 (B3) (2004) 3311–3315.
L.A. Gallardo, Multiple cross-gradient joint inversion for geospectral imaging, Geophys Res Lett 34 (19) (2007), L19301.
L.A. Gallardo, S.L. Fontes, M.A. Meju, M.P. Buonora, P.P. De Lugao, Robust geophysical integration through structure-coupled joint inversion and multispectral fusion of seismic reflection, magnetotelluric, magnetic, and gravity images: example from Santos Basin, offshore Brazil, Geophysics. 77 (2012) B237–B251.
T.L. Li, R.Z. Zhang, Y.C. Pak, et al., Multiple joint inversion of geophysical data with sub-region cross-gradient constraints, Chin J Geophys 59 (8) (2016) 2979–2988 (in Chinese).
J. Gao, H.J. Zhang, Two-dimensional joint inversion of seismic velocity and electrical resistivity using seismic travel times and full channel electrical measurements based on alternating cross-gradient structure constraint, Chin J Geophys 59 (11) (2016) 4310–4322 (in Chinese).
N.L. Bennington, H. Zhang, C.H. Thurber, Joint inversion of seismic and magnetotelluric data in the park field region of California using the normalized cross-gradient constraint, Pure Appl Geophys 172 (5) (2015) 1033–1052.
R.Z. Zhang, T.L. Li, 2D data-space joint inversion of MT, gravity, magnetic, and seismic data with cross-gradient constraints, Geophys Prospect 68 (2020) 721–731.
V. Saeed, S. Liu, A.R. Rosemary, X.Y. Hu, et al., An efficient alternating algorithm for the Lp-norm cross-gradient joint inversion of gravity and magnetic data using the 2-D fast fourier transform, IEEE Trans Geosci Rem Sens 60 (2022), 4500416.
Y.G. Li, D.W. Oldenburg, 3-D inversion of magnetic data, Geophysics 61 (2) (1996) 394–408.
Y.G. Li, D.W. Oldenburg, Joint inversion of surface and three-component borehole magnetic data, Geophysics 65 (2) (2000) 540–552.
S. Liu, S. Jin, S. Xuan, X. Liu, Three-dimensional data-space joint inversion of gravity and magnetic data with correlation-analysis constrains, Ann Geophys 65 (2022), https://doi.org/10.4401/ag-8570.
R.Z. Zhang, T.L. Li, Joint MT and gravity inversion using structural constraints: a case study from the Linjiang copper mining area, Jilin, China, Minerals 9 (2019) 407.
京公网安备11010802044758号