# Time-Varying Hydraulic Gradient Model of Paste-Like Tailings in Long-Distance Pipeline Transportation

Authors: Li Yang; Shuai Li; Xin-min Wang
Journal: Advances in Materials Science and Engineering
 (2017)

---

## Abstract

Paste-like tailings slurry (PTLS) is always simplified as a Bingham plastic fluid, leading to excessive computational errors in the calculation of the hydraulic gradient. In the case of paste-like tailings in long-distance pipeline transportation, to explore a high-precision and reliable hydraulic gradient formula, the rheological behavior of paste-like tailings slurry was analyzed, a time-varying hydraulic gradient model was constructed, and a series of laboratory shear tests were conducted. The results indicate that the PTLS shows noticeable shear-thinning characteristics in constant shear tests; the calculated hydraulic gradient declined by about 56%, from 4.44 MPa·km−1 to 1.95 MPa·km−1 within 253 s, and remained constant for the next four hours during the pipeline transportation. Comparing with the balance hydraulic gradient obtained in a semi-industrial loop test, the computational errors of those calculated by using the time-varying hydraulic gradient model, Jinchuan formula, and Shanxi formula are 15%, 78%, and 130%, respectively. Therefore, our model is a feasible and high-precision solution for the calculation of the hydraulic gradient of paste-like tailings in long-distance pipeline transportation.


---

## Body

## 1. Introduction

In the last twenty years, paste technology has been widely accepted as an environment-friendly, economical, and effective way in tailings disposal and stope backfilling [1]. Paste is dewatered tailings with little or no water bleed, which is nonsegregating in nature (Figure 1). It can be disposed on the surface with less chance of failure, lower seepage probability, and a smaller impact scope compared to a conventional impoundment [2, 3].Figure 1
Comparison of different tailings surface disposal methods.Table1 summarizes several typical examples of paste tailings disposal mines in the world [4]. Evidently, many mines’ tailings are actually disposed as “paste-like” rather than as designed paste tailings. The primary reason is that tailings are finely ground before mineral extraction, resulting in the tailings slurry being abundant in slow-settling particles such as electronegative slimes, ultrafine tailings, and beneficiation reagents [5]. Therefore, it is difficult to achieve a paste level with an underflow concentration that varies with the tailings’ sizes and properties [6].Table 1
Application examples of paste tailings disposal mines.
Mines name and locations
Application environment
Tailings concentration
Capacity (t/d)
Jinshandian Iron Mine, Hubei, China
Southern rainy regions
62% (Paste-like)
4000
Wushan Cu-Mo Mine, Inner Mongolia, China
Arid and cold regions
65% (Paste-like)
30000
Musselwhite Mine, Ontario, Canada
Cold and flat regions
65% (Paste-like)
4000
Toromocho Mine, Junin, Peru
Arid and high earthquake intensity zone
67% (Paste-like)
117000
Esperanza Mine, Coquimbo, Chile
Arid and high earthquake intensity zone
69% (Paste)
20000
Garpenberg Mine, Stockholm, Sweden
Cold regions
70% (Paste)
—In previous studies, paste-like tailings slurry (PLTS) was simplified as a Bingham plastic fluid, whose flow properties were time-independent, and the existing empirical formulas for calculating the hydraulic gradient of a paste-like flow were derived based on the time-invariant Bingham model [7]. However, an increasing number of studies have shown that the rheological characteristics of PLTS are consistent with those of a pseudo-plastic fluid, and the existing empirical formulas are optimal for calculating the error and for small applications [8]. Pornillos [9] found that the yield stress of thickened tailings slurry showed distinct time-varying behavior under a continuous shear stress. Senapati and Mishra [10] evaluated the head loss of paste-like coal combustion slurry using the power-law-fluid-based head loss model and got results with a higher precision than that with empirical formulas. Li et al. [11] calculated the hydraulic gradient of cemented paste-like backfill in long pipeline transportation and found that the design error reached as high as 40% using empirical formulas. Therefore, avoiding the simplification of PLTS as a Bingham plastic fluid and deriving a time-varying hydraulic gradient model for PLTS is of great significance.This paper discusses the disposal of paste-like tailings in the largest underground iron mine in Asia and presents the time-varying behavior of PLTS by employing groups of shear tests. A time-varying hydraulic gradient model of PLTS in long-distance pipeline transportation was derived based on the theory of flocculent structure and the Herschel–Bulkley model. The reliability of the time-varying hydraulic gradient model was analyzed and compared with a semi-industrial loop test and conventional semiempirical formulas.


