High-accuracy gravity field and gravity gradient forward modelling based on 3D vertex-centered finite-element algorithm

Abstract

Gravity anomalies generated by density non-uniformity are governed by the 3D Poisson equation. Most existing forward methods for such anomalies rely on integral techniques and cell-centered numerical approaches. Once the gravitational potential is calculated, these numerical schemes will inevitably lose high accuracy. To alleviate this problem, an accurate and efficient high-order vertex-centered finite-element scheme for simulating 3D gravity anomalies is presented. Firstly, the forward algorithm is formulated through the vertex-centered finite element with hexahedral meshes. The biconjugate gradient stabilized algorithm can solve the linear equation system combined with an incomplete LU factorization (ILU-BICSSTAB). Secondly, a high-degree Lagrange interpolating scheme is applied to achieve the first-derivate and second-derivate gravitational potential. Finally, a 3D cubic density model is used to test the accuracy of the vertex-centered finite-element algorithm, and thin vertical rectangular prisms and real example for flexibility. All numerical results indicate that our high-order vertex-centered finite-element method can provide an accurate approximation for the gravity field vector and the gravity gradient tensor. Meanwhile, compared to the cell-centered numerical algorithm, the high-order vertex-centered finite element scheme exhibits higher efficiency and accuracy in simulating 3D gravity anomalies.

Keywords

gravity anomalies / 3D Poisson equation / vertex-centered finite-element algorithm / numerical simulation / ILU-BICGSTAB iterative method

References

[1] HiramatsuY, SawadaA, KobayashiW, et al. . Gravity gradient tensor analysis to an active fault: A case study at the Togi-gawa Nangan fault, Noto Peninsula, central Japan. Earth, Planets and Space, 2019, 71(1): 107 J]

[2] BeikiM, PedersenL B. Eigenvector analysis of gravity gradient tensor to locate geologic bodies. Geophysics, 2010, 756I37-I49 J]

[3] LiW-y, LiuY-x, LiB, et al. . Hydrocarbon exploration in the South Yellow Sea based on airborne gravity, China. Journal of Earth Science, 2016, 27(4): 686-698 J]

[4] MosherC R W, FarquharsonC G. Minimum-structure borehole gravity inversion for mineral exploration: A synthetic modeling study. Geophysics, 2013, 78(2): G25-G39 J]

[5] EssaK S, DiabZ E. Source parameters estimation from gravity data using Bat algorithm with application to geothermal and volcanic activity studies. International Journal of Environmental Science and Technology, 2023, 20(4): 4167-4187 J]

[6] FediM, FerrantiL, FlorioG, et al. . Understanding the structural setting in the Southern Apennines (Italy): Insight from gravity gradient tensor. Tectonophysics, 2005, 397(1–2): 21-36 J]

[7] ElhusseinM, Abo-EzzE R, GomaaO A, et al. . Two-sided fault parameter assessment through the utilization of a particle swarm optimizer to different gravity horizontal gradients-orders with application to various cases studies. Environmental Earth Sciences, 2023, 82(21): 502 J]

[8] EssaK S, GéraudY. Parameters estimation from the gravity anomaly caused by the two-dimensional horizontal thin sheet applying the global particle swarm algorithm. Journal of Petroleum Science and Engineering, 2020, 193: 107421 J]

[9] BiswasA. Interpretation of gravity and magnetic anomaly over thin sheet-type structure using very fast simulated annealing global optimization technique. Modeling Earth Systems and Environment, 2016, 2130 J]

[10] EssaK S, Abo-EzzE R. Potential field data interpretation to detect the parameters of buried geometries by applying a nonlinear least-squares approach. Acta Geodaetica et Geophysica, 2021, 562387-406 J]

[11] PedersenL B, BastaniM. Estimating rock-vector magnetization from coincident magnetic field measurements and gravity gradient tensor. Geophysics, 2016, 81(3): 55-64 J]

[12] RimH, LiY-guo. Single-hole imaging using borehole gravity gradiometry. Geophysics, 2012, 77(5): G67-G76 J]

[13] EssaK S, Abo-EzzE R, GéraudY, et al. . A full interpretation applying a metaheuristic particle swarm for gravity data of an active mud diapir, SW Taiwan. Journal of Petroleum Science and Engineering, 2022, 215110683 J]

[14] LuG-y, CaoS-j, ZhuZ-qiang. Joint inversion of gravity and multiple components of tensor gravity data. Journal of Central South University, 2016, 23(7): 1767-1777 J]

[15] RoyA, WuL-yuan. Generalized Gauss-FFT 3D forward gravity modeling for irregular topographic mass having any 3D variable density contrast. Computers and Geosciences, 2023, 172105297 J]

[16] WonI J, BevisM. Computing the gravitational and magnetic anomalies due to a polygon: Algorithms and Fortran subroutines. Geophysics, 1987, 52(2): 232-238 J]

[17] BanerjeeB, DasG S P. Gravitational attraction of a rectangular parallelepiped. Geophysics, 1977, 42(5): 1053-1055 J]

