##plugins.themes.bootstrap3.article.main##

This is a methodology testing study that uses the case history of Vale Mining’s Brumadinho Dam disaster and data from the Barcelona investigation report. This study aims to compare the results of two constitutive models for tailings: CASM and NorSand, running in the finite element program Plaxis 2D. The same senior authors carried out a previous study on four tailings dams and adopted the NorSand model and parameters obtained from piezocone and laboratory tests. Therefore, this paper aims to check if the same methodology would predict the Brumadinho Dam failure. The results show that both methods can predict Brumadinho failure, and the triggering mechanism suggested in the Barcelona report seems to be the cause of failure.

Downloads

Download data is not yet available.

Introduction

Any major failure can lead to tragedy and an opportunity to learn from past errors and avoid similar future events. This is certainly the case of the Brumadinho Tailings Dam (Figs. 1 and 2) disaster in early 2019, which took nearly 270 lives. This dam failure, also known as Vale’s B1 Córrego do Feijão Dam, together with Samarco’s Fundão Dam failure, was a turning point in dam safety measures [1]. Since then, government dam safety authorities have urged mining companies to increase safety using technical studies, risk classification, employing engineer of records (EoR), improving instrumentation, and monitoring systems, and decommissioning unsafe tailings dams.

Fig. 1. Brumadinho Dam before failure.

Fig. 2. Cross-section of Brumadinho.

The standard GISTM [2] brought up a series of best practices for tailings management. This document recommends the use of a stress-strain approach to analyse the behaviour and safety of existing tailings dams.

Ortigao et al. [3] analysed four upstream-built major tailings dams in late 2019 and spent three years using finite element analysis with the NorSand model [4], [5]. Two of them presented high liquefaction potential (LP), which correlated with low-stress reduction factors (SRF) above two. The other two presented unsafe behaviour with high LPs and lower SRFs.

These studies were based on a series of assumptions ranging from quantity and quality of site investigation to parameter selection to be used with the NorSand constitutive model. These authors also carried out laboratory tests on shallow tailings samples.

This paper presents a methodology test in which Ortigao et al. [3] methods and assumptions were tested at Brumadinho Dam. The specific questions to be answered were:

  • Is NorSand’s constitutive model capable of predicting Brumadinho dam behaviour, like the CASM?
  • Is the characteristic value adopted by Ortigao et al. [3] for extracting the state parameter a valid methodology?

A wealth amount of data exists from Brumadinho failure published in the reports by Robertson et al. [6] and Arroyo et al. [7]. The first one blamed creep, a too-steep dam slope, and high porepressures as the leading cause of failure. On the other hand, Arroyo et al.’s [7] report was able to model the mechanism that triggered liquefaction, referred to here as the Barcelona model. They used the CASM (Clay and Sand Model) [8], [9] and wrote a dll (dynamic link library) for use with the Plaxis 2D software.

The Barcelona report led to different conclusions from the previous one: the triggering mechanism was a deep borehole carried out by Vale’s drill rig (Fig. 3) through the dam slope and soft tailings. The borehole water level was too high, and this led to hydraulic fracture, which, in turn, should have triggered the disaster.

Fig. 3. Vale’s drill-rig operating on the sloping berm of Brumadinho Dam drilling through the embankment and tailings (photo by Mr A Fanaya one week before failure).

In 2022, the Plaxis group significantly improved the software with NorSand model implementation. Therefore, the aim of these studies was:

  • Reproduce the Barcelona model [7], [10] with CASM using the same parameters and the latest Plaxis 2D version.
  • idem, but using NorSand.
  • Analyse liquefaction potential through both models and compare the results.

Clay and Sand Model (CASM)

The Barcelona report used the CASM model for reproducing the tailings’ brittle behaviour. It is not difficult to calibrate model parameters, and it has been successfully used for static liquefaction [10]–[12].

Mánica et al. [11] extended the original CASM to model viscous behaviour, as an attempt to model creep, which Robertson et al. [6] considered one of the major causes of Brumadinho failure.

The constitutive model has been compiled into a dll and implemented in the Plaxis software as a User Defined Soil Model (UDSM). The Barcelona report presents a full programme of validation of the constitutive model.

Fonseca et al. [13] present details of the extensive laboratory testing programme carried out at the University of Porto, Portugal, which enabled CASM parameter extraction.

