Abstract

For modern high earth dams, sufficient safety margin is considered in the designs of flood discharge capacity and dam crest elevation to prevent flood overtopping. However, for high earth dams which may induce catastrophic consequences, during their long operational period, extremely hazardous scenarios which could occur and threaten dam safety need to be considered. For the earth dams located in areas with intensive seismicity, there is a possible scenario that the release structures fail due to seismic landslides and gate failures caused by a severe earthquake when the flood begins to enter the reservoir. Thus, it is desirable to investigate the influence of failure duration of release structures on dam overtopping risk. Based on the Bayesian network, a methodology for overtopping risk analysis of earth dams considering effects of failure duration of release structures is proposed. The overtopping risk of the PBG earth-rockfill dam was analyzed to illustrate the methodology. The critical release structures which dominate the dam overtopping risk are identified. The dam overtopping risk is most sensitive to the failure duration of the spillway. The tolerable failure duration of the spillway is approximately 3 days, and when the failure duration of the spillway reaches 4 days, the dam overtopping risk drastically rises to an unacceptable level. The case study suggests that the proposed methodology could be helpful to analyze the influences of possible failure durations of release structures on dam overtopping risk and could facilitate preparation for emergency plans.

1. Introduction

Dams are important infrastructures for water resources management. There are more than 98000 dams in China, among which earth dams are more than 80000 [1]. Dam breaks have occurred frequently in the past 100 years and have caused huge losses of life and property [2]. Overtopping is one of the main causes of the earth dam break [35]. Overtopping risk analysis could provide useful information for decision making on reservoir operation and dam safety management.

Numerous research works on risk analysis of dam overtopping have been reported in the literature. As regards the uncertainties in overtopping risk analysis, the effects of uncertainties in flood, reservoir characteristics, outflow discharge, initial water surface level, wind velocity, dam height, and so on have been well investigated (e.g., [613]). Regarding approaches for overtopping risk analysis, a variety of approaches have been applied, such as stochastic differential equations (e.g., [14]), bivariate copula functions (e.g., [1517]), sequential uncertainty fitting method (e.g., [18]), Latin hypercube sampling (e.g., [19, 20]), maximum entropy method (e.g., [21]), and Bayesian network (e.g., [2224]).

In order to obtain hydropower, many high earth dams have been built in western China. Sufficient safety margin is considered in the designs of flood discharge capacity and dam crest elevation to prevent flood overtopping. In addition, the water level of the reservoir in flood season shall not exceed the specified limit water level. Therefore, the safety of the dams can be guaranteed even for the possible maximum flood. However, in western China, the occurrence frequency and intensity of earthquakes are high. Strong earthquakes may trigger massive rockfalls and landslides around the dams. For instance, the Wenchuan earthquake in 2008 caused more than 200000 landslides [2527]. The deposits formed by rockfalls or landslides may block release structures. And the gates of the release structures may not be opened due to excessive seismic distortion. For example, the landslides induced by the Wenchuan earthquake buried the inlet and/or outlet of the release structures of the Yuzixi, Yingxiuwan, and Taipingyi hydropower stations. The earthquake led to a temporary malfunction of the flood discharge gates of the Taipingyi, Futang, Jiangsheba, and Gengda sluices.

Therefore, although the overtopping risk of the dams under normal operating conditions is quite small, if an extremely strong earthquake in the flood season leads to failure of the release structures and the flood discharge function of the release structures cannot be restored in time, dam overtopping may occur. Thus, for the earth dams located in areas with intensive seismicity, especially in the flood season, it is necessary to consider the impact of the failure of release structures caused by seismic landslides on dam overtopping risk. The important effect of the unavailability of the spillway on dam overtopping risk has been investigated (e.g., [2831]). In fact, the flood discharge capability of the failed release structures has to be recovered as quickly as possible. Hence, the influence of the failure duration of release structures on dam overtopping risk is worth being investigated. Based on the Bayesian network, a methodology for overtopping risk analysis of earth dams considering effects of failure duration of release structures is proposed in this paper. A case study is presented to illustrate the methodology.

2. Methodology for Overtopping Risk Analysis considering Failure Duration of Release Structures

