CASE STUDIES OF DAM-BREACH GEOMETRY AND HYDRAULICS

Victor M. Ponce

Professor, Department of Civil and Environmental Engineering,
San Diego State University, San Diego, California 92182-1324

Ahmad Taher-Shamsi

Associate Professor, Amirkabir University of Technology, Tehran, Iran;
Currently on sabbatical leave at the Department of Civil and Environmental Engineering,
San Diego State University, San Diego, California 92182-1324

Ampar V. Shetty

Research Associate, Department of Civil and Environmental Engineering,
San Diego State University, San Diego, California 92182-1324


ABSTRACT

A review of dam-breach parameters and prediction methods developed from the analysis of historic embankment-dam failures is performed. A new shape factor representing dam-breach geometry is formulated. A relationship between Froude number based on peak discharge and the new shape factor is identified.


INTRODUCTION

Dams provide many benefits to society; however, floods resulting from dam failures have produced some of the most devastating disasters of the last two centuries (Wahl, 1998). Therefore, earth-dam failure has become a subject of increasing concern among dam engineers, federal, state, and local officials, and society at large.

Despite the existence of ample documentation, the prevailing mechanisms for earth-dam failure are still not clearly understood. Studies have shown that earth-dam failures can be due to many causes. Generally, these are classified into: (1) hydrologic, and (2) geotechnical. Failure usually results in the eventual development of a breach compromising a certain length of embankment. The breach width has been documented in a large number of cases. Dam-breach studies have shown that the shape and time evolution of the breach determine to a large extent the characteristics of the outflow hydrograph (Ponce, 1982).

The objective of this paper is to develop a dimensionless relationship between geometric and hydraulic parameters governing earth-dam breach failures. For this purpose, the following geometric and hydraulic parameters are defined:

  • Breach maximum depth Db
  • Breach top width Wb
  • Breach peak outflow Qp
  • Breach time-to-failure tf
  • Dam maximum height Hd
  • Dam top width Wd
  • Reservoir volume Vr
  • Reservoir characteristic length Lr

CASE STUDIES OF DAM-BREACH FAILURES

Valuable information is available from embankment failures that have been documented. Wahl (1998) identified many breach-prediction equations developed since 1984. These equations have been based on analyses of case studies comprising from 20 to 60 failures. Sources of case-study data for breached dams are numerous. However, many of the failures have taken place before the need was recognized to fully document the breach process and associated breach characteristics.

Babb and Mermel (1968) summarized over 600 dam incidents throughout the world; however, high quality, detailed information was lacking in most cases. During the 1980's, several researchers compiled databases of well documented case studies in efforts to develop predictive relations for breach peak outflows (SCS, 1981; Ponce, 1982; Singh and Snorrason, 1982; MacDonald and Langridge-Monopolis, 1984; Costa, 1985; Froehlich, 1987, 1995a, 1995b; Singh and Scarlatos, 1988).

Wahl (1998) performed a literature search to produce a single database containing all case studies cited, comprising a total of 108 embankment failures. The type, amount and quality of data available for individual case studies varied widely. There are many instances of significant discrepancies between similar data reported by different investigators.

Researchers have proposed many relations for estimating breach parameters and peak outflows. The effect of breach parameters on peak discharge and evacuation time has been reported in the literature. The peak discharge was examined by Singh and Sorranson (1984). Petrascheck and Sydler (1984) demonstrated the sensitivity of peak outflow, inundation level, and flood-arrival time to changes in breach width and breach-formation time. The shape and time evolution of the breach, along with the size of dam and reservoir, determine to a large extent the characteristics of the outflow hydrograph during an earth-dam breach (Ponce, 1982). Singh and Sorranson (1982) provided the first quantitative guidance on breach width; they plotted breach width vs dam height for twenty dam failures.

MacDonald and Langridge-Monopolis (1984) related the breach time-to-failure to the volume of eroded material. The volume of eroded material was related to the breach-formation factor, defined as the product of the outflow volume times the initial depth of water above the breach bottom.

Froehlich (1987) developed nondimensional prediction equations for estimating average breach width, average side-slope factor, and breach-formation time. Froehlich's predictions were based on the characteristics of the dam, including reservoir volume, depth of water above the breach bottom, breach depth, length of embankment at the dam crest, breach bottom width, and empirical coefficients that account for overtopping vs non-overtopping failures and the presence or absence of a core wall. Singh and Snorrason (1982; 1984) presented relations for peak outflow as a function of dam height and reservoir storage.

Froehlich (1995a) developed a regression equation for the prediction of peak outflow based on reservoir volume and hydraulic head. He used data from 22 case studies for which peak outflows were available. Harris and Wagner (1967) applied the Schoklitsch sediment-transport equation to dam-breach outflows. Lou (1981), and Ponce and Tsivoglou (1981) developed a model that linked the Meyer-Peter and Muller sediment-transport equation to the one-dimensional differential equations of unsteady water and sediment flow.

Walder and O'Connor (1997) presented a physically-based model of dam-breach formation, and used it to relate dimensionless peak outflow to a dimensionless parameter containing drop in reservoir level, volume of released water, and mean vertical erosion rate of the breach.

Ponce (1982) defined a Froude number based on breach peak outflow and related it to a shape factor defined as follows: S = (WbDb) / (WdHd). The result showed a definite trend to an inverse relationship. Ponce's results were very similar to those reported by Black (1925). However, both relations show a poor fit to the data, as depicted in Fig. 1 .

DATA   ANALYSIS

