Latin hypercube lattice sample selection strategy for correlated random hydraulic conductivity fields

[1] Given an established probabilistic description of hydraulic conductivity, realizations of the hydraulic conductivity field can be generated using a lattice sampling technique, a special case of Latin hypercube sampling, and subsequently input to a groundwater flow and transport model. Realizations of dependent variables such as contaminant concentration are obtained as the output of the model and their corresponding statistics can be calculated. Compared with the other three alternative methods for generating a random hydraulic conductivity field, the suggested method is more efficient in terms of computational effort because it needs fewer samples to achieve the same statistical accuracy.

[1]  A. Owen.  A Central Limit Theorem for Latin Hypercube Sampling , 1992 .

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

[3]  B. M. Davis.  Uses and abuses of cross-validation in geostatistics , 1987 .

[4]  R. L. Iman,et al.  An iterative algorithm to produce a positive definite correlation matrix from an approximate correlation matrix (with a program user's guide) , 1982 .

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

[6]  John R. Nuckols,et al.  Stochastic modeling of exposure and risk in a contaminated heterogeneous aquifer. 2 : Application of Latin Hypercube Sampling , 1999 .

[7]  R. Iman,et al.  A distribution-free approach to inducing rank correlation among input variables , 1982 .

[8]  G. Heuvelink,et al.  Latin Hypercube Sampling of Gaussian Random Fields , 1999, Technometrics.

[9]  H. D. Patterson.  The Errors of Lattice Sampling , 1954 .

[10]  J. P. Delhomme,et al.  Spatial variability and uncertainty in groundwater flow parameters: A geostatistical approach , 1979 .

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

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

[13]  E. Paleologos,et al.  Stochastic Methods for Flow in Porous Media, Coping With Uncertainties , 2003 .

[14]  Clayton V. Deutsch,et al.  GSLIB: Geostatistical Software Library and User's Guide , 1993 .

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

[16]  Clayton V. Deutsch,et al.  Geostatistical Software Library and User's Guide , 1998 .

[17]  S. P. Neuman,et al.  Effects of kriging and inverse modeling on conditional simulation of the Avra Valley Aquifer in southern Arizona , 1982 .

[18]  Francois Alabert,et al.  The practice of fast conditional simulations through the LU decomposition of the covariance matrix , 1987 .

[19]  M. Stein.  Large sample properties of simulations using latin hypercube sampling , 1987 .

[20]  M. D. McKay,et al.  A comparison of three methods for selecting values of input variables in the analysis of output from a computer code , 2000 .

[21]  A. Owen.  Controlling correlations in latin hypercube samples , 1994 .