Pressure Maintenance and Improving Oil Recovery by Means of Immiscible Water-Alternating-CO Processes in Thin Heavy-Oil Reservoirs

Download as pdf or txt
Download as pdf or txt
You are on page 1of 12
At a glance
Powered by AI
The document discusses using water-alternating-CO2 injection processes to improve pressure maintenance and oil recovery in thin heavy oil reservoirs. Experiments and numerical simulations were used to evaluate different well configurations and operating parameters.

Experiments using a physical model with horizontal and vertical wells were conducted to evaluate initial oil saturation, production rates, water cuts, recovery, and residual oil distributions. Numerical simulations were then used to match the experimental measurements and optimize parameters like slug size and water/CO2 ratios.

Thermal recovery techniques have issues with heat losses and may leave behind more viscous components in thin reservoirs. Steam injection requires large water consumption. Chemical flooding can have high costs from chemicals and issues with brine salinity and adsorption.

Pressure Maintenance and Improving

Oil Recovery by Means of Immiscible


Water-Alternating-CO2 Processes
in Thin Heavy-Oil Reservoirs
Sixu Zheng, SPE, and Daoyong Yang, SPE, University of Regina, Saskatchewan, Canada

Summary meaning that only a small fraction of oil can be displaced and
Techniques have been developed to experimentally and numeri- recovered, especially when low-temperature oxidation is the dom-
cally evaluate performance of water-alternating-CO2 processes in inant reaction model (Islam et al. 1989; Moore et al. 1999; Shen
thin heavy-oil reservoirs for pressure maintenance and improving 2002). Among the thermal methods proposed for heavy-oil recov-
oil recovery. Experimentally, a 3D physical model consisting of ery, by far the most successful ones are based on steam injection;
three horizontal wells and five vertical wells is used to evaluate however, large water consumption is another major constraint of
the performance of water-alternating-CO2 processes. Two well implementing the steam-based processes in certain areas. Chemi-
configurations have been designed to examine their effects on cal flooding has been considered one of the promising enhanced-
heavy-oil recovery. The corresponding initial oil saturation, oil- oil-recovery (EOR) techniques in heavy-oil reservoirs (Thomas
production rate, water cut, oil recovery, and residual-oil-satura- et al. 2001). For example, polymer flooding has been successfully
tion (ROS) distribution are examined under various operating applied in the Pelican-Lake heavy-oil field (Strauss et al. 2010).
conditions. Subsequently, numerical simulation is performed to In addition to consumption of a large amount of chemical agents,
match the experimental measurements and optimize the operating such as alkali and surfactants, which must be transported and
parameters (e.g., slug size and water/CO2 ratio). The incremental handled safely in the field, several environmental problems may
oil recoveries of 12.4 and 8.9% through three water-alternating- result from chemicals that are injected into the formation (Stoll
CO2 cycles are experimentally achieved for the aforementioned et al. 2010). Conventionally chemical flooding techniques are sen-
two well configurations, respectively. The excellent agreement sitive to brine salinity, leading to high cost or infeasibility in
between the measured and simulated cumulative oil production many cases (Flaaten et al. 2009). High levels of adsorption within
indicates that the displacement mechanisms governing water- the porous medium lead to some chemical EOR techniques asso-
alternating-CO2 processes have been numerically simulated and ciated with surfactants being uneconomical (Fletcher et al. 2012).
matched. It has been shown that water-alternating-CO2 processes Regarding heavy-oil reservoirs discovered in western Canada,
implemented with horizontal wells can be optimized to signifi- although the shallow depth makes it extremely difficult to achieve
cantly improve performance of pressure maintenance and oil re- miscible CO2 displacement, CO2 has been found to be an effec-
covery in thin heavy-oil reservoirs. Although well configuration tive agent for recovering heavy-oil resources through an immisci-
imposes a dominant impact on oil recovery, the water-alternating- ble displacement process, and satisfactory results have been
gas (WAG) ratios of 0.75 and 1.00 are found to be the optimum achieved in heavy-oil reservoirs, such as the Bati Raman field and
values for Scenarios 1 and 2, respectively. Alaska North Slope fields (Issever et al. 1993; Sahin et al. 2008;
Ning et al. 2011). It is of practical and fundamental importance
to evaluate suitability of pressure maintenance and improving
Introduction heavy-oil recovery with immiscible CO2 injection in thin pay
With escalating oil consumption and depleting conventional oil zones in which other EOR techniques are not applicable.
resources, the vast heavy-oil resources discovered in western Can- Compared with natural gas and air, CO2 still has high solubil-
ada (i.e., Alberta and Saskatchewan) have attracted much more ity in heavy oil and bitumen, though it is difficult for CO2 to
attention recently. Approximately 62% of Canada’s total heavy- achieve miscibility with heavy oil. The mechanisms associated
oil resources are located in Saskatchewan, including 3.4 billion with oil recovery include viscosity reduction, swelling effect,
m3 of proven reserves (Bowers and Drummond 1997). As for interfacial-tension reduction, and blowdown recovery (Beecher
these resources, more than 90% of the proven reserves are con- and Parkhurst 1926; Jha 1986; Rojas and Farouq Ali 1988). The
tained in reservoirs with pay zones of less than 10 m, approxi- injected CO2 not only acts as a pressure-maintenance agent but
mately 55% of which have pay zones of less than 5 m (Huang also enhances oil recovery during the pressure-maintenance pro-
et al. 1987). It is a great challenge to recover such heavy-oil cess in heavy-oil reservoirs. In recent years, the WAG process has
resources economically and effectively because of the high vis- been found to be much more effective than gas injection alone,
cosity and fast depletion of reservoir pressure. Thermal recovery because it combines the improved microscopic displacement effi-
techniques, such as cyclic steam stimulation, steamflooding, ciency of the gasflooding with improved macroscopic sweep effi-
steam-assisted gravity drainage, and in-situ combustion, appear to ciency by injecting water (Christensen et al. 2001). Although
be effective for recovering heavy oil because the oil viscosity is numerous immiscible or miscible WAG processes with use of
very sensitive to temperature. However, characteristics of thin res- CO2 (CO2-WAG) have been implemented in light- or medium-oil
ervoirs are such that thermal-based processes are ineffective or reservoirs (Asghari et al. 2007; Righi and Pascual 2007; Shi et al.
uneconomical, mainly because of excessive heat losses to adja- 2008; Chen et al. 2010a, b), few attempts have been made to eval-
cent formations or aquifers. In-situ combustion processes may uate performance of CO2-WAG processes under the well configu-
yield more viscous or immobile components in the formation, rations associated with horizontal wells for pressure maintenance
in a heavy-oil reservoir.
In general, the CO2-WAG performance is significantly
Copyright V
C 2013 Society of Petroleum Engineers
affected by reservoir heterogeneity, rock wettability, fluid proper-
This paper (SPE 157719) was accepted for presentation at the SPE Heavy Oil Conference ties, miscibility conditions, trapped gas, injection techniques, and
Canada, Calgary, 12–14 June 2012, and revised for publication. Original manuscript
received for review 7 April 2012. Revised manuscript received for review 25 September
well operational parameters. Performance of a CO2-WAG process
2012. Paper peer approved 14 November 2012. is largely affected not only by the injection parameters, including

60 February 2013 SPE Reservoir Evaluation & Engineering


