Egyptian Journal of Petroleum (2014) 23, 4551

Egyptian Petroleum Research Institute

Egyptian Journal of Petroleum

www.elsevier.com/locate/egyjp
www.sciencedirect.com

FULL LENGTH ARTICLE

Using of Picketts plot in determining the reservoir

characteristics in Abu Roash Formation, El-Razzak
Oil Field, North Western Desert, Egypt
A.A. El-Khadragy a, M.A. Ghorab b, T.F. Shazly b, M. Ramadan b,
M.Z. El-Sawy b,*
a
b

Geology Department, Faculty of Sciences, Zagazig University, Egypt

Department of Exploration, Egyptian Petroleum Research Institute, Egypt

Received 24 June 2013; accepted 16 September 2013

Available online 27 March 2014

KEYWORDS
North Western Desert;
Well logging;
Pickett plot;
Cross plot;
Porosity;
Permeability

Abstract The present work was devoted to evaluate the reservoir characteristics of the Abu Roash
Formation in the vicinity of El-Razzak Oil Field, North Western Desert, Egypt. The area of study is
bounded by latitudes 30 250 and 30 550 N and longitudes 27 500 and 28 400 E. Nine distributed
wells were utilized for this study.
A new technique has been applied through Picketts plot, to develop some of reservoir petrophysical parameters. These parameters include capillary pressure, pore throat aperture radii, height
above the free water table and bulk volume of water. This technique depends on the use of log
log plots of effective porosity versus resistivity combined with empirical relationships for calculating
the capillary pressure expressed as a function of permeability, porosity and water saturation. Also,
this technique gave the values of petrophysical exponents (m, n and a) which were used to calculate
the accurate value of water saturation in both clean and shaly rocks and then adjust estimation of
hydrocarbon saturation. The integration of these petrophysical parameters on a loglog graph of
porosity versus resistivity gives the importance for Pickett plot to be used in reservoir interpretation.
2014 Production and hosting by Elsevier B.V. on behalf of Egyptian Petroleum Research Institute.

1. Introduction
* Corresponding author.
E-mail address: [email protected] (M.Z. El-Sawy).
Peer review under responsibility of Egyptian Petroleum Research
Institute.

Production and hosting by Elsevier

The use of Pickett plots (ResistivityPorosity Relationship) as

an aid in characterizing petrophysical ow units in carbonate
reservoirs has been documented by [1] and [2]. This technique
was utilized for calculating some of reservoir petrophysical
parameters of Abu Roash Formation (G member) for nine
wells namely: (IG30-1, IJ30-1, IG33-1, IH35-2, IG34-1, IF321, IF34-2, IE34-6 and IE36-1) scattered in the oil eld as

1110-0621 2014 Production and hosting by Elsevier B.V. on behalf of Egyptian Petroleum Research Institute.
http://dx.doi.org/10.1016/j.ejpe.2014.02.007

46

A.A. El-Khadragy et al.

shown in Figs. 1 and 2. Sanyal and Ellithorpe [3] and Greengold [4] have shown that, Picketts plot should result in a
straight line with a slope equal to m.
This plot shows a useful model for putting together the
petrophysical parameters including water saturation, permeability, capillary pressure, pore throat aperture radii and
height above free water table [5]. This technique highers
Picketts plot as one of the most important plots for reservoir
evaluation, where the use of this plot enables the log analyst to
see within one single plot several keys of geological and
reservoir engineering parameters.

from Picketts plot for the G member of Abu Roash Formation are tabulated in Table 1.
Many lines are drawn at different values of porosity (90%,
60%, and 20%) equivalent to the rst drawn line.
Aguilera [5] used the following steps to show how to construct Picketts plot incorporating capillary pressure, pore
throat aperture radii and height above the free water table,
where the technique is applied on the nine wells for the G
Member of Abu Roash Formation in the El-Razzak Oil Field.
IE36-1 well is taken as an example in this paper.
2.1. Permeability

2. The theory of Picketts plot

SW I1=n

I Rt =R0 Rt =FRw

Reservoir rock must have the ability to allow petroleum uids

