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]] Hiramatsu Y, Sawada A, Kobayashi W, 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]] Beiki M, Pedersen L B. Eigenvector analysis of gravity gradient tensor to locate geologic bodies. Geophysics, 2010, 75(6): I37-I49, J]

[[3]] Li W-y, Liu Y-x, Li B, 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]] Mosher C R W, Farquharson C G. Minimum-structure borehole gravity inversion for mineral exploration: A synthetic modeling study. Geophysics, 2013, 78(2): G25-G39, J]

[[5]] Essa K S, Diab Z 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]] Fedi M, Ferranti L, Florio G, 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]] Elhussein M, Abo-Ezz E R, Gomaa O 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]] Essa K S, Géraud Y. 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]] Biswas A. 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, 2(1): 30, J]

[[10]] Essa K S, Abo-Ezz E R. Potential field data interpretation to detect the parameters of buried geometries by applying a nonlinear least-squares approach. Acta Geodaetica et Geophysica, 2021, 56(2): 387-406, J]

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

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

[[13]] Essa K S, Abo-Ezz E R, Géraud Y, 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, 215: 110683, J]

[[14]] Lu G-y, Cao S-j, Zhu Z-qiang. Joint inversion of gravity and multiple components of tensor gravity data. Journal of Central South University, 2016, 23(7): 1767-1777, J]

[[15]] Roy A, Wu L-yuan. Generalized Gauss-FFT 3D forward gravity modeling for irregular topographic mass having any 3D variable density contrast. Computers and Geosciences, 2023, 172: 105297, J]

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

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

[[18]] Johnson L R, Litehiser J 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]] Chakravarthi V, Raghuram H M, Singh S 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]] Ren Z-y, Chen C-j, Zhong Y-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]] Wu L-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]] Zhao G-d, Chen B, Chen L-w, et al.. High-accuracy 3D Fourier forward modeling of gravity field based on the Gauss-FFT technique. Journal of Applied Geophysics, 2018, 150: 294-303, J]

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

[[24]] Pan K-j, Zhang Z-h, Hu S-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]] Zhang J, Wang C Y, Shi Y-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]] Farquharson C G, Mosher C R W. Three-dimensional modelling of gravity data using finite differences. Journal of Applied Geophysics, 2009, 68(3): 417-422, J]

[[27]] Jahandari H, Bihlo A, Donzelli F. Forward modelling of gravity data on unstructured grids using an adaptive mimetic finite-difference method. Journal of Applied Geophysics, 2021, 190: 104340, J]

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

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

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

[[31]] Gharti H N, Tromp J, Zampini S. Spectral-infinite-element simulations of gravity anomalies. Geophysical Journal International, 2018, 215(2): 1098-1117, J]

[[32]] Martin R, Chevrot S, Komatitsch D, 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]] Long J-b, Farquharson C 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]] Dai S-k, Chen Q-r, Li K, 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]] Dai S-k, Zhao D-d, Wang S-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]] Liu J-x, Liu P-m, Tong X-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]] Chen J-p, Haber E, Oldenburg D W. Three-dimensional numerical modelling and inversion of magnetometric resistivity data. Geophysical Journal International, 2002, 149(3): 679-697, J]

[[38]] Pan K-j, Wang J-x, Hu S-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]] Ren Z-y, Chen C-j, Pan K-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]] Dubey C P, Tiwari V M. Computation of the gravity field and its gradient. Computers & Geosciences, 2016, 88(C): 83-96, J]

[[41]] Wang X-l, Liu J-x, Li J, et al.. Fast 3D gravity and magnetic modelling using midpoint quadrature and 2D FFT. Scientific Reports, 2023, 13: 9304, J]