Blog

Seepage stability analysis of geogrid reinforced tailings dam | Scientific Reports

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Scientific Reports volume  14, Article number: 1814 (2024 ) Cite this article Non Woven Geotextile Fabric

Seepage stability analysis of geogrid reinforced tailings dam | Scientific Reports

To investigate the influence of a geogrid-reinforced tailings dam on the seepage stability of the dam body, this paper was based on the field test of a reinforced tailings accumulation dam. The study utilized the finite element strength reduction method to simulate the stability of the main dam of the Fengshuigou tailings reservoir under different seepage conditions using ABAQUS software. Additionally, the paper discussed the impact of conventional heightening, dry beach length, and geogrid reinforcement on the position and safety factor of the saturation line of the dam body. The results showed that when the dam body was raised, the saturation line rose by 2.8–5.3 m, resulting in a decrease in the safety factor. The geogrid effectively reduced the height of the saturation line in the tailings dam. In comparison to the unreinforced condition (dam heightening), the saturation line of the tailings dam decreased by 0.9–2.8 m under the local reinforcement condition and by 3.2–12.5 m under the overall reinforcement condition. The geogrid significantly improved the stability of the tailings dam. Furthermore, under the local reinforcement condition, the safety factor of the dam increased by 3.8–5.5%, and under the overall reinforcement condition, it increased by 35.9–42.9%, when compared to the unreinforced condition. Increasing the dry beach length improved the stability of the tailings dam, and under normal working conditions, the safety of the tailings dam was much higher than under the minimum dry beach condition. These results served as a reference for the design of the dam and the new tailings reservoir, laying a foundation for the sustainable development and environmental protection of the tailings pond.

Caroline Chalumeau, Hans Agurto-Detzel, … Audrey Galve

Zhen Zhang, Min Liu, … Jinrong Su

Xuejiao Hou, Danghan Xie, … Jaap H. Nienhuis

The stability of the tailings dam was closely related to the safe operation of the whole tailings reservoir. According to a review of relevant information1, the failure of a tailings dam seriously threatened the safety, life, and property of downstream inhabitants, as well as polluted and damaged the ecological environment. For example, in 2019, the Córrego do Feijão tailings dam I in Brumadinho, Minas Gerais, broke and released a large amount of tailings, resulting in the deaths of 660 people and polluting downstream rivers2. In 1985, a dam failure in a tailings pond in Stava, northern Italy, resulted in 268 deaths and significant economic losses3. With governments’ increasing attention to the safety of tailings ponds, the overall safety level has been significantly improved. However, the problem of heavy rainfall climate has gradually become the main factor causing the tailings dam to break4. Because the permeability of the tailings dam was greatly affected by the seepage field, the location of the saturation line in the seepage field was closely related to the length of the dry beach and the rainfall intensity5.

To make an accurate evaluation of the stability of the tailings dam, similar physical model tests and numerical simulations were usually used to evaluate. For the physical model, selecting reasonably similar materials was particularly important for the test results. However, choosing similar materials often involves many subject knowledge and was more complicated. In contrast, numerical simulation has gradually become an important method to evaluate the stability of tailings dams. At present, many scholars have used the finite element method to establish two-dimensional or three-dimensional models to carry out research6,7. Lu et al.8 proposed that proper simplification and generalization of complex terrain in three-dimensional (3D) numerical calculation had little impact on the results and could meet the accuracy requirements. Based on the stochastic limit equilibrium method, Mafi et al.9 analyzed the dam’s stability with three different slopes. Dastpak et al.10 analyzed the stability of geosynthetic reinforced slopes based on non-circular certainty and randomness. Aroni Hesari et al.11 used the horizontal slice method to study the seismic internal stability of geosynthetic reinforced soil slopes. Fatehi et al.12 used the pseudo-static method to examine the stability of reinforced slopes under seismic load. Doğan and Güllü13 proposed a 3D voxel model generation method for finite element structural analysis. Wang et al.14 analyzed the stability of tailings dams under dry–wet cycles and proposed an effective calculation method for the saturation line of tailings dams under dry–wet cycles. Wang15 analyzed the seepage condition of the dam body under the current elevation and the final design elevation of the tailings pond through theoretical analysis and numerical calculation. Zhang et al.16 analyzed the influence of different dry beach and upstream-side slope ratios on the seepage stability of the tailings dam. Naeini et al.17 used SIGMA/W software to analyze the stress-pore pressure coupling. It can be seen from the above research results that the limit equilibrium method was mainly used to solve the safety factor of the tailings dam. When the limit equilibrium method was used to analyze the influence of pore water pressure on the stability of the tailings dam, the pore water pressure was treated as zero in the case of an unsaturated area, ignoring its influence. At the same time, many factors affected the stability of the tailings dam, but most were related to saturation line, dry beach length, and pore pressure. The research on seepage stability of tailings dams reinforced by geogrid was relatively weak.

