Efficient modelling of gravity effects due to topographic masses using the Gauss-FFT method
Efficient modelling of gravity effects due to topographic masses using the Gauss-FFT method
复制标题
使用高斯-FFT 方法对地形质量造成的重力效应进行有效建模
DOI：
10.1093/gji/ggw010
复制
发表时间：
2016
影响因子：
2.8
通讯作者：
Wu Leyuan
中科院分区：
文献类型：
--
作者：
Wu Leyuan
We present efficient Fourier-domain algorithms for modelling gravity effects due to topographic masses. The well-known Parker's formula originally based on the standard fast Fourier transform (FFT) algorithm is modified by applying the Gauss–FFT method instead. Numerical precision of the forward and inverse Fourier transforms embedded in Parker's formula and its extended forms are significantly improved by the Gauss–FFT method. The topographic model is composed of two major aspects, the geometry and the density. Versatile geometric representations, including the mass line model, the mass prism model, the polyhedron model and smoother topographic models interpolated from discrete data sets using high-order splines or pre-defined by analytical functions, in combination with density distributions that vary both laterally and vertically in rather arbitrary ways following exponential or general polynomial functions, now can be treated in a consistent framework by applying the Gauss–FFT method. The method presented has been numerically checked by space-domain analytical and hybrid analytical/numerical solutions already established in the literature. Synthetic and real model tests show that both the Gauss–FFT method and the standard FFT method run much faster than space-domain solutions, with the Gauss–FFT method being superior in numerical accuracy. When truncation errors are negligible, the Gauss–FFT method can provide forward results almost identical to space-domain analytical or semi-numerical solutions in much less time.
