Tailings ponds are one of the three major production facilities in metal mines. The volume of tailings increases year by year, but the storage capacity of existing tailings ponds is limited. Therefore, tailings dams must become more fine-grained and larger. The potential hazard they represent should not be underestimated. This paper reveals the causes and regional distribution patterns of 342 tailings dam failures globally from 1915 to 2021 through statistical analysis. It was found that tailings pond failures occur almost every year, with an average of 4.4 accidents/year (1947–2021). The frequency has been gradually increasing in recent years, and most tailings pond failures are directly related to heavy rainfall or earthquakes. The frequency of tailings pond failures was significantly higher in Asia (21.3%) and the Americas (57.9%), especially in China (n = 43) and the United States (n = 107). Causes of tailings pond failures differed among regions. Most tailings pond failures in Asia and Europe were related to hydroclimate, while those in South America were mainly triggered by earthquakes. This study will provide theoretical data for the pre-design as well as the safe and stable operation of global tailings ponds, which will help to prevent global tailings pond failures.

Unusual rainfall is the primary cause of the failure of the tailing dams, and overtopping is the most representative model of the tailing dam failure. The upstream tailing dam was selected as the research object to study the whole process of breach extension and the overtopping dam-failure mechanism under the full-scale rainfall condition. The results showed that the significant size grading phenomenon in the front, middle, and end of the tailing pond was obvious due to the flow separation effect, and its average particle diameter was D50. At different moments of rainfall, the height of the infiltration line at different positions of the dam body was different; at the rainfall of 3600 s, the height of the infiltration line lagged behind the height of the tailing pond, and this phenomenon from the tail of pond to the outside of the dam slope became more obvious. After the rainfall of 3600 s, the height of the infiltration line lagging behind the water level in the pond basically disappeared, and the rate of infiltration line rise kept pace with the rate of water level. The process of overtopping dam-failure experienced dam overtopping (gully erosion), formation of a multistepped small “scarp,” breach rapid expansion, formation of large “scarp,” and burst (fan-shaped formation). The width and depth of the breach showed a positive correlation, and the widening rate of the breach was 3 to 8 times of the deepening rate, especially in the middle of the dam break, widening behavior occupied the dominant factor. The shape of the dam body after failure was parabolic, and the dam body had obvious elevation changes. These results provide the theoretical guidance and engineering application value for improving the theory and early warning model of the upstream tailing dam.

The shape of tailing particles is essential factors of their macroscopic mechanical properties. Scholars have studied the effects of controllable factors, such as loading method, confining pressure, and strain rate, on the strength of tailing sand. However, research on the tailing particle structure and shape through laboratory tests has proved to be difficult due to the uncertain and discrete tailing particle distribution. Thus, the macro-mesoscopic response of heteromorphic tailing particles is rarely investigated. In this paper, the macro-mesoscopic response of heteromorphic tailing particles is studied using multisphere approximation, and numerical simulation of triaxial tests on the particles is conducted. Nonlinear evolution patterns of porosity, internal friction angle, and cohesion of heteromorphic tailing particles with the variation of angularity were investigated using the flexible boundary program developed in this study, which revealed the intrinsic relationship between the mesostructure evolution mechanism and the macroscopic engineering characteristics of heteromorphic tailing particles. The research results showed that (1) changes in angularity led to tailing particle rearrangements and, in turn, porosity changes. With increased angularity and confining pressure, particle sphericity decreased, and the deviatoric and peak stress increased accordingly. In the meantime, the softening was more significant as the peak stress was exceeded, while the cohesive force generally increased. (2) With fixed particle shape and angularity, the internal friction angle decreased slightly as the effective confining pressure increased. (3) In the shearing process, the simulated contact force chain evolution of tailing particles with different shapes was similar. The disordered contact force chains gradually undergo directional connection, i.e., the increased confining pressure reduced the number of free tailing particles and increased the number of stressed particles. (4) The triaxial stress-strain and peak stress in rigid boundary simulations under different confining pressures were slightly lower than those in the flexible boundary simulations. However, the difference was insignificant, indicating the good feasibility and reasonability of rigid boundary simulations for the macroscopic mechanical behaviors in triaxial tests. The research results could offer more direct insights into the macro-mesoscopic response and mechanical mechanisms of nonspherical particles and provide references for the simulation of tailings at the microscopic levels.

