Wettability and Production
Wettability and Production
Wettability and Production
Related terms:
Figure 3-7. Displacement of oil by water in a water-wet system. A bank of mobile oil
develops ahead of the advancing water with an average water saturation,SWB, and an oil
saturation quotient develops behind the front from 1-SWCto SORat the point of injection;
the average water saturation behind the front is Sw .
(3.6)
The mobility ratio is a single term that describes the rate and efficiency of oil
displacement by other immiscible fluids. Diminishing values of M(less than one) are
generally considered as favorable values, and increasing values (greater than one) are
considered as unfavorable values. For a waterflood to be economically successful, M
should have a value equal to or less than 10. Enhanced oil recovery processes are
designed to decrease the value of this single parameter.
Miscible displacement EOR processes that use solvents such as light petroleum
fractions, condensed gases, carbon dioxide, and alcohol to dissolve oil at the ad-
vancing front eliminate capillary forces by reducing the interfacial tension to zero.
As the miscible front expands away from the injecting well, more and more oil is
dissolved into the solvent until the solvent is no longer effective and thus making
the displacement less efficient than theoretical predictions from linear laboratory
core floods.
• Relative permeability: If both layers have the same relative permeability charac-
teristics, average water saturation will be higher in the tighter layer than in the
more permeable layer, because the average pressure is always higher in the
less permeable layer.
• Pore size: If pore size in the tight layer is smaller than that in the more perme-
able layer, then it will tend to reduce crossflow. This effect can be estimated
from capillary pressure curves.
• Reservoir geometry: Geometrical nature and extent of interlayer communica-
tion have some effect on observed field performance.
• Permeability anisotropy: In most petroleum reservoirs, vertical permeability is
significantly less than horizontal permeability.
• Reservoir n-layer system: Analysis of performance can be handled to acceptable
accuracy merely by the previously presented formulae and for (kh)t and ( h)t,
respectively.
Principles of Waterflooding
Tarek Ahmed, in Reservoir Engineering Handbook (Fifth Edition), 2019
Step Assuming EA and EV are 100%, calculate cumulative water injected at break-
3. through by applying Equation 14-42:or equivalently:
Step Select several values of injection time less than the breakthrough time, i.e., t
7. < tBT, and set:
Step Calculate the surface water-oil ratio WORs exactly at breakthrough by using
8. Equation 14-28:where fwBT is the wellbore water-cut at breakthrough (notice
that fwBT = fwf).
It should be pointed out that the surface water-oil ratio WORs as calculated by
applying the above expression is only correct when the areal sweep efficiency
EA and vertical sweep efficiency Ev are both 100%. The modification approach of
calculating WORs when EA and EV are 100% are discussed later in the chapter.
Fluid Displacement
John R. Fanchi, in Integrated Reservoir Asset Management, 2010
Exercises
12-1.Find the following values for the example shown in Figure 12.2: fractional
flow of water at the flood front fwf, water saturation at the flood front Swf, and
average water saturation behind the flood front Swbt.
12-2.Consider the following oil–water relative permeability
table.Swkrwkrow0.300.0001.0000.350.0100.5900.400.0200.3200.450.0340.1800.500.046-
0.0800.550.0680.0300.600.1280.0100.650.1660.0010.700.2000.00010.800.2400.0000At
3014.7 psia, oil viscosity = 0.594 cp and water viscosity = 0.503. Calculate and
plot the fractional flow of water.
12-3.Find the following values for the data in Exercise 12-2: fractional flow of water
at the flood front fwf, water saturation at the flood front Swf, and average water
saturation behind the flood front Swbt.
12-4.Find the mobility ratio of water to oil for the data in Exercise 12-2. Is the
mobility ratio favorable or unfavorable?
12-5.The slope of the main tangent line for the data in Exercise 12-2 is 3.145.
Suppose the water injection rate is 400 STB/D, the separation between injector
and producer is 300 feet, the cross-sectional area is 40,000 ft2, and the porosity
is 15 percent. Estimate the time to water breakthrough at the producer using 12-6.
Eq. (12.3.1). Assume the water formation volume factor is 1.01 RB/STB.
Show that Eq. (12.4.8) is a solution of Eq. (12.4.2). 12-7.
Determine the rate of finger growth of a unit mobility flood in a horizontal 12-8.
medium using Eq. (12.4.14). Hint: Set M = 1 in Eq. (12.4.14) and simplify.
Explain why the mobility ratio condition M < 1 is considered “favorable” for a 12-9.
displacement flood using Eq. (12.5.45).
Copy the file VFILL1_HM.DAT to ITEMP.DAT and run IFLO by double click- 12-10.
ing on the IFLO.EXE file on your hard drive. Select option “Y” to write
the run output to files. When the program ends, it will print “STOP”.
Close the IFLO window. You do not need to save changes. Copy the file
ITEMP.ROF to VFILL1_HM.ROF, and copy ITEMP.ARR to VFILL1_HM.ARR.
Open VFILL1_HM.ROF using a text editor. Search the file for INITIAL FLUID
VOLUMES.(a)How much oil is initially in place?(b)How much water is initially
in place?(c)How much gas is initially in place?(d)How much of the gas exists
in a free gas phase?
Run the visualization program 3DVIEW and load the file VFILL1_HM.ARR cre-
ated in Exercise 12-9. To load the file after 3DVIEW is open, click on the “File”
button and select “Open Array File.” Select the file called “VFILL1_HM.ARR”
and click on the “OK” button. Select the oil saturation attribute at the begin-
ning of the run. To select this attribute, click on the “Model” button and select
“Select Active Attribute.” From the list of options select “SO” for oil saturation.
You are looking at the side of the reservoir. To see the top of the reservoir, place
the cursor in the black field near the reservoir display, hold the left mouse
button down, and pull the mouse toward you. You should see the reservoir
image rotate. Continue rotating until you see the top of the reservoir. Sketch
the image and indicate which part of the image represents the reservoir.
In general, therefore, where higher water saturations (following water injection) have
higher water advance rates (Fig. 4.18) we have a self-sharpening system, and where
higher water saturations (following water injection) have lower water advance rates
we have a nonsharpening system (Fig. 4.18).
The applicable part of the curve is therefore that between Sw = Swf and 1 − Sor, where
Swf is the shock front advance water saturation.
To ensure mass balance we need to remove equal areas A and B in the Buckley–Lev-
erett plot to determine Swf—see Fig. 4.18.
A better way is that proposed by Welge. Integrating the saturation distribution from
x = 0 to the shock front, it can be shown that a tangent to the water fractional flow
curve will give both Swbt (the water saturation at the shock front) and the average Sw
behind the shock front (Fig. 4.19).
Figure 4.19. Welge tangent curve.
4.4.1.3.1 Steps
1. Draw the fractional flow curve as shown above.
[4.40]
or in field units:
[4.41]
where q = a constant water injection rate, and q is in bbl/day, t is in days and A is in
ft2.
Welge derived an equation that relates the average displacing fluid saturation to the
saturation at the producing end of the system:
(4.8)
where
Welge also related the cumulative water injected and the water saturation at the
producing end:
(4.9)
Thus, the reciprocal of the slope of the tangent line gives the cumulative water influx
at the time of water breakthrough. When a value of Qi and the injection rate are
known, the time to reach that stage of the flood can be computed.
For a liquid-filled, linear system, the average water saturation at breakthrough, , is:
(4.10)
(4.11)
where Swf is the water saturation at the flood front and fwf is the water cut at the
flood front. After breakthrough, water saturation is obtained from Equations 4.8 and
4.9 where, as mentioned earlier: (1) the tangent point, Sw2, represents the water
saturation at the producing end of the system, (2) the value of fw at the point of
tangency is the producing water cut, (3) the saturation at which the tangent intersects
fw = 1.0 is the average water saturation, and (4) the inverse of the slope of the tangent
line is equal to the cumulative injected fluid in pore volumes (Qi). If connate water
is mobile, appropriate corrections need to be made [44].
(4.12)
After water breakthrough, a number of saturation greater than Swf are selected; the
slope of the tangent line and average water saturation are determined for each value
of Sw chosen. Oil production after breakthrough is then determined by observing the
change in water saturation [3]:
(4.13)
The incremental oil production from Equation 4.13 can be added to the break-
through production from Equation 4.12, and the resulting total production for the
linear system can be listed as a function of Sw, time, or other parameters. If the pore
volumes in these equations are in ft3, divide by 5.615 to get barrels.
The mechanism of increased displacing fluid viscosity can be quantified using the
Buckley-Leverett (1942) theory. Figure 3.5 shows two fractional flow curves. One is
for a waterflooding case with the viscosity ratio of water to oil 0.1, the other one is
for a polymer flooding case with the viscosity ratio of polymer to oil 1. From the
fractional flow curve, we can estimate the average water saturation at breakthrough
by drawing a tangent from the connate water saturation Swc (0.2 in this case) and
intersecting the horizontal line of fw=1, and the corresponding water saturation is
the average water saturation. Here, fw is the water cut in the producing fluid. From
Figure 3.5, the average water saturation in the waterflooding case is 0.58, whereas
the average water saturation in the polymer flooding case is 0.76. The difference is
0.18. In other words, by simply increasing the viscosity of displacing fluid, the oil
recovery factor can be increased by 18% at breakthrough.
One economic impact of polymer flooding which has been less discussed is the
reduced amount of water injected and produced, compared with waterflooding.
Because polymer improves the mobility ratio and sweep efficiency, less water is
injected and less water is produced. In some situations like an offshore environment
and desert area, water and the treatment of water could be costly (Sheng, 2011).
Spontaneous imbibition
James J. Sheng, in Enhanced Oil Recovery in Shale and Tight Reservoirs, 2020
(10.2)
In the above equation, pc is the water-air capillary pressure at the front water
saturation Swf; Sw is the average water saturation behind the front according to the
Handy derivation; kw is the effective water permeability at Sw; is the porosity; A
is the flow cross-section area; and μw is the water viscosity. It is assumed that water
displaces air in a pistonlike manner. No gravity is assumed to play the role in the
process. Only the capillary force overcomes the viscous force within the imbibition
zone. As more water is imbibed, water saturation Sw is increased and kw is increased,
but pc is declined exponentially with Sw. Handy's experimental data through cores
confirmed the above linear relationship. Makhanov's (2013) experimental data also
demonstrated the above relationship, but some imbibition data showed that the
imbibition rate slowed at later time.
Waterflooding
Mohammad Ali Ahmadi, in Fundamentals of Enhanced Oil and Gas Recovery from
Conventional and Unconventional Reservoirs, 2018
Also, the following assumptions have been made for developing Buckley–Leverett
approach:
• Oil and water two-phase flow in reservoir, no gas present in the reservoir