# The Effect of Rheology of Viscoelastic Polymer on Pressure Transient Response in Near-Wellbore Regions

Authors: Jia Zhang; Shiqing Cheng; Jie Zhan; Qi Han
Journal: Geofluids
 (2021)

---

## Abstract

Viscoelastic polymer solution shows shear thinning behavior at low shear rates and shear thickening behavior at high shear rates in reservoirs. However, models that ignored shear thickening behavior were commonly employed to interpret transient pressure data derived from tested wells in viscoelastic polymer flooding systems; although, viscoelastic polymer solutions show shear thickening behavior in the near-wellbore region due to high shear rate. To better characterize the oilfield with pressure transient analysis in viscoelastic polymer flooding systems, we developed a numerical model that takes into account both shear thinning behavior and shear thickening behavior. A finite volume method was employed to discretize partially differential flow equations in a hybrid grid system including PEBI mesh and Cartesian grid, and the Newton-Raphson method was used to solve the fully implicit nonlinear system. To illustrate the significance of our model, we compared our model with a model that ignores the shear thickening behavior by graphing their solutions on log-log plots. In the flow regime of near-wellbore damage, the pressure derivative computed by our model is distinctly larger than that computed by the model ignoring shear thickening behavior. Furthermore, the effect of shear thickening behavior on pressure derivative differs from that of near-wellbore damage. We then investigated the influence of shear thickening behavior on pressure derivative with different polymer injection rates, injection rates, and permeabilities. The results can provide a benchmark to better estimate near-wellbore damage in viscoelastic polymer flooding systems. Besides, we demonstrated the applicability and accuracy of our model by interpreting transient pressure data from a field case in an oilfield with viscoelastic polymer flooding treatments.


---

## Body

## 1. Introduction

Crude oil plays an important role in contributing to the world’s energy consumption. Polymer flooding is widely employed in oilfields as a mature and economic process by increasing sweep efficiency in oil development [1–5]. In the design and adjustment of the polymer flooding process, basic information of an oilfield is the key to success. Pressure transient analysis is an important method to characterize reservoirs [6–12], which can provide valuable information including permeability, degree of near-wellbore damage (skin factor) [13–20], and wellbore storage constant.To expand pressure transient analysis to polymer flooding, we must overcome several challenges. The most important one is that polymer solution shows either or both shear thinning and shear thickening in porous media. More specifically, shear thinning behavior is that the viscosity of the viscous polymer solution decreases with the shear rate [21–24]. However, for viscoelastic polymer solutions, the viscosity decreases with the shear rate at low shear rates and increases when the shear rate increases at high shear rates [25–28]. Due to the above non-Newtonian nature of the polymer, the well test model used for the water flooding system is no longer suitable for the polymer flooding system.For the viscous polymer, a simplified single-phase model solved by the finite difference method was presented by Poollen and Jargon [29] to investigate the pressure transient behavior of a field test. Then, Ikoku and Ramey [30] and Odeh and Yang [31] developed analytical solutions of single-phase models to introduce new methods of pressure transient analysis for the viscous polymer solution. Further, Lund [32] expanded the pressure transient theory for the non-Newtonian shear-thinning fluids into composite reservoirs, in which he assumed a radial flow system consisting of non-Newtonian and Newtonian fluids with a single-phase model. To model the flow of the shear-thinning polymer solution in a vertically fractured well system, Murtha and Ertekin [33] developed a numerical simulator for single-phase flow. On this basis, Vongvuthipornchai and Raghavan [34] analyzed the pressure transient behavior of fractured wells after the injection of the viscous polymer solution. Mahani et al. [35] combined numerical simulations and analytical pressure transient computations to infer the in situ polymer rheology in a polymer flooding system. Kamal et al. [36] developed analytical solutions of single-phase well test models combining the non-Newtonian shear-thinning fluids and the multicomposite reservoir, which can be used to obtain reservoir properties in various banks and the location of the flood front. However, the above models ignored the effect of concentration gradients and saturation gradients on the pressure transient behavior during polymer flooding. To fill this gap, Zhang and Yao [37] proposed a numerical model for two-phase flow to investigate the transient pressure behavior of wells after the injection of the viscous polymer solution. Jia et al. [38] proposed a numerical model of two-phase flow considering the shear-thinning behavior of polymer solutions to study pressure transient response in polymer flooding systems.Though the above models are satisfying in the well test interpretation for the viscous polymer flooding, the influence of shear thickening behavior on pressure transient behavior is ignored in viscoelastic polymer flooding systems. To fill this gap, Yang et al. [39] and Xie et al. [40] investigated the effect of shear-thickening behavior on the pressure response of injection wells during polymer injection with single-phase models. Furthermore, Ma and McClure [41] discovered that the shear thickening behavior influences the early part of the shut-in pressure of polymer injection wells in a viscoelastic polymer flooding system with a single-phase model. However, they did not further analyze the impact of shear thickening behavior on transient pressure response in well test interpretation.Therefore, the impact of shear thickening behavior on the pressure transient response still needs to be discussed in well test interpretation. Ignoring the impact of shear thickening behavior can lead to errors in the diagnosis and interpretation of the transient pressure data, which are obtained from the falloff tests in viscoelastic polymer flooding systems. To solve this problem, a numerical model for two-phase is urgently needed to investigate the transient pressure behavior of falloff tests in viscoelastic polymer flooding systems.The objective of this paper is to better characterize the oilfield with pressure transient analysis by investigating the effect of shear thickening behavior on the pressure transient response of falloff tests in viscoelastic polymer flooding systems. Because partially hydrolyzed polyacrylamide (HPAM) is the most commonly employed viscoelastic polymer in polymer flooding, we analyzed the effect of shear thickening behavior of HPAM on the pressure transient response with a developed two-phase numerical model. To validate the applicability and accuracy of our model, we interpreted the transient pressure data obtained from a polymer injection well in an oilfield.


## 2. Model Development

To investigate the transient pressure response in a viscoelastic polymer flooding system, we developed a numerical model for two-phase flow. To solve our model more efficiently, the flow equations are discretized with a finite volume method in a hybrid grid system including PEBI mesh and Cartesian grid. Then, we applied a fully implicit Newton-Raphson method to solve the discrete flow equations.


### 2.1. Model Description

