Fast potential field modeling based on digital filtering

Article Sidebar

PDF (Español (España))
Published: Oct 1, 1994
DOI: https://doi.org/10.22201/igeof.00167169p.1994.33.4.586
Keywords:
Fast potential field modelling, FFT

Main Article Content

Uwe Koppelt
Universität Leipzig Institut für Geophysik und Geologie, Leipzig, Germany
Javier Rojas
Universität Leipzig Institut für Geophysik und Geologie, Leipzig, Germany

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. 

Article Details

How to Cite
Koppelt, U., & Rojas, J. (1994). Fast potential field modeling based on digital filtering. Geofisica Internacional, 33(4), 523–529. https://doi.org/10.22201/igeof.00167169p.1994.33.4.586
Issue
Vol. 33 No. 4 (1994): October 1, 1994
Section
Article
Creative Commons License

This work is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License.

References

BEZVODA, V., 1987. Fast algorithms in geophysics and geochemics (in czech). Doctor thesis, Charles University, Prague.

BEZVODA, V., J. HRABE and K. SEGETH, 1992. Fast algorithm for solving potential field problems. In: Vogel A., Sarwar A. K. M., Gorenflo R., Kouchev O. J. (Eds.): Theory and practice of geophysical data inversion. Vieweg, Braunschweig. DOI: https://doi.org/10.1007/978-3-322-89417-5_7

BHATTACHARYYA, B. K. and M. E. NAVOLIO, 1975. Digital convolution for computing gravity and magnetic anomalies due to arbitrary bodies. Geophysics 40, 981-992. DOI: https://doi.org/10.1190/1.1440592

BRACEWELL, R., 1965. The Fourier transform method and its application. McGraw-Hill, New York.

CHAI, Y. and W. J. HINZE, 1988. Gravity inversion of an interface above which the density contrast varies exponentially with depth. Geophysics 53, 837-846. DOI: https://doi.org/10.1190/1.1442518

ELLIOT, D. F., 1987. Handbook of digital signal processing, Engineering applications. Academic press, London.

HANSEN, R. O. and X. WANG, 1988. Simplified frequency domain expression for potential fields of arbitrary three-dimensional bodies. Geophysics 53, 365-374. DOI: https://doi.org/10.1190/1.1442470

KANASEWICH, E. R., 1981. Time sequence analysis in geophysics. University of Alberta Press, Edmonton.

MESKÓ, A., 1977. A new algorithm for the computing of gravitational attraction due to irregulary shaped bodies. Ann. Univ. Budapest 18,65-73.

NAGY, D., 1988. A short program for three-dimensional gravity modelling. Acta Geod. Geoph. Mont. Hung. 23, 449-459.

PARKER, R. L., 1973. The rapid compution of potential anomalies. Geophys. Journ. Roy. Astr. Soc. 31, 447-455. DOI: https://doi.org/10.1111/j.1365-246X.1973.tb06513.x

PEDERSEN, L. B., 1978. Wavenumber expression for potential fields from arbitrary 2-2- and 3- dimensional bodies. Geophysics 43, 626-630. DOI: https://doi.org/10.1190/1.1440841

SIDERIS, M. G., 1985. A fast Fourier transform method for computing terrain corrections. Manuscripta Geodaetica 10, 66-73. DOI: https://doi.org/10.1007/BF03655071

TALWANI, M., J. L. WORZEL and M. LANDISMAN, 1959. Rapid gravity computation for two-dimensional bodies with application to the Mendocino submarine fracture zone. J. Geophys. Res. 64, 49-59. DOI: https://doi.org/10.1029/JZ064i001p00049