Lecture15-16 HydrographsUpdated
Lecture15-16 HydrographsUpdated
Lecture15-16 HydrographsUpdated
Day-to-Day
High flows
Hydrology 10000
terms of their
frequency of Median flow
occurence 100
10
1
0 10 20 30 40 50 60 70 80 90 100
Percent$of$time$indicated$that$discharge$was$equaled$or$exceeded
Flow/duration/curve
As watersheds urbanize, peak flows increase – how do we
define natural flows and extreme events with this data?
20000
15000
10000
5000
0 1 3 5 7 9 11 13 15 17 19 21 23 25 27 29 31 33 35 37 39 41 43 45 47 49 51 53 55 57 59 61 63 65 67
Years
As watersheds urbanize, stormwater quality worsens
Permitted (pipe flows), versus non-point, unpermitted, and illegal flows
National Pollutant Discharge Elimination
System (NPDES)
• The EPA NPDES storm water program
regulates stormwater discharges associated
with industrial and construction activities and
municipal storm sewer systems
• The Joint Task Force was created by Harris
County, the City of Houston, HCFCD, and
TxDOT to work together to address stormwater
issues in the Houston Metro Area
• Best Management Practices (BMPs) are the
focus of the effort to reduce contamination in
lakes and rivers from stormwater
Primary Non-Point Sources of
Bacteria to Bayous
High Density Land Uses WW Treatment Plants
Ag / CAFO
Wildlife
am Storm Sewer
Stre
Stream Overflows
Construction Residential
& Pets
Some Models Capable of Continuous Runoff
Calculations including Storm Events (e.g., HSPF)
Figure 1-15
Time-Area Convolution Principle - Watershed
Figure 2-6(a)
Time-Area Convolution – Flow Schematic
Figure E2-2a
Time-Area Convolution – Accumulating runoff from the watershed as a function of time
Figure E2-2c
Time-Area
Method
• Watershed travel times
• Rainfall Intensities
• Resulting Hydrograph
Time Area Hydrograph
Peak Flow at Q2
• Q1 = P1 * A1
• Q2 = P2*A1 + P1*A2
• Q3 =P3*A1 + P2*A2
+ P1*A3
• And So Forth
Development of
design storm
hydrograph
Development of
watershed
hydrograph
UH developed from a single storm
Total storm hydrograph
Hydrograph minus
baseflow, rainfall minus
losses
Hydrograph adjusted
to be a 2-hr UH
Concerns with UH
Assumptions of linearity are inherent in
UH development
Linearity can be violated if…
intensity variations are large over long-
duration storms
Storage effects in watershed are important
Typically should not exceed areas of 3 – 5
mi2 in urban areas, 10 sq mi natural A
Divide the watershed into subareas if needed
Unit Hydrograph Convolution
Deriving hydrographs from multiperiod rainfall excess
or
Where
○ Qn = storm hydrograph ordinate
○ Pi = rainfall excess
○ Uj = UH ordinate
where j = n - i + 1
Unit Hydrograph Convolution
Can view this
graphically
Note that the final
hydrograph goes
to time-step 10
○ UH goes from 0-7
○ 4 Rainfall periods
of 1 time step
UH Convolution Example
Pn= [0.5, 1.0, 1.5, 0.0, 0.5] in
Un= [0, 100, 320, 450, 370, 250, 160, 90, 40, 0] cfs
Storm Hydrograph - Qn
1400
1200
1000
Outflow (cfs)
800
600
400
200
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13
Time (hr)
UH Convolution Example
Pn= [0.5, 1.0, 1.5, 0.0, 0.5 ] in
UH
Time (hr) P 1U n P 2U n P 3U n P 4U n P 5U n Qn
Storm Hydrograph - Qn
1400
0 0 0 1200
Outflow (cfs)
1 50 0 50 1000
800
2 160 100 0 260 600
3 225 320 150 0 695 400
200
4 185 450 480 0 0 1115
0
5 125 370 675 0 50 1220 0 1 2 3 4 5 6 7 8 9 10 11 12 13
Time (hr)
6 80 250 555 0 160 1045
7 45 160 375 0 225 805
8 20 90 240 0 185 535
9 0 40 135 0 125 300
10 0 60 0 80 140
11 0 0 45 45
12 0 20 20
13 0 0
Unit Hydrograph Convolution
Can reverse procedure
Multiperiod rainfall excess hydrograph à UH
Un = 0 0 0 0
[0,100,320,450,370,25 1
2
100
320
50
160
0
100 0
50
260
0,160,90,40,0] 3 450 225 320 150 0 695
4 370 185 450 480 0 0 1115
Calculate storm 5 250 125 370 675 0 50 1220
6 160 80 250 555 0 160 1045
hydrograph 7 90 45 160 375 0 225 805
8 40 20 90 240 0 185 535
9 0 0 40 135 0 125 300
10 0 60 0 80 140
11 0 0 45 45
12 0 20 20
13 0 0
Convert a 1-hr UH to a 3-hr UH
UH=[0,25,50,100,80,60,40,20,0]
Time (hr) UH (cfs) Hr 1 Hr 2 Hr 3 Sum 3-HR UH
tp = Ct*(L*Lc)0.3
= 8.6 hr
Qp = 640*Cp A/tp
= 4465 cfs
Tb ~ 4tp
= 34.4 hr
Duration = tp/5.5 hr
= 1.6 hr
SCS Method - Triangle
Early method assumed a simple triangle
hydrograph with certain parameters
Where
SCS Method - Triangle Lag Time
Lag time (tP) is most often estimated
using formula below – 3000 watersheds
Make sure units match
10 yr – 7.35 in
25 yr – 9.51 in
100 yr – 13.89
in
UH Applications
Design storm inches by duration (10-yr, 25-yr,100-yr)
Developed by the NWS for all urban areas
Duration 10 - yr 25 - yr 100 – yr
3
Hydrograph Combining/Routing
1
Hydrograph Combining/Routing
Subarea G + Subarea F = Combine PT1
11
3
Hydrograph Combining/Routing
Subarea G + Subarea F = Combine PT1
3
Hydrograph Combining/Routing
Subarea G + Subarea F = Combine PT1
3
4
Hydrograph Combining/Routing
Combine PT3 + Route to 4 + Subarea A = RESULTING HYDROG
3 Q
4
T
Channel Routing
• Simulate the movement of water through a channel
Continuity equation
• Used to predict the magnitudes, volumes, and temporal
Hydrologic Routing patterns of the flow (often a flood wave) as it translates
down a channel.
Hydraulic Routing
• 2 types of routing : hydrologic and hydraulic.
Momentum Equation
• both of these methods use some form of the continuity
equation.
56
Flood Routing – Reservoir & River
Tends to shift hydrograph in time and reduce the flow
S S
Q Q
59
Use of Manning Equation
Continuity equation Stage is also related to the
Hydrologic Routing
outflow via a relationship such as
Hydraulic Routing
Manning's equation
Momentum Equation
1.49
Q= ARh S f
2 1
3 2
n
60
Hydraulic Routing
Continuity equation • Hydraulic routing methods combine the continuity
equation with some more physical relationship
Hydrologic Routing describing the actual physics of the movement of
Hydraulic Routing the water.
• The momentum equation is the common
Momentum Equation relationship employed.
• In hydraulic routing analysis, it is intended that the
dynamics of the water or flood wave movement be
more accurately described
61
Routing Methods
• Modified Puls
Modified Puls
Kinematic Wave • Kinematic Wave
Muskingum
• Muskingum
62
Modified Puls
• The modified puls routing method is probably most often
Modified Puls applied to reservoir routing
Kinematic Wave • The method may also be applied to river routing for
Muskingum certain channel situations.
Muskingum-Cunge
• The modified puls method is also referred to as the
storage-indication method.
• The heart of the modified puls equation is found by
considering the finite difference form of the continuity
equation.
63
Modified Puls
Modified Puls I 1 + I 2 - ( O1 + O 2 = S 2 - S 1 Continuity Equation
Kinematic Wave
Muskingum 2 2 Dt
Muskingum-Cunge
Dynamic æ 2 S1 ö 2 S2
Modeling Notes I1 I2 ç
+ + - O1 ÷ = + O2 Rewritten
è Dt ø Dt
64