1887
Volume 73, Issue 2
  • E-ISSN: 1365-2478

Abstract

ABSTRACT

The main method for calculating the gravity field involves discretizing the density sources into a stack of rectangular prisms with a regular grid distribution. The analytical formulation of the gravity anomaly for a right‐angled rectangular prism is affected by depth, with the kernel function decaying as depth increases. In addition, the efficiency of the computation and the storage requirements often pose challenges. We present a fast computational method for three‐dimensional gravity forward modelling of subsurface space using rectilinear grid and apply the Block–Toeplitz Toeplitz–Block method to the rectilinear grid. The size of the upright rectangles increases with depth to offset the effect of depth. We assume that the observation points are distributed on a homogeneous grid, and the kernel sensitivity matrices exhibit a Block–Toeplitz Toeplitz–Block structure, which is symmetric. For rectilinear dissections of subsurface space in MATLAB, the logarithmic interval size is used. The rectilinear mesh can offset the effect of depth to some degree allowing gravity anomalies to decrease more quickly. For the test of a single model, the gravity anomalies decrease faster and more rapidly in the case of the rectilinear grid compared to the uniform grid. In addition to this, we performed simulations on more complex models and demonstrated that using the Block–Toeplitz Toeplitz–Block method on this basis greatly improves the computational efficiency.

Loading

Article metrics loading...

/content/journals/10.1111/1365-2478.13658
2025-01-26
2025-12-15
Loading full text...

Full text loading...

