Accelerating oil reservoir simulations using POD on the fly

Summary A reduced order model (ROM) is presented for the long-term calculation of sub-surface oil/water flows. As in several previous ROMs in the field, the Newton iterations in the full model (FM) equations, which are implicit in time, are projected onto a set of modes obtained by applying proper orthogonal decomposition (POD) to a set of snapshots computed by the FM itself. The novelty of the present ROM is that the POD modes are (i) first calculated from snapshots computed by the FM in a short initial stage, and then (ii) updated on the fly along the simulation itself, using new sets of snapshots computed by the FM in even shorter additional runs. Thus, the POD modes adapt themselves to the local dynamics along the simulation, instead of being completely calculated at the outset, which requires a computationally expensive preprocess. This strategy is robust and computationally efficient, which is tested in 10 and 30 year simulations for a realistic reservoir model taken from the SAIGUP project. This article is protected by copyright. All rights reserved.

[1]  Louis J. Durlofsky,et al.  Linearized reduced-order models for subsurface flow simulation , 2010, J. Comput. Phys..

[2]  David Ryckelynck,et al.  A robust adaptive model reduction method for damage simulations , 2011 .

[3]  José M. Vega,et al.  Adaptive POD‐based low‐dimensional modeling supported by residual estimates , 2015 .

[4]  José M. Vega,et al.  Reduced order models based on local POD plus Galerkin projection , 2010, J. Comput. Phys..

[5]  J. Jansen,et al.  Model-reduced gradient-based history matching , 2011 .

[6]  Jacob K. White,et al.  A trajectory piecewise-linear approach to model order reduction and fast simulation of nonlinear circuits and micromachined devices , 2001, IEEE/ACM International Conference on Computer Aided Design. ICCAD 2001. IEEE/ACM Digest of Technical Papers (Cat. No.01CH37281).

[7]  José M. Vega,et al.  Construction of Bifurcation Diagrams Using POD on the Fly , 2014, SIAM J. Appl. Dyn. Syst..

[8]  Margot Gerritsen,et al.  MODELING FLUID FLOW IN OIL RESERVOIRS , 2005 .

[9]  P. Beran,et al.  Reduced-order modeling: new approaches for computational physics , 2004 .

[10]  P. Sagaut,et al.  Towards an adaptive POD/SVD surrogate model for aeronautic design , 2011 .

[11]  Louis J. Durlofsky,et al.  Enhanced linearized reduced-order models for subsurface flow simulation , 2011, J. Comput. Phys..

[12]  E. Ansari.  Development of a surrogate simulator for two-phase subsurface flow simulation using trajectory piecewise linearization , 2014, Journal of Petroleum Exploration and Production Technology.

[13]  A. Velazquez,et al.  A method to generate computationally efficient reduced order models , 2009 .

[14]  Louis J. Durlofsky,et al.  Development and application of reduced‐order modeling procedures for subsurface flow simulation , 2009 .

[15]  Max D. Gunzburger,et al.  Centroidal Voronoi Tessellation-Based Reduced-Order Modeling of Complex Systems , 2006, SIAM J. Sci. Comput..

[16]  L. Sirovich.  Turbulence and the dynamics of coherent structures. I. Coherent structures , 1987 .

[17]  I. Kevrekidis,et al.  Equation-free/Galerkin-free POD-assisted computation of incompressible flows , 2005 .

[18]  A. Chatterjee.  An introduction to the proper orthogonal decomposition , 2000 .

[19]  José M. Vega,et al.  On the use of POD-based ROMs to analyze bifurcations in some dissipative systems , 2012 .

[20]  Charles-Henri Bruneau,et al.  Enablers for robust POD models , 2009, J. Comput. Phys..

[21]  V. Hurai,et al.  BASIC program for interpretation of microthermometric data from H 2 O and H 2 O-NaCl fluid inclusions , 1989 .

[22]  D. Ryckelynck,et al.  A priori hyperreduction method: an adaptive approach , 2005 .

[23]  Louis J. Durlofsky,et al.  Use of Reduced-Order Modeling Procedures for Production Optimization , 2010 .

[24]  A. Velazquez,et al.  Reduced-Order Model for Viscous Aerodynamic Flow Past an Airfoil , 2010 .

[25]  David Ryckelynck,et al.  A priori reduction method for solving the two-dimensional Burgers' equations , 2011, Appl. Math. Comput..

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

[27]  Eusebio Valero,et al.  Local POD Plus Galerkin Projection in the Unsteady Lid-Driven Cavity Problem , 2011, SIAM J. Sci. Comput..

[28]  D. Rempfer,et al.  On Low-Dimensional Galerkin Models for Fluid Flow , 2000 .

[29]  P. Holmes,et al.  The Proper Orthogonal Decomposition in the Analysis of Turbulent Flows , 1993 .

[30]  Stein Krogstad,et al.  Open-source MATLAB implementation of consistent discretisations on complex grids , 2012, Computational Geosciences.

[31]  Siep Weiland,et al.  Missing Point Estimation in Models Described by Proper Orthogonal Decomposition , 2004, IEEE Transactions on Automatic Control.

[32]  Accelerating CFD via local POD plus Galerkin projections , 2013 .

[33]  A. Velazquez,et al.  Generation of Aerodynamics Databases Using High-Order Singular Value Decomposition , 2008 .

[34]  Olav Møyner,et al.  MRST-AD - an Open-Source Framework for Rapid Prototyping and Evaluation of Reservoir Simulation Problems , 2015, ANSS 2015.