TABLE 1—COMPOSITIONAL-ANALYSIS RESULTS OF TABLE 2—COMPOSITION OF SYNTHETIC BRINE
HEAVY-OIL SAMPLE
Component Concentration (mg/L)
Carbon Carbon Carbon
Number (wt%) Number (wt%) Number (wt%) Chloride 43 677
Sodium 24 700
C1 0.00 C11 1.26 C21 2.39 Calcium 1500
C2 0.00 C12 1.54 C22 1.44 Magnesium 850
C3 0.00 C13 1.95 C23 1.98 Potassium 413
C4 0.00 C14 2.18 C24 1.67
C5 0.00 C15 2.55 C25 1.67
C6 0.00 C16 2.29 C26 1.70 21.0 C and 3500 kPa, respectively. The density of the oil sample
C7 0.00 C17 2.29 C27 1.72 is measured as 976.0 kg/m3 at 15.6 C and atmospheric pressure;
viscosity of the oil sample was 5,861 cp at 20.0 C and atmos-
C8 0.00 C18 2.60 C28 1.75
pheric pressure when measured with a Brookfield DV-IIþ Pro-
C9 0.96 C19 2.23 C29 1.50
grammable Viscometer (Brookfield Engineering Laboratories).
C10 0.99 C20 2.02 C30þ 61.33 The asphaltene content of the heavy-oil sample was determined as
14.4 wt% (n-pentane insoluble).
The compositional analysis of the oil sample is tabulated in
WAG ratio (i.e., water-slug size to CO2-slug size), injection rate, Table 1. The heavy components (C30þ) account for 61.3 wt%,
and cycle time, but also by the production parameters, including whereas no hydrocarbon component lighter than C9 exists in this
production rate and bottomhole pressure at the producer (Christen- heavy oil. Density of reservoir brine is measured as 1051.7 kg/m3
sen et al. 2001; Chen et al. 2010a, b). In particular, WAG ratio is a at 20.0 C; chloride, sodium, calcium, and magnesium ions are the
crucial factor that affects the oil recovery and has an optimum dominant components. The brine viscosity is determined as 1.03
value in a hydrocarbon reservoir (Reid and Robinson 1981; Rojas cp at 20.0 C. In this study, a synthetic brine was prepared and
and Farouq Ali 1986; Chen et al. 2010a, b; Torabi et al. 2010). An used as a proxy for reservoir brine, the compositions of which are
improved oil recovery in a miscible CO2-WAG process is provided in Table 2.
achieved as the result of an increase in sweep efficiency, and the In this study, research grade CO2 with a purity of 99.998%
highest oil recovery is achieved when the overall WAG ratio for was purchased from the Praxair, Canada. The Computer Model-
the entire reservoir reaches 0.5 (Chen et al. 2010a, b). Regarding ling Group (CMG) WinProp module (Version 2009.11, Computer
the CO2-WAG application in a heavy-oil reservoir, a high rather Modeling Group Limited, Canada) with the Peng-Robinson (PR)
than low WAG ratio is recommended, and a WAG ratio of 4:1 or equation of state (EOS) (Peng and Robinson 1976) was used to
5:1 is found to be the optimum value (Rojas and Farouq Ali 1986). calculate the densities of CO2 at different pressures and tempera-
In this study, a 3D physical model with five vertical wells and tures. Ottawa sand No. 710 (Bell & Mackenzie Company Lim-
three horizontal wells is developed and applied to experimentally ited) was used to pack the 3D physical model, and its screen
evaluate the performance of pressure maintenance and improving analysis is listed in Table 3.
oil recovery by means of CO2-WAG processes in thin heavy-oil
reservoirs. Experimentally, displacement experiments with a given
well configuration are started first with waterflooding followed by Experimental Setup
CO2-WAG processes, to maintain reservoir pressure and improve A digital image of the experimental setup is shown in Fig. 1. The
oil recovery. In addition to recording the pressure drop, the pro- whole setup can be divided into three subsystems: fluid-supply
duced oils are collected and measured, and gas-flow rate and system, 3D displacement system, and fluid-production system.
its cumulative volume are measured and logged automatically The major component of this experimental setup is the 3D dis-
throughout the experiments. Theoretically, numerical simulation placement model with dimensions of 304  304  127 mm. The
is performed to match the production history and serves as a basis wall of the 3D model has a thickness of 25.4 mm. The five-spot-
for future optimization of CO2-WAG processes. Once the simula- pattern size covers an area of 253  253 mm, which represents a
tion model is built, history matching is performed to match the ex- typical well spacing of approximately 3.2  104 m2 (8 acres) in
perimental measurements. Subsequently, sensitivity analysis was heavy-oil fields (Rojas and Farouq Ali 1986). Twenty-eight bolts
performed to examine the effect of water-slug size. are used to fasten the steel lid, and a rubber gasket is applied for
sealing the model. Three parallel horizontal wells are at the bot-
tom of the void space, and five vertical wells are distributed as a
Experimental Materials conventional five-spot well pattern in the model. This 3D dis-
The heavy-oil sample and reservoir brine were from the Lashburn placement model can be used to perform experiments under five-
area, Canada. The initial reservoir temperature and pressure were spot, line-drive, and other well patterns associated with horizontal
wells. All the wells have a diameter of 6.4 mm, around which
small holes are perforated to allow fluids to flow. The wells are
TABLE 3—TYPICAL SCREEN ANALYSIS OF OTTAWA SAND
wrapped with a 250-mesh wire-screen (Ferrier Wire Goods, Can-
NO. 710
ada) to prevent sand production during the tests. The maximum
operating pressure of the 3D displacement system is 6.0 MPa at
US Standard Sieve Size
room temperature because the rectangular model shows a tend-
Retained ency to bulge under a high internal pressure. In this study, a leak-
Mesh mm (wt%) age test was conducted up to the pressure of 5.0 MPa at room
temperature for 30 minutes before any experiments.
40 0.425 Trace Two well configurations have been designed to examine their
50 0.300 0.3 effects on heavy-oil recovery in this study. Scenario 1 includes
70 0.212 7.0 two parallel horizontal wells, one of which is injector and the
100 0.150 52.0 other is producer. Scenario 2 has four vertical injectors with one
140 0.106 32.0 horizontal producer. Both of these scenarios are illustrated in
Fig. 2. Porosities of Scenarios 1 and 2 are determined to be 33.1
200 0.075 8.0
and 33.5%, respectively. Absolute permeabilities of Scenarios 1
270 0.053 0.7
and 2 are determined to be 1,700 and 1,800 md, respectively. The

February 2013 SPE Reservoir Evaluation & Engineering 61


Syringe pump

Transfer cyclinder Data acquisition

Gas flow meter

CO2

3D model

BPR Collector

Fig. 1—Digital image of the experimental setup used in this study.

