Moving-footprint-based large-scale model decomposition method for forward modeling of gravity and gravity gradient anomalies
SHI Ze-Yu1(), ZHANG Zhi-Hou1,2(), LIU Peng-Fei1, FAN Xiang-Tai1
1. Faculty of Geosciences and Environmental Engineering, Southwest Jiaotong University, Chengdu 611756, China 2. Ministry of Education Key Laboratory of High-speed Railway Engineering, Southwest Jiaotong University, Chengdu 610031, China
The computational efficiency of the forward modeling for gravity and gravity gradient anomalies determines the feasibility of inverse modeling. It also forms the basis for the efficient building of sufficient and diverse deep learning sample data. Inspired by the application of moving-footprint—a fast forward modeling method in the aerospace electromagnetic field and based on the fast space-domain forward modeling of geometric lattice functions of grid points, the authors proposed a computation method for the forward modeling of gravity and gravity gradient anomalies by applying “moving-footprint”, aiming to further improve the speed of the forward calculation for gravity and gravity gradient anomalies. Specifically, this method selects the subspace in a certain effective range directly below an observation point in the underground half-space. The observation point anomaly approximates the total anomalies of the cuboid units in the corresponding subspace while ignoring the anomalies produced by the cuboid units outside the subspace. When the observation point moves, the corresponding subspace moves accordingly. Therefore, the large-scale underground half-space cuboid model can be decomposed into the subspace corresponding to each calculation point for the forward calculation. As shown by the results of a model test, when 32×32×15 subspace was selected in the underground half-space of a 256×256×15 rectangular parallelepiped model for calculation, the relative average error of gravity anomalies and partial gradient anomalies was less than 10% and the calculation speed was increased by 19 times. Moreover, the calculation time of 1024×1024×15 rectangular parallelepiped model is approximately 32 minutes. Compared with the existing algorithms with a bottleneck in the ultra-conventional calculations, the method proposed in this study has significant advantages regarding computation.
石泽玉, 张志厚, 刘鹏飞, 范祥泰. 重力及其梯度异常正演的Moving-footprint大尺度模型分解方法[J]. 物探与化探, 2022, 46(3): 576-584.
SHI Ze-Yu, ZHANG Zhi-Hou, LIU Peng-Fei, FAN Xiang-Tai. Moving-footprint-based large-scale model decomposition method for forward modeling of gravity and gravity gradient anomalies. Geophysical and Geochemical Exploration, 2022, 46(3): 576-584.
Yao C L, Hao T Y, Guan Z N, et al. High-speed computation and efficient storage in 3D gravity and magnetic inversion based on genetic algorithms[J]. Chinese Journal of Geophysics, 2003, 46(2): 252-258.
Yang W C, Shi Z, Hou Z Z, et al. Discrete wavelet transform and multiple decomposition of gravity anomaly[J]. Chinese Journal of Geophysics, 2001, 44(4): 534-541.
Vajda P, Vaníček P, Novák P, et al. Secondary indirect effects in gravity anomaly data inversion or interpretation[J]. Journal of Geophysical Research, 2007, 112(B6):B06411(1-11).
Li Y, Oldenburg D W. Fast inversion of large-scale magnetic data using wavelet transforms and a logarithmic barrier method[J]. Geophysical Journal International, 2003, 152(2): 251-265.
Xiong G C. Several problems of three-dimensional Fourier transform of gravity and magnetic field[J]. Chinese Journal of Geophysics, 1984, 27(1): 103-109.
Shin Y H, Choi K, Xu H. Three dimensional forward and inverse models for gravity fields based on the Fast Fourier Transform[J]. Computers & Geosciences, 2006, 32(6):27-738.
Wu L Y, Tian G. High-precision Fourier forward modeling of potential fields[J]. Geophysics, 2014, 79(5): G59-G68.
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.