to ow through its interconnected pores. The rocks ability to
conduct uids is termed as permeability which depends on its
effective porosity.
Many empirical equations are derived to estimate the relative permeability. The following equation estimates the permeability [9] in terms of irreducible water saturation, as follows:

F aUm
t

K1=2 250U3t =Swi

The water saturation in Eq. (6) is irreducible, that corresponds to the beginning of a Krw, which equals zero.
The irreducible water saturation, in turn, can be solved by
using the following equation:

The theory of this plot started with Archies formulae [6]. By

rearranging the Archie equation we get:

Eqs. (1) and (3) can be combined to yield.

Rt

aUm
t Rw I

aUm
t Rw Sw

By using the logarithm with base 10, the Eq. (4) [7] leads
log Rt m log Ut logaRW log I

This is the equation of a straight line on loglog paper. The

line has a slope of (m) which is determined manually by measuring a distance on the Rt axis (in cm) and dividing it by the
corresponding distance on the porosity axis. The intercept
when PHI = 1 is the value of aRw as shown in Fig. 3, and
by knowing the value of Rw, the value of tortuosity factor
(a) can be determined.
According to [8], the saturation exponent n which is a
function of water saturation equals the value of porosity exponent m. The results of the petrophysical exponents obtained

Figure 1

Swi Ut  Sw =Ueff

After incorporating into Eq. (4), the following formula is

obtained [10]
3
1=2
Rt aUm

t Rw 250Ut =K

n

or Rt aU3nm
Rw 250=K1=2
t

8
n

By having the logarithm of both sides, the proceeding equation becomes [11]:
n

log Rt 3n  m log Ut logaRw 250=K1=2 

Generalized sketch map of Egypt showing the location of the studied area.

10

Picketts plot in determining the reservoir characteristics

Figure 2

47

Location map showing the distribution of the studied wells.

Porosity

aRw=0.025

0.1

m=x/y=1.7
0.01
0.01

Figure 3

0.1

10

Rt (ohm.m) 100

Picketts plot for G Member of Abu Roash Formation in IE 36-1 well.

Table 1 The results of petrophysical exponents and slopes obtained from Picketts plot for the G Member of Abu Roash Formation
in the nine wells.
No.

Well name

aRw

Slope of Sw

Slope of K

Slope of Pc

Slope of r

Slope of BVW

1
2
3
4
5
6
7
8
9

IG 30-1
IJ 30-1
IG 33-1
IH 35-2
IG 34-1
IF 32-1
IF 34-2
IE 34-6
IE 36-1

0.019
0.016
0.025
0.021
0.014
0.032
0.025
0.040
0.025

0.63
1.00
1.00
1.00
1.00
1.00
0.70
1.00
1.00

1.79
1.97
1.82
2.00
2.26
1.79
1.66
1.83
1.70

0.56
0.51
0.55
0.50
0.45
0.56
0.60
0.55
0.59

7.16
7.88
6.88
8.0
9.04
7.16
6.64
7.30
6.80

3.24
3.60
3.10
3.60
4.10
3.20
3.00
3.30
3.10

3.24
3.60
3.10
3.60
4.10
3.20
3.00
3.30
3.10

0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0

48

A.A. El-Khadragy et al.

From Eq. (10), by plotting Rt vs. A on loglog coordinates,

the relation should result in a straight line with a slope equal to
(3n  m) with a constant (aRw) and permeability. Since
m = n, the slope equals (4m). The intercept of the straight
line with 100% porosity gives aRw (250/K1/2)n.
Aguilera [10] used Eq. (10) to construct parallel lines of
constant permeability. From Eq. (10), the slope of
(3n  m) was determined.
The line of constant permeability equal to (1 md) can be
drawn by calculating Rt from Eq. (9), using a Rw and porosity
for each well as shown in Fig. 4 for the G Member. The permeability ranges from 0.01 to 100 md, then by determining a
certain value of porosity, the corresponding value of Rt can
be calculated for each value of permeability.

where: porosity and water saturation are expressed as fractions

and


2:25
=0:929
14
Pc S0:8
wi  Ut
By solving Eq. (14) for Swi and inserting it into Eq. (4), this
leads to the following equation:
2:25
Rt Um

t aRw 0:0929Pc=Ut

1=0:8 n