Given how to improve the tailings dam stability, many scholars conducted relevant research, mainly including the following aspects: First, they used geosynthetics or fibers to reinforce tailings dams. It included geogrid, geotextile, geomembrane, geobag, various fibers, etc. Xiao18 used a Geo-slope to analyze the geosynthetic-reinforced Jinshandian tailings dam. The results showed that reinforcement could effectively improve the stability of the fill dam. Zhou et al.19 studied tailing-geotextile composites, analyzed their permeability and physical properties, and generalized the composite model. Li et al.20 established a stability analysis method based on limit equilibrium theory to simulate the stability of tailings dams using geotextile bags. Zheng et al.21 used basalt fiber to reinforce a tailings dam and improved its stability. Second, they used chemical reagents or biological consolidation effectively to improve the stability of the dam. Yang et al.22 used polyacrylamide to enhance the dynamic characteristics of tailings, thereby improving the stability of tailings dams. Lu et al.23 used microbially induced calcite precipitation to reinforce tailings. Sun et al.24 studied the mechanical reinforcement effect of a tailings dam slope under four conditions: goosegrass, vetiver grass, clover, and bare slope. The results showed that vetiver grass was more conducive to dam reinforcement. The third aspect was to improve the dam’s stability from the perspective of anti-seepage, mainly involving material properties and foundation treatment. In addition to using anti-seepage membranes, the ground could be improved by grouting to have low permeability, low settlement, and high bearing capacity25,26. At the same time, when the tailings dam produced a seepage channel, it could be blocked by grouting27,28. The fourth approach was to use online monitoring technology to monitor the dam in real-time. Olivier et al.29 used environmental-seismic noise to monitor the stability of tailings dams, which provided a valuable tool for the remote monitoring of the structural stability of tailings dams. In summary, many scholars have conducted a lot of research on the seepage stability and dam reinforcement of tailings dams. Although some achievements have been made, the study mainly focused on the influence of a single factor. The research on the impact of geogrid-reinforced tailings dams under multi-factor coupling still needs further discussion.

In this study, based on the previously conducted field test30 of the Fengshuigou tailings reservoir, reasonable model parameters were selected by ABAQUS finite element software, and the conventional accumulation dam and geogrid reinforced accumulation dam of the field test were simulated and calculated, respectively. By comparing the test and simulation results, the accuracy of the model parameter inversion was assessed. Based on the fluid–solid coupling theory and the finite element strength reduction method, two tailings dam models with two tailings dam elevations (current and design elevations) and two reinforcement types (geogrid local and overall reinforcement) were established. The stability of the tailings dam under three dry beach lengths (minimum dry beach, flood operation, and normal operation) was calculated. The influence of dry beach length, conventional heightening, and geogrid reinforcement on the saturation line and safety factor of the dam body was studied (see Fig. 1 for the specific process), and reference suggestions were made for the engineering application of reinforced tailings reservoirs.

Flowchart of the research process utilized in this study.

This study was based on field tests of reinforced tailings stacking dams in the Fengshuigou tailings pond30. ABAQUS simulations were performed for the field tests, a finite element model was established, an inverse trial analysis was performed, and the simulation results were compared with the field monitoring results to obtain a more reasonable calculation model and calculation parameters. On this basis, we analyzed the influence law of geogrid reinforcement on the pore pressure (saturation line), internal pressure, and deformation of stacked dams.

The field test of the reinforced accumulation dam was conducted on the dry beach surface of the No.5 auxiliary dam of the Fengshuigou tailings reservoir, and the side length of the specific test area was approximately 15 m. After excavating about 1 m shallow pits in the test area, a 4 m high accumulation dam was built, including conventional dams, and two-way geogrid reinforced dam. The outer slope ratio of the surrounding accumulation dam was 1:1, and the width of the dam crest was 3 m. The field test of the reinforced tailings dam is shown in Fig. 2a, and a 3D schematic diagram is shown in Fig. 2b. The performance parameters of the geogrid used are listed in Table 1.

Field test of reinforced tailings dam: (a) Test site diagram; (b) Three-dimensional diagram.

In this study, both the stacked and geogrid dams were modeled as solids in the finite element calculations. The Mohr–Coulomb plastic model was utilized for the stacked dam and subgrade, and the elastic model was applied for the geogrid. The embedded constraint simulated the interaction between the tailings sand and the geogrid. Considering the excessive contact surface of the model, the calculation quickly became non-convergent and held calculation time; thus, it was assumed that no relative sliding occurred between the tailing sand and the geogrid. The strength reduction method was used to calculate the stability of the dam. A different strength reduction coefficient \(F_{r}\) was assumed in the calculation, and the finite element analysis was carried out according to the strength parameters after reduction to observe whether the calculation is convergent. In the calculation process, \(F_{r}\) (using dichotomy) was continuously increased, and the strength reduction factor \(F_{r}\) at the critical failure was the slope stability safety factor \(F_{s}\) . The shear strength parameters after reduction are shown in Eqs. (1) and (2):

where \(c\) and \(\varphi\) were the shear strength of the soil; \(c_{m}\) and \(\varphi_{m}\) were the shear strength required for the actual; \(F_{r}\) was the strength reduction factor.

Only gravity loads were considered in the model calculation, without considering external factors such as overload and temperature. The cross-sectional view of the model of the unreinforced and gridded reinforced pile dam is shown in Fig. 3. The laying of the geogrid is divided into four layers: the first layer is laid on the bottom of the accumulation dam; the second, third, and fourth layers are laid every 1 m in the middle; and the top surface is not laid with reinforcement. Based on the engineering investigation report and with reference to the literature31, the model calculation parameters are listed in Table 2.

Cross section of tailings accumulation dam (unit: m): (a) Conventional dam; (b) Geogrid reinforced dam.

Simulations were performed to reduce the influence of boundary effects by expanding the boundaries in the left and lower directions of the tailings dam model. Based on data from the literature32, the distance from the left boundary to the toe of the slope of the computational model was taken as 1.5 times the slope height, and the bottom of the slope extended downward by 1 times the slope height. The boundary conditions were set as follows: the bottom boundary of the stacked dam was constrained in the x, y, and z directions, the upper and lower boundaries were constrained in the x direction only, the boundaries of both sides of the section were constrained in the y direction only, and the top of the slope and the slope surface were not constrained. The pore pressure boundary on the right side of the model was determined by the relationship between the water level in the reservoir of the surrounding pile dam and time, and the slope and foot of the slope were set as the free seepage section boundary.