## 2. Flocculent Structures (FS) and Rheological Model of PLTS

### 2.1. Preparation of PLTS

The Sijiaying Iron Mine in Hebei, China, successfully commissioned in 2007, is the largest underground iron mine in Asia to adopt a total paste-like solution for all of its tailings. About 40 million tons/year of argillized ultrafine tailings slurry will be used for stopes backfilling, and the remainder will be “stacked” on the surface. Figure2 presents a flow chart of the preparation and transportation of PLTS in Sijiaying.Figure 2
The preparation and transportation of PLTS in Sijiaying.As part of a widely accepted practice, efficient tailings-thickening equipment, such as deep-cone thickeners with a diameter of 25 m, is used in Sijiaying to obtain PLTS. After being uniformly mixed with flocculants, the sedimentation of argillized ultrafine tailings is accelerated, and the tailings concentration can be improved to 65%. With a solids content less than 300 ppm, the overflow water can be reused for flotation [12]. The tailings disposal storage, which is located in the subsidence of the Fangezhuang Coal Mine, is 30 km from Sijiaying. The paste-like underflow is transported by using piston pumps in a long-distance pipeline and then dumped and evaporated. There are several benefits cited by Sijiaying in their decision to develop the paste-like tailings disposal system: improving the dam safety in floods and earthquakes, water recycling, and decreasing the footprint of tailings disposal sites [13].


### 2.2. FS of PLTS

The tailings slurry in the highly oxidized and argillized hematite mine of Sijiaying is abundant in electronegative slimes and ultrafine tailings. The cationic polyacrylamide (CPAM) with a molecular weight of 18 million exhibited good flocculation and sedimentation after a series of flocculant selection tests. Charge neutralization and bridging constitute the main mechanisms of CPAM [14]. The negatively charged slimes and tailings are neutralized by cations after the CPAM is uniformly mixed, which weakens the strong repulsive interactions of the tailings. Consequently, the thick hydration shells get thinner. The long molecular chains from CPAM will release linearly structured reactive groups after the process of hydrolysis and polycondensation; they then experience mutual attraction by the formation of multiple colloidal bonds, capturing and gathering through bridging, and finally developing into FS [15]. FS, which are unstable, can not only be stretched or damaged by mechanical mixing, but also lapped or recombined by bridging (Figure 3) [16].Figure 3
The dynamic damage and recovery process of FS.This dynamic damage and recovery process of FS can be described as(1)dSdt=k11-S-k2Sγ˙,where S is the FS coefficient, reflecting the degree of integrity of FS; r is shear rate expressed in s−1; k1 is the recovery coefficient, which depends on the flocculant dosage; and k2 is the damage coefficient, directly proportional to r. S=0 is defined as the stage where the FS are completely destroyed, while S=1 is the stage where the FS are fully developed. When (1) = 0, it means that the damage and recovery processes have reached a dynamic balance, and the FS balance coefficient S0 can be obtained as(2)S0=k1k1+k2γ˙.After integrating (1), the following equation can be obtained:(3)S=k1+k2γ˙e-k1+k2γ˙tk1+k2γ˙=S0+k2γ˙k1+k2γ˙e-k1+k2γ˙t.


### 2.3. Rheological Model of PLTS

Non-Newtonian fluids include Bingham plastic fluid, pseudo-plastic fluid, and expansion fluid, depending on different rheological properties [17]. To validate the rheological model of PLTS, shear stresses with different shear rates were tested. The results show that the behavior of PLTS is consistent with that of a pseudo-plastic fluid, and the Herschel–Bulkley model is more applicable to PLTS. The constitutive equation of the Herschel–Bulkley model is commonly written as(4)τ=τ0+ηγ˙n,where τ is the shear stress in Pa; τ0 is the yield stress in Pa; η is the apparent viscosity in Pa·s; and n is the flow index, n<1.As shown in (3), the FS show noticeable shear-thinning characteristics under a particular shear rate, as S decreases gradually with the increase in shear time and finally reaches the balanced value S0. The time-varying behavior of PLTS can be further manifested in the change laws between the parameters of the Herschel–Bulkley model and S, which is reflected in the mathematical model given by(5)τ0=Sτ1,η=η∞+Sk3,n=Sn0,where τ1 and n0 are the initial yield stress and initial flow index, respectively; η∞ is the apparent viscosity when the FS are completely destroyed; and k3 is the viscosity coefficient of the FS. The rheological model of PLTS based on the flocculent structure theory can be obtained by integrating (4) and (5):(6)τ=Sτ1+η∞+Sk3γSn0=fτ1,η∞,n0,k1,k2,k3,γ˙,t.


