The Direct Sampling method to perform multiple‐point geostatistical simulations

Multiple‐point geostatistics is a general statistical framework to model spatial fields displaying a wide range of complex structures. In particular, it allows controlling connectivity patterns that have a critical importance for groundwater flow and transport problems. This approach involves considering data events (spatial arrangements of values) derived from a training image (TI). All data events found in the TI are usually stored in a database, which is used to retrieve conditional probabilities for the simulation. Instead, we propose to sample directly the training image for a given data event, making the database unnecessary. Our method is statistically equivalent to previous implementations, but in addition it allows extending the application of multiple‐point geostatistics to continuous variables and to multivariate problems. The method can be used for the simulation of geological heterogeneity, accounting or not for indirect observations such as geophysics. We show its applicability in the presence of complex features, nonlinear relationships between variables, and with various cases of nonstationarity. Computationally, it is fast, easy to parallelize, parsimonious in memory needs, and straightforward to implement.

[1]  G. Matheron.  Éléments pour une théorie des milieux poreux , 1967 .

[2]  R. Freeze.  A stochastic‐conceptual analysis of one‐dimensional groundwater flow in nonuniform homogeneous media , 1975 .

[3]  G. Dagan.  Comment on `Stochastic-conceptual analysis of one-dimensional groundwater flow in nonuniform homogeneous media' by R. Allan Freeze , 1976 .

[4]  D. K. Pickard,et al.  Unilateral Markov fields , 1980, Advances in Applied Probability.

[5]  A. Journel.  Nonparametric estimation of spatial distributions , 1983 .

[6]  G. Dagan.  Statistical Theory of Groundwater Flow and Transport: Pore to Laboratory, Laboratory to Formation, and Formation to Regional Scale , 1986 .

[7]  G. Dagan.  Flow and transport in porous formations , 1989 .

[8]  Joel Massmann,et al.  Hydrogeological Decision Analysis: 1. A Framework , 1990 .

[9]  L. Gelhar.  Stochastic Subsurface Hydrology , 1992 .

[10]  R. M. Srivastava,et al.  Multivariate Geostatistics: Beyond Bivariate Moments , 1993 .

[11]  Jesús Carrera,et al.  Scale effects in transmissivity , 1996 .

[12]  S. Gorelick,et al.  Heterogeneity in Sedimentary Deposits: A Review of Structure‐Imitating, Process‐Imitating, and Descriptive Approaches , 1996 .

[13]  Xian-Huan Wen,et al.  Numerical modeling of macrodispersion in heterogeneous media: a comparison of multi-Gaussian and non-multi-Gaussian models , 1998 .

[14]  E. Langsholt,et al.  Development of Three‐Dimensional Hydrostratigraphical Architecture of the Unsaturated Zone Based on Soft and Hard Data , 1998 .

[15]  J. Gómez-Hernández,et al.  To be or not to be multi-Gaussian? A reflection on stochastic hydrogeology , 1998 .

[16]  Weilin Huang,et al.  Hysteresis in the sorption and desorption of hydrophobic organic contaminants by soils and sediments: 1. A comparative analysis of experimental protocols , 1998 .

[17]  R. Grayson,et al.  Toward capturing hydrologically significant connectivity in spatial patterns , 2001 .

[18]  You‐Kuan Zhang.  Stochastic Methods for Flow in Porous Media: Coping with Uncertainties , 2001 .

[19]  A. Journel.  Combining Knowledge from Diverse Sources: An Alternative to Traditional Data Independence Hypotheses , 2002 .

[20]  S. Strebelle,et al.  Modeling of a Deepwater Turbidite Reservoir Conditional to Seismic Data Using Principal Component Analysis and Multiple-Point Geostatistics , 2003 .

[21]  Sebastien Strebelle,et al.  Conditional Simulation of Complex Geological Structures Using Multiple-Point Statistics , 2002 .

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

[23]  Charles F. Harvey,et al.  When good statistical models of aquifer heterogeneity go bad: A comparison of flow, dispersion, and mass transfer in connected and multivariate Gaussian hydraulic conductivity fields , 2003 .

[24]  J. Caers.  Interpreter's Corner—Stochastic integration of seismic data and geologic scenarios: A West Africa submarine channel saga , 2003 .

[25]  Gedeon Dagan,et al.  On application of stochastic modeling of groundwater flow and transport , 2004 .

[26]  Yuhong Liu,et al.  Multiple-point simulation integrating wells, three-dimensional seismic data, and geology , 2004 .

[27]  C.A.J. Appelo,et al.  Flow and Transport , 2004, Introductory Incompressible Fluid Mechanics.

[28]  L. Feyen,et al.  Reliable groundwater management in hydroecologically sensitive areas , 2004 .

[29]  P. Renard,et al.  Dealing with spatial heterogeneity , 2005 .

[30]  Sebastien Strebelle,et al.  Post-processing of Multiple-point Geostatistical Models to Improve Reproduction of Training Patterns , 2005 .

[31]  L. Wasserman.  All of Nonparametric Statistics , 2005 .

[32]  Jesús Carrera,et al.  On the relationship between indicators of geostatistical, flow and transport connectivity , 2005 .

[33]  Colin Daly,et al.  Higher Order Models using Entropy, Markov Random Fields and Sequential Simulation , 2005 .

[34]  Louis J. Durlofsky,et al.  Optimizing the performance of smart wells in complex reservoirs using continuously updated geological models , 2005 .

[35]  Olaf A. Cirpka,et al.  Homogenization of Richards equation in permeability fields with different connectivities , 2005 .

[36]  L. Feyen,et al.  Quantifying geological uncertainty for flow and transport modeling in multi-modal heterogeneous formations , 2006 .

[37]  A. Journel,et al.  The Necessity of a Multiple-Point Prior Model , 2007 .

[38]  Sang Joon Kim,et al.  A Mathematical Theory of Communication , 2006 .

[39]  J. Caers,et al.  The Probability Perturbation Method: A New Look at Bayesian Inverse Modeling , 2006 .

[40]  X. Sanchez‐Vila,et al.  Representative hydraulic conductivities in saturated groundwater flow , 2006 .

[41]  Paul Switzer,et al.  Filter-Based Classification of Training Image Patterns for Spatial Simulation , 2006 .

[42]  Algorithm-driven and Representation-driven Random Function : A New Formalism for Applied Geostatistics , 2007 .

[43]  J. Caers,et al.  Conditional Simulation with Patterns , 2007 .

[44]  Xavier Emery.  Using the Gibbs sampler for conditional simulation of Gaussian-based random fields , 2007, Comput. Geosci..

[45]  Philippe Renard,et al.  Stochastic Hydrogeology: What Professionals Really Need? , 2007, Ground water.

[46]  S. Strebelle,et al.  Real-time Post-Processing Method to Enhance Multiple-Point Statistics Simulation , 2007 .

[47]  J. Caers,et al.  History matching by jointly perturbing local facies proportions and their spatial distribution: Application to a North Sea reservoir , 2007 .

[48]  Ragnar Hauge,et al.  Modification of the Snesim Algorithm , 2007 .

[49]  Antoine Saucier,et al.  A patchwork approach to stochastic simulation: A route towards the analysis of morphology in multiphase systems , 2008 .

[50]  J. Caers,et al.  A Distance-based Prior Model Parameterization for Constraining Solutions of Spatial Inverse Problems , 2008 .

[51]  A. Dassargues,et al.  Application of multiple-point geostatistics on modelling groundwater flow and transport in a cross-bedded aquifer (Belgium) , 2009 .

[52]  Jesús Carrera,et al.  Application of Multiple Point Geostatistics to Non-stationary Images , 2008 .

[53]  P. Renard,et al.  Issues in characterizing heterogeneity and connectivity in non-multiGaussian media , 2008 .

[54]  L. Hu,et al.  Multiple-Point Simulations Constrained by Continuous Auxiliary Data , 2008 .

[55]  L. Y. Hu,et al.  Multiple‐point geostatistics for modeling subsurface heterogeneity: A comprehensive review , 2008 .

[56]  Peter Lehmann,et al.  Measuring the effect of structural connectivity on the water dynamics in heterogeneous porous media using speedy neutron tomography , 2008 .

[57]  Alexandre Boucher,et al.  A SGeMS code for pattern simulation of continuous and categorical variables: FILTERSIM , 2008, Comput. Geosci..

[58]  S. Gorelick,et al.  Identifying discrete geologic structures that produce anomalous hydraulic response: An inverse modeling approach , 2008 .

[59]  Philippe Renard,et al.  Reducing the impact of a desalination plant using stochastic modeling and optimization techniques , 2009 .

[60]  Carlos Llopis-Albert,et al.  Gradual conditioning of non-Gaussian transmissivity fields to flow and mass transport data: 1. Theory , 2009 .

[61]  Philippe Renard,et al.  Integrating collocated auxiliary parameters in geostatistical simulations using joint probability distributions and probability aggregation , 2009 .

[62]  Katherine A. Klise,et al.  Exploring solute transport and streamline connectivity using lidar‐based outcrop images and geostatistical representations of heterogeneity , 2009 .

[63]  Alexandre Boucher,et al.  Considering complex training images with search tree partitioning , 2009, Comput. Geosci..

[64]  S. Gorelick,et al.  Combining geologic‐process models and geostatistics for conditional simulation of 3‐D subsurface heterogeneity , 2010 .

[65]  A. Dassargues,et al.  Application of Multiple-Point Geostatistics on Modelling Groundwater Flow and Transport in a Cross-Bedded Aquifer , 2010 .

[66]  Grégoire Mariethoz,et al.  A general parallelization strategy for random path based geostatistical simulation methods , 2010, Comput. Geosci..

[67]  G. Mariéthoz,et al.  Reconstruction of Incomplete Data Sets or Images Using Direct Sampling , 2010 .

[68]  G. Mariéthoz,et al.  Bayesian inverse problem and optimization with iterative spatial resampling , 2010 .

[69]  Philippe Renard,et al.  Blocking Moving Window algorithm: Conditioning multiple‐point simulations to hydrogeological data , 2010 .

[70]  G. Mariéthoz,et al.  An Improved Parallel Multiple-point Algorithm Using a List Approach , 2011 .