Research on the Overtopping-Induced Breaching Mechanism of Tailings Dam and Its Numerical Simulation
The centrifugal model test system is utilized to reappear the breaching process of a tailings dam. The results show that there are significant differences in breaching processes between tailings dams and earth-rockfill dams. Influenced by the special physical and mechanical properties of tailings, the breaching processes of tailings dams are of short duration, and the flood flow through the breach increase rapidly. Based on the centrifugal test results, considering the special physical and mechanical properties of tailings, a mathematical model for tailings dam breaching is established. In the model, the vertical undercutting and horizontal expansion are simulated utilizing the erosion rate formula which is derived from shear stress principle of water flow; limit equilibrium method is used to analyze the breach slope stability. In addition, the model adopts iterative calculation method to simulate the whole breaching process. The analysis results show that, the proposed model is applicable in modeling the dam breaching process for tailings dam due to overtopping.
Tailings dams are the places for storing discarded tailings, which are usually formed by damming the valley. China owns the largest amount of tailings dams in the world. By the end of 2015, 8869 tailings dams were scattered all over China. Among these tailings dams, the number of upstream tailings dams is up to 1425 . Compared to other kinds of tailings dams, the upstream tailings dam has the lowest safety factor [2, 3]. Through the statistical analysis of 3500 tailings dams worldwide, scholars found that breaching risk of tailings dams was more than 10 times that of earth-rockfill dams . In recent years, extreme weather events and earthquake disasters occurred frequently, which increased the dam breaching probability of tailings dams. In 2000, the Baia Borsa gold mine tailings dam in Romania breaching caused by rainstorm and snowmelt, which released 100,000 cubic meters of tailings into the Lapus and Somes tributaries of the river Tisa, damaged the ecosystem, and killed aquatic lives . In 2011, the Kayakari tailings dam in Japan breaching due to tailings liquefaction caused by earthquake, the discharged tailings destroyed many houses in the lower reaches . In 2014, the Mount Polley mine tailings dam breaching due to dam foundation instability, about 4.5 million cubic meters of tailings flowed into Polley Lake and seriously polluted the environment . Thus, once the tailings dam breaching occurs, it will seriously threaten the surrounding ecological environment and the safety of people’s lives and property. To reduce consequences of the tailings dam breaching, many scholars have carried out studies on the breaching process of tailings dams and achieved some results [8–13]. However, these studies still use the research methods applied in earth-rockfill dams, which cannot accurately reflect the particularity of tailings dams. In addition, due to the lack of the simulation method agreed with the special physical and mechanical properties of tailings, the lose evaluation of tailings dam breaching and the proposal of emergency plan are short of reliable evidence. These weak points seriously restrict the advancement of disaster prevention and mitigation level of tailings dams. Therefore, it is necessary to study the breaching mechanism of tailings dams due to overtopping for revealing the failure law of them.
In this study, centrifugal model tests are conducted to reveal the breaching mechanism of tailings dam. Based on the centrifugal tests results, a mathematical model for tailings dam breaching is established, which can consider the special physical and mechanical properties of tailings. In the model, the vertical downcutting and lateral expansion are simulated utilizing the erosion rate formula derived from shear stress principle of water flow; limit equilibrium method is used to analyze the breach slope stability. In addition, the model adopts iterative calculation method to simulate the whole breaching process. The technical details are described in the following sections.
2. Centrifuge Model Tests of Tailings Dam Breaching
Since the breaching process of actual tailings dam is difficult to observe and reproduce , model test is an important way to study the breaching mechanism of tailings dams. According to the document retrieval results, now there is only a few small-scale flume model tests [15–19] preliminarily discussing the catastrophic mechanism of tailings dams. However, the difference between these models and actual tailings dams in the stress level is too large. Whether these models can reasonably reflect the breaching process of actual tailings dams is worthy of further study. Comparing with flume model tests, centrifugal model tests can improve the stress level of tailings dam models through changing centrifugal acceleration. Besides, centrifugal model tests can reappear the breaching process of tailings dam in short time . Therefore, it is reasonable to study the breaching mechanism of tailings dam by centrifugal model tests.
2.1. Design and Manufacture of the Test Model
The test equipment consists of the following three parts: (1) 400 gt large geotechnical centrifuge, the maximum radius of the centrifuge is 5.5 m, the effective radius is 5.0 m, the maximum centrifugal acceleration is 200 g, and the effective load is 2 t. (2) Flow control system, the system is designed to adopt the idea of rotating water ring, which can provide an experimental flow of 0.01∼0.05 m3/s (prototype 100∼500 m3/s) at 100 times of gravity acceleration, with a duration up to 20 min (prototype 33.3 h). (3) Special test box integrated data and image acquisition system, the internal effective size of the test box is 1.2 m × 0.4 m × 0.8 m (length × width × height). The integrated data and image acquisition system can capture and record the whole breaching process of the tailings dam under the condition of high centrifugal acceleration, the test equipment is shown in Figure 1. The working principle of centrifuge and the derivation of similarity criterion are referenced to related literature [21, 22]. For the centrifugal model test at N times of gravity acceleration, the similarity criterion of common physical quantities is shown in Table 1.
The physical and mechanical properties of tailings are obviously different from those of general cohesive soil and noncohesive soil . As a kind of artificial soil, tailings have the characteristics of fine grain, low particle strength, and loose and large porosity [24, 25]. Therefore, the tailings from a surrounding tailings dam are chosen as dam construction materials for the centrifugal model tests. Considering the gradation feature of tailings, the tailings are sampled from different areas and depths of the tailings dam separately. Then, all the sampled tailings are mixed together, the density of mixture tailings is 2.08 g/cm3, the water content is 16.5%, the specific gravity is 3.33, and the content of clay particles is 5%. In the centrifugal model tests, the gravity acceleration is set as 50 g, so the size ratio between model dam and prototype is 1 : 50. The height of model dam is 0.4 m and the width of dam crest is 0.3 m. Considering the storage capacity of actual tailings dam, the upstream reservoir is simulated by a water box with upper opening. The internal effective size of the water box is 0.2 m × 0.4 m × 0.5 m. In order to simulate the permeable boundary, the contact section between the water box and the dam is made of perforated plates (Figure 2). Since the complexity of the breaching mechanism of tailings dam, the centrifugal model tests were carried out two groups for different dam slopes. Taking one group as an example, in this test, the downstream slope ratio is 1 : 3.0, and the upstream inflow rate is 1.7 × 10−3 m3/s. Limited by the size of model box, the tests are mainly aimed at the breaching section of tailings dam, so the downstream slope is cut off at the height of 0.2 m. The tailings dam model is shown in Figure 3.
2.2. Test Results and Analysis
According to the results of the tailings dam breaching tests, the breaching process can be divided into four stages, and the details are shown in Figures 4(a)–4(d). Stage 1: water flows over the crest of tailings dam and scours the dam body. Under erosion of the flow, many slender gullies are formed on the downstream slope of the tailings dam, as shown in Figure 4(a). Stage 2: the slender gullies on the right side of the downstream slope gradually expand and merge to form a breach channel. After the breach channel formed, the flow begins to discharge along it. At the end of this stage, a breach is formed in the right side of tailings dam crest, as shown in Figure 4(b). Stage 3: the breach continues to grow vertically and laterally. Since the flow in the breach channel is not cross the full section, erosion of the flow is concentrating on the lower region of the breach. With the breach growing in vertical direction, the breach side-slope angle continues to increase and finally results in the breach collapsing. In this stage, the discharge flow increases gradually and reaches its peak, as shown in Figure 4(c). Stage 4: growth of the breach gradually slows down until stops. At the end of this stage, when the bed shear stress in the breach is smaller than critical stress required to initiate detachment for the material, the breach erosion stops. After the dam breaching stopping, the final breach is shown in Figure 4(d).
To simplify the calculation, the measured parameters are converted into the corresponding numerical values under the constant gravity condition according to the centrifugal acceleration. Therefore, the upstream inflow rate is 4.25 m3/s, the peak flow through the breach is 21.88 m3/s, the time to peak flow is 0.13 h after dam breaching, and the breaching time is 2.56 h (Figure 5). The final breach has an inverted trapezoidal cross section, the depth is 15.6 m, the top width is 9.0 m, and the bottom width is 5.0 m (Figure 6).
3. Numerical Model
From the above test results, it can be clearly seen that surface erosion is the main erosion mode in the breaching process of tailings dam. The development process of the breach is mainly composed of “continuous vertical downcutting and lateral expansion caused by flow erosion” and “intermittent lateral large expansion caused by collapse of the breach slope.” In both tests, the breach increases very rapidly, and the breach flow quickly reaches the peak. However, the investigation of tailings dam accidents [26, 27] show that there are the same situation in actual tailings dam breaching. The breach growth mainly depends on the erodibility of the dam material. Therefore, the special physical and mechanical properties of tailings should be fully considered when establishing a mathematical model for tailings dam breaching.
3.1. Water Balance Relationship
As described above, tailings dams are different from earth-rockfill dams, and they are used for storing discarded tailings. Tailings pumped into upstream reservoir in the form of slurry, so there are not only tailings but also water storing in the reservoir. Due to sedimentation of particles, water is mainly distributed in the upper layer of the reservoir. In the breaching process of the tailings dam, the water level in the upstream reservoir is a dynamic process. The volume change of water in the reservoir is equal to the difference between inflow and outflow, so the water balance can be described bywhere is time, is water level of reservoir, is surface area of reservoir, is inflow, is breach flow, and is flow through dam crest.
Another important factor affecting the change of reservoir water level is the geometric characteristic of reservoir. The reservoir geometric characteristic can be represented by the storage capacity curve. The curve is usually given as pair values of surface area and water level in the reservoir. In some cases, only the reservoir area information or the reservoir capacity information is known, and the storage capacity curve is not available. When this occurs, the storage capacity curve can be calculated by the formula proposed by Wu :where and are coefficients, and the value of generally ranges from 1.0 to 3.0, and is water depth in the reservoir. When the water depth and the corresponding reservoir area or reservoir capacity are known, the coefficients and can be calculated by formula (2). If only the reservoir area or the reservoir capacity is known, is assumed to be 2.0 , then other related parameters can also be calculated by formula (2).
3.2. Breach Flow Routing
Based on the results of above tests, it is known that the breach formed in tailings dam breaching usually has an inverted trapezoidal cross section. Therefore, the breach flow can be calculated by the formula of the broad-crested weir:where is the bottom width of the breach; is the slope ratio of the breach slope (horizontal/vertical); and are both corrected coefficients, refer to the research of Singh , and ; and is submergence corrected coefficient for tailwater effects on weir outflow, which can be calculated by the formula proposed by Fread  and Singh .
3.3. Dam Breaching Process
A simplified method is used to simulate the erosion process of the breach, which assumes the breach having a flat bottom. The erosion at the bottom section of the breach is calculated through the formula proposed by USDA-NRCS :where is erosion rate of the breach bottom; is erosion coefficient; is bed shear stress; and is critical shear stress of the dam material. The erosion coefficient can be measured by tests [32, 33] or calculated using the empirical formula proposed by Temple and Hanson :where is specific weight of water; is dry specific weight of the dam material; and is clay content of the dam material. By considering water content, clay content, and dry density of the tailings, the formula can reasonably reflect the effect of the physical and mechanical properties of tailings on the dam breaching process. can be calculated using Manning’s formula:where is the cross-sectional area of the breach flow; is the hydraulic radius; is the density of the water-tailing mixtures, ; is the volume concentration of tailings in mixtures; and is Manning’s roughness coefficient. The volume change of the breach consists of two parts: the change of the breach bottom elevation and the change of the breach width (Figure 7). It is assumed the relationship between and is determined bywhere is slope angle of the breach and is erosion mode of the breach and for one-sided erosion and for both-sided erosion.
Through the centrifugal model tests of tailings dam breaching, it can be found that the development process of the breach mainly consists of the vertical downcutting process and the lateral expansion process. Besides erosion of the discharge flow, the breach collapsing is another main factor resulting in lateral expansion of the breach. As described above, erosion of the flow is concentrating on lower region of the breach. With the process of erosion, the breach side-slope angle gradually increases. When the angle exceeds the critical value, the sliding force of the soil is greater than that of the sliding resistance. At this point, landslide or collapse of the breach slope makes the breach width increasing suddenly. In this model, limit equilibrium method is used in the analysis of side-slope stability (Figure 8).
In Figure 8, the sliding surface in the breach slope is assumed as a plane, so the angle between the sliding surface and the bottom is α. Then, in order to further simplify the analysis of side-slope stability, the influence of resistance at both ends of the sliding soil mass is neglected. Fd is sliding force, which reflects the component of gravity of the soil on the sliding surface, and can be calculated using the following formula:where is gravity of the tailings on the sliding surface; is specific weight of the tailings; and is sliding resistance consists of frictional resistance and cohesion, which can be calculated using the following formula:where is the internal friction angle of the tailings and is the cohesion of the tailings. When , sliding occurs in the breach slope.
3.4. Algorithm of the Numerical Model
To simulate the discharge flow and breach development at each time step, the iterative method of time steps is used in the tailings dam breaching simulation. Figure 9 shows the flow chart of the numerical method. In Figure 9, input files include initial conditions, boundary conditions, and calculation parameters. After setting the time step , the breach flow , the change of water level of reservoir , the change of the breach bottom elevation , and the change of the breach width in the first time step can be calculated through initial conditions. Therefore, in the next time step, the water level of reservoir is , the breach bottom elevation is , and the breach width is . Then, checking the breach slope stability, if , sliding occurring in the breach slope and calculating , , , and , if not, going to the next time step. Before the calculation of the next time step, checking the values of and , if , starting the calculation; if , ending the calculation and outputting , , , and of current time step. The procedure of the tailings dam breaching is listed as follows (Figure 9).
4. Numerical Model Verification
To verify the rationality of the proposed model, the centrifugal model tests of tailings dam breaching are selected. In two tests, the downstream slope ratios are 1 : 3.0 and 1 : 4.0, other parameters such as properties of the tailings, geometric characteristics of the dam, and reservoir characteristics are shown in Table 2.
The comparison of calculated results and measured data is shown in Table 3, including the peak breach flow Qp, the breach width B and b, and the time to peak flow tp. The comparisons of the calculated breach flow hydrographs and measured ones are shown in Figure 10, and the evolution of the breach widths are shown in Figure 11.
Case 1. The downstream slope ratio is 1 : 3.0. The initial breach was formed at the dam crest at 0.06 h; then the flow begins to discharge along the breach. Since erosion of the flow concentrated on the breach, the breach grows quickly. In addition, the upstream reservoir has a less storage, which results in the breach flow reaching the peak quickly. The peak breach flow calculated using the numerical model is 23.21 m3/s, and the time to peak flow is 0.124 h after dam breaching. At the end of dam breaching, the top width of the final breach is 9.21 m, and the bottom width is 4.83 m. Comparing to the measured data, the peak breach flow calculated using the numerical model has increased by 6.07%, the time to peak flow has decreased by 4.62%, the top width of the final breach has increased by 2.33%, and the bottom width of the final breach has increased by 3.4% (Table 2). From Figures 10(a) and 11(a), it can be seen that the calculated breach flow hydrographs and final breach are almost in accordance with the measured ones.
Case 2. The downstream slope ratio is 1 : 4.0. The initial breach was formed at the dam crest at 0.08 h; then with the fast growth of the breach, the breach flow reaches a peak quickly. The peak breach flow calculated using the numerical model is 19.64 m3/s, and the time to peak flow is 0.149 h after dam breaching. At the end of dam breaching, the top width of the final breach is 9.78 m, and the bottom width is 4.46 m. Comparing to the measured data, the peak breach flow calculated using the numerical model has increased by 2.99%, the time to peak flow has increased by 6.43%, the top width of the final breach has increased by 5.16%, and the bottom width of the final breach has decreased by 3.04% (Table 2). From Figures 10(b) and 11(b), it can be seen that the calculated breach flow hydrographs and final breach are almost in accordance with the measured ones.
Through the analysis of the two cases, it is proved that the proposed model is applicable in modeling the tailings dam breaching due to overtopping.
In order to analyze the influence of tailings properties on the discharge process of breach flow, an earth-rockfill dam breaching model  is selected to simulate the cases mentioned above. The selected model is a typical physical-based model, which cannot fully consider the special physical and mechanical properties of tailings. The comparison of the breach flow hydrographs calculated by two models are shown in Figure 12, and the comparison of the results calculated by two models and measured data is shown in Table 4.
Table 4 shows that the time to peak flow predicted by the selected model in two cases has larger errors, but the peak breach flow and the predicted final breach widths are close to the measured ones. When the tailings properties are not considered, the breach flow increases slowly and the time to the peak flow is much later than that of the proposed model (Figure 12). Therefore, the earth-rockfill dam breaching model cannot simulate the breaching process of the tailings dam rationally.
In this paper, the centrifugal model tests are conducted to study the breaching mechanism of tailings dam. From the model tests, it is found that the breaching process of tailings dam can be divided into four stages: Stage 1, erosion on the dam slope; Stage 2, downcutting of the breach; Stage 3, erosion on the breach slope; and Stage 4, the breach collapsing. Affected by the tailings properties and reservoir characteristics, the breaching process of tailings dam is obviously different from that of earth-rockfill dam. The tailings dam breaching has a short duration, and the breach flow increases rapidly. Based on the centrifugal test results, a mathematical model for tailings dam breaching is established, which can consider the special physical and mechanical properties of tailings. In the model, the vertical downcutting and horizontal expansion are simulated using the erosion rate formula which is derived from shear stress principle of water flow; limit equilibrium method is used to analyze the breach slope stability. In addition, the model adopts iterative calculation method to simulate the whole breaching process. To verify the rationality of the proposed model, the centrifugal model tests of tailings dam breaching are selected. The analysis results show that the peak flow, the time to peak flow, and the width of the breach calculated by the proposed model have small errors; the calculated breach flow hydrographs are almost in accordance with the measured ones. Therefore, the proposed model is applicable in modeling the dam breaching process for tailings dam due to overtopping.
The numerical calculation data used to support the findings of this study are included within the article; the chart data used to support the findings of this study are included within the supplementary information file.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
This work was financially supported by the National Key Research and Development Program of China (grant no. 2017YFC0804605); the National Natural Science Foundation of China (NSFC) (grant no. 51539006); and the Nanjing Hydraulic Research Institute (NHRI) (grant no. Y118013).
The supplementary materials include centrifugal test results and calculation results of numerical model. In the supplementary material file, t is the time, q is the breach flow, b is the bottom width of the breach, and B is the top width of the breach. Measured data are the centrifugal test results, and calculated data are the calculation results of numerical model. (Supplementary Materials)
Q. M. Li, H. Zhang, and G. Li, “Comparative analysis on safety management of tailings ponds in China and Canada,” China Mining Magazine, vol. 26, no. 1, pp. 21–24, 2017, in Chinese.View at: Google Scholar
Anon, “A Review of tailings dam failures,” in International Water Power and Dam Construction, Wilmington Business Publishing, Wilmington, UK, 2001.View at: Google Scholar
S. Azam and Q. R. Li, “Tailing dam failures: a review of the last one hundred years,” Geotechnical News, vol. 28, no. 4, pp. 50–54, 2010.View at: Google Scholar
Could the Hungarian tailings dam tragedy happen in Alberta, http://www.pembina.org, 2010.
A. Michnea and I. Gherheş, “Impact of metals on the environment due to technical accident at Aurul Baia Mare, Romania,” International Journal of Occupational Medicine and Environmental Health, vol. 14, no. 3, pp. 255–259, 2001.View at: Google Scholar
L. Liu, H. W. Zhang, D. Y. Zhong, and R. Z. Miao, “Research on tailings dam break due to overtopping,” Journal of Hydraulic Engineering, vol. 45, no. 6, pp. 675–681, 2014.View at: Google Scholar
L. T. Zhang, “Summary on the dam-break of tailing pond,” Journal of Hydraulic Engineering, vol. 39, no. 5, pp. 594–600, 2013.View at: Google Scholar
The Risk, Public Liability, & Economics of Tailings Storage Facility Failures, https://www.earthworksaction.org/files/pubsothers/BowkerChambers-RiskPublicLiability_EconomicsOfTailingsStorageFacility%20Failures-23Jul15.pdf, 2015.
G. Z. Yin, X. F. Jing, W. Z. Zhang et al., “Experiments of debris flow disaster of tailings dam-break,” Disaster Advances, vol. 4, pp. 4–9, 2011.View at: Google Scholar
X. F. Jing, G. Z. Yin, Z. A. Wei et al., “Study of tailings dam-break surges with floating slurry in model experiment in different collapse gates,” Rock & Soil Mechanics, vol. 33, no. 3, pp. 745–752, 2012.View at: Google Scholar
S. S. Chen, G. M. Xu, Q. M. Zhong, and X. Gu, “Development and application of centrifugal model test system for break of earth-rock dams,” Journal of Hydraulic Engineering, vol. 43, no. 2, pp. 241–245, 2012.View at: Google Scholar
National Energy Administration of China, Technical Specification of Geotechnical Centrifuge Model Test (DL/T 5102—2013), National Energy Administration of China, Beijing, China, 2013, in Chinese.
State Administration of Work Safety, Safety Technical Regulations for the Tailing Pond (AQ2006-2005), State Administration of Work Safety, Beijing, China, 2006, in Chinese.
Z. Z. Wu and G. D. Mei, “Statistical analysis of tailings pond accidents and cause analysis of dam failure,” China Safety Science Journal, vol. 24, no. 9, pp. 70–76, 2014, in Chinese.View at: Google Scholar
V. P. Singh, Dam Breach Modeling Technology, Kluwer Academic Publishers, Dordrecht, Netherlands, 1996.
D. L. Fread, DAMBRK: The NWS Dam Break Flood Forecasting Model, National Weather Service, Silver Spring, MD, USA, 1984.
United States Department of Agriculture, “Natural Resources Conservation Service (USDA-NRCS),” in National Engineering Handbook, NRCS, Washington, DC, USA, 1997.View at: Google Scholar
W. M. Wu, Computational River Dynamics, Taylor & Francis, London, UK, 2007.