Fifteen-stage operation of gated spillways for flood routing management through artificial reservoirs

Abstract

A procedure to identify sets of operational rules for gated spillways for optimal flood routing management of artificial reservoirs is proposed. The flood retention storage of a dam having a gated flood spillway is divided into 15 sub-storages whose surface elevations are identified as critical levels. The most suitable operation set for the downstream conditions and for the dam can be chosen from many derived operation sets. The spillway gates are operated in an optimum way for any floods from very small magnitudes to the probable maximum flood (PMF), without having to forecast the actual magnitude of the incoming flood hydrograph. Decision floods are formed by dividing the PMF into 15 sub-hydrographs by 5 and 10% increments in the ranges 5–50% and 50–100% of the PMF, respectively. Many potential spillway gate openings from closed to fully open are chosen initially. As a result of a series of routing simulations of 15 decision floods, a set of 15 gate openings is determined such that all floods from very small magnitudes to the PMF may be routed without overtopping the dam crest. Next, a few more 15-stage operation rules are determined such that the gate openings of the initial stages are decreased as their critical levels are increased stepwise, with the objective of attenuating smaller floods more effectively and releasing higher outflows for larger floods close to and including the PMF. The developed model is applied to the Catalan and Aslantas dams in Turkey, both of which serve for flood mitigation as well as hydropower generation.

Citation Haktanir, T., Citakoglu, H., and Acanal, N., 2013. Fifteen-stage operation of gated spillways for flood routing management through artificial reservoirs. Hydrological Sciences Journal, 58 (5), 1013–1031.

Editor Z.W. Kundzewicz; Associate editor A. Montanari     

Dam breach simulations with multiple breach locations under wind and wave actions

Multiple breaches of a dam resulting from wind-generated waves and wave overtopping are studied for a hypothetical long non-cohesive earthen dam with an uneven crest. Both wind speed and direction affect breach locations and outflow for a particular reservoir surface geometry. Locations on the dam with longer fetches along the wind direction are more subject to wave overtopping and breaching than other locations. Higher wind speeds lead to wave overtopping and dam breaches under larger freeboards than lower wind speeds. For a specified inflow hydrograph and spillway configuration, there exists a location at which the smallest estimated peak outflow occurs among all possible breach locations and the pool drops too quickly for additional breaches to develop. Using this location for a fuse plug or a pilot channel could minimize downstream impact, perhaps as an interim or emergency measure for a dam with inadequate spillway capacity.     

Modeling parametric uncertainty in optimal open channel design using FORM-PGSL coupled approach

Optimal design of artificial open channels is essential for the planning and management of irrigation projects. In this paper a modified formulation is presented for the comprehensive design of open channels considering the seepage loss, evaporation loss and land acquisition cost along with the lining and excavation cost. The resulting formulation is solved using a recent meta-heuristic optimization technique namely probabilistic global search Lausanne (PGSL). The uncertainty associated with channel design parameter may lead to the failure of canals (channels). The parametric uncertainty in open channel design is modeled using first order reliability method (FORM). A bi-objective optimization model is presented in the study which minimizes the cost and minimizes the probability of overtopping considering a probabilistic cost function as the objective function. A new approach is proposed to solve the model in a meta-heuristic environment following PGSL as the solution method. Also a chance constrained optimization model which considers overtopping probability constraint and channel capacity constraint simultaneously along with the objective of minimization of cost is propounded and solved using PGSL. The solutions obtained using coupled FORM-PGSL approach is encouraging and the method can be used for optimal and reliable design of artificial open channels.     

Improvement of overtopping risk evaluations using probabilistic concepts for existing dams

