Next Article in Journal
Water Absorption and Distribution in a Pultruded Unidirectional Carbon/Glass Hybrid Rod under Hydraulic Pressure and Elevated Temperatures
Next Article in Special Issue
Analysis of the Models of Motion of Aqueous Solutions of Polymers on the Basis of Their Exact Solutions
Previous Article in Journal
Selective and Efficient Arsenic Recovery from Water through Quaternary Amino-Functionalized Silica
Previous Article in Special Issue
Flow Behavior of Chain and Star Polymers and Their Mixtures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Microflow Mechanism of Oil Displacement by Viscoelastic Hydrophobically Associating Water-Soluble Polymers in Enhanced Oil Recovery

1
Key Laboratory for Enhanced Oil & Gas Recovery of the Ministry of Education, Northeast Petroleum University, Daqing 163318, China
2
McDougall School of Petroleum Engineering, University of Tulsa, Tulsa, OK 74104, USA
*
Authors to whom correspondence should be addressed.
Polymers 2018, 10(6), 628; https://doi.org/10.3390/polym10060628
Submission received: 8 May 2018 / Revised: 29 May 2018 / Accepted: 31 May 2018 / Published: 7 June 2018
(This article belongs to the Special Issue Polymer Dynamics)

Abstract

:
Polymer flooding plays an important role in enhanced oil recovery (EOR), particularly in China, where partially hydrolyzed polyacrylamide (HPAM) and hydrophobically associating water-soluble polymers (HAWP) are used in onshore and offshore reservoirs, respectively. Many researchers have highlighted the elasticity of HPAM, which can be used to improve the sweep efficiency, i.e., the ratio of the area swept by an injected fluid to the oil area. On the other hand, fewer studies exist on the elasticity of HAWP. In this study, we investigate the flow of HAWP and Xanthan solutions with identical viscosities in core experiments in terms of elasticity; results reveal that the HAWP can produce shear thickening in the core. The constitutive equation for the HAWP can be obtained using the simulation results matched with the experimental data. On the basis of these experiments, we established a two-phase flow model of a polymer and oil, including the continuity, momentum, constitutive, and phase equations. The volume-of-fluid (VOF) method was used to track the interface between the two phases. A complex pore model was established based on the glass-etched model used in the experiment. We used the OpenFOAM platform to solve the mathematical model. The saturation, pressure, and stress tensor distributions were obtained. The results show that the displacement efficiency increased as the elasticity of the polymer increased; accordingly, the elasticity can enlarge the sweep area and decrease the residual oil saturation. As the elasticity increases, the stresses (the first normal stress, second normal stress, and shear stress) increase. Finally, the results obtained in this study can be used as a guideline in polymer design, screening, and optimization in the polymer flooding oilfields.

Graphical Abstract

1. Introduction

