Abstract

The pore water pressure of tailings dam has a very great influence on the stability of tailings dam. Based on the assumption of one-dimensional consolidation and small strain, the partial differential equation of pore water pressure is deduced. The obtained differential equation can be simplified based on the parameters which are constants. According to the characteristics of the tailings dam, the pore water pressure of the tailings dam can be divided into the slope dam segment, dry beach segment, and artificial lake segment. The pore water pressure is obtained through solving the partial differential equation by separation variable method. On this basis, the dissipation and accumulation of pore water pressure of the upstream tailings dam are analyzed. The example of typical tailings is introduced to elaborate the applicability of the analytic solution. What is more, the application of pore water pressure in tailings dam is discussed. The research results have important scientific and engineering application value for the stability of tailings dam.

1. Introduction

Based on the assumption that the soil is isotropic and uniform, an external surface load is instantaneously applied and is held constant; a classical one-dimensional (1D) consolidation theory was proposed by Terzaghi [1]. In order to analyze time-dependent loading, an analytical solution to the layered consolidation problem for a general set of boundary conditions and an arbitrary load history was presented by Schiffman and Stein [2]. The 1D consolidation analytical solutions considering ramp loading were presented by Olson [3]. A simple semianalytical method to solve the problem of 1D consolidation by taking into account the varied compressibility of soil under cyclic loadings was brought up by Cai et al. [4]. A rigorous solution of the conventional Terzaghi 1D consolidation under haversine cyclic loading with any rest period was proposed by Müthing et al. [5], which is achieved using Fourier harmonic analysis for the periodic function representing the rate of imposition of excess pore water pressure. A semianalytical solution to 1D consolidation of viscoelastic unsaturated soils with a finite thickness under oedometric conditions and subjected to a sudden loading was put forward by Qin et al. [6]. A semianalytical solution to 1D consolidation equation of fractional derivative Kelvin-Voigt viscoelastic saturated soils subjected to different time-dependent loadings was presented by Wang et al. [7]. Under the condition of the increasing weight of superincumbent material and the length of the drainage path varies, a solution for the 1D consolidation of a clay layer whose thickness increases with time was proposed by Gibson [8]. An exact analytical solution of the nonhomogeneous partial differential equation governing the conventional 1D consolidation under haversine repeated loading was derived and discussed by Razouki et al. [9, 10].

In order to analyze different boundary conditions, single drainage solutions for several specific variations of the permeability and shear modulus were given by Mahmoud and Deresiewicz [11]. Several analytical solutions for the consolidation analysis of a soil layer with fairly general laws of variation of permeability and compressibility for both a single-drained condition and a double-drained condition were proposed by Zhu and Yin [12]. Tang et al. [13] propose a closed-form solution for consolidation of three-layered soil with a vertical drain system and a method to solve a convergence solution by enlarging the precision for the entire thickness of the foundation and adding controlling precision for the overall average degree of consolidation of each soil layer. The solution of Terzaghi 1D consolidation equation with general boundary conditions is proposed by Mei and Chen [14], whose solution is validated by comparing it to the classical solution. Hawlader et al. [15] develop a new constitutive model for the compressibility behavior of soft clay sediments at low effective stress level and the model is used to solve finite strain 1D consolidation with pertinent initial and boundary conditions. A solution to the consolidation equation with boundary conditions that are cyclic with time is given by Rahalt and Vuez [16]. A nonlinear theory of consolidation has been developed for an ideal normally consolidated soil by Davis and Raymond [17]. A simple calculation procedure to analyze the one-dimensional response of saturated soil layers to pore pressure variations at the boundary described by a general time-dependent function is developed by Conte and Troncone [18].

The tailings dam is an important geotechnical structure in mining engineering. For a long time, the theory of reservoir dam is applied to tailings dam without any modification. However, there are many differences between tailings dam and reservoir dam, which lead to inaccurate calculation results of the pore water pressure. According to the classical Terzaghi consolidation theory, the analytical solution of the pore water pressure is discussed in this paper.

2. The General Equation of 1D Consolidation

According to Darcy’s law, the flow drag resistance of the -direction is as follows:

The negative sign on the right side of (1) indicates that the drag resistance is opposite to the direction of flow velocity.

