Optimization algorithm for rapid 3D gravity inversion

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.