Forward modeling for 2D body gravity anomaly in the space-wavenumber mixed domain

DAI Shikun1,2 WANG Xulong1,2 ZHAO Dongdong1,2 LIU Zhiwei3 ZHANG Qianjiang4 SUN Jinfei5

(1.Key Laboratory of Metallogenic Prediction of Nonferrous Metals and Geological Environment Monitoring, Ministry of Education, Changsha, Hunan Province, China 410083)
(2.School of Geosciences and Info-physics, Central South University, Changsha, Hunan Province, China 410083)
(3.Guangdong Electric Power Design Institute, China Energy Engineering Group Co., Ltd., Guangzhou, Guangdong Province, China 510663)
(4.College of Earth Sciences, Guilin University of Technology, Guilin, Guangxi Province, China 541004)
(5.No. 1 Daqing Geophysical Exploration Company of BGP, CNPC, Daqing, Heilongjiang Province, China 163357)
【Knowledge Link】gravity anomaly

【Abstract】An efficient and high-precision forward modeling of 2D body gravity anomaly in a space-wavenumber mixed domain is proposed to calculate the gravity anomaly in the case of arbitrary density distribution. Processed by 1D Fourier transform, 2D partial differential equation governing gravity potential in the spatial domain is transformed into a group of independent 1D differential equations engaged with different wavenumbers. The method preserves the vertical direction in the spatial domain, so shallow and deep meshes for modeling can be fine and coarse respectively. By this way, both the accuracy and the efficiency of calculation are taken into account. The 1D finite element method is used to solve the transformed differential equations with different wavenumbers. Moreover, the efficiency of solving linear equations with a fixed bandwidth assembled by the finite element analysis is further improved by a chasing method. To test the proposed algorithm, a constant-density 2D model and a variable-density 2D model are designed. Numerical test results show that the proposed algorithm not only guarantees the calculation accuracy but improves the calculation efficiency compared with the conventional methods.

【Keywords】 gravity anomaly; space-wavenumber mixed domain; 2D body; forward modeling;


【Funds】 National Key Research and Development Program (2018YFC0603602) National Natural Science Foundation of China (41574127) Self-innovative Projects for Postgraduates in Central South University (2018ZZTS710) Natural Science Foundation of Guangxi Province, China (2018GXSFAA050070) Special Fund for Scientific and Technological Base and Talents in Guangxi Province, China (AD19110058)