2.1. Outline of the Overtopping Risk Analysis

In this paper, the scenario for overtopping risk analysis of earth dams is supposed as follows: when a flood begins to enter the reservoir, a strong earthquake occurs near the dam site. The earthquake may cause landslides around the dam and the resulting deposits may block the release structures, such as spillway and flood discharging tunnel. In addition, the gates of the release structures may fail to open due to the seismic action. Although emergency measures are taken to recover the discharge capability of the release structures, before discharge capability recovery, the failure duration of the respective release structure should increase the dam overtopping risk.

For simplicity, the probabilistic seismic hazard is not addressed in this study. It is assumed that a strong earthquake definitely occurs near the dam site and the intensity of the earthquake is adequately large so that it is possible to induce landslides. Additionally, since this study focuses on the effects of failure duration of release structures on dam overtopping risk and the effects of uncertainties in other affecting factors (such as flood, reservoir characteristics, initial water surface level, and dam height) have been well investigated as mentioned above, only the uncertainties in seismic landslides, gate failure, and failure duration of release structures are considered in the analysis.

Dam overtopping is the result of a series of uncertain events (event chain) with cascading causality. For example, an earthquake occurs when a flood is entering a reservoir, the earthquake causes landslides, subsequently, the landslide deposits block the spillway, the discharge capability of the spillway cannot be restored timely, and then dam overtopping happens. For practical problems, multiple overtopping event chains may exist. In this study, the Bayesian network (i.e., a directed acyclic graph composed of nodes and directed arcs) is adopted to characterize complex chains of overtopping event. In the Bayesian network, the probability of occurrence of an event represented by a node (child node) can be calculated by conditional probability formula in terms of the probabilities of occurrence of events represented by the parent nodes of the child node [32].

Suppose that an earth dam hub has M release structures and there are Ni seismic landslides that may lead to failure of the ith release structure. In addition, gate failure may also result in failure of the release structure. Figure 1 shows the Bayesian network model of overtopping risk analysis for earth dams in this study. The procedure of overtopping risk analysis is as follows:(a)According to the topographical and geological conditions as well as the layout of release structures for a specific earth dam hub, potential seismic landslides and the release structures which may be affected by the landslides are determined. Then, the Bayesian network as shown in Figure 1 is established.(b)Estimating the probabilities of seismic landslide and gate failure: The method adopted in this study is detailed in Section 2.2.(c)The failure duration of release structure related to a seismic landslide or gate failure is treated as a random variable with uniform or triangular distribution, and the relevant distribution parameters are estimated by experts. The probability distribution of failure duration of a release structure which may be affected by multiple landslides and gate failure is determined in terms of the probability distributions of failure duration associated with the landslides and gate failure, as well as the probabilities of occurrence of the combinations of the landslides and gate failure. A detailed description is given in Section 2.3.(d)Considering the influence of failure durations of the release structures, reservoir routings are implemented and the conditional probabilities of dam overtopping are determined. The details of the computing process are presented in Section 2.4(e)The probability of dam overtopping is computed with the Bayesian network, given the probabilities of seismic landslides and gate failures, the probability distributions of failure duration of the release structures, and the conditional probabilities of overtopping corresponding to the given inflow hydrograph and the failure durations of the release structures.

2.2. Estimation of Probabilities of Seismic Landslide and Gate Failure

The probabilities of seismic landslide and gate failure can be quantitatively calculated using the reliability analysis method (e.g., [3335]). Due to the difficulties in the characterization of the probability distributions for material properties and loads, as well as the complexity of the computation process, the application of probabilistic methods in practical problems is still limited.

The methods for estimation of probabilities of seismic landslide and gate failure are beyond the scope of the paper, and the results of the probability estimation do not affect the implementation of the methodology for overtopping risk analysis. For simplicity, the method presented by [36] is used to estimate the probabilities of seismic landslide and gate failure. The procedure is as follows.

First, the seismic stabilities of the slopes and the seismic resistance of the gates that may affect the release structures are ranked by expert judgment. Then, the linguistic descriptions of the likelihood of landslide and gate failure under strong earthquakes are presented. Finally, according to the prescribed probabilities corresponding to the linguistic descriptions of likelihood (as shown in Table 1), the probabilities of seismic landslide and gate failure are determined.

