J. M. Vali Samani; H. Radmehr; M. Delavar
Abstract
Introduction: The greatest part of constructed dams belongs to embankment dams and there are many examples of their failures throughout history. About one-third of the world’s dam failures have been caused by flood overtopping, which indicates that flood overtopping is an important factor affecting ...
Read More
Introduction: The greatest part of constructed dams belongs to embankment dams and there are many examples of their failures throughout history. About one-third of the world’s dam failures have been caused by flood overtopping, which indicates that flood overtopping is an important factor affecting reservoir projects’ safety. Moreover, because of a poor understanding of the randomness of floods, reservoir water levels during flood seasons are often lowered artificially in order to avoid overtopping and protect the lives and property of downstream residents. So, estimation of dam overtopping risk with regard to uncertainties is more important than achieving the dam’s safety. This study presents the procedure for risk evaluation of dam overtopping due to various uncertaintiess in inflows and reservoir initial condition.
Materials and Methods: This study aims to present a practical approach and compare the different uncertainty analysis methods in the evaluation of dam overtopping risk due to flood. For this purpose, Monte Carlo simulation and Latin hypercube sampling methods were used to calculate the overtopping risk, evaluate the uncertainty, and calculate the highest water level during different flood events. To assess these methods from a practical point of view, the Maroon dam was chosen for the case study. Figure. 1 indicates the work procedure, including three parts: 1) Identification and evaluation of effective factors on flood routing and dam overtopping, 2) Data collection and analysis for reservoir routing and uncertainty analysis, 3) Uncertainty and risk analysis.
Figure 1- Diagram of dam overtopping risk evaluation
Results and Discussion: Figure 2 shows the results of the computed overtopping risks for the Maroon Dam without considering the wind effect, for the initial water level of 504 m as an example. As it is shown in Figure. 2, the trends of the risk curves computed by the different uncertainty analysis methods are similar. As it can be seen, the risk curves computed by the LHS are slightly higher than those curves computed by the MCS method. Also as it is observed, the differences between risk values of the two methods increase in longer return periods. Variations of overtopping risk with increasing the initial water level and return period related to overtopping risk in the 2-year return period for the initial water level of 470 m are shown in Table1. The results show that elongation of return period plays a more important role in increasing the risk, than the increase of initial water level.
T Method 2→2 2→50 2→100 2→1000 2→5000 2→10000
470→470 MCS 1 5 9 23 42.36 58
470→478 2 7 15.6 37 58.34 79
470→485 5.6 13.6 28.6 55.6 85.67 112.6
470→493 10.3 32.6 54 95.6 127.34 152
470→504 40.3 83 117.3 165 200.34 224.3
470→470 LHS 1 5.34 11 25.3 43 60.3
470→478 2.3 8.6 18 39.3 60.67 84
470→485 5.3 17.3 32.6 58.3 89 114.6
470→493 13.3 37.6 57.6 97 133.34 160.3
470→504 41.6 87.3 119.6 173 205 233.3
Figure 2- Overtopping risk in the initial water level of 504 m, without considering the wind effect
Conclusions: This study applies MCS and LHS methods to analyze the uncertainty and evaluate the dam overtopping risk consideringthe uncertainties in input variables, such as quintile of flood peak discharge, initial levels of water and spill coefficients. The results show that the uncertainty of water level calculated by MCS is higher than that calculated by LHS. In addition, the overtopping risk calculated by LHS is higher than that calculated by MCS. Furthermore, the increase of inflow rate influences the variations of the overtopping risk more than the increase of the return period. In addition, evaluation of the results indicates that the overtopping risk is an important issue in the Maroon dam. So, a comprehensiverisk analysis procedure in conjunction with uncertainty gives very important information for decision makers to make better judgments in dam operation based on uncertainty in inputs.
mahsa noori; Saeed Reza Khodshenas; H. Rezaeepajand
Abstract
Introduction: Dam failure and its flooding is one of the destructive phenomena today. Therefore, estimating the peak outflow (QP) with reasonable accuracy and determining the related flood zone can reduce risks. Qp of dam failure depends on important factors such as: depth above breach (Hw), volume ...
Read More
Introduction: Dam failure and its flooding is one of the destructive phenomena today. Therefore, estimating the peak outflow (QP) with reasonable accuracy and determining the related flood zone can reduce risks. Qp of dam failure depends on important factors such as: depth above breach (Hw), volume of water above breach bottom at failure (Vw), reservoir surface area (A), storage (S) and dam height (Hd). Various researchers have proposed equations to estimate QP. They used the regression method to obtain an appropriate equation. Regression is a mathematical technique that requires initial test and diagnosis. These researchers present a new regression model for a better estimation of Qp.
Materials and Methods: The data used in this study are related to 140 broken dams in the world for 34 of which sufficient data are available for analysis. Dam failure phenomenon is a rapidly varied unsteady flow that is explained by shallow waters equations. The equations in the one-dimensional form are known as Saint-Venant equations and are based on hydrostatic pressure distribution and uniform flow under rectangular steep assumption. Although hydraulic methods to predict the dam failure flood have been developed by different software, due to the complex nature of the problem and the impossibility of considering all parameters in hydraulic analysis, statistical methods have been developed in this field. Statistical methods determine the equations that can approximate the required factors from the observed parameters. Multiple regression is a useful technique to model effective parameters in Qp, which can examine the statistical aspects of the model. This work is done by different tests, such as the model coefficients necessity test, analysis of variance table and it creates confidence intervals. Data analysis in this paper is done by SPSS 16 software. This software can provide fit model, various characteristics and related tests in the Tables.
Results and Discussion:This paper proposes a new relationship with better estimation of discharge peak (Qp) based on Hw and Vw factors. Results showed how to choose the appropriate model (fitting the model) and the initial required tests, according to the diagnostic model. And it compares the estimated error (relative efficiency) of the researchers’ models with the proposed models. The number of models can be classified to three convenient linear, multiplicative and transformed bases on Vw, Hw and Qp (nonlinear terms Qp). The best models for each of the three models were selected. Their corrected determination coefficients (Adj R2) are close together and are between 0.86 until 0.864. The relative efficiency criteria based on the root mean square error (RMSE) was used to determine the best model. This standard was also used for other researchers’ models. RMSE of the three models presented in this article is lower than that of other models (from 745 to 759). Diagnostics analysis of the three models is not possible due to the large volume, so some statistical analysis for the model 2 are presented in detail. The results are given in the following Tables. Test level has been assumed to be 5%. From the point view of hydraulics, it can be said that the final equation for Qp should be proportional to Hw 1.5. So although the model (2) has the lowest RMSE, but the model (3) of the hydraulics viewpoint seems more logical and its RMSE is not very different from the model (2), so this model can be selected as the best model. Figure 1 show diagnostics diagrams of model (3). The right Figure shows the homogeneity of residuals (follow the normal law) as a histogram. This homogeneity is confirmed by the crouch graph (center Figure). The left graph shows the stabilization of residual variance. According to the preliminary and diagnostics tests results, the model (3) has been selected. Its determination coefficient (0.864) also shows good strength.
Table 1- Top models presented in this research
Model1
Model2
Model3
,
Note:
Table 2- Statistical characteristics of the proposed models
model Adjusted R Square Durbin
Watson F VIF Std.
Residual Cook's Distance Centered Leverage
1 0.862 1.716 104.383 1.283 [-1.975 , 2.908] [ 0,0.569] [0,0.363]
2 0.860 1.744 102.545 1.283 [-1.824 , 2.834] [0,0.608] [0,0.363]
3 0.864 1.687 211.048 1 [-2.202 , 2.699] [0,0.527] [0,0.335]
Figure 1- Model 3 diagnostics pattern diagrams: histogram (right), crouch diagram (middle) the estimated residuals (left)
Conclusion: In this study, data from 140 broken dams were used to provide an appropriate model for estimating the peak outflow of dam failure. Standard statistical principles including preliminary tests, diagnostic and the efficiency of the models are the innovations of this paper. Analysis showed that the three models are competitive, and that the best of them was selected. The determined coefficient of these models was from 0.86 to 0.864 ranges. Relative efficiency was calculated by the RMSE index. The results showed that these models are more accurate than the models presented by other researchers. The model (3) was presented in this research, the best result was estimated for Qp and its error was less than the other models.