Research Article  Open Access
Ji Ho Lee, Kun Sang Lee, "Multiphase, Multicomponent Simulation for Flow and Transport during Polymer Flood under Various Wettability Conditions", Journal of Applied Mathematics, vol. 2013, Article ID 101670, 8 pages, 2013. https://doi.org/10.1155/2013/101670
Multiphase, Multicomponent Simulation for Flow and Transport during Polymer Flood under Various Wettability Conditions
Abstract
Accurate assessment of polymer flood requires the understanding of flow and transport of fluids involved in the process under different wettability of reservoirs. Because variations in relative permeability and capillary pressure induced from different wettability control the distribution and flow of fluids in the reservoirs, the performance of polymer flood depends on reservoir wettability. A multiphase, multicomponent reservoir simulator, which covers threedimensional fluid flow and mass transport, is used to investigate the effects of wettability on the flow process during polymer flood. Results of polymer flood are compared with those of waterflood to evaluate how much polymer flood improves the oil recovery and wateroil ratio. When polymer flood is applied to waterwet and oilwet reservoirs, the appearance of influence is delayed for oilwet reservoirs compared with waterwet reservoirs due to unfavorable mobility ratio. In spite of the delay, significant improvement in oil recovery is obtained for oilwet reservoirs. With respect to water production, polymer flood leads to substantial reduction for oilwet reservoirs compared with waterwet reservoirs. Moreover, application of polymer flood for oilwet reservoirs extends productive period which is longer than waterwet reservoir case.
1. Introduction
After primary and secondary oil recovery, there remains lots of oil in place in the reservoirs. To gain unrecovered oil in the reservoirs, additional flood, called enhanced oil recovery (EOR) method, is needed to be applied. Polymer flood has been known as one of the most widely used chemical EOR method. Increased viscosity of displacing fluid by adding low concentrations of water soluble and high molecular weight polymer into water produces the lower mobility ratio. Favorable mobility ratio improves the sweep efficiency so that oil recovery is increased.
Craig defined wettability as the tendency of one fluid to spread on or adhere to a solid surface in the presence of other immiscible fluids [1]. The influences of reservoir wettability on multiphase flow on porous medium and hence on oil recovery have been well known. The wettability of a rock controls the location, flow, and distribution of fluids within reservoir rocks [2], which affects the relative permeability and capillary pressure [3, 4]. Therefore, the recovery of polymer flood in oilwet reservoirs is totally different from that in waterwet reservoirs. Also, considerable hydrocarbon reserves are remained in mixedwet or oilwet reservoirs. Consequently, there is a necessity to better understand the effects of wettability on the oil recovery during polymer flood. Most of polymer flood studies have been carried out without considering wettability effect or focused on waterwet reservoirs with respect to field scale [5–7]. The investigations on the role of wettability on the various aspects of oil recovery during polymer flooding have been mainly based on the experimental works [8–10]. Therefore, there is need to conduct comprehensive numerical study on the flow and transport of aqueous and oleic phases during polymer flood under various wetting conditions on field scale.
2. Mathematical Formulation
Simulations of polymer flood were conducted with UTCHEM, which is a 3D, multicomponent, multiphase, compositional model of chemical flooding processes considering complex phase behavior, chemical and physical transformation, and heterogeneous porous media properties [11].
The basic mass conservation equation for components can be written as follows: where is the component index, is the phase index including aqueous () and oleic () phases, is the porosity, is the overall concentration of component (volume fraction), is the density of component [ML^{3}], is the number of phases, is the concentration of component in phase (volume fraction), is Darcy velocity of phase [LT^{1}], is the total source/sink term for component (volume of component per unit volume of porous media per unit time), is the dispersion tensor.
The phase flux from Darcy’s law is where is the intrinsic permeability tensor, is the vertical depth, is the relative permeability, is the viscosity, and is the specific weight of phase .
To predict the reservoir behavior under multiphase condition, it is important to understand relative permeability of a reservoir rock to each of the fluids flowing through it. Relative permeability is assumed to be solely determined by own saturation and residual saturations and investigated by experimental and analytical methods [12]. Extensive studies have been conducted resulting in representative correlations between relative permeability, saturation, and other factors [12]. In this study, multiphase relative permeabilities are modeled with Coreytype functions [11]. Coreytype relative permeability is expressed with relative permeability on residual saturation, exponent defining the curvature of relative permeability, and residual saturation determining normalized saturation. Coreytype relative permeability equation is given as follows: where is the end point of relative permeability meaning relative permeability at residual saturation, is the exponent of relative permeability of phase determining the curvature of relative permeability, represents the normalized saturation of phase calculated as follows: where is the saturation of phase and is the residual saturation of phase .
Whether oil displaces water or water displaces oil through flow channels in reservoirs, flow phenomenon for immiscible fluid phases, oil and water, is involved with the pressure difference between these phases. Brooks and Corey observed a large number of data on consolidated rock cores and analyzed them by plotting the logarithm of effective saturation versus the logarithm of capillary pressure [13]. Due to this study, linear relationship between logarithm of effective saturation and logarithm of capillary pressure has been revealed and is shown as follows: where is the effective saturation calculated with residual saturations, and are constants realized from intercept and slope, means pore size distribution, is interpreted as maximum capillary pressure, the is the capillary pressure which represents pressure difference between wetting phase pressure and nonwetting phase pressure.
Capillary pressure is strong function of saturation as presented by (5). Leverett derived capillary pressure scaled by soil permeability and porosity for homogeneous reservoirs [11, 14]. Reflected on previous relations, Brooks and Corey capillary pressuresaturation is calculated as follows: where is equivalent to , is constant, and and are the permeability and the porosity of the reservoirs.
3. Numerical Modeling
This study analyzes the effect of wettability on flow and transport of fluids during polymer flood. Reservoir wettability is implemented in the numerical model by changing relative permeability and capillary pressure curves simultaneously. The reservoir depth is 2,000?ft and initial reservoir pressure is maintained at 400?psi. Horizontal area is ?ft^{2} and vertical thickness of reservoir is 25?ft. The simulation domain consists of 10 layers and each layer is discretized into grid blocks in horizontal direction. As close to production and injection wells, grid block has smaller size comparing with blocks which are away from wells to assess pressures and saturations more accurately near wells. The model assumes that the reservoir is homogeneous so the porosity and permeability are constant everywhere as 0.2 and 300?md for horizontal direction and 30?md for vertical direction. The initial saturations of oil and water are assumed to be constant for 0.62 and 0.38, respectively. Properties of water and oil in reservoir are listed in Table 1. Water for injection is assumed to be identical as water in reservoir. Viscosity of polymeric solution including salinity and mechanical effects is considered as listed in Table 1.