initial oil saturations for Scenarios 1 and 2 are measured to be The porosity of the sandpacked model is measured by means
93.1 and 92.5%, respectively. of the imbibition method. After completion of the displacement
process, the blowdown process is carried out. Once a blowdown
Experimental Procedures recovery is completed, ROS is measured. From top to bottom, the
sandpack model is divided into four layers and marked as Layers
Before conducting the experiments, several preparations must be
1 through 4 in sequence, each of which is approximately 3.0 cm
completed, including sandpacking, porosity measurement, perme-
thick. In each layer, mixtures of sand and reservoir fluids are col-
ability determination, and initial-oil-saturation determination. As
lected from 13 locations to measure the ROS. First, the weight of
for the sand-packing process, the 3D displacement model is placed
each collected sample that contains sand, residual oil, and water is
horizontally before filling the void space with a given amount of
measured. Second, the collected sample is heated in an oven at
the Ottawa sand. Subsequently, the 3D model is shaken for a tight
100.0 C for 4 hours to ensure that water is completely evaporated.
packing with a pneumatic vibrator (NP 35, Northern Vibrator) for
Then the remaining sand/oil mixture is weighed. It should be
at least 10 hours. The lid is then removed, and more sand is added
noted again that the heavy-oil sample used in this study does not
to the formed void space in the model; the model is then covered,
contain any hydrocarbon components lighter than C9, of which
tightened, and shaken again. The same process is normally re-
the normal boiling point is 143.0 C. As such, hydrocarbon com-
peated four to five times until no void space is observed at the top
ponents will not be distilled when the heating temperature is set to
of the physical model. During the packing process, the bottom of
be 100.0 C for evaporating water. Third, toluene and kerosene are
the sandpack is packed more tightly than the top layer.
used as solvents to completely remove the residual heavy oil in
the sample. Finally, the sand is heated in the oven at 110.0 C and
then weighed when it is completely dried. Consequently, ROS is
calculated on the basis of the weight differences, porosity of the
sandpack model, oil density, and sand density. Detailed proce-
dures for the aforementioned preparations can be found elsewhere
(Zheng et al. 2011).
The sandpacked 3D physical model is first flooded with oil
until irreducible water saturation is reached. The physical model
is then left undisturbed for 24 hours to equilibrate the distribution
of fluids. For each scenario, water injection is started from the
(a) injector(s) with a constant injection rate of 2.0 cm3/min, which is
in accordance with the typical average pore velocity of 1.0 ft/D in
an oil field. Water injection is terminated once no more oil is
produced.
After the waterflooding process with a specific well configura-
tion, several cycles of water and CO2 are injected into the physi-
cal model alternatively to evaluate their performance of pressure
maintenance and improving oil recovery. During the CO2-WAG
processes, CO2 is introduced into the physical model from a CO2
cylinder, in which a gas regulator is used to control the injection
(b) pressure of 3900 kPa. A backpressure regulator (BPR) (EB1HP1,
Equilibar) is connected with the production port to maintain the
Fig. 2—Schematic of well configurations. (a) Scenario 1: one physical model at the prespecified operating pressure of 3860 kPa
horizontal injector and one horizontal producer; (b) Scenario 2: during WAG processes. Therefore, the pressure inside the 3D dis-
four vertical injectors and one horizontal producer. placement model has been maintained at approximately 3860 kPa

62 February 2013 SPE Reservoir Evaluation & Engineering


10 1.20
Calculated viscosity with GA-based correlation
Calculated swelling factor with volume-translated PR EOS
Calculated swelling factor with GA-based correlation

1.15
1

Viscosity, 103 cP

Swelling factor
1.10

0.1
1.05

0.01 1.00
0.0 0.2 0.4 0.6
Solubility of CO2 in heavy oil, mole fraction

Fig. 3—Calculated viscosity (Emera and Sarma 2008) and swelling factor (SF) (Li et al. 2011) of CO2-diluted heavy oil as a function
of CO2 solubility.

with the use of the gas regulator and BPR during the displacement Because it is not reliable to generate the EOS model for composi-
process. All the experiments are performed at the ambient temper- tional simulator without sufficient laboratory pressure/volume/
ature of 21.0 6 1.0 C. temperature data, as recommended by the CMG, genetic-algo-
Conical-bottom glass centrifuge tubes (Kimble) are used to rithm (GA) -based correlations (Emera and Sarma 2008) have
collect the produced fluids. The capacity of the tube is 100 cm3. been used to calculate CO2 solubility, swelling factor, and the vis-
There are several subdivisions, among which the largest toler- cosity of the oil/CO2 mixture in this study. The GA-based viscos-
ance is 1.00 cm3. The volume of produced liquids can be meas- ity correlation can reproduce 130 viscosity data for oil/CO2
ured accurately through the use of the tubes. A gas flowmeter mixtures, with an accuracy of 6.0%, whereas the swelling-factor
(XFM17S, Aalborg) is connected with the collector to measure correlation is able to reproduce 85 swelling-factor data for oil/
the produced gas. The gas-flow rates and gas cumulative volume CO2 mixtures, with an accuracy of 0.61% (Emera and Sarma
are logged and stored automatically in a computer at a preset time 2008). Fig. 3 shows the predicted viscosity and swelling factor of
interval. an oil/CO2 mixture, as a function of CO2 solubility through the
use of the GA-based correlations. A significant viscosity reduction
and swelling effect can be achieved with the dissolution of CO2 in
Numerical Simulation heavy oil (Beecher and Parkhurst 1926; Jha 1986). Fig. 3 also
To better understand the mechanisms of CO2-WAG processes shows the calculated swelling factor with the use of a newly
with the purpose of pressure maintenance in heavy-oil reservoirs, developed binary-interaction-parameter correlation in the vol-
numerical simulation has been conducted to history match the ume-translated PR EOS, which has been experimentally validated
waterflooding and CO2-WAG processes. As mentioned previ- with an oil sample collected from the same area (Li et al. 2011).
ously, no hydrocarbon components lighter than C9 are found in As seen in Fig. 3, there exists an excellent agreement between the
the dead heavy-oil sample in this study, whereas heavy compo- predicted swelling factor by use of the GA-based correlations and
nents account for a large proportion of the composition. The oper- the predicted swelling factor with the PR EOS model. This also
ating pressure (less than 4000 kPa) is far below the minimum confirms the reliability of the GA-based correlations used in this
miscibility pressure of 15.5 MPa at temperature of 21.0 C, which study. Because a black-oil model has also been considered to be
is calculated by the CMG WinProp module (Version 2009.11). an effective tool for simulating performance of immiscible flood-
Neither light-components extraction nor miscible flooding has ing processes (Cobanoglu 2001), the IMEX black-oil simulator
been achieved at the experimental conditions in this study. (Version 2009.11, CMG Limited, Canada) has been used in this
study. A grid system of 13  13  6 gridblocks is used in the nu-
merical model to represent the physical model, resulting in a fine-
Injector #1
gridblock size of 2.34  2.34  2.10 cm (see Fig. 4).
Injector #3 Uniform distribution of porosity and oil saturation is used in
Injector #2 10.6 the model; the absolute permeability of each layer varies with the
9.5 depth. Because of the relatively high absolute permeability, the
Injector #4 capillary pressure effect is neglected. Because the relative perme-
8.5
abilities cannot be directly measured from experiments, they are
7.4 to be determined by use of the history-matching techniques
6.4 instead (Watson et al. 1980; Eydinov et al. 2009; Li et al. 2009;
5.3 Li and Yang 2011). In principle, trapping of the nonwetting phase
4.2 (CO2) significantly affects the hysteresis in relative permeability
Producer (Spiteri and Juanes 2004). It is experimentally observed, however,
3.2
that not much CO2 is trapped in the sandpacked model during the
2.1 immiscible CO2-WAG process. It is reasonably assumed that hys-
1.1 teresis effect is not significant in this study, and thus the hysteresis
0.0 effect on relative permeability has not been taken into account. In
this study, two sets of similar relative permeability curves are
Fig. 4—3D view of the simulation model of Scenario 2. used to accurately match the experimental measurements, though

February 2013 SPE Reservoir Evaluation & Engineering 63


40.0 Cycle 2 CO2 H2O 100
Cycle 1
CO2 H2O
H2O
CO2 80
30.0
Oil recovery, %OOIP

H2O
CO2

Water cut, %
CO2 H2O 60
CO2 H2O
20.0
Cycle 1 Cycle 2 Cycle 3
40
Scenario #1
Scenario #2
10.0
Scenario #1 20
Scenario #2

0.0 0
0.0 1.0 2.0 3.0 4.0 5.0 0.0 1.0 2.0 3.0 4.0 5.0
Injected fluids, PV Injected fluids, PV

Fig. 5—Oil recovery during displacement experiments for Sce- Fig. 6—Water cut during displacement experiments for Scenar-
narios 1 and 2. ios 1 and 2.

