The Final Project Report Water Engineering
The Final Project Report Water Engineering
1. Introduction
1.1 Aims & Objectives
1.2 Objectives
2. Literature Review
3. Methodology
3.1 Selection of Site
3.2 Time of concentration and rainfall pattern generation
3.3 Unit hydro graph
3.4 Network diagram
3.5 Preparation of predevelopment hydrograph
3.6 Channel routing
3.7 Developing a sub catchment
3.8 Generating rainfall pattern for developed land
3.9 Unit hydrograph
3.10 Reservoir routing
3.11 Preparation of HEC HMS
5. Conclusion
6. References
1. Introduction
In the rapid growth of the construction industry it is a challenge to meet the council
requirement to execute a proposed project. The area of land undertaken for development in
this report is under the Muswellbrooke council in new south wales, Australia. Total area of
2.63km2 catchment was considered which was further divided into 3-sub catchments. One of
the 3 sub catchments were considered for the development, and a detailed hydrologic
analysis was carried out mainly based on pre-development and post the development.
Hyetographs were generated; unit hydrograph and direct runoff hydrographs for each of the
sub catchments were generated. The flow was routed through the channel over the different
junction and an outflow hydrograph is generated at the end. Soft wares like HEC HMS was
used to verify the results obtained from long hand calculations and to gain more accuracy.
Development activity will reduce the losses and will contribute to generation of more surface
runoff. To compensate the extra volume of runoff water, a hydraulic structure has been
designed so that the stability of the land and the structure is assured. Comparison of the pre
and post development is vital to assure that the project is safer, and to get an approval at the
end.
1.2 Objectives
According to McCuen, R.H., Wong time to concentration is defined as the time needed by water
after rainfall to reach the watershed from the most remote point. It very much depends upon the
topographical features of the land under the watershed. For designing of any hydrological
project time to concentration is a very important parameter to find out the peak discharge.
Bransby Williams equation is used to find out the time of concentration for full catchment.
However, we have used the other formula for calculating the time of concentration for
subdivided catchments.
Bransby William equation says that the minimum or maximum flood for a catchment from
produced rainfall is the time taken by the water from rainfall to reach the discharge point from
the distance point in catchment.
Cudennec, C., Fouad, Y., Sumarjo Gatot showed in their results that unit hydrograph can be
defined as the direct runoff hydrograph from 1 cm effective rainfall occurred uniform over the
entire catchment at uniform rate for specified duration.Cudennec,2004,For deriving the unit
hydrograph sum up the volume of direct runoff which is done by adding all the direct runoff
values which is further multiplied with the time interval. Synthetic unit hydrograph is defined as
the hydrograph drawn upon the assumptions of the other hydrograph. Gray,1961, SUH is based
on the catchment storage and it is made using dimension less unit hydrograph.
It was stated in research paper of Duchesne, J., The SCS triangular Unit Hydrograph is used for
generation of the catchment response for 1mm depth of rainfall excess. Unit hydrograph was
calculated for each of the sub catchments. After getting the area of total and subdivided
catchment, derive the unit hydrograph for three sub-catchments from which the 5- minute unit
hydrograph is derived for each sub-divided catchment.
In the study of Gray, D.M., S hydrograph is a kind of summation hydrograph used to adjust the
duration of unit hydrograph. For developing a S- unit hydrograph one should know the storm
duration. S hydrograph allow us to generate adjustment to existing unit hydrograph. S
hydrograph can be defined as the hydrograph generated from uniform rate occurring rainfall for
an indefinite period of time .S hydrograph is used to generate a hydrograph of the duration
which isn’t the integral multiple of the available unit hydrograph. S hydrograph is used where
principle of superposition is not implemented to calculate the unit hydrograph. The reason for
calling it S hydrograph because of the shape of hydrograph is very similar to alphabet S.
Dooge and J.C. says that Muskingum method of channel routing is used to estimate
downstream hydrograph from upstream hydrograph. Fenton, 1989, This method is used to
estimate the floods, therefore also known as Muskingum method of flood routing. Flood routing
is an important aspect of hydrology which progressively determines the timing and magnitude of
a flood wave along a river stream. Muskingum method use two parameters X and K for finding
out the discharge
K= gradient of storage
investi Kožar, I., Lozzi-Kožar, D. come to conclusion that Utilization of the mass preservation
guideline encourages the introduction of fundamental and differential details dependent on
either the water volume or the water dimension of a self-assertively molded repository.
Reservoir routing is done for estimating the water run off at any place which is very useful in
many surveys before staring any big project.
3. Methodology:
3.1 Selection of site:
The site was selected based on the Land developers request. Using the Six maps
software’s topographic view and a location was selected in the Muswellbrook shire. Using
the rule of drawing perpendicular bisectors to the topographic lines a map was drawn of the
catchment area by using the distance tool of the Six maps. Detailed hydrologic analysis is to
be carried out on the catchment selected.
Figure 1: CATCHMENT
The catchment area was calculated using the Six maps software Area tool. Additionally,
using the coordinate tool coordinates of the catchment was found.
To get council approval for the catchment selected the area is to be divided into sub
catchments. By obeying the same rule of perpendicularly bisecting the topographic lines the
catchment was divided into 3 sub catchments. Following the same steps as mentioned
above the area for each of the sub catchments was also found.
DESCRIPTION AREA
Using the equal Area slope found out by the help of Elevation Finder. The time to
concentration of the site was calculated by using Bransby Williams equation.
Knowing the Area of the catchment the Length of the main stream and the equal area slop
the time to concentration was fond out to be 58 mins. The storm duration was taken as 60
mins because in the ARR data hub storm patterns of duration 60 min is available.
Time to concentration of each of the sub catchments were also calculated using the formula
tc =0.76 A0.38
The time of concentration of the sub catchments were less than the time of concentration of
the entire Catchment. Which proved as a check for calculation.
Table 2: values approximated to the closest whole number for ease in calculation
Using the ARR data hub knowing our catchments coordinated we generated the Temporal
distribution of rainfall on our catchment for 60 min storm duration.
The council of Muswellbrook Shire requirement for developing a land with residential
building included hydrologic analysis for storm of 1 in 5-year occurrence.
From the 60 min storm temporal distribution frequent storms were taken and an average
value was taken as the percentage of precipitation of 60 min over 5 min intervals.
Knowing the duration of the storm and the frequency of the storm using BOM IFD curves the
precipitation in mm was found out. From where the storm hyetograph was plotted for the
catchment.
Rainfall
Time Rainfall depth Initial Loss Continuos Rainfall excess Rainfall excess
intensity(mm/
(min) (mm) (mm) Loss (mm/hr) Rate (mm/hr) Depth (mm)
hr)
0-5 7.86 94.29 15.16 0.00 0.00
5-10 7.30 87.57 0.00 0.00
10-15 7.03 84.33 2.50 81.83 6.82
15-20 4.46 53.48 2.50 50.98 4.25
20-25 4.38 52.61 2.50 50.11 4.18
25-30 4.55 54.55 2.50 52.05 4.34
30-35 5.17 61.99 2.50 59.49 4.96
35-40 3.59 43.08 2.50 40.58 3.38
40-45 5.28 63.42 2.50 60.92 5.08
45-50 6.35 76.19 2.50 73.69 6.14
50-55 7.75 93.03 2.50 90.53 7.54
55-60 4.69 56.27 2.50 53.77 4.48
Rainfall Depth
(mm)
9.00
8.00
7.00
Intensity mm/hr)
6.00
5.00
4.00
3.00
2.00
1.00
0.00
0-5 5-10 10-15 15-20 20-25 25-30 30-35 35-40 40-45 45-50 50-55 55-60
Time (min)
Figure 6: Hyetograph
By using the initial and continuous loss model taking account of the infiltration losses the
rainfall excess hyetograph is plotted.
Rainfall excess
Depth (mm)
8.00
7.00
6.00
Intensity (mm/hr)
5.00
4.00
3.00
2.00
1.00
0.00
0-5 5-10 10-15 15-20 20-25 25-30 30-35 35-40 40-45 45-50 50-55 55-60
Time (min)
0.300
0.250
0.200
0.150
0.100
0.050
0.000
0 10 20 30 40 50 60
(tc = 50 min)
5 min UH
0.5
0.45
0.4
0.35
0.3
0.25
0.2
0.15
0.1
0.05
0
0 10 20 30 40 50 60 70
(tc = 55 min)
0 0 0 0 0
5 0.04 0.04 0 0.04 0.44
10 0.08 0.08 0.04 0.04 0.44
15 0.12 0.12 0.08 0.04 0.44
20 0.16 0.16 0.12 0.04 0.44
25 0.2 0.2 0.16 0.04 0.44
30 0.24 0.24 0.2 0.04 0.44
35 0.28 0.28 0.24 0.04 0.44
40 0.32 0.32 0.28 0.04 0.44
45 0.36 0.36 0.32 0.04 0.44
50 0.4 0.4 0.36 0.04 0.44
55 0.44 0 0.44 0.4 0.04 0.44
60 0.4 0.04 0.44 0.44 0 0
0.2
0.15
0.1
0.05
0
0 5 10 15 20 25 30
(tc = 20 min)
Figure 12: Unit hydrograph; SUB CATCMENT 3
Unit
20 min Unit 5 min Unit
Time Hydrograp
Shift Hydrograph Shift Difference Hydrograph
(min) h
(m3/s) (m3/s)
(m3/s)
0 0 0 0 0
5 0.0497 0.0497 0 0.0497 0.1988
10 0.0994 0.0994 0.0497 0.0497 0.1988
15 0.1491 0.1491 0.0994 0.0497 0.1988
20 0.1988 0 0.1988 0.1491 0.0497 0.1988
25 0.1491 0.0497 0.1988 0.1988 0 0
Using the approach of superposition and linearity we plotted the flood hydrographs of each
sub catchments preparation for the Outflow Hydrograph
Where;
Subject to condition of C0 + C2 + C3 = 1
8.00
7.00
6.00
4.00
3.00
2.00
1.00
0.00
1 2 3 4 5 6 7 8 9 10 11 12
After
Devolupment
tc2 30 min
Area (80%) 0.27 km2
5 min UH
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0 5 10 15 20 25 30 35 40
First, we need to give the name of our project from the file after selecting new in left of screen.
After that, select the Basin model manger from components by which the system allows to draw
your sub-catchments.
In this project we have made 3 sub-catchments and named them sub-catchment 1, sub-
catchment 2 and sub-catchment 3. Sub-catchments are made with the help of sub-basin
creation tool.
After drawing the sub-catchments, from the reach creation tool reach1, reach2 and reach 3 are
made. The purpose of drawing the reach is to tell system where the sub-catchments are flowing
into.
The next step is to make junctions where the sub-catchments ends into. Therefore, for sub-
catchment 1 and subjection 2, junction A is made. From junction A the water flow goes into
junction B. the purpose of making junction is to tell system where the flow of sub-catchment is
meeting into.
For each sub-catchment, reach and junction there is data required to be filled before moving to
the other components of project. After clicking on sub-catchment, the system will ask for the
area of that catchment and in this project, transform method is selected to unit specified
hydrograph method and baseflow method to none and this is set according to the need.
Here it is very important to give downstream of every sub-catchment, reach and junction to get
accurate results from the system.
In this work, the routing is done at reach 3 with Muskingum method by giving the system values
of K and X derived from the work. Value of k is 0.1166 and weighting factor(x)=0.2
Another component which time-series data manager is created from the component tool. In
here, gage 1 is named for telling the system units, data source and time interval. Here system
ask about the starting and ending date of rainfall with starting and ending time of storm which is
8 may2019 and time duration is 1 hour in this project. Rainfall depth in mm is filled into the table
option next to time window for time-series gage.
Figure 22: Adding rainfall pattern to the catchment
Select the next component which is paired data manager. The data type for this is unit
hydrograph curves and total of three hydrographs are made named unit hydrograph 1, 2 and 3.
For each unit hydrograph the value we got for the 5-min unit hydrograph is filled into time and
discharge table for each hydrograph. For each sub-catchment, from the transform select the unit
hydrograph for their catchment respectively.
Figure 23: Adding unit-hydrographs for each sub-catchments
From meteorological data manger a new section is made named met-1. After selecting it
hydrograph system will ask to include sub-basins for which yes option is selected. Then from
the specified hydrograph for each sub-catchment the created gage1 is selected.
Figure 24: Assigning rainfall patterns to each sub-catchment
In control specification manager a new section named Control 1 is made. Again, the system will
ask for starting and ending date, starting- ending time of storm and time interval. The useful
data is filled into the required options. However, in this step the time duration is increased from
1 hour to 2:30 hour in our project.
At last from the above options, compute is selected from where create compute is opted for
simulation run.
Figure 26: Simulation runs
Project will only run if all the input data is correct. If it is okay, then a notification bar of success
will appear on the screen.
4. RESULTS AND ANALYSIS:
4.1 Pre-Dev Hydrographs:
16.000
14.000
12.000
10.000
8.000
6.000
4.000
2.000
0.000
0 20 40 60 80 100 120 140 160 180
20.000
15.000
10.000
5.000
0.000
0 20 40 60 80 100 120 140 160 180
4.500
4.000
3.500
3.000
2.500
2.000
1.500
1.000
0.500
0.000
0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95
The DRH from Sub Catchment 1 and Sub Catchment 2 meets at point A where they
combine. Summation of the DRH from the two sub catchments gives us the DRH at the
point A.
DRH at A
45.000
40.000
35.000
30.000
25.000
20.000
15.000
10.000
5.000
0.000
0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 105 110 115
40.00
35.00
30.00
25.00
20.00
15.00
10.00
5.00
0.00
0
5
10
20
30
35
40
50
55
65
75
80
85
90
95
15
25
45
60
70
0
0
0
5
5
10
11
12
13
14
15
10
11
12
13
14
15
Hydrograph at A Hydrograph at B
The routing is achieved by using the following parameters in the Muskingum equations for
channel routing.
Dt 5 min
K 7 min
x 0.2
TR 60 min
Dt <= 0.25 TR 0.25 TR 15 min
Dt <= K K 7 min
Dt >= 2 K x 2K x 2.8 min
Dt 5 min
40.00
35.00
30.00
25.00
20.00
15.00
10.00
5.00
0.00
0
5
60
70
75
80
85
90
95
10
15
20
25
30
35
40
45
50
55
65
0
5
0
5
0
5
0
5
0
5
0
5
10
10
11
11
12
12
13
13
14
14
15
15
5 min UH
12.00
10.00
8.00
6.00
4.00
2.00
0.00
0 2 4 6 8 10 12
Using the previous concepts, the post developed DRH will add with the DRH form sub
catchment 1 and point A using simple summation as following.
Post Dev Hydrograph at A
50.000
45.000
40.000
35.000
30.000
25.000
20.000
15.000
10.000
5.000
0.000
0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 105 110
Using the Muskingum method for channel routing the hydrograph at A will be routed over
channel length of 200 meters to get the routed hydrograph at B.
50
45 Chart Title
40
35
30
25
20
15
10
0
0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100105110115120125130135140145
Hydrograph at A Hydrograph at B
Figure 36: Post Dev Routed Hydrograph at B
TR 60 min
Dt <= 0.25 TR 0.25 TR 15 min
Dt <= K K 7 min
Dt >= 2 K x 2K x 2.8 min
Dt 5 min
60
70
75
80
85
90
95
10
15
20
25
30
35
40
45
50
55
65
0
5
0
5
0
5
0
5
0
5
6
10
10
11
11
12
12
13
13
14
14
14
Routed hydrogrpah Sub catchment 3 Final outflow Hydrograph
Figure 38: Post Dev Outflow Hydrograph
Pre dev & Post Dev Hydrographs
50.00
45.00
40.00
35.00
30.00
25.00
20.00
15.00
10.00
5.00
0.00
0
5
30
35
75
80
10
15
20
25
40
45
50
55
60
65
70
85
90
95
0
5
0
5
0
5
0
5
0
5
0
5
10
11
14
14
15
10
11
12
12
13
13
15
PRE dev POST dev
Figure 39: Comparison between Pre and Post development Hydrographs
This is the time series table which shows the time and date of storm. It shows the value of
precipitation, loss, excess, direct flow, base flow and total flow values.
Figure 47: the hydrograph obtained from sub catchment run off at 1
REACH 2
REACH 3
Routing is done using Muskingum method filling the value of k=0.1166 and x=0.2
JUNCTION B
45.00
40.00
35.00
30.00
25.00
20.00
15.00
10.00
5.00
0.00
0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 10 10 11 11 12 12 13 13 14 14
0 5 0 5 0 5 0 5 0 5
Dimensions of Reservoir
(Trapezoid)
side slope (m) = 1:0.2
A = 70 m
B = 90 m
Depth = 2.0 m
Death water height = .2 m
Length = 230 m
vol of pond = 36300 m3
Height (m) length (m) width (m) Area (m2) Vol (m3)
0 35 10 350 0
0.2 35.12 10.12 355.4144 70.54144
0.4 35.36 10.36 366.3296 144.3488
0.6 35.72 10.72 382.9184 224.7744
0.8 36.2 11.2 405.44 315.3434
1 36.8 11.8 434.24 419.84
1.2 37.52 12.52 469.7504 542.3942
1.4 38.36 13.36 512.4896 687.568
1.5 39.26 14.26 559.8476 804.2529
Figure 55: Analysis for stage storage curve
VOLUME
900
800
700
600
500
400
300
200
100
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6
Using a combination of orifice and weir to control the discharge at the outlet. We used 10
orifices with 1 m diameter, and 20 m length of 1 m depth weir.
2s/dt + Q
300
250
200
150
100
50
0
0 5 10 15 20 25 30 35 40 45
25.00
20.00
15.00
10.00
5.00
0.00
0 20 40 60 80 100 120 140 160
Time
The Routed hydrograph from the reservoir complied with the condition of reservoir routing. The
routed hydrograph cuts the recession limb of the post dev hydrograph and the peak discharge is
reduced as well as there is increase in time to achieve peak discharge.
4.8 Pre-dev Hydrograph & Post-dev routed Hydrograph:
Chart Title
45.00
40.00
35.00
30.00
25.00
Discharge
20.00
15.00
10.00
5.00
0.00
0
10
30
40
60
80
90
20
50
70
0
10
12
14
11
13
15
Time
The post-dev routed outflow hydrograph has reduced max Q, which satisfied the council
requirement for development of land. Which states that the Qp after development must be lower
than or equal to the pre-development hydrograph.
Pre Dev Post Dev
Outflow Routed
Hydrograph Hydrograph
TIME (m3/s) Qn (m3/s)
0 0.00 0.00
5 0.00 0.00
10 0.00 0.00
15 2.07 0.00
20 6.17 2.10
25 10.26 10.40
30 14.38 14.16
35 17.39 17.51
40 20.70 21.10
45 24.01 25.25
50 28.09 28.87
55 33.11 32.20
60 38.32 35.75
65 40.58 38.92
70 38.92 40.36
75 34.69 39.65
80 30.74 37.22
85 27.51 33.49
90 24.11 29.07
95 20.79 24.27
100 17.41 20.15
105 13.35 17.11
110 8.64 13.94
115 4.26 10.27
120 1.63 4.99
125 0.62 2.07
130 0.24 0.86
135 0.09 0.35
140 0.03 0.14
145 0.01 0.06
150 0.01 0.02
155 0.00 0.00
McCuen, R.H., Wong, S.L. and Rawls, W.J., 1984. Estimating urban time of concentration.
Journal of hydraulic Engineering, 110(7), pp.887-904.
Cudennec, C., Fouad, Y., Sumarjo Gatot, I. and Duchesne, J., 2004. A geomorphological
explanation of the unit hydrograph concept.
1961. Synthetic unit hydrographs for small watersheds. Journal of the Hydraulics Division,
87(4), pp.33-54.al Processes, 18(4), pp.603-621.
Dooge, J.C., Strupczewski, W.G. and Napiórkowski, J.J., 1982. Hydrodynamic derivation of
storage parameters of the Muskingum model. Journal of Hydrology, 54(4), pp.371-387.
Kožar, I., Lozzi-Kožar, D. and Jeričević, Ž., 2010. A note on the reservoir routing problem.
European Journal of Mechanics-B/Fluids, 29(6), pp.522-533.