Research Article  Open Access
TimeVarying Hydraulic Gradient Model of PasteLike Tailings in LongDistance Pipeline Transportation
Abstract
Pastelike 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 pastelike tailings in longdistance pipeline transportation, to explore a highprecision and reliable hydraulic gradient formula, the rheological behavior of pastelike tailings slurry was analyzed, a timevarying hydraulic gradient model was constructed, and a series of laboratory shear tests were conducted. The results indicate that the PTLS shows noticeable shearthinning 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 semiindustrial loop test, the computational errors of those calculated by using the timevarying hydraulic gradient model, Jinchuan formula, and Shanxi formula are 15%, 78%, and 130%, respectively. Therefore, our model is a feasible and highprecision solution for the calculation of the hydraulic gradient of pastelike tailings in longdistance pipeline transportation.
1. Introduction
In the last twenty years, paste technology has been widely accepted as an environmentfriendly, 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].
Table 1 summarizes several typical examples of paste tailings disposal mines in the world [4]. Evidently, many mines’ tailings are actually disposed as “pastelike” 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 slowsettling 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].

In previous studies, pastelike tailings slurry (PLTS) was simplified as a Bingham plastic fluid, whose flow properties were timeindependent, and the existing empirical formulas for calculating the hydraulic gradient of a pastelike flow were derived based on the timeinvariant Bingham model [7]. However, an increasing number of studies have shown that the rheological characteristics of PLTS are consistent with those of a pseudoplastic 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 timevarying behavior under a continuous shear stress. Senapati and Mishra [10] evaluated the head loss of pastelike coal combustion slurry using the powerlawfluidbased 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 pastelike 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 timevarying hydraulic gradient model for PLTS is of great significance.
This paper discusses the disposal of pastelike tailings in the largest underground iron mine in Asia and presents the timevarying behavior of PLTS by employing groups of shear tests. A timevarying hydraulic gradient model of PLTS in longdistance pipeline transportation was derived based on the theory of flocculent structure and the Herschel–Bulkley model. The reliability of the timevarying hydraulic gradient model was analyzed and compared with a semiindustrial 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 pastelike 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. Figure 2 presents a flow chart of the preparation and transportation of PLTS in Sijiaying.
As part of a widely accepted practice, efficient tailingsthickening equipment, such as deepcone 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 pastelike underflow is transported by using piston pumps in a longdistance pipeline and then dumped and evaporated. There are several benefits cited by Sijiaying in their decision to develop the pastelike 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].
This dynamic damage and recovery process of FS can be described aswhere is the FS coefficient, reflecting the degree of integrity of FS; is shear rate expressed in s^{−1}; is the recovery coefficient, which depends on the flocculant dosage; and is the damage coefficient, directly proportional to . is defined as the stage where the FS are completely destroyed, while 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 can be obtained as
After integrating (1), the following equation can be obtained:
2.3. Rheological Model of PLTS
NonNewtonian fluids include Bingham plastic fluid, pseudoplastic 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 pseudoplastic fluid, and the Herschel–Bulkley model is more applicable to PLTS. The constitutive equation of the Herschel–Bulkley model is commonly written aswhere is the shear stress in Pa; is the yield stress in Pa; is the apparent viscosity in Pa·s; and is the flow index, .
As shown in (3), the FS show noticeable shearthinning characteristics under a particular shear rate, as decreases gradually with the increase in shear time and finally reaches the balanced value . The timevarying behavior of PLTS can be further manifested in the change laws between the parameters of the Herschel–Bulkley model and , which is reflected in the mathematical model given bywhere and are the initial yield stress and initial flow index, respectively; is the apparent viscosity when the FS are completely destroyed; and 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):
2.4. TimeVarying Hydraulic Gradient Model of PLTS
The hydraulic gradient, or resistance loss, is the most important parameter in longdistance 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 timeinvariant Bingham model:where is the hydraulic gradient in Pa/m; is the pipe diameter in m; and is the average flow velocity in m/s.
There is usually a long pipefeeding 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 timevarying behavior, resulting in a large error in the hydraulic gradient calculation. Therefore, the hydraulic gradient formula of PLTS should be modified by using the timevarying flocculent structure theory and the Herschel–Bulkley model. The timevarying hydraulic gradient model of PLTS can then be obtained by integrating (3), (5), and (7):where and are dimensionless parameters in the simplified formula.
From (8), it can be seen that the modified , , and 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 Table 2, 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 slowsettling particles causes a large flocculant dosage and a pastelike underflow from the deepcone thickener. The particle size of the slowsettling particles is denoted by (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.


3.2. Experimental Procedure
The HAAKE VT550 rotary viscometer was used for laboratory shear tests. The immersion sensor, whose crossshaped 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 Figure 4. The shear rate is set as 20 s^{−1}, 30 s^{−1}, 40 s^{−1}, 50 s^{−1}, and 60 s^{−1}.
(a)
(b)
The PLTS of Sijiaying clearly showed shearthinning 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 Table 4. 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 leastsquares method:

To 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.
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 longdistance 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 rate was set to 400, 450, 500, and 550 m^{3}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.

It is evident from Table 5 that the PLTS shows a timevarying behavior during the longdistance pipeline transportation, with the hydraulic gradient declining gradually and finally reaching equilibrium. When the flow rate was 500 m^{3}h^{−1}, the calculated hydraulic gradient of the PLTS declined by about 56%, from 4.44 MPakm^{−1} to 1.95 MPakm^{−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 ms^{−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 timevarying behavior of the PLTS.
4.2. Reliability Analysis
The conventional semiempirical hydraulic gradient formulas, which were derived based on a twophase flow laboratory date, show a limited applicable scope [18]. For example, Durand’s formula may be suitable for the shortdistance transportation of coarse and lowconcentration 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 highconcentration rodmilling sand [19].
To evaluate the reliability of the calculated hydraulic gradient by using (8), a semiindustrial 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 semiindustrial loop test, the computational errors of those calculated by using the timevarying hydraulic gradient model, Jinchuan formula, and Shanxi formula are 15%, 78%, and 130%, respectively. Therefore, the timevarying hydraulic gradient model is more feasible and reliable than semiempirical formulas.
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 shearthinning pseudoplastic fluid rather than a timeinvariant Bingham plastic fluid. Therefore, the timevarying hydraulic gradient model of PLTS based on the theory of flocculent structure and the Herschel–Bulkley model is more scientific and reasonable.
(ii) The shearthinning 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 pastelike tailings in longdistance pipeline transportation, the timevarying hydraulic gradient model, which produces smaller errors than frequently used semiempirical formulas, has proved to be a feasible and highprecision solution for hydraulic gradient calculation.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This study was supported by a grant from the National Key R&D Program of China (2017YFC0804605).
References
 E. Yilmaz, A. Kesimal, B. Ercikdi, and R. Kaya, “The use of paste for disposal of mill tailings,” in Proceedings of the 3rd International Scientific Conference on Modern Management of Mine Producing, Geology and Environmental Protection, pp. 109–120, Varna, Bulgaria, 2003. View at: Google Scholar
 R.B.M. Verburg, “Environmental benefits associated with the use of paste for surface disposal of tailings,” in Proceedings of the 50th Canadian Geotechnical Conference of the Canadian Geotechnical Society, pp. 484–491, Ottawa, Canada, 1997. View at: Google Scholar
 S. Li, Q. Chen, and X. Wang, “Superiority of filtered tailings storage facility to conventional tailings impoundment in southern rainy Regions of China,” Sustainability (Switzerland), vol. 8, no. 11, article no. 1130, 2016. View at: Publisher Site  Google Scholar
 W. Yan, “Tailings disposal methods and applications at home and abroad,” Metal Mine, vol. 483, no. 9, pp. 1–14, 2016 (Chinese). View at: Google Scholar
 S. LI, X.M. WANG, and Q.L. ZHANG, “Dynamic experiments on flocculation and sedimentation of argillized ultrafine tailings using flyashbased magnetic coagulant,” Transactions of Nonferrous Metals Society of China (English Edition), vol. 26, no. 7, pp. 1975–1984, 2016. View at: Publisher Site  Google Scholar
 S. Li and X.M. Wang, “Flyashbased magnetic coagulant for rapid sedimentation of electronegative slimes and ultrafine tailings,” Powder Technology, vol. 303, pp. 20–26, 2016. View at: Publisher Site  Google Scholar
 X.M. Wang, J.W. Zhao, J.H. Xue, and G.F. Yu, “Features of pipe transportation of pastelike backfilling in deep mine,” Journal of Central South University of Technology (English Edition), vol. 18, no. 5, pp. 1413–1417, 2011. View at: Publisher Site  Google Scholar
 A. Haimoni and D. J. Hannant, “Hydraulic transport of solids in pipes: a simple method for the prediction of pressure drop,” Materials and Structures, vol. 26, no. 3, pp. 144–151, 1993. View at: Publisher Site  Google Scholar
 E.U. Pornillos, “A technique for measuring the reduction of yield stress of thickened tailings,” in Proceedings of the 14th International Seminar on Paste and Thickened Tailings, pp. 167–173, Australian Centre of Geomechanics, Perth, Australia, 2011. View at: Google Scholar
 P. K. Senapati and B. K. Mishra, “Design considerations for hydraulic backfilling with coal combustion products (CCPs) at high solids concentrations,” Powder Technology, vol. 229, pp. 119–125, 2012. View at: Publisher Site  Google Scholar
 S. Li, X.M. Wang, Q.L. Zhang, and S. Wang, “Timevarying behavior of pastelike superfine tailings in long selfflowing transportation,” Journal of Northeastern University (Natural Science), vol. 37, no. 7, pp. 1045–1049, 2016 (Chinese). View at: Google Scholar
 R.J. Jewell and A.B. Fourie, Paste and Thickened Tailings: A Guide, Australian Centre for Geomechanics, Perth, Australia, 2nd edition, 2006.
 F. Palkovits, “Paste thickening: considerations for backfill vs. tailings management,” Engineering and Mining Journal, vol. 212, no. 9, pp. 34–40, 2011. View at: Google Scholar
 H. E. Ries Jr. and B. L. Meyers, “Flocculation mechanism: Charge neutralization and bridging,” Science, vol. 160, no. 3835, pp. 14491450, 1968. View at: Publisher Site  Google Scholar
 E. Antunes, F. A. P. Garcia, P. Ferreira, A. Blanco, C. Negro, and M. G. Rasteiro, “Modelling PCC flocculation by bridging mechanism using population balances: Effect of polymer characteristics on flocculation,” Chemical Engineering Science, vol. 65, no. 12, pp. 3798–3807, 2010. View at: Publisher Site  Google Scholar
 E. A. Toorman, “Modelling the thixotropic behaviour of dense cohesive sediment suspensions,” Rheologica Acta, vol. 36, no. 1, pp. 56–65, 1997. View at: Publisher Site  Google Scholar
 R.P. Chhabra and J.F. Richardson, NonNewtonian Flow and Applied Rheology: Engineering Applications, ButterworthHeinemann, Oxford, UK, 2008.
 Q.L. Zhang, G.Y. Hu, and X.M. Wang, “Hydraulic calculation of gravity transportation pipeline system for backfill slurry,” Journal of Central South University of Technology, vol. 15, no. 5, pp. 645–649, 2008. View at: Publisher Site  Google Scholar
 I.T. Seitshiro, J. Katende, and H. Sato, “Slurry pipeline design of multisized solids: application of innovated models,” International Journal of the Society of Materials Engineering for Resources, vol. 20, no. 2, pp. 201–206, 2014. View at: Google Scholar
Copyright
Copyright © 2017 Li Yang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.