Bona 2001
Bona 2001
Bona 2001
Summary in contact with the oil and the location of the two fluids is the
This paper summarizes an experimental and theoretical study opposite of the water-wet case.
conducted to investigate the effect of rock wettability on the elec- Putting an electric field on the rock will produce different
trical response of water- and oil-saturated core samples in a wide effects depending on the wettability. In the water-wet case, the ions
frequency band. The aim of the study was to evaluate a technique dissolved in the water will be able to move in the film of water that
not previously applied to wettability determination as a practical coats the pores and will travel fairly long distances before being
alternative to conventional tests, and to ascertain the consequent trapped within a water/rock or water/oil interface. In an oil-wet
economic and technical benefits in reservoir appraisal. system (at irreducible water saturation), the distance traveled by
the ions, contained in spherical water droplets inside the larger
Introduction pores or in globules that extend over several pores, will have more
There are two main problems in determining reservoir wettability. or less the same diameter as these pores. The difference in the
One is that current methods are time-consuming. It takes several shape of the water inclusions creates two different electrical
weeks to carry out a wettability test using the Amott or USBM responses and allows determination of rock wettability.
technique, and only a limited number of plug samples can be ana- The key parameters are the dielectric constant, e, and the electri-
lyzed for a reservoir study. The other problem is that wettability cal conductivity, s. Both describe the interaction of the rock with the
may be altered by extracting the core from the reservoir, bringing electric field and are complex functions of the angular frequency w
it to the surface, transporting it to the laboratory, and preparing of the field. In mathematical terms, we conventionally write:
samples for analysis. This influences the quality of subsequent
measurements, particularly relative permeabilities results. No sat- e(w)=e'(w) - i e''(w) . . . . . . . . . . . . . . . . . . . . . . . . . . . . .(1a)
isfactory solution has yet been found for the first problem. and s(w)=s'(w) + i s''(w), . . . . . . . . . . . . . . . . . . . . . . . . . . .(1b)
Regarding the second, some companies use bland drilling muds
and take great care during core handling and cleaning to minimize where i= −1 is the imaginary unity (note the minus sign in the
wettability alteration; others prefer to restore the wettability by first of the two expressions). Quantities e(w) and s(w) are perfect-
aging the core in the reservoir fluid at reservoir conditions. ly equivalent, being related by the simple relationship
However, we cannot be sure in either case that representative core
s(w)=iwe0e(w), . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .(1c)
samples are obtained.
Dielectric constant measurements produce data that correlate which follows directly from Maxwell’s equations (e0=8.854´10-12
well with rock wettability. They have the advantage of taking just Farad/m is the permittivity of the vacuum). In practice, if we know
a few hours and are easy to perform, consisting of standard electri- one quantity, we can immediately calculate the other. The real part of
cal impedance tests similar to those usually performed to deter- the dielectric constant e and the imaginary part of conductivity s are
mine the Archie parameters. Furthermore, in-situ wettability can in a measure of how much energy from the electric field is stored in the
theory be determined with the dielectric results, comparing the for- rock; the imaginary part of e and the real part of s are a measure of
mation conductivity obtained from a high-frequency tool with the how dissipative the rock is in relation to the electric field. A widely
conductivity obtained from a low-frequency one. used parameter in the oil industry is the dc conductivity, sdc, which
appears, for example, in Archie’s laws; sdc is the low-frequency limit
Theory and Definitions of the real part of the complex conductivity s(w), and—according to
The effects of wettability on the electric properties of a rock can be Eqs. 1a and 1c—it can be written as we0e''(w), which provides the
divided into two categories: connectivity effects that dominate at low key to expressing Archie’s laws in terms of dielectric constant.
frequencies (below 1 kHz), and shape effects that appear at interme- Although several dielectric mechanisms contribute to the com-
diate to high frequencies (above 10 kHz). Many studies have been plex dielectric constant e, a previous study4 showed that the response
made in the past on low-frequency effects1,2 (low-frequency effects of the rock in the 100 Hz to 10 GHz frequency range can be
control the Archie parameters), while few studies have been described well by overlapping the contributions of the three terms: a
published about the intermediate- /high-frequency effects. This is real constant (e¥), a complex function of the frequency (eMW(w)) that
probably owing to the difficulty in performing wide-band, high- describes Maxwell-Wagner polarization, and a complex Power-Law
frequency electrical measurements. This study takes into considera- function (ePL(w)) associated with charge transport. Thus,
tion both the low- and intermediate- /high-frequency effects.
To understand how wettability affects the electrical response of A (M ) = A∞ + AMW (M ) + APL (M ), . . . . . . . . . . . . . . . . . . . . . . .(2a)
a rock, consider that wettability governs the geometrical distribu-
tion of the water phase in the pore space. According to W.G. where
Anderson,3 when the rock is water-wet, water tends to occupy the
small pores and contact most of the rock surface. The water will ∆A
AMW (M ) = >
. . . . . . . . . . . . . . . . . . . . . . . . . .(2b)
form a continuous phase down to very low saturations, while the 1 + (iMJ )1−=
oil will be isolated in pockets even at rather high oil saturations.
In an oil-wet system, on the other hand, the rock is preferentially N
B (iM )
and APL (M ) = . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (2c)
Copyright © 2001 Society of Petroleum Engineers
iMA0
This paper (SPE 69741) was revised for publication from SPE 50720, first presented at the In the above expressions, e¥, De, t, a, b (0£a£b£1), B, and
1999 SPE International Symposium on Oilfield Chemistry, Houston, 16–19 February.
Original manuscript received for review 24 August 1999. Revised manuscript received 6
N (0£N£0.5) are all fitting parameters. Fig. 1 shows the trend of
November 2000. Manuscript peer approved 20 November 2000. the functions eMW(w) and ePL(w) and their sum for a given choice
1E+03
1E+02
1E+01
1E+00
1E-01
1E+02 1E+03 1E+04 1E+05 1E+06 1E+07 1E+08 1E+09 1E+10
frequency, Hz
1E+03
Fig. 2—Diagram of Maxwell-Wagner (MW) polarization in a
1E+02
brine/oil/rock system at irreducible water saturation. A water-
wet condition is assumed.
1E+01
frequency, Hz
and charge transport would be a nondispersive process. However,
this never happens in real rocks. The ions are always loosely bound
(i.e., confined in zones limited by potential barriers that are difficult
Fig. 1—Frequency response of the functions eMW(w), ePL(w), and to overcome), and consequently N>0; as a result, ePL(w) has a real
their sum. Above: real components; below: imaginary compo- part that is not zero and there is an effect on the storage factors e' and
nents. Values attributed to the parameters: De=315; t=0.123 s'', leading to dispersive conduction mechanisms.
ms; a=0.26; b=1; B=4.22´10-3; and N=0.019. PL mechanisms indicate that loss and storage decrease with fre-
quency (Fig. 1). A peculiar characteristic of a PL mechanism is that
of fitting parameters in the 100 Hz to 10 GHz frequency interval. the phase does not depend on frequency. Such behavior has been
The real and imaginary parts are illustrated for both functions. As connected to hierarchical self-similarity (many-bodies effects) in
shown, charge transport is the prevalent mechanism at the low fre- the dynamics of the microscopic relaxation that originates from the
quencies, while Maxwell-Wagner polarization controls the rock dispersive nature of charge transport or to the fractal character of
response at intermediate and high frequencies. the geometric microstructure of the system. However, its origin is
not yet clear and no conclusive results exist, particularly in terms of
Maxwell-Wagner Polarization. The Maxwell-Wagner (MW) theoretical models. It is noteworthy that most of the experimental
effect5-8 is described by the function eMW (w). This dielectric mech- evidence of the PL behavior of dielectric constants comes from
anism is characteristic of heterogeneous systems like water- and studies on rocks; we refer to the literature for further details.6,9-14
oil-saturated rocks, where insulating and conductive phases coex- Multiplying Eq. 2c by iwe0 and taking the real part of the resulting
ist. Under the effect of an external electric field, the ions dissolved expression gives a real conductivity at low frequency of the form
in the water concentrate in the regions where electrical properties s'(w)=BwN cos(Np/2) [the imaginary unit i can be expressed as
undergo sudden changes (Fig. 2a). These charge concentrations exp(ip/2), thus i N=exp(iNp/2)=cos(Np/2)+isin(Np/2)]. Because
produce interfacial or MW polarization of the water inclusions. N is usually very small, the above formula is appropriate in the 100
MW polarization essentially occurs at the oil/water interfaces or on Hz to 1 kHz frequency range, and the dc conductivity of the system
pore walls. When the electric field is switched off, the ions return may be obtained with reasonable accuracy from the s' value at 100
to their initial random distribution (Fig. 2b). The time constant that Hz. In practice, the final expression is sdc= B×628Ncos(1.57N),
controls the process is termed relaxation time and is usually indi- given that w is equal to 2p times the frequency. We expect that rock
cated by t; the reciprocal of t is termed the relaxation frequency. wettability affects both parameters B and N.
The relaxation frequency for MW polarization lies between 10 kHz
and 100 MHz, depending on the type of rock, the water salinity, Materials
and the spatial distribution of the water. Measurements of the complex dielectric constant were carried out on
At low frequencies the alternating electric field is slow enough 39 samples. These included eight sintered glass filters, four Berea
for the ions to follow the field variations. Because the polarization sandstone samples, and 27 reservoir samples (both sandstone and car-
is able to develop fully, eMW'(w) assumes its highest value, and the bonate) taken from nine different fields. The samples were analyzed
loss factor eMW"(w) is directly proportional to the frequency. As the in two different wetting conditions, obtained in the following ways:
frequency increases, eMW"(w) continues to increase, but the storage • Water-wet samples. This group comprised glass filters and
factor eMW'(w) begins to decrease owing to the phase delay Berea samples with their original, strong water wettability, and
between the fluctuations of the ions and the electric field. Above reservoir samples cleaned with routine extraction (the cleaning
the relaxation frequency, both eMW"(w) and eMW'(w) drop off, as the resulted in alteration of original wettability to water-wet).
electric field is too fast to influence the diffusion of the ions, and • Oil-wet samples. This group comprised glass filters and Berea
the MW polarization disappears. samples made oil-wet through silanization, and reservoir samples
that were extracted and then silanized. A method published by N.E.
Charge Transport. Charge transport at low frequencies (<10 kHz) Takach et al.15 was used to make the samples oil-wet, which involves
is modeled with the function ePL(w). This is a Power-Law (PL) adsorbtion of bis(dimethylamino)dimethylsilane onto the solid sur-
effect, typical of a large number of disordered solids where the faces of the samples at a temperature of 290°C in a vacuum oven. Of
charge carriers are weakly bound or confined in regions of low the many methods tested (aging the samples in crude oil, fluxing
mobility.9 Electrical transport is a discontinuous process in these sys- with asphaltenes, and other silanization processes), this proved to be
tems owing to charge hopping from one low-mobility region to the only efficient one to produce a strong oil-wet condition.
All the samples analyzed were disc-shaped with a diameter of 1 MHz frequency range (all the samples), and Time Domain
or 2 in. and a thickness of between 0.65 and 0.95 cm. The 27 reser- Reflectometry (TDR) in the 25 MHz to 10 GHz frequency range
voir samples were taken from nine homogeneous core plugs (one (glass filters and samples from reservoirs K and B).
per reservoir). Every plug was sliced into two to four discs. Table
1 shows the main petrophysical characteristics of the samples ana- Measurements in the 100 Hz to 40 MHz Range. Measuring the
lyzed. The Amott wettability index is available for some of the impedance of the sample is the basis for measuring the dielectric
glass filters and Berea samples. constant in this frequency range. From the measured impedance,
The samples were saturated with brine and then centrifuged in sample area, and thickness, the dielectric constant can be calculated.
an oil environment (capillary pressure applied: 0.5 to 1 bar). A syn- An HP 4194A Impedance Analyzer was used for these measure-
thetic brine with a salinity of 75.8 g/l NaCl and a refined mineral ments. The samples were placed in a parallel golden plate capacitor
oil (Soltrol by Phillips Co.) were used. Soltrol has an ideal dielec- specifically designed for the experiment, and a standard two-contact
tric behavior, maintaining a static dielectric constant throughout method was used. The measurement is easy, rapid, and cost-
the frequency range investigated and having negligible losses.
Table 2 shows the electrical properties of the fluids used. TABLE 2—ELECTRICAL PROPERTIES OF THE FLUIDS (25°C)
100 Hz
Experimental Equipment and Procedure Static Dielectric Const. Conductivity
The complex dielectric constant of the 39 samples listed in Table 1 (dimensionless) (S/m)
was measured in the 100 Hz to 10 GHz frequency range, using two Brine 67.0 9.5
different techniques: impedance measurements in the 100 Hz to 40 Soltrol 2.015 -
HP 4194A sample
Sampling Signal
head generator
coaxial line
fringing
ZL metallic
cella field
flange
portacampioni
TDR 10 0
3
10
TDR
10 2
ε' 10 -1
tan(δ)
10 1 f [Hz]
f [Hz]
10 -2
2 3 4 5 6 7 8 9 10
10 10 10 10 10 10 10 10 10 10 2 10 3 10 4 10 5 10 6 10 7 10 8 10 9 10 10
10 1
10 6
(b) (b)
5
10
10 4 10 0
TDR
TDR
tan(δ)
10 3
2
10
1
ε" 10 -1
10
f [Hz]
0
10
f [Hz]
10 -2
10 2 10 3 10 4 10 5 10 6 10 7 10 8 10 9 10 10 10 2 10 3 10 4 10 5 10 6 10 7 10 8 10 9 10 10
Fig. 6—Real and imaginary part of the dielectric constant meas- Fig. 7—Loss tangent angle tan(d) measured for two glass4 sam-
ured on four sandstone samples. Open circles: oil-wet sample ples (above) and two sandstone samples (below). Open circles:
(K3). Open squares: oil-wet sample (K98). Solid up triangles: water- oil-wet samples glass4C and K3. Solid triangles: water-wet
wet sample (K13). Solid down triangles: water-wet sample (K99). samples glass4H and K13.
of glass4: one sample (C) was oil-wet, and the other (H) was water- ples this effect is more intense than in the oil-wet ones, as observed
wet. Both the real part (e') and the imaginary part (e") of the meas- on the glass samples.
ured dielectric constant are plotted as a function of the frequency. At higher frequencies, the MW effect becomes important: both
The water-wet sample has higher e' and e" values than the oil-wet e' and e" show a change in slope. This occurs at around 10 kHz in
across the frequency range studied, despite the fact that the water the oil-wet samples and at higher frequencies in the water-wet ones.
saturations of the two samples are similar. In the low-frequency Both e' and e" level off, but at 100 kHz (oil-wet samples) and at
range, e' and e" decrease with increasing frequency and both samples about 5 MHz (water-wet samples) there is a second change in slope.
exhibit a PL response, dominated by the function ePL(w). Charge At these frequencies the MW polarization is at its peak and the
transport is the dominant mechanism in this frequency region; response of the samples is dominated by the eMW(w) function. The
according to Eq. 2c, both ePL'(w) and ePL"(w) show a straight line of reciprocal of the frequency at which the second change in slope
negative slope (N-1) when plotted on a bilogarithmic scale, which occurs gives the relaxation time t for the MW effect: t»10 ms is
is consistent with the experimental data. The contribution of the obtained in the oil-wet case, and t»0.1 ms in the water-wet one.
charge transport ePL(w) is much greater in the water-wet sample, At very high frequencies, e' tends to a constant value (e¥),
where the water forms a continuous phase, than in the oil-wet one. which is higher for the water-wet samples, while e" goes to zero.
At higher frequencies, e' and e" exhibit a change in slope. The The oil-wet sample K3 shows a second, less evident relaxation
MW effect begins to come into play and the charge transport ePL(w) between 100 MHz and 1 GHz: e" has a relative maximum, and e'
no longer dominates the dielectric response of the samples. In the shows a corresponding decrease. We think that this second relax-
oil-wet case, e" shows a peak at about 100 MHz, whereas in the ation is also the result of an MW polarization effect, with an asso-
water-wet case the relaxation peak takes place at around 75 kHz ciated relaxation time of around 1 to 10 ns. This fact will be dealt
(this peak can be noted only by a change in slope of e", as it is still with further in the section “Discussion of Results.”
masked by the PL contribution). At very high frequencies, the stor-
age e' tends to a constant value (e¥), while the loss e" goes to zero. Loss Tangent. When the dielectric constant is drawn as a simple vec-
tor diagram, the real and imaginary components are 90° out of phase.
Dielectric Constant of Reservoir Samples. Fig. 6 shows the The vector sum forms an angle d with the real axis e'. The loss tangent
dielectric constant plots measured on the four samples from reser- is defined as tan(d)=e"/e'. It is a measure of the relative loss tendency
voir K. Two of the samples (K3 and K98) were oil-wet, and two of the rock, being the ratio of the energy lost to the energy stored.
(K13 and K99) were water-wet. Both the real part (e') and the Fig. 7 shows the plot of tan(d) vs. frequency in the 100 Hz to 10
imaginary part (e") are plotted as a function of the frequency. GHz frequency range for (a) the glass4 samples and (b) two of the K
The water-wet samples have higher e' and e" values than the oil- reservoir samples. Two regimes can be identified: a low-frequency
wet ones only in the frequency range studied with the Impedance region where the charge transport dominates, characterized by a PL
Analyzer HP4194A (<40 MHz). In the low-frequency range, both that gives tan(d)=cot(Np/2) values that are independent of the
e' and e" decrease for all samples in an almost linear way; the elec- frequency, and a high-frequency region where the effects of the trans-
tric response is dominated, as in the case of the glass samples, by port and MW polarization, easily identifiable by the relative maxi-
a PL effect ePL(w), with 0<N<0.5 (Eq. 2c). In the water-wet sam- mum, are summed. Our results show that the water-wet samples
TABLE 3—VALUES OF THE PARAMETERS OF EQS. 2A THROUGH 2C OBTAINED BY FITTING THE EXPERIMENTAL DATA
Reservoir Sample Freq. range τ σ dc fmax
Code Code Wettab investigated ∆ε ε∞ (µs) α β B N (S/m) Atan( φ) (MHz)
glass3
C o-w 100Hz-10GHz 0.5 3.2 8.1 ×10 –5 0.4 1 8.8 ×10 –8 0.02 1.1 ×10 –7 0.040 1.9 ×10 3
D o-w 0.5 3.5 1.0 ×10–4 0.4 1 1.8 ×10 –7 0.03 2.3 ×10 –7 0.037 1.8 ×10 3
G w-w 49.6 7.1 3.1 0.4 1 4.8 ×10 –6 0.085 8.1 ×10 –6 0.55 0.22
H w-w 50.2 7.5 4.3 0.4 1 3.3 ×10 –6 0.14 7.9 ×10 –6 0.62 0.079
glass4
C o-w 100Hz-10GHz 0.6 3.7 1.35 × 10 –3 0.4 1 3.8 ×10 –8 0.04 6.0 ×10 –8 0.043 150
D o-w 0.6 3.8 1.4 × 10 –3 0.4 1 9.0 ×10 –9 0.035 1.7 ×10 –8 0.041 145
G w-w 53.7 9.6 1.3 0.4 1 4.1 ×10 –5 0.032 5.1 ×10 –5 0.59 0.22
H w-w 70.7 6.9 2.1 0.4 1 8.5 ×10 –6 0.11 1.7 ×10 –5 0.66 0.37
Berea
C o-w 100Hz-40MHz 17.7 2.1 7.2 ×10 –2 0.5 1 3.2 ×10 –9 0.30 1.1 ×10 –7 0.40 10.7
D o-w 21.8 2.7 8.7 ×10 –1 0.57 1 7.8 ×10 –9 0.25 1.0 ×10 –7 0.35 1.75
G w-w 111 3.0 3.0 ×10 –1 0.5 1 1.7 ×10 –3 0.13 3.8 ×10 –3 — —
H w-w 106 3.7 2.8 ×10 –2 0.2 1 2.4 ×10 –3 0.07 2.5 ×10 –3 — —
A
2 o-w 100Hz-40MHz 80 50.8 1.8 ×10 –3 2.10 ×10 –4 1 5.49 ×10 –5 0.020 7.3 ×10 –2 — —
9 w-w 198 24.5 3.9 ×10 –3 0.13 1 2.93 ×10 –2 0.010 3.3 ×10 –2 1.85 8.5
D
6 o-w 100Hz-40MHz 248 9.4 4.78 0.42 1 2.17 ×10 –5 0.125 4.9 ×10 –5 0.85 0.34
8 w-w 110 10.5 0.131 0.18 1 6.91 ×10 –4 0.058 1.0 ×10 –3 1.205 4.5
P
10 o-w 100Hz-40MHz 173 12.5 2.750 0.26 0.78 5.42 ×10 –5 0.057 7.8 ×10 –5 0.80 0.53
16 w-w 315 14.2 0.123 0.26 1 4.22 ×10 –3 0.019 4.8 ×10 –3 1.26 8.2
V
8 o-w 100Hz-40MHz 186 12.4 4.65 0.39 1 5.52 ×10 –5 0.19 1.9 ×10 –4 0.84 2.1
16 w-w 208 15.9 21 ×10 –3 0.21 1 9.23 ×10 –3 0.105 1.8 ×10 –2 1.6 27
K
3 o-w 100Hz-10GHz 250 8.2 9.13 0.44 1 1.93 ×10 –5 0.077 3.1 ×10 –5 0.80 0.14
3.3 3.45 ×10 –4 0.2 1 0.19 280
98 o-w 100Hz-40MHz 223 5.6 6.89 0.5 1 5.42 ×10 –5 0.087 9.5 ×10 –5 — —
13 w-w 100Hz-10GHz 107 11.8 0.13 0.19 1 1.07 ×10 –3 0.097 2.0 ×10 –3 1.46 3.0
99 w-w 100Hz-40MHz 130 8.16 0.109 0.20 1 3.26 ×10 –4 0.091 5.9 ×10 –4 1.24 7.7
L
10 o-w 100Hz-40MHz 177 12.9 0.34 0.32 1 1.42 ×10 –3 0.031 1.7 ×10 –3 — —
98 o-w 88.4 7.1 0.24 0.42 1 1.21 ×10 –4 0.09 2.2 ×10 –4 0.68 3.8
17 w-w 61 13.5 29.0 ×10 –3 0.44 1 5.22 ×10 –2 0.015 5.7 ×10 –2 — —
99 w-w 72 14.1 24.0 ×10 –3 0.13 1 8.20 ×10 –4 0.103 1.6 ×10 –3 0.96 16
B
7 o-w 100Hz-10GHz 64.6 6.4 48.3 ×10 –3 0.18 0.77 7.06 ×10 –6 0.224 2.9 ×10 –5 0.67 17
2.3 0.25 ×10 –3 0.2 1 — —
11 w-w 222 4.6 14.7 ×10 –3 0.24 1 2.23 ×10 –3 0.059 4.0 ×10 –3 1.3 120
23 w-w 137 15.5 24.7 ×10 –3 0.21 1 2.01 ×10 –2 0.016 2.5 ×10 –2 — —
PL
9 o-w 100Hz-40MHz 13 23 5.2 ×10 –3 0.48 1 5.10 ×10 –2 0.014 5.8 ×10 –2 — —
11 o-w 70 6.4 9.3 ×10 –3 0.56 1 2.08 ×10 –2 0.0034 3.4 ×10 –2 — —
6 w-w 159 17.1 1.50 0.01 1 0.4 1.4 ×10 –4 0.40 — —
16 w-w 110 22.2 0.805 0.11 1 0.24 2.0 ×10 –4 0.24
— —
R
7 o-w 100Hz-40MHz 64 8.9 3.80 0 0.6 4.74 ×10 –6 0.155 1.3 ×10 –5 0.54 0.13
24.4 0.073 0.34 1 0.50 3.4
98 o-w 57.6 4.4 0.187 0.36 1 2.20 ×10 –6 0.187 7.2 ×10 –6 0.72 6.6
11 w-w 66.4 9.8 0.014 0.12 1 1.50 ×10 –4 0.151 3.96 ×10 –4 0.98 36
99 w-w 32 9.2 0.0084 0.21 1 2.01 ×10 –5 0.194 6.89 ×10 –5 — —
Note: w-w= water-wet; o-w= oil-wet