Engineering Hydrology
Engineering Hydrology
com
1
byjusexamprep.com
0 – 25 Mild
25 – 50 Moderate
>50 Severe
𝐶𝑣 2
𝑁=( )
∈
Cv = coefficient of variation.
∈ = allowable percentage error.
For a given number of rain gauge standard error ∈𝑠
𝐶𝑣
∈𝑠 =
√𝑛
iv. ESTIMATION OF MISSING RAINFALL DATA
Let N1, N2, N3, … and Nx be the normal precipitation values for station' 1 to m', and 'x' Normal
precipitation is an average rainfall value for a day.
Let P1, P2, P3, … and Px be the rainfall for station '1 to m', and 'PX' is the rainfall of station x.
Case 1: when N1, N2, … Nm differs from Nx by less than 10%, the Value of Px is given as
follows
𝑃1 + 𝑃2 + 𝑃3 + 𝑃3 +. . . +𝑃𝑚
𝑃𝑥 =
𝑚
Case 2: when one or more of N1, N2,…Nm differs from Nx by 10% or more, the Value of Px is
calculated by:
𝑁𝑥 𝑃1 𝑃2 𝑃𝑚
𝑃𝑥 = [ + +. . . . . . . + ]
𝑚 𝑁1 𝑁2 𝑁𝑚
2
byjusexamprep.com
v. AVERAGE PRECIPITATION/RAINFALL
a. ArithmeticMean/Average Method
This method is suitable if rainfall is uniformly distributed and the area is not very large.
𝑃1 + 𝑃2 + 𝑃3 +. . . . . +𝑃𝑛
𝑃𝑎𝑣𝑔 =
𝑛
This method does not give very good results and hence is not used very frequently. Any
station outside the area of consideration is not taken into account in this method.
b. Thiessen polygon/Weighted Area Method
In this method, the rainfall recorded at each station is given a weightage based on an area
closest to the station. That is why this method is also known as the weightage area method.
𝑃1 𝐴1 + 𝑃2 𝐴2 +. . . +𝑃𝑛 𝐴𝑛
𝑃𝑎𝑣𝑔 =
𝐴1 + 𝐴3 +. . . +𝐴𝑛
∑𝑛
𝑖=1 𝑃𝑖 𝐴𝑖
⇒𝑃𝑎𝑣𝑔 =
𝐴
The ratio Ai/A is called the weightage factor for each station.
This method of finding average rainfall is suitable when the area is large, and rainfall is non-
uniformly distributed. This method is superior to the arithmetical mean method.
c. Isohyetal Method
An isohyet is a line joining all the points having the same value of rainfall, and isohyetal maps
are the one which shows contours of equal rainfall magnitude.
In an isohyetal method, it is assumed that the precipitation in areas between isohyetal lines is
equal to the mean of the precipitation of at isohyetal lines.
Mathematically, various following cases are possible.
Case 1:
𝑃1 + 𝑃2 𝑃 + 𝑃2 𝑃 + 𝑃𝑛
( ) 𝐴1 + ( 2 ) 𝐴2 +. . . + ( 𝑛−1 ) 𝐴𝑛−1
𝑃𝑎𝑣𝑔 = 2 2 2
𝐴1 + 𝐴2 + 𝐴3 +. . . +𝐴𝑛−1
Case 2: If isohyets outside the considered area are given.
3
byjusexamprep.com
𝑃0 + 𝑃1 𝑃 + 𝑃2 𝑃 + 𝑃3 𝑃 + 𝑃𝑛 𝑃 + 𝑃𝑛+1
( ) 𝐴0 + ( 1 ) 𝐴1 + ( 2 ) 𝐴3 + ( 𝑛−1 ) 𝐴𝑛−1 + ( 𝑛 )𝐴
𝑃𝑎𝑣𝑔 = 2 2 2 2 2
𝐴0 + 𝐴1 + 𝐴2 . . . . . +𝐴𝑛−1 + 𝐴𝑛
Case 3:
𝑃1 + 𝑃2 𝑃 + 𝑃3 𝑃 + 𝑃𝑛
( ) 𝐴1 + ( 2 ) 𝐴2 +. . . + ( 𝑛−1 ) 𝐴𝑛−1 + 𝑃𝑛+1 𝐴𝑛
𝑃𝑎𝑣𝑔 = 2 2 2
𝐴1 + 𝐴2 + 𝐴3 +. . . +𝐴𝑛−1 + 𝐴𝑛
CHAPTER-2- EVAPORATION
i. Measurement of evaporation
a. Experimental method
𝐿𝑎𝑘𝑒 𝐸𝑣𝑎𝑝𝑜𝑟𝑎𝑡𝑖𝑜𝑛 = 𝐶𝑃 × 𝑃𝑎𝑛 𝑒𝑣𝑎𝑝𝑜𝑟𝑎𝑡𝑖𝑜𝑛
1. Class A 0.7
3. Colorado 0.78
b. Empirical Method
V
E = k m (e w −ea ) 1 + 9
16
E = Rate of evaporation per day.
Km = constant, which depends on the size of the water body.
ew = Saturated vapour pressure in mm of mercury.
ea = vapour pressure of air (ea) in mm of mercury.
V9 = Mean monthly wind velocity in km/hr at the height of about 9m from the ground
surface.
ii. Transpiration
➢ It is the process by which water leaves the body of a living plant and reaches the
atmosphere as water vapour.
➢ It can be measured by an instrument called as phytometer
➢ Following are the methods to find out the evapotranspiration.
• Lysimeter
• Field plot
• Penman's equation
• Blaney Criddle equation.
4
byjusexamprep.com
L
Mathematically, V =
t1
b. Current meter
V = aNs + b
To reduce the drag force, these are streamlined in shape. The value of this sounding weight is
given as
W = 50 v y
Here, w = weight in ‘N’
5
byjusexamprep.com
V0.2y + V0.8y
Vavg =
2
Case 2: For shallow k/a one point formula
Vavg = V0.6y
Vavg = kvs
n
Q= Qi
i=1
Q1 = A1V1
Q2 = A2 V2
Qn = An Vn
Here, A1 = w1y1
A2 = w2y2
An = wny2
Here,
2
w2
w1 + 2
w1 =
2w1
w3 + w4
w3 =
2
2
wn −1
wn +
2
Similarly, wn −1 =
2wn
6
byjusexamprep.com
b. Dilution method:
l 1 1
v= −
2 cos t1 t
Q = AV
d. Moving Boat Method:
0 + y1 y1
Q1 = A1Vf = Vb t1 vf = vb0vf t1
2 2
y1
= Vr2 sin cos t1
2
Similarly,
y + y2
Q2 = Vr2 sin t2 1
2
y + y3
Q3 = Vr2 sin t3 2
2
7
byjusexamprep.com
e. Slope Area Method :
V2 − V22
hf = (h1 − h2 ) + 1 − he
2g
hf = frictional
hc = eddies
Through experiments, it has been found that eddy head loss
V12 − V22
he = Ke
2g
hf
- Q=k
L
K = (K1 K2 K3 K 4 K5...... Kn )
1/n
CHAPTER-4- INFILTRATION
i. INFILTRATION CAPACITY
Ft = Ff + ( Fi − Ff )e( − kt )
Ft = Infiltration rate or capacity at a time 't'.
Ff = Final infiltration rate or capacity.
Fi = Initial infiltration rate or capacity.
K = Decay constant (T –1
or /s or / hr)
ii. INFILTRATION INDICES
a. 𝝓 - index
8
byjusexamprep.com
The 𝜙 - index is the average rainfall above which the rainfall volume is equal to runoff volume.
The 𝜙 𝐼𝑛𝑑𝑒𝑥 is derived from the rainfall hyetograph with the knowledge of the resulting runoff
volume. The initial losses are also considered infiltration
b. w-index
This is the average infiltration rate during the entire duration of rainfall. In calculatingthe w-
index, the initial losses are separated from total abstractions to refine the 𝜙 𝐼𝑛𝑑𝑒𝑥.
Mathematically it is defined as
P − R − ta
w − index =
Tc
9
byjusexamprep.com
Here, n is in cm
A is in km2
t is in hr
O is in m3/sec
ii. SYNTHETIC UNIT HYDROGRAPH
Snyder selected three parameters for the development of SUH. These parameters relate to the
catchment characteristics
(i) Basin time width T
(ii) Peak discharge QP
(iii) Lag time,i.e. basin lag time tp. (Snyder defined lag time as the time interval from mid pt.
of rainfall to the peak of UH (instead of centroid).
10
byjusexamprep.com
Where tp is in hr
Ct is a coefficient reflecting slope, land use, and associated storage characteristics of the basin.
Its value varies between 1.35 to 1.65, the average being 1.5
L = basin length measured along the watercourse from the basin divide to the gauging station
in km.
Lca = Distance of centroid of the catchment from the gauging point (in km)
T is in hr
Qp is in m3/s
A = Catchment area in km2
11
byjusexamprep.com
D '− tr
tpr = tp =
4
To plot the smooth synthetic unit hydrograph, US army crops of engineering gave the width of
SUH as
5.87
W50 = 1.08
Qp
A
3.35 W50
W75 = =
Qp
1.08
1.75
A
W50 and W75 are the widths of synthetic unit hydrograph in hr at 50% and 75% of Qp,
respectively. Where Qp is in m3/s and A is an area of the catchment in km2.
12
byjusexamprep.com
i. STEADY CONFINED FLOW (FULLY PENETRATING WELL)
A. Theims Theory
2KB (h2 − h1 )
Q=
r2
ln
r1
h1 + s1 = h2 + s2
Q= Rate of flow
r2, r1 = radius
B. Duipit’s theory
Further, at the edge of the zone of influence, s =0, r 2 = R and h2= H;at the well wall r1 = rw,
2KBSw
Q=
R
ln
rw
13
byjusexamprep.com
ii. STEADY UNCONFINED FLOW
A. Theim’s Theory
Q=
(
k h22 − h12 )
r
ln 2
r1
B. Duipit’s theory
Q=
(
k H2 − h2w )
R
ln
rw
R = 3000 S K
hw = depth of water in the pumping well of radius rw.
14
byjusexamprep.com
ii. EMPIRICAL FORMULAE
The empirical formulae used for estimating the flood peak are essentially regional formulae
based on statistical correlation of the observed peak and important catchment properties.
These equations are given as follows.
(a) Dicken's Equation
This equation is applicable for North and central India. Mathematically, it is given as
QP = CD A3/4
Here,
QP = m3/s
A = Area in km2
CD =dicken's constant (6 to 30)
(b) Ryve's Equation
This equation is applicable for southern parts of India. Mathematically it is given as
QP = CR A2/3
124A
QP = 123 A
A + 10.4
p = probability of exceedance.
15
byjusexamprep.com
Reliability
This is the probability that a particular flood or rainfall is never equalled or exceeded (r = 0) in
a period of 'n' years, Mathematically
⇒ Re liability = (1 − p)n
Risk
This is a probability that a particular flood or rainfall is equalled or exceeded at least once in a
period of 'n' years.
P (at least once) = 1 – p(never happen)= 1- Reliability
The most commonly used distributions are :
(a) Gumbel's distribution
(b) Log Pearson Type III distribution.
Gumbel's Method:
As per Gumbel's method
X T = X + K . n −1
X
X = Mean Value for variate = X (from annual series)
n
n = no. of yrs of record
σn-1 = Standard deviation of the sample of size n
(X − X)
2
n−1 =
n −1
K = Frequency factor
y T − yn
K=
Sn
yT = reduced variate
1 T
= (−) ln ln
t T − 1
16
byjusexamprep.com
iv. CONFIDENCE LIMIT
Confidence interval indicates the limits about the calculated value between which the true value
can be said to lie with a specific probability based on sampling errors only.
For confidence probability ' ', the confidence interval of variate XT is bounded by values X1
and X2 given by
X1 2 = XT f ( ) .Sc
where f( ) = function of confidence probability ‘α’. It can be found using the following table
Α is percent 50 68 80 90 95 99
f( ) 0.674 1.0 1.282 1.645 1.96 2.58
N −1
Sc = Probable error = b
n
b = 1 + 1.3K + 1.1 K2
y T − yn
K=
Sn
T
y T = − ln ln
T −1
n = Sample size
T = Return period
σn–1 = Standard deviation of sample
v. Rainfall run off correlation
R = ap + b
Here,
n PR − P R
a=
n P2 − ( P)2
b=
R − a P
n
vi. Flood Routing
17
byjusexamprep.com
b. Wedge storage
It is the wedge-like volume formed between the actual surface profile and the top surface of
prism storage.
c. Muskingum method
S = (I − Q)t
I + I1 Q2 + Q1
S = 2 t − t
2 2
S1 = K XI1 + (1 – x ) Q1
S2 = K XI2 + (1 – x ) Q2
****
18