Stability analysis of dam failures. Application to aurul tailing pond in Baia Mare (Romania)

La gestión de estériles de una explotación minera es un punto clave en el desarrollo económico de una actividad extractiva, y en especial, del entorno natural y social en el que se emplaza dicho proyecto. La minería de metales preciosos lleva asociada la construcción de balsas de residuos muy peligrosos, fruto de su proceso extractivo, como por ejemplo la cianuración en el caso del oro. Para un correcto funcionamiento de dichos emplazamientos es necesario escoger correctamente el método constructivo a partir de estudios de reconocimiento previos, como estudios de estabilidad geotécnica, contexto geológico de la zona, sismicidad, hidrología, etc. Así mismo, han de llevarse a cabo unas exhaustivas medidas de control y vigilancia para asegurar las condiciones de seguridad exigidas. La ruptura de la balsa de decantación de Aurul S.A. en Baia Mare (Rumania) el 30 de Enero del año 2000 ha sido escogido como caso de estudio de estabilidad de diques. ABSTRACT Tailing's management of a mining exploitation is a key point in the economical development of the extractive activity and, especially, of the natural and social environment of the site. Precious metals mining has high hazardous embankment construction associated, product of its extractive process, i.e. gold cyanidation. A correct operation of those sites makes necessary to choose a suitable construction method, based on previous studies as geotechnical stability studies, geological context of the area, seismicity, hydrology, etc. At the same time, exhaustive control and monitoring must be carried out in order to assure the required safety conditions. Aurul's decantation pond failure in Baia Mare (Romania), on 30th January 2000, has been chosen as a stability analysis case-study.

Engineering closure of an open pit gold operation in a semi-arid climate

Along with material characteristics and geometry, the climate in which a mine is located can have a dramatic effect on the appropriate options for rehabilitation. The paper outlines the setting, mining, milling and waste disposal at Kidston Gold Mine's open pit operations in the semi-arid climate of North Queensland, Australia, before focusing on the engineering aspects of the rehabilitation of Kidston. The mine took a holistic and proactive approach to rehabilitation, and was prepared to demonstrate a number of innovative approaches, which are described in the paper. Engineering issues that had to be addressed included the geotechnical stability and deformation of waste rock dumps, including a 240 m high in-pit dump: the construction and performance monitoring of a “store and release” cover over potentially acid forming mineralised waste rock; erosion from the side slopes of the waste rock dumps; the in-pit co-disposal of waste rock and thickened tailings; the geotechnical stability of the tailings dam wall; the potential for erosion of bare tailings; the water balance of the tailings dam; direct revegetation of the tailings; and the pit hydrology. The rehabilitation of the mine represents an important benchmark in mine site rehabilitation best practice, from which lessons applicable worldwide can be shared.

Beurteilung der geomechanischen Stabilität und Integrität von Endlagerbergwerken im Salzgebirge auf der Grundlage geologischer und ingenieurgeologischer Untersuchungen

Any safety assessment of a permanent repository for radioactive waste has to include an analysis of the geomechanical stability of the repository and integrity of the geological barrier. Such an analysis is based on geological and engineering geological studies of the site, on laboratory and in-situ experiments, and on numerical calculations. Central part of the safety analysis is the geomechanical modelling of the host rock. The model should simulate as closely as possible the conditions at the site and the behaviour of the rock (e.g., geology, repository geometry, initial rock stress, and constitutive models). On the basis of the geomechanical model numerical calculations are carried out using the finite-element method and an appropriate discretization of the repository and the host rock. The assessment of the repository stability and the barrier integrity is based on calculated stress and deformation and on the behaviour of the host rock measured and observed in situ. An example of the geomechanical analysis of the stability and integrity of the Bartensieben mine, a former salt mine, is presented. This mine is actually used as a repository for low level radioactive waste. The example includes all necessary steps of geological, engineering geological, and geotechnical investigations.

Evaluation of the geotechnical properties of MSW in two Brazilian landfills

The characteristics of municipal solid waste (MSW) play a key role in many aspects of waste disposal facilities and landfills. Because most of a landfill is made up of MSW, the overall stability of the landfill slopes are governed by the strength parameters and physical properties of the MSW. These parameters are also important in interactions involving the waste body and the landfill structures: cover liner, leachate and gas collection systems. On the other hand, the composition of the waste, which affects the geotechnical behavior of the MSW, is dependent on a variety of factors such as climate, disposal technology, the culture and habits of the local community. It is therefore essential that the design and stability evaluations of landfills in each region be performed based on the local conditions and the geotechnical characteristic of the MSW. The Bandeirantes Landfill, BL, in Sao Paulo and the Metropolitan Center Landfill, MCL, in Salvador, are among the biggest landfills in Brazil. These two disposal facilities have been used for the development of research involving waste mechanics in recent years. Considerable work has been made in the laboratory and in the field to evaluate parameters such as water and organic contents, composition, permeability, and shear strength. This paper shows and analyzes the results of tests performed on these two landfills. The authors believe that these results could be a good reference for certain aspects and geotechnical properties of MSW materials in countries with similar conditions. (C) 2010 Elsevier Ltd. All rights reserved.

