Breaching Parameters for Earth and Rockfill Dams

Dam risk analysis is at the heart of dam failure prevention and mitigation. In order to assess dam risk, it is essential to conduct a quantitative analysis of the process of a dam breach, which can be described by such parameters as breach geometry, breaching duration, and peak outflow rate. The main objective of this paper is to develop robust empirical formulas with physical meaning for predicting dam breaching parameters based on past dam failure data. A database of 182 earth and rockfill dam failure cases has been compiled; among these cases nearly one-half are for large dams higher than 15 m. A multiparameter nonlinear regression model is recommended to develop empirical relationships between five breaching parameters (breach depth, breach top width, average breach width, peak outflow rate, and failure time) and five selected dam and reservoir control variables (dam height, reservoir shape coefficient, dam type, failure mode, and dam erodibility). The relative importance of each control variable is evaluated. The dam erodibility is found to be the most important factor, influencing all five breaching parameters. The reservoir shape coefficient and the failure mode also play an important role in the prediction models. Two case studies are presented to show the application of the empirical models developed in this paper.     

Embankment Dam Breach Parameters and Their Uncertainties

Potential flood hazards that would be created by breached embankment dams need to be evaluated to select spillway design floods and to prepare emergency action plans. The breaches are often modeled simply, usually in the shape of a trapezoid that is defined by its final height, base width or average width, and side slopes, along with the time needed for the opening to form completely. Data collected from 74 embankment dam failures were used to develop mathematical expressions for the expected values of the final width and side slope of a trapezoidal breach along with its formation time. Information is provided that allows variances of the predicted quantities to be calculated as well. The findings of the statistical analysis were then applied in a Monte Carlo simulation to estimate the degree of uncertainty of predicted peak flows and water levels downstream from breached embankment dams.     

Time for Development of Internal Erosion and Piping in Embankment Dams

A method is presented for the approximate estimation of the time for progression of internal erosion and piping, and development of a breach leading to failure in embankment dams and their foundations. The method accounts for the nature of the soils in the dam core, the foundation, and the materials in the downstream zone of the dam.?Guidance is also provided on the detectability of internal erosion and piping, taking account of the mechanism of initiation, continuation, and progression to form a breach, for internal erosion and piping in the embankment, the foundation and from the embankment to foundation.?It is shown that in many dams which have poor internal erosion and seepage control and are constructed mainly of earthfill, the time for potential development of piping is short, and for these dams continuous monitoring of seepage or surveillance would be needed to detect the piping in time to give warning of possible failure, and to give time to attempt intervention to prevent the failure.     

Model for Predicting Floods due to Earthen Dam Breaching. I: Formulation and Evaluation

This paper proposes a dam-breach model which predicts, in a simple but physically based manner, not only the peak discharge but the whole outflow hydrograph and breach development. The following aspects are taken into account: the geometry of the embankment, the shape of the reservoir, the hydraulic characteristics of the flow through the breach and its erosive capacity, and the shape of the breach. The model needs only one calibration parameter and can be easily applied to real cases. The application of the model with a single value of the calibration parameter produced excellent results in the simulation of 12 historic earthfill dam failures, with a discharge range covering three orders of magnitude.     

Overtopping Breaching of Noncohesive Homogeneous Embankments

Homogeneous small-amplitude embankments were constructed in flumes from a range of uniform noncohesive materials and breached by overtopping flows under constant reservoir level conditions. Embankment erosion evolves from primarily vertical to predominantly lateral in nature. The breach channel initially erodes the downstream face of the embankment with an invert slope parallel to the face, the breach invert slope then progressively flattening to a terminal value by rotating about a fixed pivot point along the base of the embankment, the location of this pivot point being a function of the size of the embankment material. The breach channel is of a curved (hourglass) shape in plan. Below the water line, breach cross-section width B variation with elevation y above the breach invert is nondimensionally described by B? = 2k?y?0.5, where for the breach cross section at the embankment crest k? = ?2.82[ln(Hb?)]+0.351, and Hb is the centerline breach crest elevation. Breach discharge Qb can be nondimensionally expressed as a function of the head hb on the breach-crest centerline and the breach crest length in plan Lb using Qb? = 0.242?Lb?(hb?)1.5. All expressions presented are applicable to full-width breach sections (double the half-breach section tested). The present findings enable prediction of the development with time of breach cross section, breach longitudinal profile, eroded volumes, and breach flows. The findings can be utilized for predictions of erosion and flooding occurring as the result of embankment failure, although in an engineering sense the quantitative findings of the present work await confirmation for larger embankments.     