In a world of growing energy demands, although shale oil, coal, natural gas and other new energies continue to develop, oil remains the leading source of primary energy [1,2,3,4]. Global oil consumption increased from 0.8 million barrels per day in 2000 to approximately 100 million barrels of consumption per day in 2017. However, the increase in production has been slow, comprising only half that of consumption. In particular, the recent low oil prices have forced oil companies to enhance their oil recovery as a response to price challenges. Therefore, chemical flooding EOR techniques, including alkali/surfactant/polymer (ASP) flooding, foam flooding, and CO2 flooding, are being increasingly applied [5,6,7]. However, polymer flooding plays the most important role, particularly in the Chinese petroleum industry. In Daqing, the largest oilfield in China, polymer flooding has entered the commercial stages. The oil recovery can be increased 13% in comparison to that of water flooding [8]. The application of polymer flooding can guarantee continued supply in the future. Therefore, researchers have focused on the mechanisms of polymer flooding. In the early stage, a polymer can enhance oil recovery, increasing the viscosity and reducing water permeability, which can, accordingly, change the flow mobility ratio [9,10,11]. There are many such types of polymers, including synthetic polymers and biopolymers. To date, the Daqing oilfield has predominantly applied a synthetic polymer of HPAM. This material has specific advantages, such as the ability to tolerate high mechanical forces, a low cost, and an overall resistance to bacterial attack [10,12]. However, it is sensitive to brine salinity, making it ineffective for offshore oil due to the high salinity of sea water. The other group of polymers is the HAWP, which, in comparison to other polymers, have a smaller number of hydrophobic groups incorporated into the polymer backbone. In China, this polymer has been used in the BH offshore heavy oilfield. The high salinity and hardness of the water as well as the high viscosity of the oil require the polymer to be both salt-resistant and have good viscosity. Additionally, the polymer should have a strong anti-mechanical degradation ability as a result of the completion involving wire-wrapped and screen gravel-packed sand. Moreover, the long-term stability of the displacing fluid in the formation is important because of the large well spacing in offshore fields [13,14,15]. Therefore, HAWP can overcome the above shortcomings to satisfy these requirements. Its molecular structure is shown in Figure 1 [16]. HAWP is designed and synthesized as a small portion of a hydrophobic comonomer, which is to be incorporated into the backbone of a polymer. There are two general methods to incorporate hydrophobic moieties into water-soluble polymer chains: direct copolymerization of hydrophobic and water-soluble monomers and post-polymerization functionalization [17]. The former is the most common method of achieving associating polymer synthesis. Hydrophobically associating block copolymers are prepared by anionic polymerization [18,19,20]. Associating polymers involve acrylamide and free radical polymerization. Although other monomers have been used, arylamide is the most successful and most suited for producing high molecular weight water-soluble polymers [21,22]. This polymer is prepared by incorporating hydrophobic groups into the polymer after the polymerization process. It has potential for use in mobility control; its advantages include its constant viscosity when temperatures increase, attributable to its interchain association. It also can tolerate high salinities [23,24]. The viscosity can increase rapidly when the concentration is larger than the critical association concentration because of the network structure formed by the association of molecules, as shown in Figure 2 [25]. HAWP have been tested to be shear-stable at high flow rates. They showed a higher resistance and overall higher residual resistance factors than conventional polyacrylamides; moreover, they did not exhibit any injection problems. Therefore, they are applied in Chinese offshore reservoirs.
In many studies, researchers have predominantly concentrated on HAWP characteristics such as its molecular structure, rheology, field applications, and enhanced oil recovery [26,27,28,29,30,31]. However, fewer studies exist on the elasticity of HAWP flow in porous media. The elasticity of HPAM has attracted the attention of researchers for many years. When a polymer solution flows in a capillary tube with a constant diameter, the elasticity cannot be observed. However, in porous media, the shape and diameter of the pores and throat change rapidly, so elasticity can be obviously observed. Many investigators have studied the flow of polymer through porous media. The total pressure drop of polymer flooding is separated into the contributions of two parts, viscous contribution and elastic contribution parts; elastic parts can also enhance the oil recovery. In a 2000 study involving the HPAM and pure, viscous Xanthan solution with the same viscosity displacement oil on glass etched and core respectively, the results showed that the elasticity could enhance oil recovery and increase microsweep efficiency [32]. Xia et al. observed that, at higher elasticity, the residual oil could be mobilized because of the higher microforces that act on the different types of oil droplets in the pores, causing the oil droplet to change its shape and move. In the micro glass-etched experiment, they observed that the oil in the deadend was “pulled out” and “tripped” because of the elasticity effect [33]. Chen and Bogue proposed a model that can be used to calculate the residual oil saturation by a first normal stress difference as well as capillary number [34]. The first normal stress difference can present the elasticity of the polymer; results indicated a decrease in the residual oil saturation as the first normal stress difference increased. With the development of computers and computational methods, computational fluid dynamics (CFD) was used to study the elastic effect on the polymer displacement oil. In 2006, Yin et al. showed that the residual oil remained in different forms and that elasticity could make the residual oil move; concurrently, the study indicated that the microsweep area enlarged as elasticity increased [35]. In 2012, Yin et al. concluded from simulation results that viscoelasticity was the main property of a polymer in the increase of the microscopic sweep efficiency; the study also concluded that, with a higher Reynolds number, the microscale sweep efficiency was improved [36]. In 2013, Qi et al. simulated the elastic polymer flow in complex pores using the finite element method [37]. The pore model was derived from the high-resolution micro-CT images of a typical sandstone. The results showed that the fluid more deeply invaded the dead pores with higher elasticity; for higher Reynolds numbers, the fluid further penetrated the dead pores.
In this study, we investigate the HAWP rheology in porous media by comparing elastic HAWP and pure viscous Xanthan flow in cores; the aim of this experiment is to investigate the elasticity of the HAWP effect in porous media. On the basis of the experiment, the constitutive equation describing the rheology of HAWP is confirmed. In addition, CFD simulations were used to simulate the viscoelastic polymer displacement oil in the complex pore model. The results are significant for the optimum design of HAWP solutions in chemical flooding EOR processes.

2. Rheological Characterization of the Flow of HAWP Solutions in Porous Media

The important characteristics of a polymer solution flowing in porous media is that, when the flow velocity is low, it shows shear-thinning; shear-thickening appears when the flow velocity is high [38,39]. When single-phase flow is observed in the core, the relation of flow rate, shear rate, pressure gradient, and apparent viscosity can be expressed by the power-law Darcy Equation. The shear rate and apparent viscosity were calculated by Morais [40] in 2009 using the non-Newtonian law; they are expressed as:
γ ˙ = α Q A K ϕ
η ( γ ˙ ) = A K P Q
where γ ˙ is the apparent shear rate, 1/s; α is the non-Newtonian correction for a power-law fluid, in this paper, it is 4; Q is the flow rate, m3/s; A is the cross-sectional area of the core, m2; ϕ is the porosity; K is the permeability, m2; η ( γ ˙ ) is the apparent viscosity, Pa·s.

2.1. Experimental Solutions

In this study, APP-4, a type of HAWP synthesized by acrylamide and anionic copolymer monomer with a hydrophobic monomer unit, was employed. The molecular weight of the HAWP solutions used was 1.1 × 107 Da and varied with different concentration; the experimental water was prepared using double-distilled water, and its salinity is 9947.8 mg/L, which is the same as that of the BH oilfield actual formation water. The concentrations of the HAWP solutions are 1500 mg/L, 1750 mg/L, and 2000 mg/L. In this study, the other type of solution used was Xanthan, a pure viscous fluid, and a solution of concentration 3000 mg/L was prepared.