### 2.4. Time-Varying Hydraulic Gradient Model of PLTS

The hydraulic gradient, or resistance loss, is the most important parameter in long-distance pipeline transportation, affecting the operation consumption and reliable operation of a PLTS disposal system. In the present study, the hydraulic gradient formula of PLTS was derived based on the time-invariant Bingham model:(7)i0=163Dτ0+32vD2η,where i0 is the hydraulic gradient in Pa/m; D is the pipe diameter in m; and v is the average flow velocity in m/s.There is usually a long pipe-feeding distance before the PLTS arrives at the storage, and the continuous shearing force provided by the pipe wall can last several hours. Under the action of a continuous shearing force, the PLTS clearly shows time-varying behavior, resulting in a large error in the hydraulic gradient calculation. Therefore, the hydraulic gradient formula of PLTS should be modified by using the time-varying flocculent structure theory and the Herschel–Bulkley model. The time-varying hydraulic gradient model of PLTS can then be obtained by integrating (3), (5), and (7):(8)τ0=Sτ1=S0τ1+ae-k1+k2γ˙t,a=τ1k2γ˙k1+k2γ˙η=η∞+Sk3=η∞+S0k3+be-k1+k2γ˙t,b=k2k3γ˙k1+k2γ˙i0=163DS0τ1+32vD2η∞+S0k3+163Da+32vD2be-k1+k2γ˙t,where a and b are dimensionless parameters in the simplified formula.From (8), it can be seen that the modified τ0, η, and i0 of the PLTS are all decreasing functions of time. Under the action of a continuous shearing force, the FS are destroyed gradually, resulting in a decline in the rheological properties.


## 2.1. Preparation of PLTS

The Sijiaying Iron Mine in Hebei, China, successfully commissioned in 2007, is the largest underground iron mine in Asia to adopt a total paste-like solution for all of its tailings. About 40 million tons/year of argillized ultrafine tailings slurry will be used for stopes backfilling, and the remainder will be “stacked” on the surface. Figure2 presents a flow chart of the preparation and transportation of PLTS in Sijiaying.Figure 2
The preparation and transportation of PLTS in Sijiaying.As part of a widely accepted practice, efficient tailings-thickening equipment, such as deep-cone thickeners with a diameter of 25 m, is used in Sijiaying to obtain PLTS. After being uniformly mixed with flocculants, the sedimentation of argillized ultrafine tailings is accelerated, and the tailings concentration can be improved to 65%. With a solids content less than 300 ppm, the overflow water can be reused for flotation [12]. The tailings disposal storage, which is located in the subsidence of the Fangezhuang Coal Mine, is 30 km from Sijiaying. The paste-like underflow is transported by using piston pumps in a long-distance pipeline and then dumped and evaporated. There are several benefits cited by Sijiaying in their decision to develop the paste-like tailings disposal system: improving the dam safety in floods and earthquakes, water recycling, and decreasing the footprint of tailings disposal sites [13].


## 2.2. FS of PLTS