Case Study of the Big Bay Dam Failure: Accuracy and Comparison of Breach Predictions

The Big Bay Dam embankment failure occurred on March 12, 2004, releasing 17,500,000?m3 (14,200?acre-ft) of water. In all, 104 structures were documented as being damaged or destroyed as a result of this failure. No human lives were lost. This paper documents data gathered and analyses performed on the hydraulics of the failure. High water levels from the failure were marked and measured. A HEC-RAS unsteady flow model was developed. Using observed breach geometry, HEC-RAS provided results that agreed with the measured high water marks from ?0.02?to??0.90?m and 0.01?to?0.62?m with associated modeled flow depths ranging from 9.3?to?5.7?m (from 30?to?19?ft). A peak breach flow of 4,160?m3/s (147,000?ft3/s) was predicted at the embankment. Breach peak flow prediction equations were found to substantially underpredict the peak flow indicated by HEC-RAS for this failure. HEC-RAS modeling utilizing predicted breach geometry and formation time also underpredicted the peak flow, but by a lesser amount. The National Resources Conservation Service models WinTR-20 and TR 66 were also assessed. WinTR-20 results compared reasonably well with the high water marks for this failure. TR-66 results did not compare well.     

Model for Predicting Floods due to Earthen Dam Breaching. II: Comparison with Other Methods and Predictive Use

A sensitivity analysis of the model presented in the companion paper is made and, on the basis of the results a criterion is proposed for the choice of values to assign to the side slopes of the breach, in order to use the model for prediction. Moreover, the sensitivity analysis shows that the outflow hydrograph and its peak value depend not only on the dam height and the stored volume, but also on the vertical distribution of the water mass in the reservoir. The model is compared with some previously published methods and the disadvantages, limitations, and errors that can be made using parametric models and predictive equations are pointed out. Finally, easy to use equations interpolating the numerical results of the model are provided that predict not only the peak discharge but the whole outflow hydrograph for overtopping failures.     

Numerical Modeling of Breach Erosion of River Embankments

The process of breach erosion of river embankments depends on the interaction among flow, sediment transport, and the corresponding morphological changes. Levees often consist of noncohesive material with a wide range of grain sizes. The dam material is mainly eroded due to the transport capacity of the overtopping water. Both bed load and suspended load are of importance. For breach formation, the lateral erosion due to slope instabilities has a significant impact. A depth averaged, two-dimensional numerical model was developed to account for these processes. The sensitivity of the discharge through the breach related to different processes and material parameters was investigated and compared to experimental and field data. The results show that the most sensitive parameter of an erosion-based dike-breach simulation is the breach side-slope angle which determines the lateral erosion. The application of the described Model 2dMb to different embankment failures at the Elbe River illustrates its capability in simulating overtopping breaching.     

某尾矿库筑坝方案比较分析

结合某山谷型尾矿库工程,在满足现有政策的前提下,从形成库容及服务年限、经济性、安全性、环境影响和运行管理等方面,对中线式尾矿筑坝法和上游式尾矿筑坝法两种筑坝方案进行综合比较和分析,并得出了该尾矿库工程采用中线式尾矿筑坝法方案为最优选择的结论.所采用的分析方法可为类似工程提供思路.     

Assessment and Rehabilitation of Embankment Dams

A series of observations, studies, and analyses to be made in the field and in the office are presented to gain a proper understanding of how an embankment dam fits into its geologic setting and how it interacts with the presence of the reservoir it impounds. It is intended to provide an introduction to the engineering challenges of assessment and rehabilitation of embankments, with particular reference to a Croton Dam embankment.     

Risk Analysis for Dam Overtopping—Feitsui Reservoir as a Case Study