2.2. Core Material and Initialization

The core of this experiment was made using a quartz sand epoxy resin of diameter 25 mm and length 80 mm. The average permeability and porosity were 2000 × 10−3 µm2 and 25~30%, respectively. The characterizations were based on the properties of the BH oilfield reservoir.

2.3. Measurement Conditions

Measurements were made at 65 °C, the same temperature as the reservoir temperature. The flow rate of the pump was 0.003~10 mL/min and was increased from low to high; when the flow rate was constant, the pressure difference between the inlet and outlet was recorded. Figure 3 shows the diagram of apparent viscosity versus shear rate in the core.
In Figure 3, the hollow circle represents the results of the measurements made when HAWP solutions were flowing in the rheometer. We can observe that the viscosity decreased with increasing shear rate. This indicates the shear-thinning characteristic. However, when the single HAWP flowed in the porous media, the apparent viscosity decreased as the shear rate increase when the shear rate was low; correspondingly, it increased with shear rate increase when the shear rate was high. Therefore, the HWAP exhibited the shear-thickening phenomenon when they flowed in the core. The main reason for this occurrence is that shear effect only occurred when a fluid was flowing in the rheometer; however, when a fluid flowed in the core, elasticity was obvious because of the change in the shape and size of the pore and throat. Therefore, stretching can cause the polymer molecule to associate and the apparent viscosity to increase. Moreover, we observed that the Xanthan solution without elasticity did not exhibit the shear-thickening phenomenon in the core. This can validate that the reason mentioned above is correct. As shown in this Figure, we also determine the position where the rheology changes from shear-thinning to shear-thickening to maintain immobilization for different HAWP concentrations. It was observed that the position where shear-thickening occurred was not related to the concentration of the polymer; instead, it was potentially related to the molecular weight and structure.
Different rheologies can be shown for different polymers. To describe the rheology accurately in numerical simulations, we matched the apparent viscosity versus shear rate curve obtained for the numerical and experimental results. In Figure 4, the circle represents the point where the results were measured with a rheometer for 1750 mg/L HAWP solutions, and the different colored lines are the results of numerical simulation using different constitutive equations. The simulations were performed using a rectangular model, the ratio of which, length to width, was large enough and the flow was two-dimensional, incompressible, and isothermal. As shown in Figure 4, the deviation in the Oldroyd-B model results was minimum. Therefore, in the following numerical simulation, we used Oldroyd-B as the constitutive equation.

3. Mathematical Model

In this study, the mathematical model includes the continuity, momentum, and constitutive equations. Because this study included two phases, the polymer and oil, interface tracking was an important question; accordingly, an interface tracking equation was necessary. In this study, we adopted the VOF method.

3.1. Continuity Equation

This problem includes two continuity equations that can be written uniformly as one equation:
ρ t + ( ρ U ) = 0
where U is the velocity vector, m/s, which can be resolved into components of the velocity in the x and y directions. ρ is the density, when the point considered is in the polymer solution, it represents the density of the polymer solution. Otherwise, it represents the density of oil. In a two-phase transition zone, it is the mixture density; t is the flow time, s.

3.2. Momentum Equation

If it is a single-phase flow, the momentum equation is represented as:
( ρ U ) t + ( ρ U U ) = p + τ
In this equation, τ includes the viscoelastic polymer contribution τs, given by τs = 2ηsD, and the solvent contribution τp, which can be calculated by the constitutive equation. However, in two-phase flow problems, interfacial tension will induce a pressure jump. If the fluids are in equilibrium, fσ is the external body-force that is mechanically balanced by a pressure jump across the interface. Therefore, the momentum equation is rewritten as:
( ρ U ) t + ( ρ U U ) ( η s + η p ) ( U ) = p + τ p η p ( U ) + f σ
where ηs is the viscosity of the solvent and ηp is the viscosity of the polymer solutions. Then, fσ can be expressed as:
f σ = p = σ ( ( α | α | ) ) ( α )
where the fluid volume fraction α is defined as α = 1 for a point inside the polymer and α = 0 for a point inside the oil. Therefore, the density and dynamic viscosity in the equations of momentum are:
ρ = α ρ p + ( 1 α ) ρ o ,   μ = α μ p + ( 1 α ) μ o
where ρp is the polymer density and ρo is the oil density. They are substituted into the continuity equation, which can be reformulated in this manner:
α t + U α = 0
The equation system is closed with the constitutive equations of the two phases. The polymer rheology is described by the Oldroyd-B Equation, which is written as:
τ p + λ τ p = η 0 D
The power-law equation is used to describe the oil rheology:
η o = K γ ˙ n
where: η0 is the zero shear viscosity of the polymer solutions, Pa·s, and K is the consistency coefficient, Pa·sn+1. γ ˙ is the shear rate, s−1. These equations complete the mathematical description of the two-fluid system for laminar flow.

4. Physical Model

