Clase 10 PDF
Clase 10 PDF
• Properties
• Darcy's Law
Úo~/e~;/r;
..{lÁ- Ley ck .¿Jd rC ;¡
(1 v/ ~~ech,*f
/I¡{J~'eú- dz 4~ L¡ de ~6/(1 . '.
~6/er~~e(A.;
!Ie/e/'or'k('~dM/ 4~i:;,;to/Í~
de~~ Le .&~M /e tX;/?'7'e,,¡(
GeJ<.e4j¿>.J!-r~ Le 4 '7/€ J<3rC¡
Darcy ~\rl'.l·("11
l·¡~ . 3.
tIc,t,n"::' a,:¡~",,,:,,,,,",, l. 1". Darcy's Experiments
• Discharge is
i u ~ 1I
Proportional to
P,ly ~. ~ I h+e _ -
I
- Area J. - . L~
- Head difference
proportionality is Datum
Q __
2
-thM -e
le
------- ¡- ~hlt?JlCltv
, 'l'.:. ,,."jrm..1l1LL
~¡~·;;':~;:;;;:::=::: QOCA h¡-h2 h -ht ) O=-KA f..h
T rr
; , ..,,'l...."' . , . ~~,.'j< - ~
L O=-KA T \. c:lédre~
aN~~
http://biosystems.okstate.edu/Darcy/Englishlindex.htm
~h~o~~:
Q ~ - kll~ -- ¿uy'ke/
L - ~I Coef'Ct'eu/e '&rfo,ee~~'~c/~ :
k= (!~Kketv~~/ ¿~.~
~
2;
Darcy's Data
> "
35 o Set 1, Series I
Hydraulic Conductivity
~ 30 o Set 1, Series 2 <> • Has dimensions of velocity [LIT]
~ ~25 t::.. Set 1, Series 3 • A combined property of the medium and the fluid
.:: • Ease with which fluid moves through the medium
~ ! 20 O Set 1, Series 4
, o :.: Set 2
K =k~
~ ~15 ¡¡.
5 P = density
J1 = dynamic viscosity } Fluid properties
o r = specific weight
o 5 lO 15 20
Gro-/;:e~~!~~)
Gradient (m/m)
Q
q =-= - K
A
dh
dl
~ &~Aa¡;j/¿~cI A{/~/~
¡¿eKe ~ ~·~u~o~ ~ ~
Aa ¡;; s- de jJ dr7 ve,~;h¡ [~]
Hydraulic Conductivity
Rocks I
deposlts
Uoconsolidaled
1< k K K K
Groundwater Velocity
.......- - - - . - - ' - - -, (darey) (cm 2) (cmIS) (mis) (gaUdayIll2)
10 5 10 J 1 ' 02
[ I ,-lO"
t- • q - Specific discharge
~~tA
, 104 10 - 4 10 10 - 1
I I ¡' I 2
Discharge from a unit cross
.¡;pt ~¡d:
g ti 10" - 10-' I
1 ' 10 10"
O , >
;;; I I ¡;; ~ " 10" - 10- 6 1"' 10-' 10-'
section area of aquifer
I
~ ~ 'O ~ ~
I
I ~
103
~
ª
in 2 i5
~ ~-5
~ g~ ~
!l; ~
a I
G
~ 1
10 10- 7
10- e
1 10 -2
¡
10- 4
102
direction of flow.
~
10- 3 10- 5
É ~E § ~
~
10
o'sEro: 1?n:
~ 'O n:
I ~
'G.Q
~ I
~
~
I
I
"'''
rx: e I -~ -
.~ ,I I;""'10-3 10-11 t- l 0-6 10-6
Average velocity of fluid
:5 ~ ~
~ li
~
;;
g
I
¡ I
r 10'-4 10-"
.
i-- 1O-' 10 -'
L10
- 2
' flowing per unit cross
~ , 10 ,,4
~ (et&',4.,( M ~ .0~/'~~
c"g
:;¡ ro
r I !
! , -5
1 '- 10 8 10- 11,5 l.- 1O- ~1 10--1:3 - 10
10-'
'VGO
I ::::> , / ~
-4
(dar;) {eJ«~ (eis) (lU-IsJ !fNrKÁ~J"- ~ /t~cC/D>A N / r
rro>o
dh=(h 2 - h¡)=(1O m-12 m) =-2 m \,.)
q = -K-i-K-J-K-k
ax ay az
q= -KJ= -(lxIO- s mi' .) x (-0.02 mlm) = 2xlO- 7 mis
ah, K ah,
= - K ax l; =-
º = qA = (2x 10- 7 mis) x 50 m2 = IxlO- 5 m3/s
qx qy ay'
v=q/;=2xI0- 7 m1s/0.3 = 6.6xI0- 7 m/s Specific discharge vector points in the opposite directiqn of h
b/'O-~v---fo /¡{~~
re¡d ,
-- ~Lh<.~/;"Vc~ L~jú>-- ):;: / X 10 -s 'J<-<j
f
----------------------------------
~ ve;0~'~ v' - l / ~4 t2yec.e}t-'eo fJ-? rdc 4'/;;1'rJ'~
'ó- ,& ~ 'yek fo;-fr~r' .'
J<-iGt-€Xú $ 11! 0 =: - K 7{e'} ~ -(/x 1()-S~)xft.i2~)::
1:: V.o :; .2 >( /0 . 72 ¡f
~ (10 - /~) ::: -
dA :=: Ih<. -~/) 2~
<?-<--<
~ i.:~~r ~'e~ . _J
• We ignored kinetic energy (Iow velocity)
Circular hydraulic ~
..... . ........... ~h j ,i¿~~ f- - • We assumed laminar flow
K, conductlvlty, ¡ ! / qx ,\ . ax \ N R = pqdlO
tles~T4ui~)
, A
dlO = effective grain size diameter
CV(,~CA. K fZ1 ' \'" ""'" /~ J~ • Darcy's Law is valid for N R < 1 (maybe up to 10)
~ CD- ......,.... ......- ----......... ......./. h 1 < h 2 < h3
"~ Darcy
Law
/e~c-fo- eqb
.., L.-/~,a jJfDf
the grain size (d)
2(
k = cd = ]2
<jJ3 2 d Kozeny - Carman eq.
predicts this 180(1- <jJ)
r¡¿~ eslo k :: /er ~t:L~¡;c/'.;p-f 1'x.0'-xr:e~
f~ ~K'dX~~
f-
~
tan-1 (a)= (1/K) q • So how do we get the parameters we need for this
equation?
- We can produce the flow pattern experimentally • Consta~t head difference (h) is rx.~ .
• Hydraulic Conductivity is measured in the lab with a applied across the sample {I(JJo'c;e. .' ~re~{e
/ r~
permeameter producing a flow rate Q
A~b
- 5teady or unsteady 1-D flow
• Darcy's Law ¡:¡
- Small cylindrical sample of medium L
IQ=KAj IK= QLI
Ah
Sample P.
~ ?"Uue;l~
L--J
2 ~
IOtube = 7r1"tube {1t
I I
Initial head
ft1 1I
~ ~
- Properties same at every
point
Heterogeneous
!~um:~~umKB ~ ~_"'"
Properties different at every
_11 bo point
Isotropic
- Properties same in every
¡fuJ¡ 1)l1J¡ i'i--->
direction
Anisotropic
Ih=ho@t=O; h=ht@t=tl - Properties different in different
directions
,yeLe.C-e'~ 4 ~ ~~~~J~iv/~oI
IkIe~r%«'c4/ l /'~¿"OP;~~
;t~r1A4!Á; A ;te/c;",- JestPx~1e
J¡(~Je~o
el ~ /C'a<..Ld'o-f) A.t d C~ I ~ ./<'.r
_ //V~~ /~ "'v- ~~
L~r~~~~~
- /frfo/()?e~~
.--: A~~ ,úfY-e~4 <"u ~~
- "!50~¡-o¡k.~
~ fe~)(.dt:L I ~ d'f~fo-~r- 4¡h ... hof<-eh4 ¿~Á ~ ~ '¿;;;eCd~
~ e~Ñvf"~'~~~ ~íe e.. __ ,,(~'<;o -0¡<~
fL./' -1 ,
~~"'l-'\....T~~·tfTV'...
í
9
Layered Porous Media ~CI~~
Layered Porous Media (Flow Perpendicular to Layers)
___ ____ L
/7 !,er
':iezometric surface
(Flow Parallel to Layers) ------f i h: i
I ~~~~~ t f.h2' ~h
3 3
Q) (L¡)
¡ :~ 1~~3~ ~
Q= I
¡=1
(Q¡) = I(
¡=1
-b¡K¡ ~~) tlh¡ = (b K¡
¡ - - - - - . ' - - - - 1 -1
~ . ,., ". .~ I
[
N N
_=
K ¿31=1 (b ¡K¡) L
K= L _
1 L ' .....; .1;$[ ,',,, 1.', ,..',:',1
b lo( W )1 N (--1.)
K¡ I( )( )( )1
~kd-.o
¿1 = 1 Ll L2 L3
( )
L
k:: ~~oi'Vt:kR ~~
ID
to determine ~
61 .2 Groundwater
groundwater flow contours
direction J '1
r¡<-; vele" fet>J I';{
~~. ~\~ d ~y.e/ » ~~ Direction of
_J
,eeJ ex-r ¡tdf~;-&.
h.r4> ¡Ú.es frP / /o~s/ /ve",G}( groundwater flow
.f 60.8
~ Á\reclC4'~ hP ~ -IR a~
L-" ,
~¡-ecC)(nA.
r ",.f)
~
/
~ h ~r- fc..-b~rn..JAR.cL
kbk~;,d!tA- '
~(
column (cm)
4 ~ h .dore;:
Q - k ¡) [r lit - ft()]; ~~>'<-.@
kJe :
_ Q" ¿:~ 4 ~ 1~4f / ~ ~ .4 ~ t:4<?¿
( tJ,-}U~/~)
gradient
-- / = 2~\
¡S
Transmissivity (T)
-r A ~ ~ 6. LV / ~N
_ 6 ~ o/e~í ¿ d ~~¡&-o I
~ j ~ k (b'W) (c/,{/d f )
- f :: TtU (r/h//f)
_ 6{ ~ Tv (w) (c/~/t/f'l ~,,& .'
~ w~ ~~ ¿()~
~ ~s /e~~ /€C-oJ<-o(k'4 ~ T:.Le ~~
A ~ A-uk ~-M7SI1~ ItJ) ~? ~ é'~~,
(21 ¡2k-h'~ ~ .'
_ (Q ~ iN !c/J/r/~)
dI = 1000 m
I
r'1 8
M - - - _' - ~
-
_ _.-...J
.,'1.:"
·.. ,.. • •• • • • • tt
••
.'.• •
0CONFINED."'
.
.'
· · · O· . · 0 ··
• • O • .' • O• '.
•
•
• %
. • :AQUtFER·· ...... •
., '. b =
. . .
.. ., . 100
m ..-----:~
·. . ." . S o n d ·K· 5' • . .'. . o o. I q' .'
(1 )
7~
_ ik<-~cr;Y¡,/A L~'~ k
k = 5 o ¿fe> : cr 1 'YIA/ )"0 nA.
- eX ~ T 11) (0)
Ji
Ir
Calculation of
Transmissivity
(}J~ h ~ /rd.~){S'lU<,:riJlz'¿c/(T)
_¿ ¡t-j'-Ctxs:?vU-'~'v(¿/~ ~4! .fe¡- MeaA-4- / ,Cc' :
17
Calculation of Transmissivity
Water table
Calculation of Transmissivity
• T = (Q/w) / (dh/dl)
• T =[(5,616 m 3/d) / 5,000 m] / (0.0005 m/m)
• T = 2,246 m2/d
• The saturated hydraulic conductivity (K) can be
back-calculated as:
• K=T/b
• K = (2,246 m 2/d) / 50 m
• K =45 m/d
• Range of transmissivity in principal aquifers:
• 1 m2/d for fractured sedimentary aquifers
• 100,000 m2/d for karst cavernous limestone aquifers
/J ÍJ /J /_ ~ ~ • , . /_p(TI L/_ //
~~ A.€ ~ ./~KS~RV,c~Z. ~ ~
~~IlY¿d~ 4~~~ (k)
_ lo- J?if
~~ ~",-~;r¡- ~ '.
--6 -- So ~
- t;.! ?0-"jO.4 ~/"'-"'-''hA 'x'y~'f'¿~ ~ ts r
'~c...~ ¿Pcu<-'¡C/OS'
- i <=)1,(2 leI ~ at:<A-fe;os- ¡ec&~~~~r faclv/tlhs
_ ho/ OC/O x-t 2¡V ~ acWftttJ.f &lIer]C()HJS ~ M~
~d/1<A-tL1 ¿~'<;S Le 4.. 7~ ry
Darcy's Law Praclice Problems
,,
:,
?- M~¡ák~;~ ~ d r
Jo 1. dlr~ :¿ ~
1. 1
,,
h., L ___ _ _
==--:1
j":{
_L
K>tA
~qOY1. ~:
, .
¿ ~ .T~ --
Á,- A~
~--
3~-1
¡flhe elevalion ofHI is 35m and Ihe elevalion ofh2 is Om, whal is Ihe hydraulic gradienl iflhe
dislance from hIlo h2 is 5.6 km? (Answer in mlkm) L- S.b~
\
rl:.·
3)?7-(
,... ::- ~ . :¿:) I tIkM
Eqllation: Whut \\ ~ kllv\\.
1-(112-11 1I'L h~-h 1~ .'5m L ,.(, km
5tk'"u
:z-I/JY<U' ~ ye,dÚ'./,?{ ,ed' ~ ~~~(/e/Uc..
1 /
Answer: 1=6. 251ll 'klll
2. Find Ihe velocity oflhe waler in Ihe picture below iflhe hydraulic condllclivity is 1141ll/day.
SOllltiol1 '
IlyJr;lldl( gradJ-:nt \11 Ob
- V~~ b.~t¡~/o1J
V=114mday * .06
1- 1;;;
3- l~ JleJ~ 4i' ~ j(d~~
3. Ifthe velocity of groundwater between two wells is 3.4 m per day and the hydraulic gradient
between the two wells is .005, what is the hydraulic conductivity ofthe aquifer between the two wells
~ d'v? t>ojol "'1 3 JI ?U¿;"/ k /) ~
(assuming it is homogenous throughout)?
Sollllion :
!~~r r
4 k?-<A-<>;e~ ~ h, ~ ) J'
K(3 . ~111 ·day) ' 005 ~fr~ '" 4. ¡/¿PU~ ¡¿:: k l ' . K _ Y.ó ~ _ ~.4")o-(/~
K= 6¡Wmday --;>Q., / JI ) / : //_ • ~_ l - O, oóS ~/~
~ ~kcTlY/~ A-<.-A::./~
4. An aquifer is 2045 wide and 28 m tall. lis hydraulic gradient is .05 and its hydraulic conductivity is
145m/day. Calculate the velocity of the groundwater as well as the amount ofwater that passes
Howe\er. 1110st ofthe aquifer is not made up of Clllpty space. only 32'~o IS (porosity of 32 U,u).
,... ~C#IC~Vd;: 7 - ~r~/k
So we mu st multipl y ~15.135 m' 'da)' by .32 to get the amount 01' water that flo"s through the aquifer
él ~ -te a~ r- rfo..- r d ékCt-v- '.fo¡t>
\\'ith a porosiry of 32 0 ' 0 .
Oarcy's Law
d /70~ ¿if Zi~ ~ ~ Ú~~
In 1856 Henri Oarcy, a French engineer, discovered a relationship thal govems fluid flow through
) helZ<A
geologic material s. He detennined that flow through a sand-filled pipe was proportional to the
change in head (the height to which the water rises in a s mall tube) divided by the length ofthe
pipe (dh/dl). K (the proportionality constant) is known as the hydraulie conduetivity and is a
measure ofthe capacity for a porous medium to transmit water. It has high values for sand and
lower values for rock and clay. The relationship of discharge to head, known as Oarey 's Law,
can be written as:
Q= -KiA
{;(-=- !::¿A
Where: Q = the volume of water discharged in a given time periodo (Iength cubed/time; e. g.
~~ e/~~~ A- ¿t~Le~~
ft'/day)
K = the hydraulic conductivity ofthe aquifer matrix. (Iength/time) The negative sign is to show
flow from higher head to lower head.
'" J¡ Á /)drc¡
t :: th ~ ;na ~~ h/ ~'Ye/fte~¡;~
0 ~/pC?~ ~r~;;r~ 101ft/e/e)
1
'"
Figure illustrating Oarey's law.
'/1.= el ~ ~ ~ !eco;' ~~»Y~1rAP
With your groundwater model you ean measure the hydraulie eonductivity ofthe material in the ~ ~Á~C4~ dd~
model and then make and eheek predietions oftravel times in the modeled aquifer.
. A B S"$'J?--fe ,fd'¡10
o-CO
h. l .
; '\ waler lable ?C/
t
re
. (¡d~erMD; /
....d:=rz ~
Measure the gradient (slope) at that diseharge rateo This is the differenee in the level of water in
wells A and B divided by the distanee between those two wells. You could also use lhe
differenee between water levels at the reeharge and diseharge points divided by the distanee
between those two points. Note: the slope will be a negative number for this exereise.
Measure the eross-seetional area of lhe model aquifer. Only measure lhe par! ofthe aquifer that
is eondueting waler. The unsaturated portian and ¡he clay are not part ofthe aquifer so they
shouldn 'l be included. The height varies sorne so you will ha ve to decide how to estimate a
representative height.
~e~
Salve Darey equation for K. K= -Q/iA
For example,
----;;>
Q = 100 mi in 10 minutes = 100cm'/600sec = 0.167 cm'/sec
i = 1 em/l2 cm = 0.083
A = 2.54em x 25 .4em = 64.5em'
K = -Q/iA = (0.167 cm'/sec)/(0.083 x 64.5cm' ) = 0.03 ¡ cm/sec =3.1 x 10" cm/see
Thi s is right in the range ofhydraulie conduetivity for well-sorted sand on the lable below .
.2f
.2
Hvdraulic Cond -_ ..... .,_ .. ...... _...... ..., ".-._. . .....
4> A (k..~c'¡;.id5.,f ¿~~ ~
~~17~ /.íZO¿~J
-~~ ~ ~
Measurements can be done in any unit and then converted using the conversion table shown
below. Be sure to keep trae k of units and be consistent in using units through a particular
exercise. Hydrogeologists typically work in units of centimeters per seeond (cm/s) or feet per
-- \'-~-
Additionaloplions
Raise one end ofthe model an ineh or so, determine the new values and redo aboye calculations.
Or restrict the Oow out ofthe discharge tube by raising it up, thereby reducing both the Q and the
gradient.
How do the values for K compare? (The values shauld be very clase ifno/ iden/ieal.)
Why are they similar? (K is a proper/y al/he aquifer ma/erial and no/ dependen/ on gradien/.)
You can determine K for the confined system by using a siphon to pull water out ofthe artesian
f46;k~~
pumping well to determine Q, measuring lhe gradient and cross sectional area of the artesian
---;;>
Ú?je;;';';"I H 3¡'-'-"-
Groundwater velocity
The velocity of graundwater Oow in an aquifer system can be determined for the model system.
Velocity is a function ofthe conductivity (K), the gradient (i) and the effective porosity or
interconneetedness ofthe pore spaces (ne). Q/A would be the velocity ofthe water ifthe aquifer
was all water or moving through a pipe. Because the sand takes up space, only the area of the
pore spaee, ne, is available for the water to move through. The velocity ofthe water must
increase by l/ne times so that Q amo un! of water is still moving through the sand. The following
v ~ Q/An, ~ -Ki/ n,
Where:
V ~ velocity (Iengthltime)
i = gradient (unitless)
n, = effective porosity (use .30 for medium sand, .33 for coarse sand, .27 for fine sand) (unitless)
3 :Zf
For example,
From exercise 1
K~ 85ft/day
V= Kvne =________________
Measure the distance fram an injection well to a discharge point (or sorne other intermediate
InjecI dye and monitor time it takes for dye lO reaeh point. __________________
Injeet and measure the time it takes for the dye lO reaeh an actively pumping wel1. Measure Ihe
new gradient. Is it faster or slower than the nonpumping situation? (7his will be fasler due lO Ihe
increased gradienl. K and n, remain Ihe same.)
'1 21
What implieations does this have for water management? (PlImping will accelerale Ihe
movemenloigroundWGler and any conlaminants loward lhe pumping well. Pumping may diverl
grollndwaler Ihal would normally discharge lo a river, lake or wet/and, Ihereby reducing Ihejlow
in Ihal surface water body.)
Other oplions
Find or make a water table map (eontours of head) of your area or use the one from the aetivity
Go with lhe Flow (from lhe Groundwater Study Guide) and measure lhe gradient (dh/dl). Have
the students draw a flow path (perpendicular to head eontours) and caleulale how long it would
take for eontamination from a landfill or septie tank (piek a souree) lO reaeh a well or the river.
Choose different geologie materials and see what differenee il makes in the travel time. For eaeh
different geologic material, you'lI need to seleet an appropriate hydraulie eonduetivity and
effeetive porosity from lhe tables above. For whieh geologie material is the travel time longest?
For whieh is il shortesl?
A well that is in a high transmissivity aquifer will have little drawdown during pumping so it will
have a high speeifie eapaeity . A well in a low transmissivity aquifer will have a large amount of
drawdo\Vn during pumping so it will have a low speeifie eapaeity. A general empirieal rule
relating specific capacity to transmissivity is:
T~Sc x 2000 where T is in gpdlft and Se is in gpmlft.
The factor of 2000 eontains the conversion from gpmlft lO gpdlft and would need to be adjusted
for other units.
3" 27
This is the lower portion of a well eonstruetion report, whieh has the pump test information at
the bott~~!fi,,!_,,",,
8. lrShortliN):: Swimmir.gPool l. Ck.oo"""'~J::"'1'
:' oa.llrt.tl ~ ""u. :'.-.av.
lltuILl 1ft, ,' p, --1 fM"rI" JoI..t CJl~ QlIu. -_--_-- __ n.- _C_ CLAY 10 ...
I~
8..0 , ,'llJÍlltt
----j---t---t ~~! l;::;~;:;.~=:~::;'-'-'--'-
SANOSTONE 10 :ll
EIIDV>;'flOED .:!(I'-N,t.lLIIU:U:::P .f
,..fu. I 42
For example, using information from this well eonstruction report (WCR)
.--
hKS" ;H-UP¡'¡llcIlJ /
Es¡¡~C/~ /R (
Sc = 10 gpm/4 ft= 2.5gpm/ft
T = Se x 2000 = 2.5 x 2000 = 5000 gpdlft (use the eonversion in the table to go from gallons to
ftJ)
71-u'o3
What is the transmissivity ofthe model?
In other words ifthe model were to run all day and the discharge was the entire right side ofthe
model, 509 gallons would discharge from the modeI in one day (24 hours).
Find several welllogs from your area and estimate the T from the pump test data. Or look at the
Choose publie systems and type in eity na me, seleet water supply system by name then scroll
down and ehoose entry points . Cliek on a well number and then ehose well eonstruetion data, the
speeifie eapaeity for the well is near the bottom right of the form before the tables on casing and
geology.
~ 30
Why do the Ts differ for wells in the area? (The reason Iransmissivily differs is relaled primarily
lo Ihe hydraulic conduclivily oJlhe aquijer in Ihe immediale area around Ihe \Vell. In sand and
gravel \Vells one well mighl hil a larger more continuous seam oJhigh porosily malerial and Ihus
have o higher K and resullant higher T. In bedrock wells a majar injluence an K ond Ihus T
would be ji-ocluring or bedding planes in Ihe rack Ihol ollowJor increosedjlaw. Tronsmissivily
will olso be ajJecled by Ihe deplh Ihe \Vell penelrales inlo Ihe oquijer ond lo sorne degree Ihe well
ejJiciency.)
TCUESS Program
A computer can also ... used to estimate transmissivity. Tguess is a program written by Ken
Bradbury and E. R. Rothschild which uses basic information from a well construction report or
other pump test data to estimate the transmissivity ofthe aquifer being tested. The program is
available from the Wisconsin Geological and Natural History Survey and from the ¡ntemalional
Ground Water Modeling Center in Golden, CO. The program code is available in the March
April 1985 issue ofGroundwater.
References
The following paper is downloadable from the web and has more information on the concepts of
groundwater flow.
Basic Ground-Water Hydrology, RC Heath, USGS Water Supply Paper 2220
http://water.usgs.gov/pubs/wsp/wsp22201
--
f 31