Risk and uncertainty analysis by mathematical and statistical methods is often used to assess systematic risks and uncertainties. This research presents the procedure and application of risk and reliability analysis to dam overtopping. Annual maximum series of peak discharges of Feitsui Reservoir in northern Taiwan are used to analyze five extreme flood events with different frequencies. The highest water levels of the five extreme flood events were computed by using reservoir routing and considering seven factors subject to uncertainty. Afterward, the overtopping risk of Feitsui Dam was assessed by five uncertainty analysis methods: Rosenblueth’s point estimation method (RPEM), Harr’s point estimation method (HPEM), Monte Carlo simulation, Latin hypercube sampling, and the mean-value first-order second-moment (MFOSM) method. The results show that values of overtopping risk computed by different methods are similar. One may apply some approximated methods (MFOSM, HPEM and RPEM) to avoid the computational burden by applying sampling methods. Furthermore, the accuracy of results by approximated methods compared with that by sampling methods may differ from case to case. The selection and application of the uncertainty methods depend upon the information availability of the model parameters and model complexity. One may need to examine the model parameters and model complexity before determining appropriate methods to be used in a study.     

Discharge and Suspended Sediment Transport during Deconstruction of a Low-Head Dam

In March of 2003, the 43?m wide, 2.2?m high St. Johns Dam (Sandusky River, Ohio) was breached to lower the water level in the reservoir. In November of the same year, the dam was removed in an effort to restore aquatic habitat and connectivity in the river. During both the breach and the dam removal, high resolution time series of discharge and suspended sediment concentrations were monitored 200?m downstream of the dam. Discharge and suspended sediment during the breach were not discernible from background values. In contrast, the dam removal resulted in a peak suspended sediment concentration of 59?mg/L and a peak discharge of 33.5?m3/s, which returned to background levels of 19?mg/L and 1.5?m3/s, respectively, approximately 8?h after the removal. The floodwave during the removal attenuated by 50% at the City of Fremont, 53?km downstream, illustrating the diffusive nature of the channel and the limited risk of flooding downstream. Levels of suspended sediment and discharge during the removal were comparable to subsequent discharge events. Spatial distributions of turbidity in and upstream of the dam pool and archived turbidity data from the City of Tiffin, 13?km downstream of the dam, suggest that sediments stored in the impoundment did not statistically enhance turbidity up to 2 years after the removal. Generally, the removal had a minor impact on water quality and posed no risk to public safety or to downstream aquatic habitats.     

Embankment Dam on Liquefiable Foundation—Dynamic Behavior and Densification Remediation

Earthquake-induced liquefaction is a major concern for embankment dam safety. Many liquefaction-induced earth embankment failures or near failures have been reported around the world during various earthquakes. Such embankment damages were particularly destructive when the underlying saturated granular soils liquefied, resulting in cracking, settlement, lateral spreading, and slumping of the embankment. Through a series of four highly instrumented geotechnical centrifuge model tests, seismic behavior of a zoned embankment dam with saturated sandy soil foundation was studied under moderate earthquake conditions. The beneficial effects of foundation densification were investigated. Valuable insights into the dynamic behavior of the employed embankment–foundation systems are provided. Test results suggest that there may be an optimum depth of densification treatment beneath an earth dam beyond which the reduction of the earthquake-induced deformations are relatively minor and that relatively small and isolated zones (e.g., at depth) of loose material within a densified volume of soil may not impair the overall effectiveness of treatment and do not necessarily result in damaging displacements.     

Modeling Landslide Dambreak Flood Magnitudes: Case Study

Landslide dams typically comprise unconsolidated and poorly sorted material and are vulnerable to rapid failure and breaching, resulting in significant and sudden flood risk downstream. Hence they constitute a serious natural hazard, and rapid assessment of the likely peak flow rate is required to enable preparation of adequate mitigation strategies. To determine the relative utility and accuracy of dambreak flood forecasts, field estimates of peak outflow rates from the failure of the Poerua landslide dam in October 1999 were compared with estimates from physical laboratory modeling, empirical methods, and computer modeling. There was reasonable agreement among the field estimates, laboratory modeling, and computer modeling. Some empirical estimates were less reliable. Reasonably reliable estimates of peak outflow can be obtained from computer model routines sufficiently rapidly to be of use in an emergency management situation. The laboratory modeling demonstrated the effect of dam batter slopes and valley bed slope on peak outflow; this information could be used to refine empirical or numerical estimates of peak outflow.     

某尾矿库加高扩容坝体加固处理研究

某尾矿库加高扩容时,对土堤抛石挤淤,土堤周围出现隆起现象。为避免今后尾砂堆高过程中出现此类情况,保证堆积坝及浆砌石坝稳定安全,对堆积坝和浆砌石坝后尾矿区域进行水平加筋,以利尾矿排水固结,并与加筋形成复合坝基。经水平加筋后,效果良好,可供类似工程参考。     