2.3. Probability Distribution of Failure Duration of Release Structure

For scarce knowledge and historical data for failure duration of release structure (especially for a specific project) affected by seismic landslide and gate failure after a strong earthquake, in this study, the probability distribution of failure duration of release structure related to a landslide or gate failure is assumed to follow uniform distribution or triangular distribution [37] and the relevant distribution parameters are estimated by experts.

The uniform distribution for failure duration of release structure is as follows:where and are the lower and upper limits of the failure duration T.

The triangular distribution for failure duration of release structure is expressed aswhere is the most likely value of T.

For a release structure which may be affected by multiple landslides and gate failure, the probability distribution of failure duration of the release structure is determined in terms of the probability distributions of failure duration related to the landslides and gate failure, as well as all possible combinations of the landslides and gate failure and associated probabilities. For the case where multiple landslides and gate failure simultaneously occur, it is assumed that only the probability distribution of failure duration with the largest values of and is considered.

For instance, suppose a release structure which may be affected by two landslides (denoted by LS1 and LS2) and gate failure (denoted by GF); the probabilities of LS1, LS2, and GF are P(LS1), P(LS2), and P(GF), respectively; the probability distributions of failure duration related to LS1, LS2, and GF are fLS1(T), fLS2(T), and fGF(T), respectively; in addition, fLS1(T) and fGF(T) have the largest and smallest values of and , respectively. Thus, the probability distribution of failure duration of the release structure is determined as f(T) = fLS1(T) × P(LS1) × P(LS2) × P(GF) + fLS1(T) × P(LS1) × [1−P(LS2)] × P(GF) + fLS1(T) × P(LS1) × P(LS2) × [1−P(GF)] + fLS1(T) × P(LS1) × [1−P(LS2)] × [1−P(GF)] + fLS2(T) × [1−P(LS1)] × P(LS2) × P(GF) + fLS2(T) × [1−P(LS1)] × P(LS2) × [1−P(GF)] + fGF(T) × [1−P(LS1)] × [1−P(LS2)] × P(GF) and f(T = 0) = [1−P(LS1)] × [1−P(LS2)] × [1−P(GF)].

2.4. Reservoir Routing considering Failure Duration of Release Structure and Conditional Probability of Dam Overtopping

In the course of reservoir routing, the water level in the reservoir rises as the outflow through release structures is less than the inflow. When the water level exceeds the dam crest, the dam overtopping occurs. Given the water level-storage capacity curve, the flood hydrograph, and the discharge curve for release structures, the water level in the reservoir can be determined by solving the water balance equation [23]. In reservoir routing, the discharge of a failed release structure is set to be null until its functionality is recovered. An illustrative example of reservoir routing considering failure duration of release structure is presented as follows.

Assume that an earth dam hub has n release structures (RS1, RS2, …, RSn). When a flood begins to enter the reservoir, all the release structures fail due to a strong earthquake; the failure duration of each release structure is 0 < t1 < t2<<tn; the corresponding reservoir routing is shown in Figure 2. During 0 ≤ t < t1, all the release structures do not participate in flood discharge, so that the total discharge q = 0. During t1 ≤t < t2, only the discharge capability of the release structure RS1 is recovered; thus, q = q1, q1 is the discharge of the release structure RS1. During t2 ≤ t < t3, the discharge capability of the release structure RS2 is recovered; hence, q = q1 + q2, q2 is the discharge of the release structure RS2. When t > tn, the release structures all recover discharge capabilities; therefore, q = q1 + q2++qn. The illustrative example demonstrates the case where the failure duration of the release structure is different for each other. For the case where a part of or all of the release structures have identical failure duration, the reservoir routing process is similar. If a release structure does not fail, the corresponding failure duration is defined as zero during reservoir routing.

Assume that there are n release structures, for each of which m failure durations are considered in reservoir routing. Thus, there are mn combinations of failure durations of the release structures for reservoir routing. Given flood hydrograph and a combination of failure durations of the release structures, the water level can be determined by reservoir routing and the corresponding conditional probability of dam overtopping is expressed aswhere the failure durations of the release structures are represented by , and is the dam crest elevation.