CPE8P cells were used to mesh the two types of pile dams, and 1298 soil cells were used after dividing the cells of the conventional pile dam. Then, the four layers of geogrid reinforcement were divided into 438, 589, 759, and 759 cells. The stacked dam cells and geogrid cells are shown in Fig. 4. In the finite element simulation, we also simulated the discharge slurry time of 3 h, so the numerical simulation results showed a sharp decline that recovered within 3–5 h.

Unit mesh model of tailings accumulation dam.

A permeability coefficient test was conducted out to obtain a permeability coefficient close to the actual operating conditions of the project. The tailing material parameter affecting the measured water level was the permeability coefficient of each soil layer; therefore, the permeability coefficient was selected for stochastic inversion. The parameters involved in this permeability coefficient inversion included the permeability coefficients of the stacked dam body and geogrid. According to the report provided by the geogrid manufacturer and the related literature30, the range of the permeability coefficient and the initial value of the trial calculation were determined, as shown in Table 3.

In this section, the positive algorithm was used as the basic method for inversion parameter trial calculation, combined with ABAQUS software for field test simulation. According to the two types of monitoring data, pore and soil pressure, the actual monitoring values of pore and soil pressure were analyzed and compared with the simulated calculated values, and the distribution law was derived. The model parameters were adjusted and calculated individually until the difference between them met the accuracy requirements, thus determining the infiltration parameters that could be used.

The main field test monitoring instruments were pore water pressure monitoring instruments (seepage pressure meters) and internal soil pressure monitoring instruments (soil pressure boxes). Two seepage pressure gauges were placed horizontally at 2 m intervals in the middle cross-section of each stacking dam, and one earth pressure box was placed horizontally and vertically in the middle cross-section of the stacking dam, as shown in Fig. 5.

Deployment of field measurement instruments.

The simulated and field-monitored pore-water pressure variations over time for the two types of stockpile dams after slurry injection into the surrounding tailings dams are shown in Figs. 6 and 7, respectively. From the figure, it can be seen that the simulated and monitored values of the pore pressure inside the two types of stacked dams agreed with each other as a function of time. The variation laws were as follows: For the simulated pore pressure, the trend over time was basically the same as the trend of the water level in the reservoir of the surrounding stacked dam. That is, the simulated pore pressure increased linearly with time for 3 h before the start of the observation, gradually decreased linearly to a stable value after reaching a peak, and then tended to stabilize. For the monitored pore pressure, the pore pressure in the first 6 h of observation gradually increased with time and then tended to stabilize.

Simulation and monitoring of pore water pressure with time for conventional accumulation dams: (a) Simulation value; (b) Monitoring value.

Simulation and monitoring of pore water pressure with time for geogrid reinforced accumulation dams: (a) Simulation value; (b) Monitoring value.

Further analysis of the variation in the dam pore pressure throughout the monitoring process was presented in Table 4. The simulated peak pore water pressure in the table was compared with the measured peak, and it can be seen that the error was 0.3–6.5%. That is, the simulated peak pore water pressure of the stacked dam closely matched the measured peak. This indicates that the calculation model and parameters determined by the inverse trial calculation are reasonable and could better match the actual situation. The final stable value of pore pressure of the conventional dam was the largest, and the stable value of pore pressure of the geogrid-reinforced dam was not much different, indicating that reinforcement had a significant effect on the “resistance” of pore pressure inside the tailings dam. Reinforcement of tailings dams could effectively reduce the water level inside the dam body. This was because the reinforcement treatment of the tailings dam could make the reinforced tailings complex form a drainage prism, promoting the discharge of water inside the dam body, and reducing the pore water pressure inside the dam body.

The simulated and field-monitored values of earth pressure (lateral and vertical earth pressure) over time for the two types of stockpile dams after slurry injection into the surrounding tailings dams are shown in Figs. 8 and 9, respectively. From the figure, it can be observed that the simulated and monitored values of the soil pressure inside the two types of stockpile dams are consistent with each other as a function of time, and the variation patterns are as follows. For the simulated soil pressure, its trend with time is basically consistent with the trend of the water level in the reservoir of the surrounding pile dam. That is, the lateral soil pressure and vertical soil pressure both increased linearly with time continuously in the first 3 h of observation, reached a peak, gradually decreased linearly to a stable value, and then tended to stabilize. For monitored soil pressure, both transverse and vertical soil pressure gradually increased with time and then stabilized, but transverse pressure changed more drastically compared to vertical soil pressure before stabilization.

Simulation and monitoring of soil pressure with time for conventional accumulation dams: (a) Lateral soil pressure; (b) Vertical soil pressure.

Simulation and monitoring of soil pressure with time for geogrid reinforced accumulation dams: (a) Lateral soil pressure; (b) Vertical soil pressure.

Further analysis of the changes in the dam earth pressure throughout the monitoring process was presented in Table 5. According to the table, the simulated peak soil pressure inside the two types of dams basically matched the monitored peak with an error between 1.7 and 6.9%. This indicated that the calculation model and parameters determined by the inverse trial calculation were reasonable and can better matched the actual situation.

The possible reasons for the errors in pore and earth pressure simulation and detection values are as follows: (1) There was a deviation between the position of the output point of the simulated time field and the position of the measuring point in the test; (2) The osmometer or earth pressure box in the test was displaced; (3) The field test was affected by the environment, weather, and other factors, leading to a deviation between the simulated and actual working conditions; (4) The mesh generation was not accurate in the process of establishing the model.

After the above analysis, the numerical value calculated by the calculation model determined by the trial calculation was consistent with the actual situation. Therefore, the determination of model parameters was reasonable, and the permeability coefficients of the dam body and geogrid were determined to be 2.0 × 10−5 and 1.8 × 10−4 cm s−1, respectively.

