Ensemble size investigation in adaptive ES-MDA reservoir history matching

Author(s):  
Paulo Henrique Ranazzi ◽  
Marcio Augusto Sampaio
2019 ◽  
Vol 24 (1) ◽  
pp. 217-239
Author(s):  
Kristian Fossum ◽  
Trond Mannseth ◽  
Andreas S. Stordal

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.


Processes ◽  
2021 ◽  
Vol 9 (11) ◽  
pp. 1980
Author(s):  
Lihua Shen ◽  
Hui Liu ◽  
Zhangxin Chen

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.


SPE Journal ◽  
2019 ◽  
Vol 24 (04) ◽  
pp. 1490-1507 ◽  
Author(s):  
Sigurd Ivar Aanonsen ◽  
Svenn Tveit ◽  
Mathias Alerini

Summary This paper considers Bayesian methods to discriminate between models depending on posterior model probability. When applying ensemble-based methods for model updating or history matching, the uncertainties in the parameters are typically assumed to be univariate Gaussian random fields. In reality, however, there often might be several alternative scenarios that are possible a priori. We take that into account by applying the concepts of model likelihood and model probability and suggest a method that uses importance sampling to estimate these quantities from the prior and posterior ensembles. In particular, we focus on the problem of conditioning a dynamic reservoir-simulation model to frequent 4D-seismic data (e.g., permanent-reservoir-monitoring data) by tuning the top reservoir surface given several alternative prior interpretations with uncertainty. However, the methodology can easily be applied to similar problems, such as fault location and reservoir compartmentalization. Although the estimated posterior model probabilities will be uncertain, the ranking of models according to estimated probabilities appears to be quite robust.


2012 ◽  
Vol 17 (1) ◽  
pp. 83-97 ◽  
Author(s):  
Reza Tavakoli ◽  
Gergina Pencheva ◽  
Mary F. Wheeler ◽  
Benjamin Ganis

Sign in / Sign up

Export Citation Format

Share Document