The rate-independent model (i.e., CASM) requires eight material parameters and the viscoplastic extension requires only two additional ones (Table I).

Parameter Description Range of values
λ The slope of the CSL on the v : ln p’ plot 0.1–0.2 for clays 0.01–0.05 for sand
κ The slope of the unload-reload line on the v : ln p’ plot *
ν Poisson’s ratio 0 < ν ≤ 0.495
ϕ cs Critical state friction angle 15° < ϕcs s ≤ 48°
n Shape control of the yield surface (Fig. 5) n ≥ 1
r Spacing ratio (Fig. 4) r > 1
m The shape of the plastic potential function (Fig. 6) m > 2
Table I. CASM Parameters

Fig. 4. CASM parameters (e, r, G) [8].

Fig. 5. Influence of parameter n on the yield locus [8].

Fig. 6. Normalised plastic potential function for different m values [7].

In addition to the seven compulsory parameters, two of the four parameters below should be input, selected according to the available data. The alternatives are:

  • Γ and eini
  • Γ and po
  • Γ and ξini
  • eini and po
  • eini and ξini
  • po and ξini

where Γ is “altitude” of the CSL, given as the intersection of the CSL with the ordinate for a reference p’ref = 1 kPa, ξini is initial value of the state parameter, R is overconsolidation ratio = p’max/p’o.

There are other optional parameters to consider creep rate effects [7], not discussed in this paper. Table II presents the CASM parameters used in the Barcelona report and also in the following analyses, excluding parameters for rate effects.

Parameter Unit Fine tailings Tailings mix Coarse tailings
λ 0.056 0.05 0.046
κ 0.019 0.019 0.019
n 0.3 0.3 0.3
f cs Degrees 35.7 34.2 33.08
n 3.25 4.9 5.1
r 5.18 12.3 75
m 2.367 2.48 2.465
Γ 2.14 2.05 2.353
R 1.05 1.05 1.05
N 5 5 5
η m2 day kN−1 1.00E−09 4.00E−10 1.00E−11
Table II. CASM Parameters for the Tailings (Arroyo et al. [6])

In addition to the mechanical (stress-strain) formulation, it is necessary to account for the hydraulic behaviour of the tailings in both saturated and unsaturated states.

The saturated hydraulic behaviour of the tailings is simply characterized by a constant value of permeability. Permeability is assumed anisotropic with different values in the vertical and the horizontal directions.

Norsand Model

The NorSand [4] is a critical state model that takes into account the void ratio (Figs. 7 and 8) Table III.

Fig. 7. State parameter in NorSand model.

Fig. 8. Definition of the state parameter (Ψ) [4].

Parameter Typical range Remarks
CSL
Γ 0.9–1.4 The altitude of CSL, defined at 1 kPa
l 0.01–0.07 The slope of CSL, defined on natural logarithm
Plasticity
M tc 1.2–1.5 Critical friction ratio, triaxial compression as the reference condition
N 0.2–0.5 Volumetric coupling coefficient for inelastic stored energy
H 25–500 Plastic hardening modulus for loading, often f(ψ); as a first estimate for refinement, use H =/ λ
χtc 2–5 Relates maximum dilatancy to ψ. Triaxial compression as the reference condition
Elasticity
I r 100–600 Dimensionless shear rigidity (Gmax/p’)
ν 0.1–0.3 Poisson ratio
Table III. NorSand Model Parameters (Jefferies & Been [4])

Hardening in the NorSand model is a function of the current state. The normal consolidation line (NCL) is also a function of the current state. Thanks to these changes to the Cambridge models, NorSand can simulate the behaviour of contractive and loose sands, being suitable to represent strain softening of sandy materials, thus ideal for liquefaction simulation.

The State Parameter Extraction

The in situ piezocone test (CPTU) has been the primary tool for obtaining the initial state parameter for loose sandy tailings [14]. The authors used the same methodology used by Ortigao et al. [3], adopting Robertson’s [15] and a characteristic value corresponding to the first quartile (25% lower).

Table IV presents the NorSand parameters adopted for Brumadinho tailings. These data were derived from the Barcelona report and Fonseca et al. [13].