The tailings slurry in the highly oxidized and argillized hematite mine of Sijiaying is abundant in electronegative slimes and ultrafine tailings. The cationic polyacrylamide (CPAM) with a molecular weight of 18 million exhibited good flocculation and sedimentation after a series of flocculant selection tests. Charge neutralization and bridging constitute the main mechanisms of CPAM [14]. The negatively charged slimes and tailings are neutralized by cations after the CPAM is uniformly mixed, which weakens the strong repulsive interactions of the tailings. Consequently, the thick hydration shells get thinner. The long molecular chains from CPAM will release linearly structured reactive groups after the process of hydrolysis and polycondensation; they then experience mutual attraction by the formation of multiple colloidal bonds, capturing and gathering through bridging, and finally developing into FS [15]. FS, which are unstable, can not only be stretched or damaged by mechanical mixing, but also lapped or recombined by bridging (Figure 3) [16].Figure 3
The dynamic damage and recovery process of FS.This dynamic damage and recovery process of FS can be described as(1)dSdt=k11-S-k2Sγ˙,where S is the FS coefficient, reflecting the degree of integrity of FS; r is shear rate expressed in s−1; k1 is the recovery coefficient, which depends on the flocculant dosage; and k2 is the damage coefficient, directly proportional to r. S=0 is defined as the stage where the FS are completely destroyed, while S=1 is the stage where the FS are fully developed. When (1) = 0, it means that the damage and recovery processes have reached a dynamic balance, and the FS balance coefficient S0 can be obtained as(2)S0=k1k1+k2γ˙.After integrating (1), the following equation can be obtained:(3)S=k1+k2γ˙e-k1+k2γ˙tk1+k2γ˙=S0+k2γ˙k1+k2γ˙e-k1+k2γ˙t.


## 2.3. Rheological Model of PLTS

Non-Newtonian fluids include Bingham plastic fluid, pseudo-plastic fluid, and expansion fluid, depending on different rheological properties [17]. To validate the rheological model of PLTS, shear stresses with different shear rates were tested. The results show that the behavior of PLTS is consistent with that of a pseudo-plastic fluid, and the Herschel–Bulkley model is more applicable to PLTS. The constitutive equation of the Herschel–Bulkley model is commonly written as(4)τ=τ0+ηγ˙n,where τ is the shear stress in Pa; τ0 is the yield stress in Pa; η is the apparent viscosity in Pa·s; and n is the flow index, n<1.As shown in (3), the FS show noticeable shear-thinning characteristics under a particular shear rate, as S decreases gradually with the increase in shear time and finally reaches the balanced value S0. The time-varying behavior of PLTS can be further manifested in the change laws between the parameters of the Herschel–Bulkley model and S, which is reflected in the mathematical model given by(5)τ0=Sτ1,η=η∞+Sk3,n=Sn0,where τ1 and n0 are the initial yield stress and initial flow index, respectively; η∞ is the apparent viscosity when the FS are completely destroyed; and k3 is the viscosity coefficient of the FS. The rheological model of PLTS based on the flocculent structure theory can be obtained by integrating (4) and (5):(6)τ=Sτ1+η∞+Sk3γSn0=fτ1,η∞,n0,k1,k2,k3,γ˙,t.


## 2.4. Time-Varying Hydraulic Gradient Model of PLTS

The hydraulic gradient, or resistance loss, is the most important parameter in long-distance pipeline transportation, affecting the operation consumption and reliable operation of a PLTS disposal system. In the present study, the hydraulic gradient formula of PLTS was derived based on the time-invariant Bingham model:(7)i0=163Dτ0+32vD2η,where i0 is the hydraulic gradient in Pa/m; D is the pipe diameter in m; and v is the average flow velocity in m/s.There is usually a long pipe-feeding distance before the PLTS arrives at the storage, and the continuous shearing force provided by the pipe wall can last several hours. Under the action of a continuous shearing force, the PLTS clearly shows time-varying behavior, resulting in a large error in the hydraulic gradient calculation. Therefore, the hydraulic gradient formula of PLTS should be modified by using the time-varying flocculent structure theory and the Herschel–Bulkley model. The time-varying hydraulic gradient model of PLTS can then be obtained by integrating (3), (5), and (7):(8)τ0=Sτ1=S0τ1+ae-k1+k2γ˙t,a=τ1k2γ˙k1+k2γ˙η=η∞+Sk3=η∞+S0k3+be-k1+k2γ˙t,b=k2k3γ˙k1+k2γ˙i0=163DS0τ1+32vD2η∞+S0k3+163Da+32vD2be-k1+k2γ˙t,where a and b are dimensionless parameters in the simplified formula.From (8), it can be seen that the modified τ0, η, and i0 of the PLTS are all decreasing functions of time. Under the action of a continuous shearing force, the FS are destroyed gradually, resulting in a decline in the rheological properties.


