Forward modeling for 2D-body gravity anomaly in[JP2]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 410083, China; 2. School of Geosciences and Info-physics, Central South University, Changsha, Hunan 410083, China; 3. Guangdong Electric Power Design Institute, China Energy Engineering Group CO., Guangzhou, Guangdong 510663, China; 4. College of Earth Sciences, Guilin University of Technology, Guilin, Guangxi 541004, China; 5. No. 1 Daqing Geophysical Exploration Company of BGP, CNPC, Daqing, Heilongjiang 163357, China
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.By performing 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,thus 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.1D finite element method is used to solve the transformed differential equations with different wave numbers.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 also improves the calculation efficiency compared with the conventional methods.
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 212-D 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]
李明,王宜昌.重力垂直一次导数的计算及应用[J].石油地球物理勘探,1983,18(2):167-173.LI Ming,WANG Yichang.The computation and application of vertical first derivative of gravity[J].Oil Geophysical Prospecting,1983,18(2):167-173.
[10]
徐世浙.用有限元法计算二维重力场垂直分量及重力位二阶导数[J].石油地球物理勘探,1984,19(5):468-476.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.
[11]
姜效典,王硕儒.计算变密度任意形状地质体重力异常的B样条函数法[J].石油地球物理勘探,1990,25(3):362-373.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.
[12]
金钢燮,胡祥云,超敬来,等.复杂形体重磁异常的等参数有限元积分算法研究[J].石油地球物理勘探,2009,44(2):231-239.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.
[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]
林宝泽,肖锋,王明常.二维重力数据径向反演及应用[J].石油地球物理勘探,2018,53(2):403-409.LIN Baoze,XIAO Feng,WANG Mingchang.2D gra-vity data radial inversion[J].Oil Geophysical Prospecting,2018,53(2):403-409.
[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]
吴宣志.三度体位场波谱的正演计算[J].地球物理学报,1983,26(2):177-187.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.
[21]
柴玉璞.偏移抽样理论及其应用[M].北京:石油工业出版社,1997.
[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]
商宇航,邰振华,秦涛.基于双曲线密度模型的频率域界面反演[J].石油地球物理勘探,2018,53(4):858-864.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.
[25]
徐世浙.地球物理中的有限单元法[M].北京:科学出版社,1994.
[26]
徐士良.FORTRAN常用算法程序集[M].北京:清华大学社出版社,1995.
[27]
Blakely R J.Potential Theory in Gravity and Magnetic Applications[M].Cambridge University Press,London,1996.