Polymer Flow in Porous Media : Relevance to Enhanced Oil Recovery

Polymer flooding is one of the most successful chemical EOR (enhanced oil recovery) methods, and is primarily implemented to accelerate oil production by sweep improvement. However, additional benefits have extended the utility of polymer flooding. During the last decade, it has been evaluated for use in an increasing number of fields, both offshore and onshore. This is a consequence of (1) improved polymer properties, which extend their use to HTHS (high temperature high salinity) conditions and (2) increased understanding of flow mechanisms such as those for heavy oil mobilization. A key requirement for studying polymer performance is the control and prediction of in-situ porous medium rheology. The first part of this paper reviews recent developments in polymer flow in porous medium, with a focus on polymer in-situ rheology and injectivity. The second part of this paper reports polymer flow experiments conducted using the most widely applied polymer for EOR processes, HPAM (partially hydrolyzed polyacrylamide). The experiments addressed highrate, near-wellbore behavior (radial flow), reservoir rate steady-state flow (linear flow) and the differences observed in terms of flow conditions. In addition, the impact of oil on polymer rheology was investigated and compared to single-phase polymer flow in Bentheimer sandstone rock material. Results show that the presence of oil leads to a reduction in apparent viscosity.


Introduction
The success of polymer flooding depends on the ability of injected solutions to transport polymer molecules deep into a reservoir, thus providing enhanced mobility ratio conditions for the displacement process.In the following sections, we focus on the principal parameters that are crucial in the decision-making process for designing a satisfactory polymer flood design.
The application of polymer flooding to tertiary oil recovery may induce high injection pressures, resulting in injectivity impairment.Since the volumetric injection rate during polymer flooding is constrained by formation fracture pressure, project economics may be significantly affected.Thus, injectivity is a critical parameter and key risk factor for implementation of polymer flood projects.
A large number of injectivity studies, both theoretical and experimental, have been performed in porous media during recent decades, albeit they were mainly studies of linear cores in the absence of residual oil [1][2][3][4][5][6][7].Recently, Skauge et al. [8] performed radial injectivity experiments showing significant reduction in differential pressure compared to linear core floods.This discrepancy in polymer flow in linear cores compared to that in radial disks is partly explained by the of differing pressure conditions that occur when polymer molecules are exposed to transient and semi-transient pressure conditions in radial disks, as opposed to the steady state conditions experienced in linear core floods.In addition, they observed that the onset of apparent shear thickening occurs at significantly higher flux in radial floods.Based on these results, injectivity was suggested to be underestimated from experiments performed in linear core plugs.However, these experiments were performed in the absence of residual oil.If residual oil has a significant effect on polymer propagation in porous media, experiments performed in its absence will not be able to accurately predict polymer performance.
Experimental studies investigating the effects of residual oil on polymer propagation through porous media have been sparse, although they have generally shown decreasing levels of polymer retention in the presence of residual oil [9,10].
The polymer adsorbs to the rock surface and may also block pores due to polymer size (straining) and flow rate (hydrodynamic retention).In addition, different trapping mechanisms may take place.The polymer retention phenomena influence the flow of polymer in porous media, however, these effects are beyond the scope of this paper.The subject has been reviewed in several other books and papers, e.g., Sorbie [11] and Lake [12].
History matches performed in this study aim to highlight the injectivity of partially hydrolyzed polyacrylamides (HPAMs) in radial disks saturated with residual oil, as these conditions best mimic actual flow conditions in oil reservoirs.Results show that the presence of residual oil reduces the apparent viscosity of HPAM in flow through porous media, thus improving injectivity.These results may facilitate increased implementation of polymer EOR (enhanced oil recovery) projects, as previous projects deemed infeasible may now be economically viable.

In-Situ Rheology
Polymer viscosity as a function of shear rate is usually measured using a rheometer.During the measurement process, polymer solutions are exposed to different shear rates in a stepwise manner.For each shear rate, polymer viscosity is measured after steady state conditions are achieved; at this state, it is referred to as bulk viscosity.However, polymer molecules experience significantly different flow conditions in rheometers compared to porous media.In particular: (I) unlike rheometers, porous media exhibit an inherently complex geometry; (II) phenomena such as mechanical degradation may change rheological properties; (III) although they only demonstrate shear thinning behavior in rheometers, polymer solutions may exhibit apparent shear thickening behavior above a certain critical flow rate; (IV) due to the tortuosity of porous media and existence of several contraction-expansion channels, polymer solutions are exposed to a wide range of shear rates at each flow rate and where extensional viscosity becomes more dominant, resulting in significantly different rheology behavior compared to bulk flow.
To account for these contrasting flow conditions, in-situ viscosity has been suggested to describe the fluid flow behavior of polymer solutions in porous media.In-situ viscosity is a macroscopic parameter that can be calculated using Darcy's law for single-phase non-Newtonian fluids: It is generally measured in core flood experiments as a function of Darcy velocity.Comparison of in-situ and bulk rheology (Figure 1) shows vertical and horizontal shifts between viscosity curves.Vertical shifts may be due to phenomena such as mechanical degradation, while horizontal shifts are due to a conversion factor between in-situ shear rate and Darcy velocity, shown as α.The red line in Figure 1 shows an increase in apparent viscosity, which is due to polymer adsorption.The adsorbed layer of polymer reduces the effective pore size and blocks smaller pores, both leading to increased resistance to flow e.g., as determined by an increase in pressure at a given rate compared to a non-adsorbing situation.In contrast, a reduction in pressure (and therefore, in apparent viscosity) can be observed in the presence of depleted layers (see e.g., Sorbie [11]) which leads to slip effects.
Due to the time-consuming nature of in-situ measurements, there have been several attempts to investigate in-situ rheology, both analytically and numerically.In spite of extensive studies [13][14][15][16][17][18][19][20][21][22], limited success has been achieved to reliably relate in-situ to bulk viscosity based on polymer solution and porous media properties.Most of these models were developed based on analytical solutions of non-Newtonian flow through capillary bundles, which simplifies the complex geometry of porous media.
Colloids Interfaces 2018, 2, x FOR PEER REVIEW 3 of 27 resistance to flow e.g., as determined by an increase in pressure at a given rate compared to a nonadsorbing situation.In contrast, a reduction in pressure (and therefore, in apparent viscosity) can be observed in the presence of depleted layers (see e.g., Sorbie [11]) which leads to slip effects.Due to the time-consuming nature of in-situ measurements, there have been several attempts to investigate in-situ rheology, both analytically and numerically.In spite of extensive studies [13][14][15][16][17][18][19][20][21][22], limited success has been achieved to reliably relate in-situ to bulk viscosity based on polymer solution and porous media properties.Most of these models were developed based on analytical solutions of non-Newtonian flow through capillary bundles, which simplifies the complex geometry of porous media.In the following, the calculation procedure of in-situ viscosity is briefly explained: 1. Analytical solutions for a power-law fluid ( = ) at a given flow rate through a capillary tube with an arbitrary radius (R) can be defined by Equation (2).By comparing Equation (2) with the Poiseuille volumetric flow rate for Newtonian fluids in a tube (Equation (3)), an apparent viscosity and shear rate can be obtained from Equations ( 4) and (5), respectively.2. The analytical equation in a single tube (Equation ( 5)) can be extended to account for real porous media by using the capillary bundle approach [23][24][25].An equivalent radius of a capillary bundle model for porous media with known porosity (φ), permeability (K) and tortuosity (ψ) can be obtained by Equation (6).By calculating the Darcy velocity and substituting the equivalent radius (Equation ( 6)) into Equation ( 5), the apparent shear rate as a function of Darcy velocity can be obtained by Equation (7).In the following, the calculation procedure of in-situ viscosity is briefly explained:

1.
Analytical solutions for a power-law fluid (µ = C .γ n−1 ) at a given flow rate through a capillary tube with an arbitrary radius (R) can be defined by Equation (2).By comparing Equation (2) with the Poiseuille volumetric flow rate for Newtonian fluids in a tube (Equation ( 3)), an apparent viscosity and shear rate can be obtained from Equations ( 4) and (5), respectively.