According to the mechanical equilibrium conditions, the flow drag resistance of the -direction can also be expressed as follows:

Equations (1) and (2) can be combined as follows:

It is assumed that the water in porosity is incompressible under the consolidation process. Thus, is a constant. During the process of soil consolidation, because the porosity of soil particles is continuously compressed, the permeability coefficient of soil is decreased continuously and the permeability coefficient of soil is changed with depth, that is, . Taking the derivative of (3), one can have

According to fractional derivative rule, (4) can be expanded as follows:

Assume that the horizontal direction consolidation of the soil can be neglected; that is, the consolidation is under confined compression condition. According to the 1D consolidation condition, (6) can be deduced

Under the condition of 1D consolidation, the axial strain is equal to the volume strain of the soil. According to Pane and Schiffman [19] research, (4) and (5) can be rewritten as follows:

In the consolidation problem, it is of great significance to study the change of pore water pressure with time and position under external load. Therefore, (8) can be expanded according to pore water pressure.

Assuming that the soil profile is shown as Figure 1, the pore water pressure of a point located under load and is . According to the Terzaghi effective stress principle, the effective stress could be given as follows under consolidation:

According to Figure 1, the total stress of a point located is as follows:

Substituting (10) into (9) yields

Take the derivative of (11) with respect to and consider the relationship equals

According to the compression curve of 1D consolidation test, the volume compressibility factor can be defined as follows:

Substituting (13) into (14) gives

Substituting (15) into (8), it is noted that the last term of (8) has the following relation:

Equation (8) can be derived as

Substituting (15) into (17), (17) can be expressed as pore water pressure. Equation (18) can be obtained

Equation (18) is a general equation which reflects the 1D consolidation process. It takes into account the situation where the applied load changes with time, the soil thickness changes with time, the soil permeability changes with depth, and so on.

3. Accumulation and Dissipation of Pore Water Pressure in Tailings Dam

Tailings are mine wastes produced in the mining engineering, which are sent to tailings reservoir by pipe or flume. The tailings dam is an important part in mining engineering, which is consisted by the initial dam and fill dam. The initial dam is made by permeable rockfill generally, and the fill dam is formed by tailings. In general, the construction of a tailing dam takes many decades or even a century. According to construction method, the tailings dam can be divided into upstream tailings dam, downstream tailings dam, midline method tailings dam, and so on. The number of tailings ponds in China has reached more than 12,000 by statistics. Because the downstream of the tailings dam is residents living area or mining production area generally the social and people’s property damage a huge impact if the tailings dam is failed.

Because the upstream tailings dam has advantages of simple operation, low construction costs, the less need for coarse particles, and so on, according to statistics, 95% of the tailings dams are adopting the construction of the upstream tailings dam in China. On the other hand, the construction process of the upstream tailings dam cannot precisely control the shape of the tailings dam. The deposition structure of the tailings dam is very complex. The upstream tailings have shortcoming of long infiltration distant and poor stability. In view of the tailings slurry discharge and tailings particle deposition following the sediment mechanics, the section profile of the tailings dam has obvious regularity. The coarser the particle size is, the shorter the average distance migrates and vice versa. Generally, the deposition order of the tailings dam along the dry beach face is tailings fine sand, tailings silt, tailings sand, and tailings mud.

If particle size distribution of the tailings and the length of the dry beach face are kept constant, the interface of the tailings material should be substantially parallel to the surface of the dam slope, which is shown in Figure 2. For most of the tailings dams, because the tailings water contains more heavy metals, the initial topography and the side of tailings dam are treated as an impermeable boundary. On the other hand, the tailings dam extends very long in the direction of the reservoir area. Therefore, the consolidation of the tailings dam is a 1D problem.

As can be seen from Figure 2, the closer to inside the reservoir area is, the thicker the tailings mud is and the worse the average consolidation is. In other words, the closer to inside the reservoir area is, the greater pore water pressure of the tailings dam is and the lower shear strength of the tailings material is. Generally speaking, the construction of a tailings dam needs several decades to reach the design elevation. With the continuous production of mining, the height of the tailings dam is increasing. During this process, the accumulation and dissipation of the excess pore water pressure of the tailings material experience several different stages of development.

3.1. Tailings Thickness Increases with Time