To model the flow of the partially hydrolyzed polyacrylamide (HPAM), water, and oil in polymer flooding systems, the following basic assumptions are made:(1)
The shear thickening behavior, shear thinning behavior, convection, and adsorption of the HPAM are considered. Diffusion and dispersion are ignored due to their insignificant effect on the pressure transient response. We also consider the permeability reduction of the aqueous phase derived from the polymer adsorption and the inaccessible pore volume due to the large molecular size and flexible structure of the HPAM(2)
The reservoir is isotropic and isothermal; the initial pressure, initial water saturation, and initial polymer concentration are distributed uniformly in the reservoir. The exterior boundary of the reservoir is closed(3)
The reservoir fluids contain water and oil. The HPAM dissolves in the aqueous phase, and so we consider that the effective viscosity of the aqueous phase equals to the viscosity of polymer solution. No chemical reaction takes place during polymer flooding(4)
The relative permeability of the oil phase and aqueous phase is considered. The capillary pressure is neglected because of its small impact on the transient pressure behavior in polymer flooding systems(5)
The flow of the oil phase and aqueous phase in the reservoir follows the Darcy’s law. Fluids and rocks are kind of compressible


### 2.2. Mathematical Model

On the basis of the above assumptions, we propose the partially differential flow equations of the oil, water, and polymer.(1)∂∂tϕbαSα+∇⋅bαvα−Qα=0,α=o,w,(2)∂∂tϕfpbwSwCp+C∧p+∇⋅bwvwCp−QwCp=0,where bα is the reciprocal of the oil or aqueous formation volume factor, fp is the ratio of the accessible pore volume to the pore volume of reservoir rocks, Cp is the polymer concentration of the aqueous phase, and C∧p is the adsorbed concentration of polymer on the rock surface.Given that the relative permeability and permeability reduction are considered, the velocities of oil and aqueous can be replaced by the following equations.(3)vo=−Kkroμo∇po,vw=−KkrwμweffRk∇pw,where μweff is the effective viscosity of aqueous phase considering the dissolved polymer, and Rk is the permeability reduction factor for aqueous phase.The adsorbed concentration of polymer (C∧p) in Equation (2) is expressed by the Langmuir adsorption equation.
(4)Cp∧=C∧pmaxbpCp1+bpCp,where C∧pmax is the saturated adsorbed concentration of polymer, and bp is the constant parameter.The permeability reduction factor for the aqueous phase is described by Hou et al. [42].
(5)Rk=1+Rkmax−1bpCp1+bpCp,where Rkmax is the maximum permeability reduction factor.Although the rates of tested wells are known parameters in well tests, we still require to compute oil or aqueous flow rates of wellbores. Furthermore, the wellbore storage is taken into account in our model for pressure transient analysis, which is ignored in the model of reservoir simulation.(6)Qs=∑Qw+Qo,Qα=2πKλαhbαlnre/rw+Spwf−pα+Cbαdpwfdt,λo=kroμoλw=krwμweffRk.To calculate the effective viscosity of the aqueous phase in viscoelastic polymer flooding systems, we employ a unified viscosity model of polymer solutions proposed by Delshad et al. [43]. The model can characterize the shear thinning behavior and shear thickening behavior of the partially hydrolyzed polyacrylamide (HPAM) solution. The great advantage of Delshad’s model is that it only requires bulk rheology data and petrophysical information of the porous media.
(7)μweff=μw+μp0−μw1+λ1γ2n1−12+μmax1−exp−λ2τrγn2−1.If the shear thickening behavior of polymer solutions is neglected, Equation (7) can be simplified into the following expression.
(8)μweff=μw+μp0−μw1+λ1γ2n1−12,where γ is the shear rate of the polymer solution, μp0 is the viscosity of the polymer solution at zero shear rate, and n1, n2, and λ2 are empirical constants.In the above viscosity model,λ1, μmax, and τr are functions of the polymer concentration.
(9)λ=β1expβ2Cp,μmax=μw1+Ap1Cp+Ap2Cp2Csepsp,τr=τ0+τ1Cp,where β1, β2, Ap1, Ap2, τ0, and τ1 are constant parameters. Csep is effective salinity, and sp is constant parameters.The zero shear viscosity (μp0) in Equation (7) and Equation (8) is a function of the polymer concentration (Cp) and effective salinity (Csep).
(10)μp0=μw1+A1Cp+A2Cp2+A3Cp3Csepsp,where A1, A2, and A3 are input constants.To illustrate the influence of shear thickening behavior on polymer viscosity, we compute the viscosity of polymer solution considering or ignoring the shear thickening behavior (Figure1). The values of the input parameters for calculating the apparent viscosity are listed in Table 1.Figure 1
The viscosity of the viscoelastic polymer solution considering shear thickening behavior and shear thinning behavior and the viscosity of the polymer solution that ignores shear thickening behavior.Table 1
The input constants for calculating the apparent viscosity of the HPAM solution.
ParametersValuesParametersValuesμw, Pa·s0.00099τ0, unitless0.005989n1, unitless0.63τ1, unitless0.1562n2, unitless1.94Csep, kg/m30.33496λ2, unitless0.0088sp, unitless0.7β1, unitless0.0172A1, unitless82.53β2, unitless2.903A2, unitless51.14Ap1, unitless53.387A3, unitless60.07Ap2, unitless4.152According to Equation (7), we can calculate the viscosity of the HPAM solution as a function of the shear rate for different polymer concentrations, as illustrated in Figure 2.Figure 2
The viscosity of the partially hydrolyzed polyacrylamide (HPAM) solution as a function of the shear rate for different polymer concentrations.The shear rate (γ) in Equation (7) and Equation (8) is proposed by Wang [44, 45] based on the model of single-phase flow presented by Cannella et al. [46].
(11)γ=4vw8C′Kkrwϕ/Rk3n+14nnn−1,where C′ and n are constant input parameters.


### 2.3. Differentiation and Solution

