Multiple-point geostatistical simulation using the bunch-pasting direct sampling method

Multiple-point geostatistics has opened a new field of methodologies by which complex geological phenomena have been modeled efficiently. In this study, a modified form of direct sampling (DS) method is introduced which not only keeps the strength of DS simulation technique but also speeds it up by one or two orders of magnitude. While previous methods are based on pasting only one point at a time, here the simulation is done by pasting a bunch of nodes at a time, effectively combining the flexibility of DS with the computational advantages of patch-based methods. This bears the potential of significantly speeding up the DS method. The proposed simulation method can be used with unilateral or random simulation paths. No overlap occurs in the simulation procedure because the bunch takes the shape of the empty space around the simulated nodes. Systematic tests are carried on different training images including both categorical and continuous variables, showing that the realizations preserve the patterns existent in the training image. To illustrate the method, a Matlab implementation of the method is attached to the paper.

[1]  Pejman Tahmasebi,et al.  Multiple-point geostatistical modeling based on the cross-correlation functions , 2012, Computational Geosciences.

[2]  Hans Wackernagel,et al.  Multivariate Geostatistics: An Introduction with Applications , 1996 .

[3]  J. Chilès,et al.  Geostatistics: Modeling Spatial Uncertainty , 1999 .

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

[5]  S. De Iaco,et al.  Validation Techniques for Geological Patterns Simulations Based on Variogram and Multiple-Point Statistics , 2011 .

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

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

[8]  A. Journel,et al.  Fast FILTERSIM Simulation with Score-based Distance , 2008 .

[9]  Gregoire Mariethoz,et al.  The Direct Sampling method to perform multiple‐point geostatistical simulations , 2010 .

[10]  Marc Levoy,et al.  Fast texture synthesis using tree-structured vector quantization , 2000, SIGGRAPH.

[11]  Ragnar Hauge,et al.  Modeling of Fluvial Reservoirs with Object Models , 1998 .

[12]  Gaisheng Liu,et al.  Spatial connectivity in a highly heterogeneous aquifer: From cores to preferential flow paths , 2011 .

[13]  C. E. SHANNON,et al.  A mathematical theory of communication , 1948, MOCO.

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

[15]  Christian Lantuéjoul,et al.  Geostatistical Simulation: Models and Algorithms , 2001 .

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

[17]  Timothy C. Coburn,et al.  Geostatistics for Natural Resources Evaluation , 2000, Technometrics.

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

[19]  L. Lake,et al.  A New Approach to Shale Management in Field-Scale Models , 1984 .

[20]  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 .

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

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

[23]  J. Caers,et al.  Stochastic Simulation of Patterns Using Distance-Based Pattern Modeling , 2010 .

[24]  D. Stoyan,et al.  Stochastic Geometry and Its Applications , 1989 .

[25]  D. Stoyan,et al.  Stochastic Geometry and Its Applications , 1989 .

[26]  Julián M. Ortiz,et al.  Adapting a texture synthesis algorithm for conditional multiple point geostatistical simulation , 2011 .

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

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

[29]  Philippe Renard,et al.  Connectivity metrics for subsurface flow and transport , 2013 .

[30]  Clayton V. Deutsch,et al.  Hierarchical Object-Based Geostatistical Modeling of Fluvial Reservoirs , 1996 .

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

[32]  G. Mariéthoz,et al.  Modeling complex geological structures with elementary training images and transform‐invariant distances , 2011 .

[33]  Pejman Tahmasebi,et al.  Accelerating geostatistical simulations using graphics processing units (GPU) , 2012, Comput. Geosci..

[34]  H. Wackernagle,et al.  Multivariate geostatistics: an introduction with applications , 1998 .

[35]  W. J. Kleingeld,et al.  The Conditional Simulation of a Cox Process with Application to Deposits with Discrete Particles , 1997 .

[36]  Philippe Renard,et al.  A practical guide to performing multiple-point statistical simulations with the Direct Sampling algorithm , 2013, Comput. Geosci..