×

Investigation on principal component analysis parameterizations for history matching channelized facies models with ensemble-based data assimilation. (English) Zbl 1387.86044

Summary: Preserving plausible geological features when updating facies models is still one of the main challenges with ensemble-based history matching. This is particularly evident for fields with complex geological description (e.g., fluvial channels). There is an impressive amount of research published in the last few years about this subject. However, it appears that there is no definitive solution and both, academia and industry, are looking for practical and robust methods. Among the parameterizations traditionally investigated for history matching, the principal component analysis (PCA) of the prior covariance matrix is an efficient alternative to represent models described by two-point statistics. However, there are some recent developments extending PCA-based parameterizations for models described by multiple-point statistics. The first part of this paper presents an investigation on PCA-based schemes for parameterizing channelized facies models for history matching with ensemble-based methods. The following parameterizations are tested: standard PCA, two alternative implementations of kernel PCA and optimization-based PCA. In the second part of the paper, the optimization-based PCA is modified to allow the use of covariance localization and adapted for simultaneously adjusting the facies type and the permeability values within each facies when history matching production data with an ensemble-based method.

MSC:

86A32 Geostatistics

Software:

GSLIB
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Agbalaka CC, Oliver DS (2008) Application of the EnKF and localization to automatic history matching of facies distribution and production data. Math Geosci 40(4):353-374. doi:10.1007/s11004-008-9155-7 · Zbl 1142.86306 · doi:10.1007/s11004-008-9155-7
[2] Armstrong M, Galli A, Beucher H, Loc’h GL, Renard D, Doligez B, Eschard R, Geffroy F (2011) Plurigaussian simulations in geosciences, 2nd edn. Springer, Berlin, Heidelberg. doi:10.1007/978-3-642-19607-2 · Zbl 1393.86021
[3] Bianco A, Cominelli A, Dovera L, Nævdal G, Vallès B (2007) History matching and production forecast uncertainty by means of the ensemble Kalman filter: a real field application. In: Proceedings of the EAGE/EUROPEC Conference and Exhibition, SPE 107161. London, UK. doi:10.2118/107161-MS · Zbl 1425.86017
[4] Caers J, Zhang T (2004) Multiple-point geostatistics: a quantitative vehicle for integrating geologic analogs into multiple reservoir models. In: Integration of outcrop and modern analog data in reservoir models, vol 80. AAPG memoir, pp 383-394 · Zbl 1425.86017
[5] Chang H, Zhang D, Lu Z (2010) History matching of facies distributions with the EnKF and level set parameterization. J Comput Phys 229:8011-8030. doi:10.1016/j.jcp.2010.07.005 · Zbl 1425.86017 · doi:10.1016/j.jcp.2010.07.005
[6] Chang Y, Stordal AS, Valestran R (2015) Facies parameterization and estimation for complex reservoirs—the Brugge field. In: Proceedings of the SPE Bergen One Day Seminar, SPE 173872. Bergen, Norway. doi:10.2118/173872-MS
[7] Chen C, Gao G, Honorio J, Gelderblom P, Jimenez E, Jaakkola T (2014) Integration of principal-component-analysis and streamline information for the history matching of channelized reservoirs. In: Proceedings of the SPE Annual Technical Conference and Exhibition, SPE 170636. Amsterdam, The Netherlands. doi:10.2118/170636-MS
[8] Chen C, Gao G, Ramirez BA, Vink JC, Girardi AM (2015) Assisted history matching of channelized models using pluri-principal component analysis. In: Proceedings of the SPE Reservoir Simulation Symposium, Houston, SPE 173192. Texas, USA. doi:10.2118/173192-MS
[9] Chen Y, Oliver DS (2010) Ensemble-based closed-loop optimization applied to Brugge field. SPE Reserv Eval Eng 13(1):56-71. doi:10.2118/118926-PA · doi:10.2118/118926-PA
[10] Chen Y, Oliver DS (2014) History matching of the Norne full-field model with an iterative ensemble smoother. SPE Reservoir Eval Eng 17(2). doi:10.2118/164902-PA
[11] Cominelli A, Dovera L, Vimercati S, Nævdal G (2009) Benchmark study of ensemble Kalman filter methodology: History matching and uncertainty quantification for a deep-water oil reservoir. In: Proceedings of the International Petroleum Technology Conference, IPTC 13748. Doha, Qatar. doi:10.2523/13748-MS · Zbl 1159.86306
[12] Deutsch CV, Journel AG (1998) GSLIB: Geostatistical Software Library and User’s Guide, 2nd edn. Oxford University Press, New York
[13] Deutsch CV, Wang L (1996) Hierarchical object-based stochastic modeling of fluvial reservoirs. Math Geol 28(7):857-880. doi:10.1007/BF02066005 · doi:10.1007/BF02066005
[14] Emerick AA (2016) Analysis of the performance of ensemble-based assimilation of production and seismic data. J Pet Sci Eng 139:219-239. doi:10.1016/j.petrol.2016.01.029 · doi:10.1016/j.petrol.2016.01.029
[15] Emerick AA, Reynolds AC (2011a) Combining sensitivities and prior information for covariance localization in the ensemble Kalman filter for petroleum reservoir applications. Comput Geosci 15(2):251-269. doi:10.1007/s10596-010-9198-y · Zbl 1213.86009 · doi:10.1007/s10596-010-9198-y
[16] Emerick AA, Reynolds AC (2011b) History matching a field case using the ensemble Kalman filter with covariance localization. SPE Reserv Eval Eng 14(4):423-432. doi:10.2118/141216-PA · doi:10.2118/141216-PA
[17] Emerick AA, Reynolds AC (2013a) Ensemble smoother with multiple data assimilation. Comput Geosci 55:3-15. doi:10.1016/j.cageo.2012.03.011 · doi:10.1016/j.cageo.2012.03.011
[18] Emerick AA, Reynolds AC (2013b) History matching of production and seismic data for a real field case using the ensemble smoother with multiple data assimilation. In: Proceedings of the SPE Reservoir Simulation Symposium, SPE 163675. The Woodlands, Texas, USA. doi:10.2118/163675-MS
[19] Evensen G, Hove J, Meisingset HC, Reiso E, Seim KS, Espelid Ø (2007) Using the EnKF for assisted history matching of a North Sea reservoir model. In: Proceedings of the SPE Reservoir Simulation Symposium, SPE 106184. Houston, Texas. doi:10.2118/106184-MS
[20] Furrer R, Bengtsson T (2007) Estimation of high-dimensional prior and posterior covariance matrices in Kalman filter variants. J Multivar Anal 98(2):227-255. doi:10.1016/j.jmva.2006.08.003 · Zbl 1105.62091 · doi:10.1016/j.jmva.2006.08.003
[21] Gao G, Vink JC, Chen C, Alpak FO, Du K (2015) Enhanced reparameterization and data-integration algorithms for robust and efficient history matching of geologically complex reservoirs. In: SPE Annual Technical Conference and Exhibition, SPE 175039. Houston, Texas. doi:10.2118/175039-MS
[22] Gaspari G, Cohn SE (1999) Construction of correlation functions in two and three dimensions. Q J R Meteorolog Soc 125(554):723-757. doi:10.1002/qj.49712555417 · doi:10.1002/qj.49712555417
[23] Gavalas GR, Shah PC, Seinfeld JH (1976) Reservoir history matching by Bayesian estimation. SPE J 16(6):337-350. doi:10.2118/5740-PA · doi:10.2118/5740-PA
[24] Golub GH, van Loan CF (1989) Matrix Computations, 2nd edn. The Johns Hopkins University Press, Baltimore · Zbl 0733.65016
[25] Guardiano, FB; Srivastava, RM; Soares, A. (ed.), Multivariate geostatistics: Beyond bivariate moments, No. 5, 133-144 (1993), The Netherlands
[26] Haugen V, Nævdal G, Natvik LJ, Evensen G, Berg AM, Flornes KM (2008) History matching using the ensemble Kalman filter on a North Sea field case. SPE J 13(4):382-391. doi:10.2118/102430-PA · doi:10.2118/102430-PA
[27] Honorio J, Chen C, Gao G, Du K, Jaakkola T (2015) Integration of PCA with a novel machine learning method for reparameterization and assisted history matching geologically complex reservoirs. In: Proceedings of the SPE Annual Technical Conference and Exhibition, SPE 175038. Houston, Texas. doi:10.2118/175038-MS
[28] Houtekamer PL, Mitchell HL (1998) Data assimilation using an ensemble Kalman filter technique. Mon Weather Rev 126(3):796-811. doi:10.1175/1520-0493(1998)126<0796:DAUAEK>2.0.CO;2 · doi:10.1175/1520-0493(1998)126<0796:DAUAEK>2.0.CO;2
[29] Jafarpour B, Khodabakhshi M (2011) A probability conditioning method (PCM) for nonlinear flow data integration into multipoint statistical facies simulation. Math Geosci 43(2):133-164. doi:10.1007/s11004-011-9316-y · Zbl 1207.86010 · doi:10.1007/s11004-011-9316-y
[30] Jafarpour B, McLaughlin DB (2008) History matching with an ensemble Kalman filter and discrete cosine parameterization. Comput Geosci 12(2):227-244. doi:10.1007/s10596-008-9080-3 · Zbl 1159.86306 · doi:10.1007/s10596-008-9080-3
[31] Kwok JTY, Tsang IWH (2004) The pre-image problem in kernel methods. IEEE Trans Neural Networks 15(6). doi:10.1109/TNN.2004.837781
[32] LaVenue AM, Pickens JF (1992) Application of a coupled adjoint sensitivity and kriging approach to calibrate a groundwater flow model. Water Resour Res 28(6):1543-1569. doi:10.1029/92WR00208 · doi:10.1029/92WR00208
[33] Le DH, Younis R, Reynolds AC (2015) A history matching procedure for non-Gaussian facies based on ES-MDA. In: Proceedings of the SPE Reservoir Simulation Symposium, SPE 173233. Houston, Texas. doi:10.2118/173233-MS · Zbl 1392.86054
[34] Le DH, Emerick AA, Reynolds AC (2016) An adaptive ensemble smoother with multiple data assimilation for assisted history matching. SPE J. doi:10.2118/173214-PA (preprint)
[35] Li L, Srinivasan S, Zhou H, Gómez-Hernández JJ (2015) Two-point or multiple-point statistics? a comparison between the ensemble Kalman filtering and the ensemble pattern matching inverse methods. Adv Water Resour 86(Part B):297-310. doi:10.1016/j.advwatres.2015.05.014
[36] Liu N, Oliver DS (2005) Ensemble Kalman filter for automatic history matching of geologic facies. J Pet Sci Eng 47(3-4):147-161. doi:10.1016/j.petrol.2005.03.006 · doi:10.1016/j.petrol.2005.03.006
[37] Loc’h GL, Galli A (1997) Truncated plurigaussian method: Theoretical and practical points of view. In: Baafi EY, Schofield NA (eds) Geostatistics Wollongong ’96, vol 1. Kluwer Academic, USA, pp 211-222 · Zbl 1356.86020
[38] Lorentzen RJ, Flornes K, Nævdal G (2012) History channelized reservoirs using the ensemble Kalman filter. SPE J 17(1):137-151. doi:10.2118/143188-PA · doi:10.2118/143188-PA
[39] Mannseth T (2014) Relation between level set and truncated pluri-Gaussian methodologies for facies representation. Math Geosci 46(6):711-731. doi:10.1007/s11004-013-9507-9 · Zbl 1323.86045 · doi:10.1007/s11004-013-9507-9
[40] Mariethoz G, Caers J (2014) Multiple-point Geostatistics—Stochastic Modeling with Training Images. John Wiley & Sons, Ltd, New York
[41] 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, Holland · doi:10.1007/978-94-009-3701-7_16
[42] Moreno DL, Aanonsen SI (2011) Continuous facies updating using the ensemble Kalman filter and the level set method. Math Geosci 43(8):951-970. doi:10.1007/s11004-011-9347-4 · doi:10.1007/s11004-011-9347-4
[43] Moreno DL, Aanonsen SI, Evensen G, Skjervheim JA (2008) Channel facies estimation based on gaussian perturbations in the EnKF. In: Proceedings of the ECMOR XI—11th European Conference on the Mathematics of Oil Recovery. doi:10.3997/2214-4609.20146403
[44] Oliver DS (2014) Minimization for conditional simulation: Relationship to optimal transport. J Comput Phys 265:1-15. doi:10.1016/j.jcp.2014.01.048 · Zbl 1349.62444 · doi:10.1016/j.jcp.2014.01.048
[45] Peters L, Arts R, Brouwer G, Geel C, Cullick S, Lorentzen RJ, Chen Y, Dunlop N, Vossepoel FC, Xu R, Sarma P, Alhuthali AH, Reynolds A (2010) Results of the Brugge benchmark study for flooding optimisation and history matching. SPE Reserv Eval Eng 13(3):391-405. doi:10.2118/119094-PA · doi:10.2118/119094-PA
[46] Ping J, Zhang D (2014) History matching of channelized reservoirs with vector-based level-set parameterization. SPE J 19(3):514-529. doi:10.2118/169898-PA · doi:10.2118/169898-PA
[47] Reynolds AC, He N, Chu L, Oliver DS (1996) Reparameterization techniques for generating reservoir descriptions conditioned to variograms and well-test pressure data. SPE J 1(4):413-426. doi:10.2118/30588-PA · doi:10.2118/30588-PA
[48] Sarma P, Chen WH (2009) Generalization of the ensemble Kalman filter using kernel for non gaussian random fields. In: Proceedings of the SPE Reservoir Simulation Symposium, SPE 119177. The Woodlands, Texas. doi:10.2118/119177-MS · Zbl 1213.86009
[49] Sarma P, Durlofsky LJ, Aziz K (2008) Kernel principal component analysis for efficient differentiable parameterization of multipoint geostatistics. Math Geosci 40(1):3-32. doi:10.1007/s11004-007-9131-7 · Zbl 1144.86004 · doi:10.1007/s11004-007-9131-7
[50] Schölkopf B, Smola A, Müller KR (1998) Nonlinear component analysis as a kernel eigenvalue problem. Neural Comput 10(5):1299-1319. doi:10.1162/089976698300017467 · doi:10.1162/089976698300017467
[51] Sebacher BM, Hanea R, Heemink A (2013) A probabilistic parametrization for geological uncertainty estimation using the ensemble Kalman filter (EnKF). Comput Geosci 17(5):813-832. doi:10.1007/s10596-013-9357-z · Zbl 1393.86021 · doi:10.1007/s10596-013-9357-z
[52] Sebacher BM, Stordal AS, Hanea R (2015) Bridging multipoint statistics and truncated Gaussian fields for improved estimation of channelized reservoirs with ensemble methods. Comput Geosci 19(2):341-369. doi:10.1007/s10596-014-9466-3 · Zbl 1392.86054 · doi:10.1007/s10596-014-9466-3
[53] Skjervheim JA, Evensen G, Aanonsen SI, Ruud BO, Johansen TA (2007) Incorporating 4D seismic data in reservoir simulation models using ensemble Kalman filter. SPE J 12(3):282-292. doi:10.2118/95789-PA · doi:10.2118/95789-PA
[54] Stordal AS, Elsheikh AH (2015) Iterative ensemble smoothers in the annealed importance sampling framework. Adv Water Resour 86:231-239. doi:10.1016/j.advwatres.2015.09.030 · doi:10.1016/j.advwatres.2015.09.030
[55] Strebelle S (2002) Conditional simulation of complex geological structures using multiple-point statistics. Math Geol 34(1):1-21. doi:10.1023/A:1014009426274 · Zbl 1036.86013 · doi:10.1023/A:1014009426274
[56] Tavakoli R, Srinivasan S, Wheeler MF (2014) Rapid updating of stochastic models by use of an ensemble-filter approach. SPE J 19(3):500-513. doi:10.2118/163673-PA · doi:10.2118/163673-PA
[57] Vo HX, Durlofsky LJ (2014) A new differentiable parameterization based on principal component analysis for the low-dimensional representation of complex geological models. Math Geosci 46(7):775-813. doi:10.1007/s11004-014-9541-2 · Zbl 1323.86048 · doi:10.1007/s11004-014-9541-2
[58] Vo HX, Durlofsky LJ (2015) Data assimilation and uncertainty assessment for complex geological models using a new PCA-based parameterization. Comput Geosci 19(4):747-767. doi:10.1007/s10596-015-9483-x · Zbl 1392.86057 · doi:10.1007/s10596-015-9483-x
[59] Wang Y, Li G, Reynolds AC (2009) EnKF data assimilation with kernel PCA and discrete cosine transform parameterizations. TUPREP research report, The University of Tulsa
[60] Zhang Y, Oliver DS (2011) History matching using the ensemble Kalman filter with multiscale parameterization: A field case study. SPE J 16(2):307-317. doi:10.2118/118879-PA · doi:10.2118/118879-PA
[61] Zhang Y, Oliver DS, Chen Y, Skaug HJ (2015) Data assimilation by use of the iterative ensemble smoother for 2D facies models. SPE J 20(1):168-185. doi:10.2118/170248-PA · doi:10.2118/170248-PA
[62] Zhao H, Li G, Reynolds AC, Yao J (2013) Large-scale history matching with quadratic interpolation models. Comput Geosci 17(1):117-138. doi:10.1007/s10596-012-9320-4 · Zbl 1356.86020 · doi:10.1007/s10596-012-9320-4
[63] Zhao Y, Reynolds AC, Li G (2008) Generating facies maps by assimilating production data and seismic data with the ensemble Kalman filter. In: Proceedings of the SPE Improved Oil Recovery Symposium, SPE 113990. Tulsa, Oklahoma. doi:10.2118/113990-MS
[64] Zhou H, Gómez-Hernández JJ, Li L (2012) A pattern-search-based inverse method. Water Resour Res 48(3):W03,505. doi:10.1029/2011WR011195 · Zbl 1393.86021
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.