one set of the average relative permeability curves would provide water is able to mobilize the contacted oil that is to be displaced
a reasonable match. Regarding field applications, however, one to the producer along these paths (Mai and Kantzas 2009).
set of average relative permeability curves is used because these As shown in Fig. 5, oil recovery of the two scenarios stays
two sets of relative permeability curves are similar. close at the early stage. With a further increase of injected water
To simulate the experimental displacement process, the well PV, the well configuration associated with two horizontal wells
operational conditions defined in simulation models are set the (i.e., one horizontal injector and one horizontal producer, Scenario
same as those used in the experiments. For example, in the experi- 1) tends to yield a higher oil recovery compared with Scenario 2.
ments, the production port is connected with the BPR that is used For instance, oil recovery of Scenario 1 is approximately 1.3%
to maintain the pressure of model at the prespecified operating higher than that of Scenario 2 when 1.00 PV of water is injected.
pressure of 3860 kPa. In the simulation models, the bottomhole This can be attributed to the fact that a horizontal well is efficient
pressure of the producer is then set to be 3860 kPa. During a water- to achieve a better sweep efficiency.
injection process, the water-injection rate in the simulation model
is set at the constant value of 2.0 cm3/min, which is the same as the
CO2-WAG Processes. After waterflooding, three cycles of water
experimental water-injection rate. For the CO2-injection process, a
and CO2 injection are alternately implemented to examine per-
gas regulator connected with the CO2 cylinder is used to control
formance of pressure maintenance and improving oil recovery.
the injection pressure of 3900 kPa. Accordingly, the CO2-injection
During CO2-WAG processes, a BPR is used to maintain the pres-
pressure is set to be 3900 kPa in the simulation model.
sure of the physical model at 3860 kPa, which is higher than the
original reservoir pressure of 3500 kPa, whereas the CO2-injec-
Results and Discussion tion pressure of 3900 kPa is controlled through a gas regulator
Waterflooding. Waterflooding is started in Scenarios 1 and 2. connected with the CO2 cylinder. For Scenario 1, the same slug
The oil-recovery profiles for the two scenarios are shown in sizes (0.40 PV) of water and CO2 are used; that is, the WAG ratio
Fig. 5. After 1.00 pore volume (PV) of water injection, oil recov- is 1:1. The slug sizes of water and CO2 of Scenario 2 are 0.25 and
eries of 24.3 and 23.1% are achieved for Scenarios 1 and 2, 1.00 PV, respectively, resulting in a WAG ratio of 1:4.
respectively. A similar recovery profile for heavy-oil waterfloods Fig. 5 illustrates that the majority of the produced oil in the
can be found in the literature (Mai and Kantzas 2009). Fig. 6 CO2-WAG processes is recovered during the first cycle. After the
shows water cut in the produced liquid during displacement proc- second cycle, there is no significant increment in oil recovery.
esses of the two scenarios. The fact that the heavy-oil sample is Similar results are also documented in the literature (Sohrabi et al.
dramatically more viscous than the injected water leads to a high 2000). It is obvious that the oil recovery of Scenario 1 levels out
mobility ratio of water to heavy oil. In this way, viscous fingering during the third cycle, whereas an additional 2.0% in oil recovery
is severe during waterflooding processes, resulting in a very early is achieved during the third cycle for Scenario 2. As depicted in
water breakthrough and low oil recovery. Fig. 5, for these CO2-WAG cycles, CO2-injection processes
For light- and medium-oil waterfloods, as water is injected, oil recover a small amount of heavy oil, whereas more oil is pro-
is continuously produced until breakthrough, and a relatively high duced during the subsequent water-injection processes. Although
recovery factor can be achieved before breakthrough (Hadia et al. there is not much heavy oil produced during the water-injection
2007). As shown in Fig. 6, however, it is obvious that break- phase, the injected CO2 still contributes to oil recovery. The
through occurs very early in a heavy-oil waterflooding process, injected CO2 dissolves into heavy oil, leading to a swelling effect
and water cut increases sharply to high values at the early stage. and oil-viscosity reduction. In addition, because gasflooding leads
For these two scenarios, water cut approximates 60% at 0.15 PV to higher microscopic displacement efficiency than that of water-
of injected water, whereas water cut for Scenarios 1 and 2 is flooding, the injected CO2 changes the distribution of liquids in
measured as 95.6 and 97.5%, respectively, at 1.00 PV of injected the pores and mobilizes much more oil that is not previously con-
water. Although more heavy oil is produced after 0.20 PV of tacted by water. Subsequently, the injected water displaces the
water injection, water cut is high. This means that the injected mobilized oil to the producer. In this study, Scenario 1, with a
water may still keep contacting the in-place oil and displacing the larger water-slug size, achieves a higher oil recovery than Sce-
oil to the producer. The difference of water cut is minor, whereas nario 2, with a smaller water-slug size, when 3.0 PV of fluids is
the high water cuts indicate that much water is produced, and injected. Oil recovery of Scenario 1 reaches the plateau in a short
waterflooding is found to be inefficient at the late stage of water period, demonstrating that it reaches the production limit quickly.
injection. Figs. 5 and 6 demonstrate that considerable heavy oil is In Scenario 2, oil is still produced in the third cycle, indicating
still produced after water breakthrough. This is ascribed to the that production can be sustained for a much longer time. The bet-
fact that water flows through the continuous low-resistance path- ter performance of Scenario 1 is partially attributed to the fact
ways that are generated by water after breakthrough. The flowing that a significantly larger drainage volume can be controlled by a

64 February 2013 SPE Reservoir Evaluation & Engineering


100
Scenario #2
Scenario #1
80
Water cut, %

60

40

20

0
0.0 1.0 2.0 3.0 4.0 5.0
Blowdown recovery, %OOIP

Fig. 7—Water cut as a function of blowdown recovery for Sce- (a)


narios 1 and 2.

horizontal well compared with a vertical well because of the


larger contact area of horizontal wells, consequently leading to a
larger sweep efficiency of the CO2-WAG process. Because the
injected CO2 has a tendency to override in the pay zone, a hori-
zontal well at the bottom of the pay zone can alleviate gas over-
ride and more efficiently recover the heavy oil that is produced
through a gravity-drainage mechanism.
The water-cut profiles of the produced liquids during CO2-
WAG processes are also shown in Fig. 6. It is noticeable that water
cut decreases significantly during the CO2 injection processes, espe-
cially in the first cycle, while, in the second and third cycles, water
cut is also decreased during the CO2 injection phases. Because
more heavy oil has been displaced, whereas more water has been
injected into the sandpack model after the first WAG cycle, the dec-
rement of water cut is not as significant as that of the first cycle. Af-
ter CO2 injection, water cut increases with the further water
injection, whereas it oscillates rather than keeps an approximately
constant value during the subsequent water-injection processes. (b)

