×

Theoretical and efficient practical procedures for the generation of inflation factors for ES-MDA. (English) Zbl 1379.65033

Summary: The ensemble smoother with multiple data assimilation (ES-MDA) has proved to be a powerful assisted history-matching method. The main drawback of ES-MDA is that the inflation factors for damping the changes in model parameters have to be determined before starting the history-match. Although various authors have provided suggestions for determining the inflation factors adaptively as the history-match proceeds, these methods often result in a large number of data assimilation steps which can make ES-MDA too computationally inefficient for practical application to large-scale field problems. Here, we provide a theoretical procedure to determine exactly the minimum inflation factor at each data assimilation step that ensures the discrepancy principle is satisfied. Like previous adaptive ES-MDA methods, this method does not allow one to specify a priori the number of data assimilation steps to be done. Thus, using the exact theoretical procedure as a guide, we provide a practical efficient method for determining the inflation factors which allows one to specify a priori the number of data assimilation steps to be done with ES-MDA which still ensures that the initial inflation factor is chosen so that the discrepancy principle is approximately satisfied.

MSC:

65J22 Numerical solution to inverse problems in abstract spaces
65J20 Numerical solutions of ill-posed problems in abstract spaces; regularization
47J06 Nonlinear ill-posed problems
80A23 Inverse problems in thermodynamics and heat transfer

Software:

EnKF
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Aanonsen S I, Næ vdal G, Oliver D S, Reynolds A C and Vallés B 2009 The ensemble Kalman filter in reservoir engineering-a review SPE J.14 393-412
[2] Anderson J L 2007 Exploring the need for localization in ensemble data assimilation using hierarchical ensemble filter Physica D 230 99-111 · Zbl 1387.86039
[3] Barker J W, Cuypers M and Holden L 2001 Quantifying uncertainty in production forecasts: another look at the PUNQ-S3 problem SPE J.6 433-41
[4] Barlow R 1989 Statistics, Guide to the Use of Statistical Methods in the Physical Sciences (New York: Wiley)
[5] Bianco A, Cominelli A, Dovera L, Næ vdal G and Vallès B 2007 History matching and production forecast uncertainty by means of the ensemble Kalman filter: real field application Proc. of the EAGE/EUROPEC Conf. and Exhibition(London, UK, 11-4 June, 2007) p 107161
[6] Burgers G, van Leeuwen P and Evensen G 1998 Analysis scheme in the ensemble Kalman filter Mon. Weather Rev.126 1719-24
[7] Chen Y and Oliver D S 2009 Cross-covariance and localization for EnKF in multiphase flow data assimilation Comput. Geosci.14 579-601 · Zbl 1381.86026
[8] Chen Y and Oliver D S 2010 Ensemble-based closed-loop optimization applied to Brugge field SPE Reservoir Eval. Eng.13 56-71
[9] Chen Y and Oliver D S 2013 Levenberg-Marquardt forms of the iterative ensemble smoother for efficient history matching and uncertainty quantification Comput. Geosci.17 689-703 · Zbl 1382.65031
[10] Chen Y and Oliver D S 2014 History matching of the Norne full-field model with an iterative ensemble smoother SPE Reservoir Eval. Eng.17 244-56
[11] Devegowda D, Arroyo-Negrete E, Datta-Gupta A and Douma S G 2007 Efficient and robust reservoir model updating using ensemble Kalman filter with sensitivity-based covariance localization Proc. of the SPE Reservoir Simulation Symp.(Houston, Texas, 26-8 February, 2007) p 06144
[12] Emerick A A 2016 Analysis of the performance of ensemble-based assimilation of production and seismic data J. Pet. Sci. Eng.139 219-39
[13] Emerick A A and Reynolds A C 2011 History matching field case using the ensemble Kalman filter with covariance localization Proc. of the SPE Reservoir Simulation Symp.(The Woodlands, Texas, USA, 21-3 February) p 141216
[14] Emerick A A and Reynolds A C 2012a Combining the ensemble Kalman filter with Markov chain Monte Carlo for improved history matching and uncertainty characterization SPE J.17 639-59
[15] Emerick A A and Reynolds A C 2012b History matching time-lapse seismic data using the ensemble Kalman filter with multiple data assimilations Comput. Geosci.16 639-59
[16] Emerick A A and Reynolds A C 2013a Ensemble smoother with multiple data assimilations Comput. Geosci.55 3-15
[17] Emerick A A and Reynolds A C 2013b History-matching production and seismic data in real field case using the ensemble smoother with multiple data assimilation Proc. of the SPE Reservoir Simulation Symp.(The Woodlands, Texas, USA, 18-20 February, 2013) p 163645
[18] Emerick A A and Reynolds A C 2013c Investigation of the sampling performance of ensemble-based methods with simple reservoir model Comput. Geosci.17 325-50 · Zbl 1382.86018
[19] Evensen G 1994 Sequential data assimilation with nonlinear quasi-geostrophic model using Monte Carlo methods to forecast error statistics J. Geophys. Res.99 10143-62
[20] Evensen G 2007 Data Assimilation: the Ensemble Kalman Filter (Berlin: Springer) · Zbl 1157.86001
[21] Floris F J T, Bush M D, Cuypers M, Roggero F and Syversveen A R 2001 Methods for quantifying the uncertainty of production forecasts: comparative study Pet. Geosci.7 87-96
[22] Gao G and Reynolds A C 2006 An improved implementation of the LBFGS algorithm for automatic history matching SPE J.11 5-17
[23] Groetsch C W 1984 The Theory of Tikhonov Regularization for Fredholm Equations of the First Kind (Boston, MA: Pitman Publication) · Zbl 0545.65034
[24] Hanke M 1997 A regularizing Levenberg-Marquardt scheme, with applications to inverse groundwater filtration problems Inverse Probl.13 79-95 · Zbl 0873.65057
[25] Hanke M 2010 The regularizing Levenberg-Marquardt scheme is of optimal order J. Integral Equ. Appl.22 259-83 · Zbl 1204.65059
[26] Haugen V, Næ vdal G, Natvik L J, Evensen G, Berg A M and Flornes K M 2008 History matching using the ensemble Kalman filter on North Sea field case SPE J.13 382-91
[27] Iglesias M A 2015 Iterative regularization for ensemble data assimilation in reservoir models Comput. Geosci.19 177-212 · Zbl 1330.86015
[28] Iglesias M A and Dawson C 2013 The regularizing Levenberg-Marquardt scheme for history matching of petroleum reservoirs Comput. Geosci.17 1033-53 · Zbl 1393.86016
[29] Le D H, Emerick A A and Reynolds A C 2015 An adaptive ensemble smoother with multiple data assimilation for assisted history matching SPE J.21 2195-207
[30] Li G and Reynolds A C 2009 Iterative ensemble Kalman filters for data assimilation SPE J.14 496-505
[31] Liu N and Oliver D S 2003 Evaluation of Monte Carlo methods for assessing uncertainty SPE J.8 188-95
[32] Luo X, Stordal A S, Lorentzen R J and Næ vdal G 2015 Iterative ensemble smoother as an approximate solution to regularized minimum-average-cost problem: theory and applications SPE J.20 962-82
[33] Oliver D S, Reynolds A C and Liu N 2008 Inverse Theory for Petroleum Reservoir Characterization and History Matching (Cambridge, UK: Cambridge University Press)
[34] Peters E, Chen Y, Leeuwenburgh O and Oliver D 2013 Extended Brugge benchmark case for history matching and water flooding optimization Comput. Geosci.50 16-24
[35] Peters L et al 2010 Results of the Brugge benchmark study for flooding optimisation and history matching SPE Reservoir Eval. Eng.13 391-405
[36] Reynolds A C, Zafari M and Li G 2006 Iterative forms of the ensemble Kalman filter Proc. of 10th European Conf. on the Mathematics of Oil Recovery(Amsterdam, September 2006,)
[37] Rommelse J 2009 Data assimilation in reservoir management PhD Thesis Technical University of Delft, Delft, The Netherlands
[38] Seiler A, Aanonsen S I, Evensen G and Reivenæ s J C 2010 Structural surface uncertainty modeling and updating using the ensemble Kalman filter SPE J.15 1062-76
[39] Shirangi M G 2011 History matching production data with truncated SVD parameterization, Masters Thesis The University of Tulsa
[40] Shirangi M G 2014 History matching production data and uncertainty assessment with an efficient TSVD parameterization algorithm J. Pet. Sci. Eng.113 54-71
[41] Shirangi M G and Emerick A A 2016 An improved TSVD-based Levenberg-Marquardt algorithm for history matching and comparison with Gauss-Newton J. Pet. Sci. Eng.143 258-71
[42] Skjervheim J A, Evensen G, Aanonsen S I, Ruud B O and Johansen T A 2007 Incorporating 4D seismic data in reservoir simulation models using ensemble Kalman filter SPE J.12 282-92
[43] Tarantola A 2005 Inverse Problem Theory and Methods for Model Parameter Estimation (Philadelphia: SIAM)
[44] Tavakoli R and Reynolds A C 2010 History matching with parameterization based on the SVD of dimensionless sensitivity matrix SPE J.15 495-508
[45] Tavakoli R and Reynolds A C 2011 Monte Carlo simulation of permeability fields and reservoir performance predictions with SVD parameterization in RML compared with EnKF Comput. Geosci.15 99-116 · Zbl 1209.86014
[46] Thulin K, Li G, Aanonsen S I and Reynolds A C 2007 Estimation of initial fluid contacts by assimilation of production data with EnKF Proc. of the SPE Annual Technical Conf. and Exhibition(Anaheim, California, 11-4 November, 2007) p 109975
[47] van Leeuwen P J 1999 Comment on data assimilation using an ensemble Kalman filter technique Mon. Weather Rev.127 1374-7
[48] van Leeuwen P J and Evensen G 1996 Data assimilation and inverse methods in terms of probabilistic formulation Mon. Weather Rev.124 2898-13
[49] Vogel C R 2002 Computational Methods for Inverse Problems (Philadelphia: SIAM)
[50] Wang Y, Li G and Reynolds A C 2010 Estimation of depths of fluid contacts by history matching using iterative ensemble-Kalman smoothers SPE J.15 509-25
[51] Zafari M and Reynolds A C 2007 Assessing the uncertainty in reservoir description and performance predictions with the ensemble Kalman filter SPE J.12 382-91
[52] Zhang F, Reynolds A C and Oliver D S 2002 Evaluation of the reduction in uncertainty obtained by conditioning 3D stochastic channel to multiwell pressure data Math. Geol.34 713-40
[53] Zhang Y and Oliver D S 2011 Evaluation and error analysis: Kalman gain regularization versus covariance regularization Comput. Geosci.15 489-508 · Zbl 1254.86030
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.