2.
The analytical equation in a single tube (Equation ( 5)) can be extended to account for real porous media by using the capillary bundle approach [23][24][25].An equivalent radius of a capillary bundle model for porous media with known porosity (φ), permeability (K) and tortuosity (ψ) can be obtained by Equation (6).By calculating the Darcy velocity and substituting the equivalent radius (Equation ( 6)) into Equation ( 5), the apparent shear rate as a function of Darcy velocity can be obtained by Equation (7). .
The above expressions are considered as an analytical basis for calculating apparent viscosity in porous media.Based on Equation ( 7), a simplified linear correlation between apparent shear rate and Darcy velocity is generally suggested, i.e., Equation (8), in which the correction factor (α) is the key factor.Some proposed equations for the correction factor are summarized in Table 1.
By comparing different coefficients, different values for apparent viscosity may be obtained. .
Table 1.Summary of proposed models for correction factor (α).

Model Equation for Correction Factor (α) Description
Analytical solution 4 3n+1 4n n n−1 n is the power index in power-law region Hirasaki and Pope [26] 12 √ 150 3n+1 4n n n−1 n is the power index in power-law region Cannella et al. [16] n is the power index in power-law region, S w is water saturation, β is a constant equal to 6.
Based on the capillary bundle approach, other models were also proposed by Bird et al. [24], Christopher and Middleman [25], and Teeuw and Hesselink [15], in which the modified Blake-Kozeny model is used for power-law fluids (Equation ( 9)) and apparent viscosity is obtained using Equation (10).
Based on the discussion given by Teeuw and Hesselink [15], tortuosity has a dual effect on both shear rate and shear stress calculations.Christopher and Middleman [25] only incorporated tortuosity in shear stress calculations, while Bird et al. [24] incorporated tortuosity into the shear rate term.The various values of β chosen by different authors are summarized in Table 2.

Teeuw and Hesselink [15]
√ 2 Hirasaki and Pope [26] conducted several core flood experiments where permeability was in the range 7-23 mD, porosity in the range 18-20% and residual oil between 20% and 32%.Based on these experiments, they concluded that apparent viscosity could be calculated using the capillary bundle approach and Blake-Kozeny model as follows: where: They also included pore size distribution in their calculations: Sadowski and Bird [16] used the Ellis model to obtain viscosity from the shear rate.The following equations for apparent viscosity were suggested based on the Blake-Kozeny model and capillary bundle approach: 1 In the above expressions, µ 0 , τ 1/2 and n are Ellis model parameters that can be measured in rheometers.By applying these equations, they obtained an acceptable match between experimental and predicted results for low to medium molecular weight polymers.
In summary, none of the proposed models for non-Newtonian fluids in porous media based on the capillary bundle approach are in agreement with all experimental results.Therefore, some known limitations of the capillary bundle approach are noted as follows:

•
It neglects complex features of porous media such as tortuosity and pore size distribution.

•
It assumes unidirectional flow as it neglects interconnectivity between pores.

•
It cannot be representative for flow in an anisotropic medium due to its assumption of unique permeability along propagation direction.

•
It assumes a single radius along bundles with no variation in cross-sectional area.The contraction-expansion feature of non-Newtonian flow in porous media is of high importance, especially when studying extensional viscosity, yield stress and elasticity.

•
It is generally developed based on rheological models in which analytical solutions for velocity profiles are available (e.g., power-law and Ellis model).Analytical solutions for some models (e.g., Carreau model) are quite difficult and the equation for velocity is implicit (Equation ( 10) for the Carreau model) and needs to be solved iteratively.
Duda et al. [27] studied polymer solution rheology inside porous media and reported that experimentally measured pressure drops were greater than those predicted by capillary bundle models, especially at lower values of the Carreau power index.Based on their study, a key reason for underestimating correction factors using the capillary bundle approach is the model's failure to capture either the interconnectivity of pores or non-uniform cross-sections of pore bodies and pore throats (i.e., abrupt contractions and expansions, also known as aspect ratio).
According to the aforementioned limitations of the capillary bundle approach and lack of a universally accepted equation for calculating shear rates in porous media, the application of effective medium theory was eventually suggested.This method was able to remediate certain weaknesses in capillary bundle approach, for example, by incorporating pore interconnectivity and variation in cross-sections.Canella et al. [18] extended this method to account for power-law fluids in porous media.Core floods were conducted using xanthan in the concentration range 300-1600 ppm, rock lithology (sandstone and carbonates) in the permeability range 40-800 mD and various oil residuals (0-29%).Their general assumption was that bulk rheological properties of polymer solutions obey the power-law model, and they suggested the following equation for the relation between shear rate and Darcy velocity based on effective medium theory: .
Canella et al. achieved a satisfactory match with their experimental results by using a constant value of 6 for β, although this value far exceeds correction factors suggested by other researchers [28][29][30].Even though all published results in the literature are not covered by using this correction factor, better agreement between analytical and experimental results was obtained, such as in experiments performed by Teeuw and Hesselink [15] and Gogarty [31].
Canella et al. [18] demonstrated that apparent viscosity depends on both microscopic (connectivity, pore size distribution) and macroscopic properties (permeability, porosity) of porous media.Despite calculation improvements, neither effective medium theory nor the capillary bundle model are able to accurately estimate the correction factor.The great discrepancies in results obtained by the models described above and the wide range of correction factors suggested [17] confirm that a universally accepted model does not yet exist.Insufficiency of these models to predict in-situ viscosity may be attributed to their lack of incorporating time dependence and their use of oversimplified porous media models (e.g., capillary bundle).
To avoid over-simplification of porous media obtained by using the capillary bundle approach, pore network modelling has been suggested.In contrast to the capillary bundle approach, pore network modeling envisages porous media as interconnected bundles with idealized geometries where larger pores (pore bodies) are connected via smaller ones (pore throats).Pore network models have been used by Sorbie et al. [20] to study non-Newtonian fluids that exhibit shear thinning properties; later, several authors studied these phenomena [21,[32][33][34][35]. Using network modeling, Sorbie et al. [20] showed that in connected (2D) networks of porous media, the average shear rate in the network correlates linearly with the flow rate.This result is not obvious and indeed is rather unexpected.Thus, any formula of the form of Equation ( 8) which is linear in U, and has a "shift factor", will do well for shear thinning fluids.The paper also shows that a similar argument holds for extensional flow where the extensional rate in the porous medium correlates linearly with flow rate (U).Lopez et al. [21] applied a pore network model to study non-Newtonian fluids using the same approach as for Newtonian fluids, except that viscosity in each bundle was not assumed to be constant and was considered as a function of pressure drop.Therefore, an iterative approach was suggested to calculate pressure drop and apparent viscosity.Although they obtained satisfactory agreement between analytical and experimental results using this approach, Balhoff and Thompson [34] stated that effects of concentration were neglected, and consequently proposed a new model based on CFD calculations to include effects of concentration in calculating conductivity of pore throats.They used pore network modeling to model shear thinning polymer flow with yield stress within a sand-pack.
Zamani et al. [35] studied the effects of rock microstructures on in-situ rheology using digital rock physics and reported that microscopic properties such as aspect ratio, coordination number and tortuosity may affect deviation of in-situ from bulk rheology.
In some experiments [23,27,31,36], in-situ rheology has been reported to deviate significantly from the behavior in bulk flow, such that in-situ rheology may not be calculated directly from bulk rheology using the previously mentioned models.To achieve this, one may use these approaches assuming that either in-situ rheological properties are different from bulk rheological properties (e.g., Hejri et al. [36]) or that the relationship between apparent shear rate and Darcy velocity is non-linear (e.g., Gogarty [31]).
Calculation of in-situ rheology is a controversial subject.Until now, there has been no direct method to obtain it and, generally it has been measured by performing core floods.However, Skauge et al. [37] observed significantly different in-situ rheology for HPAM in linear compared to radial geometry.This discrepancy might be due to differing pressure regimes and flux conditions experienced by polymer solutions flowing through these inherently different flow geometries.
The problem with in-situ rheology calculations extends beyond finding the appropriate correction/shift factor.It also encompasses predicting the onset of extensional viscosity, which is treated as a separate subject in the following section.

Extensional Viscosity
Several experimental results show that, although polymer solutions (e.g., HPAM) only demonstrate shear thinning behavior in a rheometer, they may exhibit apparent shear thickening behavior above a critical shear rate in porous media (Figure 2) [23,27,31,36].Generally, polymer flow in porous media may be divided into two distinct flow regimes: shear dominant and extensional dominant flow regimes.Since apparent shear thickening occurs in the extensional flow regime, it may also be referred to as extensional viscosity.
Colloids Interfaces 2018, 2, x FOR PEER REVIEW 7 of 27 Hejri et al. [36]) or that the relationship between apparent shear rate and Darcy velocity is non-linear (e.g., Gogarty [31]).Calculation of in-situ rheology is a controversial subject.Until now, there has been no direct method to obtain it and, generally it has been measured by performing core floods.However, Skauge et al. [37] observed significantly different in-situ rheology for HPAM in linear compared to radial geometry.This discrepancy might be due to differing pressure regimes and flux conditions experienced by polymer solutions flowing through these inherently different flow geometries.
The problem with in-situ rheology calculations extends beyond finding the appropriate correction/shift factor.It also encompasses predicting the onset of extensional viscosity, which is treated as a separate subject in the following section.