In this study, the physical model used was based on the glass etched in the experiments, shown in Figure 5a. In consideration of the memory and computation time of the computer, we picked one portion from the whole physical model: the red portion, as shown in Figure 5b–d, which represents the solid skeleton grid model and pore mesh subdivided by the snappyHexMesh solver in OpenFOAM. (OpenCFD Ltd., Bracknell, UK) There were 71,613 grids; the maximum radius of the pore was 1.24 mm and the minimum radius was 0.08 mm. This model can reflect the topology and complexity of the real porous media. Elasticity can be obviously observed in this complex model.

5. Numerical Simulation

In this study, Equations (3) and (5)–(8) constitute the entire governing system of equations. They are discrete and the coupling of velocity and pressure has been addressed. The solution procedure includes: ① initialization of all the fields (α, U, τ and p), ② solving the Equation (6) so that a new α can be calculated, along with the density and viscosity, ③ solve the momentum Equation (5) so that the velocity field can be obtained. Then, on the basis of continuity Equation (3), the correction velocity and the new velocity field can be calculated. Here, the PISO (Pressure Implicit Split Operator) algorithm was used; ④ solve the constitutive Equation, so that the stress tensor can be calculated. Return to step ② until the target time is reached.
In this study, we investigated four displacement patterns: water flooding, no elastic polymer flooding, and two different elastic polymer solution floodings. The parameters are shown in Table 1.

5.1. Volume Fraction Characteristics

In this study, the phase volume fraction can represent the saturation of water. Figure 6 shows the saturation for the four patterns at 11 s. As shown in Figure 6, in comparison with the case of no elastic polymer displacement, water movement was faster and the fingers grew faster in the water displacement oil. In comparison with the water flooding with no elastic polymer flooding, the elastic polymer flooding case reveals a steady front and better sweep area. As the relaxation time increased, the sweep area was larger.
Figure 7 shows the breakthrough times of the four different patterns. It indicated that the breakthrough time of water flooding was earlier than that of no elastic polymer flooding. A comparison between the breakthrough time of no elastic polymer flooding with that of viscoelastic polymer flooding indicated that the breakthrough time of elastic polymer flooding was later. Moreover, as the elasticity increased, the front breakthrough time was later. A comparison of Figure 7c,d also indicates that the sweep area enlarged with an increase in relaxation time.
Figure 8 shows the saturation distribution at t = 28 s after the displacement fluid breakthrough. Because of poor mobility ratio, the residual oil saturation of water flooding was the largest, and the residual oil saturation decreased as the elasticity of polymer increased, as shown in Figure 8c,d.

5.2. Pressure Distribution

Figure 9 shows the pressure distribution diagram. As shown in Figure 9, in comparison to the case of water flooding, when there was no elastic polymer flooding, an additional pressure drop was induced; moreover, the elastic polymer flooding could also produce an additional pressure drop compared to the case of viscous polymer flooding. As the elasticity increased, the additional pressure drop increased. Moreover, we observed that the pressure in elastic polymer flooding was distributed more uniformly in comparison to the cases of water flooding and no elastic polymer flooding. These two flooding techniques could both produce a local high pressure area. Figure 10 shows the pressure distribution when the displacement fluid front broke through. It is shown that the pressure variation is the same as that before the breakthrough.
On the basis of the pressure distribution, the pressure drop between inlet and outlet versus time is plotted in Figure 11. As can be seen, the pressure drop increased until the displacement fluid broke in the outlet.
As shown in Figure 11, early in the displacement process, as a result of the same inflow rate and outlet conditions, their pressure drops were also identical. Subsequently, as the elasticity increased, the pressure drop increased rapidly. We also observed that, before the front breakthrough, the pressure drop increased rapidly; in contrast, after the breakthrough, the pressure increase slowed. We can conclude that the elasticity of the polymer can produce an additional pressure drop that can increase the sweep efficiency and enhance oil recovery.

5.3. Stress Tensor Distribution

In this study, the problem is the two-dimensional flow, and the stress tensor included three components. Figure 12 is the first normal stress distribution for different relaxation times. As shown in this Figure, an obvious change in the first normal stress cannot be observed because of the low relaxation time. However, local enlargement (Figure to the right of the corresponding Figure) can reveal that the normal stress tensor increases as the relaxation time increases, especially in the thin pores and throat; additionally, the shape of the pore also changes sharply where the fluid becomes more elastic, suggesting elasticity can play an important role. This is the main reason that elasticity can contribute to enhanced oil recovery. This has been discussed in many reports, particularly those pertaining to the single-phase flow problem [41,42,43,44,45].
The shear stress simulation results are shown in Figure 13. We can easily see that, as the elasticity increased, relaxation time increased; this was a result of the direction of shear stress occurring in the flow direction. The maximum of the shear stress increased obviously, especially near the walls of the pore.
Figure 14 reveals the second normal stress distribution. As shown in the local enlargement Figure, the value was larger as the relaxation time increased. This phenomenon is in agreement with the results of the first normal stress mentioned above. They both validate that the elasticity can increase the displacement efficiency.

6. Conclusions