3. Case Study

3.1. Brief Introduction of the PBG Hydropower Station

The PBG hydropower station is located on the Dadu River in Sichuan Province, China. The Dadu River is cut off by an earth-rockfill dam. A spillway, a flood discharging tunnel, and a power station diversion tunnel are arranged in the left bank, and an emptying tunnel is arranged in the right bank. The dam crest elevation is 856.00 m, and the maximum dam height is 186 m. The normal water level of the reservoir is 850.00 m, the limited water level in flood season is 841.00 m, the dead water level is 790.00 m, and the total storage capacity is 5.337 billion m3. The layout of the hydropower station is shown in Figure 3.

The annual average temperature in the Dadu River basin is 13∼18°C, and the annual rainfall is 1400∼1700 mm. The rainfall in the wet season, from May to October, accounts for 80%∼90% of the annual rainfall, and the maximum daily rainfall is more than 100 mm, even more than 400 mm. Floods in the Dadu River are mainly formed by precipitation, and the annual maximum flood peak generally occurs in flood season (from June to September). A typical flood process lasts for about 15 days, in which the flood volume in 7 days accounts for more than 53% of the total flood volume, the flood volume in 3 days accounts for more than 48% of the flood volume in 7 days, and the flood volume in 1 day accounts for more than 35% of the flood volume in 3 days. The maximum flood peak which had been observed at the dam site is 8520 m3/s.

Several large-scale active faults are distributed around the hydropower station, as shown in Figure 4. The PBG hydropower station is located in the block delimited by the faults. The historical earthquakes around the PBG hydropower station are shown in Figure 5. According to the probabilistic seismic hazard analysis, the probability of occurrence of earthquakes with a magnitude of more than 7.5 near the dam site is 0.23% in the next 50 years. And the peak value of horizontal acceleration at the dam site with an exceedance probability of 1% in the next 100 years is 0.48 g.

3.2. Bayesian Network Model for Overtopping Risk Analysis of the PBG Dam

The release structures of the PBG hydropower station include spillway, flood discharging tunnel, power station diversion tunnel, and emptying tunnel. As mentioned above, severe earthquakes are quite likely to occur in the area near the PBG hydropower station in the future. Although the excavation slopes of the release structures have been reinforced, the natural slopes above the excavation slopes may slide under the action of a strong earthquake, and the landslides may result in failures of the release structures. For risk assessment, an extreme case is assumed that the PBG hydropower station would simultaneously encounter a strong earthquake and a big flood. Considering the influences of seismic landslide and gate failure on the release structures of the PBG hydropower station and the failure durations of the release structures, a Bayesian network model for overtopping risk analysis of the earth-rockfill dam is established, as shown in Figure 6.

3.3. Model Quantification
3.3.1. Probabilities of the Seismic Landslides and Gate Failures

The method described in Section 2.2 is used to estimate the probabilities of the seismic landslides that may affect the release structures. First, the stabilities of the natural slopes above the release structures under a strong earthquake are ranked by expert judgment, then the likelihoods of seismic landslides are described with linguistic description, and finally, the probabilities of seismic landslides are determined in terms of the probabilities related to the linguistic descriptions of likelihood (see Table 1). The probabilities of gate failures induced by an earthquake are estimated similarly. Tables 2 and 3 show the ranks of slope stability and the ranks of seismic resistance of the gates, together with the linguistic descriptions of the likelihood of seismic landslide and gate failure, while the corresponding probabilities are listed in Table 4.

3.3.2. Estimation of Probability Distribution Parameters for Failure Duration of Release Structure

For comparison, the probability distribution of failure duration of release structure related to a landslide or gate failure is assumed to follow the uniform distribution and triangular distribution. The probability distribution parameters (the lower limit of failure duration), (the upper limit of failure duration), and (the most likely value of failure duration) are estimated by experts and are listed in Table 5.

3.3.3. Reservoir Routing and Determination of Conditional Probability of Dam Overtopping