Hydrologic risk analysis for dam safety relies on a series of probabilistic analyses of rainfall-runoff and flow routing models, and their associated inputs. This is a complex problem in that the probability distributions of multiple independent and derived random variables need to be estimated in order to evaluate the probability of dam overtopping. Typically, parametric density estimation methods have been applied in this setting, and the exhaustive Monte Carlo simulation (MCS) of models is used to derive some of the distributions. Often, the distributions used to model some of the random variables are inappropriate relative to the expected behaviour of these variables, and as a result, simulations of the system can lead to unrealistic values of extreme rainfall or water surface levels and hence of the probability of dam overtopping. In this paper, three major innovations are introduced to address this situation. The first is the use of nonparametric probability density estimation methods for selected variables, the second is the use of Latin Hypercube sampling to improve the efficiency of MCS driven by the multiple random variables, and the third is the use of Bootstrap resampling to determine initial water surface level. An application to the Soyang Dam in South Korea illustrates how the traditional parametric approach can lead to potentially unrealistic estimates of dam safety, while the proposed approach provides rather reasonable estimates and an assessment of their sensitivity to key parameters.     

Reliability analysis of hydraulic structures considering unit hydrograph uncertainty

Unit hydrographs (UHs), along with design rainfalls, are frequently used to determine the discharge hydrograph for design and evaluation of hydraulic structures. Due to the presence of various uncertainties in its derivation, the resulting UH is inevitably subject to uncertainty. Consequently, the performance of hydraulic structures under the design storm condition is uncertain. This paper integrates the linearly constrained Monte-Carlo simulation with the UH theory and routing techniques to evaluate the reliability of hydraulic structures. The linear constraint is considered because the water volume of each generated design direct runoff hydrograph should be equal to that of the design effective rainfall hyetograph or the water volume of each generated UH must be equal to one inch (or cm) over the watershed. For illustration, the proposed methodology is applied to evaluate the overtopping risk of a hypothetical flood detention reservoir downstream of Tong-Tou watershed in Taiwan.     

Reliability analysis of hydraulic structures considering unit hydrograph uncertainty

Unit hydrographs (UHs), along with design rainfalls, are frequently used to determine the discharge hydrograph for design and evaluation of hydraulic structures. Due to the presence of various uncertainties in its derivation, the resulting UH is inevitably subject to uncertainty. Consequently, the performance of hydraulic structures under the design storm condition is uncertain. This paper integrates the linearly constrained Monte-Carlo simulation with the UH theory and routing techniques to evaluate the reliability of hydraulic structures. The linear constraint is considered because the water volume of each generated design direct runoff hydrograph should be equal to that of the design effective rainfall hyetograph or the water volume of each generated UH must be equal to one inch (or cm) over the watershed. For illustration, the proposed methodology is applied to evaluate the overtopping risk of a hypothetical flood detention reservoir downstream of Tong-Tou watershed in Taiwan.     

A non-traditional approach to the analysis of flood hazard for dams

The traditional and still prevailing approach to characterization of flood hazards to dams is the inflow design flood (IDF). The IDF, defined either deterministically or probabilistically, is necessary for sizing a dam, its discharge facilities and reservoir storage. However, within the dam safety risk informed decision framework, the IDF does not carry much relevance, no matter how accurately it is characterized. In many cases, the probability of the reservoir inflow tells us little about the probability of dam overtopping. Typically, the reservoir inflow and its associated probability of occurrence is modified by the interplay of a number of factors (reservoir storage, reservoir operating rules and various operational faults and natural disturbances) on its way to becoming the reservoir outflow and corresponding peak level—the two parameters that represent hydrologic hazard acting upon the dam. To properly manage flood risk, it is essential to change approach to flood hazard analysis for dam safety from the currently prevailing focus on reservoir inflows and instead focus on reservoir outflows and corresponding reservoir levels. To demonstrate these points, this paper presents stochastic simulation of floods on a cascade system of three dams and shows progression from exceedance probabilities of reservoir inflow to exceedance probabilities of peak reservoir level depending on initial reservoir level, storage availability, reservoir operating rules and availability of discharge facilities on demand. The results show that the dam overtopping is more likely to be caused by a combination of a smaller flood and a system component failure than by an extreme flood on its own.     

Evaluation of dam overtopping probability induced by flood and wind