Download this article


    [1] Hubbert M K. A line-integral method of computing the gravimetric effects of two-dimensional masses [J]. Geophysics, 1948, 13 (2): 215–225.

    [2] Talwani M, Worzel J L, Landisman M. Rapid gravity computations for two-dimensional bodies with application to the Mendocino submarine fracture zone [J]. Journal of Geophysical Research, 1959, 64 (1): 49–59.

    [3] Rao D B. Modelling of sedimentary basins from gravity anomalies with variable density contrast [J]. Geophysical Journal International, 1986, 84 (1): 207–212.

    [4] Murthy I V R, Rao D B. Gravity anomalies of two-dimensional bodies of irregular cross-section with density contrast varying with depth [J]. Geophysics, 1979, 44 (9): 1525–1530.

    [5] Pan J. Gravity anomalies of irregularly shaped two-dimensional bodies with constant horizontal density gradient [J]. Geophysics, 1989, 54 (4): 528–530.

    [6] Kwok Y K. Contour integrals for gravity computation of horizontal bodies with variable density [J]. Applied Mathematical Modelling, 1991, 15 (2): 98–103.

    [7] Zhou X B. Analytic solution of the gravity anomaly of irregular 2D masses with density contrast varying as a 2D polynomial function [J]. Geophysics, 2010, 75 (2), I11–I19.

    [8] D’Urso M G. The gravity anomaly of a 2D polygonal body having density contrast given by polynomial functions [J]. Surveys in Geophysics, 2015, 36 (3): 391–425.

    [9] LI Ming, WANG Yichang. The computation and application of vertical first derivative of gravity [J]. Oil Geophysical Prospecting, 1983, 18 (2): 167–173 (in Chinese).

    [10] XU Shizhe. Calculation of vertical component of 2D gravitational field and second derivative of gravitational potential by finite element method [J]. Oil Geophysical Prospecting, 1984, 19 (5): 468–476 (in Chinese).

    [11] JIANG Xiaodian, WANG Shuoru. B spline function method for calculating gravity field of density-variable source body with arbitrary shape [J]. Oil Geophysical Prospecting, 1990, 25 (3): 362–373 (in Chinese).

    [12] KIM Kangsop, HU Xiangyun, CHO Gyonglae, et al. Study on isoparametric finite-element integral algorithm of gravity and magnetic anomaly for body with complex shape [J]. Oil Geophysical Prospecting, 2009, 44 (2): 231–239 (in Chinese).

    [13] Reeder K, Louie J, Kent G, et al. Efficient 2D finite element gravity modeling using convolution [C]. SEG Technical Program Expanded Abstracts, 2014, 33: 1254–1258.

    [14] Ren Z Y, Tang J T, Kalscheuer T, et al. Fast 3-D large-scale gravity and magnetic modeling using unstructured grids and an adaptive multilevel fast multipole method [J]. Journal of Geophysical Research: Solid Earth, 2017, 122 (1): 79–109.

    [15] LIN Baoze, XIAO Feng, WANG Mingchang. 2D gravity data radial inversion [J]. Oil Geophysical Prospecting, 2018, 53 (2): 403–409 (in Chinese).

    [16] Sharma B, Geldart L P. Analysis of gravity anomalies of two-dimensional faults using Fourier transforms [J]. Geophysical Prospecting, 1968, 16 (1): 77–93.

    [17] Rao K G C, Avasthi D N. Analysis of the Fourier spectrum of the gravity effect due to two-dimensional triangular prism [J]. Geophysical Prospecting, 1973, 21 (3): 526–542.

    [18] Bhattacharyya B K, Leu L K. Spectral analysis of gravity and magnetic anomalies due to 2-dimensional structures [J]. Geophysics, 1975, 40 (6): 993–1013.

    [19] Pedersen L B. A statistical analysis of potential fields using a vertical circular cylinder and a dike [J]. Geophysics, 1978, 43 (5): 943–953.

    [20] WU Xuanzhi. The computation of spectrum of potential field due to 3-D arbitrary bodies [J]. Chinese Journal of Geophysics, 1983, 26 (2): 177–187 (in Chinese).

    [21] Chai Yupu. Sift Sampling Theory and Its Application [M]. Petroleum Industry Press, Beijing, 1977

    [22] Tontini F C, Cocchi L, Carmisciano C. Rapid 3-D forward model of fields with application to the Palinuro Seamount gravity anomaly [J]. Journal of Geophysical Research: Solid Earth, 2009, 114 (B2): 1205–1222.

    [23] Wu L Y, Tian G. High-precision Fourier forward mode-ling of potential field [J]. Geophysics, 2014, 79 (5): G59–G68.

    [24] SHANG Yuhang, TAI Zhenhua, QIN Tao. Interface inversion in the frequency domain based on the hyperbolic density model [J]. Oil Geophysical Prospecting, 2018, 53 (4): 858–864 (in Chinese).

    [25] Xu Shizhe. Finite Element Method in Geophysics [M]. Beijing: Science Press, 1994 (in Chinese).

    [26] Xu Shiliang. FORTRAN Common Algorithm Assembly [M]. Beijing: Tsinghua University Press, 1995 (in Chinese).

    [27] Blakely R J. Potential Theory in Gravity and Magnetic Applications [M]. Cambridge University Press, London, 1996.

This Article


CN: 13-1095/TE

Vol 54, No. 06, Pages 1383-1389+1178-1179

December 2019


Article Outline



  • 0 Introduction
  • 1 Methodology
  • 2 Ordinary differential equation
  • 3 Numerical tests
  • 4 Conclusion
  • References