The main influencing factors in the construction process of a reinforced tailings dam were the water content and relative density. Therefore, the relative density and water content were used as the influencing factors to optimize the model further, and the YT1200 direct shear drawing test machine was used for the test. According to the compaction test and drying test, the relative density of tailings under natural conditions was 0.6, and the natural moisture content was 6.4%. Four working conditions with relative densities of 0.65, 0.75, 0.85, and 0.95 were selected, and the water content was natural and saturated water content (20.8%). The normal stress was determined to be 10 kPa, 20 kPa, 30 kPa, and 40 kPa through the preliminary examination to avoid the geogrid being pulled off during the test. The test results are shown in Fig. 10.

Interface characteristics of geogrid: (a) Internal friction angle; (b) Cohesion.

It was observed from the figure that when the relative density increases from 0.65 to 0.95, the curves of internal friction angle and cohesion show a monotonically increasing trend regardless of natural water content or saturated water content. With the increase of relative density, the shear strength under saturated water content was lower than that under natural water content. In summary, in the modeling of the geogrid reinforced tailings dam, the water content of tailings should have been determined by test to determine the shear strength parameters of the tailings dam. In the absence of test conditions, the interface shear strength parameters of saturated tailings below the saturation line should have been selected at a lower value, which was related to the stability and safety of the reinforced structure.

Taking the main dam of the Fengshuigou tailings reservoir as the engineering background, the current elevation of the tailings reservoir reached close to 160 m, which approached the design elevation. Currently, the mining company intended to extend its service life and planned to expand its capacity to a design elevation of 175 m. At this time, the height of the accumulation dam was 100 m, the maximum dam height reached 120 m, the calculated total storage capacity was approximately 68,443 million m3, and the tailings pond was classified as the first grade. The initial dam of the main dam of the Fengshuigou tailings reservoir was a permeable rockfill dam. The dam crest width was 4 m, the upstream slope ratio was 1:1.85, the downstream slope ratio was 1:2, the dam bottom elevation was 55 m, the dam crest elevation was 75 m, and the dam height was 20 m. There was a 4 m wide horse road at an elevation of 66 m on the outer slope of the initial dam, and the dam foundation was located on the gravel layer. The accumulation dam was constructed upstream. The height of each sub-dam was 5 m. The ratio of the inner and outer slopes of the sub dam was 1:2. The top width of the sub-dam was 15 m. The average outer slope ratio of the accumulation dam was 1:5. There was a 50 m wide platform at an elevation of 120 m.

According to the tailings pond site survey information, the tailings dam was generalized, and the heightened dam body will be extended according to the generalized layering line of the current dam body. Therefore, considering the most dangerous situation during the simulation, drainage wells and other auxiliary facilities were ignored. The simplified and generalized calculation diagram is shown in Fig. 11, and the simulation diagram is shown in Fig. 12.

Calculation diagram of the current situation of the main dam of the tailings pond (Unit: m).

The current finite element model of the main dam of the tailings pond.

In the design of the reinforcement scheme, it was assumed that33:

The tailings dam was in normal operation, and there were no unexpected natural disasters such as heavy rainfall, mountain torrents, and earthquakes in the area.

The geogrid and filler were closely meshed and embedded, there was no relative slip with the soil, and deformation coordination was maintained.

The mechanical properties of the geogrid were intact; the geogrid was free from damage, fracture, and loose splicing during construction and laying.

The tailings stockpile dam model was appropriately simplified in the finite element model, and the tailings stockpile dam was considered to have been formed by the accumulation of fine-grained tailings. The composition of the main dam profile and the physical and mechanical parameters of the geogrid are listed in Table 6.

The three sides of the tailings accumulation dam were surrounded by mountains and belonged to a valley-type tailings reservoir. A finite element numerical model was established by considering the distribution of the tailings accumulation dam on the longitudinal section of the central axis. The boundary conditions were as follows: the displacement of the dam bottom was fixed in three directions, the normal displacement of other surfaces was constrained, the upper end was the free boundary, the 55 m elevation head was set at the bottom of the primary dam, and the upstream elevation head was set according to the length of the main dam dry beach.

The Mohr–Coulomb elastoplastic model was applied to the constitutive model of rock and soil; the elastic model was selected for the reinforcement, in which the section properties of the grid and the section properties of the rock and soil are entities. The soil steady-state analysis step was adopted. Because the constitutive model of the rock and soil was the Mohr–Coulomb model, the asymmetric unsymmetric algorithm was used to analyze and calculate.

For calculations, we referred to the “code for design of tailings facilities”34 (GB 50863-2013) and “code for seismic design of structures”35 (GB 50191-2012), carried out simulation calculations according to the following conditions, and compared and analyzed the changes of saturation line position and safety factor of the main dam under different working conditions:

First, according to the current situation and future design plan of the tailings dam, and considering that the heightening affected the seepage field and stability of the tailings dam, the current elevation of 160 m and the final design elevation of 175 m were considered (See Figs. 11 and 13a);

Second, to analyze the influence of reinforcement on the state of the tailings dam, considering the reinforcement of 160–175 m elevation section and 75–175 m elevation section, geogrid was selected as the reinforcement; it was only laid in the accumulation dam, the grating laying length was 200 m, and the spacing was 2.5 m (Fig. 13b,c).

Third, considering the change in the water level in the tailings pond, the length of the dry beach was calculated by selecting the minimum length of the dry beach, flood operation, and normal operation (refer to Table 7 for the length of the dry beach at different elevations).

Three types of main dams when the final design elevation is 175 m: (a) Unreinforced; (b) Reinforcement of the 160–175 m elevation section; (c) Reinforcement of the 75–175 m elevation section.

