In most hydrocarbon reservoir development projects, geological models are fully rebuilt on a regular basis to integrate new data, in particular observations from new wells, for up-to-date forecasts. Not only this common practice is very time consuming as rebuilding models can take weeks or even months, but it also leads to major, hard-to-justify, fluctuations in reservoir volume or flow performance forecasts, especially when the modeling staff changes or a new modeling technology, workflow, or software is adopted. Rationalizing the geological model updating process is required to provide stable and reliable forecasting and make timely, well-informed, reservoir management decisions. This paper presents an innovative methodology to quickly reassess model forecasts, such as reservoir oil-in-place or oil recovery, without rebuilding any geological models provided that the new data observations are reasonably consistent with the current models. The proposed methodology uses a Bayesian framework whereby the multivariate probability joint distribution of new data predictions and forecast variables needs to be modeled. Assuming that this joint distribution is multi-Gaussian, the first step consists in computing proxies, e.g., response surfaces using experimental design, to estimate from the set of current geological models the distribution (mean and variance) of new data predictions and forecast variables as a function of the input modeling parameters (e.g., property variograms or training images, trends, histograms). Because the model stochasticity (i.e., spatial uncertainty away from wells) typically entails significant uncertainty in the prediction of new local data observations, computing the previous proxies requires generating multiple stochastic realizations for each combination of input modeling parameters. Then, using those proxies and Monte Carlo simulation, the full multivariate probability joint distribution of new data predictions and forecast variables is estimated. Plugging the actual new data values into that joint distribution finally provides new updated probabilistic distributions of the forecast variables. This new methodology is illustrated on a synthetic case study. In addition to quickly reassess reservoir volume and flow performance predictions, this new approach can be used to select new data observation types and impact maps to assess potential well locations that would optimally reduce forecasting uncertainties.
[1]
Daniel J. Mundfrom,et al.
A Monte Carlo comparison of the Type I and Type II error rates of tests of multivariate normality
,
2005
.
[2]
Clayton V. Deutsch,et al.
Impact Map for Assessment of New Delineation-Well Locations
,
2013
.
[3]
Jef Caers,et al.
Updating joint uncertainty in trend and depositional scenario for reservoir exploration and early appraisal
,
2015,
Computational Geosciences.
[4]
P. Mahalanobis.
On the generalized distance in statistics
,
1936
.
[5]
Jef Caers,et al.
Direct forecasting of subsurface flow response from non-linear dynamic data by linear least-squares in canonical functional principal component space
,
2014
.
[6]
Philippe Renard,et al.
Prediction-Focused Subsurface Modeling: Investigating the Need for Accuracy in Flow-Based Inverse Modeling
,
2015,
Mathematical Geosciences.
[7]
A. Boucher,et al.
History matching and uncertainty quantification of facies models with multiple geological interpretations
,
2013,
Computational Geosciences.
[8]
Sebastien Strebelle,et al.
Conditional Simulation of Complex Geological Structures Using Multiple-Point Statistics
,
2002
.
[9]
Anthony C. Atkinson,et al.
Optimum Experimental Designs, with SAS
,
2007
.