The maximum design flood hydrograph for the PBG hydropower station is shown in Figure 7. The flood process lasts 392 hours, and the flood peak appears on the 8th day, which is 15250 m3/s. The water level-storage capacity curve and the discharge curve for release structures are shown in Figures 8 and 9, respectively.

It is assumed that when the flood begins to enter the reservoir (the water level of the reservoir is keeping as 841.00 m in flood season), a severe earthquake occurs, and the resultant landslides and gate failures subsequently cause failures of the release structures. Considering the influence of failure durations of the release structures, the method presented in Section 2.4 is used to compute the maximum water level in the course of reservoir routing. In order to reduce the computational effort, discrete failure duration (i.e., 0d, 1d, …, 16d) is adopted in reservoir routing. Therefore, the total number of reservoir routings is mn = 174 = 83521, where m is the number of failure durations and n is the number of release structures. According to the results of 83521 reservoir routings, the conditional probabilities of dam overtopping are determined using Equation (3).

Figure 10 shows two typical results of reservoir routing: (a) the failure durations of the spillway, flood discharging tunnel, power station diversion tunnel, and emptying tunnel are 1 d, 3 d, 5 d, and 7 d, respectively, the corresponding maximum water level HR is lower than the dam crest elevation (  = 856.00 m), and thus, dam overtopping does not occur; (b) the failure durations of the spillway, flood discharging tunnel, power station diversion tunnel, and emptying tunnel are 7 d, 5 d, 3 d, and 1 d, respectively; in this case, dam overtopping occurs on the 4th day of the flood process.

4. Results and Discussion

For the case where the probability distribution of failure duration of release structure related to a landslide or gate failure is assumed to be uniformly distributed (referred to as Case I), the overtopping risk of the PBG dam is analyzed using the Bayesian network model (as shown in Figure 6) with the model quantification results presented in Section 3. The GeNIe software [38] is used to evaluate the Bayesian network model. Figure 11(a) presents the result of the overtopping risk analysis for Case I. The probability of the dam overtopping is about 10% considering the influence of failure durations of the release structures.

The intermediate nodes in Figure 11(a) show the probability distributions of failure duration of the release structures affected by multiple landslides and gate failure, which are determined in terms of the uniformly distributed failure duration related to a landslide or gate failure, as well as the probabilities of the combinations of the landslides and gate failure. It is observed that the probabilities that spillway, flood discharging tunnel, power station diversion tunnel, and emptying tunnel do not fail are 12.25%, 2.50%, 1.75%, and 12.80%, respectively. The average failure durations for each release structure are 2.62 d, 4.13 d, 4.41 d, and 9.12 d, respectively. The standard deviations of failure durations are 1.44 d, 1.46 d, 1.40 d, and 4.45 d, respectively.

A sensitivity analysis is carried out to investigate the impact of each bottom node (i.e., the respective slope and gate) in the Bayesian network model on the dam overtopping risk for Case I. With the Bayesian network model shown in Figure 11(a) as a reference point, each bottom node is assigned, in turn, a probability of 1 and 0, respectively, for being in the “failure” state, keeping the remaining nodes in their existing state and the resulting changes in the probability of dam overtopping are calculated. The sensitivity of dam overtopping risk to the bottom nodes is computed by normalizing the respective change in the probability of dam overtopping over the total change (i.e., sum of the changes in the probability of dam overtopping from each bottom node) and is presented in Figure 12(a). Since the spillway and flood discharging tunnel are major release structures for the PBG dam and the failure durations of both release structures related to seismic landslides are relatively long, the changes in failure probability of spillway slope and inlet slope of flood discharging tunnel have the greatest impacts on the dam overtopping risk. The effects of the changes in failure probability of the gates are minor as the failure durations of the release structures related to gate failure are relatively short.

