Research Article | Open Access

Volume 2014 |Article ID 876872 | https://doi.org/10.1155/2014/876872

Chuan-sheng Wu, "The Drainage Consolidation Modeling of Sand Drain in Red Mud Tailing and Analysis on the Change Law of the Pore Water Pressure", Mathematical Problems in Engineering, vol. 2014, Article ID 876872, 10 pages, 2014. https://doi.org/10.1155/2014/876872

# The Drainage Consolidation Modeling of Sand Drain in Red Mud Tailing and Analysis on the Change Law of the Pore Water Pressure

Accepted11 Jun 2014
Published22 Jul 2014

#### Abstract

In order to prevent the occurring of dam failure and leakage, sand-well drainages systems were designed and constructed in red mud tailing. It is critical to focus on the change law of the pore water pressure. The calculation model of single well drainage pore water pressure was established. The pore water pressure differential equation was deduced and the analytical solution of differential equation using Bessel function and Laplace transform was given out. The impact of parameters such as diameter , separation distance , loading rate , and coefficient of consolidation in the function on the pore water pressure is analyzed by control variable method. This research is significant and has great reference for preventing red mud tailings leakage and the follow-up studies on the tailings stability.

#### 1. Introduction

Red mud, a by-product of alumina refining, is produced in increasing quantities globally [1, 2]. Much of this red mud has traditionally been produced by methods which create a liquor of high moisture content [3, 4], which are deposited in impoundments that can reach depths of 20 m . Red mud in this tailing is infiltrated by water and in a status of saturation, resulting in a low strength and high moisture content. As a consequence, the continuous stacking of Sintering red mud into this tailing might cause the dam failure and the leak-prone of the tailing when the original designed seepage prevention rank is low, especially for the tailing built in Karst area [8, 9]. As the mineral and chemical composition of red mud is very complex, red mud and its attached fluid have strong alkalinity and corrosivity [10, 11]. It is dangerous that the failure and leakage of this tailings dam occurs and further cause serious consequences that widespread environmental pollution [12, 13].

To prevent the occurring of dam failure and leakage, sand well drainage systems were designed and constructed in this red mud tailing, the lower part of which is mainly soft Bayer red mud . With the assistance of sand well drainage systems, stacking of sintered dry red mud onto the tailings will promote the drainage of pore water and reinforcement of the basement of the tailings. Simultaneously, the pressure of pore water will change continuously as well . Consequently, it is vital to focus on the law of pore water pressure distribution and change in the lower soft red mud during the stacking of Sintering red mud. There was a great deal of reports on the application of the sand well drainage method and one-dimensional drainage-consolidation calculation for saturated soils since 1930s . However, as for the calculation on the drainage-consolidation by sand wells, the most well-known model is undoubtedly the axisymmetric radial drainage-consolidation equations of a single well, which was deduced and proposed by Barron in 1948 . This equation not only gives us the analytical solution under both the conditions of isostrain and free strain, but also has considered the fluid seepage flowing-strain interaction.

In this theoretical research field of the in soft soil foundation, although a solution to this compound problem has been given with the help of the three-dimensional parabolic equation under the 2nd type boundary condition, yet, regretfully, it has no sufficient accuracy . Therefore, it is important to give an accurate solution to this equation. This research aims at giving the analytical solution of hydrostatic pressure of the three-dimensional consolation utilizing the mathematical calculation methods such as Laplace transform. Crucially, this solution has a high applicability. This research is significantly important to the calculation of pore water pressure and the further control of leakage. In addition, it can theoretically support the follow-up studies on the tailings stability.

#### 2. Physical Model of the Drainage Consolidation System

Because of the large area and complex shape of the red mud tailings researched in this project, the theoretical calculation is difficult to be carried out without reasonable simplification and the further construction of calculation model. The sand drains are designed and sited uniformly and equally, which means that these sand drains are alike in function and effect. Then, the drainage consolidation of the whole tailing red mud can be obtained through research on each sand drain and their superposition. The schematic diagram of two well-accepted distribution modes of sand drains are shown as in Figure 1. In addition, the smear zone formed as a result of construction disturbance when the sand drain separation distance is lower than 1.0 m . It can be overlooked when the sand drain separation distance is larger than that of building and road foundation engineering . In addition, seepage resistance caused by the material type of the sand drains, for the piling rate is low and the drainage time is long. As for the piling of the red mud, it can be assumed to be continuous and the production of aluminum as well.

