Acessibilidade / Reportar erro

A simple method for assessing the probability of liquefaction of tailings dams

Abstract

In recent years, two tragic accidents in Brazilian tailings dams resulted in a lack of trust in the performance of these structures. This paper aims at proposing a new procedure for estimating the liquefaction risk of existing tailings dams, with the use of well-known probabilistic methods as a valuable tool in the decision-making process for triggering an Emergency Plan. The combination of a rigorous limit equilibrium stability method (Spencer, 1967) and the FOSM probabilistic procedure (Christian et al., 1992) has proved to adequately predict the imminent liquefaction risk of both Fundão and Feijão tailings dams, which recently collapsed in Mariana (2015) and in Brumadinho (2019), both in Minas Gerais State, Brazil. Considering the triggering mechanisms known to have occurred, the estimated probabilities of liquefaction failure at the time of the accidents were 36% for Fundão Dam, and 47% for Feijão Dam. The results show that the risk associated to the collapse of both structures was well above the acceptable level suggested by current standards.

Keywords:
Tailings dam; Liquefaction; Probabilistic analysis; Risk assessment

1. Introduction

In recent years, significant accidents have occurred in tailings dams in Brazil, raising doubts about the safety of these structures and lack of trust in their performance (Morgenstern, 2018Morgenstern, N.R. (2018). Geotechnical risk, regulation, and public policy. In The Sixth Victor de Mello Lecture, Salvador, Brazil. Brazilian Society for Soil Mecanics and Geotechnical Engineering.). The catastrophic failure of Fundão tailings dam, which occurred in Mariana, in December 2015, was one of Brazil's major environmental accidents, leaving a serious scar on the mining industry worldwide. A few years later, the world was again shocked by another tragedy, this time captured on the security video cameras showing the abrupt failure of the tailings dam at the Feijão mine, in Brumadinho, in January 2019, and the subsequent massive devastation caused by the uncontrolled flow of the saturated tailings. These two accidents were surely due to static liquefaction of the dams’ tailings.

The accidents raised questions about the methodologies for analyzing dam stability in Brazil, as well as the acceptability of the risk by society.

Concepts of risk analysis can be applied to dams, but they are not yet widely used in Brazil. Espósito & Palmier (2013)Espósito, T., & Palmier, L.R. (2013). Application of risk analysis methods on tailings dams. Soils and Rocks, 36(1), 97-117. http://doi.org/10.28927/SR.361097.
http://doi.org/10.28927/SR.361097...
applied two risk analysis methods to Brazilian dams and concluded that this type of analysis allows a better comprehension of dam behavior. At the same time, the use of probabilistic procedures has proved to be a valuable tool in assessing the safety of these structures. Regulatory guides usually specify a minimum acceptable factor of safety (FS) of 1.50 as the design criteria of dam structures. CDA (2013)CDA. (2013). Dam safety guidelines. Canadian Dam Association. indicates this value as the minimum adequate in static analyses of downstream slopes under steady state (long-term) conditions, as seen in Table 1. Australian National Committee on Large Dams (ANCOLD, 2022ANCOLD. (2022). Guidelines on risk assessment. Australian National Commitee on Large Dams.) presents, beyond the safety factor related to each case, the strength parameters that should be adopted, as displayed in Table 2.

Table 1
Minimum Factor of Safety in static analyses recommended by CDA (2013)CDA. (2013). Dam safety guidelines. Canadian Dam Association..
Table 2
Minimum safety factors recommended by ANCOLD (2022)ANCOLD. (2022). Guidelines on risk assessment. Australian National Commitee on Large Dams..

In recent years, however, the use of probabilistic analysis of dam safety has been growing rapidly, as emphasized by CDA (2013)CDA. (2013). Dam safety guidelines. Canadian Dam Association.. In Brazil, probabilistic analysis has been still limited to academic studies and ignored as a technical tool in the decision-making process of assessing the safety of existing structures.

The aim of this work is to contribute to a better understanding of the risk of failure of existing tailings dams. In this context, this paper presents a simple procedure for estimating the probability of liquefaction failure at the time of the accident, indicating that such analysis could have been of great help in anticipating these tragic occurrences.

2. Liquefaction probabilistic analysis

Liquefaction may classically occur in saturated contractive loose sand deposits. Olson & Stark (2003a)Olson, S.M., & Stark, T.D. (2003a). Yield Strength ratio and liquefaction analysis of slopes and embankments. Journal of Geotechnical and Geoenvironmental Engineering, 129(8), 727-737. http://doi.org/10.1061/(ASCE)1090-0241(2003)129:8(727).
http://doi.org/10.1061/(ASCE)1090-0241(2...
suggests that a liquefaction analysis consists of three main tasks: (1) an analysis of flow failure susceptibility, (2) a triggering analysis, and (3) a post-triggering failure stability analysis.

The first step of a liquefaction analysis of tailings deposits is to investigate if the material is loose and contractive (i.e., susceptible to flow failure). Simple methodologies for assessing the soil's susceptibility to liquefaction have been described by Olson (2001)Olson, S.M. (2001). Liquefaction analysis of level and sloping ground using field case histories and penetration resistance [Doctoral thesis]. Illinois: University of Illinois. and Robertson (2010)Robertson, P.K. (2010). Evaluation of flow liquefaction and liquified strength using the cone penetration test. Journal of Geotechnical and Geoenvironmental Engineering, 136(6), 842-853. http://doi.org/10.1061/(ASCE)GT.1943-5606.0000286.
http://doi.org/10.1061/(ASCE)GT.1943-560...
, with applications to practical cases.

