High-Efficiency Gravity Data Inversion Method Based on Locally Adaptive Unstructured Meshing

Abstract

In the 3-D density inversion calculation of gravity data, the entire subsurface space is discretized into rectangular prisms, and the density value of each prism is calculated. However, the method of dividing the entire space often results in invalid calculations. In this study, we proposed a highly efficient density inversion method using a locally adaptive unstructured mesh. In this method, the inversion scope is reduced by using the feature that the zero value of the tilt angle method corresponds to the edge of the field source. In addition, the local inversion area is divided by unstructured meshes, which can better represent irregular geological bodies and undulating terrain. In the proposed method, the size of the grid cells is changed according to the value of the tilt angle to reduce the amount of mesh, which reduces the computational complexity and improves the computational efficiency. In addition, by introducing a volume weighting function, the sensitivity of grid cells of different sizes can be balanced. Through synthetic modeling experiments, we verified that the locally adaptive unstructured mesh method can improve the efficiency and accuracy of inversion, flexibly deal with undulating terrain, and obtain the distribution features of irregular bodies. We applied this method to gravimetric data of the North Qinling region of the Shaanxi province, which clearly shows the spatial distribution information of five high-density geological bodies, and the depth range of ore bodies formed is from 240 to 1481 m according to our inversion results.