This study develops a probability-based methodology to evaluate dam overtopping probability that accounts for the uncertainties arising from wind speed and peak flood. A wind speed frequency model and flood frequency analysis, including various distribution types and uncertainties in their parameters, are presented. Furthermore, dam overtopping probabilities based on monthly maximum (MMax) series models are compared with those of the annual maximum (AMax) series models. An efficient sampling scheme, which is a combination of importance sampling (IS) and Latin Hypercube sampling (LHS) methods, is proposed to generate samples of peak flow rate and wind speed especially for rare events. Reservoir routing, which incorporates operation rules, wind setup, and run-up, is used to evaluate dam overtopping probability.     

Hydrological risk analysis of dam overtopping using bivariate statistical approach: a case study from Geheyan Reservoir,China

Hydrological risk analysis is essential and provides useful information for dam safety management and decision-making. This study presents the application of bivariate flood frequency analysis to risk analysis of dam overtopping for Geheyan Reservoir in China. The dependence between the flood peak and volume is modelled with the copula function. A Monte Carlo procedure is conducted to generate 100,000 random flood peak-volume pairs, which are subsequently transformed to corresponding design flood hydrographs (DFHs) by amplifying the selected annual maximum flood hydrographs (AMFHs). These synthetic DFHs are routed through the reservoir to obtain the frequency curve of maximum water level and assess the risk of dam overtopping. Sensitive analysis is performed to investigate the influence of different AMFH shapes and correlation coefficients of flood peak and volume on estimated overtopping risks. The results show that synthetic DFH with AMFH shape characterized by a delayed time to peak results in higher risk, and therefore highlight the importance of including a range of possible AMFH shapes in the dam risk analysis. It is also demonstrated that the overtopping risk is increased as the correlation coefficient of flood peak and volume increases and underestimated in the independence case (i.e. traditional univariate approach), while overestimated in the full dependence case. The bivariate statistical approach based on copulas can effectively capture the actual dependence between flood peak and volume, which should be preferred in the dam risk analysis practice.     

Seismic upgrade of hydraulic fill dam by buttressing

The vulnerability of hydraulic fill dams under strong earthquake shaking has long been recognized. When located in areas of high seismic hazard, seismic upgrading of these types of dams is often required to meet current dam safety standards. Selection of an appropriate design concept for seismic upgrading of such dams requires careful consideration of seismically induced deformations when the hydraulic fill is to remain as part of the dam. This paper presents a case history of the seismic upgrade of Butt Valley Dam, a hydraulic fill dam located in Northern California. The dam was strengthened to withstand the Maximum Credible Earthquake (MCE) by buttressing of its upstream and downstream slopes. The paper discusses the evaluation of alternatives to upgrade the dam, the design criteria, and the design and analysis of the seismic upgrade. It is shown that a conservative and robust design was developed based on well-established engineering principles and multiple lines of defense, and sound use of analysis procedures including finite-difference non-linear dynamic deformation analyses.     

Simulating bed evolution following the Barlin Dam (Taiwan,China) failure with implications for sediment dynamics modeling of dam removal

The failure of the Barlin Dam in Taiwan, China offers an important case study for evaluating concepts in modeling the rapid erosion and channel recovery following intentional and unplanned dam removals. We present a modeling effort that applied a 1D and quasi-2D uncoupled hydraulics and sediment model (NETSTARS) to evaluate how discretization and parameterization influence thefitofbed elevationpredic-tions to observations following dam failure. Our analysis evaluated the model sensitivity to sediment transport function, active layer thickness, and number of stream tubes used to define the cross-section. Results indicate that a) the model is more sensitive to active layer thickness and sediment transport function than to the number of stream tubes, b) development of dam removal models are likely to benefit from varying the active layer thickness in time, and c) increased lateral discretization does not appear to improve model fit in the steep and rapidly changing river environment at our site. We conclude with discussion on differences between, identifying the need for, and general use of 1D, quasi-2D, and fully 2D models in dam removal and failure analysis.     

Risk analysis for clustered check dams due to heavy rainfall

