Fast potential field modeling based on digital filtering

Abstract

In recent years many papers dealing with potential field modelling have been published. The main disadvantages of some methods using analytical expressions for potential field spectra are considered. A universal modelling algorithm for potential fields of regularly constructed discrete models is presented. It is based on the discrete convolution thorem. To insulate the output from artifacts of undesired periodicity a complete treatment of edge effects is described. In order to speed up the computation, the fast Fourier transform and its symmetries are utilized. We are dealing with a discrete convolution of two well behaved functions; more over, one of them having final support. By computing Green's function accurately, no other errors except for round-off errors can affect the data within the desired output range.