Ch. 8 - Flood Routing
Ch. 8 - Flood Routing
INTRODUCTION
The knowledge of Hydrology is essential to ensure proper functioning of hydraulic
structures like Barrages and Dams. Storm water from catchment ultimately reaches
streams. At a certain time, a peak flow is built up. For safety of a Barrage or Dam it is
necessary to plan its safe-passage over the structure. This requires flood forecasting.
Flood routing is one of techniques which may be used for flood forecasting.
Basing on the type of equations used, different techniques have been developed to route
the flood. These are Hydraulic Routing and Hydrologic Routing. Hydraulic Routing is
based on equations of continuity and momentum whereas Hydrologic Routing is based on
the equation of continuity and empirical equation.
These two techniques are used for catchment routing, reservoir routing and channel
or river routing.
RESERVOIR ROUTING
The main components of Dam are an inflow channel, storage reservoir and outflow
structures like spillways, tunnels etc. Once a flood enters a reservoir, part of it may be
stored in the reservoir and balance safely passes through or over outflow structures. The
main function of a reservoir is to store water, from which releases are made according to
water-demands on downstream of reservoir. A multipurpose hydroelectric project has
storage of water as well as generation of electricity.
The reservoirs may be small or large. An example of small reservoir is pond of a barrage.
A small reservoir has small capacity and hence water levels in barrage pond are sensitive
to quantity of inflow and time of inflow. The outflow from a small reservoir is solely a
function of pond elevation if the outflow is not controlled. In case of large reservoirs the
moderate inflow may not have large impact on reservoir elevation however large floods
need to be negotiated keeping in view the operational rules. Sometimes the water
requirement may have marked influence on releases from the reservoir. The outflow from
the reservoir also includes discharge through power house. The inflow hydrograph is
tackled in such a way that peak outflow is within safe limits. The inflow hydrograph may
be a forecast flood or an actual one.
I-O=ds/dt - - - - - - - - - - - - - - - - - - - - - - - - - - - (8.1)
Equation (8.1) shows that if inflow is assumed constant the reservoir storage is a mere
function of outflow. This condition is however a simplification and may not be met
practically. In case of large reservoirs many other factors contribute a substantial change
in storage e.g. groundwater flow. If average values of inflow and outflow are considered
for opted time interval ∆t then equation (8.1) can be written as
As mentioned above the subscripts ‘j’ and ‘j+1’ denote values of Inflow, Outflow and
Storage at beginning and end of jth time interval ∆t. The time ∆t is known as routing
period. This period should not be so large that peak of inflow hydrograph is not
intercepted.
208
Elevation vs Surface Area Realationship
124
122
120
118
Elevation (m)
116
114
112
110
108
106
104
102
40000 42000 44000 46000 48000 50000 52000 54000
Surface Area (m²)
2. The discharge capacity of overflow structure with change in water level should be
calculated. For this purpose the applicable discharge formula need to be applied.
The well known weir equation is:
Q = Cd B H3/2 - - - - - - - - - - - - - - - - - - - - - - - - - (8.3A)
Where,
Q = Total Discharge
Cd = Coefficient of Discharge
B = Width of Weir
H = Differential head over the crest of the weir neglecting velocity of approach
The coefficient of discharge depends on degree of submergence of the weir. Its value is
determined experimentally e.g. by model tests. The value of discharge can also be
determined from Gibson’s curve. Its value generally ranges from 1.6 to 2.2 having units
of m1/2/s. A mean value of 1.70 is often used in SI units.
Q = Cd A V - - - - - - - - - - - - - - - - - - - - - - - - - (8.3B)
Where,
Q = Discharge in m³/s
A = Area of Flow
V = velocity of flow
If ‘H’ is the head of water up to centre of pipe, Equation (8.3B) can be written as
209
Q = Cd A H1.5 √2g
‘A’ is cross-sectional area of the pipe, ‘H’ is water head and ‘g’ is acceleration due to
gravity. Its value is 9.81 m/sec² in SI units. The value of coefficient of discharge is
determined by model tests. If velocity of approach is appreciable the head ‘H’ must be
increased by velocity head. Equation (8.4) then becomes
For other types of outflow structures like Sluice Gates, etc. different equations are used
for calculations of discharge and can be found in books of hydraulics. Once the outflow is
determined for different reservoir elevations, a graph is plotted between storage and
outflow. A typical such graph is shown in figure (8.2). Please note that outflow is taken
along x-axis and [(2S/∆t) + O] is taken along y-axis. Where ‘S’ is storage and ‘O’ is
outflow / discharge. The quantity [(2S/∆t)+O] is called ‘Storage Indication’.
300
250
200
Outflow (m³/s)
150
100
50
-
- 500 1,000 1,500 2,000 2,500 3,000 3,500 4,000
Figure 8.2 Outflow and Storage Indication Relationship for certain reservoir
3. The inflow hydrograph should be known. It may be actual or forecasted flood. The
inflow is added for successive values to get I1+I2. The initial outflow is assumed
equal to inflow. Corresponding to this initial outflow storage indication [(2S/∆t) + O]
is found from storage indication curve. To this value double of outflow is subtracted
to get [(2S/∆t) - O]. To this value of [(2S/∆t) - O], I1+I2 is added to get next value of
[(2S/∆t) + O]. Read out next outflow from storage indication curve and repeat the
procedure till whole of inflow hydrograph is used to get outflow values.
210
4. Now the inflow and outflow hydrographs are plotted. The difference in peak of
inflow and outflow hydrograph is known as attenuation and time between two peaks
is known as reservoir lag. Please note that by routing procedure we can find expected
maximum reservoir level due to certain inflow hydrograph. This is done to ensure that
reservoir level remains within safe limits.
Example 8.1
Table 8.1 shows Inflow hydrograph entering a small reservoir having a Ogee shaped
overflow spillway. The weir has a length of 3 meters and coefficient of discharge of 2.2.
The weir crest is at elevation 104.00. When the inflow hydrograph enters the reservoir, its
water level is at 140.00 m and Storage of 40000 m³. The elevation versus storage
relationship is given in Table 8.2. Route the inflow hydrograph assuming no base flow.
0 50 9 295 18 50
1 75 10 252 19 45
2 175 11 212 20 39
3 355 12 177 21 38
4 455 13 145 22 28
5 525 14 108 23 20
6 510 15 84 24 15
7 450 16 64 25 6
8 365 17 53 26 5
211
Solution
The first step in flood routing is to establish a relationship between reservoir water
elevation and storage of the reservoir. This can be done by slicing the reservoir
horizontally, determining surface area by appropriate method e.g. by planimetring or
using computer program. The volume of each slice can be computed using formulae of
solid geometry. The calculation of volume of reservoir referred here as storage is
calculated in Table 8.3. The formula employed is ∆V = (h/3)[A1+A2+√A1A2] where h is
difference in height between two consecutive slices/contours. The calculations are shown
in table. The elevation Vs surface area is plotted in Figure 8.1.
Noting that slope of curve becomes mild at higher elevations as it should, because higher
elevations can accommodate more volume of water.
The next step is to calculate discharge using equation 8.3A with known coefficient of
discharge and crest level of spillway. The elevation is changed, net head over spillway
found and discharge calculated as shown in Table 8.4
212
Table 8.4 Discharge over Spillway
3/2
Water El. Head Above Crest, H Q = Cd B H
The routing interval is taken equal to time interval of inflow hydrograph which is 1.0
hour. The routing interval should not be that large so as to miss the peak flow
observation. We now have in hand, at a particular elevation storage and outflow value
which enables us to calculate storage indication against each elevation (Table 6.5). A
curve is drawn between storage indication and out flow as shown in Figure 8.3.
213
Storage Indication vs Outflow
600
500
400
Outflow 'm³/s'
300
200
100
0
0 200 400 600 800 1,000 1,200
The last step is to find outflow hydrograph which is done in Table 8.7. The column-3
shows sum of two consecutive inflow values. First value of outflow is considered equal
to inflow. In this example initial level of reservoir is assumed flush with spillway crest so
outflow is zero. The first value of storage is equal to reservoir volume up to spillway
crest. Both being known so [(2S/∆t) - O] is computed. The value of (I j+I j+1) is added to
[(2S/∆t)-O] to get [(2S/∆t)+O]. Corresponding to this value of storage indication outflow
is read from Figure 8.3. The outflow is subtracted twice from storage indication to give
next value of [(2S/∆t)-O]. The process is repeated to get all values of outflow till values
obtained are small enough representing steady flow.
214
Table 8.7 Calculations for Outflow hydrograph
Time Oj
Step (j) Time Inflow, I j I j+I j+1 [(2S/∆t)-O] [(2S/∆t)+O] Outflow
number (Hour) (m³/s) (m³/s) (m³/s) (m³/s) (m³/s)
1 0 50 0 50
2 1 75 125 51 125 37
3 2 175 250 41 301 130
4 3 355 530 27 571 272
5 4 455 810 13 837 412
6 5 525 980 4 993 494
7 6 510 1035 2 1,039 519
8 7 450 960 6 962 478
9 8 365 815 13 821 404
10 9 295 660 21 673 326
11 10 252 547 27 568 271
12 11 212 464 31 491 230
13 12 177 389 35 420 193
14 13 145 322 38 357 159
15 14 108 253 42 291 125
16 15 84 192 45 234 94
17 16 64 148 47 193 73
18 17 53 117 49 164 58
19 18 50 103 50 152 51
20 19 45 95 50 145 47
21 20 39 84 51 134 42
22 21 38 77 51 128 38
23 22 28 66 51 117 33
24 23 20 48 52 99 24
25 24 15 35 53 87 17
26 25 6 21 54 74 10
27 26 5 11 54 65 5
Inflow and outflow hydrographs are plotted simultaneously (Figure 8.4) and following
parameters are found:
o The difference in time on x-axis between peak of inflow and outflow hydrograph.
This is termed as Reservoir Lag.
o The difference between peak discharge of inflow and outflow hydrographs is
calculated. This is termed as attenuation.
215
600
500
400
300
200
100
0
0 5 10 15 20 25 30
The inflow should be routed such that outflow hydrograph is flattened enough having low
peak discharge. This is necessary for safe passages of flood. If abrupt rise in inflow
occurs then special consideration is made for outflow at outlet structures so that sufficient
time is available for operation. This aspect is considered while designing reservoirs.
Let,
S = Storage
I = Inflow
O = Outflow
216
Wedge
Storage
=K(I-O)
O
Prism
Storage
=KO
Figure 8.5 Prism and Wedge Storage in Channel
Then total storage ‘S’ is therefore sum of prism and wedge storage. That is:
S = K [X I + (1 - X) O] - - - - - - - - - - - - - - - - - - - - - - - (8.6)
Where ‘X’ is a dimensionless constant for certain reach or segment of channel. ‘K’ is
storage constant having dimensions of time. Both x and K are determined from inflow
and outflow hydrographs for reach under consideration. These constants vary from reach
to reach and are determined as follows.
1. The inflow and outflow hydrographs are known for the reach. Find values of (I-O)
for each time interval.
2. Find the mean and cumulative mean values of ( I - O) which is storage.
3. Assume value of ‘X’ and find the term [X I + (1-X) O] for each time interval
using assumed value of ‘X’. The storage value is already calculated against time
as explained in step 2.
4. Plot [X I + (1 - X) O] values against storage. Inspect if data plotted nearly fits a
straight line. If not assume new value of X and repeat steps 1-4.
5. The best-fit straight line corresponds to required value of ‘X’. The slope of this
straight line is our required value of ‘K’.
Now we proceed for channel flow routing once values of ‘X’ and ‘K’ are known now.
Routing mean finding outflow hydrograph for given inflow hydrograph. We rewrite
equation 8.6 at beginning and end of jth time interval ∆t.
S j = K I j + (1 - X) O j - - - - - - - - - - - - - - - - - - - - - - - (8.7)
217
S j+1 = K I j+1 + (1 - X) O j+1 - - - - - - - - - - - - - - - - - - - - (8.8)
From equation (8.2), (I j+I j+1)/2-(O j+O j+1)/2=(S j+1-S j)/∆t. Substituting (S j+1-S j) from
equation (8.9) in equation (8.2) and rearranging, we have
Where,
Co = [(∆t/K)-2X]÷[2(1-X)+(∆t/K)] - - - - - - - - - - - - - - - - (8.11)
C1 = [(∆t/K)+2X]÷[2(1-X)+(∆t/K)] - - - - - - - - - - - - - - - - (8.12)
C2 = [2(1-X)-(∆t/K)]]÷[2(1-X)+(∆t/K)] - - - - - - - - - - - - - (8.13)
Example 8.2
Table 8.8 shows Inflow and Outflow hydrographs for certain reach of a channel. Find the
Muskingum’s Coefficients ‘X’ and ‘K’ and route the given inflow hydrograph.
218
Solution
The calculations are best done in tabular form. Refer Table 8.9, wherein first column
shows time interval and second column shows known inflow hydrograph. The known
outflow hydrograph is given in column-4. The routing interval is taken as 12 hours equal
to the interval of inflow hydrograph for computational convenience. However this
interval can be reduced if more precise information is sought. The storage units are 0.50
(m³/s)-day. The change in storage within a interval is found using equation 8.9 as given in
column 6. The cumulative storage is shown in column 7.
Next an appropriate value of ‘X’ is assumed between 0.1 to 0.5 and weighted flow is
determined as XI+(1-X)O for each interval of time as shown in column 8, 12 and 14 of
Table 8.9. The cumulative weighted flow is then plotted against cumulative storage as is
219
shown in figs. 8.6, 8.7 and 8.8. The value of ‘X’ giving a best-fit straight line is the
required Muskingham’s constant X. The slope of this best fit straight line is required ‘K’
X=0.1
800.00
700.00
600.00
Storage 'S' [(m³/s)-day]
500.00
400.00
300.00
200.00
100.00
0.00
0 50 100 150 200 250 300 350 400 450 500
[XI+(1-X)O] (m³/s)
X=0.2
800.00
700.00
600.00
Storage 'S' [(m³/s)-day]
500.00
400.00
300.00
200.00
100.00
0.00
0 50 100 150 200 250 300 350 400 450 500
[XI+(1-X)O] (m³/s)
220
x=0.3
800.00
600.00
Storage 'S' [(m³/s)-day]
500.00
400.00
300.00
200.00
100.00
0.00
0 50 100 150 200 250 300 350 400 450 500
[XI+(1-X)O] (m³/s)
From known values of ‘X’ and ‘K’, constants ‘C0’, ‘C1’, ‘C2’ are determined using
equations 3.11, 3.12 and 3.13.
The outflow hydrograph is computed using equation 8.10 as shown in Table 8.10. The
initial outflow is assumed equal to inflow.
221
Table 8.10 Computations for Channel Routing
Time Constants for X=0.3 and K=1.3
Inflow (Ij) Outflow ‘O'
Interval
(m³/s) C0Ij+1 C1Ij C2Oj (m³/s)
(Hours)
(1) (2) (3) (4) (5) (6)
0 41 41.00
12 36 -4.32 22.55 23.37 41.60
24 39 -4.68 19.80 23.71 38.83
36 123 -14.76 21.45 22.13 28.82
48 342 -41.04 67.65 16.43 43.04
60 577 -69.24 188.10 24.53 143.39
72 720 -86.40 317.35 81.73 312.68
84 742 -89.04 396.00 178.23 485.19
96 675 -81.00 408.10 276.56 603.66
108 450 -54.00 371.25 344.09 661.34
120 322 -38.64 247.50 376.96 585.82
132 247 -29.64 177.10 333.92 481.38
144 190 -22.80 135.85 274.39 387.44
156 145 -17.40 104.50 220.84 307.94
168 120 -14.40 79.75 175.52 240.87
180 96 -11.52 66.00 137.30 191.78
192 79 -9.48 52.80 109.31 152.63
204 67 -8.04 43.45 87.00 122.41
216 55 -6.60 36.85 69.77 100.02
228 50 -6.00 30.25 57.01 81.26
The inflow and outflow hydrograph are shown in figure 8.9. One can find difference in peak
flows and time interval between peak flows which is one of the objectives of channel flow
routing.
Channel Routing by Muskingum Method
800
700
600
500
Flow (m³/s)
400
300
200
100
0
0 50 100 150 200 250
Time (Hours)
Fig 8.9
222
QUESTIONS
1. What is meant by the term Routing? Why we need Routing through reservoirs?
2. Explain Reservoir Routing Procedure. What assumptions are made for it?
3. What information is required prior to routing an inflow hydrograph through a
reservoir?
4. Explain Muskingum’s Method for Channel Routing.
5. What factor governs selection of suitable routing interval?
6. Take a sample data for a reservoir and route a forecast flood through it.
EXERCISE
1. A reservoir has the following elevation, discharge and storage relationships:
When the reservoir level was at 100.50 m, the following flood hydrograph entered the
reservoir.
223
2. The inflow hydrograph for a river reach is as given below:
The values of Muskingum’s Constants are X = 0.30 and K = 18 hours. Find the peak
outflow due to given inflow hydrograph.
224