×

Moving averages for Gaussian simulation in two and three dimensions. (English) Zbl 0970.86549

Summary: The square-root method provides a simple and computationally inexpensive way to generate multidimensional Gaussian random fields. It is applied by factoring the multidimensional covariance operator analytically, then sampling the factorization at discrete points to compute an array of weighted averages that can be convolved with an array of random normal deviates to generate a correlated random field. In many respects this is similar to the \(LU\) decomposition method and to the one-dimensional method of moving averages. However it has been assumed that the method of moving averages could not be used in higher dimensions, whereas direct application of the matrix decomposition approach is too expensive to be practical on large grids. In this paper, it is shown that it is possible to calculate the square root of many two- and three-dimensional covariance operators analytically so that the method of moving averages can be applied directly to the problem of multidimensional simulation. A few numerical examples of nonconditional simulation on a \(256\times 256\) grid that show the simplicity of the method are included. The method is fast and can be applied easily to nested and anisotropic variograms.

MSC:

86A32 Geostatistics
62M99 Inference from stochastic processes
65F99 Numerical linear algebra
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Alabert, F., 1987, The practice of fast conditional simulation through the LU decomposition of the covariance matrix: Math. Geology, v. 19, no. 5, p. 369–386.
[2] Armstrong, M., and Diamond, P., 1984, Testing variograms for positive-definiteness: Math. Geology, v. 16, no. 4, p. 407–421.
[3] Bracewell, R. N., 1978, The Fourier transform and its applications: McGraw-Hill Book Co., New York, 444 p. · Zbl 0502.42001
[4] Davis, M., 1987a, Production of conditional simulations via the LU decomposition of the covariance matrix: Math. Geology, v. 19, no. 2, p. 91–98.
[5] Davis, M., 1987b. Generating large stochastic simulations–the matrix polynomial approximation method: Math. Geology, v. 19, no. 2, p. 99–107.
[6] Dietrich, C. R., 1993, Computationally efficient generation of Gaussian conditional simulations over regular sample grids: Math. Geology, v. 25, no. 4, p. 439–451. · Zbl 0970.86538
[7] Dietrich, C. R., and Newsam, G. N., 1993, A fast and exact method for multidimensional Gaussian stochastic simulations: Water Resour. Res., v. 29, no. 8, p. 2861–2869.
[8] Dietrich, C. R., and Newsam, G. N., 1995, Efficient generation of conditional simulations by Chebyshev matrix polynomial approximations to the symmetric square root of the covariance matrix: Math. Geology, v. 27, no. 2, p. 207–228. · Zbl 0970.86547
[9] Gradshteyn, I. S., and Ryzhik, I. M., 1980, Table of integrals, series, and products (corrected and enlarged edition): Academic Press, New York, 1160 p. · Zbl 0521.33001
[10] Journel, A. G., 1974, Geostatistics for conditional simulation of ore bodies: Econ. Geology, v. 69, p. 673–687.
[11] Journel, A., and Huijbregts, C. J., 1978, Mining geostatistics: Academic Press, New York, 600 P.
[12] Luster, G. R., 1985, Raw materials for portland cement: applications of conditional simulation of coregionalization: unpubl. doctoral dissertation, Stanford University, 531 p.
[13] Matheron, G., 1971, The theory of regionalized variables and its applications: Tech. Rept. No. 5, Les Cahiers du Centre de Morphologie Mathématique de Fontainebleau, 211 p.
[14] Rao, C. R., 1965, Linear statistical inference and its applications: John Wiley & Sons, New York, 522 p. · Zbl 0137.36203
[15] Riesz, F., and Sz.-Nagy, B., 1955, Functional analysis: Frederick Ungar Publ. Co., New York, 468 p.
[16] Sneddon, I. N., 1951, Fourier transforms: McGraw-Hill Book Co., New York, 542 p. · Zbl 0038.26801
[17] Vanmarcke, E., 1983, Random fields: analysis and synthesis: MIT Press, Cambridge, Massachusetts, 382 p. · Zbl 1206.60051
[18] Whittle, P., 1954, On stationary processes in the plane: Biometrika, v. 41, p. 434–449. · Zbl 0058.35601
[19] Yaglom, A. M., 1962, An introduction to the theory of stationary random functions: Prentice-Hall, Englewood, Cliffs, New Jersey, 235 p. · Zbl 0121.12601
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.