Implementation of the Singular Points Method for Gravity Data by Fast Fourier Transforms

A new implementation of the singular points method is presented in this article. The FFTs algorithm is used to replace the Fourier series used in the standard method. Because FFTs are a fast algorithm, the speed of computation is accelerated, and it can easily implement 3-D computation for grid data in the frequency domain. We have described the equations of 3-D computation. It makes the singular points method more suitable for gravity data processing. After that, the 2-D computation of the new algorithm is tested on synthetic data and on real gravity profile data from the Qian-An area in Northeast China. The 3-D computation of the algorithm is applied to the gravity grid data from Dian-Nan. Many accurate results have been obtained. It has been demonstrated that the Fourier series can be replaced by FFTs in the implementation of the singular points method. The new algorithm is more suitable for data processing because it can be used to interpret profile data as well as grid data directly. As a supplement to the singular points method, this method has the potential for widespread application. It is an effective semi-automatic interpretation method for gravity data.