Gravity gradient tensor is the second derivatives of gravitational field. Compared with the traditional gravitational exploration, gravity gradient tensor can reach a better resolution. Gravity gradient tensor has 5 independent components, so they can contain more geological information. Therefore, gravity gradient tensor can be used to recover the causative bodies in the subsurface with high accuracy. Due to strong adaptability and flexibility of the unstructured grid, it can smoothly approximate the irregular and complex boundaries of anomalous source. Compared with the structured grid, the unstructured grid can provide researchers with more accurate discretization and calculation with less computation time. In order to reduce the ambiguity of the inversion, the authors chose to jointly and simultaneously invert all gravity gradient components based on tetrahedron grids with the help of the so called generalized objective function which is widely used in geophysical inverse problems. The authors applied the algorithm to each gravity tensor component at the beginning to ensure that the depth weighting function makes a difference in the inversion. It turns out that the depth weighting function works well. To explore all components of the gravity gradient tensor as much as possible, the authors described the joint-inversion in detail in this paper. The inversion results show that 3D inversion of gravity gradient tensor based on unstructured grids can obtain the position of causative bodies in the subsurface and the density distribution. Based on a comparison with rectangular grid inversion results, the authors highlighted the advantages of unstructured grid inversion. And the practicability of the method was verified through two synthetic models.
Wang C, Zhu Z Q, Lu G Y , et al. Forward simulation of analytic signals of gravity gradient tensor[J]. The Chinese Journal of Nonferrous Metals, 2013,23(9):2479-2497.
Shu Q, Zhang W Z, Zhou J X , et al. Error Study on gravity gradient tensor components transformation[J]. Geophysical and Geochemical Exploration, 2016,40(1):111-116.
Shu Q, Zhu X Y, Zhou J X , et al. Forward modeling expressions for right rectangular prism with constant density[J]. Geophysical and Geochemical Exploration, 2015,39(6):1217-1222.
[4]
董亮 . 非结构化网格生成技术研究及应用[D]. 镇江:江苏大学, 2010.
[4]
Dong L . The study and application of unstrucred mesh generation technique[D]. Zhenjiang:Jiangsu University, 2010.
Tang J T, Wang F Y . 2.5D direct current resistivity simulation based on the unstructured mesh[J]. Computing techniques for geophysical and geochemical exploration, 2008,30(5):413-418.
[7]
任政勇 . 基于非结构化网格的直流电阻率自适应有限元数值模拟[D]. 长沙:中南大学, 2007.
[7]
Ren Z Y . DC resistivity modeling by adaptive finite-element method with unstructured triangulation[D]. Changsha:Central South University, 2007.
[8]
孙丽影 . 起伏地形条件下大地电磁二维自适应有限元数值模拟[D]. 长沙:中南大学, 2013.
[8]
Sun L Y . Adaptive finite element numerical simulation of 2D MT with the topography[D]. Changsha:Central South University, 2013.
Liu C S . Three-dimension magnetotellurics adaptive edge finite-element numerical simulation based on unstructured meshes[D]. Changsha:Central South University, 2009.
Tong X Z . Research of forward using finite element method and regularized inversion using hybrid genetic algorithm in magnetocelluric sounding[D]. Changsha:Central South University, 2013.
[11]
曾思红 . 重力梯度张量正演研究及边界提取[D]. 长沙:中南大学, 2010.
[11]
Zeng S H . Gravity gradient tensor forward modeling and boundary extraction[D]. Changsha:Central South University, 2010.
Cao S J, Zhu Z Q, Lu G Y , et al. Forward modeling of full gravity gradient tensors based H-adaptive mesh refinement[J]. Progress in Geophysics, 2010,25(3):1015-1023.
[13]
Jahandari H, Farquharson C G . Forward modeling of gravity data using finite-volume and finite-element methods on unstructured grids[J]. Geophysics, 2012,78(3):G69-G80.
[14]
Okabe M . Analytical expressions for gravity anomalies due to homogeneous polyhedral bodies and translations into magnetic anomalies[J]. Geophysics, 1979,44(4):730-741.
[15]
Lelièvre P G, Farquharson C G, Hurich C A . Inversion of first-arrival seismic traveltimes without rays, implemented on unstructured grids[J]. Geophysical Journal International, 2011,185(2):749-763.
[16]
吉日嘎拉图 . 基于非结构化网格的重震联合反演[D]. 长沙:中南大学, 2014.
[16]
Ji R G L T . Joint inversion of gravity data and seismic traveltimes on unstructured grids[D]. Changsha:Central South University, 2014.
[17]
Okabe M . Finite element analysis of metal cutting mechanism[D]. Tokyo: University of Tokyo, 1974.
[18]
Guo L H, Meng X H, Shi L . 3D correlation imaging of the vertical gradient of gravity data[J]. Journal of Geophysics and Engineering, 2011(8):6-12.
[19]
Li Y G, Oldenburg D W . 3-D inversion of gravity data[J]. Geophysics, 1998,63(1):109-119.