Reservoir Flood Routing (Hydrograph Routing) 12 PDF
Reservoir Flood Routing (Hydrograph Routing) 12 PDF
Reservoir Flood Routing (Hydrograph Routing) 12 PDF
, NUST 5/3/2007
Chapter 12
Introduction
The design of dams involves the sizing of appurtenances like the spillway, outlets
and determining the storage capacity. It is a complex phenomenon and if not
professionally done, can lead to catastrophic dam failures. Of late in Zimbabwe,
dam failures are becoming common and mainly due to poor designs and
workmanship. In most cases it will be failure of the spillway to cope with the peak
discharge which may occur during the design period based on the recurrence
interval. With regard to this, the hydrology of proposed dam sites needs to be
studied. This data is necessary to carryout flood routing which will enable sizing
the spillway and other dam appurtenances.
Flood routing
A process that shows how the inflow hydrograph during peak discharge (flood
wave) can be reduced in magnitude and lengthened in time (attenuated) by use
of storage in the reach between two points. It predicts the temporal and spatial
variation of a flood wave as it traverses a river reach or reservoir. It can be used
to predict the outflow hydrograph from a catchment area. The resultant outflow
hydrograph is then used to size the spillway.
• Hydrologic
• Hydraulic
Hydrologic
It employs the equation of continuity with either an analytical or an assumed
relation between storage and discharge within the system. Hydrologic routing is
useful for flood prediction and evaluation of flood control measures.
Hydraulic routing
It uses the equation of continuity and the equation of motion (momentum). It
utilizes the partial differential equation for unsteady flow in open channels. Its
advantage is that it adequately describes the dynamics of flow than the
hydrologic technique.
Translation
Attenuation
Q
S=I-O
Outflow
Outflow
The outflow hydrograph is governed by the inflow hydrograph, weir and storage
characteristics. The procedure is summarized in the flow chart in Fig. 2.
Storage characteristics
This is the relationship between storage and height of the dam basin. A plot of
height against storage gives a capacity curve.
Weir characteristics
The shape of the weir (cill) of the spillway affects the overflow efficiency. For
small dams the spillway is normally stepped or sloped to the downstream with a
vertical face upstream. The efficiency of a weir is measured by the coefficient of
discharge. In practice the discharge coefficient varies between 1.7-1.84 (Baban,
1995). The discharge over a an un-gated spillway weir section can be described
from energy, momentum and continuity consideration by the form:
Q = C d LH x (1)
3
x =exponent, theoretically
2
Input parameters
Weir characteristics
Storage, S Outflow, Q
H
Spillway design
Weir
Inflow hydrograph
The problem of flood routing is solved by applying the hydrologic equation:
I = Q + ∆S (2)
ds
⇒ I =Q+ (3)
dt
S is often defined as the surcharge storage or the storage above the spillway
crest i.e. outflow is zero when S=0.
(I t + I t + ∆t ) (Qt + Qt + ∆t )
∆t − ∆t = S t + ∆t − S t (4)
2 2
The routing period, t selected should be sufficiently short such that the
hydrograph during the interval t to t + ∆t can be assumed as a straight line, i.e.
I t + I t + ∆t
I mean = (5)
2
(I t + I t + ∆t ) Qt Q
∆t + S t − ∆t = S t + ∆t = S t + ∆t + t + ∆t ∆t (6)
2 2 2
(I t + I t + ∆t ) S t Qt S t + ∆t Qt + ∆t
+ − = + (7)
2 ∆t 2 ∆t 2
The above equation is difficult to solve since both the S and Q can be complex
functions of H, the head over the weir, and therefore a semi-graphical method is
used. The equation is solved successively through time for a known inflow
hydrograph if the storage volume can be related to outflow.
The storage-elevation relationship is computed from the basin survey and the
stage-discharge relationship from the weir equation (Figs. 3 and 4).
If annual river discharges are available, the inflow hydrograph can be derived but
on many occasions this data is not available and the simplest method of deriving
the inflow hydrograph is by applying the Triangular Model Hydrograph and the
Rational Formula. The procedure involves computing:
storage, S
Fig. 3
Fig. 4
Example
(i) Derive the inflow and outflow hydrographs using data from previous
examples for a reservoir with a 36m long spillway and vertical sides
and plan area of 5000m2. The storage is given as S = 50 000H.
(ii) Estimate the peak outflow discharge
Solution
Calculate the maximum probable flood (previous example).
(i) MPF=65m3/s)
(ii) Rainfall intensity from equation previous equation
5
65 × 6.87
I= =45mm/hr
10
1 2 3 4
T/Tp T, hrs Q/Qp Q,
ratio m 3/ s
0 0 0 0
0.1 0.138 0.015 0.975
0.4 0.552 0.280 18.20
0.6 0.828 0.600 39.00
0.8 1.104 0.890 57.85
1.0 1.380 1.000 65.00
1.4 1.932 0.750 48.75
1.8 2.484 0.420 27.30
2.2 3.036 0.240 15.60
3.0 4.140 0.075 4.875
3.5 4.830 0.036 2.340
4.0 5.520 0.018 1.170
4.5 6.210 0.009 0.585
5.0 6.900 0.004 0.260
Procedure
70
60
50
Discharge, Q(m3/s)
40
30
20
10
0
0 1 2 3 4 5 6 7 8
-10
Time, hrs
(iii) Using the weir equation (1) calculate the flow rates, Q, associated with
height H (column 3, Table 3).
(iv) Compute column 4 and 5 in Table 3. The time interval is taken as half an
hour. Since Q is in cubic meters per second, then T must be in the same
units, T=0.5X3600=1800 seconds. Column 4 and 5 are computed by
putting the values of Q and S in the equations: S/T+Q/2 and S/T-Q/2.
(v) Plot graphs, Q against S/T+Q/2, and Q against S/T-Q/2 (Fig. 6)
1 2 3 4 5
H, m Storage Flow rate, S/T- S/T+Q/2
, Q Q/2
(m) (m3 ) ( m3)
0 0 0 0 0
0.1 15 000 2.05 7.31 9.36
0.2 30 000 5.80 13.77 19.56
0.3 45 000 10.65 19.68 30.32
0.4 6 0000 16.39 25.14 41.53
0.5 75 000 22.91 30.21 53.12
0.6 90 000 30.12 34.94 65.06
0.7 105 000 37.95 39.36 77.31
0.8 120 000 46.37 43.48 89.85
0.9 135 000 55.33 47.34 102.66
1.0 150 000 64.80 50.93 115.73
8
140
120
100
S/T-Q/2 & S/T+Q/2
80
S/T-Q/2
S/T+Q/2
60
40
20
0
0 10 20 30 40 50 60 70
Discharge, m 3/s
(iv) Determine the inflow, column 2, Table 4. Since the time interval has been
taken as 0.5hrs (routing period), then the new inflow ordinates are taken
from Fig. 5 (inflow hydrograph) for the respective half-hour intervals
(column 1 and 2, Table 4).
(v) Compute the mean inflow for each time interval, (In+In+1)/2 (column 3)
e.g. (0+18)/2 =9m3/sec
(vi) Compute the outflow, column 4 as follows. The calculations start at
interval t=0 and proceed step by step:
Step1
The outflow at the start of interval is 0m3/s. Referring to Fig. 6 when Q=O, S/T-
Q/2 =0
Step2
In Table 4, column 6, the value from step one is added to the mean inflow i.e.
9+0=9m3/s and recorded in column 6.
Step3
From Fig. 6, when S/T+Q/2 =9m3/s, the outflow value is 2m3/s (column 4)
Step4
From Fig. 6 when the outflow value is 2m3/sec, the value of S/T-Q/2 is 8, column
5. Steps 1, 2,3, and 4 are thus repeated giving the outflow hydrograph ordinates
in column 4. These ordinates are plotted against time as shown in Fig. 7. From
the outflow hydrograph the peak outflow discharge is 46m3/sec.
From the weir equation (1) the maximum height, Hmax corresponding to the
peak outflow can be determined:
2
⎡ 46 ⎤ 3
H max =⎢ ⎥ =0.796m
⎣1.8 x36 ⎦
The Hmax is the called wet free board and a dry freeboard is normally added to
the wet free board to get the crest level of the dam wall. The whole process can
be repeated again with a new spillway length instead of 36m. A larger spillway
length will automatically reduce the wet freeboard but the maximum length is
governed by site conditions.
1 2 3 4 5 6
Time, Inflow Mean Outflow, S/T-Q/2 S/T+Q/2
hrs inflow Q
(m /s) (m3/s)
3
(m3/s)
0 0 0 0
+
9 9
0.5 18 2 8
37.9 45.9
1 57.8 19 28
60.9 88.9
1.5 64 40 40
56.4 96.4
2 48.8 44 42
38.05 80.05
2.5 27.3 40 40
21.45 61.45
3 15.6 26 32
12.8 44.8
3.5 10 18 26
7.5 33.5
4 5 12 20
3.65 23.65
4.5 2.3 8 16
2.05 18.05
10
5 1.8 5 12
1.5 13.5
5.5 1.2 3 9
0.89 9.89
6 0.58 2 8
0.34 8.34
6.5 0.1 1.5 6
70
60
50
Discharge, Q(m3/s)
40
In flow
O u tflow
30
20
10
0
0 1 2 3 4 5 6 7
T im e, h rs
• prism storage
• wedge storage (Fig. 1).
11
Prism storage
The volume of water that would be stored in the reach if the flow was uniform
throughout.
Wedge storage
It is due to the non uniformity during the flood, when the inflow is different from
the outflow and is the volume of water stored due to the actual water profile and
depends on the difference of Inflow (I) and Outflow (O).
The Stream flow routing is much complicated than reservoir routing because the
volume is not simply a function of stage. This is because:
A common method of flow routing is the Muskingum method. The storage within
a reach at any given time can be expresses as:
b ⎡ mn
XI + (1 − X )O n ⎤
m
S= (1)
a ⎢⎣ ⎥⎦
12
Hydrologic routing:
ds
I −O = (3)
dt
I 1 + I 2 O1 + O2 S 2 − S1
− = (4)
2 2 ∆t
S 2 − S1 = K [X (I 2 − I 1 ) + (1 − X )(O2 − O1 )] (5)
Substituting equation 5 in 4:
I 1 + I 2 O1 + O2 K [ X (I 2 − I 1 ) + (1 − X )(O2 − O1 )]
− = (6)
2 2 ∆t
O2 = C O I 2 + C1 I 1 + C 2 O1 (7)
− KX + 0.5∆t
in which: Co = (8)
K − KX + 0.5∆T
KX + 0.5∆t
C1 = (9)
K − KX + 0.5∆T
K − KX − 0.5∆t
C2 = (10)
K − KX + 0.5∆T
13
Check: C O + C1 + C 2 = 1 (11)
• assume a value of X
• compute the expression [XI + (1 − X )O ] )
• plot values against cumulative storage to obtain the storage loop
• repeat steps above steps until a linear relationship is established
• K is the gradient of the straight line obtained
The routing period should be less than the time interval for the flood wave
through the reach: t ≅ 1 − 1 of the flood wave travel time through the
3 4
reach from the inflow hydrograph (Raghunath, 1991).. The routing period must
be expressed in equal increments just like in reservoir routing.
Example 1
The inflow and outflow hydrographs for reach of a river are given Table 1.
Determine the value of Muskingum coefficients K and X for the reach.
Table 1
Solution
• From the inflow and outflow hydrographs, a routing period t-=24 hr=1 day
is taken.
• The storage S is computed, S = I −O
S + S2
• Mean storage is computed from equation S mean = 1
2
• cumulative storage is computed
• for trail values of X=0.2, 0.25, 0.3 the values of [XI + (1 − X )O ] are
computed
14
• storage loops are plotted for S vs. [XI + (1 − X )O ] for each trial value (Fig.
2)
• the value of 0.25 approximates a straight line and hence is chosen as the
value of x
• K is the gradient of the straight line
• hence for the given reach the value of Muskingum coefficients are X=0.25
and K=0.7 day
15
Example 2
The inflow hydrograph for a stream are given in Table 4 for which the Muskingum
coefficients of K=36 hr and x=0.15. Rout the flood through the reach and
determine the outflow hydrograph. Outflow at the beginning of the flood may be
taken as the same as the inflow.
Table 4
Solution
16
From equation 7 O2 = C O I 2 + C1 I 1 + C 2 O1
KX + 0.5∆t
C1 = =0.31
K − KX + 0.5∆T
K − KX − 0.5∆t
C2 = =0.67
K − KX + 0.5∆T
• CHECK: 0.02+0.31+0.67=1
• O1 is assumed equal to 42m3/s
• O2 = 0.02 × 45 + 0.31 × 42 + 0.67 × 42 = 42.06m 3 / s
• this value of O2 becomes O1 for the next routing period and the process is
repeated till the flood is completely routed through the reach (Table 5).
The resulted outflow hydrograph is plotted (Fig. 3). The reduction in the
peak flood is 108m3/s and the lag time is 36hrs i.e. the peak flow is 84hrs
(3.5 days) after the commencement of the flood through the reach.
References
1. Mitchell T.B., (1998), Annual River Flood Frequency in Zimbabwe, the
Zimbabwe Engineer, Vol. 63 N10-1 pg.9, Thompson Publications, Zimbabwe
2. DWD, (1972), An Assessment of the Surface Water Resources of Zimbabwe
and Guidelines for Development and Planning, Ministry of Water Resources
and Development, Zimbabwe Government
3. Barban R., (1995), Design of Diversion Weirs, John Wiley, UK
4. Kabell T., (1987), Guidelines to the Design of Dams, DWD Annual
Conference, Zimbabwe Government
5. Kabell T., (1983), Reservoir Sedimentation in Zimbabwe, Ministry of Water
Resources and Development, Zimbabwe Government
17
18