C.P.hazel Groundwater Hydraulics Final 20111031
C.P.hazel Groundwater Hydraulics Final 20111031
C.P.hazel Groundwater Hydraulics Final 20111031
by
C.P. HAZEL
Second Edition
June 2009
Lectures presented by C.P.Hazel of the Irrigation and Water Supply Commission, Queensland to the Australian Water
Resources Council's Groundwater School, Adelaide in 1973. The lectures were presented initially in Imperial units at
the 1973 Groundwater School and then converted to metric units for use at the 1975 and subsequent Groundwater
Schools. They were retyped and updated in 2009.
TABLE OF CONTENTS
PREFACE .............................................................................................................................. 1
SECTION 1: PROPERTIES OF WATER AND WATER BEARING MATERIALS ..................... 3
1.1 INTRODUCTION ................................................................................................................... 3
1.2 FLUID MECHANICS ............................................................................................................... 3
1.2.1 Hydrostatics ..................................................................................................................... 3
1.2.2 Hydrodynamics ................................................................................................................ 5
1.3 SOIL MECHANICS ................................................................................................................10
1.3.1 Grain-Void Relationship ....................................................................................................10
1.3.2 Fluid Flow Properties .......................................................................................................11
1.3.3 Soil Pressures..................................................................................................................12
1.3.4 Properties of Rock Types .................................................................................................13
SECTION 2: OCCURRENCE OF GROUNDWATER ........................................................... 17
2.1 INTRODUCTION ..................................................................................................................17
2.2 ORIGIN OF GROUNDWATER....................................................................................................17
2.3 HYDROLOGIC CYCLE ............................................................................................................17
2.4 FACTORS AFFECTING THE ABSORPTION OF WATER ......................................................................19
2.5 VERTICAL DISTRIBUTION OF SUB-SURFACE WATER......................................................................19
SECTION 3: AQUIFERS ................................................................................................. 21
3.1 DEFINITIONS .....................................................................................................................21
3.2 AQUIFER FUNCTIONS ...........................................................................................................23
3.3 TYPES OF AQUIFER FORMATIONS ............................................................................................23
3.4 HYDRAULIC PROPERTIES .......................................................................................................23
SECTION 4: GROUNDWATER FLOW ............................................................................. 30
4.1 INTRODUCTION ..................................................................................................................30
4.2 DARCY’S LAW.....................................................................................................................30
4.3 HYDRAULIC CONDUCTIVITY (K) ..............................................................................................32
4.4 RELATION OF K TO PARTICLE VELOCITY ....................................................................................33
4.5 RELATION OF K TO INTRINSIC PERMEABILITY .............................................................................33
4.6 REYNOLD’S NUMBER ............................................................................................................35
4.7 RANGE OF VALIDITY OF DARCY’S LAW ......................................................................................35
4.8 GROUNDWATER FLOW RATE ..................................................................................................35
4.9 FLOW ANALOGIES ...............................................................................................................36
4.10 TYPES OF GROUNDWATER FLOW .............................................................................................36
4.11 STATES OF GROUNDWATER FLOW ...........................................................................................36
SECTION 5: BORE DISCHARGE TESTS .......................................................................... 38
5.1 INTRODUCTION ..................................................................................................................38
5.2 BACKGROUND ....................................................................................................................38
5.3 DEFINITIONS .....................................................................................................................38
5.4 FLOWING AND NON-FLOWING BORES .......................................................................................39
5.5 PLANNING A PUMPING TEST ...................................................................................................40
5.5.1 Test Design.....................................................................................................................40
5.5.2 Identify Site Constraints ...................................................................................................41
5.5.3 Purpose of the Test .........................................................................................................42
5.5.4 Specify Test Conditions ....................................................................................................42
5.5.5 Pumping rate and bore diameter ......................................................................................42
5.5.6 Bore Depth and Bore Screen ............................................................................................42
5.5.7 Observation Bores and Piezometers ..................................................................................43
5.6 MEASUREMENTS .................................................................................................................44
5.6.1 Time ...............................................................................................................................45
5.6.2 Water Levels/Heads .........................................................................................................45
5.6.3 Discharge Rate ................................................................................................................45
5.6.4 Temperature ...................................................................................................................45
5.6.5 Water Quality ..................................................................................................................46
5.7 SETUP AND INSTRUMENTATION ...............................................................................................46
5.8 DATA RECORDING AND PRESENTATION .....................................................................................51
5.8.1 Possible Corrections to Drawdown Data ............................................................................51
5.9 TESTING NON-FLOWING BORES ..............................................................................................51
5.9.1 Antecedent Conditions .....................................................................................................51
GROUNDWATER HYDRAULICS TOC i
5.9.2 Constant Discharge Test ..................................................................................................51
5.9.3 Recovery Test .................................................................................................................52
5.9.4 Constant Drawdown Test .................................................................................................52
5.9.5 Step Drawdown Test .......................................................................................................52
5.9.6 Step Drawdown Test (Extended First Step) .......................................................................53
5.9.7 Variable Discharge/Variable Drawdown Test .....................................................................54
5.9.8 Multiple Aquifer Testing ...................................................................................................54
5.9.9 Applicability of Testing Procedures ...................................................................................54
5.9.10 Pump Stoppages .............................................................................................................54
5.9.11 Slug Tests .......................................................................................................................55
5.10 TESTING FLOWING BORES .....................................................................................................56
5.10.1 Antecedent Conditions .....................................................................................................56
5.10.2 Risk in Closing Low Pressure Bores ...................................................................................56
5.10.3 Flow Recession Test (Constant Drawdown) .......................................................................56
5.10.4 Static Test (Recovery) .....................................................................................................57
5.10.5 Dynamic (Step Drawdown) Tests......................................................................................58
5.10.6 Opening Dynamic Test .....................................................................................................58
5.10.7 Closing Dynamic Test ......................................................................................................59
5.10.8 Order of Tests .................................................................................................................60
5.11 DISINFECTION....................................................................................................................60
SECTION 6: EVALUATION OF AQUIFER PROPERTIES USING OBSERVATION BORES .. 61
6.1 INTRODUCTION ..................................................................................................................61
6.2 SELECTING THE TYPE OF ANALYSIS ..........................................................................................61
6.3 CONFINED AQUIFER TEST ANALYSIS ........................................................................................62
6.3.1 Constant Discharge Tests .................................................................................................62
6.3.2 Variable Discharge Tests ..................................................................................................87
6.3.3 Other Methods ................................................................................................................87
6.4 SEMI-CONFINED AQUIFER TEST ANALYSIS .................................................................................87
6.4.1 General...........................................................................................................................87
6.4.2 Constant Discharge .........................................................................................................88
6.5 UNCONFINED AQUIFERS WITHOUT DELAYED YIELD......................................................................93
6.5.1 Constant Discharge .........................................................................................................93
6.5.2 Jacob’s Corrections for Drawdowns in Thin Unconfined Aquifers .........................................97
6.6 UNCONFINED AQUIFERS WITH DELAYED YIELD AND SEMI-UNCONFINED AQUIFERS ............................. 101
6.6.1 General......................................................................................................................... 101
6.6.2 Boulton’s Method........................................................................................................... 102
6.7 SOFTWARE ...................................................................................................................... 110
6.8 IDENTIFYING AQUIFER TYPE FROM TEST DATA ......................................................................... 110
SECTION 7: BORE PERFORMANCE TESTS................................................................... 112
7.1 INTRODUCTION ................................................................................................................ 112
7.2 EQUATION TO DRAWDOWN .................................................................................................. 112
7.3 EVALUATION OF AQUIFER PARAMETERS................................................................................... 113
7.3.1 Constant Discharge Test Analysis ................................................................................... 114
7.3.2 Variable Discharge Test Analysis .................................................................................... 114
7.4 EVALUATION OF NON-LINEAR HEAD LOSSES ............................................................................ 126
7.4.1 Drawdown Method ........................................................................................................ 126
7.4.2 Pressure Differential Method .......................................................................................... 127
7.4.3 Range of Intercepts ....................................................................................................... 128
7.4.4 Step Drawdown Test Analysis ........................................................................................ 129
7.4.5 Graphical Analysis ......................................................................................................... 130
7.4.6 Eden-Hazel Analysis....................................................................................................... 135
7.5 INTERMITTENT PUMPING TEST ANALYSIS ................................................................................ 143
7.6 EVALUATION OF LONG TERM PUMPING RATE ............................................................................ 148
7.7 SPECIFIC CAPACITY ........................................................................................................... 151
7.8 EVALUATION OF BORE EFFICIENCY......................................................................................... 152
7.8.1 When the Equation to Drawdown is Known ..................................................................... 153
7.8.2 When the Equation to Drawdown is Not Known ............................................................... 155
SECTION 8: EVALUATION OF AQUIFER PROPERTIES WITHOUT PUMPING TESTS.... 156
8.1 INTRODUCTION ................................................................................................................ 156
8.2 AREAL METHODS .............................................................................................................. 156
GROUNDWATER HYDRAULICS TOC ii
8.2.1 Numerical Analysis ........................................................................................................ 156
8.2.2 Flow-Net Analysis .......................................................................................................... 158
8.3 ESTIMATING TRANSMISSIVITY .............................................................................................. 161
8.3.1 General......................................................................................................................... 161
8.3.2 Specific Capacity of Bores .............................................................................................. 161
8.3.3 Rough Method .............................................................................................................. 163
8.3.4 Logs of Bores ................................................................................................................ 164
8.3.5 Laboratory Analysis ....................................................................................................... 164
8.4 ESTIMATING STORAGE COEFFICIENT AND SPECIFIC YIELD ........................................................... 165
8.4.1 Confined Aquifers .......................................................................................................... 165
8.4.2 Unconfined Aquifers ...................................................................................................... 165
8.4.3 Water Balance ............................................................................................................... 166
8.4.4 Barometric Efficiency ..................................................................................................... 166
8.4.5 Tidal Efficiency .............................................................................................................. 167
SECTION 9: CORRECTIONS AND EFFECTS TO BE ALLOWED FOR WHEN ANALYSING 168
9.1 GENERAL ........................................................................................................................ 168
9.2 DELAYED YIELD FROM STORAGE ........................................................................................... 168
9.3 INCREASED DRAWDOWN CAUSED BY DEWATERING .................................................................... 168
9.4 ANOMALIES IN DRAWDOWN READINGS ................................................................................... 169
9.5 PARTIAL PENETRATION ....................................................................................................... 169
9.6 ANTECEDENT CONDITIONS .................................................................................................. 170
9.7 POSSIBLE DEVELOPMENT DURING PUMPING ............................................................................. 171
9.8 PROXIMITY OF BOUNDARIES ................................................................................................ 171
9.8.1 Method of Images ......................................................................................................... 171
9.9 WATER TEMPERATURE VARIATIONS IN HOT BORES .................................................................... 178
9.10 VARIATIONS IN ATMOSPHERIC PRESSURE ................................................................................ 178
9.11 TIDAL EFFECTS................................................................................................................. 178
9.12 OTHER FACTORS TO BE CONSIDERED ..................................................................................... 179
SECTION 10: APPLICATION OF AQUIFER PROPERTIES ............................................... 180
10.1 INTRODUCTION ................................................................................................................ 180
10.2 VOLUME IN STORAGE ......................................................................................................... 180
10.3 VOLUME REMOVED FROM STORAGE ....................................................................................... 180
10.4 GROUNDWATER FLOW ........................................................................................................ 181
10.5 LEAKAGE ...................................................................................................................... 181
10.6 DRAWDOWN INTERFERENCE EFFECTS ..................................................................................... 182
10.6.1 Drawdown Within the Area of Influence .......................................................................... 182
10.6.2 Comparative Spread of Area of Influence ........................................................................ 183
10.6.3 Determination of Radius of Influence .............................................................................. 183
10.7 DRAINAGE PROBLEMS......................................................................................................... 184
10.7.1 Mine Dewatering ........................................................................................................... 184
SECTION 11: GROUNDWATER MANAGEMENT .............................................................. 190
11.1 GROUNDWATER YIELD ANALYSIS ........................................................................................... 190
11.1.1 Bore Yields.................................................................................................................... 190
11.1.2 Aquifer Yields ................................................................................................................ 193
11.2 CONTROL OF GROUNDWATER USE ......................................................................................... 196
11.3 CONJUNCTIVE USE OF GROUNDWATER AND SURFACE WATER ....................................................... 196
11.4 GROUNDWATER RECHARGE .................................................................................................. 197
11.4.1 What is Recharge? ........................................................................................................ 197
11.4.2 Definitions .................................................................................................................... 197
11.4.3 Necessity for Recharge .................................................................................................. 197
11.4.4 Natural Recharge .......................................................................................................... 198
11.4.5 Artificial or Managed Recharge ....................................................................................... 199
11.5 SEA WATER INTRUSION IN COASTAL AQUIFERS ........................................................................ 203
11.5.1 General......................................................................................................................... 203
11.5.2 Ghyben-Herzberg Concept ............................................................................................. 203
11.5.3 The Dynamic Concept .................................................................................................... 206
11.5.4 Location of the Interface ................................................................................................ 208
11.5.5 Structure of the Interface .............................................................................................. 212
11.5.6 Control of Intrusion ....................................................................................................... 212
REFERENCES .................................................................................................................... 213
GROUNDWATER HYDRAULICS TOC iii
GLOSSARY OF SYMBOLS USED ........................................................................................ 218
METRIC MULTIPLES......................................................................................................... 221
THE GREEK ALPHABET ..................................................................................................... 221
INDEX………………………………………………………………………………………………………. 222
FIGURES
Figure 1-1 Steady flow ............................................................................................................. 5
Figure 1-2 Continuity principle .................................................................................................. 5
Figure 1-3 Dynamic viscosity .................................................................................................... 7
Figure 1-4 Capillary rise ........................................................................................................... 9
Figure 1-5 Capillary rise in tubes .............................................................................................10
Figure 1-6 Pressure distribution ...............................................................................................13
Figure 1-7 Properties of pure water .........................................................................................16
Figure 2-1 Hydrologic cycle .....................................................................................................18
Figure 2-2 Vertical distribution of sub-surface water .................................................................20
Figure 3-1 Aquifer types .........................................................................................................22
Figure 3-2 Homogeneous anisotropic formation ........................................................................24
Figure 4-1 Laminar flow in a porous medium ............................................................................31
Figure 5-1 Flowing and non-flowing bores ................................................................................40
Figure 5-2 Cross section of a confined aquifer (after Kruseman and de Ridder, 1990) .................41
Figure 5-3 Cross section of an unconfined aquifer (after Kruseman and de Ridder, 1990) ...........41
Figure 5-4 Common discharge measuring devices.....................................................................48
Figure 5-5 Some common water level measuring devices ..........................................................49
Figure 5-6 Typical bore hole pump installation..........................................................................50
Figure 6-1 Steady state flow derivation – confined aquifer ........................................................62
Figure 6-2 Steady state flow example, confined aquifer ............................................................66
Figure 6-3 Non-steady state flow derivation - confined aquifer ..................................................67
Figure 6-4 Type curves for non-steady state flow in leady aquifer .............................................75
Figure 6-5 Type curve solution, confined aquifer, non-steady state, constant Q ..........................77
Figure 6-6 Modified non-steady state flow example – confined aquifer, constant Q, constant r,
varying t ................................................................................................................81
Figure 6-7 Modified non-steady state flow example – confined aquifer, constant Q, constant t,
varying r ................................................................................................................83
Figure 6-8 Steady state flow example – semi-confined aquifer ..................................................91
Figure 6-9 Type curve solution, semi-confined aquifer, non-steady state ....................................95
Figure 6-10 Unconfined aquifer, steady state flow derivation ......................................................96
Figure 6-11 Steady state flow example – unconfined aquifer .......................................................99
Figure 6-12 Unconfined aquifer, variation of S with time ........................................................... 100
Figure 6-13 Delayed yield type curves ..................................................................................... 104
Figure 6-14 Boulton’s delay index curve .................................................................................. 108
Figure 6-15 Unconfined aquifer with delayed yield, non-steady state flow example .................... 109
Figure 6-16 Typical response curves for different aquifer types ................................................. 111
Figure 7-1 Constant drawdown example - straight line solution ............................................... 121
Figure 7-2 Constant drawdown test example using Eden-Hazel method ................................... 125
Figure 7-3 Step drawdown test – graphical analysis example .................................................. 133
Figure 7-4 Step drawdown test – graphical analysis, determination of “a” and “C” .................... 134
Figure 7-5 Step drawdown test – Eden-Hazel analysis............................................................. 138
Figure 7-6 Step drawdown test – Eden-Hazel analysis, determination of “a” and “C” ................ 139
Figure 7-7 Drawdown versus discharge curves for various times of discharge .......................... 142
Figure 7-8 Intermittent pumping – “F” versus “n” curves ........................................................ 146
Figure 7-9 Intermittent pumping – “F” versus “p” curves ........................................................ 147
Figure 7-10 Determination of long term pumping rate .............................................................. 150
Figure 8-1 Numerical analysis array ....................................................................................... 156
Figure 8-2 Numerical analysis example .................................................................................. 158
Figure 8-3 Typical flow net .................................................................................................... 159
Figure 8-4 Elemental square.................................................................................................. 160
Figure 8-5 Typical specific capacity - time - discharge curves .................................................. 163
Figure 9-1 Idealised section views of a discharging well in a semi-infinite aquifer bounded by a
perennial stream, and of the equivalent hydraulic system in an infinite aquifer ......... 172
GROUNDWATER HYDRAULICS TOC iv
Figure 9-2 Generalised flow net showing stream lines and potential lines in the vicinity of a
discharging well dependent upon induced infiltration from a nearby stream ............. 173
Figure 9-3 Idealised section views of a discharging bore in a semi-infinite aquifer bounded by an
impermeable formation, and of the equivalent hydraulic system in an infinite aquifer174
Figure 9-4 Generalised flow net showing stream lines and potential lines in the vicinity of a
discharging well near an impermeable boundary .................................................... 175
Figure 9-5 Family of type curves for the solution of the modified Theis formula ........................ 177
Figure 10-1 Recharging pumped water to maintain water levels in sensitive areas...................... 187
Figure 10-2 Groundwater flow into a strip pit ........................................................................... 188
Figure 10-3 Multiple aquifer flow into a strip pit ....................................................................... 188
Figure 11-1 Constant discharge test – Callide Valley ................................................................. 192
Figure 11-2 Stable saltwater interface...................................................................................... 204
Figure 11-3 The dynamic saltwater interface ............................................................................ 207
Figure 11-4 Saltwater wedge in a confined aquifer ................................................................... 209
Figure 11-5 Saltwater wedge in an unconfined aquifer .............................................................. 211
TABLES
Colin P. Hazel
June 2009
where:
η = dynamic viscosity.
τ = intensity of shear.
dv/dy = velocity gradient in the transverse direction.
v
F
Fluid y
v( 4 R)
NR
where R is the hydraulic radius, and is equal to the ratio of the cross-sectional area to the wetted
perimeter.
Surface Tension (σ)
Surface Tension is a property associated with the free surface of any liquid or the interface between
any two non-miscible liquids.
It is well known that many insects are able to walk on the surface of liquids in apparent contradiction
to Archimedes' Principle. This property tends to suggest that there is a kind of membrane or skin that
envelopes all liquids. In fact this very nearly describes what actually the situation is.
A molecule in the interior of a fluid is acted upon by attractive forces in all directions and the vector
sum of these is zero. However, at the surface, a molecule is acted upon by a net inward cohesive
force perpendicular to the surface. Hence work is required to bring molecules to the surface.
The surface tension of a liquid is the work that must be done to bring enough molecules from inside
the liquid to form one new unit area on the surface.
The intermolecular forces also come into play when a liquid is in contact with a solid object - in this
case there is an attraction of the molecules of the solid for those of the liquid.
Hence, for example, if pure water is placed in a clean glass container it will be observed that the
water surface in contact with the glass turns up and lies flat on the glass. Thus the attraction of glass
for water is greater than that of water for itself.
On the other hand mercury placed in a glass container will be observed at the surface contact, to be
drawn away from the glass indicating that the attraction of glass for mercury is less than that of
mercury for itself.
Capillarity or the rise or fall of a liquid in a capillary tube is caused by surface tension and depends on
the relative magnitudes of the cohesion of the liquid and the adhesion of the liquid to the walls of the
containing vessel. Liquids rise in tubes they wet (adhesion > cohesion) and fall in tubes they do not
wet (cohesion > adhesion).
GROUNDWATER HYDRAULICS Page 8 of 225
α
2r p = atmospheric
hc p < atmospheric
density ρ, surface
tension σ
where:
ip = the pressure gradient.
η = the dynamic viscosity of the fluid.
k = the intrinsic permeability (square micrometre).
This may be written:
k gi
v
.....1.25
where:
i = hydraulic gradient.
This equation is of the same form as the Darcy equation:
v = Ki
Hence:
gk kg .....1.26
K
The intrinsic permeability (k) is related solely to the properties of the porous medium. The coefficient
of permeability or hydraulic conductivity (K) is related not only to the properties of the porous medium
but also to the properties of the fluid.
From extensive laboratory testing, Hazen found that the coefficient of permeability of sands in a loose
state depended on two quantities he called the effective grain size and the uniformity coefficient.
The effective grain size D10 of a sample is a grain size diameter such that 10 percent of the particles
are finer and 90 percent coarser.
Dn is defined as that diameter such that n% of the particles in a sample are finer.
The uniformity coefficient U is defined as:
D60
U
D10 .....1.27
1.3.3 Soil Pressures
Pore Water Pressure
If the pores or interstices of a porous medium are filled with water then this water is subject to the
same principles as outlined in section 1.1.1, hydrostatics.
Capillary
Capillary Water Zone
Zone of
Saturation Groundwater
Bedrock
K1 b1
b
K2 b2
y K3 b3
b
K v eff
n bm .....3.3
m 1
K vm
where:
Kveff = the effective hydraulic conductivity for vertical flow.
Kvm = the vertical hydraulic conductivity of the mth layer.
bm = the thickness of the mth layer.
b = the total thickness of the aquifer.
Equations 3.2 and 3.3 are analogous to the flow of electricity through a series of resistors. In equation
3.2 the resistors are in parallel and in equation 3.3 the resistors are in series.
The dimensions of hydraulic conductivity are Volume/time/area (Length/Time). The units to be used
are m/day or m3/day/m.
where:
b = thickness of the aquifer.
ρ = density of the fluid (for water = 1,000 kg/m3).
g = acceleration due to gravity.
θ = porosity.
β = compressibility of water, 4.8 x 10-7 kPa i.e.4.8 x 10-10 m2/N.
Ew = 1/β bulk modulus of water 2.08 x 106 kPa.
α = compressibility of soil matrix. Sand and gravel α ~ 10-8 m2/N
Es = 1/α = bulk modulus of soil matrix.
C = a dimensionless ratio, which may be considered unity in an uncemented granular
material. In a solid aquifer such as limestone having tubular solution cavities, C is apparently
equal to the porosity. The value for a sandstone ranges between these limits, depending
upon the degree of cementation.
Bear (1972, p.207), suggests that the above expression is in error and the correct expression for
porous media should be:
S b g( (1 ))
(1 )
b g( )
Ew Es .....3.6
Bear's approach considers a constant control volume with both the fluid and soil matrix being free to
move across the boundaries when the porous media is subjected to varying pressures. Jacob
considered that the control volume itself undergoes deformation.
' 1 dUb
b
i.e. Ub d .....3.8
if the hydrostatic pressure (p) is held constant.
where:
Ub = the volume of a fixed mass of porous medium.
Because groundwater is generally associated with a relatively constant external stress (σ) and a
variable hydrostatic pressure (p) another coefficient of compressibility, α b, is often defined with
respect to a unit change in p.
1 dU b
b
U b dp .....3.9
where σ is now held constant.
Two more kinds of compressibility in addition to αb have been proposed as:
1. rock (or solid) matrix compressibility, αs, which is the fractional change in volume of the
solid matrix (Us) with unit change in p;
1 dU s
s
i.e.
U s dp .....3.10
2. pore compressibility, αp, defined as the fractional change in pore volume (Up) with unit
change in p.
1 dU p
p
U p dp
i.e. .....3.11
It follows then, since:
Ub Us Up
dU b dU s dU p
dp dp dp
and since:
1 (1 )
Us (1 )U b U Us
i.e. b
and
1
Up Ub Ub Up
i.e.
1 dU b
b
U b dp
(1 ) dU s dU p
U s dp U p dp
(1 ) s p
.....3.12
From the definition of S (by Bear):
S b g( (1 )
where:
β is now written for αp and α is written for αs.
The major bracketed term is the compressibility of the porous medium, and bρg is in the form of a
stress. Since the compressibility is the inverse of the bulk modulus, then the Storage Coefficient is in
the form of a strain of the porous medium.
It can be thought of as the amount of deformation of water and matrix per unit increase in hydrostatic
pressure, or can be referred to as the elastic storage. The major bracketed term is the compressibility
of the porous medium, and bρg is in the form of a stress. Since the compressibility is the inverse of
the bulk modulus, then the Storage Coefficient S is in the form of a strain of the porous medium.
It can be thought of as the amount of deformation of water and matrix per unit increase in hydrostatic
pressure, or can be referred to simply as the elastic storage.
Specific Mass Storativity or Specific Storage (Ss)
The Specific Mass Storativity (Ss) or Specific Storage is defined as the volume of water released from
or taken into storage per unit volume of aquifer per unit change in head.
(1 )
Ss g(
Ew Es .....3.13
S
Ss
i.e. b
where:
S = Storage Coefficient.
b = thickness of aquifer.
Specific Mass Storativity has the dimensions length-1.
It has the units metres-1.
Specific Yield or Phreatic Storage Coefficient (S) (Unconfined Aquifers)
The specific yield of an aquifer is the volume of water which will drain under gravity from a unit
volume of aquifer. For a section of aquifer it is the ratio of the drainable water to the saturated
volume.
An unconfined aquifer has elastic storage properties as does a confined aquifer, but these are so small
in comparison with the non-elastic storage properties that the specific yield is normally referred to as
the Storage Coefficient of an unconfined aquifer.
It is used to determine the recoverable volume of water stored between the standing water level and
a specified dead storage level.
For an unconfined aquifer the specific yield ranges from about 0.1 to 0.3.
Specific yield has no dimensions.
Kb
B
Sy
or
T
B
Sy
.....3.18
where:
K = hydraulic conductivity of the aquifer.
b = aquifer thickness.
T = transmissivity of the aquifer.
1/α = the Boulton delay index (an empirical constant).
Sy = the specific yield after a long pumping time.
Large values of B indicate a fast drainage. If B = ∞ the yield is instantaneous with the lowering of the
water table, so the aquifer would be confined without delayed yield.
The dimensions of B are in length, and the units are metres.
Q dh
q K
A dl .....4.3
where:
q = discharge per unit area.
This statement, that the rate of flow of water through porous media is proportional to the head loss
and inversely proportional to the length of the flow path, is known universally as Darcy's Law. Darcy's
Law forms the basis for the present day knowledge of groundwater flow.
As water percolates through a permeable material, the individual water particles move along paths
which deviate erratically but only slightly from smooth curves known as "flow lines". If adjacent flow
lines are parallel the flow is said to be "linear" or "laminar".
The hydraulic principles involved in laminar flow are illustrated Figure 4-1.
In Figure 4-1, the points A and B represent the extremities of the flow lines. At each extremity a
stand pipe, known as a potentiometric tube, has been installed to indicate the level to which the water
rises at these points. The water level in the tube at B is designated as the potentiometric level at B
and the vertical distance from this level to point B is the pressure head at B. The vertical distance
between A and B represents the "position head" or difference in elevation heads above a set datum
level.
k w g
K
.....4.8
Hence:
v = -Ki .....4.9
where:
v = average discharge velocity or unit discharge.
K = coefficient of permeability, or hydraulic conductivity.
i = hydraulic gradient.
In Civil Engineering, the value K is commonly called the "coefficient of permeability". In Groundwater
Hydrology K is called the hydraulic conductivity. Equation 4.9 is commonly known as Darcy's Law.
Darcy's Law may also be written in the form:
Q = -KiA .....4.10
where:
Q = discharge rate through an area A.
i = the hydraulic gradient.
A = the cross-sectional area normal to the direction of flow.
K = the hydraulic conductivity of the material.
m3
K
m 2 day.m.m 1
which is the volumetric rate of flow through a given cross-sectional area under unit gradient. For the
average actual or particle velocity, we must know also the porosity of the medium. Thus,
dh
Q vA KA
dl .....4.11
K dh
v
dl .....4.12
where:
v = average water particle velocity.
θ = porosity, as a decimal fraction.
dh/dl = hydraulic gradient = i.
k K
Material
(millidarcys) (m/day)
Silt, sandy silts, clayey sands, till 1 - 102 8.64 x 10-4 - 8.64 x 10-2
Well sorted sands, glacial outwash 103 - 105 8.64 x 10-1 - 8.64 x 101
The channels through which the water particles travel in a mass of soil have a variable and irregular
cross-section. As a consequence, the real velocity of flow is extremely variable. However, the average
flow through such channels is governed by the same laws that determine the rate of flow through
straight capillary tubes having a uniform cross-section.
If the cross-section of the tube is circular, the velocity of flow increases with the square of the
diameter of the tube. Since the average diameter of the voids in soil at a given porosity increases
practically in proportion to the grain size D it is possible to express K as:
K = constant x D2 .....4.15
From his experiments with loose filter sands of high uniformity (uniformity coefficient not greater than
about 2), Allen Hazen obtained the empirical equation,
5.2 BACKGROUND
Generally pumping tests are carried out on bores for one of two reasons:
to determine the hydraulic characteristics of the aquifer. These aquifer characteristics,
such as Transmissivity and Storage Coefficient, are used to determine the ability of the
aquifer to store and transmit water and hence assess its response to stresses such as
recharge and discharge; or
to determine the long term pumping capability of the bore itself under sustained
pumping. This is important for the correct selection of pumping equipment for
commissioning of the bore or for the determination of the number and spacing of bores
to achieve a required objective.
Occasionally pumping tests are carried out for other reasons including:
to determine the existence and location of sub-surface boundaries which may affect
adversely or beneficially the long term pumping performance of a particular bore;
to check on the performance of a particular groundwater basin; or
to determine the radius of influence of a bore for dewatering or interference purposes.
On occasions it is not feasible or possible to carry out a pumping test. (The bore diameter may be
very small and preclude the insertion of pumping equipment or the hydraulic conductivity of the
aquifer is such that the pumping rate would be exceedingly small.) However, it may be essential that
the hydraulic conductivity be determined. In cases such as this slug tests are frequently used to
calculate the hydraulic conductivity in the near vicinity of the bore being tested. Such determinations
are representative of area in the immediate vicinity of the hole and not of the aquifer in general.
Whether flowing or non-flowing conditions are concerned, the testing process involves the removal of
water from the aquifer at a controlled rate whilst monitoring the head or change in head (water level
or potentiometric level) response over time.
5.3 DEFINITIONS
Before proceeding with the various types of tests a number of terms will be defined.
Figure 5-3 Cross section of an unconfined aquifer (after Kruseman and de Ridder,
1990)
5.5.2 Identify Site Constraints
Conditions at the site impose constraints on pumping test design which should be recorded and taken
into account prior to the start of the design process. A few examples, some of which are irrelevant if
the bore is already constructed, are:
limitation of placement (e.g. buildings, roads, railway lines or other bores);
limitation on pumping rate(s);
limitation on test duration;
local authority limitation on discharge of water;
location of the point of discharge for the pumped water to prevent recirculation;
use of an existing bore (this may limit number and location of observation bores,
problems with effects of bore storage, head losses); and
limitation on the placement of bores due to known presence of aquifer discontinuities,
presence of recharge and discharge zones.
<545 152 ID
409-954 203 ID
818-1910 254 ID
1640-3820 305 ID
2730-5450 356 OD
4360-9810 406 OD
6540-16400 508 OD
10,900-20,700 610 OD
16,400-32,700 762 OD
5.6 MEASUREMENTS
The adequate evaluation of a pumping test relies very much upon the recording of a number of sets
of measurements throughout the test. These include measurements of:
1. Time.
2. Discharge.
3. Water level or Head.
4. Temperature.
5. Water quality.
To analyse a pumping test accurately each set of measurements must include the recording time,
head and discharge. The values obtained for the aquifer characteristics cannot be of greater accuracy
than that of the basic data. Care should be taken then in the measuring and recording of time
discharge and head. For hot flowing bores measurements of temperature should be recorded as well.
Ideally, the natural fluctuation in hydraulic head of the aquifer should be known before the test
commences (e.g. hydrographs). This information can be used to correct the drawdown observed
during the test.
In coastal aquifers where the hydraulic head is affected by tidal movements, a complete hydrograph
should be obtained, including maximum and minimum levels.
For long-term tests (days), the levels of near-by surface waters and any precipitation should also be
recorded and, if the aquifer being tested is a confined, the barometric pressure should recorded as
well.
Remeasurement 24 2
Town water supply and industrial As for irrigation bores but 100 hour duration
5.11 DISINFECTION
Because of the risk cross contamination from one aquifer to another the down hole pumping
equipment should be disinfected thoroughly before and after a pumping test. This is normally done
with a chlorine solution to kill any bacteria which may be present. The details of disinfection will not
be covered in these documents but can be found in references such as Driscoll, 1988.
Separating variables:
dr 2 b
Kdh
r Q
Integrating between r1 and r2 and h1 and h2:
r2 h
dr 2 b 2
dh
r1
r Q h1
r2 2 Kb h2
log e r1
h h1
Q
r2 2 Kb
log e h2 h1
r1 Q
Hence:
2.3 log10 r2 / r1
K .....6.3
2 b h2 h1
or
2.3 log10 r2 / r1
T Kb .....6.4
2 h2 h1
Because:
h1 +s1 = h2 +s2
then:
s1 - s2 = h2 - h1
and
2.3Q log10 r2 / r1
T
2 s1 s 2 .....6.5
If readings are taken in two observation bores during a pumping test on a bore then equation 6.5 may
be used to determine the transmissivity of the aquifer. However, a straight line solution does exist
which will allow the computation of transmissivity by graphical means and incorporates the
drawdowns in any number of observation bores.
y2 y1 m( x2 x1 )
This is an equation to a straight line with slope m.
Thus a plot of drawdown “s” against “log 10r” under steady state conditions will result in a straight line
with slope.
2.3Q
2 T
If, on the plot, r2 is taken as 10r1 then log10r2 - log10r1 = 1 and
2.3Q
s2 s1
2 T
or 2.3Q
s1 s2
2 T
= Δs' .....6.7
where Δs' is the change in drawdown per log cycle from the distance-drawdown plot.
From Equation 6.7:
2.3Q
s'
2 T
and
2.3Q
T
2 s' .....6.8
The slope is negative because the drawdown, "s", decreases as distance "r" increases, or, stating in
another way, if r2 is greater than r1, in the above expression then s2 - s1 is always negative. T will
always be positive.
This is a very convenient method of determining transmissivity under Steady State conditions.
Procedure
1. On semi-logarithmic graph paper plot the drawdowns (on natural scale) against the
radial distance from the pumping bore at which the drawdowns were measured (on
logarithmic scale).
2. Fit a straight line through the plotted points.
3. From the plot read the drawdown per log cycle, Δs'.
4. Using equation 6.8 calculate T, from T = 2.3Q
2 s'
Assumptions
In deriving the Thiem equation, the following assumptions were made:
1. the aquifer is homogeneous, isotropic, and of infinite areal extent;
2. the discharging bore penetrates and receives water from the entire thickness of the
aquifer;
3. the transmissivity is constant at all places and all times;
4. discharge has continued at a constant rate for a time sufficient for the hydraulic system
to reach a steady flow condition;
5. flow lines are radial; and
6. flow is laminar.
1.20
S-2
Δ´ = 1.28m
1.40
N-2
1.60
S-1
1.80
N-1
s=1.96 m
2.00
1.00 10.00 100.00 1000.00
Distance" r" (m)
Q1 Ti1W1
h
T2 r
r .....6.10
where:
T = transmissivity. h
i = hydraulic gradient (at the inner face) = r
W = the width of flow across the section (at the inner face equals the circumference of the
circle = 2πr).
The rate of change of volume within the cylindrical shell is expressed as:
dV h
2 rdr S
dt t .....6.13
where:
S = coefficient of storage.
For unconfined aquifer (or water table) conditions, S is equivalent to the specific yield of the material
dewatered by pumping. For pressure conditions where water is drawn from storage by compression of
the aquifer, S is the storage coefficient. Specific yield and storage coefficient have been defined
previously.
Substituting the above values in equation 6.9:
2
h h h h
T 2 r T( 2
dr).2 (r dr) 2 rdr S
r r r t
2 2
h h h h h h
T 2 r T (2 r 2 rdr 2
2 dr 2 (dr) 2 2
) 2 rdr S
i.e. r r r r r t
Dividing through by 2πrTdr and neglecting differentials higher than first order:
2
h 1 h S h
2
r r r T t .....6.14
Equation 6.14 is the partial differential equation for non-steady state radial flow.
This equation may be written in cartesian co-ordinates as:
2 2 2
h h h S h
2 2 2
x y z T t
h
Note: when, 0 the entire right hand side of equation 6.14 is zero indicating that there are no
t
changes in storage in the aquifer and steady state flow occurs.
The equation for steady state radial flow is then:
2
1 h h
0 .....6.15
r r r2
e x 1 x x2 x3
1
x x 2! 3! 4!
e x x2 x3 x4
dx loge x ........
u
x 2.2! 3.3! 4.4! u
u2 u3
( 0.577216 loge u u .....)
2.2! 3.3!
and
Q u2 u3 .....6.17
s ( 0.577216 loge u u ......)
4 T 2.2! 3.3!
Q e x
s dx
4 Tu x
Q
s W (u )
4 T
where:
r 2S
u
4Tt
GROUNDWATER HYDRAULICS Page 69 of 225
The value of the series is commonly expressed as W(u), the "well function of "u".
Values of W(u) for values of u from 10-15 to 9.9 are given in tabular form in Appendix 6-A. Some of
these values are presented in graphical form in Figure 6-4. However, if the value of W(u) is required
then in most cases sufficient accuracy is obtained by calculating it using the first two terms of the
series (i.e. W(u) ≈ -0.577216 -logeu). For a given value of u, T may be determined from:
Q
T W (u )
4 s .....6.18
And S may be determined from:
S = 4Ttu/r2 .....6.19
Theis Type Curve Solution
Theis devised a simple graphical method of superposition that makes it possible to obtain solutions of
equation 6.18 and 6.19.
By re-arranging equations 6.18 and 6.19 and taking logarithms of both sides, the following
relationships are obtained:
log Q
log s logW (u ) .....6.20
4 T
and r2 log 4T
log log u
t S .....6.21
If the discharge is held constant, the bracketed parts of the equations 6.20 and 6.21 are constant,
and s is related to r2/t in the same manner as W (u) is related to u. If the values of W (u) and u were
plotted on logarithmic paper, a type curve for the relationship between s and r 2/t would result. Values
of “s” could then be plotted against r 2/t on transparent logarithmic paper to the same scale as the
type curve and would be the same shape as the type curve, but would be displaced by an amount
Q/4πT on the “s” and “W (u)” axes, and by an amount 4T/S on the " r2/t " and “u” axes.
The plot of “s” versus r2/t could be moved over the type curve, keeping the axes parallel, until its
shape was matched with a section of the type curve.
For this matching position, corresponding values of “W (u)” and “u”, "r2/t" and “s” are recorded for
any point on the graphs. The point selected is called the “Match Point”. This match point does not
have to lie on the plotted curve. In fact, for convenience of calculations, a match point is frequently
chosen so that “W (u)” = 1 and “u” = 1 and the values of s and r2/t are read.
Transmissivity and storage coefficient are then calculated from equations 6.18 and 6.19 making sure
that the units used are consistent.
If readings from only one observation bore are used then “s” could be plotted against “1/t” and “r 2”
would be introduced in equation 6.19. This eliminates the necessity for calculating many values of r2/t.
However, it is recommended that the alternative method presented below be used in practice.
Alternative Type Curve Solution of Equation
The above method involves the calculation of the reciprocal for every "time" measurement, and apart
from the extra calculations, additional sources of error are introduced because of these calculations.
A more convenient method is to plot "s" against "t" on logarithmic paper, and superimpose this on a
type curve. This can be done for each observation bore. It is also possible to plot "s" against “t/r 2” for
all observation bores.
Following the same reasoning as previously, “s” has the same relationship to “t/r 2” as the well function
“W (u)” has to “1/u”.
The type curve in this case is then a plot of W (u) against 1/u on logarithmic paper. A stylised curve of
W (u) against 1/u is given in Figure 6-4.
The uppermost curve in this family of curves is in fact the curve of "W(u)" versus “I/u” (i.e. no
leakage), or the Theis curve. If a deviation from the W (u) versus I/u curve does occur it can be
interpreted from the same type curve if the deviation is caused by leakiness.
GROUNDWATER HYDRAULICS Page 70 of 225
Procedure
Summarising the type curve solution, the procedure to be followed is:
1. Plot drawdown versus time for the observation bore, on transparent logarithmic paper
with time on the horizontal axis.
2. Fit the plotted curve to the type curve of W (u) versus 1/u (Figure 6-4).
3. Select a match-point and read off values of W (u), (1/u), s and t,
4. Calculate the transmissivity, T, and storage coefficient, S, from equations 6.18 and
6.19.
5. Before calculating storage coefficient, S, from equation 6.19, ensure that the units of T
and t are such that S will be dimensionless.
The scale of the logarithmic paper on which drawdown and time are plotted must of course be the
same as the scale of the type curve paper.
If the recorded time units are not consistent with the required time units in equation 6.19 plot the
recorded time values and make one only conversion to the match-point value when solving equation
6.19.
N\u Nx1O-15 Nx1O-14 Nx1O-13 Nx1O-12 Nx1O-11 Nx1O-10 Nx1O-9 Nx1O-8 Nx1O-7 Nx1O-6 Nx1O-5 Nx1O-4 Nx1O-3 Nx1O-2 Nx10-1 N
1.0 33.9615 31.6590 29.3564 27.0518 24.7512 22.4486 20.1460 17.8435 15.5409 13.2383 10.9357 8.6332 6.3315 4.0379 1.8229 0.2194
1.1 33.8662 31.5637 29.2611 26.9585 24.6559 22.3533 20.0507 17.7482 15.4456 13.1430 10.8404 8.5379 6.2363 3.9436 1.7371 0.1860
1.2 33.7792 31.4767 29.1741 26.8715 24.5689 22.2663 19.9637 17.6611 15.3586 13.0560 10.7534 8.4509 6.1494 3.8576 1.6595 0.1584
1.3 33.6992 31.3966 29.0940 26.7914 24.4889 22.1863 19.8837 17.5811 15.2785 12.9759 10.6734 8.3709 6.0695 3.7785 1.5889 0.1355
1.4 33.6251 31.3225 29.0199 26.7173 24.4147 22.1122 19.8096 17.5070 15.2044 12.9018 10.5993 8.2968 5.9955 3.7054 1.5241 0.1162
1.5 33.5561 31.2535 28.9509 26.6483 24.3458 22.0432 19.7406 17.4380 15.1354 12.8328 10.5303 8.2278 5.9266 3.6374 1.4645 0.1000
1.6 33.4916 31.1890 28.8864 26.5838 24.2812 21.9786 19.6760 17.3735 15.0709 12.7683 10.4657 8.1634 5.8621 3.5739 1.4092 0.0863
1.7 33.4309 31.1283 28.8258 26.5232 24.2206 21.9180 19.6154 17.3128 15.0103 12.7077 10.4051 8.1027 5.8016 3.5143 1.3578 0.07465
1.8 33.3738 31.0712 28.7686 26.4660 24.1634 21.8608 19.5583 17.2557 14.9531 12.6505 10.3479 8.0455 5.7446 3.5481 1.3080 0.06471
1.9 33.3197 31.0171 28.7145 26.4119 24.1094 21.8068 19.5042 17.2016 14.8990 12.5964 10.2939 7.9915 5.6906 3.4050 1.2649 0.05620
2.0 33.2684 30.9658 28.6632 26.3607 24.0581 21.7555 19.4529 17.1503 14.8477 12.5451 10.2426 7.9402 5.6394 3.3547 1.2227 0.04890
2.1 33.2196 30.9170 28.6145 26.3119 24.0093 21.7067 19.4041 17.1015 14.7989 12.4964 10.1938 7.8914 5.5907 3.3069 1.1829 0.04261
2.2 33.1731 30.8705 28.5679 26.2653 23.9628 21.6602 19.3576 17.0550 14.7324 12.4498 10.1473 7.8449 5.5443 3.2614 1.1454 0.03719
2.3 33.1286 30.8261 28.5235 26.2209 23.9183 21.6157 19.3131 17.0106 14.7080 12.4054 10.1028 7.8004 5.4999 3.2179 1.1099 0.03250
2.4 33.0861 30.7835 28.4809 26.1783 23.8758 21.5732 19.2706 16.9680 14.6654 12.3628 10.0603 7.7579 5.4575 3.1763 1.0762 0.02844
2.5 33.0453 30.7427 28.1401 26.1375 23.8349 21.5323 19.2298 16.9272 14.6246 12.3220 10.0194 7.7172 5.4167 3.1365 1.0443 0.02491
2.6 33.0060 30.7035 28.1009 26.0983 23.7957 21.4931 19.1905 16.8880 14.5854 12.2828 9.9802 7.6779 5.3776 3.0983 1.0139 0.02185
2.7 32.9683 30.6657 28.3631 26.0806 23.7580 21.4554 19.1528 16.8502 14.5476 12.2450 9.9425 7.6401 5.3400 3.0615 0.9849 0.01918
2.8 32.9519 30.6294 28.3268 26.0242 23.7216 21.4190 19.1164 16.8138 14.5113 12.2087 9.9061 7.6038 5.3037 3.0261 0.9573 0.01686
2.9 32.8968 30.5943 28.2917 25.9891 23.6865 21.3839 19.0813 16.7788 14.4762 12.1736 9.8710 7.5687 5.2687 2.9920 0.9309 0.01482
3.0 32.8629 30.5604 28.2578 25.9552 23.6526 21.3500 19.0474 16.7440 14.4423 12.1397 9.8371 7.5348 5.2349 2.9591 0.9057 0.01305
3.1 32.8302 30.5276 28.2250 25.9224 23.6198 21.3172 19.0146 16.7121 14.4095 12.1069 9.8043 7.5020 5.2022 2.9273 0.8815 0.01149
3.2 32.7984 30.4958 28.1932 25.8907 23.5880 21.2855 18.9829 16.6803 14.3777 12.0751 9.7726 7.4703 5.1706 2.8965 0.8583 0.01013
3.3 32.7676 30.4651 28.1625 25.8599 23.5573 21.2547 18.9521 16.6495 14.3470 12.0444 9.7418 7.4395 5.1399 2.8668 0.8361 0.008939
3.4 32.7378 30.4352 28.1326 25.8300 23.5274 21.2249 18.9223 16.6197 14.3171 12.0145 9.7120 7.4097 5.1102 2.8377 0.8147 0.007891
3.5 32.7088 30.4062 28.1036 25.8010 23.4985 21.1959 18.8933 16.5907 14.2881 11.9855 9.6830 7.3807 5.0813 2.8099 0.7942 0.006970
3.6 32.6806 30.3780 28.0755 25.7729 23.4703 21.1677 18.8651 16.5625 14.2599 11.9574 9.6548 7.3526 5.0532 2.7827 0.7745 0.006160
3.7 32.6532 30.3506 28.0481 25.7455 23.4429 21.1403 18.8377 16.5351 14.2325 11.9300 9.6274 7.3252 5.0259 2.7563 0.7554 0.005448
3.8 32.6266 30.3240 28.0214 25.7188 23.4162 21.1136 18.8110 16.5085 14.2059 11.9033 9.6007 7.2985 4.9933 2.7306 0.7371 0.004820
3.9 32.6006 30.2980 27.9954 25.6928 23.3902 21.0877 18.7851 16.4825 14.1799 11.8773 9.5748 7.2725 4.9735 2.7056 0.7194 0.004267
4.0 32.5753 30.2727 27.9701 25.6675 23.3649 21.0623 18.7598 16.4572 14.1546 11.8520 9.5495 7.2472 4.9482 2.6813 0.7024 0.003779
4.1 32.5506 30.2480 27.9454 25.6428 23.3402 21.0376 18.7351 16.4325 14.1299 11.8273 9.5248 7.2225 4.9236 2.6576 0.6859 0.003349
4.2 32.5265 30.2239 27.9213 25.6187 23.3161 21.0136 18.7110 16.4084 14.1058 11.8032 9.5007 7.1985 4.8997 2.6344 0.6700 0.002969
4.3 32.5029 30.2004 27.8978 25.5952 23.2926 20.9900 18.6874 16.3884 14.0823 11.7797 9.4771 7.1749 4.8762 2.6119 0.6546 0.002633
4.4 32.4800 30.1774 27.8748 25.5722 23.2696 20.9670 18.6644 16.3619 14.0593 11.7567 9.4541 7.1520 4.8533 2.5899 0.6397 0.002336
4.5 32.4575 30.1549 27.8523 25.5497 23.2471 20.9446 18.6420 16.3394 14.0368 11.7342 9.4317 7.1295 4.8310 2.5684 0.6253 0.002073
4.6 32.4355 30.1329 27.8303 25.5277 23.2252 20.9226 18.6200 16.3174 14.0148 11.7122 9.4097 7.1075 4.8091 2.5474 0.6114 0.001841
4.7 32.4140 30.1114 27.8088 25.5602 23.2037 20.9011 18.5985 16.2959 13.9933 11.6907 9.3882 7.0860 4.7877 2.5268 0.5979 0.001635
4.8 32.3929 30.0904 27.7878 25.4852 23.1826 20.8800 18.5774 16.2748 13.9723 11.6697 9.3671 7.0650 4.7667 2.5068 0.5848 0.001453
4.9 32.3723 30.0697 27.7672 25.4646 23.1620 20.8594 18.5568 16.2542 13.9516 11.6491 9.3465 7.0444 4.7462 2.4871 0.5721 0.001291
5.0 32.3521 30.0495 27.7470 25.4444 23.1418 20.8392 18.5366 16.2340 13.9314 11.6289 9.3263 7.0242 4.7261 2.4679 0.5598 0.001148
5.1 32.3323 30.0297 27.7271 25.4246 23.1220 20.8194 18.5168 16.2142 13.9116 11.6091 9.3065 7.0044 4.7064 2.4491 0.5478 0.001021
5.2 32.3129 30.0103 27.7077 25.4051 23.1026 20.8000 18.4974 16.1948 13.8922 11.5896 9.2871 6.9850 4.6871 2.4306 0.5362 0.0009086
5.3 32.2939 29.9913 27.6887 25.3861 23.0835 20.7809 18.4783 16.1758 13.8732 11.5706 9.2681 6.9659 4.6681 2.4126 0.5250 0.0008086
5.4 32.2752 29.9726 27.6700 25.3674 23.0648 20.7622 18.4596 16.1571 13.8545 11.5519 9.2494 6.9473 4.6495 2.3948 0.5140 0.0007198
5.5 32.2568 29.9542 27.6516 25.3491 23.0465 20.7439 18.4413 16.1387 13.8361 11.5336 9.2310 6.9289 4.6313 2.3775 0.5034 0.0006409
5.6 32.2388 29.9362 27.6336 25.3310 23.0285 20.7259 18.4233 16.1207 13.8181 11.5155 9.2130 6.9109 4.6134 2.3604 0.4930 0.0005708
5.7 32.2211 29.9185 27.6159 25.3133 23.0108 20.7082 18.4056 16.1030 13.8004 11.4978 9.1953 6.8932 4.5958 2.3437 0.4830 0.0005085
5.8 32.2037 29.9011 27.5985 25.2959 22.9934 20.6908 18.3882 16.0856 13.7830 11.4804 9.1779 6.8758 4.5785 2.3273 0.4732 0.0004532
5.9 32.1866 29.8840 27.5814 25.2789 22.9763 20.6737 18.3711 16.0685 13.7659 11.4633 9.1608 6.8588 4.5615 2.3111 0.4637 0.0004039
6.0 32.1698 29.8672 27.5646 25.2620 22.9595 20.6569 18.3543 16.0517 13.7491 11.4465 9.1440 6.8420 4.5448 2.2953 0.4544 0.0003601
6.1 32.1533 29.8507 27.5481 25.2455 22.9429 20.6403 18.3378 16.0352 13.7326 11.4300 9.1275 6.8254 4.5283 2.2797 0.4454 0.0003211
6.2 32.1370 29.8344 27.5318 25.2293 22.9267 20.6241 18.3215 16.0189 13.7163 11.4138 9.1112 6.8092 4.5122 2.2645 0.4366 0.0002864
6.3 32.1210 29.8184 27.5158 25.2133 22.9107 20.6081 18.3055 16.0029 13.7003 11.3978 9.0952 6.7932 4.4963 2.2494 0.4280 0.0002555
6.4 32.2053 29.8027 27.5001 25.1975 22.8949 20.5923 18.2898 15.9872 13.6846 11.3820 9.0795 6.7775 4.4806 2.2346 0.4197 0.0002279
6.5 32.0898 29.7872 27.4846 25.1820 22.8794 20.5768 18.2742 15.9717 13.6691 11.3665 9.0640 6.7620 4.4652 2.2201 0.4115 0.0002034
6.6 32.0745 29.7719 27.4693 25.1667 22.8641 20.5616 18.2590 15.9564 13.6538 11.3512 9.0487 6.7467 4.4501 2.2058 0.4036 0.0001816
6.7 32.0595 29.7569 27.4543 25.1517 22.8491 20.5465 18.2439 15.9414 13.6388 11.3362 9.0337 6.7317 4.4351 2.1917 0.3959 0.0001621
6.8 32.0446 29.7421 27.4395 25.1369 22.8343 20.5317 18.2291 15.9265 13.6240 11.3214 9.0189 6.7169 4.4204 2.1779 0.3883 0.0001448
6.9 32.0300 29.7275 27.4249 25.1223 22.8197 20.5171 18.2145 15.9119 13.6094 11.3608 9.0043 6.7023 4.4059 2.1632 0.3810 0.0001293
7.0 32.0156 29.7131 27.4105 25.1079 22.8053 20.5027 18.2001 15.8976 13.5950 11.2924 8.9899 6.6879 4.3916 2.1508 0.3738 0.0001135
7.1 32.0015 29.6989 27.3963 25.0937 22.7911 20.4885 18.1860 15.8834 13.5808 11.2782 8.9757 6.6737 4.3775 2.1376 0.3668 0.0001032
7.2 31.9875 29.6849 27.3823 25.0797 22.7771 20.4746 18.1720 15.8694 13.5668 11.2642 8.9617 6.6598 4.3636 2.1246 0.3599 0.00009219
7.3 31.9737 29.6711 27.3685 25.0659 22.7633 20.4603 18.1582 15.8556 13.5530 11.2504 8.9479 6.6460 4.3500 2.1118 0.3532 0.00008239
7.4 31.9601 29.6575 27.3549 25.0523 22.7497 20.4472 18.1446 15.8420 13.5394 11.2368 8.9343 6.6324 4.3364 2.0991 0.3467 0.00007364
7.5 31.9467 29.6441 27.3415 25.0389 22.7363 20.4337 18.1311 15.8286 13.5260 11.2234 8.9209 6.6190 4.3231 2.0867 0.3403 0.00006583
7.6 31.9334 29.6308 27.3282 25.0257 22.7231 20.4205 18.1179 15.8153 13.5127 11.2102 8.9076 6.6057 4.3100 2.0744 0.3341 0.00005886
7.7 31.9203 29.6173 27.3152 25.0126 22.7100 20.4074 18.1048 15.8022 13.4997 11.1971 8.8946 6.5927 4.2970 2.0623 0.3280 0.00005263
7.8 31.9074 29.6048 27.3023 24.9997 22.6971 20.3945 18.0919 15.7893 13.4868 11.1842 8.8817 6.5798 4.2842 2.0503 0.3221 0.00004707
7.9 31.8947 29.5921 27.2895 24.9869 22.6844 20.3818 18.0792 15.7766 13.4740 11.1714 8.8689 6.5671 4.2716 2.0386 0.3163 0.00004210
8.0 31.8821 29.5795 27.2769 24.9744 22.6718 20.3692 18.0666 15.7640 13.4614 11.1589 8.8563 6.5545 4.2591 2.0269 0.3106 0.00003767
8.1 31.8697 29.5671 27.2645 24.9619 22.6594 20.3568 18.0542 15.7516 13.4490 11.1464 8.8439 6.5421 4.2468 2.0155 0.3050 0.00003370
8.2 31.8574 29.5548 27.2523 24.9497 22.6471 20.3445 18.0419 15.7393 13.4367 11.1342 8.8317 6.5298 4.2346 2.0042 0.2996 0.00003015
8.3 31.8453 29.5427 27.2401 24.9375 22.6350 20.3324 18.0298 15.7272 13.4246 11.1220 8.8195 6.5177 4.2226 1.9930 0.2943 0.00002699
8.4 31.8333 29.5307 27.2282 24.9256 22.6230 20.3204 18.0178 15.7152 13.4126 11.1101 8.8076 6.5057 4.2107 1.9820 0.2891 0.00002415
8.5 31.8215 29.5189 27.2163 24.9137 22.6112 20.3086 18.0060 15.7034 13.4008 11.0982 8.7957 6.4939 4.1990 1.9711 0.2840 0.00002162
8.6 31.8098 29.5072 27.2046 24.9020 22.5695 20.2969 17.9943 15.6917 13.3891 11.0865 8.7840 6.4822 4.1874 1.9604 0.2790 0.00001936
8.7 31.7982 29.4957 27.1931 24.8905 22.5879 20.2853 17.9827 15.6801 13.3776 11.0750 8.7725 6.4707 4.1759 1.9498 0.2742 0.00001733
8.8 31.7868 29.4842 27.1816 24.8790 22.5765 20.2739 17.9713 15.6687 13.3661 11.0635 8.7610 6.4592 4.1646 1.9393 0.2694 0.00001552
8.9 31.7755 29.4729 27.1703 24.8678 22.5652 20.2626 17.9600 15.6574 13.3548 11.0523 8.7497 6.4480 4.1534 1.9290 0.2647 0.00001390
9.0 31.7643 29.4618 27.1582 24.8566 22.5540 20.2514 17.9488 15.6462 13.3437 11.0411 8.7386 6.4368 4.1423 1.9187 0.2602 0.00001245
9.1 31.7533 29.4507 27.1481 24.8455 22.5429 20.2404 17.9378 15.6352 13.3326 11.0300 8.7275 6.4258 4.1313 1.9087 0.2557 0.00001115
9.2 31.7424 29.4398 27.1372 24.8346 22.5320 20.2294 17.9268 15.6243 13.3217 11.0191 8.7166 6.4148 4.1205 1.8987 0.2513 0.000009988
9.3 31.7315 29.4290 27.1264 24.8238 22.5212 20.2186 17.9160 15.6135 13.3109 11.0083 8.7058 6.4040 4.1098 1.8888 0.2470 0.000008948
9.4 31.7208 29.4183 27.1157 24.8131 22.5105 20.2079 17.9053 15.6028 13.3002 10.9976 8.6951 6.3934 4.0992 1.8791 0.2429 0.000008018
9.5 31.7103 29.4077 27.1051 24.8025 22.4999 20.1973 17.8948 15.5922 13.2896 10.9870 8.6845 6.3828 4.0887 1.8595 0.2387 0.000007185
9.6 31.6998 29.3972 27.0946 24.7920 22.4895 20.1869 17.8843 15.5817 13.2791 10.9765 8.6740 6.3723 4.0784 1.8599 0.2347 0.000006439
9.7 31.6894 29.3868 27.0843 24.7817 22.4791 20.1765 17.8739 15.5713 13.2688 10.9662 8.6637 6.3620 4.0681 1.8505 0.2308 0.000005771
9.8 31.6792 29.3766 27.0740 24.7714 22.4688 20.1663 17.8637 15.5611 13.2585 10.9559 8.6534 6.3517 4.0579 1.8412 0.2269 0.000005173
9.9 31.6690 29.3664 27.0639 24.7613 22.4587 20.1561 17.8535 15.5509 13.2483 10.9458 8.6433 6.3416 4.0479 1.8320 0.2231 0.000004637
Time since pumping started Bore N-1 Bore N-2 Bore N-3
r = 61 m r = 122 m r = 244 m
1 Drawdown (m)
W(u)
1.0
0.1
Match Point
W(u) = 1, 1/u = 10
s = 0.17m, t = 8.2 mins
(5.7 x 10-3 days
"1/u"
0.1
1.0 10.0 100.0 1000.0
0.01
Time (mins)
1 10 100 1000
Figure 6-5 Type curve solution, confined aquifer, non-steady state, constant Q
The alternative type curve solution is shown is Figure 6-5, as a plot of drawdown versus time (in
minutes) on logarithmic paper. There is no need to convert the time units from minutes to days
before plotting. The conversion can be done prior to the calculation of S from equation 6.19. The
analysis has been carried out for one bore only (bore N-2). Plots for N-1 and N-3 could be carried out
separately or plotted as one curve as t/r 2.
Solution
From equation 6.18:
Q W (u )
T
4 s
Q = 2720 m3/day
From match point:
W(u) = 1, s = 0.17m
2720x1
T
4 0.17
= 1270 m2/day
From equation 6.19:
4Tut
S
r2
r2 = 1.5 x 104 m2
From match point:
T = 5.7 x 10-3 days, u = 0.1
S = 1.93 x 10-4
Q r 2S
s ( 0.577216 loge )
4 T 4Tt .....6.22
This may be rewritten as:
Q 4Tt
T (loge 0.562 loge 2 S )
4 s r S
2.3Q 2.25Tt
(log10 2 )
4 s r S .....6.23
If a series of drawdown measurements are taken at one observation bore at various times during a
constant discharge pump test, then all terms in equation 6.23 are constant except "s" and "t".
Equation 6.23 can be rearranged to give:
2.3Q 2.25T
s (log10 2 log10 t )
4 T r S
aQ bQ log10 t .....6.24
where:
2.3 s
b = constant
4 T Q
2.25T
a b log10 2 = constant
r S
A plot of drawdown "s" versus "log10t" will give a straight line with slope "bQ" and, provided non linear
head loss is absent, an intercept of "aQ".
“bQ” is called “Δs” and is defined as the drawdown per log cycle.
2.25Tt 0
S
r 2 100
2.25Tt 0
r2 .....6.29
t0 is determined by extending the "time-drawdown" curve back until it intersects the zero drawdown
line.
Procedure
1. Plot on semi-logarithmic graph paper, drawdown versus time, with time on the
logarithmic scale.
2. Calculate Δs, the drawdown per log cycle,
3. Transmissivity is determined using equation 6.26.
4. Storage Coefficient is determined either from equation 6.28 or 6.29.
5. In equation 6.29, "t0" is obtained by extrapolating the straight line plot back to intersect
the zero drawdown line.
Note: equation 6.26 can be used to determine transmissivity using data from either the pumped bore
or observation bore. The presence of non linear head loss does not affect the slope on a semi-log plot,
merely its location.
2.25Tt
D log10 constant
S
In practice, drawdown “s” is plotted against “r” on semi-logarithmic graph paper.
If the slope of the distance drawdown plot is called Δs' , then:
Δs' = -2Δs
= -2 x 2.3 Q/4πT
and
T = -2.3 Q/2π Δs' .....6.31
Δs’ is always negative so T is always positive.
In the same manner as was used to derive equation 6.28:
2.25Tt
S
r 10( 2 s / s ')
2
.....6.32
Transmissivity and storage coefficient are then calculated using equations 6.31 and 6.32 respectively.
Procedure
1. Plot on semi-logarithmic graph paper, drawdown in the observation bores, at time "t",
versus the distance of the observation bores from the pumping facility "r". The distance
should be plotted on the logarithmic scale.
2. Determine the drawdown per log cycle Δs'.
3. Calculate the Transmissivity T from equation 6.31.
4. Calculate the Storage Coefficient from equation 6.32.
1.20
1 10 100 1000
Time "t" (mins)
Figure 6-6 Modified non-steady state flow example – confined aquifer, constant Q, constant r, varying t
2.3Q1 2.25Tt1
s1 log10 S (discharging at Q1) .....6.33
4 T r2
and
2.3Q2 2.25Tt 2
s2 log10
4 T r 2 S (recharging at Q2) .....6.34
where:
t1 = time since pumping began.
t2 = time since pumping stopped.
To obtain a rise in water level during pumping it is not necessary to have a complete cessation of
pumping. Consider initially the case where the rate of recharge from the recharging bore is equal to
the rate of discharge from the discharging bore.
Let Q1 = Q = Q
Residual drawdown = s1 + (-s2)
2.3Q t
log10 1
4 T t2 .....6.35
0.90
0.95 m
1.10
1.30
10 100 1000 10000
Distance (m)
Figure 6-7 Modified non-steady state flow example – confined aquifer, constant Q, constant t, varying r
If the recharge to the bore, i.e. the reduction in discharge rate, is given as Q2, then the recovery in
the bore will be given by equation 6.34 and the transmissivity given by:
2.3Q2 .....6.37
T
4 s( re cov ery )
r 2S
0.05
4Tt
or
0.2Tt
S
r2
From the Modified Non-Steady State flow equation the storage coefficient is given by:
2.25Tt
S
r 210s / s
If the Modified Non-Steady State flow equations are to be applicable:
2.25Tt 0.2Tt
r 210s / s r2
i.e.
2.25 / 10(s/Δs) ≤ 0.2
or
10(s/Δs) ≥ 11.25
or
(s/Δs) ≥ 1.05 .....6.40
For the Modified Non-Steady State Flow equations to be applicable the drawdown in the observation
bore must be greater than 1.05 Δs.
This is a very quick and easy check on the applicability of the Modified Non-Steady State Flow
Equations, for radial flow in confined aquifers. Occasions will arise where the straight line
approximation is valid even where the drawdown is less than 1.05 Δs. However, this will be apparent
in the plot, and it will be permissible to use the Modified Non-Steady State Flow Equations.
2.3Q 1.12L
sm (log10 )
2 T r .....6.44
where:
sm = maximum (or steady state) drawdown at distance r from the pumping bore.
Q = discharge rate.
T = the transmissivity of the aquifer.
L = the leakage factor of semi-pervious material = √(Tc).
r = distance of observation bore from pumping bore.
Thus a plot of "s" versus "r" on semi-logarithmic paper, with "r" on the logarithmic scale, will result in
a straight line in the range where r/L is small.
In the range where r/L is large, the points fall on a curve that approaches the zero drawdown axis
asymptotically.
The slope of the straight line portion of the curve, i.e. the drawdown per log cycle of r, is expressed
as:
2.3Q
sm .....6.45
2 T
(r0 / 1.12) 2
c
T .....6.47
Procedure
1. Plot on semi-logarithmic paper "sm" versus "r" ("r" on the logarithmic scale), i.e. the
maximum or steady state drawdown measured in each observation bore against the
corresponding distance of the observation bore from the pumped bore.
2. Draw the straight line of best fit through the points which appear to fall on a straight
line.
3. Determine the slope of the line of best fit, i.e. the drawdown per log cycle.
4. Using the discharge rate Q calculate T from equation 6.45.
5. Extend the straight line to intersect the r axis and record this intercept as r o.
6. Calculate the hydraulic resistance "c" of the semi-pervious layer by substituting r and T
into equation 6.47. The leakage coefficient can be calculated from the reciprocal of the
hydraulic resistance.
Table 6-4 Data for semi-confined aquifer test analysis
Time Since Pumping Began Drawdown in metres
Example
Table 6-4 (taken from Cooper, 1963 p. C54) gives postulated drawdowns in observation bores at
various distances from a bore which discharges at a constant rate of 5 450 m 3 /day for 1,000 minutes
from a leaky artesian aquifer.
Determine transmissivity and hydraulic resistance for the aquifer utilising equation 6.45 and the
drawdown data at 1,000 minutes.
The semi-logarithmic plot and analysis are given in Figure 6-8.
0.5 c = ((r0/1.12)2 ) / T
From the plot:
Δsm = 1.57m r0 = 750m
Hence:
1
c = 353 days
Leakage Coefficient: (see note)
Leakage Coefficient =1/c
1.5 = 2.83 x 10-3 days
1.37 m NOTE: the values derived for c and 1/c are
approximate only. One assumption for the straight line
solution is r/L < 0.05.
2 From equation 3.15:
L = √(Tc)
= 670m
2.5 i.e. only r values < 0.05 x 670 (i.e. 34m) can be used.
More reliance should be placed on the Non-Steady
10 100 1000 State solution given in Figure 6.9
Distance (m)
Q L(u, v)
T
4 s
6. Compute storage coefficient from:
4Tut
S
r2 .....6.49
7. Compute hydraulic resistance from:
C = b'/K'
= r2/4Tv2 .....6.50
If required the leakage coefficient could have been calculated as:
1/c = 4Tv2/r2
8. Make sure that the units used are consistent.
Example
Using the data presented in Table 6-4 for bore 3; determine the transmissivity storage coefficient,
hydraulic resistance and leakage coefficient for the aquifer.
The plot and type curve analysis are presented in Figure 6-9.
The test data is plotted on double logarithmic paper and overlaid on Figure 6-4.
Assumptions
The assumptions made in deriving the Non-Steady State flow equations presented above are:
1 to 6 from section 6.4.2.1
7. The flow to the bore is in a non-steady state i.e. the drawdown differences with time
are not negligible nor is the hydraulic gradient constant with time.
8. The water removed from storage is discharged instantaneously with decline in head.
9. The bore diameter is very small, so the storage in the bore can be neglected.
dr K
2 hdh
r Q
K (b s 2 ) 2 (b s1 ) 2
Q .....6.56
log e (r2 / r1 )
If the right hand side is expanded and multiplied by unity in the form 2b/2b:
s12 s 22
2 Kb ( s1 ) (s 2 )
2b 2b
Q .....6.57
log e (r2 / r1 )
2 T ( s1' s 2' )
Q
loge (r2 / r1 ) .....6.58
where:
s12 s 22
s1' (s1 ) s 2' (s2 )
2b and 2b
It is now in the same form as the Thiem equation for radial flow to a bore in a confined aquifer,
equation 6.5, except that:
s12 and s22 are written instead of s1 and s2 respectively.
( s1 ) (s2 )
2b 2b
Transmissivity
Using the corrected values for drawdowns the transmissivity can be calculated in the same manner as
was used for the Thiem equation for Steady State flow confined aquifers, Section 6.2.1.1, and the
various equations for Non-Steady State flow.
Storage Coefficient
The apparent value of storage co-efficient can be determined using corrected values of drawdown and
transmissivity.
S-2
1.20
N- 2
Δs' = 1.04 m
1.40
S-1
1.60
N- 1
1.74 m
1.80
2.00
1.00 10.00 100.00 1000.00
Distance (metres)
r 2S A
uA
4 Kbt .....6.64
Kb
B .....6.68
SY
T
B .....6.69
SY
0 0 41 0.128
Table 6-5 presents data collected at an observation bore during pumping test at "Vennebulton".
The drilling log and the performance of a shallow observation bore during the test Indicated that the
test should be analysed using the Boulton Method.
The bore was pumped for 25 hours at 873 m 3/day.
The observation bore is 90 m from the discharging bore.
The logarithmic plot and match points are shown on Figure 6-15.
The left hand portion of the time-drawdown curve is superimposed on the family of Boulton type
curves and adjusted, while keeping the co-ordinate axes parallel, until a good matching position with
one of the left hand type-curve portions is found. This is the case with the left hand portion of the
curve for r/B = 0.6.
A match point (A) is chosen. This point is characterised by the following co-ordinate.
I/uA = 10
W(uA,r/B) = 1
s = 0.070 m
t = 16 min
= 1.11 x 10-2 days
Substitution of these values into equations 6.63 and 6.64 yields.
Now the right hand portion is superimposed on the right hand portion of type curve r/B = 0.6 and
again a match point (Z) is chosen. This match point is characterised by the following co-ordinate:
I/uY = 1
W(uY,r/B) = 1
s = 0.105 m
t = 250 min = 1.74 x 10-1 days
(a b log t )Q CQ n .....7.1
Generally n = 2 and the more common form of the equation to drawdown is:
(a b log t )Q CQ 2 .....7.2
where:
swt = drawdown in the discharging bore at time "t" after discharge commenced.
2.3
b
4 T = a constant in most cases. .....7.3
2.25T
a b log
rw2 S = a constant. .....7.4
C = a constant.
T = transmissivity.
Q = discharge rate.
rw = effective radius of the bore.
t = time since discharge began at discharge rate Q.
S = storage coefficient.
n = a constant (normally 2).
CQn = non-linear head loss at discharge rate Q.
The value of "b" and hence transmissivity, can be determined from a single test.
To determine "a" and "C" and fully describe the equation to drawdown in a discharging bore it is
necessary to solve at least two, and preferably more, simultaneous equations. This can be done
arithmetically or graphically.
This then requires two or more tests which are carried out at different discharge rates, or a variable
discharge test such as a step drawdown test.
Procedure
1. On transparent logarithmic paper, to the same scale as the type curve, plot Q/s w
against t/rw2 (or values of Q for a single well may be plotted against values of t).
2. Fit the plotted curve over the type curve.
3. Select a match-point and read off values of Q/sw, t/r2, G(α) and α.
4. Calculate T, from:
Q
T
2 G( )s w
2.25T
F
rw2 S = constant.
Equation 7.12 is in the form of an equation to a straight line with 1/Q as the ordinate and log t as the
abscissa. A plot of 1/Q against log t will yield a straight line with slope E. this slope is the change in
1/Q per log cycle and is called Δ1/Q.
Transmissivity
Rearranging equation 7.12 and substituting Δ l/Q for E we obtain:
2.3
T
4 s w (1 / Q) .....7.13
where:
sw = constant drawdown in the discharging bore.
Δ(1/Q) = change in 1/Q per log cycle of time.
0.000270 Hence:
T = 482 m2/day
1/Q (day/m3)
0.000275
0.000280
0.000285 0.000286
0.000290
0.000295
1 10 Time (mins) 100 1000
1 rw2 S 2
10
Tt .....7.19
sw 2 .3 2.25Tt
log 2
Q 4 T rw S
2.3
( s w / Q)
4 T
sw / Q 2.25Tt
log
( s w / Q) rw2 S .....7.20
2.25Tt
S
r 10(1 / Q ) / (1 / Q )
2
w
10 2 Tt
S
rw2
2.25
i.e. < 10-2
(1 / Q ) / (1 / Q )
10
10.00
15.00
Flow recession (Constant Draw dow n) segment of test. Each of these points w ith a constant draw dow n and
3883 m3/day
varying flow rates could have a line draw n through it parallel to the line through recovery line and be used to
20.00 determine the equation to draw dow n for the bore.
3731
3785
3709
3705
3633
3600
3458
3818
3623
3589
3523
3502
3480
3458
3546
3513
3437
3427
25.00
30.00
0.00 1000.00 2000.00 3000.00 4000.00 5000.00 6000.00 7000.00 8000.00 9000.00
i n
i 1 Qi log( t t i )
s wt
(a b log t ) CQ
Q .....7.25
which is an equation to a straight line.
If swt/Q is plotted against Q, then a straight line will result with slope C and intercept of (a + b log t)
on the swt/Q axis.
The value of "b" can be determined from the plot of drawdown versus log time on semi-logarithmic
paper.
Knowing the time "t" for which the various values of s wt were determined, "a" can be evaluated.
If the values of swt/Q were determined at t = 1, then the intercept on the swt/Q axis would be equal to
"a".
Case 2: n ≠ 2
When n ≠ 2, a plot of swt/Q versus Q will result in a curve which intersects the s wt/Q axis at
(a + b log t).
s wt
(a b log t ) CQ n 1
Q .....7.26
The value of the intercept, (a + b log t), is then subtracted from each value of s wt/Q to give:
swt - (a + b log t) = CQn-1 .....7.27
Taking log of both sides:
log (swt - (a + b log t)) = log C + (n-1) log Q
A plot of (swt/Q - (a + b log t)) versus Q on logarithmic graph paper will yield a straight line with slope
(n-1) and intercept C on the (swt/Q - (a + b log t)) axis.
The values of a, b, C and n are then readily evaluated.
Although an expression involving n ≠ 2 may occasionally appear to give a better fit to field data, it is
doubtful whether this refinement is ever justified. It is suggested that departures from a straight line
of the plot s /Q against Q are far more likely to be due to departures from the laminar flow
component from aQ than to departures of the non-linear flow component from CQ2.
s2 (a b log t )Q2 CQ 22
By subtraction:
s1 - s2 = (a + b log t) (Q1 - Q2) + C (Q12 - Q22)
Dividing by (Q1 - Q2):
s1 s2
(a b log t ) C (Q1 Q2 )
Q1 Q2 ....7.28
This is an equation to a straight line.
s1 s2
A plot of versus (Q1 + Q2) will result in a straight line with slope C and intercept of
Q1 Q2 s1 s2
(a + b log t) on the axis.
Q1 Q2
This same process can be used to compare all drawdowns at a particular time "t" with drawdowns
resulting from all other tests which have been carried out on a particular bore.
It is not necessary to use drawdowns in this analysis.
If:
h0 = the potentiometric head prior to pumping (s.w.l. or static head).
h1 = the potentiometric head at time "t" after commencement of pumping at discharge rate
Q1.
h2 = the potentiometric head at time "t" after commencement of pumping at discharge rate
Q2.
then:
s1 = h0 - h1
s2 = h0 - h2
and
s1-s2=h2-h1 .....7.29
This means then that water levels or, in the case of artesian bores, back pressures, may be inserted in
the equation instead of draw-down values.
Remember, however, that the order is reversed. That is h2 - h1 is inserted in place of s1 - s2.
Stage p1-pn
Q p Qn-Q1 Qn+Q1 sn-s1/Qn-Q1
(n) = sn-s1
1
In an important analysis it may be desirable not only to compute all differences in head loss between
Stage 1 and other stages as in Table 7-3, but to carry out this procedure using each of the stages in
turn as a base of reference.
In each case data from column 6 is plotted as abscissa against data from column 7 as ordinate. Each
plotted point should be labelled as to its source.
Sometimes it will be found that all combinations involving a certain observation are in disagreement
with those for all other combinations and this may be adequate reason for rejection of the
incompatible observation as being probably in error.
This pressure differential analysis has the advantage of giving quite a number of points but they are
plotted at a relatively high value of discharge. Where possible the drawdown method and the pressure
differential method should be used in conjunction with one another. More points will be available for
plotting and more confidence can be placed in the analysis.
7.4.3 Range of Intercepts
In some cases when plotting swt/Q versus Q a scatter of points results and difficulty is experienced in
selecting the best straight line to fit the points.
It would be desirable to have one end of the required straight line fixed - if not exactly at least within
certain limits.
The intercept is given as:
Intercept = (a + b log t) .....7.30
"b" is able to be determined from the semi-logarithmic plot as Δs/Q, and the time "t" is known. The
only unknown in the intercept is "a".
A range of intercepts can then be obtained by determining a range of values for "a".
Example
Take the case of an artesian bore:
Assume for the selected artesian bore the storage coefficient S lies between 10-3 and 10-5.
All practical values of r for a 152 mm diameter bore should lie between 10 -2 and 1 metre.
Therefore all practical values of r2S should lie between 10-3 and 10-7.
Since:
2.25T
a b log
r 2S
then:
2.25T 2.25T
b log a b log
10 3 10 7 .....7.31
b log 2.25T 103 a b log 2.25T 107 .....7.32
and "a" must fall within these limits.
Step Δs Δs/Q
Q sw 1min sw 1min/Q
(n) (m) (m/m3/day)
Step
Q sw 1min sn-s1 Qn-Q1 Qn+Q1 sn-s1/Qn-Q1
(n)
1 1135 4.59
2 1640 7.25
Step
Q sw 1min sn-s1 Qn-Q1 Qn+Q1 sn-s1/Qn-Q1
(n)
3 2180 10.6
Step
Q sw 1min sn-s1 Qn-Q1 Qn+Q1 sn-s1/Qn-Q1
(n)
4 2725 14.41
14.00 Step 4 (corrected), Q = 2725 m3/day, Δs = 1.54m The distance D3, at time t = x mins, is to be
D3 transferred to time t + x mins and added to
3
D3 the test data extended line for Step 3. In this
14.41m case
16.00 t4 = 80 mins
x = 40mins
t3 = 60 mins
t3 + x = 100 mins
18.00 Step 5 (corrected), Q = 11640 m3/day, Δs =0.82m D3 is added at 100 mins
22.00
24.00
1 10 Time (mins) 100 1000
draw dow ns step 2 as base step 3 as base step 4 as base step 5 as base
Figure 7-4 Step drawdown test – graphical analysis, determination of “a” and “C”
or
i m
s wt aQm b log(t t i ) Qi CQmn
i 1
i m
.....7.35
Since a, b and C are constants, then a plot of swt versus b log(t ti ) Qi on natural scale graph
i 1
paper will give a series of parallel straight lines (one for value of Q i) with slope "b" and intercepts of
(aQi + CQin) on the swt axis.
If these intercepts are divided by Qi, then "a” and "C" can be determined in the same manner as was
presented in sections 7.4.1 and 7.4.2.
8. For each value of "t", columns 6, 9, 12 etc are added and the sum i.e. i 1 Qi log(t ti ) is
recorded in the last column (column 13 in this case since only three steps are being
considered).
9. Column 3 is plotted against the last column; i.e. pressure (head or drawdown) is plotted
i n
against Qi log(t ti ) to give a series of parallel straight lines, one line for each value
i 1
of Q (as recorded in column 2).
10. The slope of the parallel lines is calculated and recorded. This slope is the value "b' i.e.
Δs/Q.
11. The intercepts are recorded. These intercepts record the drawdown at one time unit
(i.e. 1 minute, 1 day or whichever time unit is used in column 1). If the pressure
differential method is to be used to analyse the results, the distance between the
straight lines would be recorded.
12. The intercepts divided by the corresponding values of Q are plotted, on natural scale
graph paper, against the values of Q (or the corresponding pressure differential method
can be used).
In Figure 7-5, back pressures at the bore head, i.e. static head -drawdown, are plotted as well as
drawdowns listed in Table 7-2. This is similar to plotting negative values for drawdowns and negative
slopes result. This method of analysis is convenient when i accurate
n
values of static head, or S.W.L.,
are not known. In this case the static test crossed the zero i 1 Qi log(t ti ) line at 26.91 m indicating that
the correct static head for the bore was 26.91m and not 26.57m which was recorded at the start of
the test in Table 7-2. The cooling effect in the column of water in the bore prior to testing was then
some 0.34m. The equation to drawdown has been derived on the basis that the drawdown in each
step was 26.91m minus the back pressure during the step.
It will be noticed that the equation to drawdown is slightly different from that obtained in the Step
Drawdown Test Analysis, Figure 7-3 and Figure 7-4, when the antecedent flow recession and static
tests were ignored.
Table 7-5 Format for Eden-Hazel spreadsheet analysis
12.3m
3785
3458
3437
3883
3731
3709
3705
3633
3623
3600
3589
3546
3523
3513
3502
3458
3427
3480
Steps 1 - 18 Flow recession segment of test
0.00
0.00 1000.00 2000.00 3000.00 4000.00 5000.00 6000.00 7000.00 8000.00 9000.00
i n
i 1 Qi log( t t i )
Figure 7-6 Step drawdown test – Eden-Hazel analysis, determination of “a” and “C”
0 26.91 0
Note: the drawdown intercept at 1 minute can be obtained as the distance between the zero discharge line and the relevant
discharge line.
Base Step 20
Base Step 21
21 1640
Base Step 22
22 2180
Base Step 23
23 2725
For demonstration purposes both drawdown and pressure differential methods were used. In this case
and in most cases it is sufficient to use only the drawdown method.
The values of "a", "b" and "C" have now been determined and the equation to drawdown is:
swt = (a + b log t) Q + C Q2
= (3.1 x 10-3 + 5.34 x 10-4 log t) Q + 8.1 x 10-7 Q2
where:
Q = discharge rate in m3/day.
t = time in minutes since discharge began.
Q (log t t )
Note: time is in minutes in this equation because the time units used to calculate i i
for plotting were in minutes. If "day" time units had been used in the calculations and plotting we
would have had 1 day as the unit time intercept and the "a" term would have changed accordingly. It
would have been increased by "b log1440". The constants "b" and "C" terms in the equation to
drawdown would be unaltered and the time units in the equation to drawdown would have to be
days.)
This equation can now be used to produce drawdown versus discharge curves for various periods of
continuous pumping.
By way of illustration two curves have been drawn; one for a period of 10 5minutes (approximately
2 months) and one for a period of 106 minutes (approximately 2 years). These are given in
Figure 7-7.
where:
drawdown swt is in metres.
discharge Q is in m3/day.
30 time t is in minutes.
Note: the equation has been derived for time in minutes. If
time in days is used in the equation then the value of "a" will
have to be adjusted accordingly.
20 The plot shows the drawdown discharge relationship for times
of 2 months and 2 years continuous discharge.
10
0
0 500 1000 1500 2000 2500 3000 3500 4000 4500 5000
Discharge Rate (m3/day)
Figure 7-7 Drawdown versus discharge curves for various times of discharge
t (n 1 p) t ( n 1)
s wt s rt s w,t p
t (1 p ) t 1 .....7.39
2
The non-linear head loss (CQ ) is considered to be identical for discharge and recharge conditions.
Then:
t (n 1 p) t ( n 1)
a 2 a
s w,t (n 1 p) (n 1)(CQ 2
s ) s log t s w ,t p (n 1)(CQ s s log t
b t (1 p ) b t 1
t (n 1 p) t ( n 1)
s log t s w ,t p s log t
t (1 p ) t 1
t (n 1 p) t ( n 1)
s{ log t log t} s w,t p
t (1 p ) t 1
(n 1 p)!
s(log log(n 1)!) s w,t p
p
(n 1 p)!
s log s w ,t p
(n 1)! p .....7.40
where, for convenience the following has been written:
(n-1+p)! = (n-1+p)(n-2+p)(n-3+p).....(1+p)p
then:
s w,t F s s w,t
(n 1 p) p
.....7.41
1.2
(n 1 p)! p = 0.25
" F " log
(n 1)! p
1
0.8
F
0.6
0.4
0.2
0
1 10 100 1000 10000
Number of pumping Cycles (n)
(n 1 p)! n = 350
2.5 " F " log
(n 1)! p
n = 200
2
n = 70
F n = 30
1.5
1 n = 10
0.5
0
0 0.2 0.4 0.6 0.8 1 1.2
p
sw
QL QT
s105
.....7.42
1 2.07 2860
2 2.11
3 2.14
4 2.16
5 2.18 2860
10 2.22
20 2.28
40 2.33
60 2.37 2860
80 2.39
100 2.41
120 2.42
210 2.47
270 2.49
300 2.50
360 2.52
2.60
1 10 100 1000
Time (mins)
From Table 7-8 it can be seen that in this case the efficiency decreases with an increase in discharge
rate but increases with the duration of pumping. If continuous pumping is not required then in some
instances it may be more economical to pump at lower rates for longer periods.
7.8.2 When the Equation to Drawdown is Not Known
The above procedure is able to be used if the equation to drawdown for the bore is known. In many
instances the equation is not known and only a single pumping test is available. In such cases the
following procedure can be used.
The efficiency is determined using the same comparison as in the previous section i.e. the drawdown
in the bore operating at 100% efficiency is compared with the actual drawdowns in the pumping bore.
The theoretical drawdown at the bore for a 100% efficient bore will be:
2.3Q 2.25Tt
stw log 2
4 T rw S
and the actual drawdown in the pumping bore is taken from a plot of the test results. This analysis is
exactly the same as the previous case but it is limited to the test discharge rate. It does not give the
ability to determine the efficiency for rates other than the test rate.
Once again T can be calculated from the test data, rw is known and S can be assumed for confined
and unconfined conditions.
Driscoll uses more approximations still. Using the approximations to specific capacities given in
equation 7.45 and equation 7.46 for confined and unconfined aquifer conditions respectively, he
compares the specific capacity of a theoretical bore of 100% efficiency with the specific capacity
(Q/sw) of the pumping bore at time 1 day.
This is a very quick method and gives a good indication of the efficiency of the bore.
By way of comparison, the values of corrected drawdowns for Richmond Town Bore No. 3 as plotted
in Figure 7-3 were extrapolated to 1 day. The specific capacities were calculated for each of the
steps and the bore efficiencies estimated using Driscoll's approximation. The bore efficiencies for each
of the steps were also calculated using equation 7.50. Table 7-9 gives the efficiencies calculated
using each method.
Table 7-9 Comparison of efficiency calculations
Time Discharge Rate Efficiency (%)
(day) (m3/day) Using Using
Equation 7-50 Approximate Method
1 1135 76 72
1 1640 71 67
1 2180 66 61
1 2725 62 57
1 3270 58 54
1 0 3
y 4
a
x
a 2 S h0
h1 h2 h3 h4 4h0 .....8.4
T t
where:
Δh0 is the change in head at node (bore) 0 during the time interval Δt.
This is merely an introduction to numerical analysis.
Example
This data was presented originally in imperial units but has been converted for these notes.
Stallman tried this method successfully on several such arrays in the Arkansas Valley, Colorado, during
the winter of 1965-66. Bores 1- 4 were spaced 304m apart, so that "a" (304/√2) was equal to 216m.
From the estimated values of T and S, "a" normally is determined from the convenient empirical
relation a2S/T =10± days but in the Arkansas Valley, nearby boundaries made it necessary to use
a2S/T = 4± days. The elevations of the measuring points at each of the five wells were determined
above a convenient arbitrary datum, and the water levels, measured in metres above datum, were
obtained from automatic water level sensors. The winter data from the Nevius site near Lamar,
Colorado were plotted as shown in Figure 8-2, in which Σh, in mm, represents the left hand side of
equation 8.4, and Δh0/Δt is in mm/day.
From equation 8.4 note that the slope of the straight line, which Stallman drew form a least-squares
fit, is Σh/( Δh0/Δt) = 3.99 days.
Hence:
a2S/T = 3.99 days
The value T/S is known as the Hydraulic Diffusivity.
S was known from neutron moisture probe tests, made during periods of both high and low water
levels, to be about 0.18.
a2 = 4.7 x 104 m2
then:
T = (a2S/ Σh)/ (Δh0/Δt)
= (4.7 x 104 x 0.18) 3.99
= 2.12 x103 m2 / day
This value of T was in close agreement with the results of a nearby pumping test.
Conversely, if the value of T were known, the storage coefficient could have been determined.
Δl
Δw
h2
Q B CQ 2
(1 ) .....8.15
sw log t sw
Thus we see that Q/sw diminishes not only with time but with pumping rate Q. In unconfined aquifers
it may be necessary to adjust factor B to further account for delayed yield from storage.
The relation of specific capacity to discharge and time for a particular bore in a confined aquifer is
shown in Figure 8-5 (Jacob, 1947, Figure 5; converted to metric for these notes).
2.3Q r
T log 2 .....8.16
2 ( s1 s 2 ) r1
Let:
r1 = the effective radius of the bore, in this case 0.3 m;
r2 = radius of influence of bore, assume 300 m;
s1 = drawdown in the bore (one only drawdown, no test has been carried out); and
s2 = 0.
From this the following table was prepared, comparing average values of K for different alluvial
materials in the valley.
The same geologist who prepared Table 8-1 then carefully examined the logs of other bores and test
holes, for which no pumping tests were available. He assigned values of K to each bed of known
thickness, on the basis of the descriptive words used by the person who prepared the log. The values
of K assigned may have been equal to, or more or less than values given in the table (depending upon
cleanliness, sorting, mixing, etc), and thus necessarily involved subjective judgement; however, as
experience is gained, the geologist can estimate K and T with fair to good accuracy. The T values
from all sources also are compared carefully with the saturated-thickness map.
This method for estimating T has been used successfully in the Arkansas Valley in Colorado, in the
Arkansas Valley in Arkansas and Oklahoma (Bedinger and Emmett, 1963), in Nebraska, in California,
and can be used elsewhere.
8.3.5 Laboratory Analysis
Laboratory determinations for K of cores of consolidated rocks, such as partly to well cemented
sandstone, may be used in place of estimates. Reconstitution of disturbed samples of unconsolidated
material is not possible, however, so laboratory determinations for K generally do not give reliable
values. However, they may be very useful in indicating relative values. Table 1-2 at the end of
Section 1 gives indicative values of hydraulic conductivity for various rock types.
The above methods may also be used by the geologist in measuring exposed sections of rocks
containing water-bearing beds.
Sand, coarse 73
Clay 0.4
(Courtesy R. T. Hurr)
1 5 x 10-6 5 x 10-6
10 5 x 10-5 5 x 10-6
One may either multiply the thickness in metres by 5 x 10 -6 or interpolate between values in the first
two columns, thus for b = 300m, S = 1.5 x 10-3, etc. Values thus obtained are not absolutely correct;
they represent the minimum value of S as no adjustments have been made for porosity or for
compressibility of the aquifer. However, the values so obtained are reliable estimates.
8.4.2 Unconfined Aquifers
It is more difficult to estimate the specific yield of unconfined aquifers than the storage coefficient of
confined aquifers, but it can be done. It should be remembered that the specific yield applies only to
the material in the zone of water-table fluctuation and to the material within the cones of depression
of pumping bores. In general the specific yield ranges between 0.1 and 0.3 (10-30 percent), and long
periods of pumping may be required to drain water-bearing material.
where:
hx = the amplitude of groundwater fluctuations at distance x from shore.
h0 = the amplitude of tidal fluctuations at the shore.
x = distance of the observation bore from the shore.
t0 = the tidal period.
T = the transmissivity.
It has been further shown that the time lag t L of a given maximum or minimum water level in a bore
after high or low tide occurs can be obtained from:
tL x (t 0 S / 4 T ) .....8.23
If the storage coefficient has been estimated from the tidal efficiency then a value of transmissivity
could also be estimated from equation 8.21.
10.5 LEAKAGE
From the analysis of tests carried out on leaky aquifers the leakage coefficient, (1/c), can be
determined.
From this, the rate of leakage from an overlying aquifer to a lower aquifer can be computed by using
Darcy's Law.
Darcy's Law may be written as:
Q = -KiA
where:
K = the hydraulic conductivity.
i = the hydraulic gradient.
A = area through which the flow is occurring.
For vertical flow in the semi-confining layer or leaking layer Darcy's Law can be written as:
s
Q K| A .....10.5
b|
r12
4 *103
r22
whence:
r1
4 * 103 63.2
r2
Thus, changes in artesian head or pressure in a confined aquifer spread outward very readily from the
discharging bore, although it requires time for the drawdown to be measurable at a given distance
from a discharging bore. By contrast, in an unconfined aquifer, changes in water level occur very
slowly as gravity drainage takes place, so the cone of depression enlarges very slowly.
10.6.3 Determination of Radius of Influence
If it is assumed that we are dealing with an aquifer in porous media and that the assumptions for the
modified non-steady state equations apply then the drawdown at any point within the area of
influence at a given time t is given by:
2.3Q 2.25Tt
s log 2
4 T r S
The radius of influence at any time t is the radius at which the drawdown is zero.
For this to apply the log term must be zero i.e. the value within the log term must equal 1.
Thus at the radius of influence:
Tt
r0 1 .5 .....10.8
S
Bore1 Bore 2
A B
Bore 3 Bore 4
2.3Q 2.25Tt
sw (log10 2 )
4 s r S
The drawdown at any point in the aquifer is the summation of the individual drawdowns resulting
from each pumping bore.
The pumping bores must of course be drilled deep enough to enable such drawdowns to occur.
Table 10-1 is a spreadsheet showing the design of the borefield.
Design
Drawdown 15.00 15.00 15.00 15.00
Occasionally it is necessary to dewater a working area while still maintaining the groundwater levels
on one or more zones outside of the working area. One way of achieving this is to pump water from
the dewatering bores and then recharge the water back into the aquifer where required. This
technique is illustrated in Figure 10-1.
Figure 10-1 Recharging pumped water to maintain water levels in sensitive areas
10.7.1.4 Strip Pits
In strip pits (such as open cut mines) which intersect aquifers, water flows into the pit as parallel flow
through each of the four sides and as radial flow through the corners as shown in Figure 10-2. This
model is equivalent to the discharge from one bore, with a quarter of the bore at each corner of the
pit, and parallel flow through two sides and two ends. Equations 10.13 and 10.14 give the relevant
equations for this model for the unconfined and confined aquifer situations.
Strip Pit
Note: (1) The maximum drawdown in any aquifer is the distance from SWL to the bottom of that aquifer.
The drawdown in aquifer 2 will then be greater than that in aquifer 1.
(2) The radius of influence in any aquifer at any time is independent of discharge rate. It depends
only on T and S for that aquifer. The magnitude of drawdown within the cone of depression is dependent on
the discharge rate.
Ground Surface
Mine Discharge
SWL SWL
Aquifer 1
Aquifer 2
where:
QL = long term pumping rate.
QT = test rate.
sw = working drawdown.
s105 = test drawdown (extrapolated) to 105 minutes, or to the time required.
This rule should be applied with caution particularly in those cases where non-linear head losses are
present. It assumes, if non-linear head loss is present, that it is proportional to the discharge whereas
in fact the non-linear head loss is proportional to the discharge squared. In general it is considered
most unwise to predict a long term pumping rate beyond the constant discharge test rate unless you
are very familiar with the conditions in the area and the efficiency of the contractor.
The straight comparison of drawdowns to discharge is most useful when “ironing out” variations in
drawdown caused by relatively small variations in pumping rate during a pumping test.
The figure of 100,000 minutes is more or less an arbitrary end point which seems reasonably
convenient to use. The previous period is too short a time and the subsequent log period may be too
long a time for a normal irrigation cycle. A longer period is used of course for Town Water Supply or
industrial requirements.
Example
Table 11-1 Pumping test data
Time Drawdown Discharge Rate
(mins) (m) (m3/day)
1 2.07 2860
2 2.11
3 2.14
4 2.16
5 2.18 2860
10 2.22
20 2.28
40 2.33
60 2.37 2860
80 2.39
100 2.41
120 2.42
210 2.47
300 2.50
Table 11-1 presents test data obtained by the Queensland Irrigation and Water Supply Commission
in 1962.
Time (mins)
2.00
2.10
Drawdown (m)
2.20
2.30
2.40
2.50
2.60
1 10 100 1000
Figure 11-2.
At point A on the interface, for equilibrium to exist, the pressure on the salt water side of the interface
(pAs) must equal the pressure on the fresh water side (pAf).
From hydrostatics:
pAs = ρsgz .....11.2
and
pAf = ρfgz + ρfghf .....11.3
where:
pAs = pressure at A on the salt water side.
pAf = pressure at A on the fresh water side.
ρf = density of fresh water.
ρs = density of salt water.
hf = head of fresh water above mean sea level.
z = depth to the interface below mean sea level.
Since the pressure at A is the same on both sides of the interface:
ρsgz = ρfgz + ρfghf
or
f
z hf
( s f )
1000
= 0.5
1026 1000
= 19.2m below mean sea level
At this point the salt water interface would be expected to be 19.7m below the standing water level.
Remarks
Although the Ghyben-Herzberg concept, which is based on hydrostatics, implies no flow of the
interface, groundwater in coastal areas is invariably moving, and a dynamic concept is required.
Without fresh water flow i.e. no water level gradient, a horizontal interface would develop with the
fresh water floating on salt water.
A more correct picture of salt water intrusion is given in
Figure 11-3. Where the flow lines have a vertical component, the Ghyben-Herzberg concept gives too
small a depth to salt water. Further inland, where the flow lines are nearly horizontal, the error is
negligible.
The Ghyben-Herzberg concept is also applicable, with the same limitations, to confined aquifers,
where the potentiometric surface replaces the water table.
It should be remembered that sea water intrusion is quite natural and cannot be overcome
completely. What has to be determined is the magnitude of intrusion which will be acceptable in a
particular area. A fresh water flow to sea is necessary to stabilise the interface, and the position of the
wedge depends on the magnitude of the flow to the sea. A reduction in the magnitude of fresh water
flow will result in a movement of the toe of the wedge inland until stability is again achieved. A
smaller fresh water flow is associated with a smaller gradient which in turn results in a flatter wedge.
If the magnitude of the fresh water flow is increased, then the toe of the wedge moves seaward until
stability is achieved and a steeper wedge results. The operation of sea water intrusion is then a
management problem.
11.5.3 The Dynamic Concept
In many areas of sea water intrusion the depth of salt water computed from the Ghyben-Herzberg
relation differs markedly from observations. Often these discrepancies can be attributed to
assumptions that the head in the salt water is at mean sea level and that the fresh and salt waters are
static. An expression taking into account flows in each fluid can be derived by considering the fluid
heads.
Figure 11-3 hf and hs are potentiometric heads, related to a particular datum, in a region occupied by
fluids of densities ρf and ρs respectively. The potentiometric head is the summation of the elevation
head and the pressure head. As mean sea level has been adopted as datum level, then the following
relationships apply.
For fresh water, the potentiometric head, hf, is given by:
hf = z + p/( ρfg) .....11.5
and
hs = z + p/( ρsg) .....11.6
where:
p = pressure at A.
z = elevation at A above the datum level.
ρf = density of fresh water.
ρs = density of salt water.
hf = potentiometric head of fresh water.
hs = potentiometric head of salt water.
Since there is no flow across the interface the pressure (p) at A on the salt water side of the interface
must equal the pressure (p) at A on the freshwater side of the interface.
Equations 11.5 and 11.6 may then, at point A, be reduced to:
This then enables the depth to salt water, above datum, to be defined by the heads and densities
across the interface.
If hs = 0 then equation 11.7 reduces to equation 11.4, the Ghyben-Herzberg relation.
If hs = 0, i.e. there is no potentiometric head in the salt water then no movement of salt water occurs
and the wedge position is table.
However, if hs is negative, i.e. below mean sea level, then a hydraulic gradient exists from the sea to
the land and the wedge will move inland until stability is achieved.
If hs is positive, i.e. above mean sea level, then a hydraulic gradient exists from the land to the sea
and the wedge will move towards the sea until stability is achieved.
Example
Measurements of potentiometric heads in observation bores in a coastal aquifer reveal that at a
particular location the fresh water potentiometric head is 2m above mean sea level and the salt water
potentiometric head is 1m below mean sea level. If at this location, the densities of fresh water and
salt water are 1000 kg/m3 and 1026 kg/m3 respectively, calculate the depth of the salt water
interface.
Solution
From equation 11.7:
s f
z hs hf
s f s f
1026 1000
( 1) (2)
(1026 1000) 1026 1000
= -39.5 - 76.9 m
= -117.4m above mean sea level
The interface is then located 119.4m below the fresh water potentiometric level.
If the salt water head had been neglected the depth to the interface would be given by:
1000
z 2
(1026 1000)
= -76.9m above mean sea level
This indicates the invalidity of the Ghyben-Herzberg relation for general applications in intrusion areas.
Tidal fluctuations can cause fluctuations in potentiometric heads in coastal aquifer. Since the above
expressions are related to mean sea level, care should be taken to use the mean potentiometric head
during a tide cycle when applying these equations.
11.5.4 Location of the Interface
Take the case of a static salt water wedge where the origin is taken at the toe of the wedge.
In this static case the potentiometric head of the salt water is zero.
In practice there must be an outlet, where the aquifer is in contact with the sea, through which fresh
water flows to the sea. This outlet is shown on the following figures but the error involved in taking
the total length of the wedge, as shown on the figures to be the length of the wedge from shoreline is
considered to be negligible.
Figure 11-4.
The fresh water discharge rate (q) through a unit width of the aquifer, required to maintain a stable
wedge is given by:
dhx
q K (zx a) .....11.8
dx
From equation 11.7 (with hs = 0 and new origin):
f
z hx
( s f )
mx zx a
f
hx a
( s f )
f
i.e. dm dhx
( s f )
( s f )
and dhx dmx
f
dmx
q Km x .....11.9
dx
Integrating between x = 0 and x = L, and mx = b and mx = 0
0
L K( s f ) 2
qx 0 m x
2 f b
K( s f )
qL b2 ....11.10
2 f
where:
q = discharge rate per unit width of aquifer.
L = distance from shoreline to toe of wedge.
K = hydraulic conductivity of aquifer.
b = thickness of confined aquifer.
a = distance from top of confined aquifer to mean sea level.
Since all terms on the right hand side of equation 11.10 are constants, it may be concluded that the
length of the wedge is dependent on the unit discharge rate. If q increases, L must decrease and the
wedge moves seaward; if q decreases, L must increase and the wedge moves inland.
It will be observed also from the above equations that the shape of the wedge is parabolic.
11.5.4.2 Unconfined Aquifer
The saltwater wedge in an unconfined aquifer is illustrated in Figure 11-5.
dhx
q K (zx hx ) ....11.11
dx
From equation 12.29 (with hs =0 and new origin):
f
zx hx
( s f )
f dhx
and q Kh x (1 )
( s f ) dx
f
i.e. qdx Kh x (1 )dhx
s f
Kh02 f
qL (1 ) ....11.12
2 s f
Kz 02 s f f
qL ( ) 2 (1 ) ....11.13
2 f s f
If h0 is very small compared with z0 then ho + z0 (i.e. the saturated thickness "b") can be substituted
for z0 in equation 11.13 to give:
Kb 2
....11.14
s
qL ( s f )
2 f
It can be seen from equation 11.14 and equation 11.10 that qL for an unconfined aquifer is s times
qL for the confined aquifer.
where:
q = discharge rate per unit width of aquifer.
L = length of wedge.
K = hydraulic conductivity of the aquifer.
T = transmissivity of the aquifer.
zo = depth to bottom of aquifer below mean sea level.
b = saturated thickness.
The shape of the wedge is again parabolic. Again, since all terms on the right hand side of
equation 11.13 are constants then L is dependent on q.
11.5.5 Structure of the Interface
Since the two liquids are miscible the interface is not an ideal flow line of zero thickness but a
transition zone in which the water density varies from that of sea water to that of fresh. Field
measurements of interfaces have revealed a mixing zone ranging from a metre or so to some hundred
metres.
Some of the factors which can affect the location of the zone include tidal fluctuations, pumping and
natural recharge and discharge of fresh water. These influences cause the interface to shift continually
toward a new equilibrium position. Each movement, however, causes dispersion to occur, and a
transition zone with a salinity gradient is established.
Dispersion depends on the co-efficient of dispersion of the aquifer and the distance traversed by the
groundwater. The thickness of the transition zone at any location depends upon the co-efficient of
dispersion, the unsteady fresh water flow field, the hydraulic conductivity and the tidal pattern. One
could expect a thinner transition zone where the tidal range is small and thicker zone where the tidal
range is large.
11.5.6 Control of Intrusion
Once sea water intrusion develops in a coastal aquifer, it is not easy to control. The slow rates of
groundwater flow, the density differences between fresh and sea waters, and the flushing required
usually mean that contamination, once established, may require years to remove under natural
conditions.
Several methods have been suggested to control intrusion. These include the reduction of pumping or
modification of pumping practices, artificial recharge to create a mound parallel to the coast or the
establishment of a pumping trough parallel to the coast.
However, as stated previously control of sea water intrusion is basically a management problem. The
allowable magnitude of intrusion must be decided upon and entered as a constraint on the system
when determining annual yield.
Bear in mind that 1 m of water applied to an aquifer will have an entirely different effect on the salt
water wedge than would a 1m rise in water level.
Remarks
It should be remembered that, while a certain fresh water flow is required to maintain a stable wedge
it is not always necessary that this flow be lost to the seas. Much of it could be intercepted by a series
of collector trenches parallel to the coast and utilised inland.
r 2S
u -
4Tt
U - uniformity coefficient
- bulk volume
v - velocity (m/sec) (m/day)
- a leakage function
vs - seepage velocity (m/sec) (m/day)
V - volume (m3)
W(u) - well function of u
y - a dimensionless factor
- distance (cm) (m)
- thickness of fluid in viscosity (m)
Y0 - Bessel function
z - a dimensionless factor
α (alpha) - compressibility of soil matrix
- reciprocal of Boulton Delay Index
- function in constant drawdown analysis
Β (beta) - compressibility of water (kPa)-1
γ (gamma) - specific weight of fluid (use not recommended)
Δ(delta) - incremental value
η (eta) - dynamic viscosity (poise) (decapoise)
θ (theta) - porosity (dimensionless)
ν (nu) - kinematic viscosity (stokes)
π (pi) - 3.1416
ρ (rho) - density (kg/m3)
σ (sigma) - surface tension (N/m)
Σ (sigma) - summation
τ (tau) - intensity of shear
А α alpha A a
В β beta B b
Г γ gamma G g
Δ δ delta D d
Е ε epsilon Ĕ ĕ
Ζ ζ zeta Z z
Η η eta Ē ē
Θ θ theta Th th
Ι ι iota I i
Κ κ kappa K k
Λ λ lambda L l
Μ μ mu M m
Ν ν nu N n
Ξ ξ xi X x
Ο ο omicron Ŏ ŏ
Π π pi P p
Ρ ρ rho R r
Σ σ sigma S s
Τ τ tau T t
Υ υ upsilon Y y
Φ φ phi Ph ph
Χ χ chi Ch ch
Ψ ψ psi Ps ps
Ω ω omega Ō ō