Combining sensitivities and prior information for covariance localization in the ensemble Kalman filter for petroleum reservoir applications

Sampling errors can severely degrade the reliability of estimates of conditional means and uncertainty quantification obtained by the application of the ensemble Kalman filter (EnKF) for data assimilation. A standard recommendation for reducing the spurious correlations and loss of variance due to sampling errors is to use covariance localization. In distance-based localization, the prior (forecast) covariance matrix at each data assimilation step is replaced with the Schur product of a correlation matrix with compact support and the forecast covariance matrix. The most important decision to be made in this localization procedure is the choice of the critical length(s) used to generate this correlation matrix. Here, we give a simple argument that the appropriate choice of critical length(s) should be based both on the underlying principal correlation length(s) of the geological model and the range of the sensitivity matrices. Based on this result, we implement a procedure for covariance localization and demonstrate with a set of distinctive reservoir history-matching examples that this procedure yields improved results over the standard EnKF implementation and over covariance localization with other choices of critical length.

[1]  J. Nocedal Updating Quasi-Newton Matrices With Limited Storage , 1980 .

[2]  A. Tarantola Inverse problem theory : methods for data fitting and model parameter estimation , 1987 .

[3]  Charles R. Johnson Matrix theory and applications , 1990 .

[4]  Charles R. Johnson,et al.  Topics in matrix analysis: The Hadamard product , 1991 .

[5]  P. Houtekamer,et al.  Data Assimilation Using an Ensemble Kalman Filter Technique , 1998 .

[6]  Dean S. Oliver,et al.  Conditioning Geostatistical Models to Two-Phase Production Data , 1999 .

[7]  S. Cohn,et al.  Ooce Note Series on Global Modeling and Data Assimilation Construction of Correlation Functions in Two and Three Dimensions and Convolution Covariance Functions , 2022 .

[8]  P. Houtekamer,et al.  A Sequential Ensemble Kalman Filter for Atmospheric Data Assimilation , 2001 .

[9]  J. Whitaker,et al.  Distance-dependent filtering of background error covariance estimates in an ensemble Kalman filter , 2001 .

[10]  A. Reynolds,et al.  Optimization Algorithms for Automatic History Matching of Production Data , 2002 .

[11]  Dean S. Oliver,et al.  History Matching of Three-Phase Flow Production Data , 2003 .

[12]  Andrew C. Lorenc,et al.  The potential of the ensemble Kalman filter for NWP—a comparison with 4D‐Var , 2003 .

[13]  Albert C. Reynolds,et al.  An Improved Implementation of the LBFGS Algorithm for Automatic History Matching , 2006 .

[14]  G. Evensen Data Assimilation: The Ensemble Kalman Filter , 2006 .

[15]  Albert C. Reynolds,et al.  Assessing the Uncertainty in Reservoir Description and Performance Predictions With the Ensemble Kalman Filter , 2007 .

[16]  Tor Arne Johansen,et al.  Incorporating 4D Seismic Data in Reservoir Simulation Models Using Ensemble Kalman Filter , 2007 .

[17]  Akhil Datta-Gupta,et al.  Efficient and Robust Reservoir Model Updating Using Ensemble Kalman Filter With Sensitivity-Based Covariance Localization , 2007 .

[18]  Jeffrey L. Anderson EXPLORING THE NEED FOR LOCALIZATION IN ENSEMBLE DATA ASSIMILATION USING A HIERARCHICAL ENSEMBLE FILTER , 2007 .

[19]  A. Reynolds,et al.  Estimation of Initial Fluid Contacts by Assimilation of Production Data With EnKF , 2007 .

[20]  T. Bengtsson,et al.  Estimation of high-dimensional prior and posterior covariance matrices in Kalman filter variants , 2007 .

[21]  Ning Liu,et al.  Inverse Theory for Petroleum Reservoir Characterization and History Matching , 2008 .

[22]  Deepak Devegowda,et al.  Streamline-Assisted Ensemble Kalman Filter for Rapid and Continuous Reservoir Model Updating , 2008 .

[23]  Dean S. Oliver,et al.  THE ENSEMBLE KALMAN FILTER IN RESERVOIR ENGINEERING-A REVIEW , 2009 .

[24]  Dean S. Oliver,et al.  Ensemble-Based Closed-Loop Optimization Applied to Brugge Field , 2010 .

[25]  A. Reynolds,et al.  Estimation of Depths of Fluid Contacts by History Matching Using Iterative Ensemble-Kalman Smoothers , 2010 .

[26]  D. Oliver,et al.  Cross-covariances and localization for EnKF in multiphase flow data assimilation , 2010 .