or Rt Um2:8125n
aRw 1:0961Pc1:25 
t

n

16

Taking the logarithm of both sides of the proceeding

equation, this leads to the following form:
n

log Rt m 2:8125 n log Ut logaRw 1:0961 Pc1:25 

17

2.2. Capillary pressure

Capillary pressure is the difference in pressure between two
immiscible uids across a curved interface at equilibrium.
The following section presents the involvement of capillary
pressure (Pc) into Picketts plots by using an empirical relationship. Kwon and Pickett [12] published an empirical relationship incorporating the capillary pressure, as follows:
Pc AK=1000Ut B

11

where:
Pc is the mercuryair capillary pressure in pounds per
square inch,
B is approximately equal to 0.45,
A is a constant ranging between 151.35 and 22.91.
They based their correlation on the loglog cross plots of A
vs. water saturation and found the following equation:
A 19:5S1:7
w

12

where: Sw is the water saturation (fraction).

By inserting Eqs. (6) and (12) into Eq. (11) [5], this leads to
the following equation:
2
5 1:7
Pc 7:94319:5S1:7
w Ut Swi =250

0:45

13

Eq. (17) indicates that, a crossplot of Rt vs. Ut on a loglog

paper results in a straight line with a slope equal to
(m + 2.8125 n) with a constant (aRw) and constant capillary
pressure. The straight line cuts the 100% porosity axis at a
point equal to [aRw (1.0961 Pc  1.25)n] on the resistivity
scale. Since m = n, the slope will be equal to the values as
shown in Table 1.
Fig. 5 for the G Member exhibits a Picketts plot including
parallel lines of constant capillary pressures that range from 10
to 800 psi. First, a line of constant capillary pressure equal to
100 psi is drawn by assuming a value for t and calculating the
corresponding value of Rt from Eq. (16). The point with the
values of Rt and t is plotted and from it, a straight line is
drawn through this point with a positive slope. This straight
line corresponds to a constant capillary pressure of 100 psi.
The same procedure is followed for the other capillary
pressures of interest.
2.3. Pore throat aperture radii
Washburn [13] presented
capillary pressure (Pc in
mercury (r in dyne/cm),
and the pore aperture
equation is:

a technique for relating the mercury

dyne/cm2) to the surface tension of
contact angle of mercury in air (U)
radius (r in lm). The Washburn

Porosity

100%

90%

50%

20%

100
0.1

10
1md
0.1
0.01md
aRw=0.025

m=n=1.7

0.01
0.01

Figure 4

15

0.1

10

Rt (ohm.m)

100

Picketts plot incorporating formation permeability for G Member of Abu Roash Formation in IE 36-1 well.

Picketts plot in determining the reservoir characteristics

49
10

Porosity

100%

90%

Pc(psi)

25

50

20%

50%

100
200
400
100
800

0.1

10
1md
0.1
0.01md

aRw=0.025

m=n=1.7

0.01
0.01

Figure 5
well.

0.1

10

100

Rt (ohm.m)

Picketts plot incorporating formation permeability and capillary pressure for G Member of Abu Roash Formation in IE 36-1

Pc 2rcos h=r

18

or r 2rcos h=Pc

19

Extrapolation of the constant r straight line to 100% porosity yields the product {aRw[1.0961(2(0.147)r(cosh/r)1.25)]n}
in the resistivity scale.
If r = 480 dyne/cm and b = 140, Eq. (21) can be reduced
to the form:
n
o
n
Log Rt m 2:8125 nlog Ut aRw 1:0961108:1=r1:25 

If Pc is in pounds per square inch, the pore throat radius in

microns can be calculated from the equation:
r 20:147rcos h=Pc

20

22

Inserting Eq. (20) into Eq. (17) and taking the logarithm of
both sides [5] the equation will be:

A practical pore classication system has been represented

by [14]. The porosity is classied by the geometry of the pores
and by pore size. Pore geometry is classied as intergranular,
intercrystalline, vuggy and fractured, while pore size comprises
megaports (>10 lm), macroports (102 lm), mesoports
(20.5 lm), microports (0.50.25 lm) and nano (<0.25 lm).
Fig. 6 for the G Member shows Picketts plot discussed
previously, including parallel lines of constant pore throat
aperture radii that range from 4 (macro) to 0.25 (nano) in