## 3. Laboratory Shear Tests

### 3.1. Properties of Tailings

The properties of the tailings and the particle size are the most important parameters affecting the rheological model and FS of the PLTS. The properties of the tailings of Sijiaying are summarized in Table2, where they are classified as ultrafine, with 70% of the content smaller than 37 μm and 8.9% larger than 74 μm. As this mine is a highly oxidized and argillized hematite mine, the tailings slurry is abundant in electronegative slimes and ultrafine tailings. Since like charges repel, the strong electrostatic repulsion between the slow-settling particles causes a large flocculant dosage and a paste-like underflow from the deep-cone thickener. The particle size of the slow-settling particles is denoted by D (Table 3); the content less than 6.5 μm in size exceeds 71%, which will attract nearby polar water molecules and form thick hydration shells blocking the flocculation and sedimentation processes.Table 2
Main tailings properties of Sijiaying.
Propertiesindexes
pHvalue
Porosity/%
Bulk density/t⋅m−3
Density/t⋅m−3
Averagesize/μm
Content of<37 μm/%
Content of>74 μm/%
Value
7-8
55.66
1.40
3.09
19
70.02
8.9Table 3
Particle sizes of slow-settling particles.
D/μm
D < 6.5
18 >D > 6.5
26 >D > 18
36 >D > 26
D > 36
Content/%
71.20
15.15
5.20
3.05
5.40


### 3.2. Experimental Procedure

The HAAKE VT550 rotary viscometer was used for laboratory shear tests. The immersion sensor, whose cross-shaped rotor overcomes the slip effect of the PLTS, has higher measurement accuracy than a traditional coaxial cylinder. The control software, which is installed in the computer, can easily handle and record slight changes in the rheological parameters during shear tests. The PLTS of Sijiaying with a mass concentration of 65% was prepared in advance, rapidly stirred using an electric mixer, and kept stagnant for about 5 min. Subsequently, 900 g of PLTS was taken in a 500 mL beaker, and placed under the viscometer. The initial yield stress was defined as the shear stress when the rotor began rolling. Constant shear tests were conducted to measure the apparent viscosity and shear stress of the PLTS in five groups, with shear rates within 600 s.


### 3.3. Results and Analyses

The variation in the shear stress and apparent viscosity with the shear time are shown in Figure4. The shear rate is set as 20 s−1, 30 s−1, 40 s−1, 50 s−1, and 60 s−1.Figure 4
Varying tendency of rheological parameters with shear time: (a) shear stress; (b) apparent viscosity.
(a)
(b)The PLTS of Sijiaying clearly showed shear-thinning characteristics in the constant shear tests conducted in the laboratory. When the shear time was zero, the shear stress and apparent viscosity of the PLTS had the highest values, which meant that the FS were fully developed. Under a particular shear rate, the FS of PLTS were destroyed gradually, resulting in a rapid decline in the shear stress and apparent viscosity. The balance time was defined as the time when the shear stress and apparent viscosity acquired steady values. The calculated balance times of the five groups with different shear rates were 381 s, 324 s, 307 s, 267 s, and 230 s, respectively. The damage process was accelerated with the increase in shear rate. The greater the shear rate, the shorter the balance time.


### 3.4. Error Analysis

All the balanced rheological parameters of the PLTS in the constant shear tests conducted in the laboratory are summarized in Table4. A nonlinear equation can be obtained by substituting the data in (6). After a series of arithmetic operations, an optimal solution was obtained by using the damped least-squares method:(9)τ=fτ1,η∞,n0,k1,k2,k3,γ˙,t=f118.5,1.48,0.99,6.95×10-3,3.6×10-4,1.85,γ˙,t.Table 4
The balance rheological parameters of PLTS.
Shear rate/s−1
20
30
40
50
60
Balance apparent viscosityηb/(Pa⋅s)
2.76
2.55
2.26
2.02
1.76
Balance shear timetb/s
381
324
307
267
230
Balance shear stressτb/Pa
100.2
111.5
121.9
132.8
144.7To verify the accuracy of (9), a brief comparison of the calculated shear stress and the tested values is shown in Figure 5. With the shear rates being 20, 40, and 60 s−1, the calculated shear stresses are consistent with the tested values, implying that the rheological model of the PLTS based on the theory of flocculent structure and the Herschel–Bulkley model does have a statistical significance.Figure 5
Comparison of the calculated shear stress and the tested values.


