Application of the Kozeny-Carman Equation to

2 downloads 0 Views 387KB Size Report
Oct 29, 2001 - To cite this article: Hilmi S. Salem (2001) Application of the Kozeny-Carman Equation to. Permeability Determination for a Glacial Outwash ...
Energy Sources

ISSN: 0090-8312 (Print) 1521-0510 (Online) Journal homepage: http://www.tandfonline.com/loi/ueso19

Application of the Kozeny-Carman Equation to Permeability Determination for a Glacial Outwash Aquifer, Using Grain-size Analysis Hilmi S. Salem To cite this article: Hilmi S. Salem (2001) Application of the Kozeny-Carman Equation to Permeability Determination for a Glacial Outwash Aquifer, Using Grain-size Analysis, Energy Sources, 23:5, 461-473, DOI: 10.1080/009083101300058480 To link to this article: http://dx.doi.org/10.1080/009083101300058480

Published online: 29 Oct 2001.

Submit your article to this journal

Article views: 119

View related articles

Citing articles: 5 View citing articles

Full Terms & Conditions of access and use can be found at http://www.tandfonline.com/action/journalInformation?journalCode=ueso20 Download by: [Professor Prof. Hilmi S. Salem]

Date: 21 April 2017, At: 00:52

Energy Sources, 23:461 ± 473, 2001 Copyright # 2001 Taylor & Francis 0090-8312 /01 $12.00 ‡ .00

Application of the Kozeny-Carman Equation to Permeability Determination for a Glacial Outwash Aquifer, Using Grain-size Analysis HILMI S. SALEM Atlantic Geo-Technology Halifax, Nova Scotia, Canada Permeability (hydraulic conductivity) is a measure of the ability of a porous medium to convey a single Xuid or multiphase Xuids through its body under a certain hydraulic gradient. In this study, the permeability was determined for the grain-size fractions and the layers of a glacial outwash aquifer and its aeration zone (northern Germany) using grain-size analysis and the Kozeny-Carman equation. The aquifer and the aeration zone are composed of unconsolidated sediments characterized by lateral and vertical heterogeneities. The sediments consist, primarily, of silts, sands and gravels, with a majority of sands and a variety of grain sizes. The permeability of the aquifer and the aeration zone ranges from 2:85 £ 10 ¡4 to 4:60 £ 10¡1 m/s, corresponding to the grain-size fractions, and from 4:02 £ 10¡6 to 3:35 £ 10¡2 m/s, corresponding to the layers. These ranges of permeability agree well with those obtained experimentally from pump tests and laboratory measurements. Several empirical equations correlating between various parameters, with coeYcients of correlation ranging from 0.65 to 1.0, were obtained. Keywords glacial aquifers, grain-size distribution, Kozeny-Carman equation, permeability, porosity, speci® c surface area

Numerous fresh-water aquifers of glacial origin are present in diŒerent regions of the world, such as the northern central part of the United States of America, the southern part of Canada, and the northern part of Europe. Glacial aquifers generally consist of silts, sands, and gravels, with variable amounts of clays. The sediments of glacial aquifers have a wide range of grain size, grain shape, and grain type (mineralogy). Large quantities of the sands and gravels were deposited as outwash material, swept out from the melting glaciers by the forces of the melt-water streams and moraines in front of the glaciers. In their studies on glacial aquifers, McDonald and Wantland (1961), Frohlich (1973), Kelly (1977), Urish (1981) , and Frohlich and Kelly (1985) pointed out that glacial aquifers are greatly heterogeneous, laterally and vertically. The area of investigation, Segeberger forest, forms about 10% of the total area of the province of Schleswig-Holstein (S-H), northern Germany. The sediments of the area are of Pleistocene age, similar to all glacial deposits in northern Europe (Einsele and Schulz, 1973). The altitude of the area ranges from 25 to 65 m above sea Received 31 March 2000; accepted 11 May 2000. Address correspondence to Hilmi S. Salem, Atlantic Geo-Technology, 26 Alton Drive, Suite 307, Halifax, Nova Scotia, B3N 1L9 Canada. E-mail: hilmisalem@ canada.com

461

462

H. S. Salem

level. The annual precipitation on the area is about 800 mm/yr, distributed as follows: 500 mm evaporation, 60 mm runoŒ, and 240 mm in® ltration recharging the aquifer system. The aquifer ranges in thickness from 30 to 70 m, with a water table about 5 to 10 m deep and about 20 m underneath the areas of higher altitude. The aquifer is underlain by an aquiclude composed of boulder clays, known as ``Geschiebemergel.’’ Fluid ¯ ow through porous media may be quantitativel y evaluated by knowledge of several physical parameters and conditions, as well as several lithological attributes. These include velocity of ¯ ow; gravity acceleration; pressure drop; hydraulic gradient; permeability; porosity; speci® c surface area; density, viscosity, and temperature of ¯ uid; pore-water salinity; geometry of pores and pore channels (tortuosity); shape and size of pores and pore channels; and size, shape, type, packing, sorting, distribution, and orientation of grains (anisotropy). In the present study, emphasis is placed on determination of permeability, with respect to the porosity and the speci® c surface area, for the grain-size fractions and the layers of the S-H aquifer and the aeration zone, using grain-size analysis and the Kozeny-Carman equation.