To solve the partially differential flow equations, we discretize them in a hybrid grid system including PEBI mesh and Cartesian grid first (Figure3). The near-well zone consists of refined PEBI meshes because the shear thickening behavior of polymer solutions occurs mainly in this zone due to high shear rate. To reduce computation time, the Cartesian grids of large size are distributed far away from the wellbore. In the hybrid system, each grid or mesh contains cell (c), face (f), and node.Figure 3
A hybrid grid system including PEBI mesh and Cartesian grid.In unstructured grids, the finite volume method [47] is easily formulated, and thus we employ it to discretize the above partially differential flow equations. To make discrete equations easy to comprehend, we employ the discrete operators div and grad for the divergence and gradient operators, respectively [48, 49].Thus, the flux passing through the facef can be written as
(12)qf=−Tfgradpf,where q is the array of discrete fluxes, T is the array of discrete transmissibilities through faces, and p is the array of the pressure gradient.Following the above discrete operators, we discretize Equations (1) and 2 with upstream weighting
(13)Rαn+1,l=VcϕcbαcSαcn+1,l-ϕcbαcSαcnΔtn+1+divbαfvαfn+1,l-Qαcn+1,l=0,(14)Rpn+1,l=VcϕcfpbwcSwcCpc+C^pcn+1,l-ϕcfpbwcSwcCpc+C^pcnΔtn+1+divbwfvwfCpfn+1,l-QwcCpcn+1,l=0.Note that we denote arrays of discrete variables with boldfaced letters in the above discrete equations. In Equations (13) and 14, vαf can be substituted by the following expression.
(15)vαf=−TfλαfgradpαfThen, we discretize the well equation.
(16)Rqn+1.l=Qscn+1−∑jQαcn+1,l=0.The flow rate of the oil/aqueous phase in Equation (16) can be written as [49]
(17)Qαn+1,l=2πλαchbαclnre/rw+Sn+1,lpwfc−pαcn+1,l+Cbαcn+1,lpwfn+1,lc−pwfncΔtn+1.Coupling with Equations (13)-(17), we obtain
(18)Rn+1,lXn+1,l=Ron+1,l,Rwn+1,l,Rpn+1,l,Rqn+1,lT=0,where X is the array of primary variables, l is the iterating step, and n is the time step.To solve the above fully implicit discrete equations, the Newton-Raphson method is employed due to its simple formula in programming and fast convergence.With the above methodology, we can obtain the transient pressure of wells to investigate the influence of non-Newtonian behavior on the transient pressure behavior in viscoelastic polymer flooding systems.


## 2.1. Model Description

To model the flow of the partially hydrolyzed polyacrylamide (HPAM), water, and oil in polymer flooding systems, the following basic assumptions are made:(1)
The shear thickening behavior, shear thinning behavior, convection, and adsorption of the HPAM are considered. Diffusion and dispersion are ignored due to their insignificant effect on the pressure transient response. We also consider the permeability reduction of the aqueous phase derived from the polymer adsorption and the inaccessible pore volume due to the large molecular size and flexible structure of the HPAM(2)
The reservoir is isotropic and isothermal; the initial pressure, initial water saturation, and initial polymer concentration are distributed uniformly in the reservoir. The exterior boundary of the reservoir is closed(3)
The reservoir fluids contain water and oil. The HPAM dissolves in the aqueous phase, and so we consider that the effective viscosity of the aqueous phase equals to the viscosity of polymer solution. No chemical reaction takes place during polymer flooding(4)
The relative permeability of the oil phase and aqueous phase is considered. The capillary pressure is neglected because of its small impact on the transient pressure behavior in polymer flooding systems(5)
The flow of the oil phase and aqueous phase in the reservoir follows the Darcy’s law. Fluids and rocks are kind of compressible


## 2.2. Mathematical Model

On the basis of the above assumptions, we propose the partially differential flow equations of the oil, water, and polymer.(1)∂∂tϕbαSα+∇⋅bαvα−Qα=0,α=o,w,(2)∂∂tϕfpbwSwCp+C∧p+∇⋅bwvwCp−QwCp=0,where bα is the reciprocal of the oil or aqueous formation volume factor, fp is the ratio of the accessible pore volume to the pore volume of reservoir rocks, Cp is the polymer concentration of the aqueous phase, and C∧p is the adsorbed concentration of polymer on the rock surface.Given that the relative permeability and permeability reduction are considered, the velocities of oil and aqueous can be replaced by the following equations.(3)vo=−Kkroμo∇po,vw=−KkrwμweffRk∇pw,where μweff is the effective viscosity of aqueous phase considering the dissolved polymer, and Rk is the permeability reduction factor for aqueous phase.The adsorbed concentration of polymer (C∧p) in Equation (2) is expressed by the Langmuir adsorption equation.
(4)Cp∧=C∧pmaxbpCp1+bpCp,where C∧pmax is the saturated adsorbed concentration of polymer, and bp is the constant parameter.The permeability reduction factor for the aqueous phase is described by Hou et al. [42].
(5)Rk=1+Rkmax−1bpCp1+bpCp,where Rkmax is the maximum permeability reduction factor.Although the rates of tested wells are known parameters in well tests, we still require to compute oil or aqueous flow rates of wellbores. Furthermore, the wellbore storage is taken into account in our model for pressure transient analysis, which is ignored in the model of reservoir simulation.(6)Qs=∑Qw+Qo,Qα=2πKλαhbαlnre/rw+Spwf−pα+Cbαdpwfdt,λo=kroμoλw=krwμweffRk.To calculate the effective viscosity of the aqueous phase in viscoelastic polymer flooding systems, we employ a unified viscosity model of polymer solutions proposed by Delshad et al. [43]. The model can characterize the shear thinning behavior and shear thickening behavior of the partially hydrolyzed polyacrylamide (HPAM) solution. The great advantage of Delshad’s model is that it only requires bulk rheology data and petrophysical information of the porous media.
(7)μweff=μw+μp0−μw1+λ1γ2n1−12+μmax1−exp−λ2τrγn2−1.If the shear thickening behavior of polymer solutions is neglected, Equation (7) can be simplified into the following expression.
(8)μweff=μw+μp0−μw1+λ1γ2n1−12,where γ is the shear rate of the polymer solution, μp0 is the viscosity of the polymer solution at zero shear rate, and n1, n2, and λ2 are empirical constants.In the above viscosity model,λ1, μmax, and τr are functions of the polymer concentration.
(9)λ=β1expβ2Cp,μmax=μw1+Ap1Cp+Ap2Cp2Csepsp,τr=τ0+τ1Cp,where β1, β2, Ap1, Ap2, τ0, and τ1 are constant parameters. Csep is effective salinity, and sp is constant parameters.The zero shear viscosity (μp0) in Equation (7) and Equation (8) is a function of the polymer concentration (Cp) and effective salinity (Csep).
(10)μp0=μw1+A1Cp+A2Cp2+A3Cp3Csepsp,where A1, A2, and A3 are input constants.To illustrate the influence of shear thickening behavior on polymer viscosity, we compute the viscosity of polymer solution considering or ignoring the shear thickening behavior (Figure1). The values of the input parameters for calculating the apparent viscosity are listed in Table 1.Figure 1
The viscosity of the viscoelastic polymer solution considering shear thickening behavior and shear thinning behavior and the viscosity of the polymer solution that ignores shear thickening behavior.Table 1
The input constants for calculating the apparent viscosity of the HPAM solution.
ParametersValuesParametersValuesμw, Pa·s0.00099τ0, unitless0.005989n1, unitless0.63τ1, unitless0.1562n2, unitless1.94Csep, kg/m30.33496λ2, unitless0.0088sp, unitless0.7β1, unitless0.0172A1, unitless82.53β2, unitless2.903A2, unitless51.14Ap1, unitless53.387A3, unitless60.07Ap2, unitless4.152According to Equation (7), we can calculate the viscosity of the HPAM solution as a function of the shear rate for different polymer concentrations, as illustrated in Figure 2.Figure 2
The viscosity of the partially hydrolyzed polyacrylamide (HPAM) solution as a function of the shear rate for different polymer concentrations.The shear rate (γ) in Equation (7) and Equation (8) is proposed by Wang [44, 45] based on the model of single-phase flow presented by Cannella et al. [46].
(11)γ=4vw8C′Kkrwϕ/Rk3n+14nnn−1,where C′ and n are constant input parameters.