Blowdown Recovery. After three cycles of water and CO2 are Fig. 8—Distribution of ROS in the top layer for (a) Scenario 1
injected into the 3D physical model under each well configura- and (b) Scenario 2.
tion, pressure inside the model is decreased to atmospheric pres-
sure for blowdown recovery. As shown in Fig. 7, the blowdown ROS in the top layer can be observed directly. Digital images of
process contributes an increment of 2.8 and 4.1% in oil recovery the top layer for the two scenarios are shown in Fig. 8. In Scenario
for Scenarios 1 and 2, respectively. The blowdown recovery is 1 (see Fig. 8a), it is obvious that the color of sand along the hori-
mainly attributed to the solution-gas-drive mechanism (Jha 1986). zontal producer (Well 2) is light, demonstrating that heavy oil
CO2 gas has a good ability to dissolve into oil. During the blow- along the horizontal producer has been displaced effectively. It
down process, with pressure decline, energy stored by CO2 disso- also can be seen that the color of the bottom-left portion for Sce-
lution is then released and continues to drive the oil to the nario 1 is much darker than that of the top-right portion, indicat-
wellbore (Mangalsingh and Jagai 1996). Besides the expanded ing that a stable displacement front along the horizontal injector
gas exerting a displacement force to drive oil to producers, such a was difficult to establish. A large portion of injected fluids invades
high efficiency of foamy-oil flow under solution-gas drive in vis- the sandpack through the heel of the horizontal injector and then
cous oil is mainly the result of (1) a significant reduction in oil displaces the heavy oil to the producer, resulting in a well-
viscosity, because small gas bubbles are formed in the oil, and (2) displaced triangle area (i.e., top-right portion). This observation
simultaneous flow of continuous oil phase and discontinuous gas implies that it is possible to improve the sweep efficiency by
phase in the form of tiny bubbles (Smith 1988). changing the position of heels and toes of these two parallel hori-
It is experimentally found that little fluid is produced when zontal wells. In Scenario 2 (see Fig. 8b), the color of sand for the
pressure is decreased from 3800 kPa to 2000 kPa, whereas most top layer is very dark, indicating that much heavy oil is still at the
of the fluids are produced at pressure less than 1500 kPa. The top layer and the heavy oil at the top layer cannot be swept effi-
blowdown process is found to last for a relatively long period ciently, although approximately 5.00 PV of fluids has been
(i.e., more than 2 hours) in this study. Approximately 40 L of CO2 injected. Such observed dark-sand distribution in the 3D model
(ambient temperature and atmospheric pressure) is produced dur- corresponds to its low oil recovery with the displacement proc-
ing the blowdown process, whereas gas is still producing even esses previously applied.
though the pressure of the physical model approximates to the To examine the ROS in the physical model, ROS is measured
atmospheric pressure. This implies that CO2 storage also occurs in this study. Distribution of ROS in the sandpack for Scenarios 1
as the injected CO2 dissolves into reservoir fluids and occupies and 2 is illustrated in Figs. 9 and 10, respectively. As mentioned
the pore spaces that were previously filled by reservoir fluids. previously, from top to bottom, the sandpack is divided into four
layers, marked as Layers 1 through 4 in sequence. For each layer,
Distribution of ROS. After the blowdown process, the steel lid 13 samples were collected and measured from different locations
on the top of the 3D model is removed and the distribution of that are marked by the circle symbols in Figs. 9 and 10. The ROS

February 2013 SPE Reservoir Evaluation & Engineering 65


0.20 0.20
10.0 0.25 10.0 0.25
0.30 0.30
0.35 0.35
0.40 0.40
0.45 0.45
5.0 0.50 5.0 0.50
0.55 0.55
0.60 0.60
0.65 0.65
Y, cm

Y, cm
0.0 0.0

–5.0 –5.0

–10.0 –10.0

–10.0 –5.0 0.0 5.0 10.0 –10.0 –5.0 0.0 5.0 10.0
X, cm X, cm
(a) (b)

0.20 0.20
10.0 0.25 10.0 0.25
0.30 0.30
0.35 0.35
0.40 0.40
0.45 0.45
5.0 0.50 5.0 0.50
0.55 0.55
0.60 0.60
0.65 0.65
Y, cm

Y, cm

0.0 0.0

–5.0 –5.0

–10.0 –10.0

–10.0 –5.0 0.0 5.0 10.0 –10.0 –5.0 0.0 5.0 10.0
X, cm X, cm
(c) (d)

Fig. 9—Distribution of ROS for Scenario 1 in (a) Layer 1, (b) Layer 2, (c) Layer 3, and (d) Layer 4.

contour of each layer is then generated on the basis of these meas- tem and a gas/liquid system are determined by history matching
ured values. As shown in Fig. 9, for Layer 1, the bottom-left por- the experimental measurements (Watson et al. 1980; Spiteri and
tion shows a higher ROS, which is consistent with the dark color Juanes 2004; Eydinov et al. 2009; Li et al. 2009, 2011). More spe-
of the bottom-left portion in the digital image for Scenario 1 (see cifically, the initial inputs of the three-phase relative permeability
Fig. 8a). For Layer 2 of Scenario 1, the ROS along the horizontal curves are generated with the commonly used correlation (i.e.,
producer is low, indicating heavy oil along the producer has been Stone II correlation). Subsequently, the relative permeability
efficiently displaced. That is attributed to the fact that the horizon- curves are tuned and determined once the experimented results
tal producer controls a large reservoir area and has high sweep effi- have been history matched. Such determined relative permeability
ciency. The average ROSs of Layers 3 and 4 for Scenario 1 are curves for two scenarios are shown in Figs. 11 through 14. As
relatively low (See Fig. 9). This is ascribed to the fact that two par- can be seen from Figs. 11 and 13, water relative permeability
allel horizontal wells are at the bottom of the physical model, lead- increases slightly when the water saturation is less than 0.3. The
ing to high sweep efficiency in the bottom layers of the model. relative permeability values are low at a high water-saturation
As shown in Fig. 10 for Scenario 2, Layer 1 shows that the range. The performance of waterflooding is very sensitive to the
ROS in the top layer is high, which is consistent with the dark relative permeability values at a low water saturation range (less
color as seen in Fig. 8b, whereas the average ROS decreases from than 0.2). For gas/liquid relative permeability curves (see Figs. 12
top to bottom layers. Layer 4 illustrates low ROS; a clear dis- and 14), oil relative permeability values decrease sharply. This
placement front from four vertical injectors to the horizontal pro- implies that gas fingering is severe because of the high mobility
ducer is formed because the horizontal producer is at the bottom ratio of gas to heavy oil, resulting in an early breakthrough. The
of the model in Scenario 2. Therefore, the horizontal producer at relative permeability of oil reduces significantly as gas saturation
the bottom layer is beneficial for recovering the bottom heavy oil increases.
in reservoirs with gas-cap drive. The oil viscosity also affects the heavy-oil/water relative per-
meability curves. With an increase in oil viscosity, ROS is
increased and relative permeabilities are decreased at a higher
Numerical Simulation and Performance Optimization. water-saturation range (Wang et al. 2006; Nejad et al. 2011). Nu-
Because no direct interpretation methods have been made avail- merical simulation of the displacement processes is conducted
able to determine relative permeability for such displacement with injection rate, injection pressure, and production pressure as
experiments, the relative permeability curves for a water/oil sys- input constraints. The cumulative oil production and water cut are

66 February 2013 SPE Reservoir Evaluation & Engineering


0.20 0.20
0.25 0.25
10.0 0.30 10.0 0.30
0.35 0.35
0.40 0.40
0.45 0.45
0.50 0.50
5.0 0.55 5.0 0.55
0.60 0.60
0.65 0.65
Y, cm

Y, cm
0.0 0.0

–5.0 –5.0

–10.0 –10.0

–10.0 –5.0 0.0 5.0 10.0 –10.0 –5.0 0.0 5.0 10.0
X, cm X, cm
(a) (b)

0.20 0.20
0.25 0.25
10.0 0.30 10.0 0.30
0.35 0.35
0.40 0.40
0.45 0.45
0.50 0.50
5.0 0.55 5.0 0.55
0.60 0.60
0.65 0.65
Y, cm
Y, cm

0.0 0.0

–5.0 –5.0

–10.0 –10.0

–10.0 –5.0 0.0 5.0 10.0


–10.0 –5.0 0.0 5.0 10.0
X, cm
X, cm
(c)
(d)

Fig. 10—Distribution of ROS for Scenario 2 in (a) Layer 1, (b) Layer 2, (c) Layer 3, and (d) Layer 4.

compared with the experimental data in Figs. 15 and 16. The indicating that the waterflood becomes less efficient afterward.
good matching results show that the numerical simulation cap- After CO2 injection, however, a large amount of oil is still pro-
tured the overall behavior of both the waterflooding and CO2- duced during the subsequent water-injection process. This is
WAG processes. As can be seen in Figs. 15 and 16, oil recovery attributed to the fact that the injected CO2 with high microscopic
approaches the plateau after 1.00-PV injection (2,000 minutes), sweep efficiency can swell the heavy oil, reduce oil viscosity, and

1.0
1.0 Krg
Krog
Krw
Krow 0.8
0.8

0.6
0.6
Kr
Kr