As dam liquefaction usually occurs under undrained conditions, the key factor in computing the probability of liquefaction failure is to estimate the probability of undrained failure. Martin & McRoberts (1999)Martin, T.E., & McRoberts, E.C. (1999). Some considerations in the stability analysis of upstream tailings dams. In Tailings and Mine Waste ’99, Proceedings of the 6th International Conference on Tailings and Mine Waste (pp. 187-302), Fort Collins, Colorado. have proposed that the results of undrained analysis must be understood in the context of the potential triggers of undrained shear. These authors state that the stability margin against the shear failure of tailings dams is better expressed in terms of both the probability of undrained failure and a safety factor from a rigorous limit equilibrium method.

Martin & McRoberts (1999)Martin, T.E., & McRoberts, E.C. (1999). Some considerations in the stability analysis of upstream tailings dams. In Tailings and Mine Waste ’99, Proceedings of the 6th International Conference on Tailings and Mine Waste (pp. 187-302), Fort Collins, Colorado. suggest there are several potential triggers for undrained failure, including erosion, high impoundment rate, fast steepening of the downstream slope, seepage breakout at the dam’s face, foundation and/or embankment movement, sudden reduction of mean effective stress, severe storm runoff, blockage and failure of spillways, earthquake, and seismic deformation.

Olson & Stark (2003a)Olson, S.M., & Stark, T.D. (2003a). Yield Strength ratio and liquefaction analysis of slopes and embankments. Journal of Geotechnical and Geoenvironmental Engineering, 129(8), 727-737. http://doi.org/10.1061/(ASCE)1090-0241(2003)129:8(727).
http://doi.org/10.1061/(ASCE)1090-0241(2...
stated that the strength parameter mostly used in stability analyses is the undrained stress ratio su/σ'v, which is the ratio of undrained strength (su) to initial vertical effective stress. This stress ratio parameter offers a rational way to account for effective stress variations, state parameter (or void ratio), and grain size characteristics (Olson & Stark, 2003bOlson, S.M., & Stark, T.D. (2003b). Use of laboratory data to confirm yield and liquefied strength ratio concepts. Canadian Geotechnical Journal, 40(6), 1164-1184. http://doi.org/10.1139/t03-058.
http://doi.org/10.1139/t03-058...
).

This parameter may be defined in undrained triaxial tests in the laboratory or estimated from field piezocone (CPTu) tests results (Olson & Stark, 2002Olson, S.M., & Stark, T.D. (2002). Liquefied strength ratio from liquefaction flow failure case histories. Canadian Geotechnical Journal, 39(3), 629-647. http://doi.org/10.1139/t02-001.
http://doi.org/10.1139/t02-001...
, 2003aOlson, S.M., & Stark, T.D. (2003a). Yield Strength ratio and liquefaction analysis of slopes and embankments. Journal of Geotechnical and Geoenvironmental Engineering, 129(8), 727-737. http://doi.org/10.1061/(ASCE)1090-0241(2003)129:8(727).
http://doi.org/10.1061/(ASCE)1090-0241(2...
; Sadrekarimi, 2014Sadrekarimi, A. (2014). Effect of the mode of shear on static liquefaction analysis. Journal of Geotechnical and Geoenvironmental Engineering, 140(12), 04014069. http://doi.org/10.1061/(ASCE)GT.1943-5606.0001182.
http://doi.org/10.1061/(ASCE)GT.1943-560...
; Robertson, 2010Robertson, P.K. (2010). Evaluation of flow liquefaction and liquified strength using the cone penetration test. Journal of Geotechnical and Geoenvironmental Engineering, 136(6), 842-853. http://doi.org/10.1061/(ASCE)GT.1943-5606.0000286.
http://doi.org/10.1061/(ASCE)GT.1943-560...
).

After defining the parameters needed for the analysis, the probability of undrained failure can be obtained in several ways. In this paper, Spencer’s limit equilibrium method and the FOSM probabilistic method (First Order Second Moment, described by Christian et al., 1992Christian, J.T., Ladd, C.C., & Baecher, G. (1992). Reliability and probability in stability analysis. Journal of the Geotechnical Engineering Division, 120(2), 1071-1111. .) were adopted, for being practical and suitable in previous applications (Dell’Avanzi & Sayão, 1998Dell’Avanzi, E., & Sayão, A.S. (1998). Evaluation of a slope failure probability. In XI Brazilian Conference on Soil Mecanics and Geotechnical Engineering (COBRAMSEG) (pp. 1289-1295). Brasília: Brazilian Society for Soil Mecanics and Geotechnical Engineering (ABMS) (in Portuguese).; Sayão et al., 2012Sayão, A.S.F.J., Sandroni, S.S., Fontoura, S.A.B., & Ribeiro, R.C.H. (2012). Considerations on the probability of failure of mine slopes. Soils and Rocks, 35(1), 31-37.; Vecci & Sayão, 2019Vecci, A.N., & Sayão, A.S.F.J. (2019). Probabilistic analysis of mine slope stability. In XVI Panamerican Conference on Soil Mechanics and Geotechnical Engineering (pp. 1757-1765). IOS Press.; Yokozawa et al., 2019Yokozawa, S.Y., Assis, A.P., & Rocha, J.S. (2019). Probabilistic analysis applied to the risk assessment of dams. In 11th ICOLD European Club Symposium (pp. 1-12), Chania, Grecia. ICOLD.).

The Spencer limit equilibrium method (Spencer, 1967Spencer, E. (1967). A method of the analysis of the stability of embankments assuming parallel inter-slice forces (pp. 11-26). London: Geotechnique.) was chosen for this analysis because it is a widely used rigorous method. In this scenario, the rupture surface is free and can be either circular or non-circular, satisfying the forces and moments equations of equilibrium. The main hypothesis adopted in this method is the parallel inter-slice forces, which means that the forces between the slices can be replaced by a resultant Qi inclined δ to the horizontal (Spencer, 1967Spencer, E. (1967). A method of the analysis of the stability of embankments assuming parallel inter-slice forces (pp. 11-26). London: Geotechnique.).