In this work, we focused on the mechanism of oil displacement using viscoelastic hydrophobically associating water-soluble polymers. The polymer flowing in porous media experiments reveals that the hydrophobically associating water-soluble polymer can exhibit shear thickening because of its elasticity; however, the viscous Xanthan solution only exhibited shear thinning. We observed that the position of occurrence was the same for different concentrations.
Furthermore, we established the mathematical model for solving the two-phase flow in porous media, which was developed for etched glass and can satisfy the geometry and topology of porous media.
We demonstrated that the polymer solution can have a stable front and better sweep area in comparison to water flooding. Elasticity is helpful in enlarging the sweep area. Before front breakthrough, the pressure drop between the inlet and outlet increased rapidly; subsequently, the increase was slow. As the elasticity increased, the additional pressure increased. As the elasticity increased, the first normal stress and shear stress both increased, especially in the small pore and throat.

Author Contributions

The idea of this work belongs to Huiying Zhong, she contributed to the study design, model development, and the overall analysis. Yuanyuan Li performed the experiments and established the mathematical models. Weidong Zhang and Daizong Guo contributed the numerical simulation. The detailed results discussion was performed by Hongjun Yin and Jun Lu. Study design of this manuscript was elaborated in discussion with the co-authors. Drafting of experimental material and methods parts, tables, and figures was also done in discussion with the co-authors. The initial drafts and the revised manuscript were prepared after discussion with the co-authors. The manuscript was finalized in cooperation with the co-authors.

Funding

This work presented in this paper was financially supported by the National Natural Science Funds for Young Scholars of China (Grant no. 51604079) and the Natural Science Foundation of Heilongjiang Province (Grant no. E2017012). The Natural Science Foundation of Heilongjiang Province (Grant no. E2016015) and the University Nursing Program for Young Scholars with Creative Talents in Heilongjiang Province (Grant no. UNPYSCT-2016126) are gratefully acknowledged.

Acknowledgments