Log Rt m 2:8125 n log Ut faRw 1:096120:147

 rcos h=r1:25 n g

21

Eq. (21) indicates that, a crossplot of log Rt vs. log Ut results in a straight line with a slope equal to (m + 2.812 n)
for intervals with constant aRw, r, U and pore throat aperture
radius r.

MACRO

8 Aperture(Mm) 4
100%

90%

50%

20%

MESO

Porosity

MICRO

0.5

0.1

0.25
NANO

100
10
1md
0.1
0.01md

aRw=0.025

m=n=1.7

0.01
0.01

0.1

10

Rt (ohm.m)

100

Figure 6 Picketts plot incorporating formation permeability and pore throat aperture for G Member of Abu Roash Formation in IE
36-1 well.

50

A.A. El-Khadragy et al.

lm. The calculated slopes from Eq. (21) are recorded in

Table 1.
First, the line of constant pore throat aperture equal to
(1 lm) is drawn, where Rt is calculated by assuming a value
for t and with a constant value of aRw for each member. By
the values of Rt and Ut, a point is plotted and from it a line
is drawn with the calculated slope. The same process is applied
for the other pore radii of interest. The gure includes pore size
which is ranging between macropores and nanopores.

The lines of constant height above the free water table are
drawn through the calculation of Rt from Eq. (24) as follows:
m2:8125n

Rt Ut

faRw 1:09610:433qw  qh

hr cos h=rh cos hh 1:25 gn

26

Using aRw, qw  qh, rh, hh, r, h and any porosity for these
data, Eq. (26) reduces to the form:
m2:8125n

Rt Ut

2.4. Height above the free-water table

aRw 1:09611:4182 h1:25 

n

27

2.5. Bulk volume of water

Calculation of the height above the free water table is done

through the following equation:
Pc 0:433qw  qh hr cos h=rh cos hh 

The product of a formations water saturation (Sw) and its

porosity (U) is the bulk volume of water (BVW).

23

BVW Sw  Ut

where:

28

where:

rh is the brinehydrocarbon surface tension.

hh is the contact angle of the hydrocarbon in water.

BVW = bulk volume of water,

Sw = water saturation of uninvaded zone (Archie equation), and
Ut = total porosity.

Inserting Eq. (23) into Eq. (17) [5] leads to the form:
log Rt m 2:8125 n log Ut
logfaRw 1:09610:433qw  qh
hr cos h=rh cos hh 1:25 n g

Sanyal and Ellithorpe [3] showed that, a crossplot of porosity

vs. resistivity on loglog coordinates should result in a straight
line with a slope equal to (nm). The concept has been also documented by [15], who have termed the product (UtSw) bulk volume of water (BVW). This product is also known as the Buckles
number [16]. The concept has been extended by [17].
Fig. 8 for the G Member exhibits a Picketts plot with
(BVW). Since m equals n, the slope of (nm) is equal to zero,
and the lines of constant BVW will be vertical lines.
The deduced values of reservoir petrophysical parameters
for the nine wells are illustrated in Table 1.

24

From Eq. (24), the crossplot of Rt vs. Ut on loglog

coordinates results in a straight line with a slope equal to
(m + 2.8125 n) with constant aRw, qw, qh and height (h)
above the free water table. The straight line cuts the coordinate
of t of 100% at a point equal to:
aRw 1:09610:433qw  qh
hr cos h=rh cos hh 1:25 n 

25

Fig. 7 for the G Member shows Picketts plot presented previously, showing parallel lines of constant height above the
free water table that ranges from 50 to 500 ft. Different slopes
are calculated and recorded in Table 1.

3. Summary and conclusion

A new technique has been applied through Picketts plot, to
develop some of reservoir petrophysical parameters. These

25 height(ft)

Porosity

100%

90%

50%

20%

50
100
200
400
100

0.1

10
1md
0.1
0.01md

aRw=0.025

m=n=1.7

0.01
0.01

0.1

10 Rt (ohm.m)

100

Figure 7 Picketts plot incorporating formation permeability and high above free water table for G Member of Abu Roash Formation in
IE 36-1 well.