Based on the above analysis, this simulation was carried out under the following four working conditions:

Working condition I current main dam elevation (160 m);

Working condition II final design elevation of the main dam without reinforcement (175 m);

Working condition III local reinforcement of the main dam at the final design elevation (160–175 m);

Working condition IV overall reinforcement of the main dam at the final design elevation (75–175 m).

As shown in Fig. 14, under the four minimum dry beach working conditions of 160 m, 175 m without reinforcement, 160–175 m with geogrid, and 75–175 m with geogrid, the distribution of the saturation line when the pore pressure of the tailings dam is 0. The seepage free surface of the tailings dam is curved, and negative pore water pressure appears in the upper part of the saturation line. This area is unsaturated, which is beneficial for the stability of the dam body. The pore pressure at the bottom of the saturation line was positive and saturated. The pore water pressure of the dam body was evenly distributed, and the pore pressure increased gradually from top to bottom. The outlet point of the saturation line of the dam body was near the dam body at the foot of the upstream slope of the primary dam. There is no intersection with the dam surface of the tailings accumulation dam, and groundwater does not overflow from the dam surface.

Cloud chart of pore pressure distribution at minimum dry beach length under different working conditions: (a) Working condition I; (b) Working condition II; (c) Working condition III; (d) Working condition IV.

As shown in Fig. 15, six or seven profiles were selected from the slope of the accumulation dam for analysis. The buried depth of the saturation line from each profile is listed according to the pore pressure distribution, as shown in Fig. 16. Under the four working conditions, when the minimum dry beach was reached, the minimum buried depth of saturation line in each profile was 9.5 m, which was more than 6 m; therefore, it met the requirements of the “code for seismic design of structures” (GB50191-2012), which states that "the buried depth of saturation line on the downstream slope of grade I tailings dam should not have been less than 6 m.” As the length of the dry beach grows, the gradient of the saturation line slowed down, the path became longer and the position of saturation line in the upstream part of the dam body was significantly reduced; however, near the downstream primary dam, the location of the seepage point of the saturation line also decreased, although the change was small.

Saturation line trend and profile layout diagram with different dry beach lengths: (a) Working condition I; (b) Working conditions II, III and IV.

Burial depth of saturation line of each section of dam body under different working conditions: (a) Working condition I; (b) Working condition II; (c) Working condition III; (d) Working condition IV.

Under the four minimum dry beach conditions (current 160 m, no reinforcement 175 m, local reinforcement of geogrid 160–175 m, and overall reinforcement of geogrid 75–175 m), the incremental displacement cloud diagram of the main dam was shown in Fig. 17. It can be observed from the incremental displacement cloud diagram that the sliding surface of the main dam was circular when it was unstable, and the sliding surface appeared at the accumulation dam under the tailings dam platform. This occurred because it was easy for the site to concentrate groundwater to the location caused by displacement deformation and ultimately caused damage. Further analysis showed that with the increase in dry beach length, the depth of the sliding surface gradually became shallow, and the area of the sliding body gradually decreased.

Incremental displacement nephogram for minimum dry beach length under different working conditions: (a) Working condition I; (b) Working condition II; (c) Working condition III; (d) Working condition IV.

In combination with the positioning of the sliding surface of the displacement nephogram, the top of the platform downhill of the accumulation dam was chosen as the characteristic control point, and the sudden change in the horizontal displacement of the characteristic control point was used as the criterion for the instability of the main dam. The relationship between the displacement of the characteristic points and the strength reduction coefficient under different dry beach lengths is shown in Fig. 18. When the current main dam was 160 m, the safety factors of the main dam under different dry beach lengths (minimum dry beach, flood operation, and normal operation) were 1.97, 2.15, and 2.43, respectively; when the elevation of the unreinforced main dam was 175 m, the safety factors of the dam body under different dry beach lengths (minimum dry beach, flood operation, and normal operation) were 1.84, 1.99, and 2.05, respectively. When the geogrid was only reinforced at an elevation of 160–175 m from the main dam, the safety factors of the dam slope under different dry beach lengths (minimum dry beach, flood operation, and normal operation) were 1.91, 2.10, and 2.15, respectively. When geogrid reinforcement was carried out at an elevation of 75–175 m from the main dam, the safety factors of the dam slope under different dry beach lengths (minimum dry beach, flood operation, and normal operation) were 2.50, 2.71 and 2.93, respectively.

Variation of dam safety factor and displacement under different working conditions: (a) Working condition I; (b) Working condition II; (c) Working condition III; (d) Working condition IV.

The simulation results under the three reinforcement conditions of 160 m, 175 m without reinforcement, 160–175 m with geogrid reinforcement, and 75–175 m with geogrid reinforcement are shown in Fig. 19. According to the “code for design of tailings facilities”34, the minimum safety factors for anti-sliding stability of a tailings dam during the minimum dry beach, flood operation, and regular operation were 1.30, 1.20, and 1.10, respectively. According to the safety factor results listed in Fig. 19, under two elevation conditions, the safety factor of the main dam ranged from 1.84 to 2.93 for different dry beach lengths. Consequently, the main dam was able to operate normally under different working conditions, and the safety reserve values were large.

Comparison of calculation results of four working conditions under different dry beach lengths: (a) Minimum buried depth of saturation line; (b) Safety factor.

When the main dam increased from 160 to 175 m, the minimum buried depth of saturation line decreased, that is, the saturation line of the main dam increased by approximately 2.8–5.3 m, which is unfavorable to seepage stability. Simultaneously, the safety factor was reduced by approximately 7–15.6%. Therefore, if the main dam was raised to 175 m, the length of the dry beach should be strictly controlled, that is, the water level in the reservoir should be strictly controlled, and the monitoring of the dam saturation line had to be strengthened to improve the comprehensive drainage capacity of the dam and reduce the water level in the reservoir and the dam saturation line position as much as possible.