The first author would like to show her gratitude to the University of Tulsa for the favorable support and invaluable inspirations to summarize this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cai, J.; Yu, B.; Zou, M.; Mei, M. Fractal analysis of invasion depth of extraneous fluids in porous media. Chem. Eng. Sci. 2010, 65, 5178–5186. [Google Scholar] [CrossRef]
  2. Liu, Y.; Li, J.; Wang, Z.; Wang, S.; Dong, Y. The role of surface and subsurface integration in the development of a high-pressure and low-production gas field. Environ. Earth Sci. 2015, 73, 5891–5904. [Google Scholar] [CrossRef]
  3. Wang, Z.; Le, X.; Feng, Y.; Zhang, C. The role of matching relationship between polymer injection parameters and reservoirs in enhanced oil recovery. J. Pet. Sci. Eng. 2013, 111, 139–143. [Google Scholar] [CrossRef]
  4. Rui, Z.; Wang, X.; Zhang, Z.; Lu, J.; Chen, G.; Zhou, X.; Patil, S. A realistic and integrated model for evaluating oil sands development with steam assisted gravity drainage technology in Canada. Appl. Energy 2018, 213, 76–91. [Google Scholar] [CrossRef]
  5. Zhong, H.; Zhang, W.; Fu, J.; Lu, J.; Yin, H. The performance of polymer flooding in heterogeneous type II reservoirs–An experimental and field investigation. Energies 2017, 10, 454. [Google Scholar] [CrossRef]
  6. Zhang, R.; Yin, X.; Winterfeld, P.H.; Wu, Y.S. A fully coupled thermal-hydrological-mechanical-chemical model for CO2 geological sequestration. J. Nat. Gas Sci. Eng. 2016, 28, 280–304. [Google Scholar] [CrossRef]
  7. Wu, X.; Pu, H.; Zhu, K.; Lu, S. Formation damage mechanisms and protection technology for Nanpu nearshore tight gas reservoir. J. Pet. Sci. Eng. 2017, 158, 509–515. [Google Scholar] [CrossRef]
  8. Wang, D.M.; Seright, R.S.; Shao, Z.B.; Wang, J.M. Key aspects of project design for polymer flooding at the Daqing oil field. SPE Reserv. Eval. Eng. 2008, 12, 1117–1124. [Google Scholar] [CrossRef]
  9. Cheng, J.; Wei, J.; Song, K. Study on remaining oil distribution after polymer injection. In Proceedings of the SPE Annual Technical Conference and Exhibition SPE 133808, Florence, Italy, 19–22 September 2010. [Google Scholar]
  10. Sheng, J.; Azri, N. Status of polymer-flooding technology. J. Can. Pet. Technol. 2005, 3, 116–120. [Google Scholar] [CrossRef]
  11. Wang, Z.H.; Liu, Y.; Pang, R.S.; Qiu, X.D. A performance evaluation of glycolipid biosurfactant synthesized by the enzyme-catalyzed method. Pet. Sci. Technol. 2005, 31, 2371–2377. [Google Scholar] [CrossRef]
  12. Guo, Y.; Peng, F.; Wang, H.; Huang, F.; Meng, F.; Hui, D.; Zhou, Z. Intercalation polymerization approach for preparing graphene/polymer composites. Polymers 2018, 10, 61. [Google Scholar] [CrossRef]
  13. Zhou, W.; Zhang, J.; Han, M.; Xiang, W.T.; Feng, G.Z.; Jiang, W.; Sun, F.J.; Zhou, S.W.; Ye, Z.B. Application of hydropholicallly associating water-soluble polymer for polymer flooding in China offshore heavy oilfield. In Proceedings of the International Petroleum Technology Conference IPTC 11635, Dubai, UAE, 4–6 December 2007. [Google Scholar]
  14. Liu, Y.; Wang, Z.H.; Zhuge, X.L.; Pang, R.S. Distribution of sulfide in oil-water treatment system and field test of treatment technology in Daqing oilfield. Pet. Sci. Technol. 2014, 32, 462–469. [Google Scholar] [CrossRef]
  15. Kang, X.D.; Zhang, J. Offshore heavy oil polymer flooding test in JZW area. In Proceedings of the SPE Heavy Oil Conference-Canada, SPE 165437, Calgary, AB, Canada, 11–13 June 2013. [Google Scholar]
  16. Chassenieux, C.; Nicolai, T.; Benyahia, L. Rheology of associative polymer solutions. Curr. Opin. Colloid Interface Sci. 2011, 16, 18–26. [Google Scholar] [CrossRef]
  17. Lu, H.; Feng, Y.; Huang, Z. Association and effective hydrodynamic thickness of hydrophobically associating polyacrylamide through porous media. J. Appl. Polym. Sci. 2008, 110, 1837–1843. [Google Scholar] [CrossRef]
  18. Valint, P.L., Jr.; Bock, J. Synthesis and characterization of hydrophobically associating block polymers. Macromolecules 1988, 21, 175–179. [Google Scholar] [CrossRef]
  19. Hill, A.; Candau, F.; Selb, J. Properties of hydrophobically associating polyacrylamides: Influence of the method of synthesis. Macromolecules 1993, 26, 4521–4532. [Google Scholar] [CrossRef]
  20. El-hoshoudy, A.N.; Desouky, S.E.M.; Elkady, M.Y.; Al-Sabagh, A.M.; Betiha, M.A.; Mahmoud, S. Hydrophobically associated polymers for wettability alteration and enhanced oil recovery-Article review. Egypt. J. Pet. 2017, 26, 757–762. [Google Scholar] [CrossRef]
  21. Samakande, A.; Hartmann, P.C.; Sanderson, R.D. Synthesis and characterization of new cationic quaternary ammonium polymerizable surfactants. J. Colloid Interface Sci. 2006, 296, 316–323. [Google Scholar] [CrossRef] [PubMed]
  22. Wever, D.A.Z.; Picchioni, F.; Broekhuis, A.A. Polymers for enhanced oil recovery: A paradigm for structure-property relationship in aqueous solution. Prog. Polym. Sci. 2011, 36, 1558–1628. [Google Scholar] [CrossRef]
  23. Alvarado, V.; Manrique, E. Enhanced oil recoverry: An update review. Energies 2010, 3, 1529–1575. [Google Scholar] [CrossRef]
  24. Wang, Z.H.; Liu, Y.; Le, X.P.; Yu, H.Y. The effects and control of viscosity loss of polymer solution compounded by produced water in oilfield development. Int. J. Oil Gas Coal Technol. 2014, 7, 298–307. [Google Scholar] [CrossRef]
  25. Carro, S.; Gonzalez-Coronel, V.J.; Castillo-Tejas, J.; Maldonado-Textle, H.; Tepale, N. Rheological properties in aqueous solution for hydrophobically modified polyacrylamides prepared in inverse emulsion polymerization. Int. J. Polym. Sci. 2017, 2017, 1–13. [Google Scholar] [CrossRef]
  26. Candau, F.; Biggs, S.; Hill, A.; Selb, J. Synthesis, structure and properties of hydrophobically associating polymers. Prog. Org. Coatings 1994, 24, 11–19. [Google Scholar] [CrossRef]
  27. Nystrom, B.; Walderhaug, H.; Hansen, F.K. Dynamic crossover effects observed in solutions of a hydrophobically associating water-soluble polymer. J. Phys. Chem. 1993, 97, 7743–7752. [Google Scholar] [CrossRef]
  28. Taylora, K.C.; Nasr-Ei-Dinb, H.A. Water-soluble hydrophobically associating polymers for improved oil recovery: A literature review. J. Pet. Sci. Eng. 1998, 19, 265–280. [Google Scholar] [CrossRef]
  29. Peiffer, D.G. Hydrophobically associating polymers and their interactions with rod-like micelles. Polymer 1990, 31, 2353–2360. [Google Scholar] [CrossRef]
  30. Lin, Y.; Luo, K.; Huang, R. A study on P(AM-DMDA) hydrophobically associating water-soluble copolymer. Eur. Polym. J. 2000, 36, 1711–1715. [Google Scholar] [CrossRef]
  31. Winkler, R.G.; Elgeti, J.; Gompper, G. Active polymers-emergent conformational and dynamical properties: A brief review. J. Phys. Soc. Jpn. 2017, 86, 101014. [Google Scholar] [CrossRef]
  32. Wang, D.; Cheng, J.; Yang, Q.; Gong, W.; Li, Q.; Chen, F.M. Viscous-elastic polymer can increase microscale displacement efficiency in cores. In Proceedings of the SPE Annual Technical Conference and Exhibition, SPE 63227, Dallas, TX, USA, 1–4 October 2000. [Google Scholar]
  33. Xia, H.; Wang, D.; Wu, J.; Kong, F. Elasticity of HPAM solutions increases displacement efficiency under mixed wettability conditions. In Proceedings of the SPE Asia Oil and Gas Conference and Exhibition, SPE 88465, Perth, Australia, 18–20 October 2004. [Google Scholar]
  34. Chen, J.; Bogue, D.C. Time-dependent stress in polymer melts and review of viscoelastic theory. J. Rheol. 1972, 16, 59. [Google Scholar] [CrossRef]
  35. Yin, H.; Wang, D.; Zhong, H. Study on flow behaviors of viscoelastic polymer solution in micropore with dead end. In Proceedings of the SPE Annual Technical Conference and Exhibition SPE 101950, San Antonio, TX, USA, 24–27 September 2006. [Google Scholar]
  36. Yin, H.; Wang, D.; Zhong, H.; Meng, S. Flow characteristics of viscoelastic polymer solution in micro-pores. In Proceedings of the SPE EOR Conference at Oil and Gas West Asia, SPE 154640, Muscat, Oman, 16–18 April 2012. [Google Scholar]
  37. Qi, M.; Wegner, J.; Falco, L.; Ganzer, L. Pore-scale simulation of viscoelastic polymer flow using a stabilised finite element method. In Proceedings of the SPE Reservoir Characterization and Simulation Conference and Exhibition SPE 165987, Abu Dhabi, UAE, 16–18 September 2013. [Google Scholar]
  38. Maerker, J.M.; Sinton, S.W. Rheology resulting from shear-induced structure in associating polymer solutions. J. Rheol. 1986, 30, 77. [Google Scholar] [CrossRef]
  39. Rui, Z.; Li, C.; Peng, P.; Ling, K.; Chen, G.; Zhou, X.; Chang, H. Development of industry performance metrics for offshore oil and gas project. J. Nat. Gas Sci. Eng. 2017, 39, 44–53. [Google Scholar] [CrossRef]
  40. Morais, A.F.; Seybold, H.; Herrmann, H.J.; Andrade, J.S. Non-newtonian fluid flow through three-dimensional disordered porous media. Phys. Rev. Lett. 2009, 103, 194502. [Google Scholar] [CrossRef] [PubMed]
  41. Clarke, A.; Howe, A.M.; Mitchell, J.; Staniland, J.; Hawkes, L.A. How viscoelastic-polymer flooding enhances displacement efficiency. SPE J. 2016, 21, 675–687. [Google Scholar] [CrossRef]
  42. Buchgraber, M.; Clemens, T.; Castanier, L.M.; Kovscek, A.R. A microvisual study of the displacement of viscous oil by polymer solutions. SPE Reserv. Eval. Eng. 2011, 14, 269–280. [Google Scholar] [CrossRef]
  43. Groisman, A.; Steinberg, V. Elastic turbulence in a polymer solution flow. Nature 2000, 405, 53–55. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Liu, Y. Exact solutions to nonlinear Schrödinger equation with variable coefficients. Appl. Math. Comput. 2011, 217, 5866–5869. [Google Scholar] [CrossRef]
  45. Lopez, X.; Valvatne, P.H.; Blunt, M.J. Predictive network modeling of single-phase non-Newtonian flow in porous media. J. Colloid Interface Sci. 2003, 264, 256–265. [Google Scholar] [CrossRef]