Check dams are commonly constructed around the world for alleviating soil erosion and preventing sedimentation of downstream rivers and reservoirs.Check dams are more vulnerable to failure due to their less stringent flood control standards compared to other dams.Determining the critical precipitation that will result in overtopping of a dam is a useful approach to assessing the risk of failure on a probabilistic basis and for providing early warning in case of an emergency.However,many check dams are built in groups,spreading in several tributaries in cascade forms,comprising a complex network.Determining the critical precipitation for dam overtopping requires a knowledge of its upstream dams on whether they survived or were overtopped during the same storm,while these upstream dams in turn need the information for their upstream dams.The current paper presents an approach of decomposing the dam cluster into(1)the heading dam,(2)border dams,and(3)intermediate dams.The algorithm begins with the border dams that have no upstream dams and proceeds with upgraded maps without the previous border dams until all the dams have been checked.It is believed that this approach is applicable for small-scale check dam systems where the time lag of flood routing can be neglected.As a pilot study,the current paper presents the analytical results for the Wangmaogou Check Dam System that has 22 dams connected in series and parallel.The algorithm clearly identified 7 surviving dams,with the remaining ones being overtopped for a storm of 179.6 mm in 12 h,which is associated with a return period of one in 200 years.     

OVERVIEW OF DAM GULLY EROSION RESEARCH

Traditionally gully erosion has been identified with the dissection of the landscape in agricultural settings but it is also recognized as a prevalent erosion feature in earthen dam auxiliary spillways and embankments. Flows through earthen spillways and over dam embankments, due to large rainfall events, have the potential to erode and breach the dam or spillway and result in catastrophic releases from the reservoir. The gully erosion process in an earthen spillway or on an embankment can be characterized by stages of initiation, development, and migration of a headcut. A headcut is defmed as a near vertical drop at the upstream end of a gully. The rate of headcut migration is important in determining the breach potential of an earthen spillway and dam embankment. A research program is being conducted to examine the gully erosion processes of earthen dam auxiliary spillways and embankments. This paper describes： l ） the unique test facilities constructed to examine the dominant factors affecting the erosion of earthen spillways and embankments; 2） the observations of the erosion processes and results to date; and 3） the predictive relationships that have been developed for dam gully erosion research at the ARS Hydraulic Engineering Research Unit laboratory in Stillwater, OK.     

Artificial neural networks and fuzzy logic based control of spillway gates of dams

In this paper a new and efficient control method based on fuzzy logic is proposed for real‐time operation of spillway gates of a reservoir during any flood of any magnitude up to the probable maximum flood. Artificial neural networks are used to model the non‐linear relationship among the main variables of the reservoir under consideration. In order to demonstrate the performance of the proposed method, we simulate the control system using different probable overflow hydrographs. The results of the proposed control method are compared with the results of the conventional control methods. The results obtained from the simulations indicate that the proposed method exhibits superior performance over the conventional reservoir flood operation, with much more flexible results. Copyright © 2004 John Wiley & Sons, Ltd.     

Bayesian optimal design of an avalanche dam using a multivariate numerical avalanche model

For snow avalanches, passive defense structures are generally designed by considering high return period events. However, defining a return period turns out to be tricky as soon as different variables are simultaneously considered. This problem can be overcome by maximizing the expected economic benefit of the defense structure, but purely stochastic approaches are not possible for paths with a complex geometry in the runout zone. Therefore, in this paper, we include a multivariate numerical avalanche propagation model within a Bayesian decisional framework. The influence of a vertical dam on an avalanche flow is quantified in terms of local energy dissipation with a simple semi-empirical relation. Costs corresponding to dam construction and the damage to a building situated in the runout zone are roughly evaluated for each dam height–hazard value pair, with damage intensity depending on avalanche velocity. Special attention is given to the poor local information to be taken into account for the decision. Using a case study from the French avalanche database, the Bayesian optimal dam height is shown to be more pessimistic than the classical optimal height because of the increasing effect of parameter uncertainty. It also appears that the lack of local information is especially critical for a building exposed to the most extreme events only. The residual hazard after dam construction is analyzed and the sensitivity to the different modelling assumptions is evaluated. Finally, possible further developments of the approach are discussed.     