Unsaturated Infinite Slope Stability Considering Surface Flux Conditions

A slope stability model is derived for an infinite slope subjected to unsaturated infiltration flow above a phreatic surface. Closed form steady state solutions are derived for the matric suction and degree of saturation profiles. Soil unit weight, consistent with the degree of saturation profile, is also directly calculated and introduced into the analyzes, resulting in closed-form solutions for typical soil parameters and an infinite series solution for arbitrary soil parameters. The solutions are coupled with the infinite slope stability equations to establish a fully realized safety factor function. In general, consideration of soil suction results in higher factor of safety. The increase in shear strength due to the inclusion of soil suction is analogous to making an addition to the cohesion, which, of course, increases the factor of safety against sliding. However, for cohesive soils, the results show lower safety factors for slip surfaces approaching the phreatic surface compared to those produced by common safety factor calculations. The lower factor of safety is due to the increased soil unit weight considered in the matric suction model but not usually accounted for in practice wherein the soil is treated as dry above the phreatic surface. The developed model is verified with a published case study, correctly predicting stability under dry conditions and correctly predicting failure for a particular storm.

Review of Stability Berm Alternatives for Environmentally Senstitive Areas, 2005

Stability berms are commonly constructed where roadway embankments cross soft or unstable ground conditions. Under certain circumstances, the construction of stability berms cause unfavorable environmental impacts, either directly or indirectly, through their effect on wetlands, endangered species habitat, stream channelization, longer culvert lengths, larger right-of-way purchases, and construction access limits. Due to an ever more restrictive regulatory environment, these impacts are problematic. The result is the loss of valuable natural resources to the public, lengthy permitting review processes for the department of transportation and permitting agencies, and the additional expenditures of time and money for all parties. The purpose of this project was to review existing stability berm alternatives for potential use in environmentally sensitive areas. The project also evaluates how stabilization technologies are made feasible, desirable, and cost-effective for transportation projects and determines which alternatives afford practical solutions for avoiding and minimizing impacts to environmentally sensitive areas. An online survey of engineers at state departments of transportation was also conducted to assess the frequency and cost effectiveness of the various stabilization technologies. Geotechnical engineers that responded to the survey overwhelmingly use geosynthetic reinforcement as a suitable and cost-effective solution for stabilizing embankments and cut slopes. Alternatively, chemical stabilization and installation of lime/cement columns is rarely a remediation measure employed by state departments of transportation.

Innovative Solutions for Slope Stability Reinforcement and Characterization: Vol. I, TR-489, 2005

Soil slope instability concerning highway infrastructure is an ongoing problem in Iowa, as slope failures endanger public safety and continue to result in costly repair work. Characterization of slope failures is complicated, because the factors affecting slope stability can be difficult to discern and measure, particularly soil shear strength parameters. While in the past extensive research has been conducted on slope stability investigations and analysis, this research consists of field investigations addressing both the characterization and reinforcement of such slope failures. The current research focuses on applying an infrequently-used testing technique comprised of the Borehole Shear Test (BST). This in-situ test rapidly provides effective (i.e., drained) shear strength parameter values of soil. Using the BST device, fifteen Iowa slopes (fourteen failures and one proposed slope) were investigated and documented. Particular attention was paid to highly weathered shale and glacial till soil deposits, which have both been associated with slope failures in the southern Iowa drift region. Conventional laboratory tests including direct shear tests, triaxial compression tests, and ring shear tests were also performed on undisturbed and reconstituted soil samples to supplement BST results. The shear strength measurements were incorporated into complete evaluations of slope stability using both limit equilibrium and probabilistic analyses. The research methods and findings of these investigations are summarized in Volume 1 of this report. Research details of the independent characterization and reinforcement investigations are provided in Volumes 2 and 3, respectively. Combined, the field investigations offer guidance on identifying the factors that affect slope stability at a particular location and also on designing slope reinforcement using pile elements for cases where remedial measures are necessary. The research findings are expected to benefit civil and geotechnical engineers of government transportation agencies, consultants, and contractors dealing with slope stability, slope remediation, and geotechnical testing in Iowa.

