The Efficient 3D Gravity Focusing Density Inversion Based on Preconditioned JFNK Method under Undulating Terrain: A Case Study from Huayangchuan, Shaanxi Province, China

Since polymetallic ores show higher anomalies in gravity exploration methods, we usually obtain the position and range of ore bodies by density inversion of gravity data. The three-dimensional (3D) gravity focusing density inversion is a common interpretation method in mineral exploration, which can directly and quantitatively obtain the density distribution of subsurface targets. However, in actual cases, it is computation inefficient. We proposed the preconditioned Jacobian-free Newton-Krylov (JFNK) method to accomplish the focusing inversion. The JFNK method is an efficient algorithm in solving large sparse systems of nonlinear equations, and we further accelerate the inversion process by the preconditioned technique. In the actual area, the gravity anomalies are distributed on the naturally undulating surface. Nowadays, the gravity inversion under undulating terrain was mainly achieved by discretizing the ground into unstructured meshes, but it is complicated and time-consuming. To improve the practicality, we presented an equivalent-dimensional method that incorporates unstructured meshes with structured meshes in gravity inversion, and the horizontal size is determined by the gradient of observed gravity and terrain data. The small size meshes are adopted at the position where the terrain or gravity gradient is large. We used synthetic data with undulating-terrain to test our new method. The results indicated that the recovered model obtained by this method was similar to the inversion method of unstructured meshes, and the new method computes faster. We also applied the method to field data in Huayangchuan, Shaanxi Province. The survey area has complicated terrain conditions and contains multiple polymetallic ores. Based on the high-density characteristics of polymetallic ore bodies in the area, we calculate the field data into 3D density models of the subsurface by the preconditioned JFNK method and infer six polymetallic ores.