National Institute of Hydrology, Patna, Bihar India. National Institute of Hydrology, Roorkee, Uttarakhand India
National Institute of Hydrology, Patna, Bihar India. National Institute of Hydrology, Roorkee, Uttarakhand India
National Institute of Hydrology, Patna, Bihar India. National Institute of Hydrology, Roorkee, Uttarakhand India
Dam Break Flood Hazard Assessment: A Case Study for a Small Dam at Source Stream
of River Ganga in Uttarakhand, India
Pankaj Mani1, Rakesh Kumar2, J. P. Patra2
1
National Institute of Hydrology, Patna, Bihar India.
2
National Institute of Hydrology, Roorkee, Uttarakhand India.
Corresponding author email id: [email protected]
Abstract: Structural measures of a flood defense systems are designed for a specified degree of
protection that are defined through the exceedance probability of the flood. However, the failure of
major storage structure induces additional flooding risk in the downstream floodplains over and above
the specified degree of protection. Hence, the regulating agencies maintains the safety guidelines to
include the dam failure while assessing the comprehensive flood risk in the downstream reach. Thus,
the dam break analysis evaluates the safety level in the downstream river reach due to flood likely to
be induced from storage after the failure of dam with or without the flood generated from
meteorological conditions. The comprehensive dam break analysis includes the estimation of the
extent and magnitude of flooding, its occurrence time and duration of flooding. The specific flood
mitigation measures include preparation of emergency action plans for downstream river reach with
computation of available warning time and evacuation plans. The emergency action plans are known
to planners, local administrators and even likely to be affected population in priory. The public risk
perception helps to develop flood preparedness plan and in efficient risk management strategy. The
dam break analysis is essentially a two step procedures, (i) modelling the development of breach in
the dam section and computing the breach (flood) outflow, and (ii) routing of flood hydrograph in the
downstream river reach to compute the various flood attributes. In this paper, dam break analysis of
Dhauliganga dam located in Uttarakhand in Kumaon Lesser Himalayas has been reported. The
various scenarios of flooding due to dam breach of the concrete face rockfill dam has been carried out
and flood inundation, time of occurrences etc. has been estimated. The hydraulic model for the river
reach of about 30 km is developed in MIKE 11 using surveyed river cross sections and 10 m
resolution DEM of the study area generated using CARTODEM provided by NRSC, Hyderabad.
Three cases of flooding are simulated; (i) flooding due to PMF in the river causing dam break
condition; (ii) flooding due to PMF without dam break; and (iii) sunny day failure condition (dam
failure with nominal inflow when reservoir is full). It is observed that the time of travel of peak flood
from dam site to major settlement area at Dharchula, about 20 km downstream location is 42 minutes
for critical case of dam failure. The maximum flood level and time of travel of the peak flood at other
important locations are estimated. However, the analysis shows that no settlements/ villages area are
under inundation even for most critical flooding conditions. The extent of flood hazard for various
cases of flooding are estimated by superimposing the inundation map over Google Earth for detailed
description of inundated areas and the infrastructures likely to be affected.
Keywords: Dam break analysis, MIKE 11, Flood inundation, Flood hazard, EAP
1. Introduction
Protection of the public lives and properties from the consequences of dam failures are
important as massive population and infrastructure are vulnerable to dam break disasters. The
prior assessment of extent, magnitude and time of flooding due to failure of dam/
embankment are the important input in planning for flood preparedness measures. It is quite
difficult to conduct analysis and determine the warning time and extent of inundation at the
time of disaster. Therefore, pre-determination of these parameters is done by simulating a
hypothetical dam break/ embankment failure situation. The dam break flood is different from
the storm generated flood in the sense that it hardly gives any response time for emergency
action. The sudden and uncontrolled release of stored water, generally coinciding with the
catastrophic climatological events raises the flood magnitude to a very high level and thus the
water spreads to a wider area in the downstream causing the losses much more as compared
to normal rainfall generated floods. Further, with the increase in extremist actions over the
objects of national and economic interest, the emergency action plans are prerequisite. Such
plans consist of the maps showing the flood inundations under various scenarios of dam
failures. This helps in estimating the potential damage. The amount and extent of flooding
and occurrence of their time are also included. The threatened settlement, property, and other
infrastructures, in addition to emergency rescue plan, communication links etc. are the
important constituents of the EAP. The hypothetical dam break analysis is carried out to
predict potential flood damage and to prepare emergency action plan (EAP) in advance. Such
analysis is mandatory not only for the existing and old dams but also for the proposed dams.
(MoEF, 2015)
A complete dam break analysis involves a balanced consideration of hydrological, hydraulic,
environmental, and geotechnical and structural parameters pertaining to dam and downstream
flood plains. The description of time dependent flood wave propagation, downstream of a
breached dam is extremely complex. It is a function of site specific parameters including
reservoir characteristics and breaching characteristics of dam. Further, the movement of flood
wave across the flood plain will be governed by another range of determinants, many are
difficult to replicate in a mathematical model; e.g. dynamic variation in terrain and surfaces,
influence of land use change etc. The dam break flood simulations are carried out by
developing various scenarios leading to dam failure in addition to re-evaluating the probable
maximum flood (PMF), the most general cause of dam failure in natural circumstances. The
insufficient capacity of the spillways may be incapable to pass the updated estimate of PMF
causing dam failure due to overtop. Another likely reasons for dam failure includes
earthquake induced structural failure, equipment failure, criminal action, sabotage etc. (CWC,
2006, EAP)
With the advancement in computational fluid dynamics (CFD) and availability of new tools
and techniques like; geographical information system (GIS) and satellite images, some of the
difficulties have been successfully accounted for. Several experimental, analytical, and
numerical models have been developed to carry out dam break analysis. Computer-aided
numerical models such as DAMBRK (Fread, 1988), SMPDBK (Wetmore and Fread, 1991),
CADAMBRK (Liong et al., 1991), NWS FLDWAV (Fread, 1993), HEC RAS (USACE,
2006), BOSS DAMBRK (Kho, 2009), and MIKE 11 (DHI, 2004) have been widely used
successfully across the world due to their high computational speed and efficiency. Majority
of dam break analysis studies have been carried out using 1D model. The NWS DAMBRK
model has been used for dam break analysis of Barna dam, Madhya Pradesh, India (NIH,
1997), Ghodahoda project Odisha India (NIH, 2000), and the proposed dam on Yamuna
river, India (Lodhi and Agrawal, 2012). BOSS DAMBRK model was used to study the dam
break analysis of proposed dam on the Gerugu river Malaysia (Kho et al., 2009). HEC RAS
model has been used for dam break analysis of Oros Dam Brazil (Gee, 2008), Danjiangkou
and Yahekou dam failures in the Han river China (Minglong and Jayawardena, 2008), and
Foster Joseph Sayers Dam in Center Coutry PA, USA (Xiong, 2011). SMPDBK model has
been used for dam break analysis of Foster Joseph Sayers Dam in Center Coutry PA, USA
(Shahraki et al., 2012). MIKE 11 model has been used for dam break analysis of Bichom and
Tenga dam (Husain and Rai, 2000), Buffalo Creek Dam, North Carolina, USA (Tingsanchali
and Chinnarasri, 2001), Indra Sagar and Omkareshwar project, India (Pillai et al., 2012),
Hirakud dam, India (Mohite et al., 2014). 1D models, though simple to use and provide
information on bulk flow characteristics, fail to provide detailed information regarding the
flow field. Hence, attempts have been made to model the 2D nature of floodplain flow. A 1D
approach is used to describe breach growth and breach flow and a 2D approach is used to
predict flood propagation in the inundated areas. In this paper, a case study of failure of
Dhauliganga dam located in Uttarakhand in Kumaon Lesser Himalayas has been reported.
The project is maintained by NHPC Limited and has provided all the relevant data required
for the study. The study envisages the identification of various scenarios of flooding,
estimation of breach parameters, modelling of dam break flood and its routing in the
downstream reach to compute the maximum flood inundation and its time of occurrences.
The failure of Dhauliganga dam has been simulated in MIKE-11 model. The movement of
flood in the downstream reach is solved through Saint-Venant equations using 6-points
Abbott scheme in MIKE 11 (DHI, 2004). The breaches is modeled by means of a “dam
break” structure and its development is described by time series for breach width, crest level
and side slope. The maximum flood level and its time of occurrence and flood warning time
have been estimated at important locations. The extent of maximum inundation for various
cases of flooding has also been computed.
India
Nepal
3. Methodology
Several commercial software like NWS DAMBRK, HECRAS, and MIKE 11 models are
available for carrying out the dam break modeling. All these software essentially carry out
hydrodynamic modelling involving solution of the basic St. Venant’s partial differential
equations of conservation of mass and momentum as shown below:
i. Conservation of mass (continuity) equation
(∂Q/∂X) + ∂(A + A0) / ∂t - q = 0 (1)
where, Q = discharge; A = active flow area; A0 = inactive storage area; h = water surface
elevation; q= lateral outflow; x = distance along waterway; t = time; S f = friction slope; Sc =
expansion contraction slope and g = gravitational acceleration.
The mathematical modeling of dam break flow is a cost and time effective approach to
(approximately) solve the above governing equations. The present dam break study has been
carried out using MIKE 11 HD model developed by Danish Hydraulic Institute (DHI). The
core of the MIKE 11 system consists of the HD (hydrodynamic) module, which is capable of
simulating unsteady flows in a network of open channels. The results of a HD simulation
consist of time series of water levels and discharges. MIKE 11 hydrodynamic module is an
implicit, finite difference model for unsteady flow computation. The model can describe sub-
critical as well as supercritical flow conditions through a numerical description which is
altered according to the local flow conditions (in time and space). The Dam Break Module in
MIKE 11 simulates the outflow hydrograph resulting from the failure of a dam.
Breach Parameter Estimation - The parameters of the breach are given as a time series in
the boundary editor. The time is relative to the start of the breach. Between the specified
times the parameters are linearly interpolated. The NWS DAMBRK method simulates the
failure by two methods; Breach failure and Piping failure. The breach failure uses a weir type
equation to determine the flow through the breach where as the latter is based on an orifice
type equation. Both methods rely on the user specifying the full development of the breach
geometry as a function of time. The estimation of the breach location, size, and development
time are crucial in order to make an accurate estimate of the outflow hydrographs and
downstream inundation. Many case studies have been performed on data from historic dam
failures, leading to guidelines, regression equations, and computer modeling methodologies for
prediction of the dam breach size and time. Tony Wahl (1988) summarized a comprehensive
literature on historic dam failures data and recommends range of breach parameter for dam break
analysis. The Federal Agencies guidelines and regression equation based approach are widely
used. In fact the guideline provides the upper and lower bound of parameters values while the
range estimated from regression equations are used for parameter sensitivity analysis (TD-39,
2014). The estimated breach parameters for the dam is shown in Table-1.
Table-1: Description of breach parameters
For dam break analysis, the most Figure-2: Inflow design flood and outflow
critical flooding scenario is hydrograph at dam section.
considered when reservoir is at FRL,
i.e. at RL 1345.0 m and peak of PMF
flood enters into reservoir
considering all spillway gates to be 4421.47
image to prepare the inundation maps. The inundation map for DB flood near Dharchula is
shown in Figure-5.
Figure-5: Synoptic view of flood inundation over Google Earth for dam break flood.
downstream reach of the Dhauliganga dam are Dhauliganga power house and Dharchula
Township, at which the maximum flood elevation (msl) has been computed for various cases
of flooding. For critical DBF, the maximum flood elevation at Dhauliganga power house and
Dharchula Township are RL 1032.0 m and RL 857.3 m respectively. Similar type of flood
simulation have also been carried out for various combination of inflow design flood and dam
break condition as shown in Table-2. The peak flood level and its occurrence time are
estimated at important locations are also shown in this table. The time of occurrence of
maximum flood level is shown within bracket. The maximum flood level occurs at Power
house is RL 1031.63 m while the elevation of service bay and machine hall at Dhauliganga
power house is RL 1039.5 m. Hence, the important plant installations are above the
maximum flood level. Further, the maximum flood level is computed at Dharchula is below
the danger level, the slope stability and river bank erosion are the major issues in the region
due to very steep slope and fragile geological formation (Paudel et al., 2013). The time of
travel of peak flood from dam site to a Dhauliganga power house is 38 minute for Case-1.
Under pre project scenario, the maximum flood level due to critical condition of flooding
(Case-2) is RL 1031.36 m. In this case, the time of travel of peak flood is 4 minutes due to
steep river slope. As, the travel time is very small, an effective public warning system is
proposed to prevent any casualty.
Table-2: Maximum flood level at important locations for various cases of flooding.
Maximum Important locations
Particulars
inundation area Power house Dharchula
MIKE 11 Chainage (km) (ha) 9.75 20
Case-1 1031.63 m 857.06 m
163.27
(PMF flood and dam break) 38 min 42 min
Case-2 1031.36 m 856.87 m
156.39
(PMF flood only) 40 min 45 min
Case-3 1026.95 m 853.84 m
94.53
(only dam break) 1 hr 6 min 1 hr 12 min
Case-4 1028.90 m 855.19 m
122.64
(100 years return period flood) 56 min 1 hr 5 min
(Time in bracket shows the time of travel of peak flood from dam site to important locations)
Acknowledgment
This paper is originated from the original study sponsored by NHPC Limited. The authors
duly acknowledges the supports provided by the project authority in carrying out the study.
References
CWC EPA 2006. Guidelines for Development and Implementation of Emergency Action
Plan (EAP) for Dams, Central Water Commission, New Delhi.
DHI 2004 MIKE 11 User Manual. Danish Hydraulic Institute: Horsholm, Denmark.
Fread DL 1988. The NWS DAMBRK model-Theoretical background/ user documentation.
Hydrology Research Laboratory, Office of Hydrology, National Weather Service
(NWS), NOAA, Silver Spring, Maryland.
Fread DL 1993. NWS FLDWAV Model-The replacement of DAMBRK for dam break flood
prediction. Hydrology Research Laboratory, Office of Hydrology, National Weather
Service (NWS), NOAA, Silver Spring, Maryland.
Freed, D. L. and Lewis, J. M. 1998. NWS FLDWAV Model: Theoretical description and
User Documentation. National Weather Services, NOAA.
Gee, D., and Brunner, G. W. 2005. Dam Break Flood Routing Using HEC-RAS and NWS-
FLDWAV, Proceedings of World Water and Environmental Resources Congress
2005, Impacts of Global Climate Change.
Husain SM, Rai NN 2000. One dimensional dam break flood analysis for Kameng hydro
electric project, India. International seminar and workshop, Seinäjoki, Finland, 2-5
October.
Kho FWL, Law PL, Lai SH, Oon YW, Ngu LH, Ting HS 2009. Quantitative dam break
analysis on a reservoir earth dam. International Journal of Environmental Science and
Technology, 6 (2): 203-210.
Liong SY, Tan CC, Alwis WAM, Selvalingam S 1991. Computer-aided dam break flow
routing. Advances in Engineering Software and Workstations, 13(3):110-115.
Lodhi MS, Agrawal DK 2012. Dam-break flood simulation under various likely scenarios
and mapping using GIS: Case of a proposed dam on River Yamuna, India. Journal of
Mountain Sciences, (9): 214-220.
Minglong D, Jayawardena AW 2008. Dam-break flood analysis in mid-down stream of Han
river. http://www.icharm.pwri.go.jp/publication/pdf/2008/4119_master_course2007
/annex/11_1_daiminglong_synoposis.pdf.
MoEF, 2015. Guidance Manual for Environmental Impact Assessment and Clearance of
River Valley Projects, http://environmentclearance.nic.in/writereaddata/Form-
1A/HomeLinks/Model.htm, last accessed on 07 July 2015.
Mohite AR, Mani P, Chatterjee C and Singh R, 2014. Dam Break Analysis of the Hirakud
Dam using MIKE 11 model, National Conference on Emerging Technology Trends in
Agricultural Engineering, November 07–09, NERIST, Nirjuli (Itanagar), India.
NIH 1997. Dam break study of Barna Dam. Technical report, CS-AR-20/ 96-97. National
Institute of Hydrology, Jal Vigyan Bhawan, Roorkee, India.
NIH 2000. Dam break study of Gogahoda project, Orissa. Technical report, CS-AR-15/ 96-
97. National Institute of Hydrology, Jal Vigyan Bhawan, Roorkee, India.
Pillai BRK, Ghosh S, Gupta OP (2012) Integrated dam break analysis – Case study of Indira
Sagar and Omkareshwar projects. India Water Week 2012 – Water, Energy and Food
Security: Call for Solutions, 10-14 April 2012, New Delhi.
Shahraki A, Zadbar A, Motevalli M, Aghajani F 2012. Modeling of earth dam break with
SMPDBK-Case Study Bidekan Earth Dam. World Applied Sciences Journal, 19 (3):
376-386.
Tingsanchali T and Chinnarasri C 2001. Numerical modelling of dam failure due to flow
overtopping. Hydrological Sciences Journal, 46(1): 113-130.
USACE 2006. HEC-RAS River Analysis System, User’s Manual Version 4.0 Beta,
Hydrologic Engineering Center, US Army Corps of Engineers.
Wetmore JN, Fread DL 1991.The NWS simplified dam-break flood forecasting model.
Hydrology Research Laboratory, Office of Hydrology, National Weather Service
(NWS), NOAA, Silver Spring, Maryland.
Xiong YA 2011. Dam break analysis using HEC-RAS. Journal of Water Resource and
Protection, 3:370-379.