Critical Evaluation of the Ensemble Kalman Filter on History Matching of Geologic Facies
- Ning Liu (Chevron Corp.) | Dean S. Oliver (U. of Oklahoma)
- Document ID
- Society of Petroleum Engineers
- SPE Reservoir Evaluation & Engineering
- Publication Date
- December 2005
- Document Type
- Journal Paper
- 470 - 477
- 2005. Society of Petroleum Engineers
- 5.6.9 Production Forecasting, 5.2 Reservoir Fluid Dynamics, 1.2.3 Rock properties, 5.1 Reservoir Characterisation, 3.3 Well & Reservoir Surveillance and Monitoring, 5.1.5 Geologic Modeling, 5.3.1 Flow in Porous Media, 5.5 Reservoir Simulation, 5.5.8 History Matching
- 0 in the last 30 days
- 560 since 2007
- Show more detail
- View rights & permissions
|SPE Member Price:||USD 12.00|
|SPE Non-Member Price:||USD 35.00|
The objective of this paper is to compare the performance of the ensembleKalman filter (EnKF) to the performance of a gradient-based minimization methodfor the problem of estimation of facies boundaries in history matching. TheEnKF is a Monte Carlo method for data assimilation that uses an ensemble ofreservoir models to represent and update the covariance of variables. Inseveral published studies, it outperforms traditional history-matchingalgorithms in adaptability and efficiency.
Because of the approximate nature of the EnKF, the realizations from oneensemble tend to underestimate uncertainty, especially for problems that arehighly nonlinear. In this paper, the distributions of reservoir-modelrealizations from 20 independent ensembles are compared with the distributionsfrom 20 randomized-maximum-likelihood (RML) realizations for a 2D waterfloodmodel with one injector and four producers. RML is a gradient-based samplingmethod that generates one reservoir realization in each minimization of theobjective function. It is an approximate sampling method, but its samplingproperties are similar to the Markov-chain Monte Carlo (McMC) method on highlynonlinear problems and are relatively more efficient than McMC.
Despite the nonlinear relationship between the data (such as productionrates and facies observations) and the model variables, the EnKF was effectiveat history matching the production data. We find that the computational effortto generate 20 independent realizations was similar for the two methods,although the complexity of the code is substantially less for the EnKF.
Several questions regarding the use of the EnKF for history matching areaddressed in this paper. The most important is a comparison of the efficiencywith a gradient-based method for a history-matching problem with known faciesproperties but unknown boundary locations. Secondly, the EnKF and agradient-based method are unlikely to give identical estimates of modelvariables, so it is also important to know if one method generates betterrealizations. Finally, because there is often a desire to use thehistory-matched realizations to quantify uncertainty, it is important todetermine if one of the methods is more efficient at generating independentrealizations.
Gradient-based history matching can be performed in several ways (e.g.,assimilating data in batch or sequentially); a variety of minimizationalgorithms can be used (e.g., conjugate gradient or quasi-Newton); and severaldifferent methods for computing the gradient are available (e.g., adjoint orsensitivity equations). In this paper, we use what we believe is the mostefficient of the traditional gradient-based methods: an adjoint method tocompute the gradient of the squared data mismatch and the limited-memoryBroyden-Fletcher-Goldfarb-Shanno (LBFGS) method to compute the direction of thechange. The remaining choice is whether to incorporate all data at once orsequentially. Simultaneous, or batch, inversion of all data is clearly awell-established history-matching procedure. Although data from wells orsensors may arrive nearly continuously, the practice of updating reservoirmodels as the data arrive is not common. There are several reasons that makesequential assimilation of data difficult for large, nonlinear models: (1) thecovariance for all model variables must be updated as new data are assimilated,but the covariance matrix is very large; (2) the covariance may not be a goodmeasure of uncertainty for nonlinear problems; and (3) the sensitivity of adatum to changes in values of model variables is expensive to compute.
Bayesian updating in general is described by Woodbury. Modifying a methoddescribed by Tarantola, Oliver evaluated the possibility of using a sequentialassimilation approach for transient flow in porous media. He found that theresults from sequential assimilation could be almost as good as those frombatch assimilation if the order of the data was carefully selected. The problemwas quite small, however, and an extension to large models was impractical.Although a sequential method has the advantage of generating a sequence ofhistory-matched models that may all be useful at the time they are generated,our comparisons of efficiency will be based primarily on the effort required toassimilate all the data. If the intermediate predictions are needed (as theywould be for control of a reservoir), the comparison provided here willunderestimate the value of the sequential assimilation.
A secondary objective of history matching is often to assess the uncertaintyin the predictions of future reservoir performance or in the estimates ofreservoir properties such as permeability, porosity, or saturation. In general,uncertainty is estimated from an examination of a moderate number ofconditional simulations of the prediction or properties. Unless therealizations are generated fairly carefully and the sample is sufficientlylarge, however, the estimate of uncertainty could be quite poor. Two largecomparative studies of the ability of Monte Carlo methods to quantifyuncertainty in history matching have been carried out, one in groundwater andone in petroleum. Neither was conclusive, partly because of the small samplesize. Liu and Oliver used a smaller reservoir model (fewer variables), but amuch larger sample size. They found that the method that minimizes an objectivefunction containing a model mismatch part and a data mismatch part, with noiseadded to observations, created realizations that were distributed nearly thesame as realizations from McMC.
The EnKF is a Monte Carlo method for updating reservoir models. It solvesseveral problems with the application of the Kalman filter to large nonlinearproblems. It has been applied to reservoir flow problems with generally goodresults. There has been no examination, however, of the distribution of themembers of a single ensemble. The adequacy of the uncertainty estimate iscompletely unknown. In the first paper on the EnKF, Evensen described how theevolution of the probability density function for the model variables can beapproximated by the motion of "particles" or ensemble members in phase space.Any desired statistical quantities can be estimated from the ensemble ofpoints. When the size of the ensemble is relatively small, however, theapproximation of the covariance from the ensemble almost certainly containssubstantial errors. Houtekamer and Mitchell noted the tendency for a reductionin variance caused by "inbreeding." When the ensemble estimate is used in aKalman filter, van Leeuwen explained how nonlinearity in the covariance updaterelation causes growth in the error as additional data are assimilated. In thispaper, the comparison is made using history matching on a truncatedpluri-gaussian model for geologic facies. It provides a difficulthistory-matching problem with significant nonlinearities that make both theEnKF and the LBFGS method difficult to apply.
|File Size||944 KB||Number of Pages||8|
1. Chen, W.H. et al.: "A NewAlgorithm for Automatic History Matching," SPEJ (December 1974) 593;Trans., AIME, 257.
2. Chavent, G.M., Dupuy, M., and Lemonnier, P.: "History Matching by Use of OptimalTheory," SPEJ (February 1975) 74l; Trans., AIME, 259.
3. Makhlouf, E.M. et al.: "AGeneral History Matching Algorithm for Three-Phase, Three-Dimensional PetroleumReservoirs," SPE Advanced Technology Series (July 1993) 1, No. 2, 83.
4. Li, R., Reynolds, A.C., and Oliver, D.S.: "Sensitivity Coefficients forThree-Phase Flow History Matching," J. Cdn. Pet. Tech. (2003) 42, No. 4,70.
5. Nocedal, J.: "Updating Quasi-Newton Matrices with Limited Storage," Math.Comp. (1980) 35, 773.
6. Fletcher, R.: Practical Methods of Optimization, second edition,John Wiley & Sons, New York City (1987).
7. Deschamps, T. et al.: "The Results of Testing Six Different GradientOptimisers on Two History Matching Problems," Proc., 6th European Conference onthe Mathematics of Oil Recovery, Edinburgh, Scotland (1998).
8. Zhang, F. and Reynolds, A.C.: "Optimization Algorithms for AutomaticHistory Matching of Production Data," Proc., 8th European Conference on theMathematics of Oil Recovery, Freiberg, Germany (2002).
9. Woodbury, A.D.: "Bayesian updating revisited," Mathematical Geology(1989) 21, 285.
10. Tarantola, A.: Inverse Problem Theory: Methods for Data Fitting andModel Parameter Estimation, Elsevier, Amsterdam (1987).
11. Oliver, D.S.: "Incorporation of Transient Pressure Data into ReservoirCharacterization," In Situ (1994) 18, No. 3, 243.
12. Zimmerman, D.A. et al.: "A Comparison of Seven Geostatistically BasedInverse Approaches to Estimate Transmissivities for Modeling AdvectiveTransport by Groundwater Flow," WRR (1998) 34, No. 6, 1373.
13. Floris, F.J.T. et al.: "Methods for Quantifying the Uncertainty ofProduction Forecasts: A Comparative Study," Petroleum Geoscience (2001) 7,SUPP, 87.
14. Liu, N. and Oliver, D.S.: "Evaluation of Monte Carlo Methods forAssessing Uncertainty," SPEJ (June 2003) 188.
15. Evensen, G.: "Sequential Data Assimilation With A NonlinearQuasi-Geostrophic Model Using Monte Carlo Methods To Forecast ErrorStatistics," J. of Geophysical Research (1994) 99, Issue C5, 10143.
16. Anderson, J.L. and Anderson, S.L.: "A Monte Carlo Implementation of theNonlinear Filtering Problem to Produce Ensemble Assimilations and Forecasts,"Monthly Weather Review (1999) 127, No. 12, 2741.
17. Nævdal, G., Mannseth, T., and Vefring, E.H.: "Near-Well Reservoir MonitoringThrough Ensemble Kalman Filter," paper SPE 75235 presented at the 2002SPE/DOE Improved Oil Recovery Symposium, Tulsa, 13-17 April.
18. Nævdal, G. et al.: "Reservoir Monitoring and ContinuousModel Updating Using Ensemble Kalman Filter ," SPEJ (March 2005) 66.
19. Gu, Y. and Oliver, D.S.: "The Ensemble Kalman Filter for ContinuousUpdating of Reservoir Simulation Models," J. of Energy Resources Technology(March 2006).
20. Houtekamer, P.L. and Mitchell, H.L.: "Data Assimilation Using anEnsemble Kalman Filter Technique," Monthly Weather Review (1998) 126, No. 3,796.
21. van Leeuwen, P.J.: "Comment on ‘Data Assimilation Using an EnsembleKalman Filter Technique'," Monthly Weather Review (1999) 127, No. 6, 1374.
22. Galli, A. et al.: "The Pros and Cons of the Truncated Gaussian Method,"Geostatistical Simulations, Kluwer Academic, Dordrecht, The Netherlands (1994)217-233.
23. Lantuéjoul, C.: Geostatistical Simulation: Models and Algorithms, Springer, Berlin (2002).
24. Liu, N. and Oliver, D.S.: "Ensemble Kalman Filter for Automatic HistoryMatching of Geologic Facies," J. Pet. Sci. and Eng. (2005) 47, No. 4, 147.
25. Liu, N. and Oliver, D.S.: "Conditional Simulation of Truncated RandomFields Using Gradient Methods," Proc., 2003 IAMG Annual Meeting, Portsmouth,U.K. (2003).
26. Li, R., Reynolds, A.C., and Oliver, D.S.: "History Matching of Three-Phase FlowProduction Data," SPEJ (December 2003) 328.