A New Prediction Model for Slope Stability Analysis

The instability of river bank can result in considerable human and land losses. The Po river is the most important in Italy, characterized by main banks of significant and constantly increasing height. This study presents multilayer perceptron of artificial neural network (ANN) to construct prediction models for the stability analysis of river banks along the Po River, under various river and groundwater boundary conditions. For this aim, a number of networks of threshold logic unit are tested using different combinations of the input parameters. Factor of safety (FS), as an index of slope stability, is formulated in terms of several influencing geometrical and geotechnical parameters. In order to obtain a comprehensive geotechnical database, several cone penetration tests from the study site have been interpreted. The proposed models are developed upon stability analyses using finite element code over different representative sections of river embankments. For the validity verification, the ANN models are employed to predict the FS values of a part of the database beyond the calibration data domain. The results indicate that the proposed ANN models are effective tools for evaluating the slope stability. The ANN models notably outperform the derived multiple linear regression models.

Statistical analysis of the error associated with the simplification of the stratigraphy in geotechnical models

The uncertainties in the determination of the stratigraphic profile of natural soils is one of the main problems in geotechnics, in particular for landslide characterization and modeling. The study deals with a new approach in geotechnical modeling which relays on a stochastic generation of different soil layers distributions, following a boolean logic – the method has been thus called BoSG (Boolean Stochastic Generation). In this way, it is possible to randomize the presence of a specific material interdigitated in a uniform matrix. In the building of a geotechnical model it is generally common to discard some stratigraphic data in order to simplify the model itself, assuming that the significance of the results of the modeling procedure would not be affected. With the proposed technique it is possible to quantify the error associated with this simplification. Moreover, it could be used to determine the most significant zones where eventual further investigations and surveys would be more effective to build the geotechnical model of the slope. The commercial software FLAC was used for the 2D and 3D geotechnical model. The distribution of the materials was randomized through a specifically coded MatLab program that automatically generates text files, each of them representing a specific soil configuration. Besides, a routine was designed to automate the computation of FLAC with the different data files in order to maximize the sample number. The methodology is applied with reference to a simplified slope in 2D, a simplified slope in 3D and an actual landslide, namely the Mortisa mudslide (Cortina d’Ampezzo, BL, Italy). However, it could be extended to numerous different cases, especially for hydrogeological analysis and landslide stability assessment, in different geological and geomorphological contexts.

Slope stability analysis of the Pacaya Volcano, Guatemala, using limit equilibrium and finite element method

The Pacaya volcanic complex is part of the Central American volcanic arc, which is associated with the subduction of the Cocos tectonic plate under the Caribbean plate. Located 30 km south of Guatemala City, Pacaya is situated on the southern rim of the Amatitlan Caldera. It is the largest post-caldera volcano, and has been one of Central America’s most active volcanoes over the last 500 years. Between 400 and 2000 years B.P, the Pacaya volcano had experienced a huge collapse, which resulted in the formation of horseshoe-shaped scarp that is still visible. In the recent years, several smaller collapses have been associated with the activity of the volcano (in 1961 and 2010) affecting its northwestern flanks, which are likely to be induced by the local and regional stress changes. The similar orientation of dry and volcanic fissures and the distribution of new vents would likely explain the reactivation of the pre-existing stress configuration responsible for the old-collapse. This paper presents the first stability analysis of the Pacaya volcanic flank. The inputs for the geological and geotechnical models were defined based on the stratigraphical, lithological, structural data, and material properties obtained from field survey and lab tests. According to the mechanical characteristics, three lithotechnical units were defined: Lava, Lava-Breccia and Breccia-Lava. The Hoek and Brown’s failure criterion was applied for each lithotechnical unit and the rock mass friction angle, apparent cohesion, and strength and deformation characteristics were computed in a specified stress range. Further, the stability of the volcano was evaluated by two-dimensional analysis performed by Limit Equilibrium (LEM, ROCSCIENCE) and Finite Element Method (FEM, PHASE 2 7.0). The stability analysis mainly focused on the modern Pacaya volcano built inside the collapse amphitheatre of “Old Pacaya”. The volcanic instability was assessed based on the variability of safety factor using deterministic, sensitivity, and probabilistic analysis considering the gravitational instability and the effects of external forces such as magma pressure and seismicity as potential triggering mechanisms of lateral collapse. The preliminary results from the analysis provide two insights: first, the least stable sector is on the south-western flank of the volcano; second, the lowest safety factor value suggests that the edifice is stable under gravity alone, and the external triggering mechanism can represent a likely destabilizing factor.

Slope stability analyses based on sediment cores of the Ligurian Margin, Southern France