0.4
0.4

0.2
0.2

0.0
0.0 0.0 0.2 0.4 0.6 0.8 1.0
0.0 0.2 0.4 0.6 0.8 1.0 Sg
Sw
Fig. 12—Liquid and gas relative permeability curves for Sce-
Fig. 11—Water and oil relative permeability curves for Scenario 1. nario 1.

February 2013 SPE Reservoir Evaluation & Engineering 67


1.0 1.0
Krw Krg
Krow Krog
0.8 0.8

0.6 0.6
Kr

Kr
0.4 0.4

0.2 0.2

0.0 0.0
0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0
Sw Sg

Fig. 13—Water and oil relative permeability curves for Scenario 2. Fig. 14—Liquid and gas relative permeability curves for Sce-
nario 2.

1600 100

1400
80
Cumulative oil production, cm3

1200

1000

Water cut, %
60

800

40
600
Measured water cut
400 Simulated water cut
Measured oil production 20
Simulated oil production
200

0 0
0 1000 2000 3000 4000 5000
Time, min

Fig. 15—Measured and simulated cumulative oil production and water cut for Scenario 1.

1600 100

1400
80
Cumulative oil production, cm3

1200

1000
Water cut, %

60

800

40
600
Measured water cut
Simulated water cut
400
Measured oil production 20
Simulated oil production
200

0 0
0 1000 2000 3000 4000 5000
Time, min

Fig. 16—Measured and simulated cumulative oil production and water cut for Scenario 2.

68 February 2013 SPE Reservoir Evaluation & Engineering


1600 12

1400

Cumulative oil production, cm3


10
1200

Cumulative WOR
8
1000

800 6
600 WAG ratio: 2.00; Injected PV: 3.6
WAG ratio: 1.00; Injected PV: 3.6 4 WAG ratio: 2.00; Injected PV: 3.6
400 WAG ratio: 0.75; Injected PV: 3.5 WAG ratio: 1.00; Injected PV: 3.6
WAG ratio: 0.50; Injected PV: 3.6 WAG ratio: 0.75; Injected PV: 3.5
WAG ratio: 0.25; Injected PV: 4.0 2 WAG ratio: 0.50; Injected PV: 3.6
200 Continuous Waterflooding; Injected PV: 2.0 WAG ratio: 0.25; Injected PV: 4.0
Experiment
0 0
0 2000 4000 6000 8000 10000 0 2000 4000 6000 8000 10000
Time, min Time, min
Fig. 17—Production profiles with different water-slug sizes for Fig. 18—Cumulative WOR with different water-slug sizes for
Scenario 1. Scenario 1.

mobilize the heavy oil in the pores, resulting in a higher oil recov- ratio of 0.5. As such, WAG ratio of 0.75 leads to a relatively higher
ery with the subsequent water injection that displaces the mobi- oil recovery and a lower cumulative WOR. For Scenario 2, the
lized heavy oil. Therefore, the CO2-WAG process is a promising slug size of CO2 is set to be 1.00 PV for all cases (see Figs. 19 and
approach to maintain pressure and improve oil recovery in heavy- 20). The production profiles for different WAG ratios follow a
oil reservoirs because of the positive and composite effects of oil similar trend, whereas a higher water-slug size yields a higher oil
swelling and viscosity reduction by CO2 dissolution, high micro- recovery. The water-slug sizes of 1.00 PV and 2.00 PV (i.e., WAG
scopic sweep efficiency of CO2, and good macroscopic sweep ef- ratios of 1.00 and 2.00, respectively) are sufficient to control the
ficiency of water. mobility of CO2 to achieve high oil recovery. However, the WAG
As mentioned previously, oil is mainly produced during the ratio of 1.0 can keep a lower cumulative WOR than that of WAG
water-injection process for each CO2-WAG cycle. In this study, ratio of 2.00 when the same volume of fluids is injected (6.00 PV).
effects of water-slug size are examined by means of the history- Therefore, the WAG ratios of 0.75 and 1.00 are determined to be
matched reservoir models for Scenarios 1 and 2, respectively. In the optimum values for Scenarios 1 and 2, respectively.
Scenario 1, the slug size of CO2 is kept to 0.40 PV for each cycle, It is found from both the experimental measurements and nu-
whereas the slug size of water is changed from 0.10 to 0.80 PV merical simulation that well configuration imposes a dominant
(i.e., WAG ratio varies from 0.25 to 2.00). Although the CO2 slug impact on oil recovery. The well configuration associated with
size of Scenario 2 (1.00 PV) is larger than that of Scenario 1 (0.40 two horizontal wells tends to yield a higher oil recovery under the
PV), oil recovery of Scenario 2 is still low. For convenient com- same WAG ratio. Regarding each well configuration, WAG ratio
parison, continuous waterflooding is used as a reference case. As also affects oil recovery and cumulative WOR, although there
shown in Fig. 17, compared with other schemes, the scheme with exists an optimum value for each scenario.
a WAG ratio of 0.25 produces less oil, though 4.00 PV of fluids
has been injected. This is attributed to the fact that compared with
the CO2-slug size, the water-slug size is too small to effectively Conclusions
control the CO2 mobility in heavy-oil reservoirs. Regarding the
same CO2 slug size, a larger water-slug size yields a higher oil re- Techniques have been developed to experimentally and numeri-
covery. As shown in Fig. 18, if the water-slug size is very large, cally evaluate the performance of CO2-WAG processes under
the extra amount of water injected leads to trapping of oil and well configurations associated with horizontal wells in thin heavy-
high cumulative water/oil ratios (WORs) (Rojas and Farouq Ali oil reservoirs for pressure maintenance and improving oil recov-
1986). It should be noted that oil recovery increases when the ery. Experimentally, a 3D physical model is used to examine the
WAG ratio is increased from 0.50 to 2.00, whereas the cumulative suitability of pressure maintenance and improving oil recovery by
WOR at WAG ratio of 0.75 remains close to the value at WAG means of an immiscible CO2-WAG process in heavy-oil

12
1600

1400 10
Cumulative oil production, cm3

1200
Cumulative WOR

8
1000

800 6

600 WAG ratio: 2.00; Injected PV: 6.00


4
WAG ratio: 1.00; Injected PV: 6.00
WAG ratio: 0.75; Injected PV: 5.25 WAG ratio: 2.00; Injected PV: 6.00
400 WAG ratio: 0.50; Injected PV: 6.00 WAG ratio: 1.00; Injected PV: 6.00
WAG ratio: 0.25; Injected PV: 5.00 2 WAG ratio: 0.75; Injected PV: 5.25
Continuous Waterflooding; Injected PV: 3.00 WAG ratio: 0.50; Injected PV: 6.00
200 WAG ratio: 0.25; Injected PV: 5.00
Experiment

0 0
0 2000 4000 6000 8000 10000 0 2000 4000 6000 8000 10000
Time, min Time, min

Fig. 19—Production profiles with different water-slug sizes for Fig. 20—Cumulative WOR with different water-slug sizes for
Scenario 2. Scenario 2.

February 2013 SPE Reservoir Evaluation & Engineering 69