## 2.3. Differentiation and Solution

To solve the partially differential flow equations, we discretize them in a hybrid grid system including PEBI mesh and Cartesian grid first (Figure3). The near-well zone consists of refined PEBI meshes because the shear thickening behavior of polymer solutions occurs mainly in this zone due to high shear rate. To reduce computation time, the Cartesian grids of large size are distributed far away from the wellbore. In the hybrid system, each grid or mesh contains cell (c), face (f), and node.Figure 3
A hybrid grid system including PEBI mesh and Cartesian grid.In unstructured grids, the finite volume method [47] is easily formulated, and thus we employ it to discretize the above partially differential flow equations. To make discrete equations easy to comprehend, we employ the discrete operators div and grad for the divergence and gradient operators, respectively [48, 49].Thus, the flux passing through the facef can be written as
(12)qf=−Tfgradpf,where q is the array of discrete fluxes, T is the array of discrete transmissibilities through faces, and p is the array of the pressure gradient.Following the above discrete operators, we discretize Equations (1) and 2 with upstream weighting
(13)Rαn+1,l=VcϕcbαcSαcn+1,l-ϕcbαcSαcnΔtn+1+divbαfvαfn+1,l-Qαcn+1,l=0,(14)Rpn+1,l=VcϕcfpbwcSwcCpc+C^pcn+1,l-ϕcfpbwcSwcCpc+C^pcnΔtn+1+divbwfvwfCpfn+1,l-QwcCpcn+1,l=0.Note that we denote arrays of discrete variables with boldfaced letters in the above discrete equations. In Equations (13) and 14, vαf can be substituted by the following expression.
(15)vαf=−TfλαfgradpαfThen, we discretize the well equation.
(16)Rqn+1.l=Qscn+1−∑jQαcn+1,l=0.The flow rate of the oil/aqueous phase in Equation (16) can be written as [49]
(17)Qαn+1,l=2πλαchbαclnre/rw+Sn+1,lpwfc−pαcn+1,l+Cbαcn+1,lpwfn+1,lc−pwfncΔtn+1.Coupling with Equations (13)-(17), we obtain
(18)Rn+1,lXn+1,l=Ron+1,l,Rwn+1,l,Rpn+1,l,Rqn+1,lT=0,where X is the array of primary variables, l is the iterating step, and n is the time step.To solve the above fully implicit discrete equations, the Newton-Raphson method is employed due to its simple formula in programming and fast convergence.With the above methodology, we can obtain the transient pressure of wells to investigate the influence of non-Newtonian behavior on the transient pressure behavior in viscoelastic polymer flooding systems.


## 3. Results and Discussion

Firstly, we introduced type curves of the transient pressure response for a falloff test in a viscoelastic polymer flooding system. Then, we compared our model with a model that ignores the influence of the shear thickening behavior. After that, we crafted four scenarios to analyze the influence of shear thickening behavior on the pressure transient response (i.e., near-wellbore damage, polymer injection concentration, injection rate, and permeability).


### 3.1. Type Curve

Figure4(a) shows the type curve computed by our model, which considers the influence of the shear thickening behavior and shear thinning behavior of the partially hydrolyzed polyacrylamide (HPAM) solution. The type curve is split into five flow regimes. The first one is the flow regime of wellbore storage (I); the flow regime of near-wellbore damage (II) is significantly influenced by the shear thickening behavior of the HPAM solution. In the flow regime of near-wellbore damage, the pressure and pressure derivative computed by our model are obviously larger than those computed by the model only considering the shear thinning behavior of the HPAM solution, as shown in Figure 4(b). This phenomenon is generated from the increase of the aqueous phase viscosity near the wellbore due to shear thickening behavior of the HPAM solution. Then, the duration of the first radial flow regime (III) is ruled by the range of high concentrate polymer solution due to polymer injection. After that, the flow regime of the transitional zone (IV) is dependent on the concentration gradients generated from polymer injection. Subsequently, the system radial flow regime (V) is dominated by the original fluids in the oil field. The pressure derivative of the system radial flow regime is much smaller than that of the first radial flow regime because the fluid viscosity of the near-well zone increases due to polymer injection.Figure 4
Pressure curves of our model and a model ignoring the influence of the shear thickening behavior. (a) Pressure and pressure derivative curves computed by our model. (b) Pressure curves of our model and a model ignoring the influence of the shear thickening behavior.
(a)(b)


