Integrating production data under uncertainty by parallel interacting Markov chains on a reduced dimensional space

In oil industry and subsurface hydrology, geostatistical models are often used to represent the porosity or the permeability field. In history matching of a geostatistical reservoir model, we attempt to find multiple realizations that are conditional to dynamic data and representative of the model uncertainty space. A relevant way to simulate the conditioned realizations is by generating Monte Carlo Markov chains (MCMC). The huge dimensions (number of parameters) of the model and the computational cost of each iteration are two important pitfalls for the use of MCMC. In practice, we have to stop the chain far before it has browsed the whole support of the posterior probability density function. Furthermore, as the relationship between the production data and the random field is highly nonlinear, the posterior can be strongly multimodal and the chain may stay stuck in one of the modes. In this work, we propose a methodology to enhance the sampling properties of classical single MCMC in history matching. We first show how to reduce the dimension of the problem by using a truncated Karhunen–Loève expansion of the random field of interest and assess the number of components to be kept. Then, we show how we can improve the mixing properties of MCMC, without increasing the global computational cost, by using parallel interacting Markov Chains. Finally, we show the encouraging results obtained when applying the method to a synthetic history matching case.

[1]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[2]  M. Loève Probability theory : foundations, random sequences , 1955 .

[3]  P. C. Shah,et al.  Reservoir History Matching by Bayesian Estimation , 1976 .

[4]  C. W. Harper,et al.  A FORTRAN IV program for comparing ranking algorithms in quantitative biostratigraphy , 1984 .

[5]  R. Ghanem,et al.  Stochastic Finite Elements: A Spectral Approach , 1990 .

[6]  C. Geyer Markov Chain Monte Carlo Maximum Likelihood , 1991 .

[7]  P. Kitanidis Quasi‐Linear Geostatistical Theory for Inversing , 1995 .

[8]  Dean S. Oliver,et al.  On conditional simulation to inaccurate data , 1996 .

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

[10]  Chao Yang,et al.  ARPACK users' guide - solution of large-scale eigenvalue problems with implicitly restarted Arnoldi methods , 1998, Software, environments, tools.

[11]  L. Hu Gradual Deformation and Iterative Calibration of Gaussian-Related Stochastic Models , 2000 .

[12]  C. Robert,et al.  Controlled MCMC for Optimal Sampling , 2001 .

[13]  Yukito Iba EXTENDED ENSEMBLE MONTE CARLO , 2001 .

[14]  H. Haario,et al.  An adaptive Metropolis algorithm , 2001 .

[15]  Dean S. Oliver,et al.  Integration of production data into reservoir models , 2001, Petroleum Geoscience.

[16]  Albert Tarantola,et al.  Inverse problem theory - and methods for model parameter estimation , 2004 .

[17]  Christian P. Robert,et al.  Monte Carlo Statistical Methods , 2005, Springer Texts in Statistics.

[18]  Michael W Deem,et al.  Parallel tempering: theory, applications, and new perspectives. , 2005, Physical chemistry chemical physics : PCCP.

[19]  Yan Chen,et al.  Data assimilation for transient flow in geologic formations via ensemble Kalman filter , 2006 .

[20]  L. Durlofsky,et al.  Efficient real-time reservoir management using adjoint-based optimal control and model updating , 2006 .

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

[22]  Yalchin Efendiev,et al.  Coarse-gradient Langevin algorithms for dynamic data integration and uncertainty quantification , 2006, J. Comput. Phys..

[23]  C. Andrieu,et al.  On the ergodicity properties of some adaptive MCMC algorithms , 2006, math/0610317.

[24]  S. Kou,et al.  Equi-energy sampler with applications in statistical inference and statistical mechanics , 2005, math/0507080.

[25]  Assessing the Dimensionality of Random Fields with Karhunen-Loève Expansion , 2007 .

[26]  P. Moral,et al.  NON-LINEAR MARKOV CHAIN MONTE CARLO , 2007 .