The Ensemble Kalman Filter for Continuous Updating of Reservoir Simulation Models

2005 ◽  
Vol 128 (1) ◽  
pp. 79-87 ◽  
Author(s):  
Yaqing Gu ◽  
Dean S. Oliver

This paper reports the use of ensemble Kalman filter (EnKF) for automatic history matching. EnKF is a Monte Carlo method, in which an ensemble of reservoir state variables are generated and kept up-to-date as data are assimilated sequentially. The uncertainty of reservoir state variables is estimated from the ensemble at any time step. Two synthetic problems are selected to investigate two primary concerns with the application of the EnKF. The first concern is whether it is possible to use a Kalman filter to make corrections to state variables in a problem for which the covariance matrix almost certainly provides a poor representation of the distribution of variables. It is tested with a one-dimensional, two-phase waterflood problem. The water saturation takes large values behind the flood front, and small values ahead of the front. The saturation distribution is bimodal and is not well modeled by the mean and variance. The second concern is the representation of the covariance via a relatively small ensemble of state vectors may be inadequate. It is tested by a two-dimensional, two-phase problem. The number of ensemble members is kept the same as for the one-dimensional problem. Hence the number of ensemble members used to create the covariance matrix is far less than the number of state variables. We conclude that EnKF can provide satisfactory history matching results while requiring less computation work than traditional history matching methods.

SPE Journal ◽  
2010 ◽  
Vol 16 (02) ◽  
pp. 331-342 ◽  
Author(s):  
Hemant A. Phale ◽  
Dean S. Oliver

Summary When the ensemble Kalman filter (EnKF) is used for history matching, the resulting updates to reservoir properties sometimes exceed physical bounds, especially when the problem is highly nonlinear. Problems of this type are often encountered during history matching compositional models using the EnKF. In this paper, we illustrate the problem using an example in which the updated molar density of CO2 in some regions is observed to take negative values while molar densities of the remaining components are increased. Standard truncation schemes avoid negative values of molar densities but do not address the problem of increased molar densities of other components. The results can include a spurious increase in reservoir pressure with a subsequent inability to maintain injection. In this paper, we present a method for constrained EnKF (CEnKF), which takes into account the physical constraints on the plausible values of state variables during data assimilation. In the proposed method, inequality constraints are converted to a small number of equality constraints, which are used as virtual observations for calibrating the model parameters within plausible ranges. The CEnKF method is tested on a 2D compositional model and on a highly heterogeneous three-phase-flow reservoir model. The effect of the constraints on mass conservation is illustrated using a 1D Buckley-Leverett flow example. Results show that the CEnKF technique is able to enforce the nonnegativity constraints on molar densities and the bound constraints on saturations (all phase saturations must be between 0 and 1) and achieve a better estimation of reservoir properties than is obtained using only truncation with the EnKF.


Water ◽  
2019 ◽  
Vol 11 (7) ◽  
pp. 1520
Author(s):  
Zheng Jiang ◽  
Quanzhong Huang ◽  
Gendong Li ◽  
Guangyong Li

The parameters of water movement and solute transport models are essential for the accurate simulation of soil moisture and salinity, particularly for layered soils in field conditions. Parameter estimation can be achieved using the inverse modeling method. However, this type of method cannot fully consider the uncertainties of measurements, boundary conditions, and parameters, resulting in inaccurate estimations of parameters and predictions of state variables. The ensemble Kalman filter (EnKF) is well-suited to data assimilation and parameter prediction in Situations with large numbers of variables and uncertainties. Thus, in this study, the EnKF was used to estimate the parameters of water movement and solute transport in layered, variably saturated soils. Our results indicate that when used in conjunction with the HYDRUS-1D software (University of California Riverside, California, CA, USA) the EnKF effectively estimates parameters and predicts state variables for layered, variably saturated soils. The assimilation of factors such as the initial perturbation and ensemble size significantly affected in the simulated results. A proposed ensemble size range of 50–100 was used when applying the EnKF to the highly nonlinear hydrological models of the present study. Although the simulation results for moisture did not exhibit substantial improvement with the assimilation, the simulation of the salinity was significantly improved through the assimilation of the salinity and relative solutetransport parameters. Reducing the uncertainties in measured data can improve the goodness-of-fit in the application of the EnKF method. Sparse field condition observation data also benefited from the accurate measurement of state variables in the case of EnKF assimilation. However, the application of the EnKF algorithm for layered, variably saturated soils with hydrological models requires further study, because it is a challenging and highly nonlinear problem.


2021 ◽  
Author(s):  
Marie Turčičová ◽  
Jan Mandel ◽  
Kryštof Eben

<p>A widely popular group of data assimilation methods in meteorological and geophysical sciences is formed by filters based on Monte-Carlo approximation of the traditional Kalman filter, e.g. <span>E</span><span>nsemble Kalman filter </span><span>(EnKF)</span><span>, </span><span>E</span><span>nsemble </span><span>s</span><span>quare-root filter and others. Due to the computational cost, ensemble </span><span>size </span><span>is </span><span>usually </span><span>small </span><span>compar</span><span>ed</span><span> to the dimension of the </span><span>s</span><span>tate </span><span>vector. </span><span>Traditional </span> <span>EnKF implicitly uses the sample covariance which is</span><span> a poor estimate of the </span><span>background covariance matrix - singular and </span><span>contaminated by </span><span>spurious correlations. </span></p><p><span>W</span><span>e focus on modelling the </span><span>background </span><span>covariance matrix by means of </span><span>a linear model for its inverse. This is </span><span>particularly </span><span>useful</span> <span>in</span><span> Gauss-Markov random fields (GMRF), </span><span>where</span> <span>the inverse covariance matrix has </span><span>a banded </span><span>structure</span><span>. </span><span>The parameters of the model are estimated by the</span><span> score matching </span><span>method which </span><span>provides</span><span> estimators in a closed form</span><span>, cheap to compute</span><span>. The resulting estimate</span><span> is a key component of the </span><span>proposed </span><span>ensemble filtering algorithms. </span><span>Under the assumption that the state vector is a GMRF in every time-step, t</span><span>he Score matching filter with Gaussian resamplin</span><span>g (SMF-GR) </span><span>gives</span><span> in every time-step a consistent (in the large ensemble limit) estimator of mean and covariance matrix </span><span>of the forecast and analysis distribution</span><span>. Further, we propose a filtering method called Score matching ensemble filter (SMEF), based on regularization of the EnK</span><span>F</span><span>. Th</span><span>is</span><span> filter performs well even for non-Gaussian systems with non-linear dynamic</span><span>s</span><span>. </span><span>The performance of both filters is illustrated on a simple linear convection model and Lorenz-96.</span></p>


2006 ◽  
Author(s):  
Vibeke Eilen Jensen Haugen ◽  
Lars-Jorgen Natvik ◽  
Geir Evensen ◽  
Aina Margrethe Berg ◽  
Kristin Margrethe Flornes ◽  
...  

2013 ◽  
Vol 30 (5) ◽  
pp. 1303-1312 ◽  
Author(s):  
Xiaogu Zheng ◽  
Guocan Wu ◽  
Shupeng Zhang ◽  
Xiao Liang ◽  
Yongjiu Dai ◽  
...  

Sign in / Sign up

Export Citation Format

Share Document