Christian et al. (1992)Christian, J.T., Ladd, C.C., & Baecher, G. (1992). Reliability and probability in stability analysis. Journal of the Geotechnical Engineering Division, 120(2), 1071-1111. . and Sayão et al. (2012)Sayão, A.S.F.J., Sandroni, S.S., Fontoura, S.A.B., & Ribeiro, R.C.H. (2012). Considerations on the probability of failure of mine slopes. Soils and Rocks, 35(1), 31-37. indicate that the partial derivatives of the safety factor function in the FOSM method may be obtained by a small variation of each random variable. They suggested that each variable should be distinctly subjected to an increase of 10% on its average value while keeping all other variables fixed. In this way, the resulting variation of the safety factor is adequately assessed. In order to estimate this probability of failure by the FOSM method, the average and standard deviation values of each random variable must be previously obtained, from a proper number of reliable field or lab test results. This requirement may impose a serious difficulty to use probabilistic methods in smaller, low-budget projects.

In the tailings dam’s stability analyses herein presented, the ratio su/σ'v was considered the only random variable, due to its notorious higher relevance.

If the tailings are susceptible to liquefaction, the probability of liquefaction may be given by Equation 1:

P L = P F O S M * P t (1)

where PL is the probability of liquefaction, PFOSM is the probability obtained by the FOSM method, and Pt is the probability of occurrence of a trigger, as indicated by D’Hyppolito (2023)D’Hyppolito, L.C.B.S. (2023). Liquefaction risk analysis in tailings dams [Doctoral thesis, Pontifical Catholic University of Rio de Janeiro]. Pontifical Catholic University of Rio de Janeiro’s repository (in Portuguese)..

After careful analysis of three Brazilian tailings dams, D’Hyppolito (2023)D’Hyppolito, L.C.B.S. (2023). Liquefaction risk analysis in tailings dams [Doctoral thesis, Pontifical Catholic University of Rio de Janeiro]. Pontifical Catholic University of Rio de Janeiro’s repository (in Portuguese). concludes that the probability of trigger for liquefaction in Brazilian tailings dams may be taken 1.0, due to frequent problems, such as insufficient geotechnical investigation, fast tailings disposal, malfunction of drainage systems (Ávila et al., 2021Ávila, J.P., Sawaya, M., Sayão, A. S. F. & Ferreira, L. A. (2021). Safety of tailings dams in brazil: assessment of recent accidents. Geotecnia Journal, (152), 435-464. Special Issue celebrating the 50th Anniversary of the Geotecnia Journal. Geotechnical Societies of Portugal (SPG), Brazil (ABMS) and Spain (SEMSIG). .) and divergences with design parameters.

CDA (2013)CDA. (2013). Dam safety guidelines. Canadian Dam Association. suggests that the probability of liquefaction be obtained from Equation 1 and compared to reference values, which depend on the consequences of a potential dam failure (Table 3). This CDA classification table regulates the necessary safety criteria in terms of design parameters and proposes the risk assessment for dam safety considering the approach shown in Figure 1, with life safety risk guidelines consistent with values used in other hazardous industries and with the principle that risks should be made as low as reasonably practicable (ALARP). In Figure 1, the probability of more than N fatalities may be considered equal to the probability of liquefaction obtained from Equation 1.

Table 3
Classification of consequences of dam failures (CDA, 2013CDA. (2013). Dam safety guidelines. Canadian Dam Association.).
Figure 1
Risk assessment for dams, as a function of number of fatalities (CDA, 2013CDA. (2013). Dam safety guidelines. Canadian Dam Association.).

Recently, the Australian National Committee on Large Dams released a review of the guidelines on risk assessment (ANCOLD, 2022ANCOLD. (2022). Guidelines on risk assessment. Australian National Commitee on Large Dams.), presenting a chart to evaluate the societal risks for existing dams, as shown in the Figure 2. It is interesting to note the similarities between CDA and ANCOLD charts.

Figure 2
Societal risk assessment for existing dams (ANCOLD, 2022ANCOLD. (2022). Guidelines on risk assessment. Australian National Commitee on Large Dams.).

2.1 Fundão dam

On the afternoon of November 5th, 2015, Fundão Dam started to collapse, in Mariana, an important iron mining center in Minas Gerais State, Brazil. This dam was built by the upstream method and was 110 m high, by the time of the accident. This event has been reported in detail by Morgenstern et al. (2016)Morgenstern, N.R., Vick, S.G., Viotti, C.B., & Watts, B.D. (2016). Fundão tailings dam review panel report on the immediate causes of the failure of the Fundão dam. USDA. Retrieved in January 15, 2024, from https://www.resolutionmineeis.us/documents/fundao-2016
https://www.resolutionmineeis.us/documen...
.

Considered as one of the Brazilian’s worst environmental disasters, the failure of Fundão Dam resulted in massive damages and caused 19 deaths. Until January 2020, a total of R$ 8.2 billion (or US$ 1.8 billion) had been already disbursed in repairs and compensations, as reported by Renova (2022)Renova. (2022). Financial statements. Renova Foundation. Retrieved in January 15, 2024, from https://www.fundacaorenova.org/en/the-foundation/
https://www.fundacaorenova.org/en/the-fo...
. The number of deaths would have been much higher if the accident had occurred at night, as most people would have been sleeping in houses devastated by the massive flow of saturated tailings released by the dam’s failure.