### 3.2. Shear Thickening Behavior

To analyze the impact of shear thickening behavior on the pressure transient response of polymer injection wells during falloff tests, four parameters are analyzed. Considering that the influence of shear thickening behavior concentrates on the flow regime of near-wellbore damage, we investigated the impact of the shear thickening behavior on the pressure transient response with different skin factors. Given that the polymer concentration and shear rate have great influences on the polymer rheology (Figure2), we also analyzed the influence of the shear thickening behavior with different polymer injection concentrations, injection rates, and permeabilities. In these scenarios, a polymer injection well is situated in a 6000m×6000m×10m reservoir model with PEBI mesh and regular Cartesian grid. Considering the computing accuracy and computational time, we designate that the radius of each wellbore grid (the minimum PEBI mesh) is 1.24 m, and the size of Cartesian grid is 50m×50m×10m. The minimum time step size at the beginning of polymer injection or shut-in is 0.1162 s, and then the time step size exponentially increases to 17117 s. Polymer injection time (t1) is 374.5 d, and shut-in time (t2) of the polymer injection well is 51 d during a falloff test. The input parameters are listed in Table 1 and Table 2. The relative permeability curves are illustrated in Figure 5.Table 2
The basic parameters of the reservoir.
ParametersValuesInitial reservoir pressure, MPa20Well radius, m0.07Wellbore storage constant, m3/MPa0.1Porosity, fraction0.3Rock compressibility, MPa-10.001Initial water saturation, fraction0.65Initial oil saturation, fraction0.35Oil viscosity, mPa·s7.41Oil compressibility, MPa-10.001Oil volume factor, m3/m31Water compressibility, MPa-10.0002Water volume factor, m3/m31Initial polymer concentration, kg/m30Saturated adsorbed concentration of polymer, kg/m30.00005Constant parameter (bp), unitless2.27Maximum permeability reduction factor, unitless1.84Figure 5
The relative permeability curves.To illustrate the difference between viscoelastic polymer and viscous polymer during polymer injection, we also draw the viscosity distribution around the polymer injection well when polymer injection time is 18.38 h. The polymer injection rate is 100 m3/d, the permeability (K) is 1D, and the polymer injection concentration (Cp) is 1 kg/m3. As shown in Figure 6, the viscosity of viscoelastic polymer is significantly larger than that of viscous polymer in the near-wellbore region.Figure 6
The viscosity distribution around the polymer injection well. (a) The viscosity distribution around the polymer injection well of viscoelastic polymer. (b) The viscosity distribution around the polymer injection well of viscous polymer.
(a)(b)


#### 3.2.1. Near-Wellbore Damage

To research the influence of shear thickening behavior on the pressure transient response when degrees of near-wellbore damage are different, our model and the model that ignores the influence of shear thickening behavior are employed to compute the pressure and pressure derivatives with different skin factors (S). In this scenario, the injection rate of the tested well (q) prior to shut-in is 100 m3/d. The permeability (K) is 1D, and the polymer injection concentration (Cp) is 1 kg/m3. As illustrated in Figure 7, in the flow regime of near-wellbore damage, the pressure derivative computed by our model is always larger than that computed by the model that ignores the shear thickening behavior when the degree of near-wellbore damage changes. The duration of the flow regime of the near-wellbore damage keeps constant when considering the effect of the shear thickening behavior. However, the pressure derivative increases, and duration expands in the flow regime of near-wellbore damage as the skin factor increases. Therefore, the effect of the shear thickening behavior on the pressure transient response cannot be offset by changing the skin factor in well test interpretation.Figure 7
Effect of shear thickening behavior on pressure derivative with different skin factors during falloff tests. (S is the skin factor to represent the degree of the near-wellbore damage. Pressure derivative computed by our model and the model that ignores the shear thickening behavior of the HPAM solution. The solid line is the pressure derivative computed by our model. The dash line is the pressure derivative computed by the model that ignores the shear thickening behavior.)


#### 3.2.2. Polymer Injection Concentration

Given that shear thickening behavior is closely related to polymer concentration, we analyzed the influence of shear thickening behavior on the pressure transient behavior with different polymer injection concentrations. Our model and the model that neglects the influence of shear thickening behavior are applied to compute the pressure transient response with different polymer injection concentrations (Cp). In this scenario, the injection rate of the tested well (q) prior to shut-in is 100 m3/d. The permeability (K) is 1D, and the skin factor (S) is 0.1. Figure 8 illustrates that the duration of the flow regime affected by the shear thickening behavior increases when the polymer injection concentration increases.Figure 8
Impact of shear thickening behavior on pressure derivative with different polymer injection concentrations during falloff tests. The solid line is the pressure derivative computed by our model. The dash line is the pressure derivative computed by the model that ignores the shear thickening behavior.


#### 3.2.3. Injection Rate

Since polymer viscosity is a function of the shear rate, the impact of shear thickening behavior on the pressure transient behavior with different injection rates (q) is investigated. In this scenario, the polymer injection concentration (Cp) is 1 kg/m3, the permeability (K) is 1D, and the skin factor (S) is 0.1. During the flow regime of near-wellbore damage, the influence of the shear thickening behavior on the pressure derivative is more significant when the injection rate grows, as shown in Figure 9.Figure 9
Effect of shear thickening behavior on pressure derivative with different injection rates during falloff tests. The solid line is the pressure derivative computed by our model. The dash line is the pressure derivative computed by the model that ignores the shear thickening behavior of the HPAM solution.


#### 3.2.4. Permeability

Because the shear rate partly depends on the permeability, we also analyzed the influence of shear thickening behavior on the pressure transient behavior with different permeabilities. Our model and the model that neglects the influence of shear thickening behavior are applied to compute the pressure derivatives with different permeabilities (K). The polymer injection concentration (Cp) is 1 kg/m3. The injection rate of the tested well (q) prior to shut-in is 100 m3/d, and the skin factor (S) is 0.1. In the flow regime of near-wellbore damage, the influence of the shear thickening on the pressure derivative is more obvious when the permeability increases, as shown in Figure 10.Figure 10
Effect of shear thickening behavior on pressure derivative with different permeabilities during falloff tests. The solid line is the pressure derivative computed by our model. The dash line is the pressure derivative computed by the model that ignores the shear thickening behavior.


