3D gravity inversion based on Cauchy distribution constraint and fast proximal objective function optimization

Peng GuoMin,Liu Zhan,Yu HuiZhen,Xu KaiJun,Chen XiaoHong

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

2018-01-01

Abstract:The choice of an optimization algorithm affects computational efficiency of 3D gravity data inversion to a great extent and thereby constrains its application in the industry. In the context of complex subsurface geological settings, sharp boundaries may exist between different rock units that have large physical property contrasts. To this end, we propose a novel 3D gravity inversion method based on Cauchy distribution constraint and fast proximal objective function (FPOF) optimization. A salient characteristic of FPOF optimization is to calculate the unknown density in each subdivided unit of the inversion domain one by one at each iteration, thus reducing computational complexity and improving computational efficiency. In addition, the inclusion of Cauchy norm into objective function imposes sparseness on model parameters to be inverted and consequently facilitates a blocky density distribution. The test of synthetic data demonstrates that the proposed inversion method not only recovers more focused anomalous density bodies but also costs less run time than conventional CG optimization method. At last, the application of the proposed inversion method to real gravity data in the Western China shows that the recovered 3D density model agrees well with the known geological knowledge.