Extensional Viscosity
Several experimental results show that, although polymer solutions (e.g., HPAM) only demonstrate shear thinning behavior in a rheometer, they may exhibit apparent shear thickening behavior above a critical shear rate in porous media (Figure 2) [23,27,31,36].Generally, polymer flow in porous media may be divided into two distinct flow regimes: shear dominant and extensional dominant flow regimes.Since apparent shear thickening occurs in the extensional flow regime, it may also be referred to as extensional viscosity.Although its source is poorly understood, extensional viscosity is considered one of the principal aspects of polymer flow in porous media due to its influence on injectivity and oil mobilization.This phenomenon was suggested to be a consequence of elastic properties of polymer solutions (elongational dominated [38] or inertia-dominated flow [39]).As a result, extensional viscosity is often used interchangeably with elongational viscosity, shear thickening behavior, viscosity enhancement, dilatant behavior and viscoelasticity.Two different models are generally used to explain this phenomenon, the transient network model [40][41][42] and coil stretch model [43].We adhere to the latter of these models.
Polymer molecules may be envisaged as entangled coils, and when exposed to a flow field, two forces may arise.First, an entropic force that attempts to maintain the existing polymer coil configuration.As coil entanglement increases, higher resistance to deformation is observed.Second the drag force resulting from interactions between solvent fluid and polymer molecules.When shear rate increases beyond a critical rate, molecule configurations change abruptly from coil to stretched states.Therefore, polymer coils start to deform, resulting in anisotropy and stress differences between elongation and compression.Consequently, normal stresses and elastic properties become more dominant.Although its source is poorly understood, extensional viscosity is considered one of the principal aspects of polymer flow in porous media due to its influence on injectivity and oil mobilization.This phenomenon was suggested to be a consequence of elastic properties of polymer solutions (elongational dominated [38] or inertia-dominated flow [39]).As a result, extensional viscosity is often used interchangeably with elongational viscosity, shear thickening behavior, viscosity enhancement, dilatant behavior and viscoelasticity.Two different models are generally used to explain this phenomenon, the transient network model [40][41][42] and coil stretch model [43].We adhere to the latter of these models.
Polymer molecules may be envisaged as entangled coils, and when exposed to a flow field, two forces may arise.First, an entropic force that attempts to maintain the existing polymer coil configuration.As coil entanglement increases, higher resistance to deformation is observed.Second the drag force resulting from interactions between solvent fluid and polymer molecules.When shear rate increases beyond a critical rate, molecule configurations change abruptly from coil to stretched states.Therefore, polymer coils start to deform, resulting in anisotropy and stress differences between elongation and compression.Consequently, normal stresses and elastic properties become more dominant.
Choplin and Sabatie [40] suggested that when polymer molecules are exposed to a simple shear flow at a constant shear rate ( .γ), molecules rotate at a constant angular velocity (ω) proportional to applied shear rate, and in each rotation polymer molecules are stretched and compressed.The time between each rotation can be calculated by Equation (18).
where k is a constant of proportionality, related to viscosity.If t is higher than the Zimm relaxation time, no dilatant behavior occurs.Consequently, the critical shear rate at the onset of dilatant behavior may be calculated based on Zimm relaxation time as follows: Polymer viscosity behavior in extensional flow may be entirely different from its behavior in pure shear flow, i.e., polymer solution may show simultaneous shear thinning and extension thickening behavior.Theoretically, extensional viscosity can be calculated from Equation (21), where N 1 is normal stress difference and . is stretch rate.The relative importance of extensional viscosity and shear viscosity is defined by a dimensionless parameter known as the Trouton ratio (Equation ( 22)), initially proposed by Trouton [44].For non-Newtonian fluids (especially viscoelastic fluids), Tr can reach very large values, such as 10 3 to 10 4 (i.e., when polymer solution demonstrates shear thinning and extension thickening simultaneously).
In Figure 2, the in-situ viscosity of viscoelastic polymers is depicted in both shear and extensional flow regime.At the onset of polymer flow, the generated hydrodynamic force from fluid flow (i.e., drag force) is below the threshold value in terms of overcoming entropic forces.Therefore, polymer configuration persists in a coil shape, and viscosity remains constant and equal to the zero-shear rate viscosity (upper Newtonian plateau).As flow rate increases, polymer molecules are exposed to larger drag forces that disentangle polymer coils and aligns them along the flow direction.This coil alignment reduces resistance to flow (i.e., induces viscosity reduction) and is referred to as shear thinning.When the orientation of polymer molecules is completely aligned, they will start to stretch at increasing flow rates.A change in the deformation of polymer molecules may cause normal stress differences.At low stretch rates ( .), N 1 is very low and by increasing the stretch rate, N 1 dramatically increases.In other words, beyond the critical shear rate ( .γ c ), instead of intramolecular interaction, intermolecular interactions will develop which generate amorphous structures much larger than average polymer chain dimensions [28,45].
Within the extensional flow regime, the apparent viscosity generally reaches a maximum value, subsequently followed by a decreasing viscosity interval.This phenomenon may be interpreted as high viscoelastic stresses causing polymer rupture and chain halving, and it has been reported as being more severe in low-permeability porous media [46].As molecular rupture occurs, new molecular weight distributions emerge (larger molecular weight fractions are distorted) and viscosity behavior of the polymer may be governed by a new molecular weight distribution.
The onset of extensional viscosity-the transition point between shear-dominant and extensional dominant flow-depends on polymer, solvent, and porous media properties.The effects of polymer properties on extensional viscosity can be investigated by using special rheometers that only generate pure extensional flow [47][48][49][50][51][52][53][54][55][56].In the following, the effects of polymer, solvent and porous media properties on the onset of extensional viscosity are explained.

Polymer Concentration
Chauveteau [55] reported that the maximum relaxation time increases with polymer concentration, thus dilatant behavior commences at lower shear rates (Figure 3).He also included the effect of concentration in the expression for Zimm relaxation time, producing Equation (23).

