Data Assimilation for Strongly Nonlinear Problems by Transformed Ensemble Kalman Filter

SPE Journal ◽  
2014 ◽  
Vol 20 (01) ◽  
pp. 202-221 ◽  
Author(s):  
Qinzhuo Liao ◽  
Dongxiao Zhang

Summary The ensemble Kalman filter (EnKF) has been widely used for data assimilation. It is challenging, however, when the relation of state and observation is strongly nonlinear. For example, near the flooding front in an immiscible flow, directly updating the saturation by use of the EnKF may lead to nonphysical results. One possible solution, which may be referred to as the restarted EnKF (REnKF), is to update the static state (e.g., permeability and porosity) and rerun the forward model from the initial time to obtain the updated dynamic state (e.g., pressure and saturation). However, it may become time-consuming, especially when the number of assimilation steps is large. In this study, we develop a transformed EnKF (TEnKF), in which the state is represented by displacement as an alternative variable. The displacement is first transformed from the forecasted state, then updated, and finally transformed back to obtain the updated state. Because the relation between displacement and observation is relatively linear, this new method provides a physically meaningful updated state without resolving the forward model. The TEnKF is tested in the history matching of multiphase flow in a 1D homogeneous medium, a 2D heterogeneous reservoir, and a 3D PUNQ-S3 model. The case studies show that the TEnKF produces physical results without the oscillation problem that occurs in the traditional EnKF, whereas the computational effort is reduced compared with the REnKF.

2008 ◽  
Vol 8 (11) ◽  
pp. 2975-2983 ◽  
Author(s):  
C. Lin ◽  
Z. Wang ◽  
J. Zhu

Abstract. An Ensemble Kalman Filter (EnKF) data assimilation system was developed for a regional dust transport model. This paper applied the EnKF method to investigate modeling of severe dust storm episodes occurring in March 2002 over China based on surface observations of dust concentrations to explore the impact of the EnKF data assimilation systems on forecast improvement. A series of sensitivity experiments using our system demonstrates the ability of the advanced EnKF assimilation method using surface observed PM10 in North China to correct initial conditions, which leads to improved forecasts of dust storms. However, large errors in the forecast may arise from model errors (uncertainties in meteorological fields, dust emissions, dry deposition velocity, etc.). This result illustrates that the EnKF requires identification and correction model errors during the assimilation procedure in order to significantly improve forecasts. Results also show that the EnKF should use a large inflation parameter to obtain better model performance and forecast potential. Furthermore, the ensemble perturbations generated at the initial time should include enough ensemble spreads to represent the background error after several assimilation cycles.


SPE Journal ◽  
2011 ◽  
Vol 16 (02) ◽  
pp. 294-306 ◽  
Author(s):  
Lingzao Zeng ◽  
Haibin Chang ◽  
Dongxiao Zhang

Summary The ensemble Kalman filter (EnKF) has been used widely for data assimilation. Because the EnKF is a Monte Carlo-based method, a large ensemble size is required to reduce the sampling errors. In this study, a probabilistic collocation-based Kalman filter (PCKF) is developed to adjust the reservoir parameters to honor the production data. It combines the advantages of the EnKF for dynamic data assimilation and the polynomial chaos expansion (PCE) for efficient uncertainty quantification. In this approach, all the system parameters and states and the production data are approximated by the PCE. The PCE coefficients are solved with the probabilistic collocation method (PCM). Collocation realizations are constructed by choosing collocation point sets in the random space. The simulation for each collocation realization is solved forward in time independently by means of an existing deterministic solver, as in the EnKF method. In the analysis step, the needed covariance is approximated by the PCE coefficients. In this study, a square-root filter is employed to update the PCE coefficients. After the analysis, new collocation realizations are constructed. With the parameter collocation realizations as the inputs and the state collocation realizations as initial conditions, respectively, the simulations are forwarded to the next analysis step. Synthetic 2D water/oil examples are used to demonstrate the applicability of the PCKF in history matching. The results are compared with those from the EnKF on the basis of the same analysis. It is shown that the estimations provided by the PCKF are comparable to those obtained from the EnKF. The biggest improvement of the PCKF comes from the leading PCE approximation, with which the computational burden of the PCKF can be greatly reduced by means of a smaller number of simulation runs, and the PCKF outperforms the EnKF for a similar computational effort. When the correlation ratio is much smaller, the PCKF still provides estimations with a better accuracy for a small computational effort.


2014 ◽  
Vol 142 (8) ◽  
pp. 2915-2934 ◽  
Author(s):  
Hailing Zhang ◽  
Zhaoxia Pu

Abstract A series of numerical experiments are conducted to examine the impact of surface observations on the prediction of landfalls of Hurricane Katrina (2005), one of the deadliest disasters in U.S. history. A specific initial time (0000 UTC 25 August 2005), which led to poor prediction of Hurricane Katrina in several previous studies, is selected to begin data assimilation experiments. Quick Scatterometer (QuikSCAT) ocean surface wind vectors and surface mesonet observations are assimilated with the minimum central sea level pressure and conventional observations from NCEP into an Advanced Research version of the Weather Research and Forecasting Model (WRF) using an ensemble Kalman filter method. Impacts of data assimilation on the analyses and forecasts of Katrina’s track, landfalling time and location, intensity, structure, and rainfall are evaluated. It is found that the assimilation of QuikSCAT and mesonet surface observations can improve prediction of the hurricane track and structure through modifying low-level thermal and dynamical fields such as wind, humidity, and temperature and enhancing low-level convergence and vorticity. However, assimilation of single-level surface observations alone does not ensure reasonable intensity forecasts because of the lack of constraint on the mid- to upper troposphere. When surface observations are assimilated with other conventional data, obvious enhancements are found in the forecasts of track and intensity, realistic convection, and surface wind structures. More importantly, surface data assimilation results in significant improvements in quantitative precipitation forecasts (QPFs) during landfalls.


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.


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.


2007 ◽  
Vol 7 (6) ◽  
pp. 17511-17536
Author(s):  
C. Lin ◽  
Z. Wang ◽  
J. Zhu

Abstract. An Ensemble Kalman Filter (EnKF) data assimilation system was developed for a regional dust transport model. This paper applied the EnKF method to investigate modeling severe dust storm episodes occurred in March 2002 over China based on surface observations of dust concentrations to explore its impacts on forecast improvement. A series of sensitivity experiments using our system reveals that the EnKF is an advanced assimilation method to afford better initial conditions with surface observed PM10 in North China and lead to improved forecasts of dust storms, but forecast with large errors can be made by model errors. This result illustrates that it requires identifying and correcting model errors during the assimilation procedure in order to significantly improve forecasts. Results also show that the EnKF should use a large inflation parameter to obtain better model performance and forecast potential. Furthermore, the ensemble perturbations generated at the initial time should include enough ensemble spreads to represent the background error after several assimilation cycles.


Author(s):  
Nicolas Papadakis ◽  
Etienne Mémin ◽  
Anne Cuzol ◽  
Nicolas Gengembre

2016 ◽  
Vol 66 (8) ◽  
pp. 955-971 ◽  
Author(s):  
Stéphanie Ponsar ◽  
Patrick Luyten ◽  
Valérie Dulière

Icarus ◽  
2010 ◽  
Vol 209 (2) ◽  
pp. 470-481 ◽  
Author(s):  
Matthew J. Hoffman ◽  
Steven J. Greybush ◽  
R. John Wilson ◽  
Gyorgyi Gyarmati ◽  
Ross N. Hoffman ◽  
...  

Sign in / Sign up

Export Citation Format

Share Document