Tailings ponds are the indispensable facilities in the mine production and operation. Once the dam is destabilized and damaged, it will pose a serious threat on the life and property of the downstream population and could also potentially cause an environmental disaster. With an engineering background, this paper dynamically and numerically simulates the evolution process of tailings flow from dam failure and the influence scope of any resulting disaster in context. The evolution characteristics of leaked tailings flow are analyzed at various downstream riverbed slopes and debris blocking dam settings. In addition, parameters such as flow rate, impact force and deposition range of leaked tailings flow at downstream arrival are studied, as well as their correlations. The results indicate that the flat terrains upstream and downstream of passage zone show a relatively larger area of inundation by tailings flow. Both the maximum and final downstream inundated ranges increase with the elevating slope of downstream riverbed, and the leaked tailings are deposited mainly in the nearby villages in front of the dam and the flat terrains of the downstream passage zone. Additionally, rational establishment of debris blocking dams on the downstream side is effective in diminishing the damage of tailings flow to the downstream section. This study can also provide an important basis for the quantitative evaluation of post-disaster influence scope for tailings pond as well as for the design of dam body.

Tailings pond accidents frequently occur during an extended period, resulting in loss of life and property, wastage of resources, and environmental pollution. Relying on tailings pond engineering, this paper carried out sample particle fragmentation experiments and settling column experiments to explore the deposition distribution pattern of tailings in both horizontal and vertical directions as well as the impact of particle size distribution on the sedimentation stratification effect. The results show that the median particle size on the dry beach surface in the horizontal direction slowly decreased with the increase in the distance from the subdam. The particle size of tailings showed great fluctuations in the vertical direction, which gradually became finer with the increase in the depth overall. At the same time, saturated sedimentation experiments suggested the inconsistent variation rule with the field test, namely, coarse on the bottom and fine at the top, and the change in particle size greatly affected the tailing sedimentation stratification effect. With the increase in fine particle content in tailings, the appearance time of the water-sand interface was shortened to within 30 min, but the sedimentation and consolidation completion times were delayed to about 1400 min. The settling column results indicate that the increase in fine particle content gradually weakened the sedimentation stratification effect, and the sedimentation pattern transformed from independent sedimentation to floc-type average sedimentation, which led to the enhanced water-retaining property of the settled layer. This may lead to an increase in the saturation line and a decrease in the length of the dry beach, seriously affecting the safe operation of the tailings pond. The research results provide some theoretical guidance and basic data for analyzing the consolidation efficiency of tailings and the stability of the tailings pond.

Background: It is adverse for the safety of a tailings dam to use fine-grained tailings as the materials for a high tailings dam because of the low penetration coefficient, the slow consolidating velocity, and the bad physical mechanical property. Furthermore, with the influence of complicated geography conditions, the phreatic line will be increased enormously when encountering special conditions, which directly affect the safe operation of the tailings dam. Methods: In this study, based on the engineering, geological, and hydrogeological conditions and survey results of a tailings dam, a 210 m fine-grained tailings dam located in three gullies was selected and used to simulate the three-dimensional seepage field of a tailings dam under a steady saturated state by using the finite element software MIDAS GTS. The permeability coefficient was inverted, the seepage field of the project under different working conditions was simulated, and the position of the phreatic line was obtained. The controlled position of phreatic lines was determined by combining the seepage field with the stability requirements. Results: Back analysis could accurately reflect the actual permeability coefficient of each partition of tailings dams. Due to the multiple areas of seepage accumulation, large valley corners, and narrowing of the dam axis, the phreatic line of the shoulder region was elevated by 2~3 m compared to the surrounding area and was thereby the most critical region of the tailings dam seepage control. The stability requirements and minimum controlled position of the phreatic line requirements could be met when the controlled position of the phreatic line was 23 m. Conclusion: This study revealed the key areas and reasons why the tailings dam’s phreatic line is prone to be uplifted under complicated geography conditions. It was very critical to control the local phreatic line by adopting local horizontal seepage drainage measures or radiation wells in the key areas of the tailings dam to ensure the safety of the tailings dam. In addition to strengthening the daily monitoring of the key areas and the exfiltration facilities of the tailings dam, it is recommended to carry out determination tests of the permeability coefficient and particle size at regular intervals. The findings could provide countermeasures for seepage control.