Use of POD in control of flow through porous media

During the design of development concepts for the exploitation of oil and gas reservoirs, frequent use is made of numerical simulation of the flow of multi-phase fluids through porous rock. Recently, increased attention has been paid to systematic optimization of well positions and operating parameters (rates, pressures) over the life of the reservoir. Here we consider optimization of the displacement of oil towards production wells through the injection of water in other wells. Model-based optimal control of this “water flooding” process generally involves multiple simulations, which makes it into a time-consuming process. A potential way to address this issue is through the use of proper orthogonal decomposition (POD), We addressed the scope to speed up optimization of water-flooding a heterogeneous reservoir with multiple injectors and producers. We used an adjoint-based optimal control methodology that requires multiple passes of forward simulation of the reservoir model and backward simulation of an adjoint system of equations. We developed a nested approach in which POD was first used to reduce the state space dimensions of both the forward model and the adjoint system. After obtaining an optimized injection and production strategy using the reduced-order system, we verified the results using the original, high-order model. If necessary, we repeated the optimization cycle using new reduced-order systems based on snapshots from the verification run We tested the methodology on a reservoir model with 882 states (441 pressures, 441 saturations) and an adjoint model of 882 states (Lagrange multipliers). We obtained reduced-order models with 35-43 states only. The reduction in computing time was 52%.

[1]  Gene H. Golub,et al.  Matrix computations , 1983 .

[2]  K. Kunisch,et al.  Control of the Burgers Equation by a Reduced-Order Approach Using Proper Orthogonal Decomposition , 1999 .

[3]  Y. Yortsos,et al.  Optimization of fluid front dynamics in porous media using rate control. I. Equal mobility fluids , 2000 .

[4]  S. S. Ravindran,et al.  Reduced-Order Adaptive Controllers for Fluid Flows Using POD , 2000, J. Sci. Comput..

[5]  H. Tran,et al.  Modeling and control of physical processes using proper orthogonal decomposition , 2001 .

[6]  Jan Dirk Jansen,et al.  Generation of Low-Order Reservoir Models Using System-Theoretical Concepts , 2004 .

[7]  R. D. Prabhu,et al.  The influence of control on proper orthogonal decomposition of wall-bounded turbulent flows , 2001 .

[8]  K. Aziz,et al.  Petroleum Reservoir Simulation , 1979 .

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

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

[11]  P. Holmes,et al.  Turbulence, Coherent Structures, Dynamical Systems and Symmetry , 1996 .

[12]  Jan Dirk Jansen,et al.  Accelerating iterative solution methods using reduced‐order models as solution predictors , 2006 .

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

[14]  Arnold Heemink,et al.  Reduced models for linear groundwater flow models using empirical orthogonal functions , 2004 .

[15]  J. Jansen,et al.  Reduced-order optimal control of water flooding using proper orthogonal decomposition , 2006 .

[16]  D. W. Peaceman Fundamentals of numerical reservoir simulation , 1977 .