Calibration of channelized subsurface flow models using nested sampling and soft probabilities

Abstract A new method for calibration of channelized subsurface flow models is presented. The proposed method relies on the nested sampling algorithm and on adaptive construction of soft probability maps. Nested sampling (NS) is a Bayesian sampling algorithm for estimating the Bayesian evidence and obtaining samples from the posterior distribution of the unknown fields. NS utilizes a set of samples (active set) that evolves to high-likelihood regions. The sample evolution process is achieved by iteratively replacing the sample with the lowest likelihood within the active set by a new sample from the prior but with higher likelihood value (constrained sampling). For channelized models, drawing samples from the prior model based on a training image and only accepting the samples satisfying the likelihood constraint is computationally inefficient due to low acceptance rates. We develop an efficient constrained sampling step utilizing soft probability maps in addition to the training image (using the Tau model) to obtain samples from the prior satisfying the likelihood constraint. The soft probability map is constructed by averaging the samples within the active set and is shown to significantly increase the acceptance rate of the nested sampling algorithm. The proposed algorithm is applied for calibration of several channelized subsurface flow models. In addition, the NS algorithm is applied for prior model selection of a channelized model with different training images obtained by changing the orientation angles of a reference training image. The results show that selecting the prior model based on the data mismatch can be misleading. This highlights the need to evaluate the Bayesian evidence (estimated by the nested sampling algorithm) as a more reliable prior model selection statistics, especially when the amount of calibration data is limited.

[1]  Mary F. Wheeler,et al.  Boosting iterative stochastic ensemble method for nonlinear calibration of subsurface flow models , 2013 .

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

[3]  Haiyan Zhou,et al.  A pattern‐search‐based inverse method , 2012 .

[4]  Zhangxin Chen,et al.  Reservoir simulation : mathematical techniques in oil recovery , 2007 .

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

[6]  Yalchin Efendiev,et al.  An efficient two‐stage Markov chain Monte Carlo method for dynamic data integration , 2005 .

[7]  Behnam Jafarpour,et al.  A Probability Conditioning Method (PCM) for Nonlinear Flow Data Integration into Multipoint Statistical Facies Simulation , 2011 .

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

[9]  D. McLaughlin,et al.  A Reassessment of the Groundwater Inverse Problem , 1996 .

[10]  Knud Skou Cordua,et al.  SIPPI: A Matlab toolbox for sampling the solution to inverse problems with complex prior information: Part 1 - Methodology , 2013, Comput. Geosci..

[11]  A Blocking Markov Chain Monte Carlo Method for Inverse Stochastic Hydrogeological Modeling , 2009 .

[12]  Mary F. Wheeler,et al.  Nested sampling algorithm for subsurface flow model selection, uncertainty quantification, and nonlinear calibration , 2013 .

[13]  H. Akaike.  A new look at the statistical model identification , 1974 .

[14]  D. Parkinson,et al.  A Nested Sampling Algorithm for Cosmological Model Selection , 2005, astro-ph/0508461.

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

[16]  Soroosh Sorooshian,et al.  Dual state-parameter estimation of hydrological models using ensemble Kalman filter , 2005 .

[17]  Mary F. Wheeler,et al.  Efficient Bayesian inference of subsurface flow models using nested sampling and sparse polynomial chaos surrogates , 2014 .

[18]  Holger Class,et al.  Bayesian updating via bootstrap filtering combined with data-driven polynomial chaos expansions: methodology and application to history matching for carbon dioxide storage in geological formations , 2013, Computational Geosciences.

[19]  Behnam Jafarpour,et al.  A Bayesian mixture‐modeling approach for flow‐conditioned multiple‐point statistical facies simulation from uncertain training images , 2013 .

[20]  Christian P. Robert,et al.  The Bayesian choice : from decision-theoretic foundations to computational implementation , 2007 .

[21]  T. LaForce,et al.  Bayesian Reservoir History Matching Considering Model and Parameter Uncertainties , 2012, Mathematical Geosciences.

[22]  Geir Nævdal,et al.  Reservoir Monitoring and Continuous Model Updating Using Ensemble Kalman Filter , 2005 .

[23]  G. Schwarz.  Estimating the Dimension of a Model , 1978 .

[24]  Ionel M. Navon,et al.  Parameter estimation of subsurface flow models using iterative regularized ensemble Kalman filter , 2013, Stochastic Environmental Research and Risk Assessment.

[25]  Mary F. Wheeler,et al.  An iterative stochastic ensemble method for parameter estimation of subsurface flow models , 2013, J. Comput. Phys..

[26]  Andres Alcolea,et al.  Inverse problem in hydrogeology , 2005 .

[27]  D. Oliver,et al.  Markov chain Monte Carlo methods for conditioning a permeability field to pressure data , 1997 .

[28]  T. Hansen,et al.  Inverse problems with non-trivial priors: efficient solution through sequential Gibbs sampling , 2012, Computational Geosciences.

[29]  A. Reynolds,et al.  Uncertainty quantification of reservoir performance predictions using a stochastic optimization algorithm , 2011 .

[30]  Mary F. Wheeler,et al.  Clustered iterative stochastic ensemble method for multi-modal calibration of subsurface flow models , 2013 .

[31]  G. Mariéthoz,et al.  Analog‐based meandering channel simulation , 2014 .

[32]  Mary F. Wheeler,et al.  Hybrid nested sampling algorithm for Bayesian model selection applied to inverse subsurface flow problems , 2014, J. Comput. Phys..

[33]  J. Skilling.  Nested sampling for general Bayesian computation , 2006 .

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

[35]  Nicolas Remy.  S-GeMS: The Stanford Geostatistical Modeling Software: A Tool for New Algorithms Development , 2005 .

[36]  Mary F. Wheeler,et al.  Rapid updating of stochastic models by use of an ensemble-filter approach , 2014 .

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