The Fourier Integral Method: An efficient spectral method for simulation of random fields

E. Pardo-Igúzquiza1, M. Chica-Olmo1
1Department of Geodynamics/IAGM, Avda. Fuentenueva s/n, Facultad de Ciencias, University of Granada, Granada, Spain

Tóm tắt

The Fourier Integral Method (FIM) of spectral simulation, adapted to generate realizations of a random function in one, two, or three dimensions, is shown to be an efficient technique of non-conditional geostatistical simulation. The main contribution is the use of the fast Fourier transform for both numerical calculus of the density spectral function and as generator of random finite multidimensional sequences with imposed covariance. Results obtained with the FIM are compared with those obtained by other classic methods: Shinozuka and Jan Method in 1D and Turning Bands Method in 2D and 3D, the points for and against different methodologies are discussed. Moreover, with the FIM the simulation of nested structures, one of which can be a nugget effect and the simulation of both zonal and geometric anisotropy is straightforward. All steps taken to implement the FIM methodology are discussed.

Tài liệu tham khảo

Anderson, T. W., 1971, The Statistical Analysis of Times Series: John Wiley and Sons, New York, p. 704. Black, T. C., and Freyberg, D. L., 1990, Simulation of One-Dimensional Correlated Fields Using a Matrix-Factorization Moving Average Approach: Math. Geol., v. 22, p. 39–62. Borgman, L., Taheri, M., and Hagan, R., 1984, Three-Dimensional Frequency-Domain Simulations of Geological Variables: Proc. NATO ASI, Lake Tahoe, Sept. 1983, p. 517–541. Bracewell, R. N., 1986, The Fourier Transform and Its Applications: Mc-Graw Hill International Editions, Singapore, p. 474. Brigham, E. O., 1988, The Fast Fourier Transform and Its Applications: Prentice-Hall International Editions, London, p. 448. Brooker, P. I., 1985, Two-Dimensional Simulation by Turning Bands: Math. Geol., v. 17, p. 81–90. Chatfield, C., 1991, The Analysis of Time Series, 4th ed.: Chapman and Hall, London, p. 241. Christakos, G., 1987, Stochastic Simulation of Spatially Correlated Geo-Processes: Math. Geol., v. 19, p. 807–831. Cox, D. R., and Miller, H. D., 1968, The Theory of Stochastic Processes: Methuen & Co. Ltd., London, p. 397. Dowd, P., 1991, A Review of Recent Developments in Geostatistics: Comp. Geosci., v. 17, p. 1481–1500. Fox, G. C., 1987, An Inverse Fourier Transform Algorithm for Generating Random Signals of a Specified Spectral Form: Comp. Geosci., v. 13, p. 369–374. Hsu, H. P., and Mehra, R., 1973, Análisis de Fourier: McGraw Hill, Bogotá, p. 274. Isaaks, E. H., and Srisvastava, R. M., 1989, Applied Geostatistics: Oxford University Press, New York, p. 561. Journel, A., 1974, Geostatistics for Conditional Simulation of Ore Bodies: Econ. Geol., v. 69, p. 673–687. Journel, A., and Huijbregts, C., 1978, Mining Geostatistics: Academic Press, New York, p. 600. Mantoglou, A., 1987, Digital Simulation of Multivariate Two- and Three-Dimensional Stochastic Processes with a Spectral Turning Bands Method: Math. Geol., v. 19, p. 129–149. Mantoglou, A., and Wilson, J. L., 1982, The Turning Band Method for Simulation of Random Fields Using Line Generation by a Spectral Method: Water Res. Res., v. 18, p. 1379–1394. Matheron, G., 1973, The Intrinsic Random Functions and Their Applications: Adv. Appl. Prob., v. 5, p. 439–468. Pardo-Igúzquiza, E., 1991, Simulación Geoestadística de Variables Geológicas por Métodos Espectrales. Ph.D. Thesis, Department of Geodynamics, University of Granada, p. 412. Schrage, L., 1979, A More Portable Fortran Random Number Generator: ACM Trans. Math. Soft., v. 5, n. 2, p. 132–138. Shinozuka, M., and Jan, C. M., 1972, Digital Simulation of Random Processes and Its Applications: J. Sound Vibr., v. 25, p. 111–128. Weaver, H. J., 1989, Theory of Discrete and Continuous Fourier Analysis: John Wiley and Sons, New York, p. 307.