As referred before, the affection region of each sand drain could be equivalent to a cylinder (Figure 2). The equivalent diameter can be deduced from the equation of , where refers to the sand drain separation distance . The value of would be 1.13 when it comes to square-like sand drain distribution and 1.05 for triangle-like one. We will also give the assumption that the seepage of pore water contained in the affection cylinder region would be along the radial direction towards the sand drain and then along the vertical direction towards the underground.

#### 3. Derivation of Pore Pressure Function of the Red Mud Foundation with Sand Drain

In Barron’s theoretical research, although a solution to this compound problem has been given with the help of the three-dimensional parabolic equation under the 2nd type boundary condition, this solution is not accurate. We have tried to calculate the analytical solution of hydrostatic pressure of the three-dimensional consolation through the mathematical calculation methods such as Laplace transform.

Following Barron theory, we do the following basic assumptions for the consideration of pore water pressure under variable load change:(1)the soil skeleton is isotropic material, regardless of the efficiency of the soil mass deformation and creep;(2)soil particle and pore water can not be compressed;(3)seepage flow both in vertical and horizontal radial directions obeys Darcy’s law and have the same permeability coefficient;(4)soil particles move in the vertical direction; the applied load is a first-order function of time.

Take a cell cube from the computational domain as shown in Figure 3 and set that a unit of the direction is and a unit of the direction is . So, the amount of soil volume change during the process of drainage consolidation is equal to the amount of pore water. Consider

Consider the seepage flow balance both in horizontal and vertical directions as follows: where is the void ratio of time , and void ratio of time is . Take this into (2) and then get

And by the compression test curve, we have where is a compression coefficient. Putting (5) into (4) can deduce where is osmotic coefficient, and which can be put into (6) to get where represents the vertical coefficient of consolidation. It can be calculated from by the definition formula , where is the density of water. That means can be gotten from the data of both consolidation test and penetration test.

Change the former equation to accord with the specific boundary of sand drain:

According to the single well drainage consolidation model shown in Figure 2, the upper boundary has curtained pore pressure, the left one is water-proof, the right is sand drain, and the physical radius and action radius of which are separately and . Then, the drainage consolidation with the addiction of boundary conditions can be modeled as the following mathematical physical equation:

In this equation, is excess hydrostatic pressure; is pole diameter of the cylindrical polar coordinates, flow coordinates; is consolidation time, flowing coordinates; is coordinates for vertical flow, calculation depth; is additional load, time function of ; is the consolidation coefficient, determined by experiment; is soil drainage distance, soil thickness, which is constant.

According to the theory of soil mechanics, the three-direction water flow in the soil is believed to be turbulent flow problem and can be converted into two parts, planar radial flow (or radiation) and vertical (or linear) laminar flow problem. Based on mathematical method, the total three direction water flow problem also can be transformed to subproblem I and subproblem II.Subproblem I is Subproblem II is

Assuming the solution of subproblem I is and solution of subproblem II is , then the solution to the total question must be confirmed to the form of or

##### 3.1. The Theoretical Solution of Subproblem II

In (12), there is where is Poisson ratio of soil; is the action radius of sand drain, half of the sand drain separation distance; is the physical radius of sand drain; is coefficient of consolidation in horizontal direction.

Equation (15) is the resolution of stress with the assumption that the soil is linear elasticity. Put (15) into (12) and then it can be deduced that Equation (16) is the equilibrium equation in horizontal direction which can be rewritten as where , is engineered loading velocity. Both and are constant. Equation (17) can be further deduced as follows by Laplace transformation: Equation (18) can be solved as

Bessel solution of (19) can be calculated as Equation (20) can be transformed into Assume , and then it can be deduced to The and can be calculated by the boundary conditions in (12), where and (24) can be deduced as Equation (25) can be simplified as In this equation, , (26) can be further deduced as following by anti-Laplace transformation: Equation (27) can be written as The singularity of (28) is , and all the null points of , and the residue of singularity of is As for the root of , can be replaced by and can certified by . Then, we have In the equation, we have . The root of formulas , can be natural frequency and all are of arithmetic number. Residue, when , is where is . According to the character of Bessel function, it can be deduced that Then, the excess hydrostatic pressure can be

##### 3.2. The Theoretical Solution of Subproblem I

The kinematic equation of subproblem I is The boundary conditions are Equation (34) can be further deduced as follows by Laplace transformation: Solving the upper equations, Putting the boundary conditions can give the value that Equation (37) can be deduced as As for the , null point of is Besides, is the moving singularity of ; other points are all first grade singularity of . The it can be deduced by the expansion of residue theorem:Then,