The cross-section of Fundão Dam (Figure 3), reported by Morgenstern et al. (2016)Morgenstern, N.R., Vick, S.G., Viotti, C.B., & Watts, B.D. (2016). Fundão tailings dam review panel report on the immediate causes of the failure of the Fundão dam. USDA. Retrieved in January 15, 2024, from https://www.resolutionmineeis.us/documents/fundao-2016
https://www.resolutionmineeis.us/documen...
, has been considered in the analyses reported in this paper. Two types of slurry tailings may be noted, both delivered in separate pipelines to the Fundão impoundment. Sandy tailings, or simply Sands, are a mixture of sand-sized and finer silt particles and are relatively free-draining, but susceptible to liquefaction when loose and saturated. The slimes, on the other hand, are much finer and clay-like in nature - soft and compressible tailings, with low permeability.

Figure 3
Cross section of Fundão dam from Morgenstern et al. (2016)Morgenstern, N.R., Vick, S.G., Viotti, C.B., & Watts, B.D. (2016). Fundão tailings dam review panel report on the immediate causes of the failure of the Fundão dam. USDA. Retrieved in January 15, 2024, from https://www.resolutionmineeis.us/documents/fundao-2016
https://www.resolutionmineeis.us/documen...
.

In Figure 3, a mixture of sand and slimes can be noted, but the exact proportion of this material is not clear. A unique parameter, chosen from field piezocone (CPTu) test results, was adopted to represent the tailings.

Regarding previous studies related to similar tailings dams (Sayão et al., 2012Sayão, A.S.F.J., Sandroni, S.S., Fontoura, S.A.B., & Ribeiro, R.C.H. (2012). Considerations on the probability of failure of mine slopes. Soils and Rocks, 35(1), 31-37.; Braga, 2019Braga, V.P.C. (2019). Probabilistic analysis of a tailings dam stability [MSc’s dissertation, Pontifical Catholic University of Rio de Janeiro]. Pontifical Catholic University of Rio de Janeiro’s repository (in Portuguese). Retrieved in January 15, 2024, from https://www.dbd.puc-rio.br/pergamum/tesesabertas/1712751_2019_completo.pdf
https://www.dbd.puc-rio.br/pergamum/tese...
), it was possible to observe that the variation of the undrained strength ratio is responsible for more than 90% of the failure probability. Additionally, Duncan (2000)Duncan, J.M. (2000). Factors of safety and reliability in geotechnical engineering. Journal of Geotechnical and Geoenvironmental Engineering, 126(4), 307-316. http://doi.org/10.1061/(ASCE)1090-0241(2000)126:4(307).
http://doi.org/10.1061/(ASCE)1090-0241(2...
indicates this parameter presents the most expressive variability of all soil’s strength parameters, with a coefficient of variation reaching about 40%. Consequently, it may be suggested that adopting the undrained strength ratio as the only random variable is justified, considering the simplification of the proposed procedure, without compromising the accuracy of the result.

Regarding the in-situ tests performed on the dam, over 50 different expeditions were reported containing a wide range of tests, including field piezocone test (CPTu), standard penetration test (SPT), and permeability test. Morgenstern et al. (2016)Morgenstern, N.R., Vick, S.G., Viotti, C.B., & Watts, B.D. (2016). Fundão tailings dam review panel report on the immediate causes of the failure of the Fundão dam. USDA. Retrieved in January 15, 2024, from https://www.resolutionmineeis.us/documents/fundao-2016
https://www.resolutionmineeis.us/documen...
report that at that time there was no CPTu or SPT made from the dam’s crest until the bedrock foundation. After the failure, additional tests were performed in the Germano dam, similar to and close to Fundão dam.

The analysis of CPTu data showed that about 70% to 80% of the sandy tailings within 75 m of the dam crest were noted to be contractive, based on the liquefaction criterion proposed by Robertson (2010)Robertson, P.K. (2010). Evaluation of flow liquefaction and liquified strength using the cone penetration test. Journal of Geotechnical and Geoenvironmental Engineering, 136(6), 842-853. http://doi.org/10.1061/(ASCE)GT.1943-5606.0000286.
http://doi.org/10.1061/(ASCE)GT.1943-560...
. At greater distances from 180 m, 95% of the sandy tailings are contractive. This indicates that the hydraulically discharged Fundão sandy tailings satisfied the contractiveness requirement for liquefaction flow sliding.

Additionally, a back analysis was carried out of a failure incident that occurred in September 2014. At that time, there was a local failure of the dam. The strength parameters resulting from this back analysis may be used to compare with the parameters from CPTu tests. From this analysis, Morgenstern et al. (2016)Morgenstern, N.R., Vick, S.G., Viotti, C.B., & Watts, B.D. (2016). Fundão tailings dam review panel report on the immediate causes of the failure of the Fundão dam. USDA. Retrieved in January 15, 2024, from https://www.resolutionmineeis.us/documents/fundao-2016
https://www.resolutionmineeis.us/documen...
obtained a value of 0.22 for su/σ'vpeak.

The average results for each undrained strength ratio in Fundão tailings are shown in Table 4, obtained from the results reported by Morgenstern et al. (2016)Morgenstern, N.R., Vick, S.G., Viotti, C.B., & Watts, B.D. (2016). Fundão tailings dam review panel report on the immediate causes of the failure of the Fundão dam. USDA. Retrieved in January 15, 2024, from https://www.resolutionmineeis.us/documents/fundao-2016
https://www.resolutionmineeis.us/documen...
.

Table 4
Average values of undrained strength ratio (D’Hyppolito, 2023D’Hyppolito, L.C.B.S. (2023). Liquefaction risk analysis in tailings dams [Doctoral thesis, Pontifical Catholic University of Rio de Janeiro]. Pontifical Catholic University of Rio de Janeiro’s repository (in Portuguese).).

The present work adopted as average value the result obtained in the back analysis from the 2014 incident. Since there was a located failure in this incident, it can be considered that the undrained strength was mobilized. Furthermore, this value is similar to the CPTu results, as shown in Table 4.