Figure 11(b) shows the result of the overtopping risk analysis for Case II where the probability distribution of failure duration of release structure related to a landslide or gate failure is assumed to be triangularly distributed. The corresponding model quantification is presented in Section 3. The probability distributions of failure duration of the release structures affected by multiple landslides and gate failure are shown in the intermediate nodes in Figure 11(b). The probabilities that spillway, flood discharging tunnel, power station diversion tunnel, and emptying tunnel do not fail are 12.25%, 2.50%, 1.75%, and 12.80%, respectively. The average failure durations for each release structure are 2.63 d, 4.11 d, 4.43 d, and 9.07 d, respectively. The standard deviations of failure durations are 1.36 d, 1.28 d, 1.18 d, and 4.17 d, respectively. The probability of the dam overtopping is about 7% considering the influence of failure durations of the release structures. It is shown that the overtopping risk of the PBG dam is not sensitive to the probability distribution type of failure duration by comparing the computational results for Case I and Case II. The sensitivity analysis also reveals that the effects of the changes in failure probability of spillway slope and inlet slope of flood discharging tunnel on the dam overtopping risk are predominant and those of the gates are minor.

In order to investigate the effect of the failure duration of each individual release structure on the dam overtopping risk, the probabilities of dam overtopping corresponding to the failure duration of 0 d, 1 d, …, 16 d for each individual release structure (keeping the failure durations of the remaining release structures in their existing probability distributions) are calculated and are shown in Figure 13. It can be seen that the probabilities of dam overtopping increase with the failure duration of each individual release structure. The probabilities of dam overtopping keep unchanged when the failure duration is greater than 8 days, since it is assumed that the release structures fail due to a severe earthquake when the flood begins to enter the reservoir and the flood flow decreases after 8 days. The increase in the failure duration of the spillway affects the dam overtopping risk quite significantly. When the failure duration of the spillway is less than 4 days, the corresponding dam overtopping risk is quite small, the maximum probability of dam overtopping is 8.1 × 10−5 (Case II). However, when the failure duration of the spillway reaches 4 days, the probability of dam overtopping drastically rises to about 50%. After then, the probability of dam overtopping continues to rapidly increase with the failure duration of the spillway. The probability of dam overtopping achieves the maximum value of about 96% when the failure duration is 8 days. The effect of the increase in the failure duration of the flood discharging tunnel on the dam overtopping risk is fairly important. When the failure duration of the flood discharging tunnel reaches 4 days, the probability of dam overtopping jumps form 3.5 × 10−4 to 13.4% and remains constant until the failure duration reaches 8 days, the corresponding probability of dam overtopping sharply increases to 53.5% (Case II). The influences of the increase in the failure durations of the power station diversion tunnel and emptying tunnel on the dam overtopping risk are relatively small compared to the spillway and flood discharging tunnel.

5. Conclusions

Overtopping is one of the main causes of earth dam break. For modern high earth dams, sufficient safety margin is considered in the designs of flood discharge capacity and dam crest elevation to prevent flood overtopping. However, for the earth dams located in areas with intensive seismicity, there is a possible scenario that the release structures fail due to a severe earthquake when the flood begins to enter the reservoir. Thus, it is desirable to investigate the influence of failure duration of release structures on dam overtopping risk. For overtopping risk analysis of earth dams considering effects of failure duration of release structures, a Bayesian network model and methods for model quantification are proposed.

As a case study, the overtopping risk of the PBG earth-rockfill dam was analyzed. The spillway and flood discharging tunnel are critical release structures which dominate the dam overtopping risk. It is found that the dam overtopping risk is most sensitive to the failure duration of the spillway. The probability of dam overtopping is practically negligible when the failure duration of the spillway is less than 4 days whereas the dam overtopping risk drastically rises to an unacceptable level when the failure duration of the spillway reaches 4 days. It implies that the tolerable failure duration of the spillway is approximately 3 days and intensive efforts have to be made to timely recover the discharge capability of the spillway for eliminating overtopping risk. For earth dams which may induce catastrophic consequences, extremely hazardous scenarios (such as the one considered in this study) which could occur and threaten dam safety during a long operational period of dams need to be considered.

The case study suggests that the proposed methodology could be helpful to analyze the influences of possible failure durations of release structures on dam overtopping risk and could facilitate preparation for emergency plans. For improving the proposed methodology, future works may be focused on the scheme for effectively eliciting and combining expert estimations of the probability distribution of failure duration of release structure as well as the efficient computation of reservoir routing considering failure duration of release structure.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

Financial support provided by the National Key R&D Program of China (no. 2018YFC0407103) is gratefully acknowledged.