References

  1. Barnett, C.T. (1976) Theoretical modeling of magnetic and gravitational‐fields of an arbitrarily shaped 3‐dimensional body. Geophysics, 41(6), 1353–1364
    [Google Scholar]
  2. Bhattacharyya, B.K. 1 & Navolio, M.E. 2 (1976) A fast Fourier transform method for rapid computation of gravity and magnetic anomalies due to arbitrary bodies. Geophysical Prospecting, 24(4), 633–649
    [Google Scholar]
  3. Boulanger, O. & Chouteau, M. (2001) Constraints in 3D gravity inversion.Geophysical Prospecting, 49(2), 265–280
    [Google Scholar]
  4. Cai, Y. & Wang, C. (2005) Fast finite‐element calculation of gravity anomaly in complex geological regions. Geophysical Journal International, 162(3), 696–708. Available from: https://doi.org/10.1111/j.1365‐246X.2005.02711.x
    [Google Scholar]
  5. Cao, S. J., Zhu, Z.‐J, Lu, G. Y., Zeng, S.‐G & Guo, W. B. (2010) Forward modelling of full gravity gradient tensors based H‐adaptive mesh refinement. Progress in Geophysics, 25(3), 1015–1023. https://doi.org/10.3969/j.issn.1004‐2903.2010.03.041
    [Google Scholar]
  6. Chen, L.W. & Liu, L.B. (2019) Fast and accurate forward modelling of gravity field using prismatic grids. Geophysical Journal International, 216(2), 1062–1071
    [Google Scholar]
  7. Chen, T. & Zhang, G.B. (2018) Forward modeling of gravity anomalies based on cell mergence and parallel computing. Computers & Geosciences, 120, 1–9
    [Google Scholar]
  8. Farquharson, C.G. & Mosher, C.R. W. (2009) Three‐dimensional modelling of gravity data using finite differences. Journal of Applied Geophysics, 68(3), 417–422
    [Google Scholar]
  9. Fedi, M. & Rapolla, A. (1999) 3‐D inversion of gravity and magnetic data with depth resolution. Geophysics, 64(2), 452–460
    [Google Scholar]
  10. Feng, RUI. (1986) A computation method of potential field with three‐dimensional density and magnetization distributions. Chinese Journal of Geophysics (in Chinese), 29(04), 399–406
    [Google Scholar]
  11. Hogue, J.D., Renaut, R.A. & Vatankhah, S. (2020) A tutorial and open source software for the efficient evaluation of gravity and magnetic kernels. Computers & Geosciences, 144, 104575
    [Google Scholar]
  12. Hu, S., Pan, K., Ren, Z., Tang, J., Wang, J. & Wang, P. (2023) Fast large‐scale gravity modeling using a high‐order compact cascadic multigrid strategy[J].Geophysics, 88(6), G125–G140
    [Google Scholar]
  13. Jahandari, H. & Farquharson, C.G. (2013) Forward modeling of gravity data using finite‐volume and finite‐element methods on unstructured grids. Geophysics, 78(3), G69–G80
    [Google Scholar]
  14. Jones, B.A., Beylkin, G., Born, G.H. & Provence, R.S. (2011) A multiresolution model for small‐body gravity estimation. Celestial Mechanics and Dynamical Astronomy, 111(3), 309–335
    [Google Scholar]
  15. Li, X. & Chouteau, M. (1998) Three‐dimensional gravity modeling in all space. Surveys in Geophysics, 19(4), 339–368
    [Google Scholar]
  16. Li, Y. & Oldenburg, D W. (2003) Fast inversion of large‐scale magnetic data using wavelet transforms and a logarithmic barrier method.Geophysical Journal International, 152(2), 251–265
    [Google Scholar]
  17. Li, Y. & Oldenburg, D.W. (1998) 3‐D inversion of gravity data. Geophysics, 63(1), 109–119
    [Google Scholar]
  18. Odegard, M.E. & Berg, J.W. (1965) Gravity interpretation using the Fourier integral. Geophysics, 30(3), 424
    [Google Scholar]
  19. Plouff, D. (1976) Gravity and magnetic fields of polygonal prisms and application to magnetic terrain corrections. Geophysics, 41(4), 727–741
    [Google Scholar]
  20. Ren, Z.Y., Tang, J.T., Kalscheuer, T. & Maurer, H. (2017) Fast 3‐D large‐scale gravity and magnetic modeling using unstructured grids and an adaptive multilevel fast multipole method. Journal of Geophysical Research: Solid Earth, 122(1), 79–109
    [Google Scholar]
  21. Renaut, R.A., Hogue, J.D. & Vatankhah, S., Liu, S. (2020) A fast methodology for large‐scale focusing inversion of gravity and magnetic data using the structured model matrix and the 2‐D fast Fourier transform. Geophysical Journal International, 223(2), 1378–1397
    [Google Scholar]
  22. Tontini, F.C., Cocchi, L. & Carmisciano, C. (2009) Rapid 3‐d forward model of potential fields with application to the Palinuro seamount magnetic anomaly (southern Tyrrhenian sea, Italy). Journal of Geophysical Research: Solid Earth, 114(2), B02103–B02119
    [Google Scholar]
  23. Vogel, C.R. (2002) Computational methods for inverse problem. Philadelphia, PA: SIAM.
    [Google Scholar]
  24. Wu, L.Y. (2016) Efficient modelling of gravity effects due to topographic masses using the Gauss‐FFT method. Geophysical Journal International, 205(1), 160–178
    [Google Scholar]
  25. Wu, L.Y. (2018) Efficient modeling of gravity fields caused by sources with arbitrary geometry and arbitrary density distribution. Surveys in Geophysics, 39(3), 401–434
    [Google Scholar]
  26. Wu, L.Y. & Tian, G. (2014) High‐precision Fourier forward modeling of potential fields. Geophysics, 79(5), G59–G68
    [Google Scholar]
  27. Yin, X., Yao, C., Zheng, Y., Xu, W., Chen, G. & Yuan, X. (2023) A fast 3D gravity forward algorithm based on circular convolution. Computers & Geosciences, 172, 105309
    [Google Scholar]
  28. Zhang, Y. & Wong, Y.S. (2015) BTTB‐based numerical schemes for three‐dimensional gravity field inversion. Geophysical Journal International, 203(1), 243–256.
    [Google Scholar]
  29. Zhang, Y.M., Wang, W.Y. & Xiong, S.Q. (2020) Research on the vertical recognition ability of gravity and magnetic data of point (line) source model with given survey accuracy. Chinese Journal of Geophysics (in Chinese), 63(11), 4220–4231, https://doi.org/10.6038/cjg2020N0412
    [Google Scholar]
/content/journals/10.1111/1365-2478.13658
Loading
/content/journals/10.1111/1365-2478.13658
Loading

Data & Media loading...

  • Article Type: Research Article
Keyword(s): gravity; modelling; numerical study; potential field

Most Cited This Month Most Cited RSS feed

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error