The standard deviation from the results of CPTu tests may be regarded as evidence of the variability of the material at depth. It is important to consider the fact that various problems were reported in the tests performed in Fundão, which led to doubts about the reliability of the results and their representativeness related to all the material contained in the dam, as presented in Morgenstern et al. (2016)Morgenstern, N.R., Vick, S.G., Viotti, C.B., & Watts, B.D. (2016). Fundão tailings dam review panel report on the immediate causes of the failure of the Fundão dam. USDA. Retrieved in January 15, 2024, from https://www.resolutionmineeis.us/documents/fundao-2016
https://www.resolutionmineeis.us/documen...
. Nevertheless, the analysis of Morgenstern et al. (2016)Morgenstern, N.R., Vick, S.G., Viotti, C.B., & Watts, B.D. (2016). Fundão tailings dam review panel report on the immediate causes of the failure of the Fundão dam. USDA. Retrieved in January 15, 2024, from https://www.resolutionmineeis.us/documents/fundao-2016
https://www.resolutionmineeis.us/documen...
results shows that it can be obtained a standard deviation of 0.03 considering all the CPTu results available in the area. Briefly, the adopted values are: su/σ'vpeak=0.22 e standard deviation =0.03.

2.2 Feijão dam

At 12:28 pm local time on January 25th, 2019, the tailings dam B1 at Córrego do Feijão Iron Ore Mine (Vale, Brazil), in Brumadinho, Brazil, suffered a sudden failure, releasing 12.7 million m3 of tailings, and resulting in a catastrophic mudflow that traveled rapidly downstream. The mudflow destroyed several structures (refectory, offices, homes, and a bridge). This dam failure is unique for having high-quality video images that provide insight into the failure mechanism (Robertson et al., 2019Robertson, P.K., de Melo, L., Williams, D.J., & Wilson, G.W. (2019). Report of the expert panel on the technical causes of the failure of Feijão dam I. Retrieved in January 15, 2024, from http://www.b1technicalinvestigation.com/
http://www.b1technicalinvestigation.com/...
).

At the time of the failure, Feijão dam was 86 m high, with a 720 m long crest. More details of the dam and the causes of its failure were reported by Robertson et al. (2019)Robertson, P.K., de Melo, L., Williams, D.J., & Wilson, G.W. (2019). Report of the expert panel on the technical causes of the failure of Feijão dam I. Retrieved in January 15, 2024, from http://www.b1technicalinvestigation.com/
http://www.b1technicalinvestigation.com/...
.

The main dam tragedy ever seen in Brazil has resulted in 270 deaths, severe damage to the environment and about US$ 8 billion already spent by Vale for social and environmental repair. The cross-section of Feijão dam is presented in Figure 4. This was the highest and the most critical section, as indicated by Robertson et al. (2019)Robertson, P.K., de Melo, L., Williams, D.J., & Wilson, G.W. (2019). Report of the expert panel on the technical causes of the failure of Feijão dam I. Retrieved in January 15, 2024, from http://www.b1technicalinvestigation.com/
http://www.b1technicalinvestigation.com/...
and Polícia Federal (Minas Gerais, 2019).

Figure 4
Critical section from Feijão dam (Robertson et al., 2019Robertson, P.K., de Melo, L., Williams, D.J., & Wilson, G.W. (2019). Report of the expert panel on the technical causes of the failure of Feijão dam I. Retrieved in January 15, 2024, from http://www.b1technicalinvestigation.com/
http://www.b1technicalinvestigation.com/...
).

Like the Fundão case, there were two types of tailings, coarse and fine tailings, and the exact location/proportion of these materials is not clear. This paper adopted a unique parameter to the tailings in the analysis, and the parameter was chosen based on the CPTu tests, which reflect the real situation in the field. In Brumadinho’s case, a larger amount of test data is available, and Robertson et al. (2019)Robertson, P.K., de Melo, L., Williams, D.J., & Wilson, G.W. (2019). Report of the expert panel on the technical causes of the failure of Feijão dam I. Retrieved in January 15, 2024, from http://www.b1technicalinvestigation.com/
http://www.b1technicalinvestigation.com/...
presented a statistical analysis of the CPTu results, as shown in Figure 5. Considering the adoption of a single representative parameter, the option made was to consider the average value equivalent to the weighted average according to the frequency of each type of tailings, and the standard deviation adopted was obtained for the sand tailings due to the greater variability observed in this case. In this scenario, these are the selected parameters: su/σ'vpeak=0.27 e standard deviation =0.12. Additionally, the value adopted is consistent with another analysis described in Braga (2019)Braga, V.P.C. (2019). Probabilistic analysis of a tailings dam stability [MSc’s dissertation, Pontifical Catholic University of Rio de Janeiro]. Pontifical Catholic University of Rio de Janeiro’s repository (in Portuguese). Retrieved in January 15, 2024, from https://www.dbd.puc-rio.br/pergamum/tesesabertas/1712751_2019_completo.pdf
https://www.dbd.puc-rio.br/pergamum/tese...
, Polícia Federal (Minas Gerais, 2019)Minas Gerais. Polícia Federal. (2019). Relatório nº 1.070/2019 - SETEC/SR/PF/MG. Polícia Federal de Minas Gerais (in Portuguese). and TUV SUD (2018)TUV SUD. (2018). Revisão periódica de segurança de barragem - Mina Córrego do Feijão - barragem I - relatório técnico. Brumadinho (in Portuguese)..

Figure 5
Statistical analysis of the peak undrained shear strength ratio (Robertson et al., 2019Robertson, P.K., de Melo, L., Williams, D.J., & Wilson, G.W. (2019). Report of the expert panel on the technical causes of the failure of Feijão dam I. Retrieved in January 15, 2024, from http://www.b1technicalinvestigation.com/
http://www.b1technicalinvestigation.com/...
).

