The mathematical models of the flow of polymer foam in porous media under three injection modes are established and the relevant numerical calculation methods are given. The profiles of the liquid phase saturation, the pressure drop and the number density of the flowing HPAM foam in artificial sandstone cores with the dimensionless distance under three injection modes are numerically calculated and analyzed. The results show that, compared with the injection mode 2 and 3, HPAM foam flows in a piston-like fashion in the artificial sandstone core under the injection mode 1 and produces the biggest pressure drop. Obviously, the flood efficiency is the highest under the injection mode 1.
Injection Mode Polymer Foam Porous Media Mathematical Model Numerical Calculation1. Introduction
It can provides the necessary theoretical foundation for polymer flooding engineering design through researching the flow of the polymer foam in porous media under different injection modes. Shan, et al. [1] given an optimal foam-injection strategy in homogeneous reservoirs with a variety of foam models by using simulation. Stevens, et al. [2] presented an overview of the operating plan for CO2 foam injection and the details of the CO2 foam injection schedule and design criteria in the East Vacuum Grayburg San Andres Unit. Hou, et al. [3] compared three main injection modes, including co-injection of gas and solution, alternative injection of gas and solution, and direct injection of foam and pointed that the direct injection of foam was the most effective mode and the alternative injection of gas and solution is the poorest. Li, et al. [4] experimentally studied the use of additives to enhance foam properties and the factors affecting blocking performance such as injection model, gas liquid ratio, and alternate slug, etc. Harsenhorst, et al. [5] illustrated with a model fit by Rossen, et al. [6] to data of Persoff, et al. [7] , a fit specifically designed for a hypothetical surfactant-alternating-gas (SAG) application. The mathematical models of the flow of polymer foam in porous media under three injection modes will be established and the numerical calculation methods will be given in this paper. According to the mathematical models and the numerical calculation methods metioned above, the profiles of the liquid phase saturation, the pressure drop and the number density of the flowing HPAM foam in artificial sandstone core under three injection modes metioned above will be numerically calculated and analyzed.
2. Injection Modes
There are three injection modes for the porous media filled with polymer-surfactant aqueous solution,that is, the injection mode 1 which is the co-injection of polymer surfactant aqueous solution and air into the porous media, the injection mode 2 which is the injection of air into the porous media and the injection mode 3 which is the injection of polymer surfactant aqueous solution into the porous media.
3. Mathematical Models
According to the mass and the population balance principles, the mathematical models of polymer foam flow in porous media under different injection modes can be obtained.
3.1. Injection Mode 1
The mathematical model of the flow of polymer foam in porous media under the injection mode 1 is as follows [8] -[11] .
The mass balance equations
the population balance equations
the assistant equations
the initial conditions
the boundary conditions
where the porosity; the liquid phase mass density,; the saturation of liquid phase; the permeability of porous media ,; the liquid phase relative permeability; the apparent viscosity of liquid phase,; the pressure of liquid phase, Pa; the source-sink term of liquid phase,; the gas mass density,; the saturation of gas phase; the gas phase relative permeability; the gas viscosity,; the pressure of gas phase,; the source-sink term of gas phase,; the gas phase saturation of flowing foam; the number density of flowing foam,; the gas phase saturation of trapped foam; the number density of trapped foam,; the velocity of the flowing foam in porous media,; the foam generation coefficient,; the local interstitial velocity of flowing foam,; the local interstitial liquid velocity,; the foam coalescence coefficient,; the liquid phase saturation relevant to the limiting capillary pressure; the source-sink term for foam bubbles in units of number per unit volume per unit time,; the liquid relative permeability where; the connate liquid phase saturation; the index; the gas relative permeability where; the fraction of the gas phase flow rate; the index; the air density under standard atmosphere,; the pressure under standard atmosphere，; flow behavior index of liquid phase; the consistency coefficient of liquid phase,; the pore tortuosity; the initial gas phase viscosity,; the constant of proportionality dependent primarily on the surfactant system,; the index; [12] [13] the foam effective viscosity,; the capillary pressure, Pa; the equilibrium surface tension,; Leverett J function;, are the bubble trapped parameters; is the initial liquid phase pressure,; the injection velocity of liquid phase,; the injection velocity of gas phase,; the length of porous media,; the backpressure,.
3.2. Injection Mode 2
Under the injection mode 2, where the Formula (19) of the mathematical model of the flow of the polymer foam in porous media under the injection mode 1 changes as follow
the mathematical model of the flow of the polymer foam in porous media consists of Formula (1)-(18), (25), (20)-(24).
3.3. Injection Mode 3
Under the injection mode 3, the mathematical model of the flow of the polymer foam in porous media consists of Formula (1), (13), (15), (19) and (21).
4. Numerical Calculation Methods
The mathematical models of the flow of polymer foam in porous media under three injection modes are numerically calculated by the finite difference method [9] [14] [15] .
For the Injection mode 1 and the injection mode 2, the mass balance Equations (1) and (2) are solved by IMPES (Implicit Pressure Explicit Saturation) method, and the population balance Equation (3) is solved by TVD (Total Variation Diminish) three order method. For the Injection mode 3, the mass balance Equation (1) is solved by IM (Implicit Method).
5. Numerical Calculations and Analyses
The flow of the foam of the 0.015 wt% HPAM, 0.200 wt% FL-605 and 0.830 wt% NaCl HPAM solution and the air under different injection modes in the artificial sandstone core are numerically calculated and analyzed.
5.1. Basic Data
The basic data of the three injection modes are shows as Table 1 [8] [9] .
5.2. Injection Mode 1
Under the injection mode 1, the profiles of the liquid phase saturation, the pressure drop, the number density of the flowing HPAM foam in the artificial sandstone core show as Figures 1-3.
In Figure 1, under the injection mode 1, the liquid phase saturation of HPAM foam at the flow front appears to be steep and sharp, the liquid phase saturation of the upstream of the flow front is about 0.38, and 1.00 of the downstream. The gas phase breaks through at about dimensionless time 0.80 PV. The liquid phase saturation at steady state is about 0.38. HPAM foam flows in a piston-like fashion in the artificial sandstone core. After HPAM foam flow front passes a particular location of the core, the liquid phase saturation changes very little.
In Figure 2, under the injection mode 1, the pressure drop of the upstream of HPAM foam flow front decreases along with the dimensionless distance x/L increasing, the pressure drop of the downstream of HPAM foam flow front is zero. The pressure drop doesn’t increase up to 948.99 kPa along with dimensionless time PV elapsing until about 2 PV.
In Figure 3, under the injection mode 1, there is the flow front of HPAM foam in the artificial sandstone core. The number density of flowing foam reaches the peak value at the flow front. The steady state profile of the
Basic data
Injection mode 1
Injection mode 2
Injection mode 3
1000.09
1000.09
1000.09
n
0.9316
0.9316
0.9316
0.1359
0.1359
0.1359
11.03
11.03
0.00
0.30
0.30
0.30
0.66
0.66
0.66
0.26
0.26
0.26
2.30
2.30
2.30
0.70
0.70
0.00
1.00
1.00
0.00
3.01
3.01
0.00
3.03
3.03
0.00
0.29
0.29
0.29
0.26
0.26
0.00
0.33
0.33
0.00
0.18
0.18
0.00
0.01
0.01
0.01
6.00
6.00
0.00
0.33
0.33
0.00
4.80
4.80
4.80
4.80
4.80
4.80
0.10
0.10
0.00
0.046
0.000
0.046
0.431
0.431
0.000
10.08
10.08
0.00
57.84
57.84
0.00
Profiles of the liquid phase saturation <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-6201857x147.png" xlink:type="simple"/></inline-formula> under the injection mode 1
Profiles of the pressure drop <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-6201857x149.png" xlink:type="simple"/></inline-formula> under the injection mode 1
Profiles of the number density <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-6201857x151.png" xlink:type="simple"/></inline-formula> of flowing foam under the injection mode 1
number density of flowing foam appears at 2 PV, after which the number density of flowing foam decreases from the inlet to the outlet.
5.3. Injection Mode 2
Under the injection mode 2, the profiles of the liquid phase saturation, the pressure drop, the number density of the flowing HPAM foam in the artificial sandstone core show as Figures 4-6.
In Figure 4, under the injection mode 2, the liquid phase saturation of the upstream of HPAM foam flow front is about 0.32, the liquid phase saturation at the flow front appears to be steep and sharp and the liquid phase saturation of the downstream of the flow front is 1.00, which is because the gas phase flow and the foam flow appear in the upstream of the flow front, the foam flow and the liquid phase flow appear in the upstream of the flow front under the injection mode 2.
In Figure 5, under the injection mode 2, the pressure drop of the upstream of HPAM foam flow front is about 27 kPa, the pressure drop at the flow front changes sharply and the pressure drop of the
Profiles of the liquid phase saturation <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-6201857x164.png" xlink:type="simple"/></inline-formula> under the injection mode 2
Profiles of the pressure drop <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-6201857x166.png" xlink:type="simple"/></inline-formula> under the injection mode 2
Profiles of the number density <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-6201857x168.png" xlink:type="simple"/></inline-formula> of flowing foam under injection mode 2
downstream of flow front decreases slowly, which is due to the fine foam texture of flow front under injection mode 2.
In Figure 6, under the injection mode 2, the number density of the flowing HPAM foam reaches the peak value sharply only at the flow front, and then the number density of the flowing HPAM foam gets to be zero, which is because no sufficient liquid phase supplement leads the foam coalescence velocity at the flow front accelerating and the number density of the flowing HPAM foam decreasing sharply under the injection mode 2.
5.4. Injection Mode 3
Under the injection mode 3, the profiles of the liquid phase saturation, the pressure drop, the number density of the flowing HPAM aqueous solution in the artificial sandstone core show as Figures 7-9.
In Figure 7, under the injection mode 3, the flow of the HPAM aqueous solution in the artificial sandstone core is a single phase fluid flow. Obviously, the liquid phase saturation is 1 constantly.
In Figure 8, under the injection mode 3, the pressure drop of the flowing HPAM aqueous solution changes linearly.
Profiles of the liquid phase saturation <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-6201857x178.png" xlink:type="simple"/></inline-formula> under the injection mode 3
Profiles of the pressure drop <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-6201857x180.png" xlink:type="simple"/></inline-formula> under the injection mode 3
Profiles of the number density <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/4-6201857x182.png" xlink:type="simple"/></inline-formula> of flowing foam under the injection mode 3
In Figure 9, under the injection mode 3, the number density of the flowing HPAM aqueous solution is zero obviously.
5.5. Comparisons of Injection Modes
Among the three injection modes mentioned above, the pressure drop of the flowing HPAM foam under the injection mode 1 is the biggest, which is because there is the highest number density of the flowing HPAM foam under the injection mode 1. Obviously, the flood efficiency of the injection 1 is the highest.
6. Conclusions
1) The mathematical models of the flow of the polymer foam in porous media under the three injection modes are established and the relevant numerical calculation methods are given.
2) The profiles of the liquid phase saturation, the pressure drop and the number density of the flowing HPAM foam with the dimensionless distance in the artificial sandstone core under the three injection modes are numerically calculated and analyzed.
3) Compared with the injection mode 2 and 3, HPAM foam flows in a piston-like fashion in the artificial sandstone core under the injection mode 1 and shows the biggest pressure drop. Obviously, flood efficiency is highest under the injection mode 1.
4) Based on the research results mentioned above, it can choose the better injection mode to reach the best flood efficiency in the polymer flooding engineering design according to the actual situation. Especially, in the injection mode 1, it also needs the future study to continue working and researching for what factor influencing the flood efficiency most, such as the polymer type, concentration and so on.
Cite this paper
WeiZhao,HaiqingCui,KeliangWang, (2015) Polymer Foam Flow in Porous Media under Different Injection Modes. Energy and Power Engineering,07,365-374. doi: 10.4236/epe.2015.78034
ReferencesShan, D. and Rossen, W.R. (2002) Optimal Injection Strategies for Foam IOR. SPE Journal, 9, 1-17.Steven, J.E., Harpole, K.J., Zornes, D.R. and Martin, F.D. (1992) CO2 Foam Field Verification Pilot Test at EVGSAU: Phase II-Foam Injection Design and Operation Plan. SPE Annual Technical Conference and Exhibition, Washington, DC, 4-7 October 1992, 115-128.Hou,Q.F., Zhu, Y.Y., Luo, Y.S. and Weng, R. (2012) Studies on Foam Flooding EOR Technique for Daqing Reservoirs after Polymer Flooding. Proceedings of the SPE Improved Oil Recovery Symposium, Tulsa, 14-18 April 2012, 1-11. http//dx.doi.org/10.2118/151955-msLi, Z.Q., Zhou, G.H. and Zhou, Z.M. (2006) The Feasibility Studies of Polymer Foam Flooding in Gudao Oilfield. Proceedings of the SPE Asia Pacific Oil & Gas Conference and Exhibition, Adelaide, 11-13 September, 1-8.
http//dx.doi.org/10.2118/101189-msHarsenhorst, R.M., Dharma, A.S., Andrianov, A. and Rossen, W.R. (2013) Extension of a Simple Model for Vertical Sweep in Foam SAG Displacement. Proceedings of the EAGE Annual Conference & Exhibition Incorporating SPE Europec, London,10-13 June 2013, 1-11. http//dx.doi.org/10.2118/164891-msBoeije, C.S. and Rossen, W.R. (2013) Fitting Foam Simulation Model Parameters to Data. Proceedings of the IOR 2013—17th European Symposium on Improved Oil Recovery, St. Petersburg, 16-18 April 2013, 1-16.
http//dx.doi.org/10.3997/2214-4609.20142604Persoff, P., Radke, C.J., Pruess, K, Benson, S.M. and Witherspoon, P.A. (1991) A Laboratory Investigation of Foam Flow in Sandstone at Elevated Pressure. SPE Reservoir Engineering, 6, 365-372. http//dx.doi.org/10.2118/18781-PAKovseck, A.R. (1994) Foam displacement in Porous Media: Experiment and Mechanistic Prediction by the Population Method. University of California at Berkely, Berkely.Zhao, W. (2015) Polymer Foam Transient Flow in the Porous Media. Northeast Petroleum University, Daqing.Bertin, H.J., Quintard, M.Y. and Castainier, L.M. (1998) Development of a Bubble-Population Correlation for Foam-Flow Modeling in Porous Media. SPE Journal, 3, 356-362. http//dx.doi.org/10.2118/52596-paGarrouch, A.A. (1999) A Viscoelastic Model for Polymer Flow in Reservoir Rocks. Proceedings of the SPE Asia Pacific Oil and Gas Conference and Exhibition, Jakarta, 20-22 April 1999, 1-10.
http//dx.doi.org/10.2118/54379-msFalls, A.H., Musters, J.J. and Ratulowski, J. (1989) The Apparent Viscosity of Foams in Homogeneous Bead Packs. SPE Reservoir Engineering, 4, 155-164.Hirasaki, G.J. and Lawson, J.B. (1985) Mechanisms of Foam Flow in Porous Media: Apparent Viscosity in Smooth Capillaries. Society of Petroleum Engineers Journal, 25, 176-190. http//dx.doi.org/10.2118/12129-paAziz, K. and Settari, A. (1979) Petroleum Reservoir Simulation. Applied Science Publishers, London, 135-139.Leonard, B.P. (1984) Third Order Upwinding as a Rational Basis for Computational Fluid Dynamics. Proceedings of the Computational Techniques and Applications: CTAC-83, Sydney, 1983, 106-120.