The simulation was continued over 1,000 days. Injection designs for water and polymeric solution with conditions of constant rate are assumed to be identical for both waterwet and oilwet reservoirs cases to access the performance of polymer flood quantitatively. To prevent high injection pressure due to increased viscosity of polymeric solution, injection of water and polymeric solution was operated in three steps as given in Figure 1. The flow rate of solution was constant at each step. In the first step, with 0.1% of polymer, injection was operated at low rate as 800?ft^{3}/day to prevent fractures on the injection well until 180 days. In the next step, the concentration of polymer was reduced to 0.05%, whereas the flow rate of injection was increased to 1,000?ft^{3}/day by 360 days. In the last step, from 360 to 1,000 days, only water was injected through well by the end of operation. The reservoir fluids were recovered from the production well constrained with at 200?psi.
In order to analyze the effect of wettability on the flow and transport of fluids during polymer flood, comparisons were made for results from simulations implementing different relative permeability curves and capillary pressure curves. Relative permeability and capillary pressure are the most important factors that show the different performances of polymer flood applied to waterwet and oilwet reservoirs. To model wettability effect, data for different relative permeabilities and capillary pressure curves are attained from several studies. Anderson derived an equation of Coreytype functions by curve fitting the data measured by Morrow et al. [15–17]. Table 1 lists fitted relative permeability parameters and Figure 2(a) shows relative permeability curves generated with the parameters.
(a)
(b)
Calculation of capillary pressure is based on Brooks and Corey equation as previously explained. Capillary pressure end point, , and capillary pressure exponent, , for waterwet and oilwet reservoirs are listed in Table 2 [17]. Capillary pressures are also calculated with the data as shown in Figure 2(b).

