Prediction under uncertainty in reservoir modeling

Reservoir simulation is routinely employed in the prediction of reservoir performance under different depletion and operating scenarios. Usually, a single history-matched model, conditioned to production data, is obtained. The model is then used to forecast future production profiles. Because the history match is non-unique, the forecast production profiles are therefore uncertain, although this uncertainty is not usually quantified. This paper presents a new approach for generating uncertain reservoir performance predictions and quantifying the uncertainty associated with forecasting future performance. Firstly, multiple reservoir realizations are generated using a new stochastic algorithm. This involves adaptively sampling the model parameter space using an algorithm, which biases the sampling towards regions of good fit. Using the complete ensemble of models generated, the posterior distribution is resampled in order to quantify the uncertainty associated with forecasting reservoir performance in a Bayesian framework. The strength of the method in performance prediction is demonstrated by using an upscaled model to history match fine scale data. The maximum likelihood model is then used in forecasting the fine grid performance, and the uncertainty associated with the predictions is quantified. It is demonstrated that the maximum likelihood model is highly accurate in reservoir performance prediction. D 2004 Published by Elsevier B.V.

[1]  A. Tarantola Inverse problem theory : methods for data fitting and model parameter estimation , 1987 .

[2]  S Hou,et al.  Prediction of Oil Production With Confidence Intervals , 2001 .

[3]  Michael Andrew Christie,et al.  Uncertainty reduction in reservoir modelling , 2002 .

[4]  Michael Andrew Christie,et al.  Tenth SPE Comparative Solution Project: a comparison of upscaling techniques , 2001 .

[5]  G. L. Chierici,et al.  Novel relations for drainage and imbibition relative permeabilities , 1984 .

[6]  M. Sambridge Geophysical inversion with a neighbourhood algorithm—I. Searching a parameter space , 1999 .

[7]  Michael Andrew Christie,et al.  Multiple history-matched models for Teal South , 2002 .

[8]  Franz Aurenhammer,et al.  Voronoi diagrams—a survey of a fundamental geometric data structure , 1991, CSUR.

[9]  Jonathan Carter,et al.  Improved Reservoir Characterization through Evolutionary Computation , 2000 .

[10]  R.C.M. Portellaand,et al.  Use of Automatic History Matching and Geostatistical Simulation to Improve Production Forecast , 1999 .

[11]  J. Glimm,et al.  A Probability Model For Errors in the Numerical Solutions of a Partial Differential Equation , 1999 .

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

[13]  D. T. Rian,et al.  Optimization Methods for History Matching of Complex Reservoirs , 2001 .

[14]  Kenneth M. Hanson,et al.  Uncertainty assessment for reconstructions based on deformable geometry , 1997, Int. J. Imaging Syst. Technol..

[15]  F. Roggero Direct Selection of Stochastic Model Realizations Constrained to Historical Data , 1997 .

[16]  M. Sambridge Geophysical inversion with a neighbourhood algorithm—II. Appraising the ensemble , 1999 .