Picketts plot in determining the reservoir characteristics

0.20

Porosity

100%

90%

51
0.10 0.06

50%

0.04

0.02 BVW

20%

100
0.1

10
1md
0.1
0.01md

aRw=0.025

m=n=1.7

0.01
0.01

0.1

10

Rt (ohm.m)

100

Figure 8 Picketts plot incorporating formation permeability and bulk volume of water for G Member of Abu Roash Formation in IE
36-1 well.

parameters include capillary pressure, pore throat aperture

radii, height above the free water table and bulk volume of
water. This technique depends on the use of loglog plots of
effective porosity versus resistivity combined with empirical
relationships for calculating the capillary pressure expressed
as a function of permeability, porosity and water saturation.
The result plot is a straight line, where the slope of this line will
be controlled by the value of m. The intersection of the line passing through the data point and the 100% porosity will be a Rw.
The integration of these petrophysical parameters on a log
log graph of porosity versus resistivity gives the importance for
Pickett plot to be used in reservoir interpretation.
References
[1] A.J. Martin, S.T. Solomon, D.J. Hartmann, Characterization of
petrophysical ow units in carbonate reservoirs, AAPG Bull. 81
(5) (1997) 734759.
[2] R. Aguilera, M.S. Aguilera, The integration of capillary pressures
and Pickett plots for determination of ow units and reservoir
containers, SPE Reservoir Eval. Eng. 5 (6) (2002) 465471.
[3] S.K. Sanyal, J.E. Ellithorpe, A generalized resistivityporosity
crossplot concept, Society of petroleum engineers California.
Regional meeting, SPE paper 7145, 1978, 8pp.
[4] G.E. Greengold, The graphical representation of bulk volume of
water on the Pickett crossplot, The Log Analyst 27 (3) (1986) 125.
[5] R. Aguilera, Incorporating capillary pressure, pore throat
aperture radii, height above free-water table and Winland r35
values on Pickett plots, AAPG Bull. 86 (4) (2001) 605620.
[6] G.E. Archie, The electrical resistivity logs as an aid in
determining some reservoir characteristics, Trans. AIME 146
(1942) 5467.

[7] G.R. Pickett, A review of current techniques for determination

of water saturation from log, J. Petrol. Technol. 18 (1966) 1425
1433.
[8] G.R. Pickett, Pattern recognition as a means of formation
evaluation, in: Society of Professional Well Log Analysts 14th
Annual Logging Symposium Transactions Paper A, 1973, pp.
A1A21.
[9] R.L. Morris, W.P. Biggs, Using log-derived values of water
saturation and porosity, in: Society of Professional Well Log
Analysts Annual Logging Symposium, 1967, 26pp.
[10] R. Aguilera, A new approach for analysis of the nuclear
magnetic log, resistivity log combination, J. Can. Pet. Technol.
29 (1) (1990) 6771.
[11] R. Aguilera, Naturally Fractured Reservoirs, PennWell Books,
Tulsa, Oklahoma, 1995.
[12] B.S. Kwon, G.R. Pickett, A new pore structure model and pore
structure interrelationship, in: Society of Professional Well Log
Analysts, 16th Annual Logging Symposium, 1975, 14pp.
[13] E.W. Washburn, Note on a method of determining the
distribution of pore sizes in a porous material, Proc. Natl.
Acad. Sci. 7 (4) (1921) 115116.
[14] A.J. Martin, S.T. Solomon, D.J. Hartmann, Characterization of
petrophysical ow units in carbonate reservoirs, AAPG Bull. 83
(7) (1999) 11641173.
[15] J.H. Doveton, W.J. Guy, W.L. Watney, G.C. Bohling, S. Ullah,
D. Adkins-Heljeson, Log analysis of petrofacies and ow units
with microcomputer spreadsheet software; Kansas Geological
Survey, University of Kansas, Lawrence, Kansas, 1996.
[16] R.S. Buckles, Correlating and averaging connate water
saturation data, J. Can. Petrol. Technol. 5 (1965) 4252.
[17] R. Aguilera, Extension of Pickett plots for the analysis of
shaly formation by well logs, The Log Analyst 31 (6)
(1990) 304313.