A comparison was made between Ponce's data and the Dam Safety Office database (Wahl, 1998). It was determined that only 25 cases reported by Wahl are useful for this study; ostensibly, these are the same as those reported by Ponce (1982). A review of this work indicates that it was collected directly from original sources. For this reason, it is preferred for use in this study.

The following peak-outflow Froude number was used in this study:

Eq. 1(Eq. 1)

in which g= gravitational acceleration.

A new shape factor Sf is defined as follows:

Eq. 2(Eq. 2)

An analysis of the relation between peak-outflow Froude number and the new shape factor Sf, based on all the available data, led to the graphs shown in Fig. 2 and Fig. 3. The relationship is:

Eq. 3(Eq. 3)

with a coefficient of determination r2= 0.75.

According to Fig. 3, there are two limiting trends:

  1. When Sf approaches 0,  F increases above 1.  For F = 1,  the shape factor Sf = 0.00269.  The ratio Vr /(Wd Hd) can be taken as the characteristic length of the reservoir Lr.  Therefore, when Lr is greater than or equal to 371Wb,   F is greater than or equal to 1.

  2. When Sf approaches 1,  F approaches 0.  In other words, when Lr equal to WbF approaches 0. In this condition, there is no flow that would create a dangerous flood. 

Given Eq. 3, the peak-outflow Froude number can be calculated as a function of the new shape factor.  For example, in the case of Teton Dam failure, with Vr= 355 million m3, Hd= 93 m, Wd= 945 m, and assuming Db= Hd, a sensitivity analyses is performed to determine peak outflows for a possible range of breach widths. The calculations are presented in dimensionless form in Fig. 4, in which Wba is the assumed value of breach width and Wbd is the documented value of breach width; Qpc is the calculated value of peak outflow and Qpd is the documented value of peak outflow. It is seen that peak outflow is not too sensitive to breach width, which enhances its prediction.

CONCLUSIONS

Based on the available data, there seems to be a good correlation between the peak-outflow Froude number F and the new shape factor Sf. In addition, the reservoir characteristic length Lr is an indicator of the quantity of water which is responsible for the magnitude of peak outflows during a dam breach.

REFERENCES

Babb, A.O. and T.W. Mermel. 1968. "Catalog of Dam Disaster, Failures and Accidents," Bureau of Reclamation, Washington, DC.

Black, E. B. 1925. "Partial Failure of Earth Dam at Horton,Kansas,"Engineering News Record, 95(2). 58-60.

Costa, J. E. 1985. "Flood from Dam Failures," U.S. Geological Survey Open-file Report 85-560, Denver, Colorado, 54 p.

Froehlich, D. C. 1987. "Embankment-Dam Breach Parameters," Hydraulic Engineering Proceedings of the 1987 ASCE National Conference on Hydraulic Engineering, Williamsburg, Virginia, August 3-7, 1987, 570-575.

Froehlich, D. C. 1995(a). "Peak Outflow from Breached Embankment Dam," Water Resources Engineering , Proceedings of the 1995 ASCE Conference on Water Resources Engineering, San Antonio, Texas, August 14-18, 1995, 887-891.

Froehlich, D. C. 1995(b). "Embankment Dam Breach Parameters Revisited," Journal of Water Resources Planning and Management, 121(1), 90-97.

Harris, G. W. and D. A. Wagner. 1967. "Outflow from Breached Earth Dams,"University of Utah,Salt Lake City, Utah.

Lou, W. C. 1981. "Mathematical Modelling of Earth Dam Breaches," Thesis,Presented to Colorado State University," at Fort Collins, Colorado, in partial fulfillment of the requirements for the Degree of Doctor of Philosophy.

MacDonald, T. C. and J. L. Monopolis. 1984, "Breaching Characteristics of Dam Failures," Journal of Hydraulic Engineering, 110 (5), 567-586.

Petrascheck, A. W. and P. A. Sydler. 1984. "Routing of Dam Break Flood," International Water Power and Dam Construction, 36,29-32.

Ponce, V. M. 1982. "Documented Cases of Earth Dam Breaches," San Diego State University, SDSU Civil Engineering Series, No. 82149.

Ponce, V. M. and A. J. Tsivoglou. 1981. "Modeling Gradual Dam Breaches", Journal of the Hydraulics Division, ASCE, 107(7), 829-838.

Singh, K. P. and A. Snorrason. 1982. "Sensitivity of Outflow Peaks and Flood Stage to the Selection of Dam Breach Parameters and Simulation Models," Journal of Hydrology, 68, 295-310.

Singh, K. P. and A. Snorrason, 1984, "Sensitivity of Outflow Peaks and Flood Stage to the Selection of Dam Breach Parameters and Simulation Models," SWS Contract Report 288, Illinois Department of Energy and Natural Resources, State Water Survey Division, Surface Water Section at the University of Illinois, 179 p.

Singh, K. P. and P. D. Scarlatos. 1985."Breach Erosion of Earthfill Dams and Flood Routing,"BEED Model, Research Report, Army Research Office, Battelle, Research Triangle Park, North Carolina, 131p.

Soil Conservation Service. 1981. "Simplified Dam-Breach Routing Procedure", Technical Release No.66 (Rev.1), December 1981, 39p.

Wadler, J. S. and J. E. O'Connor. 1997. "Methods for Predicting Peak Discharge of Floods Caused by Failure of Natural and Constructed Earth Dams," Water Resources Research , 33(10), October 1997, 12p.

Wahl, T. L. 1998. "Prediction of Embankment Dam Breach Parameters," DSO-98-004, Dam Safety Research Report, U.S. Department of the Interior Bureau of Reclamation, Dam Safety Office.


http://ponce.sdsu.edu/dambreachgeometry.html 020703