Influence of the polymer properties and numerical schemes on tertiary oil recovery processes

Abstract Chemical Enhanced Oil Recovery (EOR) processes comprise a number of flooding techniques aimed at increasing the operational life of mature oilfields. Among these, polymer flooding is one of the most developed; its functionality is to increment the aqueous viscosity, avoiding the formation of viscous fingering. Reservoir simulators consider this influence as well as other physical properties (e.g., adsorption, permeability reduction). However, the polymer degradation is usually not considered even though it plays a critical role in the viscosity. In this paper this mechanism is analyzed and coupled with the previously mentioned physical phenomena in order to present a complete study of their influence in the EOR process. Moreover, since a fully second-order accuracy scheme is used along with a Total Variation Diminishing (TVD) flux-limiting function, the influence of the latter on the recovery factor is also discussed. Results showed that the negative effect of the polymer adsorption was the most relevant physical phenomenon in terms of the oil recovery. Furthermore, the analysis of the discretization of the differential equations showed that traditional, linear first-order schemes created numerical diffusion affecting negatively the macroscopic sweeping efficiency, which disappeared when TVD techniques were used. Reservoir simulators allow determining the desired designing properties for future polymers in relationship with the characteristics of the oilfield to be exploited.

[1]  Mustapha El Ossmani,et al.  A sequential approach for numerical simulation of two-phase multicomponent flow with reactive transport in porous media , 2017, Math. Comput. Simul..

[2]  L. Dake.  Fundamentals of Reservoir Engineering , 1983 .

[3]  Hassan Smaoui,et al.  Flux-limiting techniques for simulation of pollutant transport in porous media: Application to groundwater management , 2008, Math. Comput. Model..

[4]  Weidong Zhang,et al.  Study on Mechanism of Viscoelastic Polymer Transient Flow in Porous Media , 2017 .

[5]  Gary A. Pope,et al.  Mechanistic Interpretation and Utilization of Viscoelastic Behavior of Polymer Solutions for Improved Polymer-Flood Efficiency , 2008 .

[6]  Kamy Sepehrnoori,et al.  Comparison of an Improved Compositional Micellar/Polymer Simulator With Laboratory Corefloods , 1987 .

[7]  Kamy Sepehrnoori,et al.  Description of an Improved Compositional Micellar/Polymer Simulator , 1987 .

[8]  A. Kamyabi,et al.  Simulation of two generalised Newtonian fluids flow in micropore with dead end , 2011 .

[9]  L. Lake,et al.  Enhanced Oil Recovery , 2017 .

[10]  Yuanle Ma,et al.  Computational methods for multiphase flows in porous media , 2007, Math. Comput..

[11]  Taha Sochi,et al.  Flow of non‐newtonian fluids in porous media , 2009, 0907.2399.

[12]  Huiqing Liu,et al.  Mechanistic Simulation Studies on Viscous-Elastic Polymer Flooding in Petroleum Reservoirs , 2013 .

[13]  J. Bear.  Dynamics of Fluids in Porous Media , 1975 .

[14]  Z. Zhang,et al.  Microscopic Roles of “Viscoelasticity” in HPMA polymer flooding for EOR , 2011 .

[15]  M. Kamyabi,et al.  Effects of viscoelastic polymer solutions on the oil extraction from dead ends , 2013 .

[16]  Pietro Tesi,et al.  Numerical modeling of a compositional flow for chemical EOR and its stability analysis , 2017 .

[17]  Jacob Bear,et al.  Macroscopic modelling of transport phenomena in porous media. 2: Applications to mass, momentum and energy transport , 1986 .

[18]  A. Skauge,et al.  Mechanisms of Non-Newtonian Polymer Flow Through Porous Media Using Navier–Stokes Approach , 2015 .

[19]  Numerical Aspects of the Convection–Dispersion Equation , 2014 .

[20]  Q. Tran,et al.  An analysis of physical models and numerical schemes for polymer flooding simulations , 2016, Computational Geosciences.

[21]  T. S. Mykkeltvedt,et al.  Fully implicit higher-order schemes applied to polymer flooding , 2016, Computational Geosciences.

[22]  Kamy Sepehrnoori,et al.  Application of higher-order flux-limited methods in compositional simulation , 1994 .

[23]  G. B. Savioli,et al.  On The Numerical Modeling For Surfactant Flooding Of Oil Reservoirs. , 2002 .

[24]  Pablo Druetta,et al.  Numerical Modeling and Validation of a Novel 2D Compositional Flooding Simulator Using a Second-Order TVD Scheme , 2018, Energies.

[25]  Taha Sochi,et al.  Non-Newtonian flow in porous media , 2010 .

[26]  J. Lai,et al.  An upstream flux‐splitting finite‐volume scheme for 2D shallow water equations , 2005 .

[27]  Linsong Cheng,et al.  Flow Behavior of Viscoelastic Polymer Solution in Porous Media , 2015 .

[28]  Olav Møyner,et al.  Multiscale Simulation of Polymer Flooding with Shear Effects , 2016, Transport in Porous Media.

[29]  A. Jannelli,et al.  Second Order Positive Schemes by means of Flux Limiters for the Advection Equation , 2009 .

[30]  R. Seright,et al.  Injectivity Characteristics of EOR Polymers , 2009 .

[31]  D. Kuzmin,et al.  High-resolution FEM-TVD schemes based on a fully multidimensional flux limiter , 2004 .

[32]  R. Seright.  How Much Polymer Should Be Injected During a Polymer Flood? Review of Previous and Current Practices , 2017 .

[33]  Pietro Tesi,et al.  Methods in Oil Recovery Processes and Reservoir Simulation , 2016 .

[34]  A. Stavland,et al.  A model for non-Newtonian flow in porous media at different flow regimes , 2016, Computational Geosciences.

[35]  P. Tyvand,et al.  Onset of Thermal Convection in a Vertical Porous Cylinder with a Partly Conducting and Partly Penetrative Cylinder Wall , 2011 .

[36]  D. W. Peaceman.  Interpretation of Well-Block Pressures in Numerical Reservoir Simulation(includes associated paper 6988 ) , 1978 .

[37]  Nariman Fathi Najafabadi,et al.  Modeling chemical EOR processes using IMPEC and fully IMPLICIT reservoir simulators , 2009 .

[38]  X. Yue,et al.  Micro-mechanisms of residual oil mobilization by viscoelastic fluids , 2008 .