Polymer Concentration
Chauveteau [55] reported that the maximum relaxation time increases with polymer concentration, thus dilatant behavior commences at lower shear rates (Figure 3).He also included the effect of concentration in the expression for Zimm relaxation time, producing Equation ( 23).The effect of concentration on extensional viscosity was also investigated by Lewandowska [56].In contrast to Chauveteau, he reported that dilatant behavior commences at higher shear rates with increasing polymer concentration.He attributed this observation to the higher degree of entanglement as the concentration increases, thus increasing the extent of the shear thinning region.
Briscoe et al. [57] could not identify a consistent trend between polymer concentration and onset of extensional viscosity.They assumed that only a narrow region of polymer concentrations is able to generate apparent shear thickening behavior.Below a critical concentration limit, defined as the critical overlap concentration (C*), few polymer chains are able to form transient networks.At concentrations above C*, the extent of shear thinning may increase and, consequently, the onset of apparent shear thickening may be delayed.This effect was also studied by Dupuis et al. [58], where they observed that the onset of dilatant behavior decreased with polymer concentration.However, rheological behavior above the critical shear rate deviated among different concentration ranges (low: 30-60 ppm; medium: 120-240 ppm; and high: 480-960 ppm).Jiang et al. [59] also confirmed scattered data for the onset of extensional viscosity as function of polymer concentration.Clarke et al. [60] reported that the onset of extensional viscosity is independent of concentration and only depends on molecular weight.

Molecular Weight
The lengths of polymer chains increase with molecular weight, resulting in higher inter-and intramolecular entanglement.Thus, the extent of the shear thinning region increases and, consequently, delay the onset of dilatant behavior [56].However, this explanation directly contradicts the expression for the Zimm relaxation time (Equation ( 23)), where the latter increases with molecular weight and causes critical shear rate to occur at a lower shear rate.
Jiang et al. [59] also studied the effects of molecular weight on the onset of extensional viscosity.They concluded that relaxation time increases with molecular weight, thus the onset of extensional The effect of concentration on extensional viscosity was also investigated by Lewandowska [56].In contrast to Chauveteau, he reported that dilatant behavior commences at higher shear rates with increasing polymer concentration.He attributed this observation to the higher degree of entanglement as the concentration increases, thus increasing the extent of the shear thinning region.
Briscoe et al. [57] could not identify a consistent trend between polymer concentration and onset of extensional viscosity.They assumed that only a narrow region of polymer concentrations is able to generate apparent shear thickening behavior.Below a critical concentration limit, defined as the critical overlap concentration (C*), few polymer chains are able to form transient networks.At concentrations above C*, the extent of shear thinning may increase and, consequently, the onset of apparent shear thickening may be delayed.This effect was also studied by Dupuis et al. [58], where they observed that the onset of dilatant behavior decreased with polymer concentration.However, rheological behavior above the critical shear rate deviated among different concentration ranges (low: 30-60 ppm; medium: 120-240 ppm; and high: 480-960 ppm).Jiang et al. [59] also confirmed scattered data for the onset of extensional viscosity as function of polymer concentration.Clarke et al. [60] reported that the onset of extensional viscosity is independent of concentration and only depends on molecular weight.

Molecular Weight
The lengths of polymer chains increase with molecular weight, resulting in higher inter-and intramolecular entanglement.Thus, the extent of the shear thinning region increases and, consequently, delay the onset of dilatant behavior [56].However, this explanation directly contradicts the expression for the Zimm relaxation time (Equation ( 23)), where the latter increases with molecular weight and causes critical shear rate to occur at a lower shear rate.
Jiang et al. [59] also studied the effects of molecular weight on the onset of extensional viscosity.They concluded that relaxation time increases with molecular weight, thus the onset of extensional viscosity occurs at lower shear rates.In addition, they observed that this trend was not valid above a critical molecular weight.Clarke et al. [60] proposed the following correlation for the dependency of the onset of extensional viscosity on polymer molecular weight:

Salinity Effect
The effect of salinity on polymer rheology may be crucial in some reservoir conditions [11,61,62], and depends on polymer type.For typical EOR polymers (e.g., xanthan, HPAM, or generally non-hydrolyzed polymers), increasing salinity generally reduces coil gyration and hydrodynamic radius.Due to the repulsion between ionic groups in HPAM solutions, increasing salinity compresses the electrical double layer on molecular chains and electrostatic repulsion decreases.In the case of HPAM, the reaction mechanism varies for different metal ions i.e., either monovalent (Na + ) or divalent (Ca 2+ ) cations.In the monovalent case, it may suppress the charge effect and reduce the hydrodynamic radius.In the divalent case, reactions between cations (i.e., Ca 2+ ) can play the role of cross-linkers and influence the conformation and rheological properties of HPAM.In both cases, larger shear rates are required to uncoil polymers and the apparent shear thickening commences at larger shear rates [57,58,63].

Degree of Hydrolysis
When HPAM is dissolved in water, electrostatic repulsion forces cause polymer molecules to expand easily and the shear thinning region is shortened.Therefore, as the degree of hydrolysis increases, the onset of apparent shear thickening decreases [56].

Pressure and Temperature Effect
Although polymers are considered incompressible fluids, they do exhibit some degree of compressibility.Thus, pressure may have an impact on viscosity.By increasing pressure, the free volume between polymer molecules decreases and Brownian motion of polymer chains is inhibited, consequently resulting in viscosity increase of polymer solution.Experimental results [64] indicate that the onset of extensional viscosity decrease significantly with pressure.
The effects of temperature on polymer rheology has also been studied extensively [57,59,65,66] and results show that the critical shear rate and onset of dilatant behavior are retarded with increasing temperature.This behavior may have the following two explanations.Firstly, polymer relaxation time and solvent viscosity should both decrease with increasing temperature, based on Equation (23).Secondly, solvent quality decreases with temperature.By decreasing solvent quality, coil size is reduced, and to compensate for this reduction, a larger shear rate is needed to uncoil and elongate the polymer.Therefore, the onset of extensional viscosity occurs at higher shear rates.