3. Results

3.1 Fundão dam failure

As presented, the average elected for su/σ'vpeakwas 0.22, with 0.05 for standard deviation. The other parameters were adopted equal to Morgenstern et al. (2016)Morgenstern, N.R., Vick, S.G., Viotti, C.B., & Watts, B.D. (2016). Fundão tailings dam review panel report on the immediate causes of the failure of the Fundão dam. USDA. Retrieved in January 15, 2024, from https://www.resolutionmineeis.us/documents/fundao-2016
https://www.resolutionmineeis.us/documen...
because they have a small impact on the results. The stability analysis performed for the average parameters using the Spencer method is presented in Figure 6. Table 5 exhibits the parameters adopted for the analysis.

Figure 6
Stability analysis of Fundão Dam.
Table 5
Parameters adopted in Fundão Dam analysis (D’Hyppolito, 2023D’Hyppolito, L.C.B.S. (2023). Liquefaction risk analysis in tailings dams [Doctoral thesis, Pontifical Catholic University of Rio de Janeiro]. Pontifical Catholic University of Rio de Janeiro’s repository (in Portuguese).).

The combination of the FOSM probabilistic method and the Spencer limit equilibrium method, considering that the random variable has a normal (Gaussian) distribution, resulting in 0.37 for the Reliability Index β (Christian et al., 1992Christian, J.T., Ladd, C.C., & Baecher, G. (1992). Reliability and probability in stability analysis. Journal of the Geotechnical Engineering Division, 120(2), 1071-1111. .) and the probability of failure PFOSM=36%.

The calculated probability is quite high and inconsistent with any kind of engineering construction. Based on the results, the dam was not stable and could be subjected to an undrained collapse because of a trigger.

Regarding the occurrence of triggers, according to Morgenstern et al. (2016)Morgenstern, N.R., Vick, S.G., Viotti, C.B., & Watts, B.D. (2016). Fundão tailings dam review panel report on the immediate causes of the failure of the Fundão dam. USDA. Retrieved in January 15, 2024, from https://www.resolutionmineeis.us/documents/fundao-2016
https://www.resolutionmineeis.us/documen...
, the dam was showing (i) serious issues related to drainage; (ii) disobedience to the width of the beach designed, leading to a fine waste settlement in unforeseen places; (iii) small displacements through the instrumentation; and (iv) quick loading to heighten the structure. Those problems are known for their possibility to trigger the tailings liquefaction, as already commented.

Considering that the tailings were subjected to liquefaction and the triggers occurred, the Equation 1 shows that: Pliquefaction=36%. The failure consequences can be classified as Extreme, according to Table 2 (CDA 2013CDA. (2013). Dam safety guidelines. Canadian Dam Association.), due to the resultant environmental impact, which was expressive and unreturnable.

The risk associated to the Fundão dam was way beyond the acceptable. In this scenario, measures should have been taken to mitigate this risk through probability reduction and consequence control.

3.2 Feijão dam failure

The stability analysis performed by the Spencer method with average parameters is presented in Figure 7. Table 6 exhibits the parameters adopted in the analysis. The average selected for su/σ'vpeak was 0.27 and 0.12 for standard deviation. These other parameters were adopted as indicated by Morgenstern et al. (2016)Morgenstern, N.R., Vick, S.G., Viotti, C.B., & Watts, B.D. (2016). Fundão tailings dam review panel report on the immediate causes of the failure of the Fundão dam. USDA. Retrieved in January 15, 2024, from https://www.resolutionmineeis.us/documents/fundao-2016
https://www.resolutionmineeis.us/documen...
because they have a small impact on the results.

Figure 7
Feijão dam stability analysis with average parameters, obtaining a safety factor of 1.03.
Table 6
Parameters adopted in Feijão dam analysis (D’Hyppolito, 2023D’Hyppolito, L.C.B.S. (2023). Liquefaction risk analysis in tailings dams [Doctoral thesis, Pontifical Catholic University of Rio de Janeiro]. Pontifical Catholic University of Rio de Janeiro’s repository (in Portuguese).).

The combination of the FOSM probabilistic method and the Spencer limit equilibrium method, considering that the random variable has a normal (Gaussian) distribution, achieves 0.07 for the reliability index β (Christian et al., 1992Christian, J.T., Ladd, C.C., & Baecher, G. (1992). Reliability and probability in stability analysis. Journal of the Geotechnical Engineering Division, 120(2), 1071-1111. .) and the failure probability is 47%.

As presented by Robertson et al. (2019), (i)Robertson, P.K., de Melo, L., Williams, D.J., & Wilson, G.W. (2019). Report of the expert panel on the technical causes of the failure of Feijão dam I. Retrieved in January 15, 2024, from http://www.b1technicalinvestigation.com/
http://www.b1technicalinvestigation.com/...
the dam was presenting significant issues related to drainage; (ii) there were layers of fine tailings in areas supposed to have tailings with larger particle size; (iii) the foundation data and information were very limited; and (iv) slow displacements were observed, consistent with a possible creep, through one of the instruments (InSar). In this scenario, the liquefaction triggers were present.

Considering the tailings were subjected to liquefaction and the triggers occurred, the use of Equation 1 results in Pliquefaction=47%. The failure consequences can be classified as Extreme, according to Table 2 (CDA, 2013CDA. (2013). Dam safety guidelines. Canadian Dam Association.), due to the associated numbers of deaths.

The risk associated to the Feijão dam was not satisfactory. In this scenario, actions should have been taken to mitigate this risk by reducing the probability of failure and controlling the consequences.

4. Discussion

