Three-dimensional joint inversion of gravity and magnetic data using Gramian constraints and Gauss-Newton method

Kong RuiJin,Hu XiangYun,Cai HongZhu

DOI: https://doi.org/10.6038/cjg2022Q0394

2023-01-01

Abstract:This paper proposes an approach for the joint inversion of gravity and magnetic data using Gramian constraints to address the non-uniqueness of the inverse problem. The Gramian constraints can be constructed either directly from the model or from the spatial gradient of the model. By utilizing the Gramian operator, the correlation between the anomalous density and magnetic susceptibility can be enhanced. The Gauss-Newton method is employed to minimize the objective functional, and the equation for model update is simplified to reduce computational costs. To ensure the inversion results align well with geological models, we restrict the boundary of the model parameters using the logarithmic method. Additionally, we introduce a novel step-length calculation strategy specifically tailored for the joint inversion approach. To validate our algorithm, we conducted experiments using three synthetic models. Furthermore, we applied the method to airborne gravity and magnetic surveys in the McFaulds Lake area, Ontario, Canada. The inversion results and cross-plots demonstrate a stronger relationship between the physical property parameters or gradient distributions of anomalous density and susceptibility obtained through joint inversion compared to those obtained through separate inversions. The joint inversion algorithm effectively produces inversion results that accurately fit the data and satisfy the coupling assumption. Consequently, we successfully validate the effectiveness and applicability of the developed joint inversion algorithm using practical data.