reservoirs. Two scenarios are designed to evaluate the effect of Flaaten, A.K., Nguyen, Q.P., Pope, G.A., et al. 2009. A Systematic Labo-
well configurations on pressure maintenance and oil recovery ratory Approach to Low-Cost, High-Performance Chemical Flooding.
with CO2-WAG injection. The well configuration with two paral- SPE Res Eval & Eng 12 (5): 713–723. http://dx.doi.org/10.2118/
lel horizontal wells (i.e., one is injector and the other is producer) 113469-PA.
yields a higher oil recovery with less injected fluids. After three Fletcher, P., Cobos, S., Jaska, C., et al. 2012. Improving Heavy Oil Recov-
CO2-WAG cycles, 37.6 and 32.5% oil recoveries are achieved for ery Using an Enhanced Polymer System. Paper SPE 154045 presented
Scenarios 1 and 2, respectively. Incremental oil recovery of 12.4 at SPE Improved Oil Recovery Symposium, Tulsa, Oklahoma, 14–18
and 8.9% through three CO2-WAG cycles was experimentally April. http://dx.doi.org/10.2118/154045-MS.
achieved for the aforementioned two well configurations, respec- Hadia, N. Chaudhari, L., Mitra, S.K., et al. 2007. Experimental Investiga-
tively. This is mainly because gasflooding leads to high micro- tion of Use of Horizontal Wells in Waterflooding. J. Pet. Sci. Eng. 56
scopic displacement efficiency, whereas waterflooding contributes (4): 303–310. http://dx.doi.org/10.1016/j.petrol.2006.10.004.
to large macroscopic displacement efficiency. Well configurations Huang, S.S., Pappas, E.S., and Jha, K.N. 1987. The Carbon Dioxide Im-
with horizontal wells lead to a better sweep efficiency. Theoreti- miscible Recovery Process and Its Potential for Saskatchewan Reser-
cally, simulation techniques have been developed to history match voirs. Paper SS-87-1 presented at the Technical Meeting/Petroleum
the experimental measurements. The excellent agreement between Conference of the South Saskatchewan Section, Regina, Saskatche-
the measured and simulated cumulative oil production indicates wan, Canada, 6–8 October. http://dx.doi.org/10.2118/SS-87-1.
that the displacement mechanisms governing water-alternating- Islam, M.R., Chakma, A., and Farouq Ali, S.M. 1989. State-of-the-Art of
CO2 processes have been numerically simulated and matched. It In-Situ Combustion Modeling and Operations. Paper SPE 18755 pre-
can be found from both the experimental measurements and nu- sented at SPE California Regional Meeting, Bakersfield, California,
merical simulation that well configuration has a dominant impact 5–7 April. http://dx.doi.org/10.2118/18755-MS.
on oil recovery. Displacement processes are sensitive to the rela- Issever, K., Pamir, A.N., and Tirek, A. 1993. Performance of a Heavy-Oil
tive permeability values at low water- and gas-saturation ranges. Field Under CO2 Injection, Bati Raman, Turkey. SPE Res Eng 8 (4):
Both experimental and simulation results indicate that the slug 256–260. http://dx.doi.org/10.2118/20883-PA.
size of water during CO2-WAG processes shows a positive effect Jha, K.N. 1986. A Laboratory Study of Heavy Oil Recovery With Carbon
on oil recovery. The well configuration associated with two hori- Dioxide. J. Cdn. Pet. Tech. 25 (2): 54–63. http://dx.doi.org/10.2118/
zontal wells tends to yield a higher oil recovery under the same 86-02-03.
WAG ratio. Regarding each well configuration, WAG ratio affects Li, H., Chen, S., Yang, D., et al. 2009. Estimation of Relative Permeability
oil recovery and cumulative WOR, and there exists an optimum by Assisted History Matching Using the Ensemble Kalman Filter
value for each scenario. The WAG ratios of 0.75 and 1.00 are the Method. J. Cdn. Pet. Tech. 51 (3): 205–214. http://dx.doi.org/10.2118/
optimum values for Scenarios 1 and 2, respectively. 156027-PA.
Li, H., and Yang, D. 2011. Estimation of Multiple Petrophysical Parame-
Acknowledgments ters for the PUNQ-S3 Model Using Ensemble-Based History Match-
The authors acknowledge a Discovery Grant from the Natural Sci- ing. Paper SPE 143583 presented at 2011 SPE EUROPEC/EAGE
ences and Engineering Research Council (NSERC) of Canada and Annual Conference and Exhibition, Vienna, Austria, 23–26 May.
an innovation fund from the Petroleum Technology Research http://dx.doi.org/10.2118/143583-MS.
Centre (PTRC) to D. Yang. Li, H., Zheng, S., and Yang, D. 2011. Enhanced Swelling Effect and Vis-
cosity Reduction of Solvents-CO2-Heavy Oil System. Paper SPE
150168 presented at the SPE Heavy Oil Conference and Exhibition,
References Kuwait City, Kuwait, 12–14 December. http://dx.doi.org/10.2118/
150168-MS.
Asghari, K., Dong, M., Shire, J., et al. 2007. Development of a Correlation
Between Performance of CO2 Flooding and the Past Performance of Mai, A., and Kantzas, A. 2009. Heavy Oil Waterflooding: Effects of Flow
Waterflooding in Weyburn Oil Field. SPE Prod & Oper 22 (2): Rate and Oil Viscosity. J. Cdn. Pet. Tech. 48 (3): 42–51. http://
260–264. http://dx.doi.org/10.2118/99789-PA. dx.doi.org/10.2118/09-03-42.
Beecher, C.E., and Parkhurst, I.P. 1926. Effect of Dissolved Gas Upon the Mangalsingh, D., and Jagai, T. 1996. A Laboratory Investigation of the
Viscosity and Surface Tension of Crude Oil. Pet. Trans. G-26 (1): Carbon Dioxide Immiscible Process. Paper SPE 36134 presented at
51–69. http://dx.doi.org/10.2118/926051-G. SPE Latin America/Caribbean Petroleum Engineering Conference,
Bowers, B., and Drummond, K.J. 1997. Conventional Crude Resources of Port-of-Spain, Trinidad, 23–26 April. http://dx.doi.org/10.2118/36134-
the Western Canada Sedimentary Basin. J. Cdn. Pet. Tech. 36 (2): MS.
56–63. http://dx.doi.org/10.2118/97-02-05. Moore, R.G., Laureshen, C.J., Ursenbach, M.G., et al. 1999. A Canadian
Chen, S., Li, H., and Yang, D. 2010a. Optimization of Production Per- Perspective on In Situ Combustion. J. Cdn. Pet. Tech. 38 (13): 1–8.
formance in a CO2 Flooding Reservoir under Uncertainty. J. Cdn. Pet. http://dx.doi.org/10.2118/99-13-35.
Tech. 49 (2): 71–78. http://dx.doi.org/10.2118/133374-PA. Nejad, K.S., Berg, E.A., and Ringen, J.K. 2011. Effect of Oil Viscosity on
Chen, S., Li, H., Yang, D., et al. 2010b. Optimal Parametric Design for Water/Oil Relative Permeability. Paper SCA 2011-12 presented at
Water-Alternating-Gas (WAG) Process in a CO2-Miscible Flooding International Symposium of the Society of Core Analysts, Austin,
Reservoir. J. Cdn. Pet. Tech. 49 (10): 75–82. http://dx.doi.org/ Texas, 18–21 September.
10.2118/141650-PA. Ning, S., Jhaveri, B., Jia, N., et al. 2011. Viscosity Reduction EOR With
Christensen, J.R., Stenby, E.H., and Skauge, A. 2001. Review of WAG CO2 and Enriched CO2 to Improve Recovery of Alaska North Slope
Field Experience. SPE Res Eval & Eng 4 (2): 97–106. http:// Viscous Oils. Paper SPE 144358 presented at SPE Western North
dx.doi.org/10.2118/71203-PA. American Region Meeting. Anchorage, Alaska, 7–11 May. http://
Cobanoglu, M. 2001. A Numerical Study to Evaluate the Use of WAG as dx.doi.org/10.2118/144358-MS.
an EOR Method for Oil Production Improvement at B. Kozluca Field, Peng, D.Y., and Robinson, D.B. 1976. A New Two-Constant Equation of
Turkey. Paper SPE 72127 presented at the SPE Asia Pacific Improved State. Ind. Eng. Chem. Fundam. 15 (1): 58–64. http://dx.doi.org/
Oil Recovery Conference, Kuala Lumpur, Malaysia, 6–9 October. 10.1021/i160057a011.
http://dx.doi.org/10.2118/72127-MS. Reid, T.B., and Robinson, H.J. 1981. Lick Creek Meakin Sand Unit Im-
Emera, M.K., and Sarma, H.K. 2008. A Genetic Algorithm-Based Model miscible CO2 Waterflood Project. J. Pet. Tech. 33 (9): 1723–1729.
to Predict CO2-Oil Physical Properties for Dead and Live Oil. J. Cdn. http://dx.doi.org/10.2118/9795-PA
Pet. Tech. 47 (2): 52–61. Righi, E.F., and Pascual, M. 2007. Water-Alternating-Gas Pilot in the
Eydinov, D., Gao, G., Li, G., et al. 2009. Simultaneous Estimation of Rel- Largest Oil Field in Argentina: Chihuido de la Sierra Negra, Neuquen
ative Permeability and Porosity/Permeability Fields by History Match- Basin. Paper SPE 108031 presented at Latin American and Caribbean
ing Production Data. J. Cdn. Pet. Tech. 48 (12): 13–25. http:// Petroleum Engineering Conference, Buenos Aires, Argentina, 15–18
dx.doi.org/10.2118/132159-PA. April. http://dx.doi.org/10.2118/108031-MS.