Figure 1. Molecular structure of HAWP: x = 30–100, y = 0–50, and z = 0.01–1.
Figure 1. Molecular structure of HAWP: x = 30–100, y = 0–50, and z = 0.01–1.
Polymers 10 00628 g001
Figure 2. Schematic diagram of HAWP association (C: concentration; C*: critical concentration): (a) C < C*, intramolecular association; (b) C < C*, intermolecular association; and (c) C > C*, the formation of a network structure by the association of molecules.
Figure 2. Schematic diagram of HAWP association (C: concentration; C*: critical concentration): (a) C < C*, intramolecular association; (b) C < C*, intermolecular association; and (c) C > C*, the formation of a network structure by the association of molecules.
Polymers 10 00628 g002
Figure 3. Apparent viscosity vs. shear rate for different polymer solutions.
Figure 3. Apparent viscosity vs. shear rate for different polymer solutions.
Polymers 10 00628 g003
Figure 4. Apparent viscosity matching with shear rate diagram.
Figure 4. Apparent viscosity matching with shear rate diagram.
Polymers 10 00628 g004
Figure 5. Mesh generation process for the physical simulation model: (a) the glass-etched model used in the microexperiment; (b) the black and white digitized model; (c) STL skeleton grid generated by solidworks; and (d) pore meshes generated by OpenFOAM.
Figure 5. Mesh generation process for the physical simulation model: (a) the glass-etched model used in the microexperiment; (b) the black and white digitized model; (c) STL skeleton grid generated by solidworks; and (d) pore meshes generated by OpenFOAM.
Polymers 10 00628 g005
Figure 6. Volume fraction distribution (t = 11 s): (a) water flooding; (b) polymer flooding (no elasticity); (c) polymer flooding (λ = 0.01); and (d) polymer flooding (λ = 0.03).
Figure 6. Volume fraction distribution (t = 11 s): (a) water flooding; (b) polymer flooding (no elasticity); (c) polymer flooding (λ = 0.01); and (d) polymer flooding (λ = 0.03).
Polymers 10 00628 g006
Figure 7. Displacement fluid breakthrough time: (a) 15 s water flooding; (b) 16 s no elastic polymer flooding; (c) 19 s viscoelastic polymer flooding (λ = 0.01); and (d) 21 s viscoelastic polymer flooding (λ = 0.03).
Figure 7. Displacement fluid breakthrough time: (a) 15 s water flooding; (b) 16 s no elastic polymer flooding; (c) 19 s viscoelastic polymer flooding (λ = 0.01); and (d) 21 s viscoelastic polymer flooding (λ = 0.03).
Polymers 10 00628 g007
Figure 8. Saturation distribution after breakthrough (t = 28 s): (a) water flooding; (b) no elastic polymer flooding; (c) viscoelastic polymer flooding (λ = 0.01); and (d) viscoelastic polymer flooding (λ = 0.03).
Figure 8. Saturation distribution after breakthrough (t = 28 s): (a) water flooding; (b) no elastic polymer flooding; (c) viscoelastic polymer flooding (λ = 0.01); and (d) viscoelastic polymer flooding (λ = 0.03).
Polymers 10 00628 g008aPolymers 10 00628 g008b
Figure 9. Volume fraction distribution (t = 11 s): (a) water flooding; (b) no elastic polymer flooding; (c) viscoelastic polymer flooding (λ = 0.01); and (d) viscoelastic polymer flooding (λ = 0.03).
Figure 9. Volume fraction distribution (t = 11 s): (a) water flooding; (b) no elastic polymer flooding; (c) viscoelastic polymer flooding (λ = 0.01); and (d) viscoelastic polymer flooding (λ = 0.03).
Polymers 10 00628 g009
Figure 10. Displacement fluid front breakthrough times: (a) 15 s; (b) 16 s; (c) 19 s; and (d) 21 s.
Figure 10. Displacement fluid front breakthrough times: (a) 15 s; (b) 16 s; (c) 19 s; and (d) 21 s.
Polymers 10 00628 g010
Figure 11. Displacement pressure drop vs. time for different elastic polymer solutions.
Figure 11. Displacement pressure drop vs. time for different elastic polymer solutions.
Polymers 10 00628 g011
Figure 12. First normal stress (τxx) distribution: (a) λ = 0.01 s; (b) λ = 0.05 s.
Figure 12. First normal stress (τxx) distribution: (a) λ = 0.01 s; (b) λ = 0.05 s.
Polymers 10 00628 g012
Figure 13. Shear stress (τxy) distribution: (a) λ = 0.01 s; (b) λ = 0.05 s.
Figure 13. Shear stress (τxy) distribution: (a) λ = 0.01 s; (b) λ = 0.05 s.
Polymers 10 00628 g013
Figure 14. Shear stress (τxy) distribution: (a) λ = 0.01 s; (b) λ = 0.05 s.
Figure 14. Shear stress (τxy) distribution: (a) λ = 0.01 s; (b) λ = 0.05 s.
Polymers 10 00628 g014
Table 1. Parameters of the four computational cases considered in this study.
Table 1. Parameters of the four computational cases considered in this study.
CaseFlow Rate [m3/s]OilDisplacement FluidInterfacial Tension [mN/m]
Density [kg/m3]Viscosity [mPa·s]Density [kg/m3]Viscosity [mPa·s]Relaxation Time [s]
Waterflooding1 × 10−10860910001-4.8
No elastic polymer flooding1 × 10−10860990012-4.8
Elastic polymer flooding (λ = 0.01)1 × 10−108609900ηsηp0.014.8
111
Elastic polymer flooding (λ = 0.03)1 × 10−1086099001110.034.8

