ESS 454 Hydrogeology: Instructor: Michael Brown Brown@ess - Washington.edu
ESS 454 Hydrogeology: Instructor: Michael Brown Brown@ess - Washington.edu
Hydrogeology
Module 4
Flow to Wells
• Preliminaries, Radial Flow and Well Function
• Non-dimensional Variables, Theis “Type” curve, and
Cooper-Jacob Analysis
• Aquifer boundaries, Recharge, Thiem equation
• Other “Type” curves
• Well Testing
• Last Comments
initial/final conditions
Goal here:
Need an entire course 1. Understand the basic principles
devoted to “Wells and 2. Apply a small number of well testing methods
Well Testing”
Module Four Outline
• Flow to Wells
– Qualitative behavior
– Radial coordinates
– Theis non-equilibrium solution
– Aquifer boundaries and recharge
– Steady-state flow (Thiem Equation)
• “Type” curves and Dimensionless variables
• Well testing
– Pump testing
– Slug testing
Concepts and Vocabulary
• Radial flow, Steady-state flow, transient flow, non-equilibrium
• Cone of Depression
• Diffusion/Darcy Eqns. in radial coordinates
– Theis equation, well function
– Theim equation
• Dimensionless variables
• Forward vs Inverse Problem
• Theis Matching curves
• Jacob-Cooper method
• Specific Capacity
• Slug tests
• Log h vs t
– Hvorslev falling head method
• H/H0 vs log t
– Cooper-Bredehoeft-Papadopulos method
• Interference, hydrologic boundaries
• Borehole storage
• Skin effects
• Dimensionality
• Ambient flow, flow logging, packer testing
Module Learning Goals
Assumptions
1. Aquifer bounded on bottom, horizontal and infinite, isotropic and homogeneous
2. Initially horizontal potentiometric surface, all change due to pumping
3. Fully penetrating and screened wells of infinitesimal radius
4. 100% efficient – drawdown in well bore is equal to drawdown in aquifer
5. Radial horizontal Darcy flow with constant viscosity and density
Equations in axial coordinates
Cartesian Coordinates: x, y, z
Axial Coordinates: r, q, z r q
Flow
For a cylinder of radius through
r and height b : surface of
area 2prb
Equations in axial coordinates
Darcy’s Law:
dh
Diffusion Equation: Q = KA
Area of cylinder dr
dh T 2 Leakage:
K 'Dh dh
= Ñ h+ e e= Water infiltrating through
= K 2 p rb
dt S b' confining layer with
dr
properties K’ and b’ and no
storage.
dh
= 2p Tr
Need to write in axial coordinates
dr
with no q or z dependences
ho r Gradient
needed to
induce flow
Wanted: ho-h
Drawdown as function of Drawdown must increase
distance and time to maintain gradient
h(r,t)
Radial flow
Confined Aquifer
Theis Equation
His solution (in 1935) to Diffusion equation for radial flow to well
subject to appropriate boundary conditions and initial condition:
Story: Charles Theis went to his mathematician friend C. I. Lubin who gave him
h= ho for all r at t=0 the solution to this problem but then refused to be a co-author on the paper
because Lubin thought his contribution was trivial. Similar problems in heat flow
for all time at r=infinity had been solved in the 19th Century by Fourier and were given by Carlslaw in
1921
dh
Q = -2p Tr
dr r®rw
Important step: use a non- r2S For u=1, this was the
dimensional variable that u= definition of characteristic
includes both r and t 4Tt time and length
Solutions to the diffusion equation depend only on the ratio of r2 to t!
Q ¥
e-a Q
ho - h=
4p T
ò da = W(u) W(u) is the
“Well Function”
u a 4p T
No analytic solution
Theis Equation
Need values of W for different values of the dimensionless variable u
2. Calculate “numerically”
o Matlab® command is W=quad(@(x)exp(-x)/x, u,10);
How much drawdown at well screen (r=0.5’) after 24 hours? 6.2x10-8 16.0
How much drawdown 100’ away after 24 hours? 2.5x10-3 5.4
Cone of Depression
Continues to go
down
After 1000
After 1Days
30 Days
After ofPumping
Dayof
of Pumping
Pumping
Q
Dh=2.4 feet
T=
4p W(ho - h)
Match point
at u=1 and 4Ttu
W=1 S=
r2
time=4.1 minutes
Semilog Plot of “Non-
equilibrium” Theis equation
After initial time, drawdown increases with log(time)
Ideas:
Linear drawdown
1. At early time water is 2T
delivered to well from Initial non- T
linear curve
“elastic storage” Double T -> slope
then linear with decreases to half
head does not go down
log(time)
much
Larger intercept for larger
Log time
storage Intercept time increases with S
2. After elastic storage is
Delivery from
depleted water has to elastic storage De
flow to well liv
er
yf
Double S and
ro
Head decreases to m
flo
w
intercept changes but
maintain an adequate slope stays the same
hydraulic gradient
Rate of decrease is
inversely proportional to T
Cooper-Jacob Method
Theis Well function in series expansion
2 3
Theseutermsu u4negligible
become
W(u) @ -0.5772 - ln(u) + u- + - +
as time
4 goes18
on 96
Theis equation -Q é æ r 2 S öù -Q é æ r 2 S öù Q
(ho - h) = ê0.5772 + ln ç ÷ú = êln(1.78) + ln ç ÷ú + ln ( t)
for large t 4p T ë è 4Tt øû 4p T ë è 4T øû 4 p T
Q æ 4T ö Qpe
Head decreases linearly with log(time) = constant
ln ç ÷ + lo ln ( t)
4p T è 1.78r 2 S ø 4sp T
– slope is inversely proportional to T
2.3Q Conversion to
– constant is proportional to S = constant+ log(t) base 10 log
4p T
Cooper-Jacob Method
Works for “late-time” drawdown data
Solve inverse problem: Given drawdown vs time data for a
well pumped at rate Q, what are
Using equations from the aquifer properties T and S?
previous slide
intercept
to
2.3Q Calculate T Fit line through linear range of data
T= from Q and Dh
Need to clearly see “linear” behavior
Line defined by slope and intercept
4pDh Not
acc
ept
abl
e
Slope =Dh/1
2.25Tto Dh for 1 log unit
S= Need T, to and
r2 r to calculate S
1 log unit
Summary
• Recognize causes for departure of well drawdown data from the Theis
“non-equilibrium” formula
• Be able to explain why a pressure head is necessary to recover water from a
confined aquifer
• Be able to explain how recharge is enhanced by pumping
• Be able to qualitatively show how drawdown vs time deviates from Theis
curves in the case of leakage, recharge and barrier boundaries
• Be able to use diffusion time scaling to estimate the distance to an aquifer
boundary
• Understand how to use the Thiem equation to determine T for a confined
aquifer or K for an unconfined aquifer
• Understand what Specific Capacity is and how to determine it.
When Theis Assumptions Fail
D
h Aquifer above Aquitard
K 'Dh
e=
b' Increased flow through aquitard
Confined Aquifer
Flow to Well in Confined Aquifer with
Recharge Boundary
surface
ho: Initial potentiometric surface
Lake
o well
e ad t
fixed h
t from
d ie n
G ra
Confined Aquifer
Flow to Well –Transition to
Steady State Behavior
Q dr
dh= Rearrange h1 h2
2p T r r1
r2 Confined Aquifer
h2 r2
Q dr
ò dh =
2p T
ò r
Integrate both sides
h1 r1
Determine T from drawdown
at two distances
h2 - h1 =
Q r Q r
ln 2 Result T= ln 2 In Steady-state – no
2p T r1 2p (h2 - h1 ) r1 dependence on S
Steady-State Flow
Thiem Equation – Unconfined Aquifer
When hydraulic head does not change with time
db
Q = 2p Kbr Darcy’s Law in radial coordinates surface
dr
Q dr
bdb= Rearrange b1 b2
2p K r r1
r2
b2 r2
Q dr
ò bdb =
2p K
ò r
Integrate both sides
b1 r1
Determine K from drawdown
at two distances
Q r Q r2
b22 - b12 = ln 2 Result K= ln In Steady-state – no
p K r1 p (b22 - b12) r1 dependence on S
Specific Capacity (driller’s term)
1. Pump well for at least several hours – likely not in steady-state
Q=21*.134*60*24
= 4.1x103 ft3/day
Static head is 15’ below surface
Specific capacity of:
=4.1x103/8=500 ft2/day
Pumped at 21 gallons/minute K is about 100 ft/day
for 2 hours
(typical “good”
Drawdown of 8’
sand/gravel value)
The End: Breakdown of Theis assumptions and steady-state behavior
Q
ho - h= W(u, r B) Larger r/B -> smaller steady-state drawdown
4p T
Drawdown reaches
“steady-state” when
Large K’
Tb' recharge balances flow
B= makes r/B
K' large “Type Curves” to determine T, S, and r/B
Other Type-curves –
Given without Derivations
2 Difference between r 2 Sy
rK uB =
Later flow from
G= 2 v vertical and horizontal
gravity draining
b Kh conductivity is 4Tt
important
Flow in Unconfined Aquifer
1. Elastic Storage
Flow from elastic storage
Dimensionless number
Goes from 1 to 0
Plot: H(t)/Ho vs log(Tt/rc2) Call it z
rc surface
1
Initial head
Head returns
slug
to initial state
Increased head
0 causes radial
b
0.01 0.1 1.0 10.0 flow into aquifer
z=Tt/rc2
rs
Cooper-Bredehoeft-PapadopolosTest
rs= 1.0’
rc= 0.5’ .8
Tt
z= .6
rc2
zrc2
H/Ho
1*.25 ft2
T= =
1 1 day
t 21* * .4
60 24
rs2
m = 10-2
m= 2S .2
rc
Log scale
K only determined
1
.8 2 æ Le ö
.7 r ln ç ÷
casing
è Rø
.6 K=
.5 H/Ho=.37 2Let37
H/Ho
.4
.3
t37
Le
Screen
.
1 2 4 6 8 10
minutes
Linear scale Partially Penetrating OK
R high K material
The End: Well Testing
– Borehole storage
– Skin effect
– Aquifer Storativity
– Aquifer Transmissivity
– Recharge/barrier boundaries
Well interference
Greater drawdown
Hydraulic head is measure of energy Smaller hydraulic gradient
Energy is a scalar and is additive Reduced flow to wells
Just add drawdown for each well to get total drawdown Flow divide between wells
Boundary and Dimension Effects
2-D
1-D
3-D
Reservoir geometry
Network/Flow geometry
Discussion of ways to deal with these “real-world” situations is beyond the scope of this class
Last Comments on well testing