An Iterative Ensemble Kalman Filter for Multiphase Fluid Flow Data Assimilation
- Yaqing Gu (U. of Oklahoma) | Dean S. Oliver (U. of Oklahoma)
- Document ID
- Society of Petroleum Engineers
- SPE Journal
- Publication Date
- November 2007
- Document Type
- Journal Paper
- 438 - 446
- 2007. Society of Petroleum Engineers
- 5.1 Reservoir Characterisation, 5.5 Reservoir Simulation, 5.1.1 Exploration, Development, Structural Geology, 1.6.9 Coring, Fishing, 5.5.8 History Matching, 3.3 Well & Reservoir Surveillance and Monitoring, 5.1.5 Geologic Modeling, 5.3.2 Multiphase Flow, 5.6.5 Tracers, 5.3.1 Flow in Porous Media, 7.2.2 Risk Management Systems
- 1 in the last 30 days
- 837 since 2007
- Show more detail
- View rights & permissions
|SPE Member Price:||USD 10.00|
|SPE Non-Member Price:||USD 30.00|
The dynamical equations for multiphase flow in porous media are highly nonlinear and the number of variables required to characterize the medium is usually large, often two or more variables per simulator gridblock. Neither the extended Kalman filter nor the ensemble Kalman filter is suitable for assimilating data or for characterizing uncertainty for this type of problem. Although the ensemble Kalman filter handles the nonlinear dynamics correctly during the forecast step, it sometimes fails badly in the analysis (or updating) of saturations.
This paper focuses on the use of an iterative ensemble Kalman filter for data assimilation in nonlinear problems, especially of the type related to multiphase .ow in porous media. Two issues are key: (1) iteration to enforce constraints and (2) ensuring that the resulting ensemble is representative of the conditional pdf (i.e., that the uncertainty quantification is correct). The new algorithm is compared to the ensemble Kalman filter on several highly nonlinear example problems, and shown to be superior in the prediction of uncertainty.
For linear problems, the Kalman filter is optimal for assimilating measurements to continuously update the estimate of state variables. Kalman filters have occasionally been applied to the problem of estimating values of petroleum reservoir variables (Eisenmann et al. 1994; Corser et al. 2000), but they are most appropriate when the problems are characterized by a small number of variables and when the variables to be estimated are linearly related to the observations. Most data assimilation problems in petroleum reservoir engineering are highly nonlinear and are characterized by many variables, often two or more variables per simulator gridblock.
The problem of weather forecasting is in many respects similar to the problem of predicting future petroleum reservoir performance. The economic impact of inaccurate predictions is substantial in both cases, as is the difficulty of assimilating very large data sets and updating very large numerical models. One method that has been recently developed for assimilating data in weather forecasting is ensemble Kalman filtering (Evensen 1994; Houtekamer and Mitchell 1998; Anderson and Anderson 1999; Hamill et al. 2000; Houtekamer and Mitchell 2001; Evensen 2003). It has been demonstrated to be useful for weather prediction over the North Atlantic. The method is now beginning to be applied for data assimilation in groundwater hydrology (Reichle et al. 2002; Chen and Zhang 2006) and in petroleum engineering (Nævdal et al. 2002, 2005; Gu and Oliver 2005; Liu and Oliver 2005a; Wen and Chen 2006, 2007; Zafari and Reynolds 2007; Gao et al. 2006; Lorentzen et al. 2005; Skjervheim et al. 2007; Dong et al. 2006), but the applications to state variables whose density functions are bimodal has proved problematic (Gu and Oliver 2006).
For applications to nonlinear assimilation problems, it is useful to think of the ensemble Kalman filter as a least squares method that obtains an averaged gradient for minimization not from a variational approach but from an empirical correlation between model variables (Anderson 2003; Zafari et al. 2006). In addition to providing a mean estimate of the variables, a Monte Carlo estimate of uncertainty can be obtained directly from the variability in the ensemble.
|File Size||770 KB||Number of Pages||9|
Anderson, J.L. 2001. An Ensemble Adjustment Kalman Filter for DataAssimilation. Monthly Weather Review 129 (12): 2884-2903.
Anderson, J.L. 2003. A Local Least Squares Framework for Ensemble Filtering.Monthly Weather Review 131 (4): 634-642.
Anderson, J.L. and Anderson, S.L. 1999. A Monte Carlo Implementation of theNonlinear Filtering Problem to Produce Ensemble Assimilations and Forecasts.Monthly Weather Review 127 (12): 2741-2758.
Bell, B.M. and Cathey, F.W. 1993. The iterated Kalman Filter Update as aGauss-Newton Method. IEEE Transactions on Automatic Control 38(2): 294-297.
Bishop, C.H., Etherton, B.J., and Majumdar, S.J. 2001. Adaptive SamplingWith the Ensemble Transform Kalman Filter. Part I: Theoretical Aspects. Mon.Wea. Rev. 129: 420-436.
Chen, Y. and Zhang, D. 2006. Data Assimilation for Transient Flow inGeologic Formations via Ensemble Kalman Filter. Advances in WaterResources 29 (8): 1107-1122.
Corser, G.P., Harmse, J.E., Corser, B.A., Weiss, M.W., and Whitlow, G.L.2000. Field Test Results for aReal-Time Intelligent Drilling Monitor. Paper SPE 59227 presented at theIADC/SPE Drilling Conference, New Orleans, 23-25 February. DOI:10.2118/59227-MS.
Dennis, J.E. and Schnabel, R.B. 1983. Numerical Methods for UnconstrainedOptimization and Nonlinear Equations. Prentice-Hall, Englewood Cliffs.
Dong, Y., Gu, Y., and Oliver, D.S. 2006. Sequential Assimilation of 4DSeismic Data for Reservoir Description Using the Ensemble Kalman Filter.Journal of Petroleum Science and Engineering 53 (1-2): 83-99.
Eisenmann, P., Gounot, M.-T., Juchereau, B., and Whittaker, S.J., 1994. Improved Rxo Measurements ThroughSemi-Active Focusing. Paper SPE 28437 presented at the SPE Annual TechnicalConference and Exhibition, New Orleans, 25-28 September. DOI:10.2118/28437-MS.
Evensen, G. 1994. Sequential Data Assimilation With a NonlinearQuasi-geostrophic Model Using Monte Carlo Methods to Forecast Error Statistics.Journal of Geophysical Research 99 (C5): 10143-10162.
Evensen, G. 2003. The Ensemble Kalman Filter: Theoretical Formulation andPractical Implementation. Ocean Dynamics 53: 343-367.
Evensen, G., Hove, J., Meisingset, H.C., Reiso, E., Seim, K.S., and Espelid,O. 2007. Using the EnKF forAssisted History Matching of a North Sea Reservoir Model . Paper SPE 106184presented at the SPE Reservoir Simulation Symposium, Houston, 26-28 February.DOI: 10.2118/106184-MS.
Gao, G., Zafari, M., and Reynolds, A.C. 2006. Quantifying Uncertainty for thePUNQ-S3 Problem in a Bayesian Setting with RML and EnKF. SPEJ11 (4): 506-515. SPE-93324-PA. DOI: 10.2118/93324-PA.
Gavalas, G.R., Shah, P.C., and Seinfeld, J.H. 1976. Reservoir History Matching by BayesianEstimation. SPEJ 16 (6): 337- 350. SPE-5740-PA. DOI:10.2118/5740-PA.
Golub, G.H. and van Loan, C.F. 1989. Matrix Computations. The JohnsHopkins University Press, Baltimore, second edition.
Gu, Y. 2006. History Matching Production Data Using the Ensemble KalmanFilter. PhD dissertation, University of Oklahoma.
Gu, Y. and Oliver, D.S. 2005. History Matching of the PUNQS3Reservoir Model Using the Ensemble Kalman Filter. SPEJ 10(2): 51-65. SPE-89942-PA. DOI: 10.2118/89942-PA.
Gu, Y. and Oliver, D.S. 2006. The Ensemble Kalman Filter for ContinuousUpdating of Reservoir Simulation Models. Journal of Energy ResourcesTechnology 128 (1): 79-87.
Hamill, T.M., Snyder, C., Baumhefner, D.P., Toth, Z., and Mullen, S.L. 2000.Ensemble Forecasting in the Short to Medium Range: Report From a Workshop.Bull. Amer. Meteor. Soc. 81: 2653-2664.
Houtekamer, P.L. and Mitchell, H.L. 1998. Data Assimilation Using anEnsemble Kalman Filter Technique. Monthly Weather Review 126 (3):796-811.
Houtekamer, P.L. and Mitchell, H.L. 2001. A Sequential Ensemble KalmanFilter for Atmospheric Data Assimilation. Monthly Weather Review129 (1): 123-137.
Jazwinski, A.H. 1970. Stochastic Processes and Filtering Theory.Academic Press, New York.
Journel, A. and Huijbregts, C.J. 1978. Mining Geostatistics. AcademicPress, New York. 600 p.
Kitanidis, P.K. 1995. Quasi-Linear Geostatistical Theory for Inversing.Water Resour. Res. 31 (10): 2411-2419.
Li, R., Reynolds, A.C., and Oliver, D.S. 2003. History Matching of Three-Phase FlowProduction Data. SPEJ 8 (4): 328-340. SPE-87336-PA. DOI:10.2118/87336-PA.
Liu, N. and Oliver, D.S. 2003. Evaluation of Monte Carlo Methods forAssessing Uncertainty. SPEJ 8 (2): 188 195. SPE-84936- PA.DOI: 10.2118/84936-PA.
Liu, N. and Oliver, D.S. 2005a. Critical Evaluation of the EnsembleKalman Filter on History Matching of Geologic Facies. SPEREE8 (6): 470-477. SPE-92867-PA. DOI: 10.2118/92867-PA.
Liu, N. and Oliver, D.S. 2005b. Ensemble Kalman Filter for Automatic HistoryMatching of Geologic Facies. Journal of Petroleum Science andEngineering 47 (3-4): 147-161.
Lorentzen, R.J., Nævdal, G., Vàlles, B., Berg, A.M., and A.-A. Grimstad2005. Analysis of the EnsembleKalman Filter for Estimation of Permeability and Porosity in ReservoirModels. Paper SPE 96375 presented at the SPE Annual Technical Conferenceand Exhibition, Dallas, 9-12 October. DOI: 10.2118/96375-MS.
Nævdal, G., Johnsen, L.M., Aanonsen, S.I., and Vefring, E.H. 2005. Reservoir Monitoring and ContinuousModel Updating Using Ensemble Kalman Filter. SPEJ 10 (1):66-74. SPE-84372-PA. DOI: 10.2118/84372-PA.
Nævdal, G., Mannseth, T., and Vefring, E.H. 2002. Near-Well Reservoir MonitoringThrough Ensemble Kalman Filter. Paper SPE 75235 presented at the SPE/DOEImproved Oil Recovery Symposium, Tulsa, 13-17 April. DOI: 10.2118/75235-MS.
Oliver, D.S., He, N., and Reynolds, A.C. 1996. Conditioning PermeabilityFields to Pressure Data. In European Conference for the Mathematics of OilRecovery, V 1-11.
Reichle, R.H., McLaughlin, D.B., and Entekhabi, D. 2002. Hydrologic DataAssimilation With the Ensemble Kalman Filter. Monthly Weather Review130 (1): 103-114.
Ripley, B.D. 1987. Stochastic Simulation. John Wiley & Sons, NewYork.
Skjervheim, J.-A., Evensen, G., Aanonsen, S.I., Ruud, B.O., and Johansen,T.A. 2007. Incorporating 4DSeismic Data in Reservoir Simulation Models Using Ensemble Kalman Filter.SPEJ 12 (3): 282-292. SPE-95789-PA. DOI: 10.2118/95789-PA.
Tan, T.B. and Kalogerakis, N. 1992. A Three-Dimensional Three-PhaseAutomatic History Matching Model: Reliability of Parameter Estimates.Journal of Canadian Petroleum Technology 31 (3): 34-41.
Tarantola, A. 1987. Inverse Problem Theory: Methods for Data Fitting andModel Parameter Estimation. Elsevier, Amsterdam, The Netherlands.
Thulin, K., Li, G., Aanonsen, S.I., and Reynolds, A.C. 2007. Estimation ofInitial Fluid Contacts by Assimilation of Production Data With EnKF. Paper tobe presented at the SPE Annual Technical Conference and Exhibition, Anaheim,California, 11-14 November.
Tjølsen, C.B., Damsleth, E., and Bu, T. 1994. The Effect of StochasticRelative Permeabilities in Reservoir Simulation. Journal of PetroleumScience and Engineering 10: 273-290.
Wen, X.-H. and Chen, W.H. 2006. Real-time Reservoir Model UpdatingUsing Ensemble Kalman Filter. SPEJ 11 (4): 431-442.SPE-92991-PA. DOI: 10.2118/92991-PA.
Wen, X.-H. and Chen, W.H. 2007. Some Practical Issues on RealtimeReservoir Model Updating Using Ensemble Kalman Filter. SPEJ12 (2): 156-166. SPE-111571-PA. DOI: 10.2118/111571-PA.
Zafari, M., Li, G., and Reynolds, A.C. 2006. Iterative Forms of the EnsembleKalman Filter. In Proceedings of the 10th European Conference on theMathematics of Oil Recovery—Amsterdam A030.
Zafari, M. and Reynolds, A.C. 2007. Assessing the Uncertainty inReservoir Description and Performance Predictions With the Ensemble KalmanFilter. SPEJ 12 (3): 378-387. SPE-95750-PA. DOI:10.2118/95750-PA.
Zupanski, M. 2005. Maximum Likelihood Ensemble Filter: Theoretical Aspects.Monthly Weather Review 133 (6): 1710-1726.