Porous Media Properties
In addition to polymer properties, porous media may also significantly influence the generation of extensional flow, as shown by several experimental [25] and numerical studies [67].Due to variation in cross-sectional area along its propagation path, polymer molecules are forced to accelerate and decelerate.Consequently, they will experience both stretch and shear flow in porous media, and above a critical flow rate, extensional flow will dominate shear flow.
To envisage polymer flow in porous media, the latter may be considered as a simplified contraction-expansion channel.As polymer molecules enter contractions, they will be compressed and stretched.If the flow is below a critical velocity, deformed polymer molecules have sufficient time to return to their original state.Therefore, when polymer solutions enter subsequent contractions, no stress is stored and no additional resistance to flow is observed.However, if polymer relaxation time is high and polymer molecules are not able to return to their equilibrium state between contractions, stress will be stored and accumulated, thus resulting in steep increases in pressure drop and apparent viscosity.This phenomenon can be interpreted as a memory effect of polymer molecules.
Due to the inherent nature of porous media, polymer molecules are sheared near the wall and elongated at the flow axis.Therefore, molecular momentum is transferred by both tangential and normal stress components in porous media.Seeing that polymer molecules are able to rotate in pore space, molecules are not strained and effective viscosity is only controlled by shear.In contrast, if molecules are exposed to strain for sufficient time, molecule deformation plays a major role and effective viscosity will be defined by strain [25,[67][68][69][70][71][72].
To predict the onset of extensional viscosity in porous media, the dimensionless Deborah number is defined as a ratio between the characteristic relaxation time of a fluid (θ f ) and characteristic time of porous media (θ p ), considered as the average time to travel from one pore body to another (Equation ( 25)).In other words, the Deborah number may be interpreted as the ratio between elastic and viscous forces.Based on this expression, the Deborah number is zero for Newtonian fluids and infinity for Hookean elastic solids.
Polymer solutions may have a wide range of molecular weights leading to a large number of relaxation times.Many researchers have used the longest relaxation time as representative of θ f .However, this may cause the overestimation of Deborah numbers at the onset of extensional viscosity.Relaxation times may also be calculated from normal stress differences [73].
Some experimental observations revealed that the onset of extensional viscosity occurs when N De is larger than 0.5 [74].However, the Deborah number is not constant in different experiments and a wide range of values has been reported.Marshall and Metzner [73] reported a Deborah number of 0.1 at the onset of extensional viscosity, while Chauveteau [55] reported a relatively high Deborah number of 10.This wide range of reported Deborah numbers at the onset of extensional viscosity is due to difficulties in calculating stretch rates in porous media.To support this idea, Heemskerk et al. [75] reported that by using different polymer types in the same rock sample, critical Deborah numbers (N De ) were identical.However, when the same polymer was used in different rock samples, the critical N De varied between 1 and 2. They concluded that measured relaxation times from experimental results can be used to practically define the onset of extensional viscosity, but they acknowledged that equations for calculating stretch rate are not able to capture the exact N De at the onset of extensional viscosity.Zamani et al. [67] proposed that to obtain a more accurate estimation of the critical N De , the stretch rate distribution at the pore scale is required.Metzner et al. [76] concluded that the critical Deborah number might only be used as a first estimation of the onset of extensional viscosity.In Table 3, some suggested equations for the calculation of Deborah number are summarized.
They used the inverse of the shear rate for θ p .U w is the Darcy velocity, k rw is the water relative permeability, S w is water saturation and .γ c is a constant equal to 3.97C, where C is an empirical correlation factor to account for the difference between an equivalent capillary model and real porous media Hirasaki and Pope [26] Haas and Durst [78] Heemskerk et al. [75] 1 Several experimental results [68,79] show that the Deborah number alone is not sufficient to predict the onset of extensional viscosity.As an explanation, Ranjbar et al. [80] stated that the onset of extensional viscosity highly depends on the elastic properties of polymer solutions and relaxation time alone cannot capture viscoelastic properties.Experimental results reported by Garrouch and Gharbi [79] support this idea.They investigated two different polymer solutions (xanthan and HPAM) in Berea and sand-packs.Calculated Deborah numbers for these two completely and inherently different polymer solutions inside sand-packs were (surprisingly) identical.While xanthan consists of rigid, rod-like molecules that do not show extensional viscosity, HPAM consists of flexible and elastic chain-structured molecules.
Zamani et al. [67] numerically studied the effect of porous media on the onset of extensional viscosity by using real images of porous media obtained from digital rock physics.They confirmed that microscopic features of porous media had significant impact on the onset of extensional viscosity.Furthermore, by increasing the aspect ratio and inaccessible pore volume and decreasing the coordination number, extensional viscosity occurred at lower shear rates, in agreement with several experimental results [55,68,81].
Skauge et al. [37] reported that in radial flow, the onset of extensional viscosity occurred at higher shear rates than at typical core flooding.Since radial flow is more representative of real field conditions, results obtained from radial disks should be more accurate as laboratory data for field implementation.
Briefly summarized, at low shear rates where the amplitude of the elastic component is negligible, flow is controlled by shear forces.In contrast, above a critical shear rate, flow is extensional and governed by elastic forces.Therefore, the response of polymer solutions to imposed stress may be expressed as the sum of shear and elastic components: The viscosity of polymer solutions under shear flow can be described by empirical equations such as the power-law and Carreau models.To describe viscosity under elongational flow, several models have been suggested, and some of them are summarized in Table 4.
Table 4. Proposed models for calculation of elongational viscosity.

Model Equation Description
Hirasaki and Pope [26] µ el = µ sh τ r is the characteristic relaxation time and can be calculated by dynamic frequency sweep test in the laboratory.Some empirical correlations are also proposed for dependency of different parameters on polymer concentration

Stavland et al. [62]
m is a non-zero tuning parameter which is known as the elongation exponent and depends on the molecular weight and demonstrates linear correlation with [µ] C p .α in the listed formulation is considered 2.5

Injectivity
Polymer injectivity is a crucial factor governing the economics of polymer flooding projects and its accurate estimation is a prerequisite in terms of optimizing the upper-limit injection rate [82].Injection well pressure may increase due to one of the following causes: (1) oil bank formation, (2) in-situ polymer viscosity (especially shear thickening due to viscoelasticity) and (3) different types of retention, which cause permeability reduction.
The highest pressure drops observed during polymer flooding are located in the vicinity of the injection wellbore due to dramatic variations in flow rate.Therefore, it is important to include non-Newtonian effects of polymer solutions to accurately predict polymer injectivity.Although both HPAM and xanthan demonstrate shear thinning behavior at low to moderate shear rates, HPAM exhibits apparent shear thickening above a critical flow rate due to its inherent viscoelastic nature.For field applications, injection rates in the vicinity of the injection well may easily exceed the onset of extensional viscosity, and injectivity will then dramatically decrease.In contrast to HPAM, xanthan shows exclusively shear thinning behavior and will attain its highest value of injectivity in the near-wellbore region.
Injectivity investigations at the lab scale are required before implementing field applications, and effects of polymer solution properties, in-situ rheology, temperature, pH, level of retention and the nature of porous media should be accurately measured [83,84].Furthermore, if screening criteria for polymer type are disregarded, polymer entrapment in narrow pore throats can have significant effects on its injection rate.The salinity of solutions can also affect polymer solubility, resulting in filter cake formation near injection wells or precipitation of polymer molecules in the reservoir.Inaccurate measurement of in-situ rheology and especially the onset of extensional viscosity may lead to either an underestimation or overestimation of injectivity.In some polymer flooding projects, measured injectivity may differ significantly from the simulation or analytical forecast.These unexpected injectivities may be due to the occurrence of mechanical degradation [82,85,86], induced fractures [87][88][89], or even inaccurate analytical models for calculating in-situ rheology and predicting extensional viscosity.

Radial In-Situ Rheology
Injectivity (I) may be defined as the ratio of volumetric injection rate, Q, to the pressure drop, ∆P, associated with polymer propagation between injection well and producer [1]: As previously mentioned, formation fracture pressure may constrain the value of volumetric injection rate.Due to its significant effect on project economics, accurate determination of differential pressure, and hence injectivity, at a given injection rate is essential.To achieve this, all factors affecting differential pressure during polymer flooding must be quantified.Darcy's law for radial flow may be expressed in terms of differential pressure as follows: where µ app is apparent viscosity, h is disk thickness, k e,i is effective permeability to polymer solution, r e is disk radius and r w is injection well radius.In this paper, the ratio of resistance factor (RF) to residual resistance factor (RRF) is used to represent apparent viscosity of polymer solutions propagating through porous media, thus isolating its viscous behavior, i.e., where the resistance factor (RF) represents the pressure increase of polymer relative to brine and the residual resistance factor (RRF) is defined as the ratio of pressure before and after polymer injection (i.e., pressure caused by irreversible permeability reduction induced by retention mechanisms).Due to their inherent viscoelastic behavior in porous media, synthetic polymers (e.g., HPAM) will exhibit shear-dependent apparent viscosity.Although the common consensus on apparent shear thickening as a phenomenon is accepted, its viscosifying magnitude is still an ongoing topic of debate in scientific communities.
Accurate polymer rheology estimation is a prerequisite for reasonable injectivity estimates due to the proportionality between apparent viscosity and differential pressure.In linear core floods where steady-state pressure conditions exist, polymer flux will remain constant from inlet to outlet, rendering rheology estimation a straightforward task.However, in radial flow, polymer flux is gradually reduced as it propagates from injection well to producer, therefore attaining a range of viscosities rather than one specific value.Since the degree of mechanical degradation generally increases with injection rate, discrepancies in polymer rheology obtained from different injection rates may transpire.Instead of possessing one definite rheology, polymers propagating through radial disks will exhibit both shear-dependent and history-dependent viscosity behavior, thus increasing the complexity of rheology estimation in radial compared to linear models.To date, no correction factor has been suggested to account for this dual nature phenomenon.Even when mechanical degradation is excluded, i.e., when injected and effluent viscosities are approximately equal, this dual nature phenomenon persists, and is suggested to be attributed to non-equilibrium pressure conditions experienced in radial flow and inherent history-dependent nature of polymer molecules.
In addition, synthetic polymers are susceptible to mechanical degradation at high flux, typically in the near-wellbore region, which will impart an irreversible viscosity reduction due to polymer molecule fragmentation.Mechanical degradation induces a pressure drop that improves injectivity.However, since it disrupts the carefully selected viscous properties of the polymer solutions by a non-reversible viscosity decrease, mechanical degradation is not a sought-after phenomenon in polymer flooding.A remediation measure to reduce mechanical degradation is to pre-shear the polymer before injection.Pre-shearing removes the high molecular weight part of the molecular weight distribution, which is believed to be most susceptible to mechanical degradation [6].Mechanical degradation may also be minimized by shifting to a lower molecular weight polymer.However, this would require higher amounts of polymer to obtain the same concentration, thus potentially influencing polymer project economics.
As mentioned, in radial geometry, high flux causing mechanical degradation occurs principally in the near-wellbore region, as opposed to linear geometry where this high flux persists throughout the entire propagation distance.Therefore, the time that polymer is exposed to high shear is short in radial transient flow pattern, as opposed to that of a steady-state linear core flood, [34].Based on this time-differing condition between linear and radial flow, it was suggested by Skauge et al. that polymer is degraded to a lesser extent in radial compared to linear flood when injected at the same volumetric flow rate [33].
In summary, there are two principal factors governing injectivity during polymer flooding in linear geometry: (1) viscoelasticity of polymer that induces large injection pressures mainly due to apparent shear thickening behavior at high flux; and (2) mechanical degradation in the near-wellbore region, which causes an entrance pressure drop [1].In radial disks, two additional factors should be included: (3) non-equilibrium pressure conditions due to kinetic effects; and (4) memory-effects of polymer molecules in non-constant velocity fields.
Polymer: Flopaam 3630S, 30% hydrolyzed, MW = 18 million Da.Brine: Relatively low salinity with a low content of divalent ions.Brine composition by ions is given in Table 5. Linear core floods: Core data are summarized in Table 6.All experiments were performed at room temperature and pressure.
Radial core floods: Bentheimer disks were prepared by coating with epoxy resin, vacuuming and saturating with brine.One disk was then drained with an extra heavy oil and aged for 3 weeks at 50 • C to a non-water-wet state.The crude had an initial viscosity of about 7000 cP.The extra heavy oil used for drainage and aging, was then exchanged with a flooding oil of 210 cP.Both experiments were performed at room temperature and pressure.Core data are given in Table 7.The pressure ports were located in the injection and production wells and at radii 0.5, 1.0, 1.5, 2.0, 3.0, 4.0 and 5.0 cm for the disk without oil and at radii 1.1, 2.0, and 5.0 cm for the disk containing oil.
The Bentheimer cores show a pore-throat distribution function similar to other outcrop sandstone material, Figure 4.All porous media have local pore-size variation, involving continuous contraction and expansion of pore-scale transport.Linear core floods: Core data are summarized in Table 6.All experiments were performed at room temperature and pressure.
Radial core floods: Bentheimer disks were prepared by coating with epoxy resin, vacuuming and saturating with brine.One disk was then drained with an extra heavy oil and aged for 3 weeks at 50 °C to a non-water-wet state.The crude had an initial viscosity of about 7000 cP.The extra heavy oil used for drainage and aging, was then exchanged with a flooding oil of 210 cP.Both experiments were performed at room temperature and pressure.Core data are given in Table 7.The pressure ports were located in the injection and production wells and at radii 0.5, 1.0, 1.5, 2.0, 3.0, 4.0 and 5.0 cm for the disk without oil and at radii 1.1, 2.0, and 5.0 cm for the disk containing oil.
The Bentheimer cores show a pore-throat distribution function similar to other outcrop sandstone material, Figure 4.All porous media have local pore-size variation, involving continuous contraction and expansion of pore-scale transport.Simulation: The experimental set-up enabled detailed monitoring of pressure by internal pressure ports located at various distance from the injection well.Differential pressure as function of radial distance was history matched using the STARS simulator, developed by Computer Modeling Group (CMG).The simulation model encompassed a radial grid with 360 sectors, each consisting of 150 grid block cells in radial direction, where the grid block cell size is 1 mm.Porous media

Simulation:
The experimental set-up enabled detailed monitoring of pressure by internal pressure ports located at various distance from the injection well.Differential pressure as function of radial distance was history matched using the STARS simulator, developed by Computer Modeling Group (CMG).The simulation model encompassed a radial grid with 360 sectors, each consisting of 150 grid block cells in radial direction, where the grid block cell size is 1 mm.Porous media permeability (tuning parameter) was obtained by history matching water floods prior to polymer flooding.Local permeability variation improved the history match compared to analytical solution (Darcy's law for radial flow).Permeability data obtained from water floods were used in subsequent polymer floods to isolate the effects of polymer apparent viscosity on differential pressure.In polymer floods, as the permeability obtained from the precursor water flood was held constant, apparent viscosity could be quantitatively investigated as a function of velocity and was used as the tuning parameter to history match differential pressure.The STARS simulation tool can include both shear thinning and thickening behavior of viscoelastic fluids.
Due to the inherent grid averaging calculation method of the simulation tool, the velocity in the first grid block after the injection well was below its analytical value.Because of a rapid velocity decrease with distance in radial models, this phenomenon was addressed by decreasing the injection well radius, thus effectively parallel shifting the position of the first grid block towards the injection well until the correct velocity was attained.This was a necessary step, since the tuning parameter is apparent viscosity as a function of velocity.

Polymer In-Situ Rheology in Linear Cores
Four Bentheimer outcrop cores were used to study polymer in-situ rheology in linear systems.Petro-physical properties of core samples as well as properties of polymer solutions are given in Table 6.Two experiments were carried out to examine the effect of polymer concentration on in-situ rheology of the polymer solution.Partially-hydrolyzed Flopaam 3630S at 500 ppm and 1500 ppm was injected into the cores and the in-situ rheology of the polymer solutions was measured.The two concentrations were chosen to give viscosities representative of the upper and lower limit of what would be economically viable for polymer flooding in an oil field.Both concentrations are above the polymer critical overlap concentration, C*.The results are presented in Table 6 and Figure 5.The bulk viscosity of 1500 ppm 3630S is about 34 cP which is about 5 times that of 500 ppm 3630S.Comparing in-situ rheology of 500 ppm and 1500 ppm 3630S shows that the onset and degree of apparent shear thickening behavior are fairly similar for both concentrations.This is in line with observations by Skauge et al. [8] and Clarke et al. [60] that the onset of extensional viscosity is independent of polymer concentration and only depends on polymer molecular weight.It is noted that this is generally only true for C* < C < C lim , where C lim is the economic limit for polymer concentration, typically between 1500 and 2500 ppm.Table 6 and Figure 5 show that the magnitude of resistance factor (RF) and residual resistance factor (RRF) are about 4 and 3 times higher for 1500 ppm compared to 500 ppm, respectively.This implies that polymer injectivity is a function of polymer concentration, and better injectivity is achieved with lower polymer concentrations.
A series of experiments was also performed to study the effect of the presence of residual oil on polymer in-situ rheology.In these experiments, Bentheimer cores at residual oil saturations of about 22% and different initial wettability states were flooded with polymer and the in-situ rheology behavior was compared to that of single-phase polymer injection in absence of residual oil.Prior to polymer injection, the cores containing oil were water flooded to residual oil saturation.At the end of the water flood, the flow rates were increased to generate pressures higher than that expected for the subsequent polymer flood.This was performed in order to avoid oil mobilization during the polymer flood and, indeed, no oil production was observed during the subsequent polymer flood.The results are presented and compared in Figure 6.As this figure shows, the onset of apparent shear thickening is not affected by the presence of residual oil or the wettability state of the cores.However, the slope of apparent shear thickening and magnitude of resistance factor is significantly affected by oil presence in the cores.That is, although onset of apparent shear thickening is independent of oil presence in porous media and its wettability condition, the results show that the degree of apparent shear thickening is lower when oil is present in the porous media.the subsequent polymer flood.This was performed in order to avoid oil mobilization during the polymer flood and, indeed, no oil production was observed during the subsequent polymer flood.
The results are presented and compared in Figure 6.As this figure shows, the onset of apparent shear thickening is not affected by the presence of residual oil or the wettability state of the cores.However, the slope of apparent shear thickening and magnitude of resistance factor is significantly affected by oil presence in the cores.That is, although onset of apparent shear thickening is independent of oil presence in porous media and its wettability condition, the results show that the degree of apparent shear thickening is lower when oil is present in the porous media.It is important to note that a lower resistance factor in the presence of oil is achieved while porous media is partially occupied by residual oil.Therefore, unlike the single-phase system, in which the pore volume (assuming no inaccessible pore volume) is available for polymer flow, only PV*(1-Sor) is available for polymer flow in two-phase system.This influences and reduces permeability and therefore an even higher resistance factor is expected in the presence of oil.However, the results do not show such an effect, and a lower resistance factor and polymer injectivity is observed with the presence of oil in porous media, which supports the significance of the positive effect of oil on polymer injectivity.The effluent polymer viscosity is reduced by 18% compared to the injected polymer solution for the single-phase, water-wet case, while there is no reduction in effluent viscosity for the two-phase experiments (water-wet and non-water-wet).Shiran and Skauge [90] studied wettability using the same crude oil for aging and found that intermediate wettability was achieved.The end-point water relative permeability confirms that a similar condition was obtained.
Polymer injection in cores with residual oil results in a lower resistance factor which means better polymer injectivity.Furthermore, the resistance factor in the aged core with the non-water-wet It is important to note that a lower resistance factor in the presence of oil is achieved while porous media is partially occupied by residual oil.Therefore, unlike the single-phase system, in which the pore volume (assuming no inaccessible pore volume) is available for polymer flow, only PV*(1-Sor) is available for polymer flow in two-phase system.This influences and reduces permeability and therefore an even higher resistance factor is expected in the presence of oil.However, the results do not show such an effect, and a lower resistance factor and polymer injectivity is observed with the presence of oil in porous media, which supports the significance of the positive effect of oil on polymer injectivity.The effluent polymer viscosity is reduced by 18% compared to the injected polymer solution for the single-phase, water-wet case, while there is no reduction in effluent viscosity for the two-phase experiments (water-wet and non-water-wet).Shiran and Skauge [90] studied wettability using the same crude oil for aging and found that intermediate wettability was achieved.The end-point water relative permeability confirms that a similar condition was obtained.
Polymer injection in cores with residual oil results in a lower resistance factor which means better polymer injectivity.Furthermore, the resistance factor in the aged core with the non-water-wet state is lower than the resistance factor in the water-wet core.The lower resistance factor in the presence of oil could be attributed to lower adsorption/retention of polymer molecules on rock surface, as reported by Broseta et al. [10].The rock surface in the presence of an oil film, and especially in less water-wet conditions is partially covered by crude oil polar components during flooding.Therefore, in comparison to single-phase systems, the rock surface has fewer adsorption sites to adsorb polymer molecules.The analysis of reduced apparent viscosity in the presence of oil, assumes that end-point water relative permeability remains constant for polymer as it does for water.The RRF measured with brine after the polymer injection is assumed constant for all rate variation of polymer flow.Under these assumptions a lower resistance factor and better polymer injectivity is expected.

Polymer In-Situ Rheology in Radial Flow
Recently, polymer injectivity was analyzed by matching field injectivity tests [5,6,91,92].In addition to history matching, modification of equations to incorporate fractures and polymer degradation in the near-wellbore zone were reported.The laboratory experiment simplifies the analysis as additional complications like fractures and strong heterogeneity can be avoided.
In earlier studies of radial flow experiments, Skauge et al. [37] used local pressure taps as a function of radial distance from the injection well to derive in-situ rheology.These experiments demonstrated both shear thinning and strong apparent shear thickening behavior.
Two radial flow experiments were performed on circular Bentheimer sandstone disks of 1.6 and 2.6 D permeability with 30 cm diameter and 2.2 cm thickness, see Table 7.The first experiment was performed on a disk that was drained with crude oil and aged to non-water-wet conditions.The second experiment was performed in the absence of oil on a water-wet disk.For the first experiment, the disk was flooded extensively with brine to reach residual oil saturation, Sorw = 0.22.Bump rates were applied to avoid oil mobilization by viscous forces during the subsequent polymer flood.The polymer flood was performed by first saturating the disk with polymer at a low rate to avoid mechanical degradation due to shearing.Thereafter, rate variations were performed to determine in-situ rheology of the polymer.A brine flush was performed between concentration slugs to remove non-adsorbed polymer.
Concentrations of 800 and 2000 ppm were chosen to represent lower and upper boundaries of the semi-dilute region.The second experiment included the same steps, except for water flooding to Sorw.In this case, the water flood was performed to obtain a pressure reference for the subsequent polymer flood.No oil production was detected during polymer floods.
Differential pressure was measured by internal pressure ports located at different radii from the injection well.The 800-ppm HPAM solution was injected in the presence of residual oil at flow rates of 2.2 and 2.8 mL/min, and in the absence of residual oil at 2.0 and 4.0 mL/min.Differential pressure decay as a function of radial distance from injector is shown in Figure 7.The pressure transition zone from semi-steady-state to steady-state is extended compared to the case without oil.Most notable is the difference in pressure in the injection well.While differential pressures measured from internal pressure ports are higher for the two-phase system (as expected), well injection pressure is significantly lower in the presence of residual oil.Taking the pressure ratio of pressure ports at ~1 cm from injection well as a reference, injection pressure should be 5-6 times higher for the disk with oil, compared to the one without.Instead, the injection pressure is 25% lower.There may be several reasons for this observed result.One reason may be that the presence of oil reduces the effective pore volume, thereby leading to higher flow velocities for the polymer in the near-well region.This would subsequently lead to higher effective shear forces on the polymer, producing mechanical degradation.If mechanical degradation occurs, it has only a minor effect on the shear viscosity.The shear viscosity is 15.1 mPas for the effluent sample taken at 2.0 mL/min, while it was 16.0 mPas for the injected solution (measured at 22 • C, 10 1/s).However, as discussed in Section 2.2, it is the extensional viscosity that is the determining factor for high pressures in near-well region.Changes in extensional viscosity are intrinsically hard to measure and were not performed here.It is still possible that the increase in shear forces for the case with residual oil lead to a reduction in extensional viscosity but not for the case without oil where the effective pore volume was larger.The two other reasons are related to the wetting state of the porous media.If the oil is located in smaller pores, polymer flow is diverted to larger pores where it flows at higher velocities (higher flux).Since the velocity increase takes place in larger pores, only minor degradation would be expected.A third reason may be that porous media is fractionally oil-wet and that there is a difference in the slip conditions for the water-wet and the oil-wet surfaces.This may reduce effective shear for the oil-wet surfaces leading to reduced mechanical degradation.Although there have been speculations on the "lubricating" effect of oil-wet surfaces, no clear evidence of the effect on apparent viscosity or injectivity for core material have been shown to date.It is not possible to differentiate between the three phenomena based on the pressure data alone.Each of the polymer floods were history matched using STARS (CMG).The measured differential pressures as a function of distance from injection well were used as history match parameters, while polymer apparent viscosity was used as a tuning variable.History matches and polymer rheology from both experiments for 800-ppm HPAM floods are shown in Figures 7 and 8, respectively.It is evident from Figure 8 that the polymer rheology is significantly influenced by the presence of residual oil.In terms of absolute values, the apparent viscosity is between a factor of 5 and 10, and it is higher in the absence, compared to the presence of residual oil.Furthermore, the onset of apparent shear thickening shifts to lower velocities in the presence of residual oil.This occurrence is suggested to result from reduced propagation cross-section caused by the residual oil saturation.When flow channels in porous media become narrower, the extensional flow regime is reached at a lower flux, and HPAM exhibits viscoelastic behavior at an earlier stage, thus the onset of apparent shear thickening commences at a lower flux.The effect of shifting the onset of apparent shear thickening to a lower flux may be detrimental for injectivity.However, since the apparent shear thickening seems to be much more extensive in the absence of residual oil, the rheology shows that overall injectivity is significantly improved in presence of residual oil.Each of the polymer floods were history matched using STARS (CMG).The measured differential pressures as a function of distance from injection well were used as history match History matches and polymer rheology in the presence and absence of residual oil for 2000 ppm floods are shown in Figures 9 and 10, respectively.In order to evaluate the influence of polymer concentration on in-situ rheology, 2000 ppm HPAM was injected in both disks.The differential pressures are shown in Figure 9.In this case, the injection rates were Q = 2.0 and 5.0 mL/min for the disk with no oil, and Q = 1.4 and 1.6 mL/min for the disk with oil.These data show the same trend as for the 800 ppm injection: strong reduction in injection well pressure in the presence of residual oil and extension of the transition zone.
In accordance with the 800 ppm solution, polymer viscosity was significantly higher in the absence compared to presence of residual oil, and ranged between a factor of 6 and 16 in their joint velocity interval, Figure 10.In addition, the 2000 ppm solution also showed a decrease in the onset of apparent shear thickening in the presence of residual oil, consistent with the lower concentration solution investigated.Similar to the 800 ppm solution, apparent shear thickening is observed to be much more extensive in absence of residual oil, thus improved injectivity in the presence of residual oil is further corroborated.
saturation.When flow channels in porous media become narrower, the extensional flow regime is reached at a lower flux, and HPAM exhibits viscoelastic behavior at an earlier stage, thus the onset of apparent shear thickening commences at a lower flux.The effect of shifting the onset of apparent shear thickening to a lower flux may be detrimental for injectivity.However, since the apparent shear thickening seems to be much more extensive in the absence of residual oil, the rheology shows that overall injectivity is significantly improved in presence of residual oil.History matches and polymer rheology in the presence and absence of residual oil for 2000 ppm floods are shown in Figures 9 and 10, respectively.In order to evaluate the influence of polymer concentration on in-situ rheology, 2000 ppm HPAM was injected in both disks.The differential pressures are shown in Figure 9.In this case, the injection rates were Q = 2.0 and 5.0 mL/min for the disk with no oil, and Q = 1.4 and 1.6 mL/min for the disk with oil.These data show the same trend as for the 800 ppm injection: strong reduction in injection well pressure in the presence of residual oil and extension of the transition zone.In accordance with the 800 ppm solution, polymer viscosity was significantly higher in the absence compared to presence of residual oil, and ranged between a factor of 6 and 16 in their joint velocity interval, Figure 10.In addition, the 2000 ppm solution also showed a decrease in the onset of apparent shear thickening in the presence of residual oil, consistent with the lower concentration solution investigated.Similar to the 800 ppm solution, apparent shear thickening is observed to be much more extensive in absence of residual oil, thus improved injectivity in the presence of residual oil is further corroborated.

Conclusions
A review of polymer flow in a porous medium was presented.The available EOR analytical models we evaluated have limitations in accurately describing flow of polymer at high shear rates, e.g., near injector, and this leads to underestimating or overestimating of polymer injectivity.
The experimental results presented expand our insight into polymer flow in a porous medium.Shear thinning behavior may be present in core floods while bulk rheology is predominant from rheometer measurements.Linear polymer flow experiments are dominated by apparent shear thickening which is not measured in standard rheometers.The extensional viscosity, which is the main cause of the apparent shear thickening behavior, occurs at flow velocities strongly influenced by the porous media.
Linear core floods are commonly used for evaluating polymer in-situ rheology and injectivity, but they suffer from steady-state conditions throughout the core as opposed to the well injection situation where both pressure and shear forces are nonlinear gradients.
In the linear core floods, the onset of apparent shear thickening is independent of polymer concentration, when polymer type, brine composition and porous media are held constant.It is also independent of the presence of oil and wettability for the three cases evaluated here.
Radial flow injections show more complex in-situ rheology.The in-situ rheology shows a much higher degree of apparent shear thickening in the presence of oil.This may be due to restrictions in the pore space available.In the absence of oil, high concentration polymer (2000 ppm) showed shear thinning behavior.The onset of apparent shear thickening was shifted to higher flow velocities.There is a need for further development of numerical models that incorporate memory effects and possible kinetic effects for high polymer flow rates in the near-well region.
Both linear and radial experiments confirm lower apparent viscosity when oil is present in the porous medium.This conclusion is based on the assumption that brine end-point relative permeability is unchanged for polymer injection compared to two-phase flow by water injection.No extra oil was produced during polymer injection and this support the lowering of in-situ polymer viscosity in the presence of oil.

Figure 1 .
Figure 1.Schematic comparison of in-situ and bulk rheology.

Figure 1 .
Figure 1.Schematic comparison of in-situ and bulk rheology.

Figure 2 .
Figure 2. Schematic illustration of apparent viscosity in porous media.

Figure 2 .
Figure 2. Schematic illustration of apparent viscosity in porous media.

Figure 3 .
Figure 3.Effect of polymer concentration on the onset of extensional viscosity in a model with 45 successive constrictions.Reproduced with permission from [55].

Figure 3 .
Figure 3.Effect of polymer concentration on the onset of extensional viscosity in a model with 45 successive constrictions.Reproduced with permission from [55].

[ 1 −
N De ] Masuda et al. [77] µ elas = µ sh C c (N De ) mc where C c and m c are constant and relate to pore geometry Delshad et al. [61]

Figure 4 .
Figure 4. Mercury injection derived pore throat distribution for Bentheimer core material used in the polymer flow experiments.

Figure 4 .
Figure 4. Mercury injection derived pore throat distribution for Bentheimer core material used in the polymer flow experiments.

Figure 6 .
Figure 6.Resistance factor versus interstitial velocity of 500 ppm pre-filtered Flopaam 3630S partially hydrolyzed polyacrylamide (HPAM) polymer dissolved in 1 wt% NaCl brine, single-phase polymer flow and polymer flow at residual oil saturation.

Figure 6 .
Figure 6.Resistance factor versus interstitial velocity of 500 ppm pre-filtered Flopaam 3630S partially hydrolyzed polyacrylamide (HPAM) polymer dissolved in 1 wt% NaCl brine, single-phase polymer flow and polymer flow at residual oil saturation.

Figure 7 .
Figure 7. Differential pressure profiles for 800-ppm HPAM floods in the presence and absence of residual oil in radial geometry as a function of distance from injector to producer for four flow rates.

Figure 7 .
Figure 7. Differential pressure profiles for 800-ppm HPAM floods in the presence and absence of residual oil in radial geometry as a function of distance from injector to producer for four flow rates.

Figure 7 .
Figure 7. Differential pressure profiles for 800-ppm HPAM floods in the presence and absence of residual oil in radial geometry as a function of distance from injector to producer for four flow rates.

Figure 8 .
Figure 8. Apparent viscosity from history match of differential pressure for 800 ppm HPAM in presence and absence of residual oil in radial geometry.

Figure 8 .
Figure 8. Apparent viscosity from history match of differential pressure for 800 ppm HPAM in presence and absence of residual oil in radial geometry.

Figure 9 .
Figure 9. Differential pressure profiles for 2000 ppm 3630S HPAM floods in presence and absence of residual oil in radial geometry as a function of distance from injector to producer for four flow rates.

Figure 9 .
Figure 9. Differential pressure profiles for 2000 ppm 3630S HPAM floods in presence and absence of residual oil in radial geometry as a function of distance from injector to producer for four flow rates.Colloids Interfaces 2018, 2, x FOR PEER REVIEW 21 of 27

Figure 10 .
Figure 10.Apparent viscosity from the history match of differential pressure for 2000 ppm 3630S HPAM in the presence and absence of residual oil in radial geometry.

Figure 10 .
Figure 10.Apparent viscosity from the history match of differential pressure for 2000 ppm 3630S HPAM in the presence and absence of residual oil in radial geometry.

Table 3 .
Proposed equations for Deborah number calculation.

Table 6 .
Core data for linear core floods.

Table 7 .
Core data for radial core floods.

Table 6 .
Core data for linear core floods.

Table 7 .
Core data for radial core floods.