Joint Nonlinear Inversion of Full Tensor Gravity Gradiometry Data and Its Parallel Algorithm

Abstract

Geophysical joint inversion is more frequently applied to deep crust probes. For the potential field data, it requires the introduction of large-scale observed data and sensitivity matrices. Massive matrix-vector multiplications occur during iterations, and the obtained data cannot be effectively interpreted by merely using desktop computers. To improve the resolution and the computing ability of inversion, we here propose the parallel joint nonlinear inversion of full tensor gravity gradiometry data. As the inversion is affected by linear searches, it is associated with certain classical computing performance issues. Hence, we addressed the memory and efficiency limitations, which are caused by very large calculation volumes. Also, we performed a quantitative and comprehensive feasibility analysis of parallel computing. Then, we identified the main factor influencing the inversion performance and clarified the correspondence between the cell number and the memory. A parallel inversion solution was proposed via graphics processing unit (GPU) based on the sensitivity matrix compression. The data tests demonstrated that inversion has antinoise property and that it can obtain accurate underground density distributions. Also, the parallel solution was found to be suitable for inverting cells at the million cell scale and greater because of its ability of acceleration and matrix compression. A design pattern was applied for gravity or magnetic anomaly inversion of 100 <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$\times 100\times20$ </tex-math></inline-formula> cells, and the run time was less than 1 min. Overall, we believe that the proposed solution can help implement massive potential field data inversions and promote the application of the parallel technique in other inversion research.