Considering the actual situation, after geogrid reinforcement of the main dam from 160 to 175 m, the height of the saturation line decreased by 0–3 m, and the safety factor increased by 0–5.5% under three dry beach lengths. In the ideal state, geogrid reinforcement was carried out on the entire accumulation dam. Under the three dry beach lengths, the maximum decrease of the saturation line height reached 12.5 m, and the maximum increase of the safety factor reached 0.4 times. It could be seen that only the reinforcement at the elevation of 160–175 m of the accumulation dam could reduce the saturation line and improve the safety factor of the dam body, but the effect was not obvious. If the saturation line of the main dam was significantly reduced after the overall reinforcement of the accumulation dam, the harm of seepage could be greatly reduced, and the stability could be significantly improved. Therefore, this study could be referred to in practical engineering. After the construction of the initial dam of the tailings reservoir was completed, geogrid reinforcement was carried out on the accumulation dam.

In this paper, the stability of the tailings dam was studied using the finite element strength reduction method. To further verify the reliability of the simulation results, the limit equilibrium method was used, and the safety factor of the tailings dam under different working conditions was analyzed again by ABAQUS. The specific results are shown in Table 8. It can be seen from the table that the difference in the safety factor of the dam calculated by the two different methods is small. This is consistent with the research results of Cheng et al.36; that is, the calculation results of the limit equilibrium method were slightly smaller than the strength reduction method, and the errors of the two methods were within 5%. This shows that the strength reduction method used in this paper, combined with the results of ABAQUS simulation calculation, is reliable.

Compared with other research results on the safety factor of geogrid-reinforced slopes, Zhang et al.37 used FLAC3D software to compare the stability of geogrid-reinforced and unreinforced dam bodies. The results showed that the safety factor increased from 1.69 to 2.37, 40.70%. In the case of using the strength reduction method in this paper, the safety factor under the minimum dry beach condition increased from 1.84 to 2.50, from 1.99 to 2.71 under flood operation conditions, and from 2.05 to 2.93 under normal operation conditions, achieving 35.87%, 36.18% and 42.93% growth, respectively. In the case of the limit equilibrium method, the safety factor increased from 1.75 to 2.39 under the minimum dry beach condition, from 1.90 to 2.60 under the flood operation condition, and from 1.97 to 2.83 under the normal operation condition. The corresponding increases are 36.57%, 36.84%, and 43.65%, respectively. The results demonstrate that the calculation results of the two methods are consistent with the research results of Zhang et al., especially under normal operating conditions. This also indicates that the geogrid significantly improved the seepage stability of the tailings dam.

In this paper, Fengshuigou tailings reservoir was taken as an example. The finite element strength reduction method was employed to simulate and calculate using ABAQUS. The location and safety factor of the saturation line under different working conditions were compared and analyzed. The influence of dam height, geogrid reinforcement, and dry beach length on the stability of the tailings dam under seepage was discussed, providing reference for similar tailings dam projects.

The problem of seepage stability led to the destruction of the tailings dam structure and eventually resulted in a dam break. Referring to the relevant data38,39,40, the probability of dam failure accidents due to flow sliding under various dam construction processes was shown in Table 9. It was observed from the table that, under the upstream dam construction process, the failure rate of tailings dam flow sliding was the highest, reaching 88%, significantly greater than other dam construction processes. As is known, a tailings dam break was a serious environmental disaster with potential long-term negative impacts on the surrounding environment and human society1. It caused significant losses to the local economy and led to casualties and disappearances of nearby residents. Therefore, the seepage stability of the tailings dam could not be ignored. The reinforcement of geogrid played an active role in improving the seepage stability of the tailings dam. It enhanced the shear strength of the dam, reduced the height of the saturation line, improved the safety factor, and thus reduced the risk of dam failure.

From the modeling perspective, this paper only used the finite element strength reduction method for simulation and analysis, which had certain limitations. In the future, it could consider combining the spatial analysis ability of geographic information system (GIS)41 or use the soft computing technology of gene expression programming (GEP)42,43 for alternating modeling research. In summary, these factors would impact the stability of the tailings dam, which required more in-depth research in the future.

From the actual situation of the tailings dam, the tailings dam was generalized in the modeling process, which did not reflect the distribution of tailings dam materials in the actual state. The simulation effect had a specific error compared with the basic shape of the tailings dam, resulting in the limitation of the numerical simulation. In this study, when analyzing the influence of heightening and reinforcement on the seepage stability of tailings dams, only the effect of dry beach length, heightening, and geogrid reinforcement on the stability of tailings dams under the condition of fluid–solid coupling was studied. When considering the influence of geogrid reinforcement on the stability of tailings dams, only one grid laying method was assessed, and the effect of reinforcement spacing, reinforcement position, and reinforcement length was not considered.

From the point of view of the influence of consolidation on seepage, it was essential to consider the influence of consolidation on seepage in seepage analysis, especially when dealing with fine-grained soil and other materials with certain consolidation properties. Consolidation of tailings may occurred during loading or standing, mainly reflected in the following two aspects: First, the change of permeability coefficient. In the process of consolidation, the pore structure of the tailings changed, thus affecting the permeability of the tailings. The second was the change in pore water pressure. Consolidation resulted in the change of pore water pressure in tailings, which affected the driving force of seepage. In the seepage analysis of reinforced tailings dam, this paper used the finite element method to compare the simulated values of pore water pressure and earth pressure with the monitoring values, reasonably determined the permeability coefficient, and ensured the model’s accuracy. In the future, the corresponding consolidation model could be introduced to further improve the accuracy further.

