Engineering Hydrology Lecture Notes 2010
Engineering Hydrology Lecture Notes 2010
LECTURE NOTES
CHAPTER ONE
1 METEOROLOGICAL AND HYDROLOGICAL DATA
ANALYSIS
1
To represent the actual processes taking place on the
environment and better estimation of quantities out of it
hydrological studies require extensive analysis of
meteorological
hydrological
Errors due to splashing and evaporation usually are small and can
be neglected.
4
Systematic errors due to the mechanics of operation of the
instrument can be minimized by installing a non recording
type gauge adjacent to each recording gauge to assure that at
least the total precipitation is measured.
5
b) Estimation of Missing Data
Several approaches are used to estimate the missing values.
where annual precipitation value b/n neighboring gauges differs
by less than 10% use
Station Average method,
Px=1/m(P1+P2+…Pm)
Px = The missing precipitation record
P1, P2 , …, Pm = Precipitation records at the neighboring stations
M = Number of neighboring stations
6
where annual precipitation value b/n neighboring gauges
differs by greater than 10% the ff methods will use.
Normal Ratio method
Px=Nx/M(P1/N1+P2/N2+…+Pm/Nm)
Nx = Annual-average precipitation at the gage with missing
values
N1 , N2 , …, Nm = Annual average precipitation at
neighboring gauges
Inverse Distance method weights the annual
average values only by their distances, dm, (no
need of information about average annual
precipitation at the gauges)
7
D=∑dm^-b
Px=1/D(∑dm^-b Nx)
10
1.3 Areal Estimation
Several approaches have been devised for
estimating areal precipitation from point
measurements.
1. The Arithmetic mean:- uses the mean of
precipitation record from all gauges in a
catchment.
2. Thiessen polygon(introduces a weighting factor
on rational basis).
3.Isohyetal method(The method requires the
preparation of the isohytal map of the catchment
11
from a network of gauging stations).
1.4 Hydrological Data
12
Chapter Two
Rainfall-Runoff Relationships
2.1 Hydrological Models
The two classical types of hydrological models are
1. The deterministic
2. And stochastic model
13
14
2.2 Deterministic Hydrological Models
Deterministic models permit only one outcome from a
simulation with one set of inputs and parameter values.
The first of this kind of model was the Rational Method published
by16 the Irish engineer Thomas James
The model was a single simple equation often used for drainage
design for small suburban and urban watersheds.
18
2.2.3 Distributed Process Description Based Models
20
2.4 Rational Method
One of the most commonly used for the calculation of peak
flow from small areas is the rational formula given as:
21
Assumptions inherent in the Rational Formula are as
follows:
The peak flow occurs when the entire watershed is
contributing to the flow
The rainfall intensity is the same over the entire drainage
area
The rainfall intensity is uniform over a time duration
equal to the time of concentration, tc . the time of
concentration is the time required for water to travel from
the hydraulically most remote point of the basin to the point
of interest
The frequency of the computed peak flow is the same as
that of the rainfall intensity, i.e., the 10-yr rainfall intensity
is assumed to produce the 10-yr peak flow
The coefficient of runoff is the same for all storms of all
recurrence probabilities
22
2.4.1 Runoff Coefficient
The ground cover and a host of other hydrologic abstractions
considerably affect the coefficient.
The Values of C vary from 0.05 for flat sandy areas to 0.95 for
impervious urban surfaces,
If the basin contains varying amount of different land cover or
other abstractions, a coefficient can be calculated through areal
weighing as shown in following equation
23
24
2.4.2 Rainfall intensity
Rainfall intensity,
duration and frequency
curves are necessary to
use the Rational method.
25
2.4.3 Time of Concentration
There are a number of methods that can be used to estimate
time of concentration (tc), some of which are intended to
calculate the flow velocity within individual segments of
the flow path (e.g. shallow concentrated flow, open channel
flow, etc.)
Sheet flow travel time
26
Shallow concentrated flow velocity
27
For a circular pipe flowing full, the hydraulic radius is one-
28
For small natural catchments,
29
2.5. SCS Curve Number Method
30
31
32
2.6. Time-Area Method
The time – area method of obtaining runoff or discharge from rainfall can
be considered as an extension and improvement of the rational method. The
peak discharge Qp is the sum of flow – contributions from subdivisions of
the catchment defined by time contours (called isochrones), which are lines
of equal flow – time to the river section where Q p is required.
As the assumption for the rational method, the whole
catchment is taken to be contributing to the flow after T
equals to Tc.
Hence the peak flow contributed from the whole catchment
after Tc of the commencement of rain is:
Q
p =
∑ i(n−k ) ∆A(k )
k -1
Where n, the number of incremental areas between
successive isochrones, is given by Tc/∆T, and k is a counter.
33
The unrealistic assumption made in the rational method of uniform
rainfall intensity over the whole catchment and during the whole of T c
is avoided in the time – area method, where the catchment
contributions are subdivided in time. The varying intensities within a
storm are averaged over discrete periods according to the isochrones
time interval selected. Hence, in deriving a flood peak for design
purposes, a design storm with a critical sequence of intensities can be
used for the maximum intensities applied to the contributing areas of
the catchment that have most rapid runoff. However, when such
differences within a catchment are considered, there arises difficulty
in determining Tc, the time after the commencement of the storm
34
when, by definition, Qp occurs.
2.7 STREAM FLOW HYDROGRAPH
A hydrograph is a graphical plot of discharge of a natural stream
or river versus time.
The hydrograph is a result of…….
By definition, the volume of water under an effective rainfall
hyetograph is equal to the volume of surface runoff.
Characteristic parts of hydrograph
the rising limb,
the crest segment
and the falling limb or depletion curve.
35
With reference to the following figure the effective rainfall
hyetograph consisting of a single block of rainfall with duration
D (T is also used in the lecture note alternatively)
36
2.7.1 Hydrograph Analysis
One of the major tasks of the hydrograph analysis is to produce
rainfall-runoff relationships for a catchment area,
37
Depending upon the unit of time involved, we
have:
Annual hydrographs showing the variation of daily or
weekly or 10 days daily mean flows over a year.
Monthly hydrograph
Seasonal hydrographs representing the variation of the
discharge in a particular season such as the monsoon
season or dry season
Flood hydrographs or hydrographs due to storm
representing stream flow due to a storm over a catchment
38
Two Basic component of hydrograph
Surface flow
Bas flow
39
The length of the delay before the river rises depends on the
wetness of the catchment before the storm and on the intensity of
the rainfall itself.
The volume of surface runoff, represented by the area under the
hydrograph minus the base flow,
The boundary between surface runoff and base flow is difficult
to define and depends very much on the
geological structure
and composition of the catchment.
Permeable aquifers, such as limestone and
sandstone strata,
but impervious clays and built-up areas
The base flow levels are also affected by the general climatic
state
40 of the area:
The main aims of the engineering hydrologist
are to quantify the various components of the
hydrograph, by analyzing past events, in order
to relate effective rainfall to surface runoff,
and thereby to be able to estimate and design
for future events.
41
2.7.2 Factors affecting flood hydrograph
42
2.7.3 Effective Rainfall
The portion of rainfall that finds its way into a river is
known as the effective rainfall, the rest being lost in
evaporation, detention on the vegetation and ground
surface or retention in the soil.
43
For the purposes of correlating direct runoff hydrograph
(DRH) with the rainfall, which produces the flow, it is
necessary to obtain the effective rainfall hydrograph
(hyetograph) (ERH) which can be obtained by deducting
the losses from the total rain.
At the beginning of a storm there could be considerable
interception of the rainfall and initial wetting of surfaces
before the rainfall become ‘effective’ to form surface
runoff.
The loss-rate is dependent on the state of the catchment
before the storm. But it is difficult to assess
44
quantitatively.
The two simplified methods of determining the effective
rainfall are:
I. The φ-index method
II. The initial and continuing loss method.
I. The φ-index method: this method assumes a constant
loss rate of φ-mm from the beginning of the rainfall
event. This amount accounts for interception,
evaporation loss and surface detention in pools and
hollows.
II. Initial and continuing loss rate method: In this method
all the rainfall up to the time of rise of the hydrograph
is considered lost, and there is a continuing loss-rate at
45
same level after words.
The total runoff consists of direct runoff and the base
flow. For hydrograph analysis the base flow has to be
separated from the total runoff.
There are several methods of base flow separation. Some
of them that are in common use are:
Straight-line method
Method-I, Method-II, Method-III
49
The definition of unit hydrograph implies the following.
51
The total hydrograph of direct runoff due to n successive amounts of
effective rainfall (for instance R1 and R2) is equal to the sum of the
n successive hydrographs produced by the effective rainfall (the
latter lagged by T h on the former). This is known as Law of
Superposition.
The third property of TUH assumes that the effective rainfall-
surface runoff relationship does not change with time, i.e., the same
TUH always occurs whenever the unit of effective rainfall in T h is
applied on the catchment. Using this time invariance assumption,
once a TUH has been derived for a catchment area, it could be used
to represent the response of the catchment whenever required.
52
The derivation of the unit hydrograph of a catchment from
single storms proceeds in the following stages:
54
4. The ordinates of the surface runoff ydrograph are divided by
the runoff depth d due to the ordinates of the unit hydrograph.
55
Changing of the Duration of the UH
By superposition
58
This S-Curve is due to a D-h unit hydrograph. It
has an initial steep portion and reaches a
maximum equilibrium discharge at a time equal to
the time base of the 1st unit hydrograph. This
magnitude can be calculated by
Where,
Qs = the maximum rate at which an ER
intensity of 1/T can drain out of the catchment
of area, A (km2)
T1 = unit storm in hours
59
Consider two D-h S-curves A and B displaced
by T-h if the ordinates of B are subtracted from
that of A, the resulting curve is a DRH
produced by a rainfall excess of duration T-h
and magnitude [1/D]*T cm , hence if the ordinate
difference of A and B ie, (SA-SB) are divide by T/D, the
resulting ordinates a hydrograph due to an ER of 1cm
and T-h ie, T-h unit hydrograph
60
2.9 Applications of Unit Hydrograph
As the UH establishes a relationship between the DRH and
ERH for a catchment, they are of immense value in the
study of the hydrology of a catchment.
They are of great use in:
The development of flood hydrograph for extreme rainfall
magnitudes for use in design of hydraulic structures
Extension of flood-flow records based on rainfall records
The development of flood forecasting and warning systems
61
2.10 Synthetic Unit Hydrographs
To develop unit hydrographs to a catchment, detailed
information about the rainfall and the resulting flood
hydrograph are needed. However, such information might
be available only at a few locations and in a majority of
catchments, especially those, which are at remote
locations; the data could normally be scarce. In order to
construct UH for such areas, empirical equations of
regional validity, which relate the important hydrograph
characteristics to the basin characteristics are of most
important. Unit hydrographs derived from such
relationships are known as synthetic unit hydrographs.
62
2.10.1 Snyder’s method
The most important characteristics of a basin affecting a
hydrograph due to a given storm is basin lag. Actually basin
lag (also known as lag time). Its value is determined essentially
on the physical features of the catchment, such as size, length,
stream density and vegetation.
For simplicity, Snyder has used a somewhat different definition
of basin lag (denoted by tp) in his methodology. This tp is
practically of the same order of magnitude as TL and in this
section the term basin lag is used to denote Snyder’s tp.
The first of the Snyder’s equation relates the basin lag tp.
Defined as the time interval from the mid point of the unit
rainfall excess to the peak of the unit hydrograph (Figure
below) to the basin characteristics as:
63 tp = Ct (LLc)0.3
Where,
tp in hours
L = basin length measured along the watercourse from the basin divide to the
gauging station in km.
Lc = distance along the main watercourse from the gauging station to the point
opposite (or nearest) the watershed centroid in km
Ct = a regional constant representing watershed slope and storage
The value of Ct in Snyder’s study ranged from 1.35 to 1.65. However, studies by
many investigators have shown that Ct depends upon the region under study and
wide variations with the value of Ct ranging from 0.3 to 6.0 have been reported.
64
Elements of synthetic unit hydrograph
65
CtL and n are basin constants. (n= 0.38 and CtL = 1.715, 1.03, 0.50
for mountainous, foot-hill and valley drainages of USA) Standard
duration of effective rainfall, tr (in hours)
66
Where t’p = basin lag in hours for an effective duration of tR.
Therefore Qp,
Note that when tR = tr implies QP = Qps The time base of unit
hydrograph is given by Snyder as:
67
Where W50 = width of unit hydrograph in hour at 50% peak
discharge W75 = width of unit hydrograph in hour at 75% peak
discharge q = Qp/A = peak discharge per unit catchment area in
m3/s/km2 Since the coefficients Ct and Cp vary from region to
region, in practical applications it is advisable that the value of these
coefficients are determined from known unit hydrographs of
meteorologically homogeneous catchments and then used in the basin
under study. This way Snyder’s equations are of use in scaling the
hydrograph information from one catchment to another similar
c
68
catchment.
2.11 UH from a complex storm
In nature storms are most likely occurring with changing
intensities over their total duration.
69
At various time intervals 1D, 2D, 3D, … from the start of the
ERH, let the ordinates of the unit hydrograph be u1, u2, u3, …
and the ordinates of the composite DRH be Q1, Q2, Q3,….
Then
Q1 = R1u1
Q2=R1u2 + R2u1
Q3= R1u3 +R2u2 + R3u1
Q4= R1u4 + R2u3 +R3u2
Q5= R1u5 + R2u4 +R3u3
From the above equation U1, U2, U3… can be determined
70
• The first of the Snyder’s equation relates the basin lag tp.
Defined as the time interval from the mid point of the unit
rainfall excess to the peak of the unit hydrograph
tp = Ct (LLc)0.3
• Where,
• tp in hours
• L = basin length measured along the
watercourse from the basin divide to the
gauging station in km.
• Lc = distance along the main watercourse
from the gauging station to the point opposite
(or nearest) the watershed centroid in km
• Ct = a regional constant representing watershed
slope and storage
71
72
→Important relationships:
Basin lag tp
CtL and n are basin constants. (n= 0.38 and CtL = 1.715,
1.03, 0.50 for mountainous, foot-hill and valley
drainages respectively
Standard duration of effective rainfall, tr (in
hours)
Where
W50 = width of unit hydrograph in hour at 50%
peak discharge
W75 = width of unit hydrograph in hour at
75% peak discharge
75
q = Qp/A = peak discharge per unit catchment
3 2
2.12 Instantaneous unit Hydrograph (IUH)
76
The limiting case of a unit hydrograph of zero
duration is known as instantaneous unit
hydrograph (IUH).
3. u(t) → = 0 for t → ∞;
77
2.13 Dimensionless Unit Hydrograph
78
2.14 Hydrology of Ungauged Catchments
81
Utilizing the records of average annual
precipitation input to the basins at measured
streams nearby or having similar hydrologic
characteristics, a runoff coefficient is estimated for
the drainage basin being studied.
The product of this coefficient and the computed
normal annual precipitation input to the basin and
the basin area can be used to calculate the average
annual discharge as:
82