Submarine slope failures of various types and sizes are common along the tectonic and seismically active Ligurian margin, northwestern Mediterranean Sea, primarily because of seismicity up to ~M6, rapid sediment deposition in the Var fluvial system, and steepness of the continental slope (average 11°). We present geophysical, sedimentological and geotechnical results of two distinct slides in water depth >1,500 m: one located on the flank of the Upper Var Valley called Western Slide (WS), another located at the base of continental slope called Eastern Slide (ES). WS is a superficial slide characterized by a slope angle of ~4.6° and shallow scar (~30 m) whereas ES is a deep-seated slide with a lower slope angle (~3°) and deep scar (~100 m). Both areas mainly comprise clayey silt with intermediate plasticity, low water content (30-75 %) and underconsolidation to strong overconsolidation. Upslope undeformed sediments have low undrained shear strength (0-20 kPa) increasing gradually with depth, whereas an abrupt increase in strength up to 200 kPa occurs at a depth of ~3.6 m in the headwall of WS and ~1.0 m in the headwall of ES. These boundaries are interpreted as earlier failure planes that have been covered by hemipelagite or talus from upslope after landslide emplacement. Infinite slope stability analyses indicate both sites are stable under static conditions; however, slope failure may occur in undrained earthquake condition. Peak earthquake acceleration from 0.09 g on WS and 0.12 g on ES, i.e. M5-5.3 earthquakes on the spot, would be required to induce slope instability. Different failure styles include rapid sedimentation on steep canyon flanks with undercutting causing superficial slides in the west and an earthquake on the adjacent Marcel fault to trigger a deep-seated slide in the east.

General Report on Main Session 4: Embankments and dams, slope stability and landslides

The General Reporter presents the papers from the Authors, along with some personal contributions on the subjects discussed. Embankments are classified by their use. Different kinds of slope failure and remedial measures are dealt with, as well as investigations for material characterisation and selection.

Study of the foundations, bearing soil and stability of the minaret-tower of Cordobas Mosque

This research in Cordoba's mosque tower main objective was to analyze and characterize the foundations and the underlying soil, calculating the stability of the monument as well as the settlement and deformations performed, using traditional calculation methods and also by finite elements, and to determine differences between both, as well as the stability factor of the Minaret Tower. The works done to study the soil, were drill bores and dynamic penetration tests, classification of samples by size, Atterberg limits, physical and chemical analysis, showing the typical geotechnical composition of the Guadalquivir valley: an alluvial material composed by sand, gravels and silt clays. To study the foundations, inclined bore samples were extracted with 10-65º, showing calcarenite stone ashlars and lime concrete alternating with stone and brick masonry.

A linearization-based system reliability method with application to geotechnical analysis