Tailings
Course Mixed Fine
γunsat (kN/m3) 22 22 22
γsat (kN/m3) 27 27 27
einitial 1.2 1.2 1.2
Gref (kPa) 40000 40000 40000
pref (kPa) 100 100 400
nG 0 0 0
ν 0.3 0.3 0.3
Γ 1.27 1.23 1.2
λe 0.04 0.046 0.053
Mtc 1.4 1.42 1.4
N 0.3 0.3 0.3
χtc 6 6 6
H0 140 130 140
950 950 950
R 1 1 1
S 0 0,25 0
Ψ0 0.05 0.12 0.16
Table IV. NorSand Parameters for Brumadinho Tailings

Liquefaction Potential

The liquefaction potential (LP) is a parameter that indicates the proximity of the CSL (critical state line), also called the flow liquefaction line (Fig. 9). LP is given by the ratio:

Fig. 9. Definition of liquefaction potential (LP) on the stress path.

LP=η/Mwhere η is q/p’.

LP values close to or above 0.7 indicate a high liquefaction potential and the possibility of the stress path turning to the left due to porepressure rise and drastically reducing strength.

Strength Reduction Method

The strength reduction method (SRM) [16], [17] is a well-known technique in which strength parameters are step by step reduced until large displacements take place in the model. The reduction factor is taken as the strength reduction factor (SRF) value. It works well on simple constitutive models like MC and HSM (Plaxis [16]) but has not been implemented in Plaxis for the NorSand model because a more complex model like NorSand may be influenced by other parameters.

For these SRM analyses, the Authors tested the GeoStudio program, which enables SRM calculation.

Brumadinho Geometry and FE Mesh

This investigation uses the same cross-section (Fig. 10) as the Barcelona report and the same materials. Fig. 11 presents the CPTU location on the same cross-section, Fig. 12, the tailings distribution, and Fig. 13 the finite element mesh used for Plaxis 2D analysis.

Fig. 10. Brumadinho Dam analysis cross-section.

Fig. 11. CPTU’s location on the analysed cross-section.

Fig. 12. Distribution of tailings in the cross-section [6].

Fig. 13. Finite element mesh for Plaxis 2D analysis.

Dam Construction Simulation

The dam construction was simulated in 33 stages. The first is just to apply a K0 = 0.5 to the foundation, and the second is to apply initial stresses and porepressures. From stage 3 onwards, all stages, including the dyke construction, followed by tailings filling. Tailings consolidation was applied at each stage, using the MC elastoplastic model with Mohr-Coulomb failure criteria.

The last phase, numbered 33, was when the dam no longer received tailings until failure occurred. At this moment, the tailings model was replaced by the CASM or NorSand models, but using the SSC Soft Soil Creep (Plaxis [18]) for the back of the reservoir, assuming that the elapsed time was not enough for the tailings to change its behaviour.

At this stage, all displacements were set to nil.

Flow Conditions

The flow conditions during construction were simulated by applying a hydraulic head, as recommended by Arroyo et al. [7], as there is not a lot of information during this phase of the structure. At the end of phase 33, the closed flow condition was then replaced by open borders to allow free flow, and the observed water level was applied. The dam was subjected to a precipitation of 1500 mm/year during construction and, after its inactivity, to 1100 mm/year, the same as in the Barcelona report.

Analysis Results

Fig. 14 compares the results of predicted total displacements at the end of constructions, and the maximum results seem to agree, although their location is quite different.

Fig. 14. Total displacements comparison between CASM and NorSand models

Fig. 15 presents the LP obtained through the NorSand model and shows very high values above 0.7 before failure. This is the key parameter to analyse liquefaction. As commented before, Plaxis does not use SRM for the NorSand model on the ground that a complex model, with many parameters, may not yield trustworthy results if only two strength parameters are reduced.

Fig. 15. Liquefaction potential (LP), NorSand model.

Fig. 16 presents the state parameter values at the end of construction, showing a very contractive tailings behaviour.

Fig. 16. State parameter, NorSand model (top) and CASM (bottom).

Fig. 17 shows the plastic points and liquefied points, as predicted by the NorSand model in Plaxis.

Fig. 17. Plastic points, NorSand model.

Liquefaction Triggering Mechanism

Following the Barcelona report, this investigation also analyses the same trigger mechanism by increasing porepressures to 28 m head at the same location (Fig. 18).

Fig. 18. Simulation of water pressure increase of 28 m head in the borehole.

Fig. 19 shows the results for both constitutive models. On the left, it shows the initial conditions just before hydraulic fracturing; on the right the moment the porepressures were increased, leading to high shear strains and liquefaction.

