scholarly journals Parallel Implementation of the Deterministic Ensemble Kalman Filter for Reservoir History Matching

Processes ◽  
2021 ◽  
Vol 9 (11) ◽  
pp. 1980
Author(s):  
Lihua Shen ◽  
Hui Liu ◽  
Zhangxin Chen
Keyword(s):  
Kalman Filter ◽  
Ensemble Kalman Filter ◽  
Message Passing ◽  
History Matching ◽  
Message Passing Interface ◽  
Parallel Implementation ◽  
Computational Time ◽  
Ensemble Size ◽  
Three Phase ◽  
Reservoir History Matching

In this paper, the deterministic ensemble Kalman filter is implemented with a parallel technique of the message passing interface based on our in-house black oil simulator. The implementation is separated into two cases: (1) the ensemble size is greater than the processor number and (2) the ensemble size is smaller than or equal to the processor number. Numerical experiments for estimations of three-phase relative permeabilities represented by power-law models with both known endpoints and unknown endpoints are presented. It is shown that with known endpoints, good estimations can be obtained. With unknown endpoints, good estimations can still be obtained using more observations and a larger ensemble size. Computational time is reported to show that the run time is greatly reduced with more CPU cores. The MPI speedup is over 70% for a small ensemble size and 77% for a large ensemble size with up to 640 CPU cores.

Use of Phase Streamlines for Covariance Localization in Ensemble Kalman Filter for Three-Phase History Matching

2012 ◽  
Vol 15 (03) ◽  
pp. 273-289 ◽  
Author(s):  
Shingo Watanabe ◽  
Akhil Datta-Gupta
Keyword(s):  
Kalman Filter ◽  
Gas Phase ◽  
Ensemble Kalman Filter ◽  
History Matching ◽  
Model Parameters ◽  
Ensemble Size ◽  
Three Phase ◽  
Cross Covariance ◽  
The Cross ◽  
Water Cut

Summary The ensemble Kalman filter (EnKF) has gained increased popularity for history matching and continuous reservoir-model updating. It is a sequential Monte Carlo approach that works with an ensemble of reservoir models. Specifically, the method uses cross covariance between measurements and model parameters estimated from the ensemble. For practical field applications, the ensemble size needs to be kept small for computational efficiency. However, this leads to poor approximations of the cross covariance and can cause loss of geologic realism from unrealistic model updates outside the region of the data influence and/or loss of variance leading to ensemble collapse. A common approach to remedy the situation is to limit the influence of the data through covariance localization. In this paper, we show that for three-phase-flow conditions, the region of covariance localization strongly depends on the underlying flow dynamics as well as on the particular data type that is being assimilated, for example, water cut or gas/oil ratio (GOR). This makes the traditional distance-based localizations suboptimal and, often, ineffective. Instead, we propose the use of water- and gas-phase streamlines as a means for covariance localization for water-cut- and GOR-data assimilation. The phase streamlines can be computed on the basis of individual-phase velocities which are readily available after flow simulation. Unlike the total velocity streamlines, phase streamlines can be discontinuous. We show that the discontinuities in water-phase and gas-phase streamlines naturally define the region of influence for water-cut and GOR data and provide a flow-relevant covariance localization during EnKF updating. We first demonstrate the validity of the proposed localization approach using a waterflood example in a quarter-five-spot pattern. Specifically, we compare the phase streamline trajectories with cross-covariance maps computed using an ensemble size of 2,000 for both water-cut and GOR data. The results show a close correspondence between the time evolution of phase streamlines and the cross-covariance maps of water-cut and GOR data. A benchmark uncertainty quantification (the PUNQ-S3) (Carter 2007) model application shows that our proposed localization outperforms the distance-based localization method. The updated models show improved forecasts while preserving geologic realism.

scholarly journals Assessment of multilevel ensemble-based data assimilation for reservoir history matching

2019 ◽  
Vol 24 (1) ◽  
pp. 217-239
Author(s):  
Kristian Fossum ◽  
Trond Mannseth ◽  
Andreas S. Stordal
Keyword(s):  
Kalman Filter ◽  
Data Assimilation ◽  
Ensemble Kalman Filter ◽  
History Matching ◽  
Computational Cost ◽  
Ensemble Size ◽  
Matching Problems ◽  
Reservoir Models ◽  
Reservoir History Matching ◽  
Computational Resources

AbstractMultilevel ensemble-based data assimilation (DA) as an alternative to standard (single-level) ensemble-based DA for reservoir history matching problems is considered. Restricted computational resources currently limit the ensemble size to about 100 for field-scale cases, resulting in large sampling errors if no measures are taken to prevent it. With multilevel methods, the computational resources are spread over models with different accuracy and computational cost, enabling a substantially increased total ensemble size. Hence, reduced numerical accuracy is partially traded for increased statistical accuracy. A novel multilevel DA method, the multilevel hybrid ensemble Kalman filter (MLHEnKF) is proposed. Both the expected and the true efficiency of a previously published multilevel method, the multilevel ensemble Kalman filter (MLEnKF), and the MLHEnKF are assessed for a toy model and two reservoir models. A multilevel sequence of approximations is introduced for all models. This is achieved via spatial grid coarsening and simple upscaling for the reservoir models, and via a designed synthetic sequence for the toy model. For all models, the finest discretization level is assumed to correspond to the exact model. The results obtained show that, despite its good theoretical properties, MLEnKF does not perform well for the reservoir history matching problems considered. We also show that this is probably caused by the assumptions underlying its theoretical properties not being fulfilled for the multilevel reservoir models considered. The performance of MLHEnKF, which is designed to handle restricted computational resources well, is quite good. Furthermore, the toy model is utilized to set up a case where the assumptions underlying the theoretical properties of MLEnKF are fulfilled. On that case, MLEnKF performs very well and clearly better than MLHEnKF.

scholarly journals Parameters Estimation and Prediction of Water Movement and Solute Transport in Layered, Variably Saturated Soils Using the Ensemble Kalman Filter

Water ◽  
2019 ◽  
Vol 11 (7) ◽  
pp. 1520
Author(s):  
Zheng Jiang ◽  
Quanzhong Huang ◽  
Gendong Li ◽  
Guangyong Li
Keyword(s):  
Kalman Filter ◽  
Solute Transport ◽  
Ensemble Kalman Filter ◽  
Water Movement ◽  
Parameters Estimation ◽  
Saturated Soils ◽  
Hydrological Models ◽  
State Variables ◽  
Ensemble Size ◽  
Highly Nonlinear

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.

Sign in / Sign up

Export Citation Format

Share Document