Estimation of design flood using EWT and ENE metrics and uncertainty analysis under non-stationary conditions

Concepts of Expected Waiting Time (EWT) and Expected Number of Exceedances (ENE) have been presented in much literature for estimating the Design Flood (DF) under non-stationary conditions. The parameters of the EWT and ENE are generally no less than four, which inevitably leads to the uncertainty of the DF estimation. In this paper, the Bayesian method is proposed to analyze the impact of parameter estimation uncertainty on the EWT- and ENE-based estimation of the DF and Corresponding Design Reliability (CDR). In addition, a comparison analysis between the EWT and ENE is conducted in terms of the DF and CDR with or without a consideration being given to the impact of parameter uncertainty. In the case of giving no consideration to the impact of parameter uncertainty, the experiment results indicate that the EWT-based estimations are less than that of ENE in terms of DF and CDR in the case of a decreasing trend. While in the case of an increasing trend, the EWT-based estimations are bigger than that of ENE. In the case of considering the impact of parameter uncertainty, results in the case study show that the distribution of the EWT-based estimations of DF and CDR are left shifted compared to that of the ENE. Overall, the EWT-based estimations are significantly different from that of ENE in terms of DF and CDR. Therefore, it is necessary and open for further discussions about which metric will be optimal between the EWT and ENE for estimating the DF under non-stationarity.     

Hybrid method to update the design flood shape for spillway operation in a dam: case study of the Huites Dam,Mexico

To design and review the operation of spillways, it is necessary to estimate design hydrographs, considering their peak flow, shape and volume. A hybrid method is proposed that combines the shape of the design hydrograph obtained with the UNAM Institute of Engineering Method (UNAMIIM) with the peak flow and volume calculated from a bivariate method. This hybrid method is applied to historical data of the Huites Dam, Sinaloa, Mexico. The goal is to estimate return periods for the maximum discharge flows (that account for the damage caused downstream) and the maximum levels reached in the dam (measure of the hydrological dam safety) corresponding to a given spillway and its management policy. Therefore, to validate the method, the results obtained by the flood routing of the 50-year hydrograph are compared with those obtained by the flood routing of the three largest historical floods. Both maximum flow and elevation were in the range of values observed within 37.5–75 years corresponding to the length of the historical record.     

The Upper Limit Size of Reservoir-Induced Earthquakes

We showed the relation between the magnitude of induced earthquake and the reservoir storage and dam height based on the global catalog from 1967 to 1989 compiled by Ding Yuanzhang （1989）. By multiplying reservoir storage with dam height, we introduced a new parameter named EE. We found that the cases with specific EE and magnitude do not exceed a limit. Based on the discussion of its physics, we called EE the equivalent energy. We considered this limit as the upper limit of magnitude for reservoir-induced earthquakes. The result was proved by the recent cases occurring in China. This size limitation can be used as a helpful consideration for reservoir design.     

Geobag stepped spillway for check dams: A pilot study

Check dams are widely used worldwide for the soil and water conservation. Many of them have no spillways, resulting in frequent dam breach failures during heavy rainfall events. This paper proposes a new geobag stepped spillway for small check dams. The structure of the spillway primarily consists of a stepped chute created from geobags filled with local soil compacted and consolidated with a small amount of cement. The developed composite geosynthetic material exhibits water erosion resistance,...     

基于遗传算法优化神经网络权值的大坝结构损伤识别

针对传统 BP 神经网络存在着容易陷入局部极小点、训练时间太长等缺点,本文采用基于浮点编码的遗传算法,对 BP 神经网络的初值空间进行了遗传优化。用基于浮点编码的遗传算法来优化 BP 神经网络的权值,可得到最佳初始权值矩阵,并按误差前向反馈算法,沿负梯度搜索进行网络学习。文中以混凝土重力坝结构作为算例,用结构的模态频率变化作为网络的输入向量,结构的损伤位置作为输出向量,对网络进行了训练。仿真结果表明：遗传 BP 神经网络的收敛和诊断能力优于传统 BP 神经网络,可有效地运用到大坝结构的健康诊断与损伤识别中。