Directional hydrostratigraphic units simulation using MCP algorithm

Understanding the geological uncertainty of hydrostratigraphic models is important for risk assessment in hydrogeology. An important feature of sedimentary deposits is the directional ordering of hydrostratigraphic units (HSU). Geostatistical simulation methods propose efficient algorithm for assessing HSU uncertainty. Among different geostatistical methods to simulate categorical data, Bayesian maximum entropy method (BME) and its simplified version Markov-type categorical prediction (MCP) present interesting features. In particular, the zero-forcing property of BME and MCP can provide a valuable constrain on directional properties. We illustrate the ability of MCP to simulate vertically ordered units. A regional hydrostratigraphic system with 11 HSU and different abundances is used. The transitional deterministic model of this system presents lateral variations and vertical ordering. The set of 66 (11 × 12/2) bivariate probability functions is directly calculated on the deterministic model with fast Fourier transform. Despite the trends present in the deterministic model, MCP is unbiased for the HSU proportions in the non-conditional case. In the conditional cases, MCP proved robust to datasets over-representing some HSU. The inter-realizations variability is shown to closely follow the amount and quality of data provided. Our results with different conditioning datasets show that MCP replicates adequately the directional units arrangement. Thus, MCP appears to be a practical method for generating stochastic models in a 3D hydrostratigraphic context.

[1]  Pejman Tahmasebi,et al.  Multiple-point geostatistical simulation using enriched pattern databases , 2015, Stochastic Environmental Research and Risk Assessment.

[2]  Chuanrong Zhang,et al.  A random-path markov chain algorithm for simulating categorical soil variables from random point samples , 2007 .

[3]  Thomas Tran,et al.  Improving variogram reproduction on dense simulation grids , 1994 .

[4]  Chuanrong Zhang,et al.  Two-dimensional Markov chain simulation of soil type spatial distribution , 2004 .

[5]  D. Marcotte.  Fast variogram computation with FFT , 1996 .

[6]  T. G. Orton,et al.  Accounting for the uncertainty in the local mean in spatial prediction by Bayesian Maximum Entropy , 2007 .

[7]  Patrick Bogaert,et al.  Estimating soil properties from thematic soil maps: the Bayesian maximum entropy approach , 2002 .

[8]  D. D'Or,et al.  Application of the BME approach to soil texture mapping , 2001 .

[9]  Alexandre Boucher,et al.  Simulation of Geological Contacts from Interpreted Geological Model Using Multiple-Point Statistics , 2014, Mathematical Geosciences.

[10]  Marc L. Serre,et al.  Modern geostatistics: computational BME analysis in the light of uncertain physical knowledge – the Equus Beds study , 1999 .

[11]  P. Bogaert.  Spatial prediction of categorical variables: the Bayesian maximum entropy approach , 2002 .

[12]  R. A. Leibler,et al.  On Information and Sufficiency , 1951 .

[13]  R. Froidevaux,et al.  An efficient maximum entropy approach for categorical variable prediction , 2011 .

[14]  Wenpeng Li,et al.  A review of regional groundwater flow modeling , 2011 .

[15]  Olivier Barreteau,et al.  Integrated Groundwater Management: Concepts, Approaches and Challenges , 2016 .

[16]  E. Frind,et al.  On the use of mean groundwater age, life expectancy and capture probability for defining aquifer vulnerability and time-of-travel zones for source water protection. , 2012, Journal of contaminant hydrology.

[17]  Sam Savage,et al.  Accounting for Uncertainty , 2002 .

[18]  Characterization of High Order Correlation for Enhanced Indicator Simulation , 2001 .

[19]  Patrick Bogaert,et al.  Spatial prediction of categorical variables with the Bayesian Maximum Entropy approach: the Ooypolder case study , 2004 .

[20]  Margaret Armstrong,et al.  Plurigaussian Simulations in Geosciences , 2014 .

[21]  Jens Christian Refsgaard,et al.  Review of strategies for handling geological uncertainty in groundwater flow and transport modeling , 2012 .

[22]  Dimitri D'Or.  Spatial prediction of soil properties : the Bayesian Maximum Entropy approach , 2003 .

[23]  George Christakos,et al.  Random Field Models in Earth Sciences , 1992 .

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

[25]  A. Jakeman,et al.  Integrated Groundwater Management : An Overview of Concepts and Challenges 1 , 2018 .

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

[27]  Carlo Giupponi.  Integrated Groundwater Management: Concepts, Approaches and Challenges , 2017, Environ. Model. Softw..

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

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

[30]  Weidong Li,et al.  Markov Chain Random Fields for Estimation of Categorical Variables , 2007 .

[31]  Cass T. Miller,et al.  A BME solution of the inverse problem for saturated groundwater flow , 2003 .

[32]  O. Dubrule,et al.  Modelling Asymmetrical Facies Successions Using Pluri-Gaussian Simulations , 2017 .

[33]  Omid Asghari,et al.  Multiple-point geostatistical simulation of dykes: application at Sungun porphyry copper system, Iran , 2014, Stochastic Environmental Research and Risk Assessment.

[34]  Benoit Noetinger,et al.  The FFT Moving Average (FFT-MA) Generator: An Efficient Numerical Method for Generating and Conditioning Gaussian Simulations , 2000 .