## 3.1. Properties of Tailings

The properties of the tailings and the particle size are the most important parameters affecting the rheological model and FS of the PLTS. The properties of the tailings of Sijiaying are summarized in Table2, where they are classified as ultrafine, with 70% of the content smaller than 37 μm and 8.9% larger than 74 μm. As this mine is a highly oxidized and argillized hematite mine, the tailings slurry is abundant in electronegative slimes and ultrafine tailings. Since like charges repel, the strong electrostatic repulsion between the slow-settling particles causes a large flocculant dosage and a paste-like underflow from the deep-cone thickener. The particle size of the slow-settling particles is denoted by D (Table 3); the content less than 6.5 μm in size exceeds 71%, which will attract nearby polar water molecules and form thick hydration shells blocking the flocculation and sedimentation processes.Table 2
Main tailings properties of Sijiaying.
Propertiesindexes
pHvalue
Porosity/%
Bulk density/t⋅m−3
Density/t⋅m−3
Averagesize/μm
Content of<37 μm/%
Content of>74 μm/%
Value
7-8
55.66
1.40
3.09
19
70.02
8.9Table 3
Particle sizes of slow-settling particles.
D/μm
D < 6.5
18 >D > 6.5
26 >D > 18
36 >D > 26
D > 36
Content/%
71.20
15.15
5.20
3.05
5.40


## 3.2. Experimental Procedure

The HAAKE VT550 rotary viscometer was used for laboratory shear tests. The immersion sensor, whose cross-shaped rotor overcomes the slip effect of the PLTS, has higher measurement accuracy than a traditional coaxial cylinder. The control software, which is installed in the computer, can easily handle and record slight changes in the rheological parameters during shear tests. The PLTS of Sijiaying with a mass concentration of 65% was prepared in advance, rapidly stirred using an electric mixer, and kept stagnant for about 5 min. Subsequently, 900 g of PLTS was taken in a 500 mL beaker, and placed under the viscometer. The initial yield stress was defined as the shear stress when the rotor began rolling. Constant shear tests were conducted to measure the apparent viscosity and shear stress of the PLTS in five groups, with shear rates within 600 s.


## 3.3. Results and Analyses

The variation in the shear stress and apparent viscosity with the shear time are shown in Figure4. The shear rate is set as 20 s−1, 30 s−1, 40 s−1, 50 s−1, and 60 s−1.Figure 4
Varying tendency of rheological parameters with shear time: (a) shear stress; (b) apparent viscosity.
(a)
(b)The PLTS of Sijiaying clearly showed shear-thinning characteristics in the constant shear tests conducted in the laboratory. When the shear time was zero, the shear stress and apparent viscosity of the PLTS had the highest values, which meant that the FS were fully developed. Under a particular shear rate, the FS of PLTS were destroyed gradually, resulting in a rapid decline in the shear stress and apparent viscosity. The balance time was defined as the time when the shear stress and apparent viscosity acquired steady values. The calculated balance times of the five groups with different shear rates were 381 s, 324 s, 307 s, 267 s, and 230 s, respectively. The damage process was accelerated with the increase in shear rate. The greater the shear rate, the shorter the balance time.


## 3.4. Error Analysis

All the balanced rheological parameters of the PLTS in the constant shear tests conducted in the laboratory are summarized in Table4. A nonlinear equation can be obtained by substituting the data in (6). After a series of arithmetic operations, an optimal solution was obtained by using the damped least-squares method:(9)τ=fτ1,η∞,n0,k1,k2,k3,γ˙,t=f118.5,1.48,0.99,6.95×10-3,3.6×10-4,1.85,γ˙,t.Table 4
The balance rheological parameters of PLTS.
Shear rate/s−1
20
30
40
50
60
Balance apparent viscosityηb/(Pa⋅s)
2.76
2.55
2.26
2.02
1.76
Balance shear timetb/s
381
324
307
267
230
Balance shear stressτb/Pa
100.2
111.5
121.9
132.8
144.7To verify the accuracy of (9), a brief comparison of the calculated shear stress and the tested values is shown in Figure 5. With the shear rates being 20, 40, and 60 s−1, the calculated shear stresses are consistent with the tested values, implying that the rheological model of the PLTS based on the theory of flocculent structure and the Herschel–Bulkley model does have a statistical significance.Figure 5
Comparison of the calculated shear stress and the tested values.