Brazilian dams have specific problems demanding special caution, such as drainage systems (Ávila et al., 2021Ávila, J.P., Sawaya, M., Sayão, A. S. F. & Ferreira, L. A. (2021). Safety of tailings dams in brazil: assessment of recent accidents. Geotecnia Journal, (152), 435-464. Special Issue celebrating the 50th Anniversary of the Geotecnia Journal. Geotechnical Societies of Portugal (SPG), Brazil (ABMS) and Spain (SEMSIG). .), quick loading, and lack of attention to the design conception elements, like the beach width (Martin et al., 2002Martin, T.E., McRoberts, E.C., & Davies, M.P. (2002). A tale of four upstream tailings dams. In Proceedings of Tailings Dams 2002, Las Vegas, Nevada, USA. Association of State Dam Safety Officials.). These problems may be seen frequently in tailing dams and open the chances for the occurrence of triggers.

Considering the information presented in this study, the most important element in the liquefaction risk analysis is the undrained probability of failure. The estimated probabilities in the two cases reviewed herein were -36% at the time of the Fundão dam (Mariana) failure, which collapsed in 2016, and 47% at the time of the Feijão dam (Brumadinho) failure, which collapsed in 2019.

These two cases discussed in this paper, Fundão and Feijão dams could be considered to present very high levels of associated consequences. For the Feijão dam, which caused over 200 deaths, the acceptance limit as suggested by CDA (2013)CDA. (2013). Dam safety guidelines. Canadian Dam Association. would be 10-7 and tolerability of 10-5. In Fundão dam, in which much lower fatality level (19 deaths), the acceptance limit as suggested by CDA (2013)CDA. (2013). Dam safety guidelines. Canadian Dam Association. would be 10-6 and tolerability of 10-4. The limits of tolerability proposed by ANCOLD (2022)ANCOLD. (2022). Guidelines on risk assessment. Australian National Commitee on Large Dams. are similar. However, it was not the worst case. As previously mentioned, the number of deaths in the Fundão case could have been much higher, and similar to the Brumadinho case, had the accident occurred at night. The computed results for liquefaction failure probability were much higher than these reference values and would be unacceptable for any kind of engineering construction.

5. Conclusions

This paper presents a simple procedure for estimating the risk associated to tailings dams, using simple well-known methods. The practical application of the proposed procedure was illustrated by two recent cases of dam failure, suggesting that the risk associated with the collapsed structures was higher than the acceptable level indicated by current standards. The results of the analysis presented in this paper indicate that in the case of the Fundão Dam, the estimated probability of failure was 36%, while in the case of the Feijão Dam, it was 47%. In both scenarios, the probabilities suggest a level of risk deemed intolerable for dam structures.

The probability of failure may be estimated with the procedure suggested in this paper and may be considered in the decision-making procedures, regarding the need for intervention and mitigation of the risk of accident and its consequences by triggering an Emergency Plan.

List of symbols and abbreviations

ALARP As low as reasonably practicable

ANCOLD Australian National Committee on Large Dams

CNEN Brazilian Commission of Nuclear Energy

CDA Canadian Dam Association

CPTu Field Piezocone Test

FS Factor of safety

FOSM First-Order Second Moment

PL Probability of liquefaction

PFOSM Probability obtained by the FOSM method

Pt Probability of occurrence of a trigger

PUC Pontifical Catholic University of Rio de Janeiro

su/σ'v Undrained stress ratio

su/σ'vpeak Peak undrained stress ratio

β Reliability index

ϕ Friction angle

Data availability

The datasets generated and analyzed in the course of the current study are available from the corresponding author upon request.

Acknowledgements

The authors thank PUC-Rio and the Brazilian Commission of Nuclear Energy (CNEN) for the financial support during these studies. Contributions from Dr Sandro Sandroni on earlier research on the use of probabilistic procedures at PUC-Rio are also appreciated.

  • Discussion open until February 28, 2025.

