General joint conditional simulations using a fast fourier transform method

A procedure for generating joint statistically homogeneous random fields is examined. The method is based on the spectral representation theorem. It handles large fields easily and is both rapid and flexible. Algorithm development and examples are presented. The procedure is adapted further to include the possibility of generating fields that are jointly conditioned on data from two related fields.

[1]  T. W. Anderson.  An Introduction to Multivariate Statistical Analysis , 1959 .

[2]  James R. Carr,et al.  COSIM: a FORTRAN IV program for coconditional simulation , 1983 .

[3]  M. Stein,et al.  A Bayesian analysis of kriging , 1993 .

[4]  Allan L. Gutjahr,et al.  Cross‐correlated random field generation with the direct Fourier Transform Method , 1993 .

[5]  A. G. Journel,et al.  Indicator principal component kriging , 1991 .

[6]  Andre G. Journel,et al.  Joint simulation of multiple variables with a Markov-type coregionalization model , 1994 .

[7]  M. Chica-Olmo,et al.  Spectral simulation of multivariable stationary random functions using covariance fourier transforms , 1994 .

[8]  M. W. Davis,et al.  Production of conditional simulations via the LU triangular decomposition of the covariance matrix , 1987, Mathematical Geology.

[9]  L. H. Koopmans.  The spectral analysis of time series , 1974 .

[10]  R. Ababou,et al.  Implementation of the three‐dimensional turning bands random field generator , 1989 .

[11]  M. W. Davis.  Generating large stochastic simulations—The matrix polynomial approximation method , 1987 .

[12]  Aristotelis Mantoglou,et al.  Digital simulation of multivariate two- and three-dimensional stochastic processes with a spectral turning bands method , 1987, Mathematical Geology.

[13]  A. Yaglom.  Correlation Theory of Stationary and Related Random Functions I: Basic Results , 1987 .

[14]  G. Newsam,et al.  Efficient generation of conditional simulations by chebyshev matrix polynomial approximations to the symmetric square root of the covariance matrix , 1995 .

[15]  Probability of plume interception using conditional simulation of hydraulic head and inverse modeling , 1991 .

[16]  A. Mantoglou,et al.  The Turning Bands Method for simulation of random fields using line generation by a spectral method , 1982 .

[17]  Leon E. Borgman,et al.  Three-Dimensional, Frequency-Domain Simulations of Geological Variables , 1984 .

[18]  G. Matheron.  The intrinsic random functions and their applications , 1973, Advances in Applied Probability.

[19]  Allan L. Gutjahr,et al.  Joint Conditional Simulations and Flow Modeling , 1994 .

[20]  M. Chica-Olmo,et al.  The Fourier Integral Method: An efficient spectral method for simulation of random fields , 1993 .

[21]  J. Gómez-Hernández,et al.  Joint Sequential Simulation of MultiGaussian Fields , 1993 .

[22]  M. G. Marietta,et al.  Results from a Comparison of Geostatistical Inverse Techniques for Groundwater Flow , 1996 .

[23]  Sean Hatch,et al.  Joint conditional simulations and the spectral approach for flow modeling , 1994 .

[24]  Donald E. Myers.  Vector Conditional Simulation , 1989 .