Fig. 19. Resulting shear strains in both CASM and NorSand models due to hydraulic fracturing.

Conclusions

Both CASM and NorSand models predicted the failure mechanisms well, following the footsteps of the Barcelona report.

The LP is the key parameter to predict liquefaction, but only NorSand runs smoothly in Plaxis 2D [18] and yields results of very high LP taking place just before failure.

Both models also agree that the hydraulic failure triggered liquefaction.

Finally, the methodology of selecting a characteristic value for the state parameter seemed to lead to acceptable results.

References

  1. Morgenstern NR, Vick SG, Viotti CB, Watts BD. Fundão Tailings Dam Review Panel: Report on the Immediate Causes of the Failure of the Fundão Dam. New York: Report prepared for Cleary Gottlieb Steen & Hamilton, LLP; 2016. Available from: http://fundaoinvestigation.com/thepanel-report.
     Google Scholar
  2. GISTM. Global Industry Standard on Tailings Management. Global Tailings Review.org; 2020, pp. 39.
     Google Scholar
  3. Ortigao A, Sieira ACC, Moraes J, Santos F, Soares CM. Liquefaction analysis of four high tailings dams, forthcoming publication submitted to soil and rocks. 2023.
     Google Scholar
  4. Jefferies MG. NorSand: a simple critical state model for sand. Géotechnique. 1993;43:91–103.
     Google Scholar
  5. Jefferies M, Been K. Soil Liquefaction: A Critical State Approach. 2nd ed. Spoon; 2016, pp. 690.
     Google Scholar
  6. Robertson PK, Melo L, Williams DJ, & Wilson GW. Panel report on the causes of Córrego do Feijão, (Brumadinho Dam). 2019. Available from: http://www.b1technicalinvestigation.com.
     Google Scholar
  7. Arroyo M, Gens A, Mánica MA, Tarragó R, Monforte L, Lantada N, et al. Computational analyses of Dam I failure at the Córrego de Feijão Mine in Brumadinho Final Report. CIMNE International Centre for Numerical Methods in Engineering, University of Barcelona; 2021.
     Google Scholar
  8. Yu HS. CASM: A unified state parameter model for clay and sand. Int J Numer Anal Methods Geomech. 1998;22(8):621–53.
     Google Scholar
  9. Yu HS, Zhuang PZ, Mo PQ. A unified critical state model for geo-materials with an application to tunnelling. J Rock Mech Geotech Eng. 2019;11(3):464–80.
     Google Scholar
  10. Arroyo M, Gens A. La rotura de la presa de relaves de Brumad- inho: la montaña de cristal. XI Simposio Nacional de Ingeniería Geotécnica Mieres (Asturias), 24 al 27 de mayo de 2022. Español.
     Google Scholar
  11. Mánica M, Arroyo M, Gens A, Monforte L. Application of a critical state model to the Merriespruit tailings dam failure. Proceedings of the Institution of Civil Engineers-Geotechnical Engineering, pp. 1–15, 2021.
     Google Scholar
  12. Gens AS. Hydraulic fills with a special focus on liquefaction. Proceedings of the XVII ECSMGE-2019: Geotechnical Engineering foundation of the future. International Society for Soil Mechanics and Geotechnical Engineering (ISSMGE), pp. 1–31, 2019.
     Google Scholar
  13. Fonseca AV, Cordeiro D, Molina-Gómez F, Besenzon D, Fonseca A, Ferreira C. The mechanics of iron tailings from laboratory tests on reconstituted samples collected in post-mortem Dam I in Brumadinho. Soils and Rocks. 2022;45(2):e2022001122.
     Google Scholar
  14. Robertson PK, Cabal K. Guide do Cone Penetration Testing. 7th ed. Greg Drilling; 2022.
     Google Scholar
  15. Robertson PK, Cabal K. Cone penetration test (CPT)-based soil behaviour type (SBT) classification system—an update. Can Geotech J. 2016;53:1910–27.
     Google Scholar
  16. Dawson EM, Roth WH, Drescher A. Slope stability analysis by strength reduction. Géotechnique. 1999;49(6):835–40.
     Google Scholar
  17. Matsui T, San KC. Finite element slope stability analysis by shear strength reduction technique. Soils Found. 1992;32(1):59–70.
     Google Scholar
  18. Plaxis. Plaxis 2D—tutorial and reference manual. 2022.
     Google Scholar


Most read articles by the same author(s)