## 4. Hydraulic Gradient Calculation and Reliability Analysis

### 4.1. Hydraulic Gradient Calculation

The PLTS of Sijiaying with a mass concentration of 65% is transported by piston pumps in a 30 km long pipeline and then dumped and evaporated. Since the capacity of PLTS disposal is about 30 million tons/year, transporting such a large amount of PLTS in such a long-distance pipeline is proving to be particularly challenging and difficult. Therefore, an accurately calculated hydraulic gradient is especially important in reducing the operation consumption and improving the reliability of the PLTS disposal system. The flow rateQ was set to 400, 450, 500, and 550 m3·h−1, and the diameter of the pipeline was 300 mm. The hydraulic gradient of the PLTS was calculated by using (8) and is given in Table 5.Table 5
Calculated hydraulic gradient of PLTS.
Flow rateQ/(m3⋅h−1)
Flow velocityv/(m⋅s−1)
Shear rater/s-1
Dimensionless parameters
Balance timetb/s
Maximum ofhydraulic gradient imax/(MPa⋅km−1)
Balance hydraulic gradientib/(MPa⋅km−1)
a
b
400
1.57
41.9
81.15
1.27
286
3.97
1.82
450
1.77
47.2
84.09
1.31
270
4.20
1.88
500
1.97
52.4
86.61
1.35
253
4.44
1.95
550
2.16
57.7
88.78
1.39
241
4.67
2.02It is evident from Table5 that the PLTS shows a time-varying behavior during the long-distance pipeline transportation, with the hydraulic gradient declining gradually and finally reaching equilibrium. When the flow rate was 500 m3·h−1, the calculated hydraulic gradient of the PLTS declined by about 56%, from 4.44 MPa·km−1 to 1.95 MPa·km−1 within 253 s. With an increase in the flow rate, the balance time becomes shorter while the balance hydraulic gradient becomes larger. With an average flow velocity of 1.97 m·s−1, it takes approximately four hours to transport the PLTS to the storage. Therefore, it is more systematic and reasonable to use the balance hydraulic gradient for the pipeline design and pump selection. Besides, it is advisable to stir the PLTS thoroughly before transportation, to destroy the FS completely, and reduce the hydraulic gradient effectively, based on the time-varying behavior of the PLTS.


### 4.2. Reliability Analysis

The conventional semiempirical hydraulic gradient formulas, which were derived based on a two-phase flow laboratory date, show a limited applicable scope [18]. For example, Durand’s formula may be suitable for the short-distance transportation of coarse and low-concentration crude tailings. As a widely accepted formula in China, the Shanxi Shuiliyuan formula is suitable for coarse particles, and the Jinchuan formula gives a high precision in calculations related to high-concentration rod-milling sand [19].To evaluate the reliability of the calculated hydraulic gradient by using (8), a semi-industrial loop test was conducted. The hydraulic gradients of the PLTS were calculated by recording the press values of the observation points (Figure 6). A comparison of the hydraulic gradients, calculated using semiempirical formulas, loop test, and (8), is shown in Figure 7. Comparing with the balance hydraulic gradient obtained in the semi-industrial loop test, the computational errors of those calculated by using the time-varying hydraulic gradient model, Jinchuan formula, and Shanxi formula are 15%, 78%, and 130%, respectively. Therefore, the time-varying hydraulic gradient model is more feasible and reliable than semiempirical formulas.Figure 6
Layout of the loop system.Figure 7
Comparison of calculated hydraulic gradient using different methods.


## 4.1. Hydraulic Gradient Calculation

