In gravity-anomaly-based prospecting, the computational and memory requirements for practical numerical modeling are potentially enormous. Achieving an efficient and precise inversion for gravity anomaly imaging over ...In gravity-anomaly-based prospecting, the computational and memory requirements for practical numerical modeling are potentially enormous. Achieving an efficient and precise inversion for gravity anomaly imaging over large-scale and complex terrain requires additional methods. To this end, we have proposed a new topography-capable By performing a two-dimensional Fourier transform in the horizontal directions, threedimensional partial differential equations in the spatial domain were transformed into a group of independent, one-dimensional differential equations engaged with different wave numbers. These independent differential equations are highly parallel across different wave numbers. differential equations with different wave numbers, and the efficiency of solving fixedbandwidth linear equations was further improved by a chasing method. In a synthetic test, a prism model was used to verify the accuracy and reliability of the proposed algorithm by comparing the numerical solution with the analytical solution. We studied the computational precision and efficiency with and without topography using different Fourier transform methods. The results showed that the Guass-FFT method has higher numerical precision, while the standard FFT method is superior, in terms of computation time, for inversion and quantitative interpretation under complicated terrain.展开更多
基金supported by the Natural Science Foundation of China(No.41574127)the China Postdoctoral Science Foundation(No.2017M622608)the project for the independent exploration of graduate students at Central South University(No.2017zzts008)
文摘In gravity-anomaly-based prospecting, the computational and memory requirements for practical numerical modeling are potentially enormous. Achieving an efficient and precise inversion for gravity anomaly imaging over large-scale and complex terrain requires additional methods. To this end, we have proposed a new topography-capable By performing a two-dimensional Fourier transform in the horizontal directions, threedimensional partial differential equations in the spatial domain were transformed into a group of independent, one-dimensional differential equations engaged with different wave numbers. These independent differential equations are highly parallel across different wave numbers. differential equations with different wave numbers, and the efficiency of solving fixedbandwidth linear equations was further improved by a chasing method. In a synthetic test, a prism model was used to verify the accuracy and reliability of the proposed algorithm by comparing the numerical solution with the analytical solution. We studied the computational precision and efficiency with and without topography using different Fourier transform methods. The results showed that the Guass-FFT method has higher numerical precision, while the standard FFT method is superior, in terms of computation time, for inversion and quantitative interpretation under complicated terrain.