Then, we can put (33) and (42) into (14) of total question to get the final solution. Consider

This is the final expression of the excess hydrostatic pressure changing along the value of pole diameter of the cylindrical polar coordinates , consolidation time , and calculation depth .

#### 4. Analysis on the Change Law of the Pore Water Pressure

From the derivation of former part of the research, the function of the excess hydrostatic pressure is influenced by multiple variables, whose variation characteristic cannot be vividly expressed in a three-dimensional space coordinate system. For the purpose of understanding the variation rule of with the time and space, we use MATLAB as a tool, by examining a parameter and fixing other parameters, to get the affects way of each variable effects on the value of . Based on the previous engineering experience and based on the survey and design of tailings and actual production, scope of calculation parameter is as shown in Table 1.

 Parameter and its extremum Diameter ofsand drain (m) Separation distanceof sand drain (m) Loading rate (kPa/y) Max 0.35 5.0 45 Min 0.10 1.0 25

Control variate method is employed to conclude the influence of various parameters on the law of drainage. Concretely, what this means is to study the effect of target variable from the four variables, diameter , separation distance , loading rate , and coefficient of consolidation , and fix the other three ones as initial value.

On the basis of scope of each calculation parameter (Table 1), the initial values are set as Table 2.

 Diameter of sand drain Separation distance of sand drain Calculation depth Loading rate Coefficient of consolidation Action radius 0.30 m 3.0 m 20.0 m 39 kPa/y 0.002 1.575 m

#### 5. Conclusion

This research has established the calculation model of single well drainage pore water pressure, deduced pore water pressure differential equation based on the three-dimensional consolidation, and get the analytical solution of differential equation using Bessel function and Laplace transform. The impact of parameters in the function on the pore water pressure is analyzed by control variable method. Changes of excess pore water pressure on four variables such as diameter , separation distance , loading rate , and coefficient of consolidation are analyzed and proved to be nonlinear. The calculation results help to further set up a depository stability evaluation model to prevent loading exceeding design capacity. And together with inhomogeneous deformation data of the field, it is promising to create warning systems for site managers based merely on the field logging data of pore pressure data. Therefore, this research has a guiding significance on preventing red mud tailings leakage of continuous loading pile exceeding design capacity.

#### Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