Theory For permeability …k† determination, several equations have been theoretically, experimentally, and empirically developed (e.g., Slichter, 1899; Van Terzaghi, 1923; Krumbein and Monk, 1942; Wyllie and Rose, 1950; Loudon, 1952; Rose, 1957; Harr, 1962; Chilingar et al., 1963; Davis and De Wiest, 1966; Berg, 1970; Todd, 1980; De Marsily, 1986; Katz and Thompson, 1986; Johnson et al., 1987; Salem, 1994; Salem and Chilingarian, 1999). Those equations are dependent, one way or another, on three fundamental equations: Darcy’s equation (Darcy’s law), Poiseuille’s equation, and the Kozeny-Carman equation. The temporal rate of change of the position of a ¯ uid particle is represented by means of the particle velocity. Fluid velocity is one of the major factors in the mechanism of ¯ uid ¯ ow through a medium because it is aŒected by the dynamic and kinematic relationships among the ¯ uid, medium and ¯ ow. The apparent velocity of ¯ uid ¯ ow …¾ap †, in cm/s, also known as ® lter (microscopic) velocity, is de® ned as the mean velocity of a ¯ uid, ¯ owing with a rate q, in cm3 /s, through the entire space of a porous medium, having a cross-sectional area A, in cm2 , and permeability k, in cm/s, i.e., ¾ap ˆ

q A

…1†

Equation (1) is one of the simplest forms of Darcy’s law. The actual velocity of ¯ uid …¾ac †, in cm/s, also known as distance (macroscopic) velocity is obtained by dividing ¾ap by the fractional porosity …¿† of the medium. For a ¯ uid with viscosity · (in cP) ¯ owing with a rate q (in cm3 /s) under a pressure gradient ¢P (in atm/cm), normal to cross-sectional area A (in cm2 ) of a medium with permeability k (in Darcy), the apparent velocity …¾ap †, in cm/s, is (Darcy’s law): ³ ´ q k …2† ¾ap ˆ ˆ ¢P A ·

The Kozeny-Carman Equation to Permeability Determination

463

Using the same units given in Equation (2), k (in Darcy) for the same medium, having length L (in cm), is kˆ

q·L A¢P

…3†

For a ¯ uid with viscosity · (in Poise) ¯ owing through a number of capillary tubes nc , each with capillary length Lc and capillary radius rc (both in cm), where ¢P is in dyne/cm 2 , the general form of the Poiseuille’s equation for capillary ¯ ow is !³ ´ nc ºr4c ¢P qˆ …4† 8· Lc For linear ¯ ow and capillary ¯ ow, the Darcy’s law and the Poiseuille’ s equation are quite similar (Tiab and Donaldson, 1996). Using the same units given in Equations (2)± (4), q of a ¯ uid passing through capillary tubes with a total cross-sectional area At (in cm 2 ) is ³ ´³ ´ kAt ¢P qˆ …5† Lc · Based on the equations of Darcy and Poiseuille, Kozeny (1927) derived an equation that was later modi® ed by Carman (1937, 1938) and was thus known as the Kozeny-Carman equation. The Kozeny-Carman equation is a good-working hypothesis for determination of permeability or speci® c surface area (Loudon, 1952; Wyllie and Splanger, 1952; Rose, 1957; Kezdi, 1974; De Marsily, 1986; Salem and Chilingarian, 1999). The Kozeny-Carman equation and its derivatives can be expressed as follows: ¾ap ˆ where ¾ap 5 ¢P 5 g 5 ¿ 5 Kcc 5 ss 5 · 5 Lc 5

¢Pg¿3 Kcc s2s ·Lc

…6†

apparent velocity of ¯ uid ¯ ow (cm/s); pressure gradient in the direction of ¯ uid ¯ ow (g/cm2 ); gravity acceleration (ˆ 981 cm/s 2 ); porosity (fractional); Kozeny-Carman coe cient (dimensionless); speci® c surface area (cm ¡1 ); viscosity of ¯ uid (Poise ˆ g/cm.s); length of capillary tubes (cm).

Substituting · by f¸»g, where ¸ is ¯ uid’s kinematic viscosity (in cm2 /s) and » is ¯ uid’s density (in g/cm3 ), and then substituting f¢P=»Lc g by the hydraulic gradient, I, dimensionless, Equation (6) is ¾ap ˆ

Ig¿3 ¸Kcc s2s

…7†

It is important to mention that the symbol I is used for the hydraulic gradient instead of i because i is used later for the grain-size fractions. Substituting f¿ap =Ig by k, then k (in cm/s) is

464

H. S. Salem g¿3 ¸Kcc s2s



…8†

Carman (1937, 1938) substituted ¿3 by f¿3 =…1 ¡ ¿†2 g, dimensionless, known as porosity factor, which is a measure of texture of the medium, relating k (in cm/s) to grain size; consequently, Equation (8) is kˆ

(

¿3 …1 ¡ ¿†2

)» ³

g …¸Kcc s2s

´¼

…9†

For k, in cm 2 , the Kozeny-Carman equation can be simply expressed as (Wyllie and Gregory, 1955; Chilingar et al., 1963; Pfannkuch, 1969; Salem and Chilingarian, 1999) kˆ

¿

…10†

Kcc s2s

The Kozeny-Carman coe cient, Kcc , is de® ned as the product of tortuosity …½† and shape factor …Shf †, both are dimensionless, i.e., fKcc ˆ ½Shf g. The shape factor is a measure of the shape of the grains, pores, and pore channels in a porous medium. The tortuosity and shape factor re¯ ect the geometry of the cross-sectional area of the pore channels normal to the ¯ ow direction. For unconsolidated sediments, Carman (1937, 1938) assigned a value of 2.0 for ½ and a value of 2.5 for Shf , which both results in a value of 5.0 for Kcc . The speci® c surface area, ss , is variably de® ned as the interstitial surface area of the pores and pore channels for each unit of bulk volume, grain volume, pore volume, in cm2 /cm 3 (ˆ cm¡1 ), or for a unit of weight of a material, in cm2 /g. For a medium, with grain diameter D, in cm, and porosity ¿, in fraction, Carman (1937) gave the following equation for ss : ss ˆ

6…1 ¡ ¿† D

…11†

If porosity and grain size are known, then Equation (11) can be used to obtain the speci® c surface area of fractions of sediments …ssf , in cm¡1 †. For a layer consisting of diŒerent fractions (X i to X n), with various grain sizes and various speci® c surface areas, the speci® c surface area of all sediment fractions in that layer (ssl , in cm ¡1 ) can be obtained as (Loudon, 1952) ssl ˆ L

(

n X iˆ1

X i ssi

)

…12†