## 3.1. Type Curve

Figure4(a) shows the type curve computed by our model, which considers the influence of the shear thickening behavior and shear thinning behavior of the partially hydrolyzed polyacrylamide (HPAM) solution. The type curve is split into five flow regimes. The first one is the flow regime of wellbore storage (I); the flow regime of near-wellbore damage (II) is significantly influenced by the shear thickening behavior of the HPAM solution. In the flow regime of near-wellbore damage, the pressure and pressure derivative computed by our model are obviously larger than those computed by the model only considering the shear thinning behavior of the HPAM solution, as shown in Figure 4(b). This phenomenon is generated from the increase of the aqueous phase viscosity near the wellbore due to shear thickening behavior of the HPAM solution. Then, the duration of the first radial flow regime (III) is ruled by the range of high concentrate polymer solution due to polymer injection. After that, the flow regime of the transitional zone (IV) is dependent on the concentration gradients generated from polymer injection. Subsequently, the system radial flow regime (V) is dominated by the original fluids in the oil field. The pressure derivative of the system radial flow regime is much smaller than that of the first radial flow regime because the fluid viscosity of the near-well zone increases due to polymer injection.Figure 4
Pressure curves of our model and a model ignoring the influence of the shear thickening behavior. (a) Pressure and pressure derivative curves computed by our model. (b) Pressure curves of our model and a model ignoring the influence of the shear thickening behavior.
(a)(b)


## 3.2. Shear Thickening Behavior

To analyze the impact of shear thickening behavior on the pressure transient response of polymer injection wells during falloff tests, four parameters are analyzed. Considering that the influence of shear thickening behavior concentrates on the flow regime of near-wellbore damage, we investigated the impact of the shear thickening behavior on the pressure transient response with different skin factors. Given that the polymer concentration and shear rate have great influences on the polymer rheology (Figure2), we also analyzed the influence of the shear thickening behavior with different polymer injection concentrations, injection rates, and permeabilities. In these scenarios, a polymer injection well is situated in a 6000m×6000m×10m reservoir model with PEBI mesh and regular Cartesian grid. Considering the computing accuracy and computational time, we designate that the radius of each wellbore grid (the minimum PEBI mesh) is 1.24 m, and the size of Cartesian grid is 50m×50m×10m. The minimum time step size at the beginning of polymer injection or shut-in is 0.1162 s, and then the time step size exponentially increases to 17117 s. Polymer injection time (t1) is 374.5 d, and shut-in time (t2) of the polymer injection well is 51 d during a falloff test. The input parameters are listed in Table 1 and Table 2. The relative permeability curves are illustrated in Figure 5.Table 2
The basic parameters of the reservoir.
ParametersValuesInitial reservoir pressure, MPa20Well radius, m0.07Wellbore storage constant, m3/MPa0.1Porosity, fraction0.3Rock compressibility, MPa-10.001Initial water saturation, fraction0.65Initial oil saturation, fraction0.35Oil viscosity, mPa·s7.41Oil compressibility, MPa-10.001Oil volume factor, m3/m31Water compressibility, MPa-10.0002Water volume factor, m3/m31Initial polymer concentration, kg/m30Saturated adsorbed concentration of polymer, kg/m30.00005Constant parameter (bp), unitless2.27Maximum permeability reduction factor, unitless1.84Figure 5
The relative permeability curves.To illustrate the difference between viscoelastic polymer and viscous polymer during polymer injection, we also draw the viscosity distribution around the polymer injection well when polymer injection time is 18.38 h. The polymer injection rate is 100 m3/d, the permeability (K) is 1D, and the polymer injection concentration (Cp) is 1 kg/m3. As shown in Figure 6, the viscosity of viscoelastic polymer is significantly larger than that of viscous polymer in the near-wellbore region.Figure 6
The viscosity distribution around the polymer injection well. (a) The viscosity distribution around the polymer injection well of viscoelastic polymer. (b) The viscosity distribution around the polymer injection well of viscous polymer.
(a)(b)


### 3.2.1. Near-Wellbore Damage

To research the influence of shear thickening behavior on the pressure transient response when degrees of near-wellbore damage are different, our model and the model that ignores the influence of shear thickening behavior are employed to compute the pressure and pressure derivatives with different skin factors (S). In this scenario, the injection rate of the tested well (q) prior to shut-in is 100 m3/d. The permeability (K) is 1D, and the polymer injection concentration (Cp) is 1 kg/m3. As illustrated in Figure 7, in the flow regime of near-wellbore damage, the pressure derivative computed by our model is always larger than that computed by the model that ignores the shear thickening behavior when the degree of near-wellbore damage changes. The duration of the flow regime of the near-wellbore damage keeps constant when considering the effect of the shear thickening behavior. However, the pressure derivative increases, and duration expands in the flow regime of near-wellbore damage as the skin factor increases. Therefore, the effect of the shear thickening behavior on the pressure transient response cannot be offset by changing the skin factor in well test interpretation.Figure 7
Effect of shear thickening behavior on pressure derivative with different skin factors during falloff tests. (S is the skin factor to represent the degree of the near-wellbore damage. Pressure derivative computed by our model and the model that ignores the shear thickening behavior of the HPAM solution. The solid line is the pressure derivative computed by our model. The dash line is the pressure derivative computed by the model that ignores the shear thickening behavior.)


### 3.2.2. Polymer Injection Concentration

Given that shear thickening behavior is closely related to polymer concentration, we analyzed the influence of shear thickening behavior on the pressure transient behavior with different polymer injection concentrations. Our model and the model that neglects the influence of shear thickening behavior are applied to compute the pressure transient response with different polymer injection concentrations (Cp). In this scenario, the injection rate of the tested well (q) prior to shut-in is 100 m3/d. The permeability (K) is 1D, and the skin factor (S) is 0.1. Figure 8 illustrates that the duration of the flow regime affected by the shear thickening behavior increases when the polymer injection concentration increases.Figure 8
Impact of shear thickening behavior on pressure derivative with different polymer injection concentrations during falloff tests. The solid line is the pressure derivative computed by our model. The dash line is the pressure derivative computed by the model that ignores the shear thickening behavior.


### 3.2.3. Injection Rate