1. M. Gräfe and C. Klauber, “Bauxite residue issues: IV. Old obstacles and new pathways for in situ residue bioremediation,” Hydrometallurgy, vol. 108, no. 1-2, pp. 46–59, 2011. View at: Publisher Site | Google Scholar
2. U.S. Geological Survey, “Bauxite and alumina statistics and information,” in Historical Statistics for Mineral and Material Commodities in the United States, D. T. Kelly and G. R. Matos, Eds., U.S. Geological Survey, 2008. View at: Google Scholar
3. M. Gräfe, G. Power, and C. Klauber, “Bauxite residue issues: III. Alkalinity and associated chemistry,” Hydrometallurgy, vol. 108, no. 1-2, pp. 60–79, 2011. View at: Publisher Site | Google Scholar
4. J. Wehr, I. Fulton, and N. Menzies, “Revegetation strategies for bauxite refinery residue: a case study of Alcan Gove in Northern Territory, Australia,” Environmental Management, vol. 37, no. 3, pp. 297–306, 2006. View at: Publisher Site | Google Scholar
5. C. Klauber, M. Gräfe, and G. Power, “Bauxite residue issues: II. options for residue utilization,” Hydrometallurgy, vol. 108, no. 1-2, pp. 11–32, 2011. View at: Publisher Site | Google Scholar
6. R. Bott, T. Langeloh, and J. Hahn, “Re-usage of dry bauxite residue,” in Proceedings of the 7th International Alumina Quality Workshop, A. McKinnon, Ed., pp. 236–241, AQW Inc., Perth, Australia, 2005. View at: Google Scholar
7. C. Brunori, C. Cremisini, P. Massanisso, V. Pinto, and L. Torricelli, “Reuse of a treated red mud bauxite waste: Studies on environmental compatibility,” Journal of Hazardous Materials, vol. 117, no. 1, pp. 55–63, 2005. View at: Publisher Site | Google Scholar
8. W. Liu, J. Yang, and B. Xiao, “Review on treatment and utilization of bauxite residues in China,” International Journal of Mineral Processing, vol. 93, no. 3-4, pp. 220–231, 2009. View at: Publisher Site | Google Scholar
9. P. Renforth, W. M. Mayes, A. P. Jarvis, I. T. Burke, D. A. C. Manning, and K. Gruiz, “Contaminant mobility and carbon sequestration downstream of the Ajka (Hungary) red mud spill: the effects of gypsum dosing,” Science of the Total Environment, vol. 421-422, pp. 253–259, 2012. View at: Publisher Site | Google Scholar
10. P. Wang and D. Y. Liu, “Physical and chemical properties of sintering red mud and bayer red mud and the implications for beneficial utilization,” Materials, vol. 5, no. 10, pp. 1800–1810, 2012. View at: Publisher Site | Google Scholar
11. C. Wu and D. Liu, “Mineral phase and physical properties of red mud calcined at different temperatures,” Journal of Nanomaterials, vol. 2012, Article ID 628592, 6 pages, 2012. View at: Publisher Site | Google Scholar
12. A. Gelencsér, N. Kováts, B. Turóczi et al., “The red mud accident in Ajka (Hungary): characterization and potential health effects of fugitive dust,” Environmental Science and Technology, vol. 45, no. 4, pp. 1608–1615, 2011. View at: Publisher Site | Google Scholar
13. W. M. Mayes, A. P. Jarvis, I. T. Burke et al., “Dispersal and attenuation of trace contaminants downstream of the ajka bauxite residue (red mud) depository failure, Hungary,” Environmental Science and Technology, vol. 45, no. 12, pp. 5147–5155, 2011. View at: Publisher Site | Google Scholar
14. J. Paul Rinehimer, J. Thomson, and C. Chris Chickadel, “Thermal observations of drainage from a mud flat,” Continental Shelf Research, vol. 60, pp. S1250–S135, 2013. View at: Publisher Site | Google Scholar
15. Q. Bi, J. Jia, and Z. Huang, “Numerical analysis for consolidation drainage scheme of red mud disposal,” Journal of North China Institute of Water Conservancy and Hydroelectric Power, vol. 33, pp. 87–90, 2012. View at: Google Scholar
16. K. Terzaghi, Theoretical Soil Mechanics, John Wiley & Sons, New York, NY, USA, 1943.
17. A. Onoue, “Consolidation by vertical drains taking well resistance and smear into consideration,” Soils and Foundations, vol. 28, no. 4, pp. 165–174, 1988. View at: Publisher Site | Google Scholar
18. P. Orleach, Technique to evaluate the field performance of vertical drains [M.S. thesis], Massachusetts Institute of Technology, Cambridge, Mass, USA, 1983.
19. X. W. Tang and K. Onitsuka, “Consolidation of ground with partially penetrated vertical drains,” Geotechnical Engineering, vol. 29, no. 2, pp. 209–231, 1998. View at: Google Scholar
20. R. A. Barron, “Consolidation of fine-grained soils by drain wells,” Transactions of the American Society of Civil Engineers, vol. 113, no. 2346, pp. 718–754, 1948. View at: Google Scholar
21. H. Yoshikuni and H. Nakanado, “Consolidation of fine-grained soils by drain well with filter permeabiiity,” Soil and Foundation, vol. 14, pp. 35–46, 1974. View at: Publisher Site | Google Scholar
22. G. X. Zeng and K. H. Xie, “New development of the vertical drain theories,” in Proceedings of the 12th International Conference on Soil Mechanics and Foundation Engineering, vol. 2, pp. 1435–1438, Rio de Janeiro, Brazil, August 1989. View at: Google Scholar
23. G. Zhu and J. Yin, “Consolidation analysis of soil with vertical and horizontal drainage under ramp loading considering smear effects,” Geotextiles and Geomembranes, vol. 22, no. 1-2, pp. 63–74, 2004. View at: Publisher Site | Google Scholar
24. W. Li, “Numerical simulation of sand drains foundation with a consideration of permeability coefficient changing along the radial direction,” Subgrade Engineering, vol. 6, pp. 70–73, 2012. View at: Google Scholar
25. L. Ding, D. Ling, L. Wang, and Y. Chen, “Application of composite element method in sand drain ground finite element analysis,” Chinese Journal of Computational Mechanics, vol. 21, no. 1, pp. 13–20, 2004. View at: Google Scholar
26. K. Zhao and Y. Zheng, “Shear-resistance strength of red mud after consolidation and drainage,” Guizhou Geology, vol. 13, pp. 280–286, 1996. View at: Google Scholar