The PLTS of Sijiaying with a mass concentration of 65% is transported by piston pumps in a 30 km long pipeline and then dumped and evaporated. Since the capacity of PLTS disposal is about 30 million tons/year, transporting such a large amount of PLTS in such a long-distance pipeline is proving to be particularly challenging and difficult. Therefore, an accurately calculated hydraulic gradient is especially important in reducing the operation consumption and improving the reliability of the PLTS disposal system. The flow rateQ was set to 400, 450, 500, and 550 m3·h−1, and the diameter of the pipeline was 300 mm. The hydraulic gradient of the PLTS was calculated by using (8) and is given in Table 5.Table 5
Calculated hydraulic gradient of PLTS.
Flow rateQ/(m3⋅h−1)
Flow velocityv/(m⋅s−1)
Shear rater/s-1
Dimensionless parameters
Balance timetb/s
Maximum ofhydraulic gradient imax/(MPa⋅km−1)
Balance hydraulic gradientib/(MPa⋅km−1)
a
b
400
1.57
41.9
81.15
1.27
286
3.97
1.82
450
1.77
47.2
84.09
1.31
270
4.20
1.88
500
1.97
52.4
86.61
1.35
253
4.44
1.95
550
2.16
57.7
88.78
1.39
241
4.67
2.02It is evident from Table5 that the PLTS shows a time-varying behavior during the long-distance pipeline transportation, with the hydraulic gradient declining gradually and finally reaching equilibrium. When the flow rate was 500 m3·h−1, the calculated hydraulic gradient of the PLTS declined by about 56%, from 4.44 MPa·km−1 to 1.95 MPa·km−1 within 253 s. With an increase in the flow rate, the balance time becomes shorter while the balance hydraulic gradient becomes larger. With an average flow velocity of 1.97 m·s−1, it takes approximately four hours to transport the PLTS to the storage. Therefore, it is more systematic and reasonable to use the balance hydraulic gradient for the pipeline design and pump selection. Besides, it is advisable to stir the PLTS thoroughly before transportation, to destroy the FS completely, and reduce the hydraulic gradient effectively, based on the time-varying behavior of the PLTS.


## 4.2. Reliability Analysis

The conventional semiempirical hydraulic gradient formulas, which were derived based on a two-phase flow laboratory date, show a limited applicable scope [18]. For example, Durand’s formula may be suitable for the short-distance transportation of coarse and low-concentration crude tailings. As a widely accepted formula in China, the Shanxi Shuiliyuan formula is suitable for coarse particles, and the Jinchuan formula gives a high precision in calculations related to high-concentration rod-milling sand [19].To evaluate the reliability of the calculated hydraulic gradient by using (8), a semi-industrial loop test was conducted. The hydraulic gradients of the PLTS were calculated by recording the press values of the observation points (Figure 6). A comparison of the hydraulic gradients, calculated using semiempirical formulas, loop test, and (8), is shown in Figure 7. Comparing with the balance hydraulic gradient obtained in the semi-industrial loop test, the computational errors of those calculated by using the time-varying hydraulic gradient model, Jinchuan formula, and Shanxi formula are 15%, 78%, and 130%, respectively. Therefore, the time-varying hydraulic gradient model is more feasible and reliable than semiempirical formulas.Figure 6
Layout of the loop system.Figure 7
Comparison of calculated hydraulic gradient using different methods.


## 5. Conclusions

The conclusions are made as follows:(i) FS of PLTS can not only be stretched or damaged by mechanical mixing, but also lapped and recombined through flocculant bridging. This dynamic damage and recovery process of the FS is the core of the flocculent structure theory. Laboratory shear tests have validated that PLTS is consistent with a shear-thinning pseudo-plastic fluid rather than a time-invariant Bingham plastic fluid. Therefore, the time-varying hydraulic gradient model of PLTS based on the theory of flocculent structure and the Herschel–Bulkley model is more scientific and reasonable.(ii) The shear-thinning characteristics of PLTS are determined by the FS growth level. The shear stress, apparent viscosity, and hydraulic gradient had the highest values when the FS were fully developed and then declined gradually and reached a dynamic balance with the damaging process of FS. Accordingly, additional mechanical agitation of the PLTS before transportation can completely destroy the FS and reduce the hydraulic gradient effectively.(iii) In the case of paste-like tailings in long-distance pipeline transportation, the time-varying hydraulic gradient model, which produces smaller errors than frequently used semiempirical formulas, has proved to be a feasible and high-precision solution for hydraulic gradient calculation.



---

*Source: 5276431-2017-08-23.xml*