Share and Cite

MDPI and ACS Style

Zhong, H.; Li, Y.; Zhang, W.; Yin, H.; Lu, J.; Guo, D. Microflow Mechanism of Oil Displacement by Viscoelastic Hydrophobically Associating Water-Soluble Polymers in Enhanced Oil Recovery. Polymers 2018, 10, 628. https://doi.org/10.3390/polym10060628

AMA Style

Zhong H, Li Y, Zhang W, Yin H, Lu J, Guo D. Microflow Mechanism of Oil Displacement by Viscoelastic Hydrophobically Associating Water-Soluble Polymers in Enhanced Oil Recovery. Polymers. 2018; 10(6):628. https://doi.org/10.3390/polym10060628

Chicago/Turabian Style

Zhong, Huiying, Yuanyuan Li, Weidong Zhang, Hongjun Yin, Jun Lu, and Daizong Guo. 2018. "Microflow Mechanism of Oil Displacement by Viscoelastic Hydrophobically Associating Water-Soluble Polymers in Enhanced Oil Recovery" Polymers 10, no. 6: 628. https://doi.org/10.3390/polym10060628

APA Style

Zhong, H., Li, Y., Zhang, W., Yin, H., Lu, J., & Guo, D. (2018). Microflow Mechanism of Oil Displacement by Viscoelastic Hydrophobically Associating Water-Soluble Polymers in Enhanced Oil Recovery. Polymers, 10(6), 628. https://doi.org/10.3390/polym10060628

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop