Optimization algorithm for rapid 3D gravity inversion*
Jing Lei, Yao Chang-Li, Yang Ya-Bin, Xu Meng-Long, Zhang Guang-Zhi, and Ji Ruo-Ye
1.School of geophysics and information technology, China University of Geosciences (Beijing), Beijing 100083, China.
2.The National Center for Geological Exploration Technology, Langfang 065000, China.
3.Institute of Geophysical and Geochemical Exploration, Chinese Academy of Geological Sciences, Langfang 065000,China.
Abstract:
The practical application of 3D inversion of gravity data requires a lot of computation time and storage space. To solve this problem, we present an integrated optimization algorithm with the following components: (1) targeting high accuracy in the space domain and fast computation in the wavenumber domain, we design a fast 3D forward algorithm with high precision; and (2) taking advantage of the symmetry of the inversion matrix, the main calculation in gravity conjugate gradient inversion is decomposed into two forward calculations, thus optimizing the computational efficiency of 3D gravity inversion We verify the calculation accuracy and efficiency of the optimization algorithm by testing various grid-number models through numerical simulation experiments.