In this stage, the thickness of the tailings increases with time. The tailings material produces great compressive deformation under self-weight effect. The pore water pressure produced by the self-weight pressure at the early stage is partially dissipated. Nevertheless, the pore water pressure caused by self-weight pressure is increased with the tailings thickness at the later stage, which leads to the increase of the pore water pressure with time. Therefore, the net pore water pressure is accumulated during this process.

3.2. Load on the Tailings Layer Increases with Time

The self-deposition of the tailings mud continues to extend to interior of the tailings reservoir area in this stage. The self-deposition of the tailings layer has ended where the end of the tailings reservoir is near. It overlapped with particles coarse tailings and tailings mud. The tailings mud has to bear a growing load with time. The tailings mud continues to consolidate under the combined action of self-weight and additional load. At this stage, the pore pressure produced by the self-weight pressure dissipates a little part, and the more pore water pressure produced by additional load accumulates. The net pore water pressure tends to accumulate in this process.

3.3. The Additional Load of the Tailings Material Remains Constant

For a specific part of the tailings dam, the upper boundary of the tailings dam has reached the designed elevation at a certain moment. Therefore, the additional load will remain unchanged from this moment. During this stage, the tailings mud will continue to consolidate under the combined action of self-weight and constant additional load. Thus, the pore water pressure will dissipate within the tailings dam. Obviously, the different parts of tailings mud have different consolidation stages. There is no additional load at the thickest tailings mud in the tailings reservoir until it is closed. Therefore, the tailings mud is always in the first stage and will never enter the second and third stages. From the above qualitative analysis, the pore water pressure, which is located at the slope dam segment, reaches the maximum value when the additional load just stops growing. The pore water pressure, which is located at dry beach segment and artificial lake segment, achieves the maximum value when the tailings dam reaches the maximum height. The numerical values of the pore water pressure at different locations and moments can be obtained by solving the partial differential equations of consolidation problems.

4. Analytical Solutions of the Pore Water Pressure

It is necessary to accurately calculate the accumulation and dissipation of the pore water pressure in tailings dam, of which a solution of 1D consolidation problem can be simplified. During the consolidation process, the soil parameters of tailings material (such as the change of bulk density, permeability coefficient, and consolidation coefficient) are changed with the consolidation process. The change law can be determined through many experiments. Under normal circumstances, it is difficult to obtain the analytical solution when the change law of soil parameters is considered. Therefore, the exact solution of the problem can only be depended on numerical calculation method.

As the horizontal length of the tailings reservoir is usually far greater than the vertical thickness (i.e., the horizontal length is 10 times more than the vertical thickness), the drainage consolidation effect of the horizontal direction can be ignored. Therefore, the vertical direction of drainage is needed to consider. Then, it is only a 1D consolidation problem. Based on the assumption of small-strain and constant of soil parameters, such as , , , , and are constant, (18) is degenerated into a constant coefficient nonhomogeneous parabolic partial differential equationwhere

The above three different stages are specifically discussed as follows.

4.1. The Thickness Mud of Tailings Mud Increases with Time

The overburden load of the tailings mud at this stage is zero, that is, . The consolidation equation yields

If the bottom of the tailings mud is impermeable, the coordinates origin is taken as the impermeable bottom. The boundary conditions of the problem can be given as follows:

The initial condition of the problem is given by

According to Gibson’s [20] study, (25) is chosen as the solution of (21)

It can be verified that (25) satisfies the boundary condition (22). The choice of function needs to satisfy (23). Substituting (25) into (23) gives

For any function , the partial differential equation (26) has no numerical solution. With regard to , the following parameters can be introduced to transform the equation

Equation (26) can be transformed into

Equation (28) can be obtained through the Laplace change:

Thence

Combining the solutions of (25) and (30), the solution of (21) is obtained

It can be verified that (31) satisfies the initial condition of (24).

4.2. The Additional Load of Tailings Mud Increases with Time

During this stage, the thickness of tailings mud not only does not increase with time but also gradually decreases with the increase of effective stress. For the sake of simplicity, assuming that the thickness of tailings mud is constant, the basic consolidation differential (19) can be simplified as

The initial condition of the problem is given by

The boundary conditions of the problem can be given as

The following three questions are called problem , problem , and problem , respectively:

:

:

:

If the solution of problem , problem , and problem is , , and respectively, it can be directly proved through substitution method:

The solution of belongs to the Terzaghi classical 1D consolidation problem. The analytical solution can be obtained.where

is the solution of (31). Substituting (31) into (40) by Matlab software, can be obtained. Substituting it into (39), the numerical solution can be given.

For solution to the problem , Schiffman [21] acquires the following series solutions:

When the tailings mud is below the groundwater level, the total stress of the tailings mud at the depth can be expressed as follows:

The ratio of pore water pressure to total stress can be easily determined as follows:

4.3. The Additional Load on the Tailings Mud Remains Constant

In this stage, the basic differential equation and the definite condition are problem :

:

The following questions can be called problem :

:

If the solution of problem is and the solution of problem is , the solution of problem can be easily demonstrated is as follows:

The definite condition of problem and problem is exactly the same in form, except that the function is different from the function .

Suppose that

Thenwhere

According to the orthogonal rule of the solution, the following relationship can be obtained:

Equation (48) can be simplified through (50)

Substituting (51) into (47) yields

For the third stage, the analytical solution of pore water pressure is obtained.

Similar to (42), the ratio of pore water pressure to total stress at a certain depth of tailings mud at this stage is as follows:

During this stage, does not change with time, and it can take the following value:

5. Example

The simplified calculation section of the representative section of a tailings dam is shown in Figure 3. The length of the tailings dam is 368 m. Therefore, the problem can be simplified as 1D consolidation. The slime layer below the level is assumed for impervious boundary. The rising speed of the dam height is 1.75 m/y. The saturation density of gravel is 22.0 kN/m3, the saturation density of the tailings fine sand is 19.6 kN/m3, the saturation density of the tailings sand is 19.1 kN/m3, and the saturation density of the tailings mud is 18.8 kN/m3. The permeability coefficient of gravel is 3 × 10−3 m/s, the permeability coefficient of the tailings fine sand is 5 × 10-6 m/s, the permeability coefficient of the tailings sand is 2 × 10-6 m/s, and the permeability coefficient of the tailings mud is 1.5 × 10-8 m/s. The consolidation coefficient of gravel is 3.8 × 10-1 cm2/s, the consolidation coefficient of the tailings fine sand is 5.4 × 10-3 cm2/s, the consolidation coefficient of the tailings sand is 2.8 × 10-3 cm2/s, and the consolidation coefficient of the tailings mud is 1.2 × 10-5 cm2/s. Please determine the distribution of pore water pressure when the tailings dam reaches maximum height.

The above problem can be solved through (31), (46), and (52). The calculation can be accomplished through the MATLAB software. Firstly, a coordinate system can be set up. Then, the corresponding coordinates are taken out of the interval 1 m. The region is determined according to the coordinates. Then, the pore water pressure can be calculated through the corresponding equation. Finally, the contour map can be drawn through a set of array. The results of the calculation are shown in Figure 4.

6. Discussion

Due to the similarity between the tailings dam and the reservoir dam in the geotechnical structure, many scholars directly introduce the calculation methods of the reservoir dam, whose theory is relatively mature, into the tailings dam. Considering the difference on construction cycle, construction materials between the tailings dam and reservoir dam, it makes the calculation results inconsistent with the actual.

Based on the Terzaghi consolidation theory of 1D, the tailings dam is divided into the slope dam segment, the dry beach segment, and the artificial lake segment. The solutions of the pore water pressure are derived, respectively. The analysis shows that the additional load of the slope dam segment is unchanged, which can be calculated using (46). The analysis indicates that tailings thickness of the artificial lake segment increases with time during construction period, which can be calculated by (31). The situation of the dry beach section is slightly complicated. The additional load of the dry beach section which is under the slope dam segment is basically the same, which can be calculated by (52). The additional load on the other sections of the dry beach segment increases with time, which can be calculated using (31).

The theoretical derivation is based on 1D consolidation theory. It is only considering consolidation in the vertical direction. Since the horizontal scale of most tailings dams is much larger than the vertical direction, it has little influence to ignore the drainage of horizontal. From the point of engineering view, it is conservative to the stability of the tailings dam. Considering drainage of the horizontal, the pore water pressure will be lessened and the safety factor of the tailings dam will be greater.