Since polymer viscosity is a function of the shear rate, the impact of shear thickening behavior on the pressure transient behavior with different injection rates (q) is investigated. In this scenario, the polymer injection concentration (Cp) is 1 kg/m3, the permeability (K) is 1D, and the skin factor (S) is 0.1. During the flow regime of near-wellbore damage, the influence of the shear thickening behavior on the pressure derivative is more significant when the injection rate grows, as shown in Figure 9.Figure 9
Effect of shear thickening behavior on pressure derivative with different injection rates during falloff tests. The solid line is the pressure derivative computed by our model. The dash line is the pressure derivative computed by the model that ignores the shear thickening behavior of the HPAM solution.


### 3.2.4. Permeability

Because the shear rate partly depends on the permeability, we also analyzed the influence of shear thickening behavior on the pressure transient behavior with different permeabilities. Our model and the model that neglects the influence of shear thickening behavior are applied to compute the pressure derivatives with different permeabilities (K). The polymer injection concentration (Cp) is 1 kg/m3. The injection rate of the tested well (q) prior to shut-in is 100 m3/d, and the skin factor (S) is 0.1. In the flow regime of near-wellbore damage, the influence of the shear thickening on the pressure derivative is more obvious when the permeability increases, as shown in Figure 10.Figure 10
Effect of shear thickening behavior on pressure derivative with different permeabilities during falloff tests. The solid line is the pressure derivative computed by our model. The dash line is the pressure derivative computed by the model that ignores the shear thickening behavior.


## 3.2.1. Near-Wellbore Damage

To research the influence of shear thickening behavior on the pressure transient response when degrees of near-wellbore damage are different, our model and the model that ignores the influence of shear thickening behavior are employed to compute the pressure and pressure derivatives with different skin factors (S). In this scenario, the injection rate of the tested well (q) prior to shut-in is 100 m3/d. The permeability (K) is 1D, and the polymer injection concentration (Cp) is 1 kg/m3. As illustrated in Figure 7, in the flow regime of near-wellbore damage, the pressure derivative computed by our model is always larger than that computed by the model that ignores the shear thickening behavior when the degree of near-wellbore damage changes. The duration of the flow regime of the near-wellbore damage keeps constant when considering the effect of the shear thickening behavior. However, the pressure derivative increases, and duration expands in the flow regime of near-wellbore damage as the skin factor increases. Therefore, the effect of the shear thickening behavior on the pressure transient response cannot be offset by changing the skin factor in well test interpretation.Figure 7
Effect of shear thickening behavior on pressure derivative with different skin factors during falloff tests. (S is the skin factor to represent the degree of the near-wellbore damage. Pressure derivative computed by our model and the model that ignores the shear thickening behavior of the HPAM solution. The solid line is the pressure derivative computed by our model. The dash line is the pressure derivative computed by the model that ignores the shear thickening behavior.)


## 3.2.2. Polymer Injection Concentration

Given that shear thickening behavior is closely related to polymer concentration, we analyzed the influence of shear thickening behavior on the pressure transient behavior with different polymer injection concentrations. Our model and the model that neglects the influence of shear thickening behavior are applied to compute the pressure transient response with different polymer injection concentrations (Cp). In this scenario, the injection rate of the tested well (q) prior to shut-in is 100 m3/d. The permeability (K) is 1D, and the skin factor (S) is 0.1. Figure 8 illustrates that the duration of the flow regime affected by the shear thickening behavior increases when the polymer injection concentration increases.Figure 8
Impact of shear thickening behavior on pressure derivative with different polymer injection concentrations during falloff tests. The solid line is the pressure derivative computed by our model. The dash line is the pressure derivative computed by the model that ignores the shear thickening behavior.


## 3.2.3. Injection Rate

Since polymer viscosity is a function of the shear rate, the impact of shear thickening behavior on the pressure transient behavior with different injection rates (q) is investigated. In this scenario, the polymer injection concentration (Cp) is 1 kg/m3, the permeability (K) is 1D, and the skin factor (S) is 0.1. During the flow regime of near-wellbore damage, the influence of the shear thickening behavior on the pressure derivative is more significant when the injection rate grows, as shown in Figure 9.Figure 9
Effect of shear thickening behavior on pressure derivative with different injection rates during falloff tests. The solid line is the pressure derivative computed by our model. The dash line is the pressure derivative computed by the model that ignores the shear thickening behavior of the HPAM solution.


## 3.2.4. Permeability

Because the shear rate partly depends on the permeability, we also analyzed the influence of shear thickening behavior on the pressure transient behavior with different permeabilities. Our model and the model that neglects the influence of shear thickening behavior are applied to compute the pressure derivatives with different permeabilities (K). The polymer injection concentration (Cp) is 1 kg/m3. The injection rate of the tested well (q) prior to shut-in is 100 m3/d, and the skin factor (S) is 0.1. In the flow regime of near-wellbore damage, the influence of the shear thickening on the pressure derivative is more obvious when the permeability increases, as shown in Figure 10.Figure 10
Effect of shear thickening behavior on pressure derivative with different permeabilities during falloff tests. The solid line is the pressure derivative computed by our model. The dash line is the pressure derivative computed by the model that ignores the shear thickening behavior.


## 4. Field Application

An appropriate model for history matching is the key to attain satisfactory interpretation results for a well test on the basis of geologic information and features of the type curve. However, the shear thickening behavior of viscoelastic polymers was usually ignored by researchers. Now, it has been demonstrated that the shear thickening behavior of viscoelastic polymer solutions cannot be ignored in pressure transient analysis. To validate the accuracy and practicability of our model in viscoelastic polymer flooding systems, we interpreted the transient pressure data obtained from a tested injection well in a reservoir with HPAM polymer flooding, in which we employed our model and the mode that ignores the shear thickening behavior.


### 4.1. Basic Information of a Falloff Test in a Reservoir with HPAM Polymer Flooding