where L 5 dimensionless coe cient ranging from 1.1 (for rounded grains) to 1.4 (for angular grains); ssi 5 speci® c surface area of fraction i, in cm¡1 (ˆ 6=Dmi †; Dmi 5 mean diameter of the grains in fraction i (in cm).

The Kozeny-Carman Equation to Permeability Determination

465

Methodology The permeability, k, was determined for the grain-size fractions and the layers of the S-H aquifer and its aeration zone with respect to the porosity, ¿, and the speci® c surface area, ss . That was achieved by using grain-size analysis (GSA) of the sediments obtained from 6 wells at depth (Z) ranging from 1 to 32 m and by applying the Kozeny-Carman equation. The fractions of silts, ® ne, medium, and coarse sands and gravels, as well as a small amount of clays were recognized at various depths. The weight percentage (WP) and the grain-size distribution (GSD), both obtained from GSA, were plotted on semi-log paper to construct the accumulative curves that represent the various fractions of the sediments. The accumulative curves were then used to obtain the grain size (Gs , in mm) as D10 , D50 , and D90 , which are, respectively, de® ned as the size of the grains (in terms of diameter D) regressed at WP of 10, 50, and 90% . The uniformity coe cient …Uc †, dimensionless, was obtained as the ratio of D regressed at 90% to D regressed at 10% , i.e., fUc ˆ D90 =D10 g. For each accumulative curve, the values of Uc and D50 (in mm) were introduced in the following empirical equations, given by Urish (1981) for similar heterogeneous glacial aquifers, to obtain the porosity, in % , as maximum …¿max † and minimum …¿min †, from which the average porosity …¿av † was obtained as f…¿max ‡ ¿min †=2g log ¿max ˆ 1:62563 ¡ 0:08653 log D50 ¡ 0:03636 log U c

…13†

log ¿min ˆ 1:53902 ¡ 0:18968 log D50 ¡ 0:08201 log Uc

…14†

The speci® c surface area of the grain-size predominant fractions of the grain sizes of the sediments in each layer (ssf , in cm ¡1 ) was obtained in accord with Equation (11). The speci® c surface area of all the fractions in each layer (ssl , in cm ¡1 ) was obtained in accord with Equation (12). For this purpose, one requires Gs regressed at D10 , D50 , and D90 and the corresponding WP (both read from the accumulative curves), as well as a value for L, which was assumed as 1.3, because of the variations in the grain shape of the sediments. The values of porosity and speci® c surface area were then used in the KozenyCarman equation (Equation 9) to determine the permeability …kf † and …k1 †, in cm/s, for the fractions and the layers, respectively. A value of 0.01 cm 2 /s was used for ¸ and a value of 5.0 (dimensionless) was used for Kcc . The values of D10 , D50 , and D90 (in mm), U c (dimensionless), ¿av (in % ), ssf and ssl (in cm ¡1 ), kf and k1 (in m/s), corresponding to the various layers penetrated by the six wells at diŒerent depths …Z†, are given in Table 1. Empirical equations [Equations (15)± (32)], along with their coe cients of correlation …Rc †, ranging from 0.65 to 1.0, are given in Table 2. These 18 equations represent relationships correlating among the various parameters determined in this study. Examples of the relationships are given in Figures 1± 4.

Results and Discussion The grain-size analysis showed that the sediments exhibit a range of Gs from µ 63 mm (clays and silts) to ¶ 2 £ 104 mm (coarse gravels). The analysis also showed that the amount of clays is very small, agreeing with the results obtained by Schroeter (1983), who indicated that the clay fraction does not exceed 1.5% of the

466

1± 5 77 722

3.5± 4 76.2 77

5.5± 6.5 710.5

1± 6 78 714 722

1± 6 712.5

1± 9 710 719 721 732

B

C

D

E

F

Z (m)

A

WELL

85 20 102 83 95

230 160

16 145 180 130

12 215

15 10 80

100 600 125

D10 (mm)

440 925 390 290 500

600 380

220 1,800 3,500 600

110 800

200 35 440

590 3,200 525

D50 (mm)

6,750 30,000 2,200 775 7,000

3,800 1,400

1,550 4,600 4,400 4,600

575 17,500

525 182 2,600

3,800 5,400 2,700

D90 (mm)

79.4 1,500 21.6 9.3 73.7

16.5 8.8

96.9 31.7 24.4 35.4

47.9 81.4

35.0 18.2 32.5

38.0 9.0 21.6

Uc

33.5 26.0 36.6 39.9 33.1

35.1 38.6

39.5 28.1 25.0 33.7

45.0 29.6

42.3 50.8 35.2

33.3 29.3 35.2

¿av (% )

91 48 98 125 81

65 97

166 24 13 67

300 53

173 844 89

68 13 75

ssf (cm ¡1 †

153 1,133 212 431 148

226 153

1,896 284 173 575

2,878 158

2,285 7,100 190

172 49 218

ssl (cm¡1 )

1.20 £ 10¡2 3.00 £ 10¡2 1.19 £ 10¡2 8.42 £ 10¡3 1.51 £ 10¡2

2.52 £ 10¡2 1.31 £ 10¡2

4.71 £ 10¡3 1.33 £ 10¡1 3.93 £ 10¡1 2.27 £ 10¡2

1.78 £ 10¡3 2.96 £ 10¡3

4.22 £ 10¡3 5.37 £ 10¡5 2.53 £ 10¡3 7.02 £ 10¡4 4.44 £ 10¡3

2.08 £ 10¡3 5.27 £ 10¡3

3.57 £ 10¡5 9.49 £ 10¡4 2.18 £ 10¡3 3.02 £ 10¡4

1.94 £ 10¡5 3.30 £ 10¡3

2.76 £ 10¡5 4.02 £ 10¡6 2.98 £ 10¡3

3.33 £ 10¡3 3.35 £ 10¡2 2.23 £ 10¡3

2.13 £ 10¡2 4.60 £ 10¡1 1.94 £ 10¡2 4.80 £ 10¡3 2.85 £ 10¡4 1.36 £ 10¡2

k1 (m/s)

kf (m/s)

Table 1 Grain size of sediments (diameter of grains; D10 , D50 , and D90 ), in mm, regressed at weight of 10, 50, and 90% ; uniformity coe cient …Uc †, dimensionless; average porosity …¿av †, in % ; speci® c surface area of the predominant grain-size fractions of sediments in a layer …ssf †, in cm¡1 ; speci® c surface area of all grain-size fractions of sediments in a layer …ssl †, in cm¡1 ; permeability of the predominant grain-size fractions of sediments in a layer …kf †, in m/s; and permeability of all grain-size fractions of sediments in a layer …k1 †, in m/s. The parameters were determined from grain-size analysis for samples taken from 6 wells penetrating the aeration zone and the aquifer (S-H, northern Germany) at depth (Z) ranging from 1 to 32 m

The Kozeny-Carman Equation to Permeability Determination

467

Table 2 Empirical equations [Equations (15)± (32); listed as indicated in the text] with coe cients of correlation (Rc , ranging from 0.65 to 1.0) correlating among grain size of sediments (diameter of grains; D10 , D50 , and D90 ), in mm, regressed at weight of 10, 50, and 90% ; average porosity …¿av †, in% ; speci® c surface area of the predominant grain-size fractions of sediments in a layer …ssf †, in cm¡1 ; speci® c surface area of all grain-size fractions of sediments in a layer …ssl †, in cm¡1 ; permeability of grain-size fractions …kf †, in m/s; and permeability of layers …k1 †, in m/s. These equations are based on grain-size analysis for samples obtained from 6 wells at depth (Z) ranging from 1 to 32 m, penetrating the aeration zone and the aquifer (S-H northern Germany). The equations represent grain size ranging from µ 63 mm to ¶ 2 £ 104 mm, which correspond to clays and silts (µ 63), ® ne sands (63± 200), medium sands (200± 630), coarse sands (630± 2000), ® ne gravels (2000± 6300), medium gravels …6:3 £ 103 –2 £ 104 †, and coarse gravels …¶ 2 £ 104 †. Equations (21), (26), (28) and (32) correspond, respectively, to Figures 1± 4 Eq. # (15) (16) (17) (18) (19) (20) (21) (22) (23) (24) (25) (26) (27) (28) (29) (30) (31) (32)

Equation …¡0:095548 † ¿av ˆ 52:525 D10

…¡0:157640 † ¿av ˆ 92:052 D50 …¡0:130970 † ¿av ˆ 97:993 D90 …¡1:0499† ssl ˆ 36299 D10 † ssl ˆ 1:03 £ 105 D…¡0:90457 50 …¡0:63162 † ssl ˆ 57009 D90 …4:9434† ssf ˆ 1:9334 £ 10¡6 ¿av …4:4436† ssl ˆ 5:4557 £ 10 ¡5 ¿av …1:1724† kf ˆ 1:04 £ 10¡4 D10 kf ˆ 1:03 £ 10¡6 D…1:5673† 50 …1:9423† k1 ˆ 1:589 £ 10¡7 D10 …¡0:096014 † ¿av ˆ 23:436 kf † ssf ˆ 10:155 k…¡0:28482 1 …¡0:57796 † ssl ˆ 35:937 kf kf ˆ 32:143 s…¡1:7268† sf …¡1:8536† k1 ˆ 44:284 ssl ssl ˆ 4:9474 s…0:99285† sf …1:0031† k1 ˆ 4:3827 £ 10¡2 kf

Nature of Correlation

Rc

Inverse Correlation

0.65

Inverse Correlation

0.95

Inverse Correlation

0.92

Inverse Correlation

0.90

Inverse Correlation

0.96

Inverse Correlation

0.91

Direct Correlation

0.96

Direct Correlation

0.92

Direct Correlation

0.77

Direct Correlation

0.99

Direct Correlation

0.97

Inverse Correlation

0.92

Inverse Correlation

0.85

Inverse Correlation

0.96

Inverse Correlation

1.00

Inverse Correlation

0.99

Direct Correlation

0.96

Direct Correlation

0.72

Figure

Fig. 1

Fig. 2 Fig. 3

Fig. 4

sediments. The D10 (Table 1) covers a Gs range of 10± 600 mm (median ˆ 100 mm), corresponding to silts, ® ne and medium sands, and the small clay portion. The D50 (Table 1) covers a Gs range of 35± 3500 mm (median ˆ 500 mm), corresponding to silts, ® ne, medium, and coarse sands, and ® ne gravels. The D90 (Table 1) covers a Gs range from 182 to 3 £ 104 mm (median ˆ 3800 mm), corresponding to ® ne, medium, and

468

H. S. Salem

coarse-grained sands and gravels. The median values (100, 500, and 3800 mm) correspond, respectively, to ® ne sands, medium sands, and ® ne gravels. The sediments within the aquifer exhibit a general range of Gs from 100 to 6300 mm, corresponding to ® ne, medium, and coarse sands, as well as ® ne gravels. The Gs range of the aquifer coincides with the median values of Gs (100, 500, 3800). These results indicate that the values of Gs (regressed at D10 , D50 , and D90 ) cover together the whole range of the size of the grains present in the aquifer and the aeration zone. The obtained physical parameters (Table 1) have the following ranges and median values (given in brackets): ¿av ˆ 25± 51% (35.1% ), corresponding to a range of ¿max of 33± 55% and a range of ¿min of 16± 42% ; ssf 5 13± 844 cm¡1 (81 cm¡1 ); ssl ˆ 49± 7100 cm¡1 (218 cm¡1 ); kf ˆ 2:85 £ 10¡4 –4:6 £ 10¡1 m/s (1:15 £ 10¡2 m/s); k1 ˆ 4:02£ 10¡6 –3:35 £ 10¡2 m/s (2:18 £ 10¡3 m/s). The wide ranges of these parameters can be contributed to the high degree of lateral and vertical heterogeneitie s characterizing the sediments of the aquifer and the aeration zone. It is important to mention that the median values are used instead of the mean values, because the median values are more representative of the wide variations of the parameters obtained than the mean values. In addition, a median value (representing the value of the middle variable of a class of data and corresponding to the 50% point on an accumulative frequency range of values) is not sensitive to extreme values, as in the case of the mean value. Regarding the porosity in particular, the diŒerence between the mean value (35.3% ) and the median value (35.1% ) is negligible because ¿ has no extreme values, similar to the other parameters. The interrelationships among the various parameters determined for the fractions and the layers [Table 2; Equations (15)± (32); Figures 1± 4] show the presence of a strong interconnection between Gs (represented by D10 , D50 , and D90 ), ¿, ss , and k. Table 2 shows that an increase in ¿ is associated with a decrease in Gs [Equations (15)± (17)]; an increase in ss is associated with a decrease in Gs (Equations (18)± (20) and with an increase in ¿ [Equations (21)± (22); Figure 1]; and an increase in k is associated with an increase in Gs [Equations (23)± (25)] and with a decrease in ¿ [Equation (26); Figure 2], as well as with a decrease in ss [Equations (27)± (30);

Figure 1. Average porosity …¿av †, in % , versus speci® c surface area of the predominant grainsize fractions of sediments in a layer …ssf †, in cm¡1 ; 19 readings representing 6 wells penetrating the aeration zone and the aquifer (S-H, northern Germany).

The Kozeny-Carman Equation to Permeability Determination

469

Figure 2. Permeability of the predominant grain-size fractions of sediments in a layer …kf †, in m/s, versus average porosity …¿av †, in % ; 19 readings representing 6 wells penetrating the aeration zone and the aquifer (S-H, northern Germany).

Figure 3. Permeability of the predominant grain-size fractions of sediments in a layer …kf †, in m/s, versus speci® c surface area of all grain-size fractions of sediments in a layer …ssl †, in cm¡1 ; 19 readings representing 6 wells penetrating the aeration zone and the aquifer (S-H, northern Germany).

Figure 3]. The speci® c surface area of the predominant grain-size fractions …ssf † and that of all the fractions together in a layer …ssl † show a progressive increase with each other [Equation 31]. Also, the permeability of fractions …kf † and that of layers …k1 † increase progressively with each other [Equation (32); Figure 4]. The fact that ssf < ssl and kf > k1 (Table 1) suggests that an increase in the number of fractions in a layer, particularly those of a ® ner grain size, leads to an increase in the speci® c surface area of the sediments, which results in a lower permeability of that layer. The uniformity coe cient …Uc ˆ D90 =D10 † does not show strong correlations with the various parameters obtained because it re¯ ects a wide range of grain sizes, indicating lack of uniformity of the sediments. The range of Uc [º 9–97, with a single reading of 1500 (median ˆ 32.5); Table 1] agrees well with the ranges of heterogeneous glacial deposits given by other researchers.

470

H. S. Salem

Figure 4. Permeability of the predominant grain-size fractions of sediments in a layer …kf †, in m/s, versus permeability of all grain-size fractions of sediments in a layer …k1 †, in m/s; 19 readings representing 6 wells penetrating the aeration zone and the aquifer (S-H, northern Germany).

The speci® c surface area is a sensitive parameter to the variations of Gs . The results indicate that an increase in the value of ss is related to an increase in the amount of ® ne-grained sediments (associated with an increase in ¿ and a decrease in k). Also, a decrease in the value of ss is related to an increase in the amount of coarse-grained sediments (associated with a decrease in ¿ and an increase in k). It is important to mention that if clays were present in a greater amount, the values of ss would be much higher than those obtained, and thus the values of k would be much lower than those obtained. Davis and De Wiest (1966) presented k values for clays, ranging from 10¡10 to 10¡8 m/s (0.01± 1 mD). For hydrocarbon reservoirs, composed of shaly (clayey) sandstones, Salem and Chilingarian (1999) obtained values of ss of up to º 4:8 £ 104 cm ¡1 …k ˆ 1:9mD†, corresponding to clays …Gs ˆ 1:2 mm†. The variations of ss aŒect strongly the variations of k, because ss contributes to the nature of the grain-to-grai n contact and the size of the pores, which both play a considerable role in the mechanism of ¯ uid ¯ ow through the pore channels (capillaries).

Conclusions Permeability of aquifers can be determined empirically, theoretically or experimentally from laboratory measurements, tracer tests, or pump tests. For the glacial aquifer in S-H (northern Germany), the permeability, along with the speci® c surface area and the porosity, was determined from grain-size analysis using the KozenyCarman equation. The aquifer is composed of unconsolidated, heterogeneous sediments that consist of a variety of grain sizes and diŒerent mineralogical compositions. Variations of the parameters that control the hydraulic ¯ ow through the aquifer are strongly aŒected by the grain-size distribution of the sediments. The nature of the relationships between the various parameters is primarily dependent on the grain-size distribution and variations of the grain size. A decrease in the grain size leads to an increase in the speci® c surface area, which results in a decrease in the

The Kozeny-Carman Equation to Permeability Determination

471

permeability. The high degree of heterogeneity of the aquifer contributes to the wide variations of the physical properties of the sediments, re¯ ected in the water-¯ ow nonuniformity. The results obtained agree well with those obtained experimentally for the same aquifer or other aquifers composed of glacial deposits. The obtained empirical equations, correlating among several parameters, can be successfully applied to similar deposits.

Nomenclature A At D Dm D10 D50 D90 Gs I Kcc L Lc Rc Shf Uc X i¡n X Y Z a g k kf kl nc q rc ss ssi–n ssf ssl L · ¸ ¿ ¿av ¿max ¿min

Cross-sectional area of a porous medium, normal to ¯ uid ¯ ow (cm2 ) Total cross-sectional area of capillary tubes in a porous medium (cm2 ) Diameter of grains (cm, mm, mm) Mean diameter of grains (cm, mm, mm) Size of grains regressed at 10% of weight (cm, mm, mm) Size of grains regressed at 50% of weight (cm, mm, mm) Size of grains regressed at 90% of weight (cm, mm, mm) Grain size (cm, mm, mm) Hydraulic gradient (dimensionless) Kozeny-Carman coe cient (dimensionless) Length of medium (cm) Length of capillary (cm) Correlation coe cient (dimensionless) Shape factor of grains, pores, and pore channels (dimensionless) Uniformity coe cient (dimensionless) Weight percentage of grain-size fractions (i± n) in a layer (% ) Variable of the X-axis in the empirical equations of Figures 1± 4 (various dimensions) Variable of the Y-axis in the empirical equations of Figures 1± 4 (various dimensions) Depth (m) Empirical-equation coe cient in Figure 1 (dimensionless) Gravity acceleration (ˆ 981 cm/s2 ) Permeability (cm/s, m/s, mD, Darcy, cm2 ) Permeability of the predominant grain-size fractions in a layer (m/s) Permeability of a layer (m/s) Number of capillaries in a medium Rate of ¯ uid ¯ ow (cm3 /s) Radius of capillary (cm) Speci® c surface area (cm¡1 , cm 2 /g) Speci® c surface area of grain-size fractions (i± n) in a layer (cm ¡1 ) Speci® c surface area of the predominant grain-size fractions in a layer (cm ¡1 ) Speci® c surface area of all grain-size fractions of sediments in a layer (cm ¡1 ) Coe cient of roundness and angularity of grains (dimensionless) Viscosity of ¯ uid (Centipoise ``cP,’’ Poise, g/cm.s) Kinematic viscosity of ¯ uid (cm2 /s) Porosity (fraction or % ) Average porosity (% ) Maximum porosity (% ) Minimum porosity (% )

472 » ½ ¾ac ¾ap ¢P GSA GSD S-H WP

H. S. Salem Density of ¯ uid (g/m3 ) Tortuosity (dimensionless) Actual velocity of ¯ uid (cm/s) Apparent velocity of ¯ uid (cm/s) Pressure drop (atm/cm, dyne/cm2 , g/cm 2 ) Grain-size analysis Grain-size distribution Schleswig-Holstein Weight percentage (% )

References Berg, R. R. 1970. Method for determining permeability from reservoir rock properties. Trans. Gulf Coast Assoc. Geol. Soc., 20:303± 317. Carman, P. C. 1937. Fluid ¯ ow through granular beds. Trans. Inst. Chem. Engs., 15:150± 156. Carman, P. C. 1938. The determination of the speci® c surface of powders I. J. Soc. Chem. Indus., 57:225± 234. Chilingar, G. V., R. Main, and A. Sinnokrot. 1963. Relationship between porosity, permeability, and surface areas of sediments. J. Sed. Petrol., 33:759± 765. Davis, S. N., and R. J. M. De Wiest. 1966. Hydrogeology, New York: John Wiley & Sons, 463 pp. De Marsily, G. 1986. Quantitative Hydrogeology: Groundwate r Hydrology for Engineers, London: Academic Press, 440 pp. Einsele, G., and H. D. Schulz. 1973. Ueber den Grundwasserhaushalt im norddeutschen Flachland. Teil I: Grundwasserneubildung bewaldeter und waldfreier Sander¯ aechen Schleswig-Holstein. Bes Mitt. Dt. Gewaesserkdl. Jb., 36, 72 pp. Frohlich, R. K. 1973. Detection of fresh water aquifers in the glacial deposits of Northwestern Missouri by geoelectrical method. Bull. Water Resour. 9:723± 734. Frohlich, R. K., and W. E. Kelly. 1985. The relation between hydraulic transmissivity and transverse resistance in a complicated aquifer of glacial outwash deposits. J. Hydrol. 79:215± 229. Harr, M. E. 1962. Groundwater and Seepage, New York: McGraw-Hill, 315 pp. Johnson, D. L., J. Koplik, and R. Dashen. 1987. Theory of dynamic permeability in ¯ uidsaturated porous media. J. Fluid Mech. 176:379± 402. Katz, A. J., and A. H. Thompson. 1986. Quantitative prediction of permeability in porous rocks. Phys. Rev. B. 34:8179± 8181. Kelly, W. E. 1977. Geoelectric sounding for estimating aquifer hydraulic conductivity. Ground Water, 15:420± 425. Kezdi, A. 1974. Handbook of Soil Mechanics, Vol. 1, Amsterdam: Elsevier, 294 pp. Kozeny, J. 1927. Ueber kapillare Leitung des Wassers im Boden. Sitzungsber. Akad. Wiss., Wien, 136a, 271± 306. Krumbein, W. C., and G. D. Monk. 1942. Permeability as a function of the size parameters of unconsolidated sand. Trans. AIME 151:153± 163. Loudon, A. G. 1952. The computation of permeability from simple soil tests. Geotechnique, 3:165± 182. McDonald, H. R., and D. Wantland. 1961. Geophysical procedures in groundwater study. Trans. Am. Soc. Civ. Engs. 126:122± 135. Pfannkuch, H. O. 1969. On the correlation of electrical conductivity properties of porous systems with viscous ¯ ow transport coe cients. Intern. Assoc. Hydrau. Res., 1st Intern. Symp. Found. Trans. Phen. Porous Media 23:496± 505. Rose, W. D. 1957. Fluid ¯ ow in petroleum reservoirs. IÐ The Kozeny paradox. Circ. Ill. St. Geol. Surv., 236:1± 8.

The Kozeny-Carman Equation to Permeability Determination

473

Salem, H. S. 1994. The electric and hydraulic anisotropic behavior of the Jeanne d’Arc Basin reservoirs. J. Pet. Sci. Eng. 12:49± 66. Salem, H. S., and G. V. Chilingarian. 1999. Determination of speci® c surface area and mean grain size from well log-data and their in¯ uence on the physical behavior of oŒshore reservoirs. J. Pet. Sci. Eng. 22:241± 252. Schroeter, J. 1983. Der EinXuss von Textur- und Struktureigenschaften poroeser Medien auf die Dispersivitaet. Ph.D. dissertation, Kiel University, F.R. Germany, 232 pp. Slichter, S. C. 1899. Theoretical investigation of the motion of groundwater, U.S. Geol. Surv., NTIS, Spring® eld, VA, 19th Annual Report Part II, pp. 295± 384. Tiab, D., and E. C. Donaldson. 1996. Petrophysics. Theory and Practice of Measuring Reservoir Rock and Fluid Transport Properties, Houston, TX: Gulf Pub. 706 pp. Todd, D. K. 1980, Groundwate r Hydrology, New York: John Wiley and Sons. 535 pp. Urish, D. W. 1981. Electric resistivity-hydraulic conductivity relationships in glacial outwash aquifers. Water Resour. Res. 17:1401± 1408. Van Terzaghi, K. 1923. Die Berechnung der DurchlaessigkeitsziŒer des Tones aus dem Verlauf der hydrodynamischen Spannungserscheinungen. Sitzungsber. Akad. Wiss., Wien, 132:105± 113. Wyllie, M. R. J., and A. R. Gregory. 1955. Fluid ¯ ow through unconsolidated porous aggregates, eŒects of porosity and particle shape on Kozeny-Carman constant. Indus. Eng. Chem. 47:1379± 1388. Wyllie, M. R. J., and W. D. Rose. 1950. Application of the Kozeny equation to consolidated porous media. Nature 165 (June):972. Wyllie, M. R. J., and M. B. Splanger. 1952. Application of electrical resistivity measurements to problem of ¯ uid ¯ ow in porous media. Bull. Am. Assoc. Pet. Geol. 36:359± 403.