Green Earth
Green Earth
Green Earth
Flow Routing Q
Q
Procedure to
determine the flow
hydrograph at a t
point on a Q
watershed from a
known hydrograph
upstream t
As the hydrograph Q
travels, it
attenuates t
gets delayed
Why route flows?
dS
t I O
dt t t
Hydrologic river routing
(Muskingum Method)
Wedge storage in
reach Advancing
Flood
S Prism KQ Wave Q
I>Q I
S Wedge KX ( I Q)
I Q
K = travel time of peak through the reach
X = weight on inflow versus outflow (0 ≤ X ≤ 0.5) Q Q
X = 0 Reservoir, storage depends on outflow,
no wedge
X = 0.0 - 0.3 Natural stream
Receding
Flood I Q
S KQ KX ( I Q) Wave
Q>I
S K [ XI (1 X )Q] QI
I I
Continuity Equation in
Difference Form
Referring to figure, the continuity
equation in difference form can be
expressed as DS S2 S1 _I O_ (I1 I2 ) (O1 O2 )
Dt t t 2 2
2 1
Derivation of Muskingum Routing Equation
• By Muskingum Model,
at t = t2, S2 = K [X I2 + (1 - X)O2]
at t = t1, S1 = K [X I1 + (1 - X)O1]
• Substituting S1, S2 into the continuity equation and after some
algebraic manipulations, one has
O2 = Co I2 + C1 I1 + C2 O1
• Replacing subscript 2 by t +1 and 1 by t, the Muskingum routing
equation is
Ot+1 = Co It+1 + C1 It + C2 Ot, for t = 1, 2, …
D KX 0.5Dt
where C KX 0.5 t ; C ; C2 = 1 – Co – C1
o K KX 0.5Dt 1 K KX 0.5Dt
Routing
Muskingum Routing
Equation
Q2 C0 I 2 C1 I1 C2Q1
Qt 1 C0 I t 1 C1 I t C2Qt
C0 = (– Kx + 0.5Dt) / D
C1 = (Kx + 0.5Dt) / D
C2 = (K – Kx – 0.5Dt) / D
D = (K – Kx + 0.5Dt)
Graphical Method:
Referring to the Muskingum Model, find X
such that the plot of XIt+ (1-X)Ot (m3/s) vs
St (m3/s.h) behaves almost nearly as a single
value curve. The assume value of x lies
between 0 and 0.3.
The corresponding slope is K.
Example 8.4: Estimating the
value of x and K.
Try and error to get the nearly straight
line graph.
Muskingum Routing
Procedure
Given (knowns): O1; I1, I2, …; Dt; K; X
Find (unknowns): O2, O3, O4, …
Procedure:
(a) Calculate Co, C1, and C2
(b) Apply Ot+1 = Co It+1 + C1 It + C2 Ot
starting from t=1, 2, … recursively.
Example 8.5
Given K and x.
Initial outflow, Q also given.
Solution:
Calculate Co, C1, and C2
C0 = (– Kx + 0.5Dt)/ D
C1 = (Kx + 0.5Dt)/ D
C2 = (K – Kx – 0.5Dt)/ D
D = (K – Kx + 0.5Dt)
Solution:
Route the following flood hydrograph through
a river reach for which K=12.0hr and X=0.20.
At the start of the inflow flood, the outflow
flood, the outflow discharge is 10 m3/s.
Time 0 6 12 18 24 30 36 42 48 54
(hr)
Inflow 10 20 50 60 55 45 35 27 20 15
(m3/s)
Reservoir Routing
dS
I Q
dt
Inflow and Outflow
I1 + I2 – Q1 + Q2 S2 – S1
=
2 2 Dt
Inflow & Outflow Day 3
= change in storage / time
S3 S2
I 2 I 3 / 2 Q2 Q3 / 2 dt
Repeat for each day in progression
Determining Storage
• Evaluate surface area at several different depths
• Use available topographic maps or GIS based DEM
sources (digital elevation map)
• Outflow Q can be computed as function of depth for
either pipes, orifices, or weirs or combinations
Combined
S
Pipe
Q
Comparisons:
River vs.
Reservoir
Routing
River Reach
Flood Control
Structural Measures
Non-structural Methods
Structural Measures
Storage and detention reservoir
Flood ways (new channel)
Levees (flood embankment)
Channel Improvement
Reservoirs
Reservoirs reduce flooding by temporarily
storing flood waters behind dams or in
storage or detention basins.
Reservoirs lower flood heights by holding
back, or detaining, runoff before it can flow
downstream.
Flood waters are detained until the flood has
subsided, then the water in the reservoir or
detention basin is released or pumped out
slowly at a rate that the river can
accommodate downstream.
Timah Tasoh Dam
Reservoirs
Flood control reservoirs are most commonly
built for one of two purposes. Large
reservoirs are constructed to protect property
from existing flood problems.
Smaller reservoirs, or detention basins are
built to protect property from the impacts of
new development (i.e., more runoff).
Think!