The basic parameters of a falloff test in a reservoir with HPAM polymer flooding are illustrated in Table3. The relative permeability curves are shown in Figure 5. The input parameters of the HPAM solution are exhibited in Tables 1 and 2. The polymer flooding started in October 2014. The date of falloff test is October 14, 2015. The injection time prior to shut-in is 177 d. The effective shut-in time is 71.46 hr.Table 3
The basic parameters of a field case.
ParametersValuesPorosity, fraction0.174Rock compressibility, MPa-11.46×10−3Thickness, m13Initial water saturation, fraction0.649Water viscosity, mPa·s0.99Oil viscosity, mPa·s7.41Water compressibility, MPa-11.5×10−4Oil compressibility, MPa-12×10−3Water volume factor, m3/m31.001Oil volume factor, m3/m31.174Well radius, m0.07Injection rate of tested well, m3/d60Polymer injection concentration, kg/m30.503Polymer concentration of adjacent production well, kg/m30


### 4.2. Well Test Interpretation

In this section, our model and the model that ignores the influence of shear thickening behavior are employed to interpret the transient pressure data derived from a falloff test. We outline our matching methodology in the following procedure.Step 1.
Input the parameters listed in Section4.1 and transient pressure data from the tested well into the computation procedure on the basis of our model.Step 2.
To reduce the matching time, the limitations of certain unknown parameters can be set according to the geologic and dynamic information of the field case.Step 3.
Regulate the unknown variables based on the results obtained from the above discussion results. More specifically, the features of the flow regimes of the near-wellbore damage and transitional zone are significant and identifiable, which can be used to assist with parameters adjustment in history matching.As illustrated in Figure11(a), the pressure and pressure derivative computed by our model agree well with the transient pressure data obtained from the tested well, which validates the applicability and accuracy of our model. The interpretation results are listed in Table 4. However, the pressure curves computed by the model ignoring the shear thickening behavior match the pressure data poorly in the flow regime of near-wellbore damage, as shown in Figure 11(b). Furthermore, Table 4 shows that the skin factor interpreted by the model ignoring shear thickening is obviously larger than that interpreted by our model, which can lead to a false impression on the near-wellbore damage.Figure 11
History matching results computed by our model and the model that ignores the influence of shear thickening behavior. (a) History matching using our model that considers the shear thickening behavior and shear thinning behavior of the HPAM solution. (b) History matching using the model that ignores the influence of shear thickening behavior.
(a)(b)Table 4
Interpretation results.
Interpreting parametersValues interpreted by our modelValues interpreted by the model ignoring shear thickeningWellbore storage constant, m3/MPa0.9170.917Skin factor, unitless-0.412.06Permeability,D0.420.42


## 4.1. Basic Information of a Falloff Test in a Reservoir with HPAM Polymer Flooding

The basic parameters of a falloff test in a reservoir with HPAM polymer flooding are illustrated in Table3. The relative permeability curves are shown in Figure 5. The input parameters of the HPAM solution are exhibited in Tables 1 and 2. The polymer flooding started in October 2014. The date of falloff test is October 14, 2015. The injection time prior to shut-in is 177 d. The effective shut-in time is 71.46 hr.Table 3
The basic parameters of a field case.
ParametersValuesPorosity, fraction0.174Rock compressibility, MPa-11.46×10−3Thickness, m13Initial water saturation, fraction0.649Water viscosity, mPa·s0.99Oil viscosity, mPa·s7.41Water compressibility, MPa-11.5×10−4Oil compressibility, MPa-12×10−3Water volume factor, m3/m31.001Oil volume factor, m3/m31.174Well radius, m0.07Injection rate of tested well, m3/d60Polymer injection concentration, kg/m30.503Polymer concentration of adjacent production well, kg/m30


## 4.2. Well Test Interpretation

In this section, our model and the model that ignores the influence of shear thickening behavior are employed to interpret the transient pressure data derived from a falloff test. We outline our matching methodology in the following procedure.Step 1.
Input the parameters listed in Section4.1 and transient pressure data from the tested well into the computation procedure on the basis of our model.Step 2.
To reduce the matching time, the limitations of certain unknown parameters can be set according to the geologic and dynamic information of the field case.Step 3.
Regulate the unknown variables based on the results obtained from the above discussion results. More specifically, the features of the flow regimes of the near-wellbore damage and transitional zone are significant and identifiable, which can be used to assist with parameters adjustment in history matching.As illustrated in Figure11(a), the pressure and pressure derivative computed by our model agree well with the transient pressure data obtained from the tested well, which validates the applicability and accuracy of our model. The interpretation results are listed in Table 4. However, the pressure curves computed by the model ignoring the shear thickening behavior match the pressure data poorly in the flow regime of near-wellbore damage, as shown in Figure 11(b). Furthermore, Table 4 shows that the skin factor interpreted by the model ignoring shear thickening is obviously larger than that interpreted by our model, which can lead to a false impression on the near-wellbore damage.Figure 11
History matching results computed by our model and the model that ignores the influence of shear thickening behavior. (a) History matching using our model that considers the shear thickening behavior and shear thinning behavior of the HPAM solution. (b) History matching using the model that ignores the influence of shear thickening behavior.
(a)(b)Table 4
Interpretation results.
Interpreting parametersValues interpreted by our modelValues interpreted by the model ignoring shear thickeningWellbore storage constant, m3/MPa0.9170.917Skin factor, unitless-0.412.06Permeability,D0.420.42


## 5. Conclusions

Although viscoelastic polymer solutions show shear thickening behavior in the near-wellbore region due to high shear rate, the effect of shear thickening behavior on the pressure transient response was commonly ignored in well test interpretation. To better characterize the oilfields with pressure transient analysis in viscoelastic polymer flooding systems, we developed a numerical mode that takes into account both shear thickening behavior and shear thinning behavior. The finite volume method and fully implicit Newton-Raphson method were employed to solve the model in a hybrid grid system. The investigation of our model led to the following conclusions:(1)
The influence of shear thickening behavior on pressure derivative is different from that of near-wellbore damage. In the flow regime of near-wellbore damage, the pressure derivative computed by our model is obviously larger than that computed by the model ignoring shear thickening behavior. The duration of the flow regime of the near-wellbore damage remains constant when considering the effect of the shear thickening behavior(2)
The flow regime of near-wellbore damage lasts longer when the polymer injection concentration increases. The influence of the shear thickening on the pressure derivative is more significant when the injection rate or permeability increases(3)
The skin factor interpreted by the model that ignores the effect of shear thickening behavior is much larger than that interpreted by our model, which can lead to false impression on the degree of near-wellbore damage



---

*Source: 5568336-2021-06-08.xml*