70 February 2013 SPE Reservoir Evaluation & Engineering


Rojas, G.A., and Farouq Ali, S.M. 1986. Scaled Model Studies of Carbon CO2 Flooding and WAG Processes on Recovery of Heavy Oils. Paper
Dioxide/Brine Injection Strategies for Heavy Oil Recovery From Thin CSUG 138188 presented at Canadian Unconventional Resources and
Formations. J. Cdn. Pet. Tech. 25 (1): 85–94. http://dx.doi.org/ International Petroleum Conference, Calgary, Alberta, Canada, 19–21
10.2118/86-01-07. October. http://dx.doi.org/10.2118/138188-MS.
Rojas, G.A., and Farouq Ali, S.M. 1988. Dynamics of Subcritical CO2/ Wang, J., Dong, M., and Asghari, K. 2006. Effect of Oil Viscosity on
Brine Floods for Heavy-Oil Recovery. SPE Res Eng. 3 (1): 35–44. Heavy Oil-Water Relative Permeability Curves. Paper SPE 99763 pre-
http://dx.doi.org/10.2118/13598-PA. sented at SPE/DOE Symposium on Improved Oil Recovery, Tulsa,
Sahin, S., Kalfa, U., and Celebioglu, D. 2008. Bati Raman Field Immisci- Oklahoma, 22–26 April. http://dx.doi.org/10.2118/99763-MS.
ble CO2 Application-Status Quo and Future Plans. SPE Res Eval & Watson, A.T., Seinfeld, J.H., Gavalas, G.R., et al. 1980. History Matching
Eng 11 (4): 778–791. http://dx.doi.org/10.2118/106575-PA. in Two-Phase Petroleum Reservoirs. SPE J. 20 (6): 521–532. http://
Shen, C. 2002. Limitations and Potentials of In-Situ Combustion Processes dx.doi.org/10.2118/8250-PA.
for Heavy Oil Reservoirs. Paper SPE 2002-217 presented at Canadian Zheng, S., Li, H., and Yang, D. 2011. Pressure Maintenance and Improv-
International Petroleum Conference, Calgary, Alberta, Canada, 11–13 ing Oil Recovery with CO2 Injection in Heavy Oil Reservoirs. Paper
June. http://dx.doi.org/10.2118/2002-217. SPE 150169 presented at the SPE Heavy Oil Conference and Exhibi-
Shi, W., Corwith, J., Bouchard, A., et al. 2008. Kuparuk MWAG Project tion, Kuwait City, Kuwait, 12–14 December. http://dx.doi.org/10.2118/
After 20 Years. Paper SPE 113933 presented at SPE/DOE Symposium 150169-MS.
on Improved Oil Recovery, Tulsa, Oklahoma, 20–23 April. http://
dx.doi.org/10.2118/113933-MS. Sixu Zheng is currently a PhD candidate in petroleum systems
Smith, G.E. 1988. Fluid Flow and Sand Production in Heavy-Oil Reser- engineering in the Faculty of Engineering and Applied Sci-
voirs Under Solution-Gas Drive. SPE Prod Eng 3 (2): 169–180. http:// ence at the University of Regina, Canada. His major research
dx.doi.org/10.2118/15094-PA. interests include CO2 EOR and storage, reservoir modeling
Sohrabi, M., Henderson, G.D., Tehrani, D.H., et al. 2000. Visualisation of and simulation, hybrid steam/solvent processes for heavy-oil
Oil Recovery by Water Alternating Gas (WAG) Injection Using High recovery, and production optimization. Zheng holds a BSc
Pressure Micromodels-Water-Wet System. Paper SPE 63000 presented degree in petroleum engineering from the China University of
at SPE Annual Technical Conference and Exhibition, Dallas, Texas, Petroleum and an MASc degree in petroleum systems engi-
neering from the University of Regina. He is a member of SPE.
1–4 October. http://dx.doi.org/10.2118/63000-MS.
Spiteri, E.J., and Juanes, R. 2004. Impact of Relative Permeability Hyster- Daoyong Yang is a professor in and the Chair of the Petroleum
esis on the Numerical Simulation of WAG Injection. Paper SPE 89921 Systems Engineering Program in the Faculty of Engineering
presented at the SPE Annual Technical Conference and Exhibition, and Applied Science at the University of Regina. He holds BSc
Houston, Texas, 26–29 September. http://dx.doi.org/10.2118/89921- and a PhD degrees in petroleum engineering from the China
University of Petroleum and a PhD degree in petroleum sys-
MS.
tems engineering from the University of Regina. Previously,
Stoll, W.M., al Shureqi, H., Finol, J., et al. 2010. Alkaline-Surfactant- Yang worked for 3 years as a petroleum engineer in the Petro-
Polymer Flood: From the Laboratory to the Field. Paper SPE 129164 China TuHa Oilfield Company (CNPC). He also worked 4 years
presented at SPE EOR Conference at Oil & Gas West Asia, Muscat, as a reservoir engineer in the Software Development and In-
Oman, 11–13 April. http://dx.doi.org/10.2118/129164-MS. formation Centre of Petroleum Engineering, CNPC. Yang’s
Strauss, J.P., Alexander, D.M., Al-azri, N., et al. 2010. EOR Development major research areas include CO2 EOR and storage, formation
Screening of a Heterogeneous Heavy Oil Field-Challenges and Solu- evaluation, reservoir modeling and simulation, phase behav-
tions. Paper SPE 129157 presented at the SPE EOR Conference at Oil ior, production optimization, artificial-lift methods, interfacial
& Gas West Asia, Muscat, Oman, 11–13 April. http://dx.doi.org/ interactions in EOR processes, and heavy-oil recovery. He has
authored or coauthored more than 80 refereed-journal
10.2118/129157-MS.
papers and conference papers. Yang was the recipient of the
Thomas, S., Farouq Ali, S.M., Scoular, J.R., et al. 2001. Chemical Meth- SPE Canada Region Regional Formation Evaluation Award in
ods for Heavy Oil Recovery. J. Cdn. Pet. Tech. 40 (3): 56–61. http:// 2011. He is an active member of SPE, American Chemical So-
dx.doi.org/10.2118/01-03-05. ciety, and Association of Professional Engineers and Geo-
Torabi, F., Jamaloei, B.Y., Zarivnyy, O., et al. 2010. Effect of Oil Viscos- scientists of Saskatchewan and is a registered professional
ity, Permeability and Injection Rate on Performance of Waterflooding, engineer in Saskatchewan, Canada.

February 2013 SPE Reservoir Evaluation & Engineering 71

You might also like