Los análisis de fiabilidad representan una herramienta adecuada para contemplar las incertidumbres inherentes que existen en los parámetros geotécnicos. En esta Tesis Doctoral se desarrolla una metodología basada en una linealización sencilla, que emplea aproximaciones de primer o segundo orden, para evaluar eficientemente la fiabilidad del sistema en los problemas geotécnicos. En primer lugar, se emplean diferentes métodos para analizar la fiabilidad de dos aspectos propios del diseño de los túneles: la estabilidad del frente y el comportamiento del sostenimiento. Se aplican varias metodologías de fiabilidad — el Método de Fiabilidad de Primer Orden (FORM), el Método de Fiabilidad de Segundo Orden (SORM) y el Muestreo por Importancia (IS). Los resultados muestran que los tipos de distribución y las estructuras de correlación consideradas para todas las variables aleatorias tienen una influencia significativa en los resultados de fiabilidad, lo cual remarca la importancia de una adecuada caracterización de las incertidumbres geotécnicas en las aplicaciones prácticas. Los resultados también muestran que tanto el FORM como el SORM pueden emplearse para estimar la fiabilidad del sostenimiento de un túnel y que el SORM puede mejorar el FORM con un esfuerzo computacional adicional aceptable. Posteriormente, se desarrolla una metodología de linealización para evaluar la fiabilidad del sistema en los problemas geotécnicos. Esta metodología solamente necesita la información proporcionada por el FORM: el vector de índices de fiabilidad de las funciones de estado límite (LSFs) que componen el sistema y su matriz de correlación. Se analizan dos problemas geotécnicos comunes —la estabilidad de un talud en un suelo estratificado y un túnel circular excavado en roca— para demostrar la sencillez, precisión y eficiencia del procedimiento propuesto. Asimismo, se reflejan las ventajas de la metodología de linealización con respecto a las herramientas computacionales alternativas. Igualmente se muestra que, en el caso de que resulte necesario, se puede emplear el SORM —que aproxima la verdadera LSF mejor que el FORM— para calcular estimaciones más precisas de la fiabilidad del sistema. Finalmente, se presenta una nueva metodología que emplea Algoritmos Genéticos para identificar, de manera precisa, las superficies de deslizamiento representativas (RSSs) de taludes en suelos estratificados, las cuales se emplean posteriormente para estimar la fiabilidad del sistema, empleando la metodología de linealización propuesta. Se adoptan tres taludes en suelos estratificados característicos para demostrar la eficiencia, precisión y robustez del procedimiento propuesto y se discuten las ventajas del mismo con respecto a otros métodos alternativos. Los resultados muestran que la metodología propuesta da estimaciones de fiabilidad que mejoran los resultados previamente publicados, enfatizando la importancia de hallar buenas RSSs —y, especialmente, adecuadas (desde un punto de vista probabilístico) superficies de deslizamiento críticas que podrían ser no-circulares— para obtener estimaciones acertadas de la fiabilidad de taludes en suelos. Reliability analyses provide an adequate tool to consider the inherent uncertainties that exist in geotechnical parameters. This dissertation develops a simple linearization-based approach, that uses first or second order approximations, to efficiently evaluate the system reliability of geotechnical problems. First, reliability methods are employed to analyze the reliability of two tunnel design aspects: face stability and performance of support systems. Several reliability approaches —the first order reliability method (FORM), the second order reliability method (SORM), the response surface method (RSM) and importance sampling (IS)— are employed, with results showing that the assumed distribution types and correlation structures for all random variables have a significant effect on the reliability results. This emphasizes the importance of an adequate characterization of geotechnical uncertainties for practical applications. Results also show that both FORM and SORM can be used to estimate the reliability of tunnel-support systems; and that SORM can outperform FORM with an acceptable additional computational effort. A linearization approach is then developed to evaluate the system reliability of series geotechnical problems. The approach only needs information provided by FORM: the vector of reliability indices of the limit state functions (LSFs) composing the system, and their correlation matrix. Two common geotechnical problems —the stability of a slope in layered soil and a circular tunnel in rock— are employed to demonstrate the simplicity, accuracy and efficiency of the suggested procedure. Advantages of the linearization approach with respect to alternative computational tools are discussed. It is also found that, if necessary, SORM —that approximates the true LSF better than FORM— can be employed to compute better estimations of the system’s reliability. Finally, a new approach using Genetic Algorithms (GAs) is presented to identify the fully specified representative slip surfaces (RSSs) of layered soil slopes, and such RSSs are then employed to estimate the system reliability of slopes, using our proposed linearization approach. Three typical benchmark-slopes with layered soils are adopted to demonstrate the efficiency, accuracy and robustness of the suggested procedure, and advantages of the proposed method with respect to alternative methods are discussed. Results show that the proposed approach provides reliability estimates that improve previously published results, emphasizing the importance of finding good RSSs —and, especially, good (probabilistic) critical slip surfaces that might be non-circular— to obtain good estimations of the reliability of soil slope systems.

Cone penetration testing to assess slope stability in the 1979 Nice landslide area (Ligurian Margin, SE France)

In the landslide-prone area near the Nice international airport, southeastern France, an interdisciplinary approach is applied to develop realistic lithological/geometrical profiles and geotechnical/strength sub-seafloor models. Such models are indispensable for slope stability assessments using limit equilibrium or finite element methods. Regression analyses, based on the undrained shear strength (su) of intact gassy sediments are used to generate a sub-seafloor strength model based on 37 short dynamic and eight long static piezocone penetration tests, and laboratory experiments on one Calypso piston and 10 gravity cores. Significant strength variations were detected when comparing measurements from the shelf and the shelf break, with a significant drop in su to 5.5 kPa being interpreted as a weak zone at a depth between 6.5 and 8.5 m below seafloor (mbsf). Here, a 10% reduction of the in situ total unit weight compared to the surrounding sediments is found to coincide with coarse-grained layers that turn into a weak zone and detachment plane for former and present-day gravitational, retrogressive slide events, as seen in 2D chirp profiles. The combination of high-resolution chirp profiles and comprehensive geotechnical information allows us to compute enhanced 2D finite element slope stability analysis with undrained sediment response compared to previous 2D numerical and 3D limit equilibrium assessments. Those models suggest that significant portions (detachment planes at 20 m or even 55 mbsf) of the Quaternary delta and slope apron deposits may be mobilized. Given that factors of safety are equal or less than 1 when further considering the effect of free gas, a high risk for a landslide event of considerable size off Nice international airport is identified