Closed-loop reservoir management on the Brugge test case

This paper proposes an augmented Lagrangian method for production optimization in which the cost function to be maximized is defined as an augmented Lagrangian function consisting of the net present value (NPV) and all the equality and inequality constraints except the bound constraints. The bound constraints are dealt with using a trust-region gradient projection method. The paper also presents a way to eliminate the need to convert the inequality constraints to equality constraints with slack variables in the augmented Lagrangian function, which greatly reduces the size of the optimization problem when the number of inequality constraints is large. The proposed method is tested in the context of closed-loop reservoir management benchmark problem based on the Brugge reservoir setup by TNO. In the test, we used the ensemble Kalman filter (EnKF) with covariance localization for data assimilation. Production optimization is done on the updated ensemble mean model from EnKF. The production optimization resulted in a substantial increase in the NPV for the expected reservoir life compared to the base case with reactive control.

[1]  C. W. Harper,et al.  A FORTRAN IV program for comparing ranking algorithms in quantitative biostratigraphy , 1984 .

[2]  Charles R. Johnson,et al.  Matrix analysis , 1985, Statistical Inference for Engineers and Data Scientists.

[3]  Harald Asheim,et al.  Maximization of Water Sweep Efficiency by Controlling Production and Injection Rates , 1988 .

[4]  P. Toint,et al.  Lancelot: A FORTRAN Package for Large-Scale Nonlinear Optimization (Release A) , 1992 .

[5]  G. Evensen.  Sequential data assimilation with a nonlinear quasi‐geostrophic model using Monte Carlo methods to forecast error statistics , 1994 .

[6]  Jorge Nocedal,et al.  A Limited Memory Algorithm for Bound Constrained Optimization , 1995, SIAM J. Sci. Comput..

[7]  I. Zakirov,et al.  Optimizing Reservoir Performance by Automatic Allocation of Well Rates , 1996 .

[8]  Jorge Nocedal,et al.  Algorithm 778: L-BFGS-B: Fortran subroutines for large-scale bound-constrained optimization , 1997, TOMS.

[9]  Stephen J. Wright,et al.  Numerical Optimization , 2018, Fundamental Statistical Inference.

[10]  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 .

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

[12]  Jan Dirk Jansen,et al.  Dynamic Optimization of Water Flooding with Smart Wells Using Optimal Control Theory , 2002 .

[13]  P. Houtekamer,et al.  Ensemble size, balance, and model-error representation in an ensemble Kalman filter , 2002 .

[14]  Jan Dirk Jansen,et al.  Dynamic Optimization of Waterflooding With Smart Wells Using Optimal Control Theory , 2004 .

[15]  Jan Dirk Jansen,et al.  Improved Reservoir Management Through Optimal Control and Continuous Model Updating , 2004 .

[16]  Louis J. Durlofsky,et al.  Implementation of Adjoint Solution for Optimal Control of Smart Wells , 2005 .

[17]  Geir Nævdal,et al.  Reservoir Monitoring and Continuous Model Updating Using Ensemble Kalman Filter , 2005 .

[18]  J. Jansen,et al.  Closed-loop reservoir management , 2005 .

[19]  L. Durlofsky,et al.  Efficient real-time reservoir management using adjoint-based optimal control and model updating , 2006 .

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

[21]  A. Cominelli,et al.  Production Optimization under Constraints Using Adjoint Gradients , 2006 .

[22]  Akhil Datta-Gupta,et al.  Optimal Waterflood Management Using Rate Control , 2006 .

[23]  Xuecheng Tai,et al.  Efficient Optimization of Production from Smart Wells Based on the Augmented Lagrangian Method , 2006 .

[24]  Rolf Johan Lorentzen,et al.  A New Approach For Dynamic Optimization Of Water Flooding Problems , 2006 .

[25]  P.J.P. Egberts,et al.  Optimal waterflood design using the adjoint method , 2007 .

[26]  Albert C. Reynolds,et al.  Characterization of the Measurement Error in Time-Lapse Seismic Data and Production Data with an EM Algorithm , 2007 .

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

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

[29]  Louis J. Durlofsky,et al.  Production Optimization With Adjoint Models Under Nonlinear Control-State Path Inequality Constraints , 2008 .

[30]  Gijs van Essen,et al.  Robust Waterflooding Optimization of Multiple Geological Scenarios , 2009 .

[31]  Roald Brouwer,et al.  Closed Loop Reservoir Management , 2009 .

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

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

[34]  Dongxiao Zhang,et al.  Efficient Ensemble-Based Closed-Loop Production Optimization , 2009 .

[35]  Geert Brouwer,et al.  Results of the Brugge Benchmark Study for Flooding Optimisation and History Matching , 2010 .

[36]  Chunhong Wang,et al.  Production Optimization in Closed-Loop Reservoir Management , 2009 .