The seepage stability of the tailings dam could not be ignored. In this paper, ABAQUS software was used to simulate and analyze the influence of heightening and geogrid reinforcement on the stability of the tailings dam under different dry beach lengths. The following conclusions were drawn:

Geogrid played a positive role in reducing the saturation line of a tailings dam. In the case of dam heightening, the saturation line rose by 2.8–5.3 m. Compared with the unreinforced condition (dam heightening), the saturation line of the tailings dam under the local reinforcement condition was reduced by 0.9–2.8 m, and the saturation line of the tailings dam under the overall reinforcement condition was reduced by 3.2–12.5 m.

Geogrid also played a significant role in improving the stability of the tailings dam. In the case of heightening the dam, the safety factor was reduced by 7–15.6%. Compared with the unreinforced condition, the increase rate of the safety factor of the dam body under the local reinforcement condition was between 3.8 and 5.5%, and the increase rate of the safety factor under the overall reinforcement condition was between 35.9 and 42.9%.

The increase in dry beach length improved the stability of the tailings dam. The safety of the tailings dam under normal conditions was much higher than that under minimum dry beach conditions. Under the condition of the same dry beach length, the geogrid effectively improved the stability of the dam body compared with the unreinforced condition.

The fitting degree of the finite element model was good, and the error between the simulation and the experimental value was within 6.9%. This study provided a reference for the later heightening design of the main dam and the design of the new tailings reservoir, which had great reference significance for the design and operation of the reinforced tailings dam project.

All data generated or analysed during this study are included in this published article.

Halabi , ALM , Siacara , AT , Sakano , VK , Pileggi , RG & Futai , MM Tailings dam failures: A historical analysis of the risk .J. Fail.Anal.Prev.22(2), 464–477 (2022).

Rose, R. L., Mugi, S. R. & Saleh, J. H. Accident investigation and lessons not learned: AcciMap analysis of successive tailings dam collapses in Brazil. Reliab. Eng. Syst. Saf. 236, 109308 (2023).

Lyu , Z. , Chai , J. , Xu , Z. , Qin , Y. & Cao , J. A comprehensive review on reasons for tailings dam failures based on case history .Adv.Civil Eng.https://doi.org/10.1155/2019/4159306 (2019).

Zheng, B. et al. Risk evolution study of tailings dam failures disaster based on DEMATEL-MISM. Front. Earth Sci. 10, 924 (2022).

Lin, S. Q. et al. Regional distribution and causes of global mine tailings dam failures. Metals 12(6), 905 (2022).

Güllü, H. & Özel, F. Microtremor measurements and 3D dynamic soil–structure interaction analysis for a historical masonry arch bridge under the effects of near-and far-fault earthquakes. Environ. Earth Sci. 79, 1–29 (2020).

Tolun, M., Emirler, B., Yildiz, A. & Güllü, H. Dynamic response of a single pile embedded in sand including the effect of resonance. Period. Polytech. Civil Eng. 64(4), 1038–1050 (2020).

Lu, M. L. & Cui, L. Analysis of factors affecting seepage field of tailings dam. China Saf. Sci. J. 14(6), 17–20 (2004).

Mafi, R., Javankhoshdel, S., Cami, B., Jamshidi Chenari, R. & Gandomi, A. H. Surface altering optimisation in slope stability analysis with non-circular failure for random limit equilibrium method. Georisk: Assess. Manag. Risk. Eng. Syst. Geohazards 15(4), 260–286 (2021).

Dastpak, P., Jamshidi Chenari, R., Cami, B. & Javankhoshdel, S. Noncircular deterministic and stochastic slope stability analyses and design of simple geosynthetic-reinforced soil slopes. Int. J. Geomech. 21(9), 04021155 (2021).

Aroni Hesari, S., Javankhoshdel, S., Payan, M. & Jamshidi Chenari, R. Pseudo-static internal stability analysis of geosynthetic-reinforced earth slopes using horizontal slices method. Geomech. Geoeng. 17(5), 1417–1442 (2022).

Fatehi, M., Hosseinpour, I., Jamshidi Chenari, R., Payan, M. & Javankhoshdel, S. Deterministic seismic stability analysis of reinforced slopes using pseudo-static approach. Iran. J. Sci. Technol. Trans. Civil Eng. 47(2), 1025–1040 (2023).

Doğan, S. & Güllü, H. Multiple methods for voxel modeling and finite element analysis for man-made caves in soft rock of Gaziantep. Bull. Eng. Geol. Environ. 81, 1–20 (2022).

Wang, X. G., Zhan, H. B., Wang, J. D. & Li, P. The stability of tailings dams under dry-wet cycles: A case study in Luonan, China. Water 10(8), 1048 (2018).

Wang, Q. Research on Safety of Fine Tailings Dam under the Influence of Saturation Line Distribution Characteristics Jin-Shan-Dian (North China University Of Technology, 2013).

Zhang, C. et al. Numerical simulation of seepage and stability of tailings dams: A case study in Lixi, China. Water 12(3), 742 (2020).

Naeini, M. & Akhtarpour, A. Numerical analysis of seismic stability of a high centerline tailings dam. Soil Dyn. Earthq. Eng. 107, 179–194 (2018).

Xiao, H. The Stability Analysis and Reinforcement of Xi-Ye-Shan Tailings Dam in Jin-Shan-Dian (Wuhan University Of Science And Technology, 2010).

Zhou, Y. X. et al. Comprehensive analysis of the stability of tailings-geotextile composite: Iron Mine Tailings Dam in Gushan, Anhui, China. Front. Earth Sci. 10, 931714 (2022).

Li, Q. Y., Ma, G. W. & Lu, Y. L. An experimental and theoretical study on the tailings dam with geotextile bags. Sustainability 15(6), 4768 (2023).