Model Uncertainty in Finite-Element Analysis: Bayesian Finite Elements

In this paper, probabilistic models for structural analysis are put forward, with particular emphasis on model uncertainty. Context is provided by the finite-element method and the need for probabilistic prediction of structural performance in contemporary engineering. Sources of model uncertainty are identified and modeled. A Bayesian approach is suggested for the assessment of new model parameters within the element formulations. The expressions are formulated by means of numerical “sensors” that influence the model uncertainty, such as element distortion and degree of nonlinearity. An assessment procedure is proposed to identify the sensors that are most suitable to capture model uncertainty. This paper presents the general methodology and specific implementations for a general-purpose structural element. Two numerical examples are presented to demonstrate the methodology and its implications for probabilistic prediction of structural response.     

Fully Coupled Analysis of Failure and Remediation of Lower San Fernando Dam

The extent of flow deformation in an embankment dam is determined by the driving forces and the residual strength of the soil, as well as by the kinematic constraints. The soil conditions of berm and buttress, as well as of foundation, are also critical factors affecting seismic performance of an embankment dam. A careful examination of these factors is necessary when proposing remedial measures to a seismically deficient dam. This paper presents a set of fully coupled finite element analyses of the response of the well-known lower San Fernando Dam during the 1971 earthquake. A critical state model incorporating the concept of state-dependent dilatancy was employed to describe soil behavior over the full range of loading conditions encountered. The results show clearly that a flow slide occurred on the upstream side, and indicate that a downstream flow slide would occur, too, if the downstream berm had not been constructed before the event. The analyses show also that the addition of an upstream berm could effectively prevent the upstream flow slide.     

Two-Phase Flow Characteristics of Stepped Spillways

An experimental study on a large model flume with fiber-optical instrumentation indicated that minimum Reynolds and Weber numbers of about 105 and 100, respectively, are required for viscosity and surface tension effects to become negligible compared to gravitational and inertial forces expressed by Froude similitude. Both the location of and the flow depth at the inception point of air entrainment can be expressed as functions of a so-called roughness Froude number containing the unit discharge, step height and chute angle. The depth-averaged air concentration is found to depend only on a normalized vertical distance from the spillway crest and the chute angle for chute slopes ranging from embankment to gravity dam spillways. Air concentration profiles can be expressed by an air bubble diffusion model. The pseudobottom air concentration allows the assessment of the cavitation risk of stepped chutes and is approximated by a regression function. Finally, a new velocity distribution function is presented consisting of a power law up to 80% of the characteristic nondimensional mixture depth, and a constant value above.     

热轧钢板桩的发展和应用前景

钢板桩是一种新型环保建筑钢材,广泛应用于码头、堤防护岸、挡土墙、船坞、围堰等工程。钢板桩具有高强度、轻型、隔水性好、耐久性强、环保效果显著,施工简单,较强的救灾抢险功能等优点。分析了我国热轧钢板的发展,工程应用、标准以及市场前景。     

Characterization of Spectral Analysis of Surface Waves Shear Wave Velocity Measurement Uncertainty

The spectral analysis of surface waves (SASW) method is a nondestructive test for characterization of the variation with depth of the shear modulus of soils. While the testing procedure is well developed, only one preliminary study has investigated measurement uncertainty associated with SASW, and the methods utilized to quantify measurement uncertainty were prohibitive to routine assessment. Knowledge of this uncertainty, and ability to include its assessment in routine testing, would allow for inclusion of SASW results in reliability-based design and in assessment of the spatial variability of shear modulus. In this study, a large sample of test data was collected from two test sites. Characteristic statistics, statistical distribution, and measurement uncertainty were determined for each phase of SASW. Using the empirical statistical properties and measurement uncertainty results as validation criteria, an analytically based uncertainty assessment system was developed. Phase angle, inverse phase angle, and phase velocity data typically display a coefficient of variation (COV) of 2%, and the COV for combined phase velocity data is typically 1.5%. The COV for shear wave velocity is typically between 5 and 10%, and thus the inversion appears to magnify measurement uncertainty. Phase angle, inverse phase angle, phase velocity, and combined phase velocity data are normally distributed. Shear wave velocity samples at a given depth are generally normally distributed. Using a small sample of experimental data and the analytically based process developed in this study, the measurement uncertainty of SASW test results can be assessed as part of routine testing.