It is assumed that the deformation of the tailings is small deformation during the consolidation process. If the actual tailings are loose relatively, the deformation of the tailings is large deformation. He et al. [22] studied the pore water pressure of saturated soils using the updated Lagrangian formulation of large strain method (ULM), total Lagrangian formulation of large strain method (TLM), and small strain method (SSM). The results point out that the pore water pressure of the above three methods is exactly the same when the strain is less than 2.4%. When the strain is greater than 2.4%, the pore water pressure of the ULM and SSM still keeps the same, and the pore water pressure of the TLM is obviously less than the ULM and SSM. When the strain reaches 7.5%, the pore water pressure of the ULM and SMM began to bifurcate. When the deformation of the tailings dam is examined for a relatively short time, the consolidation of prophase tailings dam has been completed under previous load. The consolidation deformation of the whole tailings dam is relatively small, and its deformation conforms to the small deformation assumption. On the other hand, it is difficult to calculate and solve partial differential equations using large deformation theory. The differences between large deformation and small deformation need to be further studied.

It is supposed that the mechanical parameters such as permeability coefficient and consolidation coefficient are constant during the consolidation process. Previous studies [23] have shown that the permeability coefficient of rock and soil decreases with the increase of deformation during the consolidation process. It is believed that the value of decides whether it is necessary to take under consideration the effect of nonlinear property by Zhuang et al. [24]. The results show that pore water pressure calculated depends on the ratio of loading intensity to the initial effective vertical stress subjected to time-dependent loading by Conte and Troncone [25]. How the permeability coefficient of tailings material changes during the consolidation process is less studied by relevant scholars. How the permeability coefficient of tailings material changes during the consolidation process is less studied by relevant scholars. From the aspect of engineering view, the larger the pore water pressure of actual is, the smaller the safety factor of tailings dam is. Therefore, it is dangerous for the tailings dam. The law of permeability coefficient of tailing material during consolidation process needs to be further studied.

7. Conclusion

Tailings dam is a very important geotechnical structure of mine engineering. The calculation of pore water pressure has a great impact on the safety factor of tailings dam slope. How to accurately estimate pore water pressure is very difficult. Based on the assumption of 1D consolidation and small strain of tailings material, a general equation of the pore water pressure is proposed. According to dissipation and accumulation characteristics of the pore water pressure in the tailings dam, the tailings dam can be divided into the slope dam segment, the dry beach segment, and the artificial lake segment. The analytic solutions of the corresponding segment are obtained through solving the partial differential equation, which has great significance to the stability of the tailings dam.

Notations

:Shear strength of the tailings material
:Maximum increases in vertical total stress as a function of depth
:A parameter introduced to transform an equation
:The rate of the additional load on the tailings mud,
:Coordinate of direction
:A parameter introduced to transform an equation
:Porosity of the tailings material
:Counters, 1, 2, 3,…
:A function of variable
:A function introduced into solving partial differential equation
:The rate of the thickness of the tailings mud,
:A function introduced to transform an equation
:Stress of the tailings material
:Coordinate of direction
:The pore water pressure
:Coefficient of coefficient
:The height of static water level
:Thickness of tailings mud
:Porosity ratio
:Increment of stress
:Termination pore water pressure as a function of depth
:Initial pore water pressure as a function of depth
:Termination time
:Initial time
:Coefficient of consolidation
:Saturation
:Coefficients to be determined
:Coefficients to be determined
:Initial total stress
:Effective vertical stress
:Strain in the direction
:Bulk density of water
:Saturated bulk density of tailings dam
:Buoyancy unit weight of tailings material
:The actual velocity along the flow direction in the tailings dam
:Water pressure
:Coefficient of volume compressibility
:Hydraulic gradient in the direction
:The total height of tailings mud plus static water level
:The drag resistance force on the pore wall of a unit volume in the direction.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This paper is supported by the Open Research Fund of State Key Laboratory of Geomechanics and Geotechnical Engineering, Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, (Grant no. Z013009), the National Key Research and Development Program of China (Project no. 2017YFC0804601), the National Natural Science Foundation of China (Grant nos. 51764020; 51741410; 51234004), and the Natural Science Foundation of Yunnan Province (Grant no. 2015FB130). The authors would like to thank them for providing the financial support for conducting this research.