COMPUTATIONALLY EFFICIENT GENERATION OF GAUSSIAN CONDITIONAL SIMULATIONS OVER REGULAR SAMPLE GRIDS

被引:8
作者
DIETRICH, CR
机构
[1] Faculty of Environmental Sciences, Griffith University, Nathan, 4111, Queensland
来源
MATHEMATICAL GEOLOGY | 1993年 / 25卷 / 04期
关键词
BLOCK TOEPLITZ STRUCTURE; CHOLESKY FACTORIZATION; GEOSTATISTICS; MONTE-CARLO SIMULATIONS; SPATIAL RANDOM FIELD;
D O I
10.1007/BF00894778
中图分类号
P [天文学、地球科学];
学科分类号
07 ;
摘要
The generation over two-dimensional grids of normally distributed random fields conditioned on available data is often required in reservoir modeling and mining investigations. Such fields can be obtained from application of turning band or spectral methods. However, both methods have limitations. First, they are only asymptotically exact in that the ensemble of realizations has the correlation structure required only if enough harmonics are used in the spectral method, or enough lines are generated in the turning bands approach. Moreover, the spectral method requires fine tuning of process parameters. As for the turning bands method, it is essentially restricted to processes with stationary and radially symmetric correlation functions. Another approach, which has the advantage of being general and exact, is to use a Cholesky factorization of the covariance matrix representing grid points correlation. For fields of large size, however, the Cholesky factorization can be computationally prohibitive. In this paper, we show that if the data are stationary and generated over a grid with regular mesh, the structure of the data covariance matrix can be exploited to significantly reduce the overall computational burden of conditional simulations based on matrix factorization techniques. A feature of this approach is its computational simplicity and suitability to parallel implementation.
引用
收藏
页码:439 / 451
页数:13
相关论文
共 19 条
[1]   BLOCK TOEPLITZ MATRIX INVERSION [J].
AKAIKE, H .
SIAM JOURNAL ON APPLIED MATHEMATICS, 1973, 24 (02) :234-241
[2]   THE PRACTICE OF FAST CONDITIONAL SIMULATIONS THROUGH THE LU DECOMPOSITION OF THE COVARIANCE-MATRIX [J].
ALABERT, F .
MATHEMATICAL GEOLOGY, 1987, 19 (05) :369-386
[3]  
[Anonymous], 2003, MULTIVARIATE STAT AN
[4]  
[Anonymous], 1978, MINING GEOSTATISTICS
[5]   NUMERICAL SOLUTION OF LINEAR EQUATIONS WITH TOEPLITZ AND VECTOR TOEPLITZ MATRICES [J].
BAREISS, EH .
NUMERISCHE MATHEMATIK, 1969, 13 (05) :404-&
[6]   SIMULATION OF ONE-DIMENSIONAL CORRELATED FIELDS USING A MATRIX-FACTORIZATION MOVING AVERAGE APPROACH [J].
BLACK, TC ;
FREYBERG, DL .
MATHEMATICAL GEOLOGY, 1990, 22 (01) :39-62
[7]  
Borgman L, 1984, GEOSTATISTICS NATURA, P517
[8]  
DAVIS MW, 1987, MATH GEOL, V19, P91
[9]   COMPUTATIONALLY EFFICIENT CHOLESKY FACTORIZATION OF A COVARIANCE-MATRIX WITH BLOCK TOEPLITZ STRUCTURE [J].
DIETRICH, CR .
JOURNAL OF STATISTICAL COMPUTATION AND SIMULATION, 1993, 45 (3-4) :203-218
[10]   CONDITIONAL SIMULATION OF INTRINSIC RANDOM FUNCTIONS OF ORDER-K [J].
DIMITRAKOPOULOS, R .
MATHEMATICAL GEOLOGY, 1990, 22 (03) :361-380