3D joint inversion of gravity and magnetic data using Cross gradient constraint

Authors

1 Faculty of Mining, Petroleum and Geophysics, Shahrood University of Technology, Shahrood, Iran

2 Associate Professor, Faculty of Mining, Petroleum and Geophysics, Shahrood University of Technology, Shahrood, Iran.

3 Malayer university, Malayer, Iran

Abstract

Potential field methods, specifically gravity and magnetic method are fundamental geophysical techniques in mineral exploration. The interpretation of this data heavily relies on inversion modeling.
Inversion of potential field data to estimate subsurface model parameters is a numerically challenging and inherently non-unique problem. An infinite number of models can adequately satisfy the observed data. A crucial approach to mitigate this non-uniqueness is incorporating geological constraints relevant to the survey area.
While individual inversions of geophysical datasets often provide satisfactory reconstructions of subsurface anomaly locations, integrating information from multiple geophysical methods is essential for achieving a better recovery and subsequently a more robust interpretation of subsurface structures. In many practical scenarios, single-method geophysical surveys fail to provide unambiguous interpretations. Thus, the combination of geophysical methods is imperative. Joint inversion, whether for identical or different physical parameters, is an effective technique to achieve enhanced resolution. Consequently, a comprehensive geophysical interpretation based on joint inversion has gained widespread application in recent years.
This paper investigates the three-dimensional (3D) joint inversion of gravity and magnetic data by utilizing a cross-gradient function to enforce structural similarity. The performance of the joint inversion algorithm is first validated using synthetic model data. Finally, the algorithm is applied to real field data to demonstrate its practical efficacy

Keywords


Bennington, N. L., Zhang, H., Thurber, C. H., & Bedrosian, P. A. (2015). Joint inversion of seismic and magnetotelluric data in the Parkfield Region of California using the normalized cross-gradient constraint. Pure and Applied Geophysics, 172, 1033-1052.
Candansayar, M. E., & Tezkan, B. (2008). Two‐dimensional joint inversion of radiomagnetotelluric and direct current resistivity data. Geophysical Prospecting, 56(5), 737-749.
Demirci, İ., Dikmen, Ü., & Candansayar, M. E. (2018). Two-dimensional joint inversion of Magnetotelluric and local earthquake data: Discussion on the contribution to the solution of deep subsurface structures. Physics of the Earth and Planetary Interiors, 275, 56-68.
Gallardo, L. A., & Meju, M. A. (2003). Characterization of heterogeneous near‐surface materials by joint 2D inversion of dc resistivity and seismic data. Geophysical Research Letters, 30(13).
Gallardo, L. A., & Meju, M. A. (2004). Joint two‐dimensional DC resistivity and seismic travel time inversion with cross‐gradients constraints. Journal of Geophysical Research: Solid Earth, 109(B3).
Gallardo, L. A., & Meju, M. A. (2007). Joint two-dimensional cross-gradient imaging of magnetotelluric and seismic traveltime data for structural and lithological classification. Geophysical Journal International, 169(3), 1261-1272.
Gallardo, L. A., & Meju, M. A. (2011). Structure‐coupled multiphysics imaging in geophysical sciences. Reviews of Geophysics, 49(1).
Gross, L. (2019). Weighted cross-gradient function for joint inversion with the application to regional 3-D gravity and magnetic anomalies. Geophysical Journal International, 217(3), 2035-2046.
Hansen, P. C. (2010). Discrete inverse problems: insight and algorithms. Society for Industrial and Applied Mathematics.
Hu, W., Abubakar, A., & Habashy, T. M. (2009). Joint electromagnetic and seismic inversion using structural constraints. Geophysics, 74(6), R99-R109.
Joulidehsar, F., Moradzadeh, A., & Doulati Ardejani, F. (2018). An improved 3D joint inversion method of potential field data using cross-gradient constraint and LSQR method. Pure and Applied Geophysics, 175, 4389-4409.
Menichetti, V., & Guillen, A. (1983). Simultaneous interactive magnetic and gravity inversion. Geophysical Prospecting, 31(6), 929-944.
Molodtsov, D. M., Troyan, V. N., Roslov, Y. V., & Zerilli, A. (2013). Joint inversion of seismic traveltimes and magnetotelluric data with a directed structural constraint. Geophysical Prospecting, 61(6), 1218-1228.
Moorkamp, M., Heincke, B., Jegen, M., Roberts, A. W., & Hobbs, R. W. (2011). A framework for 3-D joint inversion of MT, gravity and seismic refraction data. Geophysical Journal International, 184(1), 477-493.
Oldenburg, D. W., & Li, Y. (2005). Inversion for applied geophysics: A tutorial. Near-surface geophysics, 89-150.
Oliver-Ocaño, F. M., Gallardo, L. A., Romo-Jones, J. M., & Pérez-Flores, M. A. (2019). Structure of the Cerro Prieto Pull-apart basin from joint inversion of gravity, magnetic and magnetotelluric data. Journal of Applied Geophysics, 170, 103835.
Saunders, J. H., Herwanger, J. V., Pain, C. C., Worthington, M. H., & De Oliveira, C. R. E. (2005). Constrained resistivity inversion using seismic data. Geophysical Journal International, 160(3), 785-796.
Shi, B., Yu, P., Zhao, C., Zhang, L., & Yang, H. (2018). Linear correlation constrained joint inversion using squared cosine similarity of regional residual model vectors. Geophysical Journal International, 215(2), 1291-1307.
Tavakoli, M., Kalateh, A. N., Rezaie, M., Gross, L., & Fedi, M. (2021). Sequential joint inversion of gravity and magnetic data via the cross‐gradient constraint. Geophysical Prospecting, 69(7), 1542-1559.
Vatankhah, S., Renaut, R. A., Huang, X., Mickus, K., & Gharloghi, M. (2022). Large-scale focusing joint inversion of gravity and magnetic data with Gramian constraint. Geophysical Journal International, 230(3), 1585-1611.
Vernant, P., Masson, F., Bayer, R., & Paul, A. (2002). Sequential inversion of local earthquake traveltimes and gravity anomaly—The example of the western Alps. Geophysical Journal International, 150(1), 79-90.
Wang, K. P., Tan, H. D., & Wang, T. (2017). 2D joint inversion of CSAMT and magnetic data based on cross-gradient theory. Applied Geophysics, 14, 279-290.
Wang, X., Fang, J., & Hsu, H. (2011). Three-dimensional crustal density distribution beneath North China by sequential inversion of local earthquake traveltimes and gravity anomaly. Earthquake Science, 24(2), 135-141.
Zhang, Y., & Wang, Y. (2020). Three-dimensional gravity-magnetic cross-gradient joint inversion based on structural coupling and a fast gradient method. Journal of Computational Mathematics, 758-777.
Zhou, J., Xiu, C., & Zhang, X. (2015). Simultaneous structure-coupled joint inversion of gravity and magnetic data based on a damped least-squares technique. International Journal of Geosciences, 6(02), 172.