Streamline-Assisted Ensemble Kalman Filter for Rapid and Continuous Reservoir Model Updating
- Elkin Arroyo (Texas A&M University) | Deepak Devegowda (U. of Oklahoma) | Akhil Datta-Gupta (Texas A&M University) | Jonggeun Choe (Seoul Natl. U.)
- Document ID
- Society of Petroleum Engineers
- SPE Reservoir Evaluation & Engineering
- Publication Date
- December 2008
- Document Type
- Journal Paper
- 1,046 - 1,060
- 2008. Society of Petroleum Engineers
- 4.3.4 Scale, 3.3 Well & Reservoir Surveillance and Monitoring, 5.1.2 Faults and Fracture Characterisation, 5.4.1 Waterflooding, 5.5.7 Streamline Simulation, 7.2.2 Risk Management Systems, 5.1 Reservoir Characterisation, 7.6.2 Data Integration, 5.1.5 Geologic Modeling, 5.5 Reservoir Simulation, 2.3 Completion Monitoring Systems/Intelligent Wells, 5.5.8 History Matching, 5.6.4 Drillstem/Well Testing, 5.3.1 Flow in Porous Media, 5.3.2 Multiphase Flow, 5.1.9 Four-Dimensional and Four-Component Seismic, 5.6.9 Production Forecasting
- 3 in the last 30 days
- 925 since 2007
- Show more detail
- View rights & permissions
|SPE Member Price:||USD 12.00|
|SPE Non-Member Price:||USD 35.00|
The use of the ensemble Kalman filter (EnKF) is a promising approach for data assimilation and assessment of uncertainties during reservoir characterization and performance forecasting. It provides a relatively straightforward approach to incorporating diverse data types, including production and/or time-lapse seismic data. Unlike traditional sensitivity-based history matching methods, the EnKF relies on a cross-covariance matrix computed from an ensemble of reservoir models to relate reservoir properties to production data. For practical field applications, we need to keep the ensemble size small for computational efficiency. However, this leads to poor approximations of the cross-covariance and, often, loss of geologic realism through parameter overshoots, in particular by introducing localized patches of low and high permeabilities. Because the EnKF estimates are "optimal" only for Gaussian variables and linear dynamics, these difficulties are compounded by the strong nonlinearity of the multiphase history matching problems and for non-Gaussian prior models. Specifically, the updated parameter distribution tends to become multi-Gaussian with loss of connectivities of extreme values, such as high permeability channels and low permeability barriers, which are of special significance during reservoir characterization.
We propose a novel approach to overcome some of these limitations by conditioning the cross-covariance matrix using information gleaned from streamline trajectories. Our streamline-assisted EnKF is analogous to the conventional assisted history matching, whereby the streamline trajectories are used to identify gridblocks contributing to the production response of a specific well. We then use these gridblocks only to compute the cross-covariance matrix and eliminate the influence of unrelated or distant observations and spurious correlations. We show that the streamline-assisted EnKF is an efficient and robust approach for history matching and continuous reservoir model updating. We illustrate the power and utility of our approach using both synthetic and field applications.
Proper characterization of the reservoir and the assessment of uncertainty are crucial aspects of any optimal reservoir development plan and management strategy. To achieve this goal, it is necessary to reconcile geological models to the dynamic response of the reservoir through history matching. The topic of history matching has been of great interest and an area of active research in the oil industry (Datta-Gupta and King 2007; Emanuel and Milliken 1998; Oliver et al. 2001). The past decade has seen some significant developments in assisted and automatic history matching of high-resolution reservoir models and associated uncertainty quantification. Many of these techniques involve computation of sensitivities that relate changes in production response at a well to a change in reservoir parameters. Techniques of automatic history matching that typically do not use parameter sensitivities or gradient of the misfit function are stochastic algorithms such as Markov Chain Monte Carlo (MCMC), simulated annealing and genetic algorithms (Ma et al. 2008; Sen et al. 2005). A relatively recent and promising addition to this class of techniques is the use of ensemble Kalman Filters (EnKF) for data assimilation (Gu and Oliver 2005, 2006; Naevdal et al. 2005; Gao et al. 2006; Skjervheim et al. 2007; Dong et al. 2006). It is a Monte-Carlo approach that works with an ensemble of reservoir models. Specifically, the method utilizes cross-covariances between measurements and model parameters computed directly from the ensemble members to sequentially update the reservoir models.
A major advantage of the EnKF is that it can be readily linked to any existing reservoir simulator. The ability to assimilate diverse data types and the ease of implementation have resulted in considerable interest in the approach. Moreover, EnKF uses a sequential updating technique; that is, the reservoir data is assimilated as and when it becomes available. The EnKF can assimilate the latest production data without re-running the simulator from the initial conditions. These characteristics make it particularly well-suited for continuous model updating. The increased application of downhole monitors, intelligent well systems, and permanent sensors to continuously record pressure, well rates, and temperature has provided a further boost to the sequential model updating through EnKF.
In spite of all its favorable properties, the current implementation of EnKF approach comes with its own share of challenges. A key requirement in history matching is that the final model should honor the available geological information and retain geologic realism. It has been shown that the EnKF works well when the prior distribution of parameters is Gaussian; however, the estimates are suboptimal for non-Gaussian distributions. Over a sequence of many updates, multimodal permeability distributions tend to transform to Gaussian distribution. During geologic model updating, this can lead to a loss of structure and connectivity of the extremes in the permeability field. This has serious implications in the fluid flow because of the influence of high-permeability channels and low-permeability barriers. Although there are some variants of the Kalman filter that work with non-Gaussian distributions, such as the Gaussian summation approximation, the implementation on an ensemble framework tend to be very expensive (Anderson and Moore 1979).
In the past few years, we have seen several applications of the EnKF for field-scale history matching, including some recent papers that attempt to deal with some of the challenges pertaining to its use (Gu and Oliver 2005, 2006; Naevdal et al. 2005; Gao et al. 2006; Skjervheim et al. 2007; Dong et al. 2006). In particular, localized overshooting of permeabilities has been reported, resulting in loss of geologic continuity. This is aggravated by the strong non-linearity inherent in multiphase flow simulations.
Another common difficulty experienced when using the EnKF is filter divergence. The effect of filter divergence is such that the distribution produced by the filter drifts away from the truth. Filter divergence normally occurs because the prior probability distribution becomes too narrow (loss of variance) and the observations have progressively less impact on the model updates.
One common approach to deal with filter divergence is to add some (white) noise to the prior ensemble to "inflate" its distribution and enhance the impact of new observations. Other problems and limitations of the EnKF, particularly for nonlinear problems and non-Gaussian parameter distributions, can be partly controlled using a large ensemble. However, for practical field applications, the ensemble size needs to be kept relatively small for computational efficiency.
This paper describes an approach to address many of the currently reported difficulties in the use of the EnKF applied to reservoir history matching. The unique feature of our proposed approach is that the final models that constitute the ensemble tend to retain the geological information that went into building them initially. Over a sequence of many updates, our approach tends to preserve the shape of the initial permeability distribution and consequently retains key geological features. Our approach greatly decreases the severity of the overshooting problem reported in earlier implementations of the EnKF. Moreover, it allows the use of smaller ensemble size, while providing results comparable or better than the standard EnKF.
The paper is organized as follows. First, we briefly review the major steps of the EnKF and the additional streamline-based conditioning of the cross-covariance proposed here. We also illustrate these steps using a synthetic example. Next, we discuss the underlying mathematical formulation in detail. We then demonstrate the power and practical utility of the approach using the benchmark PUNQ-S3 synthetic example (Gu and Oliver 2005) and a field example. Finally, an analysis of the scalability and speed-up factor for the parallel implementation of our code is given.
|File Size||19 MB||Number of Pages||15|
Anderson, B.D. and Moore, J.B. 1979. Optimal Filtering. EnglewoodCliffs, New Jersey: Prentice-Hall.
Barker, J.W., Cuypers, M., and Holden, L. 2001. Quantifying Uncertainty in ProductionForecasts: Another Look at the PUNQ-S3 Problem. SPEJ 6 (4):433-441. SPE-74707-PA. DOI: 10.2118/74707-PA.
Carter, J. PUNQ-S3 Test Case. http://www3.imperial.ac.uk/earthscienceandengineering/research/perm/punq-s3model/.
Datta-Gupta, A. and King, M.J. 2007. Streamline Simulation: Theory andPractice. Textbook Series, SPE, Richardson, Texas 11.
Deutsch, C.V. and Journel, A. 1992. GSLIB Geostatistical Software Libraryand User's Guide. Oxford, UK: Oxford University Press.
Devegowda, D. 2008. Streamline Assisted Ensemble Kalman Filter--Formulationand Field Application. PhD dissertation, Texas A&M University.
Dong, Y., Gu, Y., and Oliver, D.S. 2006. Sequential assimilationof 4D seismic data for reservoir description using ensemble Kalman filter.Journal of Petroleum Science and Engineering 53 (1-2): 83-99.DOI: 10.1016/j.petrol.2006.03.028.
Emanuel, A.S. and Milliken, W.J. 1998. History Matching Finite DifferenceModels With 3D Streamlines. Paper SPE 49000 presented at the SPE AnnualTechnical Conference and Exhibition, New Orleans, 27-30 September. DOI:10.2118/49000-MS.
Evensen, G. 1994. Sequential data assimilation with a nonlinearquasi-geostrophic model using Monte Carlo methods to forecast error statistics.J. of Geophysical Research 99 (C5): 10143-10162.
Evensen, G. 2003. TheEnsemble Kalman Filter: Theoretical formulation and practicalimplementation. Ocean Dynamics 53 (4): 343-367. DOI:10.1007/s10236-003-0036-9.
Floris, F.J.T., Bush, M.D., Cuypers, M., Roggero, F., and Syversveen, A-R.2001. Methods for quantifying the uncertainty of production forecasts: Acomparative study. Petroleum Geoscience 7 (Supplement, 1 May):87-96.
FrontSIM Reference Manual. 2005. Houston: Schlumberger InformationSystems.
Gao, G., Zafari, M., and Reynolds, A.C. 2006. Quantifying Uncertainty for thePUNQ-S3 Problem in a Bayesian Setting With RML EnKF. SPEJ 11(4): 506-515. SPE-93324-PA. DOI: 10.2118/93324-PA.
Gaspari, G. and Cohn, S.E. 1999. Construction of correlationfunctions in two and three dimensions (DAO Office Note 96-03R1). Quart.J. Roy. Meteor. Soc. 125 (554): 723-757. DOI:10.1002/qj.49712555417.
Gu, Y. and Oliver, D.S. 2005. History Matching of the PUNQ-S3Reservoir Model Using the Ensemble Kalman Filter. SPEJ 10(2): 217-224. SPE-89942-PA. DOI: 10.2118/89942-PA.
Gu, Y. and Oliver, D.S. 2006. The ensemble Kalman filter forcontinuous updating of reservoir simulation models. J. of EnergyResources Technology 128 (1): 79-87. DOI: 10.1115/1.2134735.
Hamill, T.M., Whitaker, J.S., and Snyder, C. 2001. Distance-dependent filtering of background error covariance estimate in anensemble Kalman filter. Monthly Weather Review 129(11): 2776-2790. DOI:10.1175/1520-0493(2001)129<2776:DDFOBE>2.0.CO;2.
He, Z. and Datta-Gupta, A. 2001. Streamline-Based Production DataIntegration Under Changing Field Conditions. Paper SPE 71333 presented atthe SPE Annual Technical Conference and Exhibition, New Orleans, 30 September-3October. DOI: 10.2118/71333-MS.
Houtekamer, P.L. and Mitchell, H.L. 2001. A sequential ensemble Kalman filter for atmospheric data assimilation.Monthly Weather Review 129 (1): 123-137. DOI:10.1175/1520-0493(2001)129<0123:ASEKFF>2.0.CO;2.
Jimenez, E., Sabir, K., Datta-Gupta, A., and King, M.J. 2007. Spatial Error and Convergence inStreamline Simulation. SPEREE 10 (3): 221-232. SPE-92873-PA.DOI: 10.2118/92873-PA.
Keppenne, C.L. 2000. Data assimilation into a primitive-equation model with a parallel ensembleKalman filter. Monthly Weather Review 128 (6): 1971-1981.DOI: 10.1175/1520-0493(2000)128<1971:DAIAPE>2.0.CO;2.
Ma, X., Al-Harbi, M., Datta-Gupta, A., and Efendiev, Y. 2008. An Efficient Two-Stage SamplingMethod for Uncertainty Quantification in History Matching GeologicalModels. SPEJ 13 (10): 77-87. SPE-102476-PA. DOI:10.2118/102476-PA.
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.
Oliver, D.S., Reynolds, A.C., Bi, Z., and Abacioglu, Y. 2001. Integration ofproduction data into reservoir models. Petroleum Geoscience 7(Supplement, 1 May): 65-73.
Sen, M.K., Datta-Gupta, A., Stoffa, P.L., Lake, L.W., and Pope, G.A. 1995.Stochastic Reservoir ModelingUsing Simulated Annealing and Genetic Algorithms. SPEFE 10(1): 49-55. SPE-24754-PA. DOI: 10.2118/24754-PA.
Skjervheim, J.-A., Evensen, G., Aanonsen, S.I., and Johansen, T.A. 2007. Incorporating 4D Seismic Data inReservoir Simulation Model Using Ensemble Kalman Filter. SPEJ12 (3): 282-292. SPE-95789-PA. DOI: 10.2118/95789-PA.
Vasco, D.W., Peterson, J.E. and Majer, E.L. 1996. Nonuniqueness in Travel-TimeTomography: Ensemble Inference and Cluster Analysis. Geophysics61(4):1209-1227. DOI:10.1190/1.1444040.