4. Results and Discussion
Based on the simulation results, comparisons were made between the performances of polymer flood and waterflood applied to waterwet and oilwet reservoirs. The results were presented with cumulative oil recovery, wateroil ratio, and oil saturation distribution of the fifth (middle) layer. To decide whether application of polymer flood for oilwet reservoirs is effective or not, water cut was also analyzed.
Figures 3(a) and 3(b) present the cumulative oil recovery and wateroil ratio which was obtained from the application of waterflood and polymer flood to waterwet reservoirs. As can be seen, the cumulative oil recovery by polymer flood is 0.58, which is considerably higher than the oil recovery by waterflood, 0.40, at the end of production. In terms of wateroil ratio, polymer flood results in lower wateroil ratio from 70 to 850 days, not from the initiation of injection. Even though polymeric solution is injected with a start of simulation, it takes time, 70 days in this model, to manifest the effect of polymer flood so that no improvement can be seen until 70 days. After 850 days, polymer flood already recovers almost movable oil in reservoirs, so that significant water production and lower increment of oil recovery are obtained. These improvements, higher oil recovery and lower wateroil ratio, during 70 days to 850 days result from increased viscosity of polymeric solution leading to lower or favorable mobility ratio. Furthermore, productive period is calculated to prove the influence of polymer flood applied to waterwet reservoir compared with waterflood case. The life span of production is determined with the assumption that producer is valid until it produces 90% of water cut. Considering results shown in Figure 5(a), polymer flood for waterwet reservoir sustains productive period of 160 days longer than that of waterflood case.
(a)
(b)
In attempts to study the effectiveness of polymer flood in oilwet reservoirs, relative permeability and capillary pressure curves are set to be as Figures 2(a) and 2(b). Figures 4(a) and 4(b) compare the cumulative oil recovery and wateroil ratio between numerical simulations of polymer flood and waterflood for oilwet reservoir. Application of polymer flood to oilwet reservoir increases cumulative recovery as much as 0.36 which is higher than the result of waterflood case, 0.08, and decreases wateroil ratio significantly. According to the calculation of productive period, while production by waterflood seems invalid for whole production period due to higher water cut as much as 90%, polymer flood extends the life span of the well as much as 315 days as shown Figure 5(b). From the analysis of the average reservoir pressure in application of polymer flood shown in Figure 6, pressure profile for oilwet reservoir has been maintained to be lower than waterwet reservoir. Due to the high viscosity of polymeric solution, pressure starts to be increased initially regardless of wettability. After postflush as waterflood is applied, it is going to be decreased.
(a)
(b)
(a)
(b)
Despite of these effective performances of polymer flood for oilwet reservoirs, utilization of polymer flood has not widely fulfilled. As shown in Figures 3 and 4, lower recovery and higher wateroil ratio are observed than those of waterwet reservoirs when both waterflood and polymer flood are applied to oilwet reservoirs. In agreement with previous studies, the performance of polymer flood in oilwet reservoir seems less effective than that in waterwet reservoir. These results could draw the conclusion that the application of polymer flood to oilwet seems not very effective and may not be recommended because it still shows low performance. However, not only the absolute value of oil recovery and wateroil ratio but also the improvement should be analyzed to consider efficiency of polymer flood. Therefore, additional analysis was repeated with respect to improvement by polymer flood. For waterwet reservoirs, polymer flood increases cumulative oil recovery up to 45% than that of waterflood from 0.40 to 0.58. Polymer flood also reduces wateroil ratio at least 50% less than the waterflood case. If there is enough capability for waterflood to produce high performance, these improvements attained by polymer flood may not be so significant. Whereas, for the oilwet reservoirs, the application of polymer flood leads to substantial improvement of oil recovery and extensive reduction of wateroil ratio compared with waterflood case due to low performance of waterflood. As shown in Figure 5(b), when waterflood is applied to oilwet reservoirs, it already exceeds 90% water cut at the early stage and it increases up to 99%. Because of high water production, waterflood should be suspended for oilwet reservoir. However, polymer flood reduces wateroil ratio as much as maximum 90% compared with waterflood case, so that it makes oilwet reservoirs as fruitful. Moreover, polymer flood leads to increment of oil recovery as much as 351% at the end of production as shown in Figure 4. Additionally, analysis on improvements of water cut by polymer flood for different wettability is calculated and shown in Figure 7. Even though polymer flood for waterwet reservoirs represents the maximum reduction of water cut as much as 22%, it results in reduction of as much as 80% for oilwet reservoirs. Also, it sustains longer productive period as much as 315 days compared with that of waterwet reservoirs, 160 days, as shown in Figures 5 and 7. Therefore, the application of polymer flood results in more significant improvement in the case of oilwet reservoir than in the case of waterwet reservoir.
Figure 8 compares remained oil saturation in the fifth layer in application of waterflood and polymer flood to waterwet reservoirs at 1,000 days or 1.98 pore volumes injected. Due to the favorable mobility ratio by polymer flood, relatively higher contrast of oil saturation between swept and unswept regions exits in reservoir. On the other hand, Figure 9 shows simulated oil saturation in the same condition like Figure 8, except for wettability as oilwet reservoir. The overall oil saturation of polymerflooded waterwet reservoirs is lower than that of waterflooded scheme and is almost close to residual oil saturation. For oilwet reservoirs, remained oil saturation is still higher than residual oil saturation after polymer flood is applied. Nevertheless, the contrast of oil saturation distribution between before and after polymer flood is applied is clearly seen in oilwet reservoirs rather than waterwet reservoirs.
(a)
(b)
(a)
(b)
5. Conclusions
According to multiphase, multicomponent simulation, polymer flood allows the mobility ratio between aqueous and oleic fluids to be favorable, which increases oil recovery and decreases wateroil ratio in both waterwet and oilwet reservoir conditions. The efficiency of polymer flood is remarkably affected by reservoir wettability. Performance of polymer flood for waterwet condition seems better than that for oilwet condition. Because polymer flood is kind of modified waterflood, polymeric solution can displace oil out of the pore easily in waterwet scheme so that it sweeps lots of mobile oil and leaves just a little oil including residual oil saturation. This mechanism results in higher recovery for waterwet reservoir than oilwet one. It could mislead us that polymer flood for oilwet reservoirs is not effective as much as that for waterwet reservoirs. If anything, application of polymer flood to oilwet reservoirs clearly shows higher improvement of oil recovery, wateroil ratio, and productive period than that of waterwet reservoirs. These results demonstrate that oilwet reservoirs seem good candidates for the application of polymer flood because the technique is very effective in terms of the improvement of performance. Therefore, reliable evaluation of polymer flood should take into account wettability of reservoirs.
Acknowledgments
This work was supported by the Energy Efficiency and Resources of the Korea Institute of Energy Technology Evaluation and Planning (KETEP) Grant funded by the Korean government, Ministry of Trade, Industry and Energy (20122010300020).
References
 F. F. Craig, The Reservoir Engineering Aspects of Waterflooding, vol. 3 of SPE Monograph Series, Society of Petroleum Engineers, 1971.
 W. G. Anderson, “Wettability literature survey—part 1: rock/oil/brine interactions and the effects of core handling on wettability,” Journal of Petroleum Technology, vol. 38, no. 11, pp. 1125–1144, 1986. View at: Google Scholar
 W. G. Anderson, “Wettability literature survey—part 5: the effects of wettability on relative permeability,” Journal of Petroleum Technology, vol. 39, no. 11, pp. 1453–1468, 1987. View at: Google Scholar
 W. G. Anderson, “Wettability literature survey—part 4: effects of wettability on capillary pressure,” Journal of Petroleum Technology, vol. 39, no. 10, pp. 1283–1300, 1987. View at: Google Scholar
 Y. Du and L. Guan, “Fieldscale polymer flooding: lessons learnt and experiences gained during past 40 years,” in Proceedings of the SPE International Petroleum Conference, pp. 383–388, Puebla, Mexico, November 2004, SPE 91787MS. View at: Google Scholar
 X. Zhang, W. Guan, and F. Pan, “Polymer flooding for middle and low permeability sandstone reservoirs,” in Proceedings of the the 6th International Petroleum Technology Conference, Beijing, China, March 2013, SPE 16536MS. View at: Google Scholar
 Y. Wang, Y. Pang, Z. Shao et al., “The polymer flooding technique applied at high water cut stage in daqing oilfield,” in Proceedings of the the North Africa Technical Conference & Exhibition, Cairo, Egypt, April 2013, SPE 164595. View at: Google Scholar
 H. Xia, D. Wang, J. Wu, and F. Kong, “Elasticity of HPAM solutions increases displacement efficiency under mixed wettability conditions,” in Proceedings of the SPE Asia Pacific Oil and Gas Conference and Exhibition, pp. 103–110, Perth, Australia, October 2004, SPE 88456. View at: Google Scholar
 H. Dong, Y. Hong, W. Rui, and D. Fan, “The effect of wettability on oil recovery of alkaline/surfactant/polymer flooding,” in Proceedings of the SPE Annual Technical Conference and Exhibition, pp. 2396–2403, San Antonio, Tex, USA, September 2006, SPE 102564. View at: Google Scholar
 Y. Zhu, Y. Zhang, Q. Hou, H. Yuan, and G. Jian, “Effect of main factors on oil recovery of surfactantpolymer flooding,” in Proceedings of the 6th International Petroleum Technology Conference (IPTC '13), Beijing, China, March 2013, IPTC16433. View at: Google Scholar
 Center for Petroleum and Geosystems Engineering, UTCHEM9.0, a ThreeDimensional Chemical Flood Simulator, University of Texas at Austin, Austin, Tex, USA, 2000.
 A. T. Corey, “The interrelation between gas and oil relative permeabilities,” Production Monthly, vol. 19, no. 1, pp. 38–41, 1954. View at: Google Scholar
 R. H. Brooks and A. T. Corey, “Properties of porous media affecting fluid flow,” Journal of the Irrigation and Drainage Division, vol. 92, no. 2, pp. 61–90, 1966. View at: Google Scholar
 M. C. Leverett, “Capillary behavior on porous solids,” Transactions of the AIME, vol. 142, no. 1, pp. 152–169, 1941. View at: Google Scholar
 G. A. Anderson, Simulation of chemical flood enhanced oil recovery processes including the effects of reservoir wettability [M.S. thesis], University of Texas at Austin, Austin, Tex, USA, 2006.
 N. R. Morrow, P. J. Cram, and F. G. McCaffery, “Displacement studies in dolomite with wettability control by octanoic acid,” SPE Journal, vol. 13, no. 4, pp. 221–232, 1973. View at: Google Scholar
 M. Delshad, N. F. Najafabadi, G. A. Anderson, G. A. Pope, and K. Sepehrnoori, “Modeling wettability alteration in naturally fractured reservoirs,” in Proceedings of the 15th SPEDOE Improved Oil Recovery Symposium, pp. 1132–1140, Tulsa, Okla, USA, April 2006, SPE 100081. View at: Google Scholar
Copyright
Copyright © 2013 Ji Ho Lee and Kun Sang Lee. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.