Numerical Calculation of Transient Flow of Polymer Foam in Porous Media ()
1. Introduction
The researches about the flow of aqueous foam in porous media are published more in literatures and reports. Kovseck [1] established the mathematical model of the transient flow of the aqueous foam in porous media based on the population balance equation and gave the relevant numerical calculation method; Bertin [2] [3] presented the mathematical model of transient flow of the aqueous foam in porous media based on the bubble correlation equation and compared the numerical calculated results of transient flow of the aqueous foam based on the bubble correlation equation in porous media with those based on the population balance equation. And the relevant research about the transient flow of polymer foam in porous media is limited in the experiment research in most cases now, due to its non-Newtonian property of liquid phase causing the flow of polymer foam more complicated than that of the aqueous foam in porous media. Sydansk [4] derived polymer-enhanced foams through the laboratory experiments, which were shown to exhibit, as compared with their counterpart conventional foams, the substantially higher effective viscosities and the much improved stability; Safouane, et al. [5] studied the drainage of foams made from Newtonian and non-Newtonian solutions of different viscosities and showed that foams made of solutions containing long flexible Poly Ethylene Oxide molecules counter-intuitively drained faster than foams made of Newtonian solutions of the same viscosity.
By using the finite difference method, one dimension transient flow of the polymer foam in porous media will be numerically calculated. In this paper, the profiles of the liquid phase saturation, the pressure drop and the number density of one dimension transient flow of HPAM foam with the dimensionless time in artificial sandstone core will be analyzed.
2. Mathematical Model
2.1. Hypothesis Conditions
1) The porous media is incompressible;
2) Polymer foam liquid phase is considered as the incompressible power-law fluid. The concentration of surfactant is constant. Gas phase obeys the ideal gas law;
3) Extensive Darcy law is adaptable, foam only influence on the gas phase relative permeability and viscocity;
4) Bubbles of the flow of the polymer foam in porous media are in the same size;
5) The effects of convection, diffusion, adsorption and gravitation of the polymer foam liquid phase in porous media are neglected;
6) The foam relative permeability is characterized by the gas relative permeability model.
2.2. Mass Balance Equation
The liquid and gas phase mass balance equations of one dimension transient flow of the polymer foam in porous media are [1] [6] [7]
(1)
(2)
where
(3)
(4)
(5)
(6)
(7)
where
is 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 local interstitial velocity of flowing foam,
;
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 number density of flowing foam,
;
the local interstitial velocity of flowing foam,
;
the index;
the foam effective viscosity,
.
2.3. Population Balance Equation
The population balance equation of one dimension transient flow of the polymer foam in porous media is [1]
(8)
(9)
(10)
(11)
(12)
(13)
where
is the gas phase saturation 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 liquid velocity,
;
the foam coalescence coefficient,
;
the liquid phase saturation relavant to the limiting capillary pressure;
the source-sink term for foam bubbles in units of number per unit volume per unit time,
;
the fraction of the gas phase trapped rate;
the relative permeability of the foam that is flowing;
the maximum fraction of the gas phase that trapped,
trapped parameter.
2.4. Assistant Equation
2.4.1. Relationship between Pressure
(14)
(15)
where [3]
is the capillary pressure, Pa;
the equilibrium surface tension,
;
the Leverett J fuction.
2.4.2. Relationship between Saturation
(16)
2.4.3. Relationship between the Number Density of Flowing Foam and the Number Density of Trapped Foam
(17)
where
,
are the bubble trapped parameters [1] .
2.5. Definite Condition
2.5.1. Initial Condition
(18)
(19)
(20)
(21)
(22)
(23)
where
is the initial liquid phase pressure,
.
2.5.2. Boundary Condition
(24)
(25)
(26)
(27)
(28)
(29)
where
is the injection velocity of liquid phase,
;
the injection velocity of gas phase,
;
the length of porous media,
,
the backpressure,
.
The mass balance Equations (1) and (2), the bubble population balance Equation (8), the assistant equations (14), (16), (17), the definite condition (18)-(29) mentioned above consist of the mathematical model of one dimension transient flow of the polymer foam in porous media.
It is generally considered that the polymer foam generation coefficient
and the coalescence coefficient
of the transient flow of the polymer foam in porous media are influenced by the physical properties of the polymer foam liquid phase, the polymer foam gas phase and the porous media, which needed to be determined by experiments.
3. Numerical Calculation Method
The mathematical model mentioned above is numerically calculated by using the finite difference method.
The mass balance Equations (1) and (2) of the mathematical model mentioned above are solved by using IMPES (Implicit Pressure Explicit Saturation) method. And the bubble population balance Equation (8) is solved by using TVD (Total Variation Diminish) three order methods [6] [8] [9] .
4. Experiment Research
Through the laboratory experiments, with measuring the time-average pressure drop
of one dimension transient flow of the HPAM foam in the artificial sandstone core, the generation coefficient
and the coalescence coefficient
of the mathematical model mentioned above in the artificial sandstone core are determined.
4.1. Experiment Meter
The main experiment meters include: DV-II + Pro viscometer, produced by American Brookfield Company, the measurement error is ±1%; TX
-500C
Spining interfacial tension meter, produced by American Bowing Company, the measurement error is 0.001 mm.
4.2. Experiment Material
HPAM: foam liquid phase is the aqueous solution of 0.045 wt% HPAM, 0.200 wt% FL-605 and 0.830 wt% NaCl, gas phase is air. The core is the homemade artificial sandstone core of 0.045 m × 0.045 m × 0.30 m.
4.3. Experiment Programme
The experiment programme shows in Figure 1.
“1” D
-250L
liquid pump; “2” D
-250L
gas pump; “3” PMQ-2 foam performance measurement apparatus; “4” the liquid can; “5” the gas can; “6” digital manometer; “7” back-pressure valve; “8” the homemade artificial sandstone; “9” digital manometer; “10” cylinder.
Shut the back-pressure valve “
7”
of the outlet of the artificial sandstone core “8” of PMQ-2 foam performance measurement apparatus “3” and set the backpressure; Start D
-250L
liquid pump “
1”
and the liquid can “4” and set the liquid injection velocity
; start D
-250L
gas pump “2” and the gas can “5” and set the air injection velocity
; measure the time-average pressure drop
of the homemade artificial sandstone core “8” continuouslly by using the digital manometer “9” and “6” of “
3”
.
4.4. Experimental Result
4.4.1. Experimental Parameters
The experimental parameters of one dimension transient flow of the HPAM foam of 0.045 wt% HPAM, 0.200 wt% FL-605 and 0.083 wt% NaCl in the artificial sandstone core show as Table 1 [1] .
4.4.2. Determination of the Generation and Coalescence Coefficient
By using the mathematical model, the relevant numerical calculated method, Table 1 mentioned above and B (Bisection) iterated algorithm, the generation coefficient
and the coalescence coefficient
of one dimension transient flow of HPAM foam of 0.045 wt% HPAM, 0.200 wt% FL-605 surfactant and 0.830 wt% NaCl in the artificial sandstone core can be determined as
,
[1] [6] .
5. Numerical Calculation Results
By using the mathematical model mentioned above, the relevant numerical calculated method, Table 1 and the generation coefficient
and the coalescence coefficient
, the profiles of the liquid phase saturation
, the pressure drop
, the number density
of one dimension transient flow of HPAM foam of 0.045 wt% HPAM, 0.200 wt% FL-605 surfactant and 0.830 wt% NaCl in the artificial sandstone core with the di- mensionless time PV are numerical calculated and analyzed.
5.1. Profiles of Liquid Phase Saturation
The profiles of the liquid phase saturation
of the one dimension transient flow of HPAM foam of 0.045 wt% HPAM, 0.200 wt% FL-605 surfactant and 0.830 wt% NaCl in the artificial sandstone core with the di- mensionless time PV show as Figure 2.
![]()
Figure 2. Profiles of liquid phase saturation
of HPAM foam with the dimensionless time PV.
In Figure 2, the liquid phase saturation
of HPAM foam of 0.045 wt% HPAM, 0.200 wt% FL-605 surfactant and 0.830 wt% NaCl at the flow front appears to be steep and sharp, the liquid phase saturation
of the upstream of the flow front is about 0.42, and 1.00 of the downstream. The gas phase breaks through at about dimensionless time 0.70 PV. The liquid phase saturation
at steady state is about 0.42. HPAM foam flows in a piston-like fashion in the artificial sandstone core. After HPAM foam flow front passing a particular location of the core, the liquid phase saturation
changes very little.
5.2. Profiles of Pressure Drop
The profiles of the pressure drop
of the one dimension transient flow of HPAM foam of 0.045 wt% HPAM, 0.200 wt% FL-605 surfactant and 0.830 wt% NaCl in the artificial sandstone core with the di- mensionless time PV show as Figure 3.
![]()
Figure 3. Profiles of pressure drop
of HPAM foam with the dimensionless time PV.
In Figure 3, the pressure drop
of the flow of HPAM foam of 0.045 wt% HPAM, 0.200 wt% FL-605 surfactant and 0.830 wt% NaCl decreases away from the core entrance with the dimensionless distance x/L. The pressure drop
of the downstream of the flow front is close to 0.00 kPa. The pressure drop
doesn’t increase up to 1130.21 kPa along with dimensionless time PV elapsing until about 2 PV. The steady pressure drop
is 1130.21 kPa.
5.3. Profiles of Number Density of Flowing HPAM Foam
The profiles of the number density
of the one dimension transient flow of HPAM foam of 0.045 wt% HPAM, 0.200 wt% FL-605 surfactant and 0.830 wt% NaCl in the artificial sandstone core with the di- mensionless time PV show as Figure 4.
![]()
Figure 4. Profiles of the number density
of HPAM foam with the dimensionless time PV.
In Figure 4, there is the flow front of the flow of HPAM foam of 0.045 wt% HPAM, 0.200 wt% FL-605 surfactant and 0.830 wt% NaCl in the artificial sandstone core. The number density
of flowing HPAM foam reaches the peak value at the flow front. The steady state profile of the number density
of flowing HPAM foam appears at 2 PV, after which the number density
of flowing HPAM foam decreases from the inlet to the outlet.
6. Experiment Verification
By using the mathematical model, the relevant numerical calculated method, Table 1 and the generation coefficient
, the coalescence coefficient
mentioned above, the curves of the calculated pressure drop
and the experiment pressure drop
of one dimension transient flow of HPAM foam of 0.045 wt% HPAM, 0.200 wt% FL-605 and 0.830 wt% NaCl with dimensionless time PV in the artificial sandstone core are compared as Figure 5.
![]()
Figure 5. Comparison between the curves of the numerical calculated pressure drop
and the experiment pressure drop
of HPAM foam.
In Figure 5, the curves of the calculated pressure drop
curve and the experiment pressure drop
of one dimension transient flow of HPAM foam of 0.045 wt% HPAM, 0.200 wt% FL-605 and 0.830 wt% NaCl in artificial sandstone core are basically coincided, which show that the mathematical model of one dimension transient flow of polymer foam in porous media and the relevant numerical calculation method in the paper are correct [10] .
7. Conclusions
1) Based on the mathematical model of one dimension transient flow of the polymer foam in porous media, the numerical calculation method of the flow mentioned above by using the finite difference method is given;
2) Through the experiments, the generation coefficient and the coalescence coefficient of the transient flow of HPAM foam of 0.045 wt% HPAM, 0.200 wt% FL-605 surfactant and 0.830 wt% NaCl in the artificial sandstone core are determined;
3) The profiles of the liquid phase saturation ,the pressure drop and the number density of one dimension transient flow of HPAM foam of 0.045 wt% HPAM, 0.200 wt% FL-605 and 0.830 wt% NaCl with the dimen- sionless time in the artificial sandstone core are numerically calculated and analyzed.
4) Through the experiments, the mathematical model of one dimension transient flow of polymer foam in porous media and the relevant numerical calculation method are verified.
5) Based on the research results mentioned above, it needs the future study to continue working and researching for the polymer concentration influencing on the transient flow of the polymer foam in porous media.