Reservoir Saturation Tool
Reservoir Saturation Tool
Reservoir Saturation Tool
Chapter overview
Saturation through casing
Porosity through casing
Lithology through casing
Permeability and pressure
Advances in production logging
Hydraulic integrity
Permanent monitoring systems:
Extending the life of a reservoir
E VA L U AT I O N A N D M O N I T O R I N G O F E X I S T I N G W E L L S
CHAPTER OVERVIEW
S AT U R AT I O N T H R O U G H C A S I N G
E VA L U AT I O N A N D M O N I T O R I N G O F E X I S T I N G W E L L S
(3)
Charts (Schlumberger Log Interpretation
Charts, 1997) are available to determine Sw
and Sh. Sw is directly related to the salinity
of the water, its value ranging from 20 capture units (c.u.) for fresh water to more
than 120 c.u. for salt-saturated water. Sh
typically varies between 18 c.u. and 22 c.u.
for oil and remains below 17 c.u. for gas.
The quantity f(Sw-Sh) in the denominator
of equation (3) shows that PNC logging
requires a contrast between Sw and Sh. Hence,
the water must be salty and the salinity must
be known to evaluate Sw. There are no hard
limits, but it is generally considered that when
f is less than 10 p.u. and water salinity less
than 30,000 ppm, the contrast is insufficient.
These conditions are usually not met in
reservoirs where production is maintained by
water injection programs. The resulting formation water, a mixture of injection water and
connate water, has a variable and often low
salinity. Also, most Venezuelan reservoirs have
connate water salinities outside this range.
6 3
YC
K + K f(1-Sw) + KC3(1-yw)
= C1 C2
KOX1 + KOX2fSw + KOX3yw
YOX
(4)
S AT U R AT I O N T H R O U G H C A S I N G
Figure 6.1
Counts
Hydrogen
Oil
Carbon
Oxygen
Water
Energy (Me V)
RST spectra obtained in an oil tank and a water tank. The broad energy windows
used across the region of the predominant C and O peaks are indicated.
Figure 6.2
VUOI-COR
ELAN effective porosity
Depth 10
(ft)
(p.u.)
VUOI-WIN
40
20
VUOI-ALFA
(p.u.)
70
3650
3675
3700
3725
Alpha processing results in an accurate volume of oil with low statistical variations.
Alpha processing
Spectral processing is hampered by
large statistical uncertainties and requires
very low logging speeds to ensure that
enough gamma rays are measured.
Another approach uses broad energy
windows in the gamma ray spectrum across
the region of the predominant C or O peaks
(Fig. 6.1). This method gives better statistical
precision but the interpretation is prone to
bias. The windows do not only represent C
and O, but are strongly influenced by other
elements present in the formation and by the
gamma radiation induced by the tool itself.
Consequently, the ratio from the C and O
windows is very dependent on environmental
effects, especially borehole effects such as
washouts and cement integrity.
To improve the statistical precision without losing accuracy, a method called Alpha
Processing combines the windows and spectral processing to obtain an accuracy comparable to the C/O yields ratio and a precision comparable to the C/O windows ratio.
Figure 6.2 shows for one well the
formation porosity with the volume of oil
computed from the yields (VUOI-COR) and
windows (VUOI-WIN) ratios, and the same
volume resulting from alpha processing
(VUOI-ALFA). VUOI-WIN exhibits substantially fewer statistical variations than
6
E VA L U AT I O N A N D M O N I T O R I N G O F E X I S T I N G W E L L S
Oil in borehole
rh = 0.74 g cm3
Gravel pack f = 35%
Sw (gravel pack) = Sw (for
.7
.6
Sw = 0 %
.5
.4
.3
Sw = 100 %
.2
.1
.0
0
.1
Porosity
(a)
S AT U R AT I O N T H R O U G H C A S I N G
Figure 6.4
Oil in borehole
rh = 0.74 g cm3
Gravel pack f = 35%
Sw (gravel pack) = Sw (formation)
.7
.6
Sw = 0 %
.5
.4
.3
Sw = 100 %
.2
.1
.6
Sw = 0 %
.5
.4
.3
Sw = 100 %
.2
.1
.0
.0
0
.1
.2
Porosity
(a)
.3
.1
.2
.3
Porosity
(b)
E VA L U AT I O N A N D M O N I T O R I N G O F E X I S T I N G W E L L S
Figure 6.5
Sigma borehole
Depth
(ft)
250
(c.u.)
GR
(gAPI)
Sigma formation
0
100 .2
Deep induction
(ohm-m) 2000
Near COR
60
(c.u.)
.1
60
RST porosity
(p.u.)
.1
.4
Far COR
.4
3500
3600
3700
S AT U R AT I O N T H R O U G H C A S I N G
Figure 6.6
Water
Moved oil
Sw RST
100 (p.u.)
Sw OH
RST porosity
50 (p.u.)
0 100 (p.u.)
Depth
(ft)
Oil
Water
Water
Oil
Oil
0 Fluid analysis
from RST
0 50 (p.u.)
Sand
Bound water
Fluid analysis
from OH logs
0 50
(p.u.)
Clay
0 0
Volume analysis
(p.u.)
100
8
7
6
3500
3600
3700
Figure 6.7
1400
Fluid conduit
Theoretical performance
Predicted performance
assuming 30% water cut
Pressure (psig)
1200
1000
800
600
400
200
0
0
100
200
300
400
,
,,
@@
A
B
B
,
,
,,
,
@
A
C
BB
,
@
,
,
A
,,
@@
,
,
BB
C
A
,,
@@
,
BB
C
A
,,
@@
,
,
,
,
,,
,
@
A
BB
,
@
,
,
E VA L U AT I O N A N D M O N I T O R I N G O F E X I S T I N G W E L L S
Figure 6.8
Water
Moved oil
Water
Perfs.
Depth
(ft)
Sw RST
100 (p.u.)
Sw OH
Permeability
0
(md) 200 100 (p.u.)
Oil
0
0 25
(p.u.)
2
15,300
15,500
9
15,600
10
11
12
13
15,700
Sand
Oil
Fluid analysis
from RST
15,400
Oil
Water
Bound water
Fluid analysis
from OH logs
0 25
(p.u.)
0 0
Clay
Volume analysis
(p.u.)
100
Conclusion
A crucial piece of information for
reservoir monitoring is water saturation.
By combining and improving two techniques
to obtain this information, the RST is a breakthrough in cased-hole logging. Sigma logging
can quickly provide Sw in formations
with known and high water salinity.
C/O logging is the alternative in reservoirs
with low formation water salinity or under
water injection.
AB@,,,B@,
CBB
RST
S AT U R AT I O N T H R O U G H C A S I N G
RST
Figure 6.9
Figure
6.10a
Inelastic scattering
Nucleus
Fast
neutron
6.10b
Thermalization
Fast
neutron
Nucleus
6.10c
Neutron capture
Slow
neutron
Nucleus
6.10d
Excited
nucleus
Time (msec)
Fig. 6.11
100
80
C
60
Net inelastic =
A bB
40
Net inelastic
B
20
A
Counts
Counts
Neutron
generator
Counts
Far
Near
detector detector
Counts
Acquisition
cartridge
Neutron
burst
Telemetry,
casing collar
locator
(CCL)
Gamma ray
(optional)
Neutron
burst
Energy
0
Energy
g-ray
Activation
n + 16O 16N + p 16O* + b 16O + g with half life = 7.1s
6 10
E VA L U AT I O N A N D M O N I T O R I N G O F E X I S T I N G W E L L S
RST
Fig.
6.12
Counts
Oxygen
Silicon
Tool background
Calcium
Iron
Carbon
Energy (MeV)
Flowing
signal
Stationary
signal
Background
signal
Fast flowing signal
Count rate
Time
Long burst
Fig. 6.14
Counts
Short burst
500
1000
Time (msec)
Fig. 6.13
6 11
1500
Various methods are used to obtain porosity in old wells. If an adequate modern openhole logging suite is available, it is usually
sufficient to rely on the openhole porosity.
However, in many cases the openhole logging
suite is insufficient or non-existent. In other
cases, the openhole porosity is no longer accurate; for example, carbonates that have been
acidized. Nuclear and acoustic measurements
must then be used to determine the porosity.
The most commonly used technique is a
neutron-type porosity measurement. Such
measurements are relatively insensitive to the
casing and the cement, and can be characterized so as to remove these effects. The
most convenient is a through-tubing measurement, such as from the 11116-in. or 212-in. RST .
The
next
section
describes
the
characterization and processing of this
measurement and shows an example in a
complicated completion. Density tools are
also used. However, in general, the effects of
casing and cement are significant, and unless
the cement bond is excellent, the results are
considered unreliable.
Figure 6.15
Borehole
salinity
0 kppm
Count rate
Count rate
Borehole
salinity
0 kppm
200 kppm
200 kppm
Long
burst
Short
burst
50
100
150
200
600
1000
1400
Time (ms)
Time (ms)
a.) Typical signal decay from the short burst of the RST near
b.) Typical signal decay from the long burst of the RST far
6 12
E VA L U AT I O N A N D M O N I T O R I N G O F E X I S T I N G W E L L S
Figure 6.16
Input
Time decay spectra
Correction to spectra
Counting loss corrections
Background adaptive filtering
Background subtraction
SBNA
SFFA
TRAT
etc.
Environmental
parameters
Database
Borehole size
Casing size weight
Lithology
Outputs
Borehole salinity Cross section
Porosity
Corrected near and far sigma
Sigma formation
BSAL SIBF
TPHI
SFNC SFFC
SIGM
Openhole
6.4
2.8755
3.5
12.7
Porosity2
Formation
flush3
(kppm NaCl)
Borehole
fluid4
(kppm NaCl)
LS
LS
LS
Z*M H
Z*M H
Z*M H
0 70 140 210
0 70 140 210
0 70 140 210
6
6
6
4.5
5.0
Openhole
10.5
18
LSD
LSD
LSD
Z*M H
Z*M H
Z*M H
0 70 140 210
0 70 140 210
0 70 140 210
8
8
8
5.5
7.0
Openhole
15.5
32
LSD
LSD
LSD
Z*M H
Z*M H
Z*M H
0 70 140 210
0 70 140 210
0 70 140 210
10
10
10
5.5
7.0
7.625
15.5
32
26.4
LSD
LSD
LSD
Z*M H
Z*M H
Z*M H
0 70 140 210
0 70 140 210
0 70 140 210
11
12
12
7.625
9.625
9.625
26.4
32.3
53.5
LS
LS
LS
Z*M H
Z*M H
Z*M H
0 70 140 210
0 70 140 210
0 70 140 210
6 13
Figure 6.17
60
50
40
30
20
10
-1.5 -0.5 0.5 1.5
Deviation from assigned sigma (c.u.)
0
0
10
20
30
40
50
60
Figure 6.18
Water
Gas
Depth
(ft)
CALI
30 (in.)
CALI
0 (in.) 30
50
TPHI
(p.u.)
50
NPHI
(p.u.)
1.85
RHOB
(g / cm3)
2.65
50
Moved Oil
Oil
Water
Sand
Gas
Bound water
Oil
Clay
Fluid analysis
(p.u.)
Volumetric analysis
0
(p.u.)
100
5800
B
C
5900
E
6400
Example of porosity and C/O measurements from an RST run inside a 3 1/2-in.
processed (see page 64) to improve precision, using SFFC as the accurate measurement and SFFA as the more precise one.
The accuracy of the processing has been
checked against the database itself. Figure
6.17 shows the measured sigma versus. the
assigned sigma for 900 points where Sf < Sb,
which is the normal situation. The average
absolute errors for the RST-A tool are 0.22 c.u.
for SIGM, 0.9 p.u. for TPHI, and 5.5 kppm for
BSAL. For the RST-B tool the corresponding
errors are 0.2 for SIGM, 1.7 p.u. for TPHI and
6.9 kppm for BSAL. The overall accuracy of
the measurement has also been checked
against the industry standard calibration
pits in Houston and in the EUROPA facility
in Aberdeen.
Figure 6.18 shows an example of a well
logged from inside a 312-in. tubing lying within a 958-in. casing. The well produces from a
lower formation and has a packer set at
6963 ft, approximately 500 ft below the
interval logged. Since the well is deviated 20,
the tubing is expected to lie against the side
of the casing. There is a gas lift mandrel in the
interval logged at 5933 ft, which could
possibly cause some stand-off. However, this
is opposite a shale and the effects on the RST
outputs are not significant. There was water
in the tubing and the tubing-casing annulus at
the time of logging.
Such a dual-string completion had not
been included in the database at the time of
logging. However, it could be approximated
by taking the total thickness of both casings
as the casing thickness. The results show
that this is a good approximation, with good
agreement between the openhole neutron
(NPHI) and the RST porosity (TPHI).
Track 4 shows the evaluation of the openhole logs. Track 3 shows the openhole fluid
analysis compared to the fluid analysis from
the RST C/O measurement. As in Figs. 6.6 and
6.8, moved oil is the difference between the oil
seen by the RST and the openhole logs. Track
2 shows NPHI, density and TPHI.
6 14
E VA L U AT I O N A N D M O N I T O R I N G O F E X I S T I N G W E L L S
Figure 6.19
GR
0
(gAPI)
S / T Projection
Depth
150
(ft)
40
(ms / ft)
Sonic Waveform
240
500
(ms )
1700
S / T line
150
S / T Plane : P & S (Rec.)
3660
200
15,000
175
150
125
100
75
50
500
Time (ms )
S / T line
150
S / T Plane : P & S (Rec.)
3660
15,100
175
150
125
100
75
50
500
Time (ms )
S / T line
150
S / T Plane : P & S (Rec.)
3660
15,200
225
200
175
150
125
100
75
50
500
Time (ms )
15,300
6 15
Figure 6.20
CALI
(in)
Vp Vs
SP
20
(mV)
GR
(gAPI)
0.1
120
150
16
Depth
(ft)
LLS
(ohm-m)
1000
1.825
LLD
0.1
(ohm-m)
1000
113
RHOB
(g / cm3)
DTCO
(ms / ft)
2.65
55
15,000
15,100
15,200
15,300
Zone D (5900 ft) was found to be gasfilled when the well was drilled. However,
TPHI shows no gas effect, and the COR
indicates oil. This is an interesting indication
of the movement of fluids in the reservoir.
Zone E (6375 ft) shows no change in
fluid content between the openhole logs
and the RST log. NPHI and TPHI overlay
almost exactly, as they do in the water zone
below, except opposite the shalier intervals
where NPHI reads higher.
Overall, this example confirms that the
RST porosity gives good results through two
strings in 1214-in. hole in these conditions of
25 p.u. sands and eccentered tubing. It also
shows the benefit of monitoring old wells.
The comparison of RST porosity and saturation with the original openhole results not
Acoustic measurements
through casing
First arrivals have been used for many
years to determine compressional slowness
through casing. However, success depends
on the cementation being good enough to
allow the generation and transmission of
strong formation signals. Also, unless there is
a near-perfect cement bond, the formation
slowness and casing slowness must be
significantly different.
These limitations can be overcome by
recording and processing full waveforms.
Two special techniques are recommended to
improve the results. First, the data should be
acquired at a lower transmitter frequency
than in openhole. Second, the processing
should include a wave separation technique
to remove the casing signal.
The amplitude of the compressional signal
increases with frequency because of the effects
of borehole resonance. However, in cased
holes more casing propagation modes are
excited at higher frequency, especially when
the cement bond is poor. Therefore, the frequency must be lowered, but without reducing
the amplitude of the compressional signal below measurable limits. This limit depends on
formation slowness, since the compressional
amplitude also increases with formation
slowness because of lower contrast (better
coupling) between borehole and formation.
Therefore, the recommended transmitter frequency is 7.5 kHz in a hard formation, which
is reduced to 5 kHz in medium and soft formations, and 3 kHz in very soft formations
(Dtc > 130 msec/ft). Other parameters, such as
filter range, sampling interval and length
must be chosen accordingly.
6 16
E VA L U AT I O N A N D M O N I T O R I N G O F E X I S T I N G W E L L S
6 17
Figure 6.21
Clay wt%
100
Well 1
Well 2
Well 3
Well 4
Well 5
Well 6
Well 7
Well 8
Well 9
Well 10
Well 11
Well 12
50
Clay wt%
100
50
Clay wt%
100
50
0
0
100
200
Gamma ray
Synthetic GR (computed
from Th, U and K concentrations) plotted against total
clay (kaolinite, illite, smectite,
chlorite and glauconite)
measured on the same
sample for 12 data sets.
Wells 3 and 7 are from
Lake Maracaibo.
100
200
Gamma ray
100
200
Gamma ray
100
200
Gamma ray
6 18
E VA L U AT I O N A N D M O N I T O R I N G O F E X I S T I N G W E L L S
Figure 6.22
Clay (wt%)
100
50
0
0
10
20
Thorium (ppm)
10
Uranium (ppm)
2.5
Potassium (wt%)
Clay (wt%)
100
50
0
0
10
20
Aluminum (wt%)
Titanium (wt%)
10
Gadolinium (ppm)
Clay (wt%)
100
50
0
0
25
50
Silicon (wt%)
15
30
Iron (wt%)
20
40
Calcium (wt%)
Figure 6.23
100
Clay wt%
Well 1
Well 2
Well 3
Well 4
50
Clay wt%
100
Well 5
Well 6
Well 7
Well 8
Well 9
Well 10
Well 11
Well 12
50
Clay wt%
100
50
0
0
50
100 0
Estimated clay
50
100 0
Estimated clay
50
100 0
Estimated clay
6 19
GR spectroscopy
The same database was used to examine
other estimators of clay. Figure 6.22 shows
the relationship of various elements to clay
in Well 3. Of the three elements thorium,
uranium and potassium (Th, U and K) that
make up the natural GR signal, only K shows
a better correlation than the total GR.
Aluminum (Al) correlates well, and silica (Si)
anti-correlates well except for a few points
that have a high iron (Fe) percentage.
In general, the most reliable estimator
is found to be Al. However, it is a complicated
measurement to make in a borehole, currently
requiring a chemical source, two spectrometers and a measurement of formation capture cross section. K is not generally as
reliable as in the example in Fig. 6.22.
Therefore, an alternative estimator is used
which is based on the observation that as the
clay percentage increases, the amounts of
several elements decreases. This was already
noted for Si in Fig. 6.22. Similar anticorrelations can be observed with calcium
(Ca) in carbonate intervals, and with Fe where
siderite or pyrite are present. After converting
from elemental to oxide concentrations and
optimizing the coefficients for the data from
Wells 1-10, the following algorithm has been
found:
50
100
Estimated clay
Figure 6.24
100
50
0
0
50
100
(6)
Figure 6.25
10,000
10,100
10,200
Depth (ft)
10,300
10,400
10,500
10,600
10,700
10,800
0.1 0 0.1
Silicon
0.1 0 0.1
Calcium
0.2 0.4
Iron + .14Al
0.1
Sulfur
0.1 0.1
0.1
0.2 0.4
Titanium Gadolinium
Capture yields from the RST far detector in a Lake Maracaibo well, unsmoothed
and not corrected for casing or cement.
6 20
E VA L U AT I O N A N D M O N I T O R I N G O F E X I S T I N G W E L L S
Figure 6.26
350
300
Frequency
250
200
150
100
50
0
0
0.1
0.2
0.3
0.4
0.5
Fe yield
Figure 6.27
300
250
Frequency
200
150
100
50
0
0.1 0.05
0.05
0.1
0.15
0.2
0.25
Ca yield
capture Ca yield.
10,000
10,100
10,200
10,300
Depth (ft)
Wi = F Yi /Si
10,400
10,500
10,600
10,700
10,800
0.1 0 0.1
Silicon
0.1 0 0.1
Calcium
0.2
0.4
Iron + .14Al
0.1
Sulfur
0.1 0.1
0.1
Titanium
Capture yields from the far detector, after casing and cement correction.
0.2 0.4
Gadolinium
(7)
Figure 6.28
6 21
(8)
Figure 6.29
10,100
10,200
Depth (ft)
10,300
10,400
10,500
10,600
10,700
0
50 0
Silicon
(wt%)
20
Calcium
(wt%)
40
10
20
Iron + .14Al
(wt%)
10
Sulfur
(wt%)
20
Titanium
(wt%)
20
40
Gadolinium
(ppm)
Elemental Concentration Logs computed with SpectroLith are compared with cores.
Elements to minerals
Equations (5) and (6) provided a means
for obtaining minerals from elements, based
on core data. There are two complications
using log data. First, magnesium (Mg) is not
normally measured by spectroscopy logging
tools, although recent work shows this can be
done. However, as discussed above, the sum
(Ca + 1.455Mg) is independent of whether
there is calcite or dolomite, which is the same
as saying that this sum equals the Ca
measured by logs. Therefore, the log Ca can
be used to estimate the total carbonate.
Second, the optimization for clay percentage
needs to allow for the fact that the
Fe measurement is contaminated by Al.
Equations (5) and (6) are therefore modified
as follows:
Clay (log) =1.91 ( 100 - SiO2
- CaCO3 - 1.99Fe)
(9)
(10)
(11)
In the case of sandstones rich in feldspar or mica, such as Wells 11 and 12, the
alternative clay computation is as follows:
Clay2(log)=
-18.5+3.34(100-SiO2 -CaCO3 -1.99Fe)
(12)
6 22
E VA L U AT I O N A N D M O N I T O R I N G O F E X I S T I N G W E L L S
Figure 6.30
10,100
10,200
Depth (ft)
10,300
10,400
10,500
10,600
10,700
0
50
Clay (wt%)
100
50
100
Carbonate (wt%)
50
100
Quartz-Feld-Mica (wt%)
Figure 6.31
100
90
GR
Spectrolith
80
Clay (wt%)
70
60
50
40
30
20
10
0
X650
X700
X750
X800
X850
Depth (ft)
Clay weight percentage in a Orinoco Belt well from SpectroLith, GR and cores.
The SpectroLith value clearly agrees more closely with the cores.
6 23
X900
(13)
The SpectroLith processing was automatic. As can be seen, the two log estimates
agree well between X855 ft and X885 ft.
Elsewhere, and in spite of the scaling factor,
the GR gives too high a clay percentage by
about 10 to 20 wt%. Furthermore, it shows
excursions, such as at X760 ft and X815 ft,
which would be interpreted as sealing shales,
but which are not shales at all. SpectroLith
gives a better picture of the reservoir.
Summary
In shaly sandstone intervals, it is important to know the percentage of clay to estimate
porosity, saturation and permeability. Very
often the clay percentage is estimated from the
GR. However, it has been shown that the GR
is unreliable as a quantitative estimator of clay.
This is especially true in old wells where core
and other logs are not often available to
calibrate and help interpret the response.
An alternative technique is to use
capture GR spectroscopy yields from an RST
tool. Using the SpectroLith procedure, these
yields can be turned into elemental and
mineralogical concentrations. Based on
examples from Lake Maracaibo and the
Orinoco Belt, it has been shown that more
accurate clay percentages can be obtained by
this technique. This improvement was
significant in enhancing the reservoir
evaluation, so that potentially good zones
were not falsely judged to be poor quality
reservoir, and so that facies with high GR
were not treated as reservoir barriers.
(14)
(15)
(16)
(17)
6 24
E VA L U AT I O N A N D M O N I T O R I N G O F E X I S T I N G W E L L S
Pressure
One of the simplest and best established
methods of obtaining layer pressure is to
make openhole wireline formation tests in a
nearby new well. The key to success is a
good correlation between the old and new
well, and a good knowledge of sealing faults
and other factors that might cause a
discontinuity in pressure. Typically in old
wells the correlation is based on a good
geological model, and on GR and resistivity
logs. Here again the SpectroLith lithology
can improve results. Chemostratigraphy,
which is based on the elemental volumes
from SpectroLith, can help make correlations
that are better than those from GR and
resistivities alone (Herron, 1988).
Figure 6.32
10,500
10,600
10,700
Depth (ft)
10,800
10,900
11,000
11,100
Shell
core data
Maraven
core data
11,200
0
50
Clay (wt %)
6 25
100
0.1
10
100
1000
10,000
Permeability (md)
Example
The example shown in Fig. 6.33 is from
Bloque IV in the Maracaibo Basin. VLD-1112
is a new well and VLD-775 is an old well in
the same block. An extensive study of cores
and logs in this area led to a new definition
of layers for the reservoir, shown as layers I
to XII on the left. This layering was strongly
confirmed by the MDT* Modular Formation
Dynamics Tester pressures recorded in VLD1112. Some layers are clearly much more
depleted than others. The permeabilities
were estimated with the K-L method using
the openhole mineral volumes.
The old well, VLD-775, is about 1 km
SE of VLD-1112 and was drilled in 1978.
After producing about 5 MMBbl of oil it
recently showed a large increase in gas cut.
The fluids in the different layers at the time of
the original openhole and a recent RST are
summarized on the right of Fig 6.33. In this
case, the openhole logging suite was
sufficient to make a reasonable estimate of
porosity and lithology. The permeability in
track 2 was derived from the mineralogy and
the K-L model. The correlation between the
two wells is based on the new layer model,
Laguna
(gAPI)
-20
GR
120
f =19%
k =31 md
XI
f =18%
k = 4 md
f =19%
k =152 md
10,500
IX
10700
f =16%
k= 3 md
f= 16%
k = 2 md
VI
10,600
VII
VI
V
f= 22%
k =1498 md
IV
IV
10800
III
II
,
XI
VIII
VII
10000
XII
VIII
KINT (md)
1
10600
f =23%
k = 306 md
IX
md
ft
2 79
-80
10,400
XII
Lagunillas I nferior
SP
Hist + RST 4 94
GR
0
Layer
MDT Pressure (psig)
Average
1300
1050
Porosity and
K (md)
120 Depth
Permeability
(ft) 1.00
10000.
Layer
Sand
Layer
Member
2 79
6 94
12 91
VLD 775
VLD 1112
Openhole 12 78
Figure 6.33
f 24%
k 1585 md
10,700
III
II
Open perforations
Oil
Sequence boundary
Flooding surface
Other marker
Gas
Closed perforations
Water
Indirect estimates of dynamic properties in an old well (VLD-775). Permeability is from openhole measurements and pressure by correlation
with the new well VLD-1112, (after Bryant, et al, 1996).
6 26
(18)
(19)
where:
q = flow rate, hp = heavy phase
lp = light phase, y = phase hold-up,
vs = slippage velocity,
qs = slippage flow rate defined as
yhp ylp A Vs
In a bubble flow regime; i.e., when one
phase is continuous while the other is
distributed as bubbles within it, then it is
assumed that the spinner can determine the
total flow rate qt. The hold-up of the phases
has to be determined from measurements,
while the slippage velocity is found from
laboratory-established correlations.
6 27
(20)
(21)
(22)
Figure 6.34
Probe
Probe
Flow
6 28
E VA L U AT I O N A N D M O N I T O R I N G O F E X I S T I N G W E L L S
Figure 6.35
Water
Moved oil
Average BUB count
0
Depth
(ft)
(cps)
Fluid density
.5
Hold up
image 0
(g / cm3)
Sand
Spinner velocity
(rps)
1.5 0
Average holdup
(%)
Oil
Perforated
zones
500
20
Fluid temperature
1 177
(oF)
183 0
Water phase
Clay
Oil phase
Bound water
Flow rate
(B / D)
Combined model
750
(V / V)
B
C
6 29
Figure 6.36
Counts sec
600
400
200
10
20
30
40
50
Time (sec)
6 30
E VA L U AT I O N A N D M O N I T O R I N G O F E X I S T I N G W E L L S
Figure 6.37
LL3
ILM
ILD
0.02
( ohm-m )
200
GR
(gAPI)
Perforations
200
Depth
(ft)
11,800
Water velocity up
(ft min)
20
CPS cycle
Time (sec)
11,900
12,000
12,200
CPS cycle
12,100
Time (sec)
12,300
An example of a WFL
showing water flowing down
from an upper zone, and not
up from a lower zone.
6 31
Figure 6.38
Oil
6000
Water
Oil
1500
Water
Oil
600
Water
80
89
90
91
Figure 6.39
Combinable Production
Reservoir Saturation Tool
Logging Tool
Oil hold-up
Pressure & temperature Gas indicator
Gamma Ray
Detector
Dual DEFT
Spinner
Water velocity
Water hold-up index
Water flow-rate index
6 32
E VA L U AT I O N A N D M O N I T O R I N G O F E X I S T I N G W E L L S
Figure 6.40
0.5
Flow rate (B D)
6000
0.57
0.52
0.51
0.49
0.74
0.63
0.48
0.40
0.87
0.68
0.34
0.21
1500
0.5
600
1
80
89
90
91
The cross-sectional hold-up image maps are shown that were acquired with the
Dual-DEFT of the Flagship string during the same flow-loop experiment that is
shown in Fig. 6.38.
Figure 6.41
Marker signal
10
20
30
40
50
60
70
80
90
Time (sec)
Oil-miscible marker
Casing
Oil
Water
RST
An oil-soluble gadolinium-based marker is injected into the oil phase. The arrival of
the marker at the RST is detected by a perturbation of the capture cross section of
the borehole fluid. Time of flight is used to compute phase velocity. A water-soluble
marker can be ejected downward for water velocity.
6 33
Spinner measurement
Figure 6.42
500
400
300
200
100
14,000 B D
(5-in. liner)
Oil
Water
0
0
100
200
300
400
500
Auxiliary measurements
6 34
E VA L U AT I O N A N D M O N I T O R I N G O F E X I S T I N G W E L L S
6 35
Figure 6.43
Oil
Water
Normalized flow rate
Flow profile
2000
(B D)
Q T ratio
Computed oil flow rate
(B D)
WFL station water
flow rate
(B D)
0
7378.0ft
Hold-up image
(DFH)
Above 0.94
0.88 - 0.94
0.82 - 0.88
0.76 - 0.82
0.71 - 0.76
0.65 - 0.71
0.59 - 0.65
0.53 - 0.59
0.47 - 0.53
0.41 - 0.47
0.35 - 0.41
0.29 - 0.35
0.24 - 0.29
0.18 - 0.24
0.12 - 0.18
0.06 - 0.12
Below 0.06
10400
10300
10200
10000
9900
9800
9700
9600
9500
9400
9300
9200
9100
9000
8900
8800
8700
7383.0ft
Borehole gas
7378.0ft
Hold-up image
(DFH)
7383.0ft
Borehole gas
Above 0.94
0.88 - 0.94
0.82 - 0.88
0.76 - 0.82
0.71 - 0.76
0.65 - 0.71
0.59 - 0.65
0.53 - 0.59
0.47 - 0.53
0.41 - 0.47
0.35 - 0.41
0.29 - 0.35
0.24 - 0.29
0.18 - 0.24
0.12 - 0.18
0.06 - 0.12
Below 0.06
A Flagship log run in a horizontal well in Trinidad. The hold-up images of water and oil for flow and shut-in conditions are shown
together with gas hold-up indication. The interpreted flow profile of water and oil is shown. The cross-sectional hold-up images at
9200 ft and 9000 ft indicate a difference in oil and water hold-up due to deviation changes only.
6 36
HYDRAULIC INTEGRITY
CCC
BBB
AAA
,,,
@@@
,,,
CCC
BBB
AAA
,,,
@@@
,,,
CCC
BBB
AAA
,,,
@@@
,,,
,
,,
@@
A
BB
C
,
,,
@@
A
BB
C
Figure 6.44
Formation
Cement
Casing
Acoustic
beam
transducer,
b) casing resonance
Transducer
and reflections,
Rotation
a)
Metal
plate
c) received signal
with large first
dependent signals.
b)
Transducer
Zmud
Vmud
Zsteel
Vsteel
a
Mud
Casing
Zcem
Cement
Formation
Cement
Water
tc
Amplitude
c)
50
60
70
Time (sec)
6 37
80
90
HYDRAULIC INTEGRITY
Figure 6.45
4
Setting slurry
cement + mud
t
gh
tc
em
en
Slurr
Li
Impedance (MRAY)
Default
threshold
WBM
max
Fluids
Diesel
1.2
1.4
1.6
1.8
2.0
2.2
2.4
6 38
E VA L U AT I O N A N D M O N I T O R I N G O F E X I S T I N G W E L L S
Figure 6.46
Fluid compensated
CBL amplitude
0
(mv)
Acoustic
impedance
Max
50
Transit time
400
(ms)
Gamma ray
(gAPI)
100
Acoustic
impedance
Min
200
Liquid
Depth
(ft)
Avg
Gas
(MRAY)
Cement 0
10
Min
10
200
VDL variable
density Max
(ms)
1200
17,650
17,700
17,750
6 39
HYDRAULIC INTEGRITY
Figure 6.47
-500.0000
Depth
(ft)
(gAPI)
150
0.0760
Internal radii
minus average
(IRBK)
External
radius
average
(ERAV)
0
(in.)
External
radius
average
(ERAV)
150 0
(in.)
150
7553
7560.5
6 40
E VA L U AT I O N A N D M O N I T O R I N G O F E X I S T I N G W E L L S
Figure 6.48
200
Depth
(ft)
Travel time
(Ms)
1000
GR
(gAPI)
CBL amplitude
100
(mv)
Average map
0
100
(mv)
Casing attenuation
150 40
0
(mv)
0
20
Cement map
17,700
17,750
17,800
Figure 6.49
Surface sensors
and controls
Single-phase flow rate
Valve and chok e control
Pressure measurements
Sand detection
Field office
Production reports
Alarms
Problem identification
Integrated
applications
Downhole sensors
FloWatcher* monitors (flow
rate and density)
PumpWatcher* gauges
(ESP control)
PressureWatch* gauges
(pressure and
temperature)
Host server
and database
District office
Production engineering
Head office
Field performance
6 42
E VA L U AT I O N A N D M O N I T O R I N G O F E X I S T I N G W E L L S
Table 6.2
Main Applications
Description
Restricted access
Highly deviated wells
Pumping wells
Reservoir management
Production optimization
Figure 6.50
Section A
Cable
protector
Cable
Gauge
mandrel
Section A
Pressure and
temperature
gauge
Exploded view of
metal-to-metal seal
Wellhead
Wellhead
outlet
outlet
Interference testing
Reservoir pressure control
Transient well testing
History matching
As the technology evolves further, monitoring will be combined with downhole control
(intelligent completions), allowing even
greater control of production, which will find
particular applications in the new generation
of multilateral wells.
Several reservoirs in Venezuela already
benefit from permanent monitoring installations, and there are many more that could do
so. As seen in Chapter 2, Venezuela is distinguished worldwide by the large variety of its
reservoirs in terms of rock and fluid properties, drive mechanisms, depth, well
completions and production practices. There
are two broad groups of applications for
permanent monitoring systems (Table 6.2):
improved reservoir management and production optimization due to the availability
of continuous, real-time downhole data; and
difficult field or well conditions, in which the
problems of well intervention prevent the
recording of downhole data at all. The
following sections describe briefly the
technique, and show examples of applications in Venezuela and elsewhere.
6 43
Figure 6.51
Wellhead
outlet
Cable
pentrator
Control
line
Electrical
cable
Upper
Filter
Flapper
valve
Digital
permanent
gauge
Packer
Lower
6 44
E VA L U AT I O N A N D M O N I T O R I N G O F E X I S T I N G W E L L S
Figure 6.52
Well A
4900
Pressure pulses
4800
Pressure (psi)
4700
4600
4500
4400
4300
4200
4100
4000
Well B
Pressure (psi)
5060
5040
5020
Effect of
short pulses
5000
Effect of
long pulses
4980
4960
13
14
15
16
17
18
19
20
21
22
23
24
25
June
6 45
Reservoir management
Complex fields require careful reservoir
management. Two North Sea fields, Gullfaks
and Veslefrikk, illustrate how permanent
monitoring can help.
In Veslefrikk, commingled production
and injection was planned from the Brent and
Intra Dunlin Sand reservoirs to reduce total
investment. Here gauges are used in dedicated wells to monitor the two reservoirs
independently. Control is obtained by selective perforation in producers and downhole
chokes in injectors. A carefully planned data
acquisition program during the initial
production phase provided information about
reservoir properties, production potential and
well behavior. In addition, two of the largest
uncertainties were partially resolved: the
degree of communication across the main
arcuate fault, and the vertical transmissibility
between the Lower and Middle Brent through
the low-quality Rannoch sandstone.
Pressure (kPa)
Figure 6.53
28,000
27,000
3000
2500
1000
1500
2000
2500
3000
3500
4000
4500
Time (hr)
Production optimization
Careful optimization of the producing
conditions is essential in reservoirs with oil
near the bubble point, and with near-critical
fluids. More generally in high producing
wells, small improvements in the operating
point can significantly increase the production rate. These conditions exist in the
deep Oligocene fields of North Monagas, and
in the Center Lake area. Permanent pressure
data makes it easier to regularly monitor and
adjust the producing conditions. For example,
in the reservoir shown in Fig. 6.52, the bubble
point of the oil in one of the zones is 3700 psi
and the initial formation pressure is 4600 psi.
So drawdown had to be less than 900 psi to
sustain gas-free production. High skin factor
in the first well meant that as large a drawdown as possible would be needed for adequate production, introducing a further complication. However, the pressure could be
carefully monitored and production optimized to maintain reservoir pressure at
around 40 psi above bubble point.
The state-of-the-art
Permanent monitoring systems are an
important factor in a long-term effort to
improve recovery factors by continuous
monitoring and control of the key processes
in the reservoir. At the present time, most
systems allow continuous real-time measurements of pressure and temperature at the
bottom of the tubing. These measurements
alone offer significant advantages in three
areas: reservoir management, production
optimization, and where well conditions or
logistics make normal wireline interventions
too difficult or expensive.
Future directions
Current research is directed at extending
permanent in-situ installations to include
arrays with a large number of diverse
sensors (Fig. 6.54). Arrays of sensors
distributed at many points in the reservoir
will allow time-varying distributions of
reservoir properties, such as saturation and
pressure, to be continuously monitored.
These measurements can eventually be used
to continuously update the reservoir model
for optimum reservoir management. Figure
6.54 is a conceptual idea of how this might
eventually be done.
One direction is the development of
sensors in direct contact with the formation.
6 46
E VA L U AT I O N A N D M O N I T O R I N G O F E X I S T I N G W E L L S
Figure 6.54
I
A1
A2
Q1
P5
A3
P6
S1 S2 S3 S4 S5
Q3
R1 R2 R3 R4 R5
A4
P7
Q4
P1
Q2
P2
C1
C2
C3
P3
P4
T1
T2
T3
T4
CONCLUSION
6 47
REFERENCES
Saturation
Roscoe, B.A., Stoller, C., Adolph, R.A., Cheeseborough, J.C., Hall, J.S.,
McKeon, D.C., Pittman, D., Seeman, B., Thomas, S.R., 1991, A new throughtubing oil saturation measurement system, Paper SPE 21413, presented at
the International Arctic Technology Conference, Anchorage, May 29-31, and
the Middle East Oil Show and Conference, Bahrain, November 1991.
Bryant, I.D., Baygun, B., Coll, C., Cordova, P., Gamero de Villaroel, H.,
Hernandez, L., Herron, M.M., Jimenez, Z., Leon, K., Matteson, A.,
Ramamoorthy, R., Rondon, L., Stewart, L., and Tariq, S.M., 1996, Integration
of old and new measurements to optimize redevelopment of the Lower
Lagunillas reservoir of Bloque IV, Lake Maracaibo, Venezuela, SPE 36096,
Presented at the Fourth Latin American and Caribbean Petroleum
Conference, Port-0f-Spain, Trinidad, April.
Johnson, D.L., Koplik, J., Schwartz, L.M., 1986, New pore-size parameter
characterizing transport in porous media, Phys. Rev. Let, 57, p 2564-2567.
Stoller, C., Scott, H. D., Plasek, R. E., Lucas, A. J., and Adolph, R. A., 1993,
Field tests of a slim carbon/oxygen tool for reservoir saturation monitoring,
SPE Paper 25375, Presented at the Asia Pacific Oil & Gas Conference &
Exhibition, Singapore, February 8-10.
Halford, F. R., MacKay, S., Barnett, S., and Petler, J. S., 1996, Production logging measurement of distributed local phase holdup, SPE Paper 35556 A,
Society of Petroleum Engineers Annual Technical Conference and Exhibition
held in Denver, Colorado, U.S.A., October 6-9.
Roscoe, B., Lenn, C., Jones, T. G. J., and Whittaker, C., 1996, Measurement
of the oil and water flow rates in a horizontal well using chemical markers and
a pulsed-neutron tool, SPE Paper 36563, Society of Petroleum Engineers
Annual Technical Conference and Exhibition held in Denver, Colorado, U.S.A.,
October 6-9.
Locke, J. and Butler, J., 1993, Characterization of rock formations for the
improved calibration of nuclear logging tools, Paper R, Transactions of the
15th European Formation Evaluation Symposium, Stavanger, Norway, May.
Plasek, R. E., Adolph, R. A., Stoller, C., Willis, D. J., Bordon, E. E., and Portal,
M. G., 1995, Improved pulsed neutron capture logging with slim carbon-oxygen tools, SPE Paper 30598, Society of Petroleum Engineers Annual Technical
Conference and Exhibition held in Dallas, Texas, U.S.A., October 22-25.
Hydraulic Integrity
Hayman, A. J., Hutin, R., and Wright, P. V., 1991, High resolution cementation
and corrosion imaging by ultrasound 1191, Paper KK, Transactions of the
Society of Professional Well Log Analysts 32nd Annual Symposium, Midland,
Texas, U.S.A., October.
Hayman, A. J., Parent, P., Rouault, G., Zurquiyah, S., Verges, P., Liang, K.,
Stanke, F. E., and Herve, P., 1995, Developments in corrosion logging using
ultrasonic imaging, Paper W, Transactions of the Society of Professional Well
Log Analysts 36th Annual Logging Symposium, Paris, France, June.
Silva, M. G. P., Miranda, C. R., and Vincente, R., 1996, Optimization of
cementing and acidizing operations is the key to minimize acid damage to the
cemented annulus, SPE Paper 36112, Presented at the 4th Latin American
Petroleum Engineering Conference, Port of Spain, April.
Permanent monitoring
Oilfield Review, Winter 1995, p. 32-47.
The Search for Oil and Gas in Latin America & the Caribbean, No. 5, 1997,
p. 28-41.
6 48