×

Investigation of the sampling performance of ensemble-based methods with a simple reservoir model. (English) Zbl 1382.86018

Summary: The application of the ensemble Kalman filter (EnKF) for history matching petroleum reservoir models has been the subject of intense investigation during the past 10 years. Unfortunately, EnKF often fails to provide reasonable data matches for highly nonlinear problems. This fact motivated the development of several iterative ensemble-based methods in the last few years. However, there exists no study comparing the performance of these methods in the literature, especially in terms of their ability to quantify uncertainty correctly. In this paper, we compare the performance of nine ensemble-based methods in terms of the quality of the data matches, quantification of uncertainty, and computational cost. For this purpose, we use a small but highly nonlinear reservoir model so that we can generate the reference posterior distribution of reservoir properties using a very long chain generated by a Markov chain Monte Carlo sampling algorithm. We also consider one adjoint-based implementation of the randomized maximum likelihood method in the comparisons.

MSC:

86A32 Geostatistics
86-08 Computational methods for problems pertaining to geophysics
76S05 Flows in porous media; filtration; seepage

Software:

L-BFGS
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Aanonsen, S.I., Nævdal, G., Oliver, D.S., Reynolds, A.C., Vallés, B.: Review of ensemble Kalman filter in petroleum engineering. Soc. Pet. Eng. J. 14(3), 393-412 (2009)
[2] Anderson, J.L., Anderson, S.L.: A Monte Carlo implementation of the nonlinear filtering problem to produce ensemble assimilations and forecasts. Mon. Weather Rev. 127(12), 2741-2758 (1999)
[3] Bonet-Cunha, L., Oliver, D.S., Rednar, R.A., Reynolds, A.C.: A hybrid Markov chain Monte Carlo method for generating permeability fields conditioned to multiwell pressure data and prior information. Soc. Pet. Eng. J. 11(3), 261-271 (1998)
[4] Burgers, G., van Leeuwen, P., Evensen, G.: Analysis scheme in the ensemble Kalman filter. Mon. Weather Rev. 126(6), 1719-1724 (1998)
[5] Chavent, G.M., Dupuy, M., Lemonnier, P.: History matching by use of optimal control theory. Soc. Pet. Eng. J. 15(1), 74-86 (1975)
[6] Chen, W.H., Gavalas, G.R., Seinfeld, J.H., Wasserman, M.L.: A new algorithm for automatic history matching. Soc. Pet. Eng. J. 14(6), 593-608 (1974)
[7] Chen, Y., Oliver, D.: Ensemble randomized maximum likelihood method as an iterative ensemble smoother. Math. Geosci. 44(1), 1-26 (2011)
[8] Marsily, G.; Lavedan, G.; Boucher, M.; Fasanino, G.; Verly, G. (ed.); David, M. (ed.); Journel, AG (ed.); Marechal, A. (ed.), Interpretation of interference tests in a well field using geostatistical techniques to fit the permeability distribution in a reservoir model, 831-849 (1984), Dordrecht
[9] Emerick, A.A., Reynolds, A.C.: Combining the ensemble Kalman filter with Markov chain Monte Carlo for improved history matching and uncertainty characterization. Soc. Pet. Eng. J. 17(2), 418-440 (2012a)
[10] Emerick, A.A., Reynolds, A.C.: Ensemble smoother with multiple data assimilation. Computers & Geosciences (2012b). doi:10.1016/j.cageo.2012.03.011
[11] Emerick, A.A., Reynolds, A.C.: History matching time-lapse seismic data using the ensemble Kalman filter with multiple data assimilations. Comput. Geosci 16(3), 639-659 (2012c)
[12] Evensen, G.: Sequential data assimilation with a nonlinear quasi-geostrophic model using Monte Carlo methods to forecast error statistics. J. Geophys. Res. 99(C5), 10,143-10,162 (1994)
[13] Evensen, G.: Sampling strategies and square root analysis schemes for the EnKF. Ocean Dyn. 54(6), 539-560 (2004)
[14] Evensen, G.: The ensemble Kalman filter for combined state and parameter estimation. In: IEEE Control Systems Magazine, pp. 83-104 (2009) · Zbl 1395.93534
[15] Gao, G., Reynolds, A.C.: An improved implementation of the LBFGS algorithm for automatic history matching. Soc. Pet. Eng. J. 11(1), 5-17 (2006)
[16] Gao, G., Zafari, M., Reynolds, A.C.: Quantifying uncertainty for the PUNQ-S3 problem in a Bayesian setting with RML and EnKF. Soc. Pet. Eng. J. 11(4), 506-515 (2006)
[17] Gelman, A., Roberts, G.O., Gilks, W.R.: Efficient Metropolis jumping rules. Bayesian Stat. Oxford University Press 5(5), 599-607 (1996)
[18] Gu, Y., Oliver, D.S.: An iterative ensemble Kalman filter for multiphase fluid flow data assimilation. Soc. Pet. Eng. J. 12(4), 438-446 (2007)
[19] Hastings, W.K.: Monte Carlo sampling methods using Markov chains and their applications. Biometrika 57(1), 97-109 (1970) · Zbl 0219.65008
[20] Houtekamer, P.L., Mitchell, H.L.: Data assimilation using an ensemble Kalman filter technique. Mon. Weather Rev. 126(3), 796-811 (1998)
[21] Houtekamer, P.L., Mitchell, H.L.: A sequential ensemble Kalman filter for atmospheric data assimilation. Mon. Weather Rev. 129(1), 123-137 (2001)
[22] Hu, L.Y.: Gradual deformation and iterative calibration of Gaussian-related stochastic models. Math. Geol. 32(1), 87-108 (2000)
[23] Kitanidis, P.K.: Quasi-linear geostatistical theory for inversing. Water Resour. Res. 31(10), 2411-2419 (1995)
[24] Krymskaya, M.V., Hanea, R.G., Verlaan, M.: An iterative ensemble Kalman filter for reservoir engineering applications. Comput. Geosci 13(2), 235-244 (2009) · Zbl 1171.86306
[25] LaVenue, A.M., Pickens, J.F.: Application of a coupled adjoint sensitivity and kriging approach to calibrate a groundwater flow model. Water Resour. Res. 28(6), 1543-1569 (1992)
[26] LaVenue, A.M., RamaRao, B.S., de Marsily, G., Marietta, M.G.: Pilot point methodology for automated calibration of an ensemble of conditionally simulated transmissivity fields, 2. Application, Water Resour. Res. 31(3), 495-516 (1995)
[27] Li, G., Reynolds, A.C.: Iterative ensemble Kalman filters for data assimilation. Soc. Pet. Eng. J. 14(3), 496-505 (2009)
[28] Li, R., Reynolds, A.C., Oliver, D.S.: History matching of three-phase flow production data. Soc. Pet. Eng. J. 8(4), 328-340 (2003)
[29] Liu, D., Nocedal, J.: On the limited memory BFGS method for large scale optimization. Math. Program. 45, 503-528 (1989) · Zbl 0696.90048
[30] Liu, N.: Assessment of Uncertainty Assessment Methods. Master’s thesis, The University of Tulsa, Tulsa, Oklahoma (2002)
[31] Liu, N., Oliver, D.S.: Evaluation of Monte Carlo methods for assessing uncertainty. Soc. Pet. Eng. J. 8(2), 188-195 (2003)
[32] Lorentzen, R.J., Nævdal, G.: An iterative ensemble Kalman filter. IEEE Trans. Automat. Contr. 56(8), 1990-1995 (2011) · Zbl 1368.93721
[33] Lorentzen, R.J., Nævdal, G., Vallès, B., Berg, A.M., Grimstad, A.-A.: Analysis of the ensemble Kalman filter for estimation of permeability and porosity in reservoir models. In: Proceedings of the SPE Annual Technical Conference and Exhibition, Dallas, TX, SPE 96375, 9-12 October 2005
[34] Metropolis, N., Rosenbluth, A.W., Rosenbluth, M.N., Teller, A.H., Teller, E.: Equations of state calculations by fast computing machines. J. Chem. Phys. 21, 1087-1092 (1953) · Zbl 1431.65006
[35] Moradkhani, H., Sorooshian, S., Gupta, H.V., Mouser, P.R.: Dual state-parameter estimation of hydrological models using ensemble Kalman filter. Adv. Water Resour. 28(2), 301-305 (2005)
[36] Nocedal, J.: Updating quasi-Newton matrices with limited storage. Math. Comput. 35(151), 773-782 (1980) · Zbl 0464.65037
[37] Oliver, D.S., Cunha, L.B., Reynolds, A.C.: Markov chain Monte Carlo methods for conditioning a permeability field to pressure data. Math. Geol. 29(1), 61-91 (1997)
[38] Oliver, D.S., He, N., Reynolds, A.C.: Conditioning permeability fields to pressure data. In: Proceedings of the European Conference for the Mathematics of Oil Recovery (1996)
[39] Oliver, D.S., Reynolds, A.C., Liu, N.: Inverse Theory for Petroleum Reservoir Characterization and History Matching. Cambridge University Press, Cambridge (2008)
[40] Peaceman, D.W.: Interpretation of well-block pressures in numerical reservoir simulation with non-square grid blocks and anisotropic permeability. Soc. Pet. Eng. J. 23(6), 531-543 (1983)
[41] Reynolds, AC; He, N.; Oliver, DS; Schatzinger, RA (ed.); Jordan, JF (ed.), Reducing uncertainty in geostatistical description with well testing pressure data, 149-162 (1999), Tulsa
[42] Reynolds, A.C., Zafari, M., Li, G.: Iterative forms of the ensemble Kalman filter. In: Proceedings of 10th European Conference on the Mathematics of Oil Recovery, Amsterdam, 4-7 September 2006
[43] Roberts, G.O., Gelman, A., Gilks, W.R.: Weak convergence and optimal scaling of random walk Metropolis algorithm. Ann. Appl. Probab. 7(1), 110-120 (1997) · Zbl 0876.60015
[44] Roberts, G.O., Rosenthal, J.S.: Optimal scaling of various Metropolis-Hastings algorithms. Stat. Sci. 16(4), 351-367 (2001) · Zbl 1127.65305
[45] Roggero, F., Hu, L.Y.: Gradual deformation of continuous geostatistical models for history matching. In: Proceedings of the SPE Annual Technical Conference and Exhibition, New Orleans, LA, SPE 49004, 27-30 September 1998
[46] Sakov, P., Bertino, L.: Relation between two common localisation methods for the EnKF. Comput. Geosci 15(2), 225-237(2011) · Zbl 1213.62150
[47] Sakov, P., Oke, P.R.: A deterministic formulation of the ensemble Kalman filter: an alternative to ensemble square root filters. Tellus, Ser. A Dyn. Meteorol. Oceanogr. 60(2), 361-371 (2008a)
[48] Sakov, P., Oke, P.R.: Implications of the form of the ensemble transformation in the ensemble square root filters. Mon. Weather Rev. 136, 1042-1053 (2008b)
[49] Skjervheim, J.-A., Evensen, G., Aanonsen, S.I., Ruud, B.O., Johansen, T.A.: Incorporating 4D seismic data in reservoir simulation models using ensemble Kalman filter. Soc. Pet. Eng. J. 12(3), 282-292 (2007)
[50] Sun, A.Y., Morris, A., Mohanty, S.: Comparison of deterministic ensemble Kalman filters for assimilating hydrogeological data. Adv. Water Resour. 32, 280-292 (2009)
[51] Tavakoli, R., Reynolds, A.A.: Monte Carlo simulation of permeability fields and reservoir performance predictions with SVD parameterization in RML compared with EnKF. Comput. Geosci. 15(1), 99-116 (2011) · Zbl 1209.86014
[52] Thulin, K., Li, G., Aanonsen, S.I., Reynolds, A.C.: Estimation of initial fluid contacts by assimilation of production data with EnKF. In: Proceedings of the SPE Annual Technical Conference and Exhibition, Anaheim, CA, SPE 109975, 11-14 November 2007
[53] Thulin, K., Nævdal, G., Skaug, H.J., Aanonsen, S.I.: Quantifying Monte Carlo uncertainty in the ensemble Kalman field. Soc. Pet. Eng. J. 16(1), 172-182 (2011)
[54] Tippett, M.K., Anderson, J.L., Bishop, C.H., Hamill, T.M., Whitaker, J.S.: Ensemble square-root filters. Mon. Weather Rev. 131, 1485-1490 (2003)
[55] van Leeuwen, P.J., Evensen, G.: Data assimilation and inverse methods in terms of a probabilistic formulation. Mon. Weather Rev. 124, 2898-2913 (1996)
[56] Wang, Y., Li, G., Reynolds, A.C.: Estimation of depths of fluid contacts by history matching using iterative ensemble-Kalman smoothers. Soc. Pet. Eng. J. 15(2) (2010)
[57] Wen, X.-H., Chen, W.H.: Real-time reservoir model updating using ensemble Kalman filter with confirming option. Soc. Pet. Eng. J. 11(4), 431-442 (2006)
[58] Wu, Z., Reynolds, A.C., Oliver, D.S.: Conditioning geostatistical models to two-phase production data. Soc. Pet. Eng. J. 3(2), 142-155 (1999)
[59] Zafari, M., Reynolds, A.C.: Assessing the uncertainty in reservoir description and performance predictions with the ensemble Kalman filter. Soc. Pet. Eng. J. 12(3), 382-391 (2007)
[60] Zhang, F., Reynolds, A.C.: Optimization algorithms for automatic history matching of production data. In: Proceedings of 8th European Conference on the Mathematics of Oil Recovery, Freiberg, Germany, 3-6 September 2002
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.