[18] JohnsonL R, LitehiserJ J. A method for computing the gravitational attraction of three-dimensional bodies in a spherical or ellipsoidal Earth. Journal of Geophysical Research, 1972, 77(35): 6999-7009 J]

[19] ChakravarthiV, RaghuramH M, SinghS B. 3-D forward gravity modeling of basement interfaces above which the density contrast varies continuously with depth. Computers & Geosciences, 2002, 28(1): 53-57 J]

[20] RenZ-y, ChenC-j, ZhongY-y, et al. . Recursive analytical formulae of gravitational fields and gradient tensors for polyhedral bodies with polynomial density contrasts of arbitrary non-negative integer orders. Surveys in Geophysics, 2020, 41(4): 695-722 J]

[21] WuL-yuan. Efficient modelling of gravity effects due to topographic masses using the Gauss-Fft method. Geophysical Journal International, 2016, 205(1): 160-178 J]

[22] ZhaoG-d, ChenB, ChenL-w, et al. . High-accuracy 3D Fourier forward modeling of gravity field based on the Gauss-FFT technique. Journal of Applied Geophysics, 2018, 150294-303 J]

[23] WuL-y, TianGang. High-precision Fourier forward modeling of potential fields. Geophysics, 2014, 79(5): G59-G68 J]

[24] PanK-j, ZhangZ-h, HuS-g, et al. . Three-dimensional forward modelling of gravity field vector and its gradient tensor using the compact difference schemes. Geophysical Journal International, 2021, 224(2): 1272-1286 J]

[25] ZhangJ, WangC Y, ShiY-l, et al. . Three-dimensional crustal structure in central Taiwan from gravity inversion with a parallel genetic algorithm. Geophysics, 2004, 69(4): 917-924 J]

[26] FarquharsonC G, MosherC R W. Three-dimensional modelling of gravity data using finite differences. Journal of Applied Geophysics, 2009, 68(3): 417-422 J]

[27] JahandariH, BihloA, DonzelliF. Forward modelling of gravity data on unstructured grids using an adaptive mimetic finite-difference method. Journal of Applied Geophysics, 2021, 190104340 J]

[28] CaiY-e, WangC Y. Fast finite-element calculation of gravity anomaly in complex geological regions. Geophysical Journal International, 2005, 162(3): 696-708 J]

[29] MayD A, KnepleyM G. Optimal, scalable forward models for computing gravity anomalies. Geophysical Journal International, 2011, 187(1): 161-177 J]

[30] JahandariH, FarquharsonC G. Forward modeling of gravity data using finite-volume and finite-element methods on unstructured grids. Geophysics, 2013, 78(3): G69-G80 J]

[31] GhartiH N, TrompJ, ZampiniS. Spectral-infinite-element simulations of gravity anomalies. Geophysical Journal International, 2018, 215(2): 1098-1117 J]

[32] MartinR, ChevrotS, KomatitschD, et al. . A highorder 3-D spectral-element method for the forward modelling and inversion of gravimetric data—Application to the western Pyrenees. Geophysical Journal International, 2017, 209(1): 406-424[J]

[33] LongJ-b, FarquharsonC G. Three-dimensional forward modelling of gravity data using mesh-free methods with radial basis functions and unstructured nodes. Geophysical Journal International, 2019, 217(3): 1577-1601[J]

[34] DaiS-k, ChenQ-r, LiK, et al. . The forward modeling of 3D gravity and magnetic potential fields in space-wavenumber domains based on an integral method. Geophysics, 2022, 87(3): G83-G96 J]

[35] DaiS-k, ZhaoD-d, WangS-g, et al. . Three-dimensional numerical modeling of gravity and magnetic anomaly in a mixed space-wavenumber domain. Geophysics, 2019, 84(4): G41-G54 J]

[36] LiuJ-x, LiuP-m, TongX-zhong. Three-dimensional land FD-CSEM forward modeling using edge finite-element method. Journal of Central South University, 2018, 25(1): 131-140 J]

[37] ChenJ-p, HaberE, OldenburgD W. Three-dimensional numerical modelling and inversion of magnetometric resistivity data. Geophysical Journal International, 2002, 149(3): 679-697 J]

[38] PanK-j, WangJ-x, HuS-g, et al. . An efficient cascadic multigrid solver for 3-D magnetotelluric forward modelling problems using potentials. Geophysical Journal International, 2022, 230(3): 1834-1851 J]

[39] RenZ-y, ChenC-j, PanK-j, et al. . Gravity anomalies of arbitrary 3D polyhedral bodies with horizontal and vertical mass contrasts. Surveys in Geophysics, 2017, 38(2): 479-502 J]

[40] DubeyC P, TiwariV M. Computation of the gravity field and its gradient. Computers & Geosciences, 2016, 88(C): 83-96 J]

[41] WangX-l, LiuJ-x, LiJ, et al. . Fast 3D gravity and magnetic modelling using midpoint quadrature and 2D FFT. Scientific Reports, 2023, 139304 J]