References

  • ANCOLD. (2022). Guidelines on risk assessment. Australian National Commitee on Large Dams.
  • Ávila, J.P., Sawaya, M., Sayão, A. S. F. & Ferreira, L. A. (2021). Safety of tailings dams in brazil: assessment of recent accidents. Geotecnia Journal, (152), 435-464. Special Issue celebrating the 50th Anniversary of the Geotecnia Journal. Geotechnical Societies of Portugal (SPG), Brazil (ABMS) and Spain (SEMSIG). .
  • Braga, V.P.C. (2019). Probabilistic analysis of a tailings dam stability [MSc’s dissertation, Pontifical Catholic University of Rio de Janeiro]. Pontifical Catholic University of Rio de Janeiro’s repository (in Portuguese). Retrieved in January 15, 2024, from https://www.dbd.puc-rio.br/pergamum/tesesabertas/1712751_2019_completo.pdf
    » https://www.dbd.puc-rio.br/pergamum/tesesabertas/1712751_2019_completo.pdf
  • CDA. (2013). Dam safety guidelines. Canadian Dam Association.
  • Christian, J.T., Ladd, C.C., & Baecher, G. (1992). Reliability and probability in stability analysis. Journal of the Geotechnical Engineering Division, 120(2), 1071-1111. .
  • Dell’Avanzi, E., & Sayão, A.S. (1998). Evaluation of a slope failure probability. In XI Brazilian Conference on Soil Mecanics and Geotechnical Engineering (COBRAMSEG) (pp. 1289-1295). Brasília: Brazilian Society for Soil Mecanics and Geotechnical Engineering (ABMS) (in Portuguese).
  • D’Hyppolito, L.C.B.S. (2023). Liquefaction risk analysis in tailings dams [Doctoral thesis, Pontifical Catholic University of Rio de Janeiro]. Pontifical Catholic University of Rio de Janeiro’s repository (in Portuguese).
  • Duncan, J.M. (2000). Factors of safety and reliability in geotechnical engineering. Journal of Geotechnical and Geoenvironmental Engineering, 126(4), 307-316. http://doi.org/10.1061/(ASCE)1090-0241(2000)126:4(307)
    » http://doi.org/10.1061/(ASCE)1090-0241(2000)126:4(307)
  • Espósito, T., & Palmier, L.R. (2013). Application of risk analysis methods on tailings dams. Soils and Rocks, 36(1), 97-117. http://doi.org/10.28927/SR.361097
    » http://doi.org/10.28927/SR.361097
  • Martin, T.E., & McRoberts, E.C. (1999). Some considerations in the stability analysis of upstream tailings dams. In Tailings and Mine Waste ’99, Proceedings of the 6th International Conference on Tailings and Mine Waste (pp. 187-302), Fort Collins, Colorado.
  • Martin, T.E., McRoberts, E.C., & Davies, M.P. (2002). A tale of four upstream tailings dams. In Proceedings of Tailings Dams 2002, Las Vegas, Nevada, USA. Association of State Dam Safety Officials.
  • Minas Gerais. Polícia Federal. (2019). Relatório nº 1.070/2019 - SETEC/SR/PF/MG Polícia Federal de Minas Gerais (in Portuguese).
  • Morgenstern, N.R. (2018). Geotechnical risk, regulation, and public policy. In The Sixth Victor de Mello Lecture, Salvador, Brazil. Brazilian Society for Soil Mecanics and Geotechnical Engineering.
  • Morgenstern, N.R., Vick, S.G., Viotti, C.B., & Watts, B.D. (2016). Fundão tailings dam review panel report on the immediate causes of the failure of the Fundão dam. USDA. Retrieved in January 15, 2024, from https://www.resolutionmineeis.us/documents/fundao-2016
    » https://www.resolutionmineeis.us/documents/fundao-2016
  • Olson, S.M., & Stark, T.D. (2002). Liquefied strength ratio from liquefaction flow failure case histories. Canadian Geotechnical Journal, 39(3), 629-647. http://doi.org/10.1139/t02-001
    » http://doi.org/10.1139/t02-001
  • Olson, S.M., & Stark, T.D. (2003a). Yield Strength ratio and liquefaction analysis of slopes and embankments. Journal of Geotechnical and Geoenvironmental Engineering, 129(8), 727-737. http://doi.org/10.1061/(ASCE)1090-0241(2003)129:8(727)
    » http://doi.org/10.1061/(ASCE)1090-0241(2003)129:8(727)
  • Olson, S.M., & Stark, T.D. (2003b). Use of laboratory data to confirm yield and liquefied strength ratio concepts. Canadian Geotechnical Journal, 40(6), 1164-1184. http://doi.org/10.1139/t03-058
    » http://doi.org/10.1139/t03-058
  • Olson, S.M. (2001). Liquefaction analysis of level and sloping ground using field case histories and penetration resistance [Doctoral thesis]. Illinois: University of Illinois.
  • Renova. (2022). Financial statements. Renova Foundation. Retrieved in January 15, 2024, from https://www.fundacaorenova.org/en/the-foundation/
    » https://www.fundacaorenova.org/en/the-foundation/
  • Robertson, P.K. (2010). Evaluation of flow liquefaction and liquified strength using the cone penetration test. Journal of Geotechnical and Geoenvironmental Engineering, 136(6), 842-853. http://doi.org/10.1061/(ASCE)GT.1943-5606.0000286
    » http://doi.org/10.1061/(ASCE)GT.1943-5606.0000286
  • Robertson, P.K., de Melo, L., Williams, D.J., & Wilson, G.W. (2019). Report of the expert panel on the technical causes of the failure of Feijão dam I. Retrieved in January 15, 2024, from http://www.b1technicalinvestigation.com/
    » http://www.b1technicalinvestigation.com/
  • Sadrekarimi, A. (2014). Effect of the mode of shear on static liquefaction analysis. Journal of Geotechnical and Geoenvironmental Engineering, 140(12), 04014069. http://doi.org/10.1061/(ASCE)GT.1943-5606.0001182
    » http://doi.org/10.1061/(ASCE)GT.1943-5606.0001182
  • Sayão, A.S.F.J., Sandroni, S.S., Fontoura, S.A.B., & Ribeiro, R.C.H. (2012). Considerations on the probability of failure of mine slopes. Soils and Rocks, 35(1), 31-37.
  • Spencer, E. (1967). A method of the analysis of the stability of embankments assuming parallel inter-slice forces (pp. 11-26). London: Geotechnique.
  • TUV SUD. (2018). Revisão periódica de segurança de barragem - Mina Córrego do Feijão - barragem I - relatório técnico. Brumadinho (in Portuguese).
  • Vecci, A.N., & Sayão, A.S.F.J. (2019). Probabilistic analysis of mine slope stability. In XVI Panamerican Conference on Soil Mechanics and Geotechnical Engineering (pp. 1757-1765). IOS Press.
  • Yokozawa, S.Y., Assis, A.P., & Rocha, J.S. (2019). Probabilistic analysis applied to the risk assessment of dams. In 11th ICOLD European Club Symposium (pp. 1-12), Chania, Grecia. ICOLD.

Publication Dates

  • Publication in this collection
    16 Aug 2024
  • Date of issue
    2024

History

  • Received
    15 Jan 2024
  • Accepted
    14 May 2024
Associação Brasileira de Mecânica dos Solos Av. Queiroz Filho, 1700 - Torre A, Sala 106, Cep: 05319-000, São Paulo - SP - Brasil, Tel: (11) 3833-0023 - São Paulo - SP - Brazil
E-mail: secretariat@soilsandrocks.com