Performance Analysis of Wells With Downhole Water Loop Installation For Water Coning Control
Performance Analysis of Wells With Downhole Water Loop Installation For Water Coning Control
0.8 6000
0.7
5000
Water Cut
0.6
Pwftop, psi
0.5 4000
Qtop = 200 bpd
0.4
Qtop = 400 bpd
0.3 3000
Qtop = 2000 bpd
0.2 Qtop = 3000 bpd Qtop = 200 bpd Qtop = 2000 bpd
2000
Qtop = 5000 bpd Qtop = 400 bpd Qtop = 3000 bpd
0.1
Qtop = 5000 bpd
0.0 1000
0 200 400 600 800 1000 1200 0 100 200 300 400 500 600
FIGURE 4: Water coning stabilization process with equilibrium water FIGURE 5: Stabilization of bottomhole flowing pressure.
cut.
show that the model may reach equilibrium condition after a con-
roughly in geometric progression. There are 10 layers in both the siderably long time period. We believe that the model properly
oil zone and aquifer in a vertical direction; the grid size is equally describes pressure interference between the three completions –
distributed in this direction. To simulate steady-state flow condi- the effect that entirely controls DWL system performance.
tions, the produced fluids from both completions are reinjected at
the outermost gridblock into the oil and water zones, respectively The reservoir and well data and fluid properties used in this
[Inikori(22)]. The reservoir model would maintain constant volume study are shown in Tables 1, 2 and 3 [Arslan(24)]. The data reveals
of fluid in place andv constant reservoir pressure, thus simulating a combination of high porosity and low-permeability – the prop-
a stabilized well production. The aquifer model is a coaxial cyl- erties of an oilfield in China [Song et al.(25)] where the porosity in
inder of large size connected to the bottom of the oil pay zone. the reservoir is as high as 0.3, while the permeability is only 80
After an initial transient period of developing water cone, the mD. High porosity with low-permeability is caused by the high
system attains steady state [Inikori(22), Arslan et al.(23), Arslan(24) clay content in the sandstone. Hu and Dietrich(26) also reported a
and Hernandez(1)]. The duration of the transient period depends similar field case in Canada. Neasham(27) gave a concrete analysis
on system properties (permeability, size, etc.) and the set limit for on the dispersed clay in sandstone reservoirs and its effect on fluid-
flowing pressure stability (convergence tolerance). Figures 4 and 5 flow properties. It is not uncommon to have reservoirs with high
U.S. Conversion SI
Data Field Unit Water Oil Gas Rock Factor Unit
Pref psi 6000 6000 6000 6000 × 6.894 757 E+00 kpa
B rb/stb 1.02 1.26 × 1.000 E+00 rm3/sm3
C 1/psi 0.000003 0.000015 0.000004 × 1.450 377 E-01 1/kpa
µ cp 0.4 0.8 × 1.000 E−03 Pa·s
ρsc lbm/ ft3 68 48 0.0702 × 1.601 846 E+01 kg/m3
Pb psi 1000 × 6.894 757 E+00 kpa
Rs at Pb Mscf/stb 0.379 × 1.801 175 2 E+0 m3/m3
porosity and low-permeability. For example, this effect is widely Table 3: Relative permeability and capillary pressure.
spread across the Norphlet formation in the southern coast of the
U.S. [Dixon et al.(28)]. The reason for choosing the data in this Case Sw krw kro Pc (psi) Pc (kpa)
study is because water coning is more severe in these kinds of res- 0.20 0.0000 1.0000 1.428 9.846
ervoirs [Joshi(29)]. However, we believe the choice of data should 0.30 0.0002 0.8374 1.250 8.618
not significantly affect findings of the study. 0.40 0.0039 0.6328 1.071 7.384
A 0.55 0.0366 0.3337 0.803 5.536
Approach 0.70 0.1526 0.1187 0.536 3.696
0.80 0.3164 0.0391 0.357 2.461
Our main objective was to develop an understanding of the ef- 1.00 1.0000 0.0000 0.000 0.000
fects of DWL well completion variables on well system perfor-
mance. Hence, three operational parameters have been selected for 0.20 0.0000 1.0000 1.428 9.846
the simulation study: top (production) rate, drainage-injection rate 0.30 0.0029 0.6944 1.250 8.618
(bottom rate) and the D/I spacing. In order to see individual effects, B 0.40 0.0128 0.4444 1.071 7.384
we would set the value of one parameter and vary the other two pa- 0.55 0.0605 0.1736 0.803 5.536
0.70 0.1898 0.0278 0.536 3.696
rameters to observe their relative effects. During each simulation
0.80 0.3547 0.0000 0.357 2.461
run, we would bring the system to equilibrium by stabilizing the 1.00 1.0000 0.0000 0.000 0.000
flow so both the water cut (WC) and flowing pressure (pwf) would
become constant. We also assumed zero skin damage and a com-
plete integrity of the cement sheath [Arslan(24)]. damaged. (Even minimal oil content in injection water would de-
With the high-speed reservoir simulator, Eclipse 2007TM, which posit a residual oil-saturated skin zone around the injection com-
was used in the study, a typical simulation run would be as follows: pletion, thus reducing permeability to water and injectivity of the
fix the drainage rate and change the top rate and distance. (For ex- completion.) Therefore, it is important to avoid oil in the injection
ample, for drainage rate 3200 B/D, we would change the top rate water [Pang and Sharma(30) and Abou-Sayed et al.(3)].
to 480 B/D, 640 B/D and 800 B/D, with D/I spacings for each top Let’s consider the oil content limit in the injection water. Figure
rate at 0 ft, 50 ft, 100 ft, 200 ft, 300 ft, 400 ft and 480 ft. The total 9 shows that oil content in injection water changes with D/I spacing
number of runs for this scenario is 1×3×7=21.) Typical simulation and injection rate. When the top rate is high and the bottom rate is
scenarios used in this study are shown in Table 4. low or the D/I spacing is short, there will be little oil in the injec-
tion water. These results correspond to the drainage completion
located 10 ft below the OWC. Also note that if the drainage com-
pletion is placed 20 ft below the OWC, there will be no oil in the
Results and Observations injection water.
In conventional wells, operated at a constant liquid produc- To produce oil at a higher rate, while avoiding oil in injec-
tion rate, the increased water cut leads to less oil produced – a tion water, an Inflow Performance Window (IPW) can be de-
common problem in production practices. In contrast, all simula- veloped for different combinations of the top and bottom rates,
tions of DWL wells consistently show lower water cut and higher as shown in Figure 10. In the figure, there are three plots of Qbot
oil rates, as shown in Figures 6 and 7. The results can be summa- vs. Qtop, each of them being an “envelope” representing different
rized as follows: D/I spacing: 10 ft, 100 ft and 480 ft. The top boundary of each
• Larger spacing between water drainage and injection com- envelope represents the maximum oil-free water drainage rate,
pletions (D/I spacing) dramatically improves well perfor- while the lower boundary corresponds to the maximum top rate
mance by reducing water cut and increasing oil rate. of oil with no water. (Points inside the envelopes correspond to a
• There is an optimum value of D/I spacing (≈200 ft in this
case), above which there is no further improvement. Table 4: Simulation scenarios.
• Critical rate increases with increased D/I spacing.
Bottom Rate Top Rate D/I Spacing
• At optimum D/I spacing, the effect of water drainage on well
(bpd) (bpd) (ft)
performance is strong; a 60% and 100% increase of water
drainage would increase oil production by 55% and 80%, 480 0,50,100,200,300,400,480
respectively. 800 640 0,50,100,200,300,400,480
• For each DWL system, there is a combination of the top 800 0,50,100,200,300,400,480
(production) rate, bottom (drainage-injection) rate and D/I
spacing that would result in water-free oil production. The 480 0,50,100,200,300,400,480
observation is consistent with previous studies of DWS sys- 2000 640 0,50,100,200,300,400,480
800 0,50,100,200,300,400,480
tems [Shirman(2), Siemek and Stopa(11), Arslan(24), Johns et
al.(12), Hernandez(1) and Utama(13)].
480 0,50,100,200,300,400,480
• For a given D/I spacing, an effective increase of oil pro- 640 0,50,100,200,300,400,480
3200
duction requires synchronized increases in production and 800 0,50,100,200,300,400,480
drainage rates, as shown in Figure 8. A sole increasing of the
top production rate is not effective. 480 0,50,100,200,300,400,480
Oil-free water looping is an important condition for DWL. If 4000 640 0,50,100,200,300,400,480
the drainage-injection water contains some oil, not only will the 800 0,50,100,200,300,400,480
aquifer be polluted, but also the injection completion will be
June 2010, Volume 49, No. 6 41
1.0 1.0
1.0 Qtop = 480 bpd Qtop = 480 bpd Qtop = 480 bpd
0.9 0.9
0.9 Qtop = 640 bpd Qtop = 640 bpd Qtop = 640 bpd
0.8 0.8
0.8 Qtop = 800 bpd Qtop = 800 bpd Qtop = 800 bpd
0.7 0.7
0.7 (Qbot = 2000 bpd) (Qbot = 3200 bpd) (Qbot = 4000 bpd)
Water Cut
Water Cut
Water Cut
0.6 0.6
0.6
0.5 0.5
0.5
0.4 0.4
0.4
0.3 0.3
0.3
0.2 0.2
0.2
0.1 0.1
0.1
0.0 0.0
0.0
0 100 200 300 400 500 600 0 100 200 300 400 500 600 0 100 200 300 400 500 600
FIGURE 6: Water cut changes with top rate and D/I spacing at various bottom rate.
combination of the rates that the oil is water-free and drainage water A typical nodal analysis involves solving the system of two
contains no oil – i.e., segregated production.) It is evident that with the equations, TPR and IPR, for the optimum values of flowing
injection point approaching the water drainage completion, the pressure and production rate. Golan and Whitson(34) and Beggs(35)
performance window of the DWL system moves into the area of presented a wide selection of available calculation models and ap-
larger water rates and smaller oil production. plied the nodal analysis method to production oil and gas wells, as
It is important to see that the DWL water re-injection tech- well as modifying the method for injection wells.
nology would work even for a close distance between drainage Arslan et al.(23) and Arslan(24) applied the nodal analysis method
and injection completions (small D/I spacing). This means that the
to DWS wells. They pointed out that oil wells would be generally
DWL system does not require deep drilling or deviated wells to in-
produced at their highest possible rate to maximize the cash flow.
ject produced water. Moreover, the system would not reduce the
water-drive ability of the aquifer (reservoir pressure) while giving However, the maximum allowable pressure drawdown at the two
a low-water-cut oil production. From the previously described re- completions of the DWS well was a limiting factor caused by prac-
sults we can see that there exists a critical (minimum) D/I spacing tical considerations, such as well integrity, sand control, gas lib-
in the DWL system. At this distance, the oil rate is at the maximum eration, etc. The nodal analysis for DWS wells was based on the
and there is no oil in the injection water. following operational principles:
• Produce at maximum possible top rate (economic goal).
• Keep the pressure drawdown at or below the maximum al-
System (Nodal) Analysis of DWL Well lowable pressure drawdown for both completions (comple-
tion limit).
Nodal analysis is a well-known method for assessing perfor-
mance of petroleum wells as production systems – from the outer 1.0
(Qbot = 2000 bpd)
boundary of the reservoir to the sand face, across the perforations 0.9
and completion section to the tubing intake, up the tubing string 0.8
with its restrictions and down-hole safety valves, past the surface 0.7
choke, the flow lines and separators. The name “nodal” was pro-
Water Cut
0.6 D/I = 50 ft
posed by Gilbert(31), who proposed considering an oil well as a se-
0.5 D/I = 100 ft
ries of hydraulic components bracketed by appropriately placed
0.4 D/I = 200 ft
nodes. Later, Brown and Beggs(32) and Mach et al.(33) systemati- D/I = 300 ft
cally studied Gilbert’s idea and further improved the system anal- 0.3
D/I = 400 ft
ysis method. They clearly defined the concept by reducing the 0.2
D/I = 480 ft
complex production system to three basic correlations describing: 0.1
400
300
Top Oil Rate, bpd
350
250
300
200 250
150 200
150
100
D/I = 10 ft D/I = 200 ft D/I = 10 ft D/I = 200 ft
100
D/I = 300 ft D/I = 480 ft D/I = 480 ft
50 D/I = 50 ft D/I = 50 ft D/I = 300 ft
50
D/I = 100 ft D/I = 400 ft D/I = 100 ft D/I = 400 ft
0 0
0 500 1000 1500 2000 0 500 1000 1500 2000
FIGURE 8: Top completion oil rate changes with top completion liquid rate at various D/I spacing and bottom rate.
42 Journal of Canadian Petroleum Technology
2500 2500 2500
0 0 0
0 100 200 300 400 500 600 0 100 200 300 400 500 600 0 100 200 300 400 500 600
D/I Spacing, feet D/I Spacing, feet D/I Spacing, feet
FIGURE 9: Oil content in injection water depends upon D/I spacing and the combination of production and drainage rates.
• Control the water drainage rate, such that the operational the DWS top completion. A considerable improvement in water
point (top and bottom rates) is below the “flip-flop” line to control and well productivity could be achieved by using DWS in
avoid reverse coning (environmental limit). a heavy oil reservoir.
• Select submersible pumps’ rating, such that that the both Because of the similarity between DWS and DWL wells, the
completions inflow could be lifted to surface (TPR limit). DWS nodal analysis method could be extended to DWL wells by
As the nodal analysis only gives well performance at a given evaluating the pressure vs. flow rate relationship for the top com-
time, the calculated production rates change with time. For a se- pletion for a range of bottom rates and D/I spacings. As an addi-
ries of timesteps, the model was used to predict the DWS well’s tional constraint, oil should be eliminated from the water drainage
recovery over the well’s life. At each timestep, oil depletion (incre- so the injected water is oil-free. This constraint is added to the
mental recovery) was estimated and the new oil pay thickness was DWL well performance plots, shown in Figures 11 and 12 for
updated from material balance. Future rates were predicted succes- drainage-injection rates 4000 B/D and 7200 BWPD, respectively.
sively in this manner. By using this method, it was concluded that In the plots, several upward lines represent TPR with variable
DWS technology could increase oil production rates more than water cut – each line for a different D/I spacing. The numbered in-
two-fold when compared with conventional wells. tercepts (1, 2, etc.) indicate the rate/pressure conditions of natural
Qin and Wojtanowicz(36) applied the DWS nodal analysis flow at the top completion. It is clear that DWL increases produc-
method to heavy oil reservoirs with bottomwater. They considered tion rate and reduces water cut compared to a conventional well.
a variable rate-dependent water cut at stabilized well production The results also show that larger D/I spacing improves DWL per-
and various ways to reduce the water cut with DWS. They found formance until D/I optimum value is reached.
out that the operational range of production rates with variable In Figures 11 and 12, lines AB and A’B’ represent the limit of
water cut was limited for heavy oil. Substantial value of pressure oil-free water drainage. For each D/I spacing and drainage rate,
drawdown needed for heavy oil inflow was diminished after the there is a minimum production rate (Qtop) needed to prevent oil
water breakthrough. Therefore, DWS bottomwater drainage was drainage. Figure 11 shows a case when the oil-free drainage limit
instrumental in stimulating inflow of heavy oil without water to is met for all values of D/I spacing. When the bottom rate (Qbot)
is 4000 B/D and the oil-free injection limitation line AB is below
2000 the IPR line, this means that reservoir energy is enough to sustain
1800 natural flow rates. In this case (the 850 B/D rate with 0.75 water
1600 cut), the DWL well’s production rate is 1050 B/D with 0.56 water
1400
cut for the 50 ft D/I spacing. Moreover, the DWL production rate
would increase to 1280 B/D with water cut reduced to 0.31 for the
1200
480 ft D/I spacing. For the two values of D/I spacing, the oil rates
Qbot, bpd
TPRDWLD/I = 50 ft
Pwf, psi
TPRDWLD/I = 50 ft
4500 1 – WC = 0.75 TPRDWLD/I = 100 ft
4500 1 – WC = 0.75 TPRDWLD/I = 100 ft
2 – WC = 0.50
2 – WC = 0.60 TPRDWLD/I = 200 ft
TPRDWLD/I = 200 ft 3 – WC = 0.43
3 – WC = 0.56 4000
4000 4 – WC = 0.36 TPRDWLD/I = 300 ft
4 – WC = 0.47 TPRDWLD/I = 300 ft
5 – WC = 0.40 TPRDWLD/I = 400 ft