Zheng, B. B., Zhang, D. M., Liu, W. S., Yang, Y. H. & Yang, H. Use of basalt fiber-reinforced tailings for improving the stability of tailings dam. Materials 12(8), 1306 (2019).

Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

Yang, Y. H. et al. Experimental study on dynamic behavior of polyacrylamide-reinforced tailings. Environ. Sci. Pollut. Res. 30(16), 47274–47288 (2023).

Lu, T. et al. Experimental Investigation of sample preparation and grouting technology on microbially reinforced tailings. Construct. Build. Mater. 312, 125458 (2021).

Sun, Y., Gu, X. W. & Xu, X. C. Ecological restoration and mechanical reinforcement effect of slope of tailings reservoir. Environ. Earth Sci. 80, 1–12 (2021).

Güllü, H., Yetim, M. E. & Bacak Güllü, E. Effect of using nano-silica on the rheological, fresh and strength characteristics of cement-based grout for grouting columns. J. Build. Eng. 76, 107100 (2023).

Güllü, H., Yetim, M. E. & Bacak Güllü, E. On the rheological, fresh and strength effects of using nano-silica added geopolymer grout for grouting columns. Eur. J. Environ. Civil Eng. https://doi.org/10.1080/19648189.2023.2245867 (2023).

Güllü, H., Al Nuaimi, M. M. & Aytek, A. Rheological and strength performances of cold-bonded geopolymer made from limestone dust and bottom ash for grouting and deep mixing. Bull. Eng. Geol. Environ. 80, 1103–1123 (2021).

Güllü, H. Comparison of rheological models for jet grout cement mixtures with various stabilizers. Construct. Build. Mater. 127, 220–236 (2016).

Olivier, G., Brenguier, F., de Wit, T. & Lynch, R. Monitoring the stability of tailings dam walls with ambient seismic noise. Lead. Edge 36(4), 350a1-350a6 (2017).

Du , CB , Liang , LD , Yi , F. & Niu , B. Effects of geosynthetic reinforcement on tailings accumulation dams .Water 13(21), 2986 (2021).

Zheng, J. J., Cao, W. Z., Zhou, Y. J. & Jiang, J. G. Pull-out test study of interface behavior between triaxial geogrid and soil. Rock Soil Mech. 38(2), 317–324 (2017).

Lu, Z. M., Chen, C. X., Zuo, B. C. & Huang, C. C. Experimentation research on factors influencing stability of anti-dip layered slope. Rock Soil Mech. 04, 629–632 (2006).

Lei, X. L. Sen Sitivity Study Of Geogrid Reinforcement Mine-Tailings Dam Based on Fem (University Of South China, 2017).

Code for design of tailings facilities (GB50863-2013). Published by China Planning Press. (2013).

Code for seismic design of special structures (GB50191-2012). Published by China Planning Press. (2012).

Cheng, Y. M., Lansivaara, T. & Wei, W. B. Two-dimensional slope stability analysis by limit equilibrium and strength reduction methods. Comput. Geotechn. 34(3), 137–150 (2007).

Zhang, H. Z., Wang, L. G. & Feng, M. S. Stability study of the geogrids reinforced tailings. Adv. Mater. Res. 327, 1–5 (2011).

Islam, K. & Murakami, S. Global-scale impact analysis of mine tailings dam failures: 1915–2020. Glob. Environ. Change 70, 102361 (2021).

Hancock, G. R. & Coulthard, T. J. Tailings dams: Assessing the long-term erosional stability of valley fill designs. Sci. Total Environ. 849, 157692 (2022).

Article  ADS  CAS  PubMed  Google Scholar 

Yu, D. Y., Tang, L. Y., Ye, F. & Chen, C. C. A virtual geographic environment for dynamic simulation and analysis of tailings dam failure. Int. J. Digit. Earth 14(9), 1194–1212 (2021).

Shen, H., Klapperich, H., Abbas, S. M. & Ibrahim, A. Slope stability analysis based on the integration of GIS and numerical simulation. Autom. Construct. 26, 46–53 (2012).

Güllü, H. A new prediction method for the rheological behavior of grout with bottom ash for jet grouting columns. Soils Found. 57(3), 384–396 (2017).

Güllü, H. A novel approach to prediction of rheological characteristics of jet grout cement mixtures via genetic expression programming. Neural Comput. Appl. 28, 407–420 (2017).

The authors would like to thank the reviewers for their constructive comments that improved the paper.

The research was supported by the youth fund of Liaoning Provincial Department of Education (LJKQZ2021153), Liaoning Provincial Department of Education Scientific Research General Project (LJ2020JCL037) and the discipline innovation team of Liaoning Technical University (LNTU20TD-12).

College of Civil Engineering, Liaoning Technical University, Fuxin, 123000, China

Changbo Du, Han Tao & Fu Yi

Beijing Jingneng Geological Engineering Co., Ltd, Beijing, 102300, China

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

C.D. and F.Y. participated in the design of the work, methodology data interpretation, and analysis for the work; carried out the statistical analyses; and drafted the manuscript. H.T. designed the study; participatedin data interpretation, analysis for the work, and methodology; carried out the statistical analyses; and draftedthe manuscript.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Du , C. , Tao , H. & Yi , F. Seepage stability analysis of geogrid reinforced tailings dams .Sci Rep 14, 1814 (2024).https://doi.org/10.1038/s41598-024-52472-y

DOI: https://doi.org/10.1038/s41598-024-52472-y

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Scientific Reports (Sci Rep) ISSN 2045-2322 (online)

Seepage stability analysis of geogrid reinforced tailings dam | Scientific Reports

Woven Geotextile Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.