This study determined the magnitude of the hydraulic anisotropy ratio (rk = kh/kv) in undisturbed tailings samples with varying fines contents. Constant-head permeability tests were performed under different consolidation stress levels and hydraulic gradients, evaluating the influence of effective stress, void ratio, and fines content on the compression index (Cc), the coefficient of consolidation (Cv), and permeability (k). To obtain rk, the experimental method of cutting samples in different directions was implemented. This method has been applied to obtain rk in various materials such as clays and rocks, but it is not known to have been used for tailings. The literature reports few rk values for tailings. According to different researchers, rk in tailings can vary from 2 for fairly uniform dams to 100 in transition zones, due to the intercalation of fine and coarse particles. In this study, rk values ranging from 1 to 10 were obtained, identifying fines content as the main factor controlling hydraulic anisotropy. Furthermore, three correlations were defined between fines content and initial and final void ratios.

Cc

compression index

Cv

coefficient of consolidation

e

void ratio

ef

final void ratio

ei

initial void ratio

emax

maximum void ratio

emin

minimum void ratio

en

natural void ratio

Fc

fines content

Fchs

fines content in the horizontal sample

Fcvs

fines content in the vertical sample

Gs

specific gravity

h

hydraulic head

i

hydraulic gradient

k

permeability coefficient in the direction perpendicular to the stratification plane

k

permeability coefficient parallel to the stratification plane

kh

permeability in the horizontal direction

kv

permeability in the vertical direction

rk

hydraulic anisotropy ratio

U

discharge velocity

β

Skempton’s parameter

ΔFc

difference in the fines content

σ

effective stress

ω

natural water content

The hydraulic conductivity or permeability of a porous medium is a measure of the ease with which a fluid flows through the medium; the magnitude of the permeability is associated with the structure of the voids in the medium (Gu et al., 2022). Permeability is one of the properties that present great natural variability in soils, so it is common to find differences of up to several orders of magnitude among the permeabilities of the materials of the same geological formation (López-Acosta, 2017).

In general, soils exhibit variations in their physical and mechanical properties that can be classified in two ways: heterogeneity, which is the variation of a parameter over space within a geological formation; and anisotropy, which is the vector variation of a property at a given point within a geological formation (Clavaud et al., 2008). Soil heterogeneity results from the natural geological processes by which soil is deposited and formed as a geomechanical material (Jamshidi Chenari and Behfar, 2017). In contrast, anisotropy is a natural outcome of the soil’s depositional history and stress history, and can be attributed to both inherent anisotropy and induced anisotropy. Inherent anisotropy refers to the preferential orientation of the soil fabric established during particle deposition and sedimentation, whereas induced anisotropy refers to changes in particle arrangement and orientation caused by the application of external loads (Jamshidi Chenari and Behfar, 2017; Jamshidi Chenari and Mahigir, 2014).

According to Witt and Brauns (1983), there are three main causes of hydraulic anisotropy: macro-stratification, micro-stratification, and the orientation of non-spherical particles. Macro-stratification is caused by changes in materials and depositional conditions over geological time, leading to successive layers with different grain-size distributions, textures, or lithologies. Layers with higher permeability will convey most of the flow under a flow regime parallel to the stratification, whereas layers with lower permeability will govern the flow regime perpendicular to the stratification, thus creating anisotropy (Nordquist, 2015). The microstructure comprises the soil skeleton, interparticle connections, particle contacts, and pore morphology (Xu et al., 2020). Hydraulic anisotropy is therefore largely controlled by the soil’s macro and microstructure. Particle orientation is typically the result of compression, where particles align perpendicular to the principal stress. Several researchers have recognised the significant influence of particle orientation on the magnitude of directional permeability. At the particle scale, hydraulic anisotropy develops due to the tortuosity of the flow path within the soil (Scholes et al., 2007).

The hydraulic anisotropy ratio (rk) is commonly defined as the ratio of the permeability coefficient parallel to the stratification plane k (horizontal flow) and the permeability coefficient in the direction perpendicular to the stratification plane k (vertical flow); therefore (Bolton et al., 2000; Chan and Kenney, 1973; Scholes et al., 2007):

1

Since the horizontal permeability is typically greater than the vertical one, the hydraulic anisotropy ratio is usually of a magnitude larger than unity; however, vertical permeability values greater than the horizontal permeability (anisotropy smaller than unity) have been found at shallow depths. Micro-cracks induced by stresses are likely to generate preferential values of vertical permeability (Bolton et al., 2000; Chapuis and Gill, 1989).

Tailings are fine solid residues generated from the beneficiation of minerals by grinding processes in a mineral concentration plant (Orozco, 2010). Their origin begins with the extraction of the mineral, followed by the stages of crushing, grinding, concentration (flotation, leaching, or other methods), and ends with the deposition of the waste in a semi-liquid state, which is transported through pipelines to the deposit (Blight, 2010; Candelaria et al., 2022; ICOLD, 2001; Vick, 1990). During discharge, particles settle according to their size and specific gravity: larger and heavier particles settle more rapidly and closer to the discharge point, while smaller and lighter particles are carried to more distal zones. This depositional mechanism gives tailings deposits a stratified, heterogeneous and anisotropic structure, with layers of different permeabilities and compressibilities that directly influence their hydraulic and mechanical behaviour (Kossoff et al., 2014; Villavicencio et al., 2014).

Permeability is one of the most important properties in the design of tailings dams, as it governs both the consolidation of the material and the water flow conditions within and through the deposit (Bussiere, 2004). Tailings generally exhibit very low permeability values and often display significant anisotropic conditions (López-Acosta and Mendoza-Promotor, 2016). Other factors on which hydraulic anisotropy in tailings dams depends are the milling process, the possible interruption of the disposal site operation, and climatic conditions (Perlea and Botea, 1984).

In the literature, many rk values are found for clays, sands, and rocks, but few results are available for tailings. Due to their layered nature, tailings dams exhibit considerable variation in permeability between the parallel (kh) and the perpendicular (kv) directions with respect to stratification (Abadjiev, 1976). In the case of quite uniform sandy beaches that are not sufficiently compacted, and for silty areas deposited underwater, the hydraulic anisotropy ratio is in the range from 2 to 10 (Saad, 2008; Vick, 1990). Transition zones between clean coarse sands and silts show higher anisotropy factors, due to the intermediate layer of finer and coarser particles; for tailings dams where discharge procedures are not well controlled, the anisotropy ratio can be as high as 100 or more (Vick, 1990; Witt and Schönhardt, 2004). The hydraulic anisotropy measured by different researchers in tailings is presented in Table 1.

Table 1.

Hydraulic anisotropy ratio in tailings

Referencerk
Saad (2008) 0.1–0.8
Pettibone and Kealy (1971) 5 to 10
Vermeulen (2001)  
Kealy and Busch (1979)10
Vick (1990) 2–10
Van der Berg (1995)7–22
Wagener et al. (1998)25
McPhail and Wagner (1989)100
Oliva-González et al. (2024) 10
ICOLD (2021)10–20

The hydraulic anisotropy ratio within the tailings dam affects the seepage pattern and influences the position of the water table; due to anisotropy, the phreatic level at the flow discharge point shows a significant increase compared with the isotropic case. Nelson et al. (1977) observed that a small rise in the water table can result in zones of saturation at the face of the reservoir and can affect the overall stability of the reservoir, so that as the accuracy required in water flow models increases, the need to take hydraulic anisotropy into account increases.

Due to the particular characteristics of tailings and the limited information available in the literature regarding the hydraulic anisotropy ratio in this type of material, the main objective of this investigation was to determine in the laboratory the magnitude of the hydraulic anisotropy ratio in undisturbed tailings samples with different fines contents. To this end, constant-head permeability tests were performed, and the results obtained were compared with those reported by other researchers. As part of the experimental programme, the compression index (Cc) and coefficients of consolidation (Cv) were also determined, evaluating the influence of parameters such as fines content, permeability, and void ratio.

The soil studied in this work corresponds to undisturbed tailings samples recovered from three tailings dams located in Mexico.

Table 2 shows a summary of the properties of the tailings studied: natural water content w, [ASTM D2216-19 (ASTM, 2019)]; specific gravity Gs, [ASTM D854-23 (ASTM, 2023)]; fines content, Fc [ASTM D1140-17 (ASTM, 2017b)], and the natural void ratio en.

Table 2.

Properties of tested tailings

Sample#Orientationω: %GsFc: %en
1Vertical43.363.0898.961.34
 Horizontal40.93 99.391.29
2Vertical40.993.2098.901.44
 Horizontal40.96 99.321.35
3Vertical42.773.0699.391.36
 Horizontal46.18 99.601.5
4Vertical12.072.8364.550.65
 Horizontal17.74 37.420.56
5Vertical27.912.8799.000.84
 Horizontal28.02 98.580.85
6Vertical12.912.9036.400.81
 Horizontal14.58 18.320.77
7Vertical17.152.8842.390.57
 Horizontal14.28 32.770.71
8Vertical24.563.2887.820.84
 Horizontal21.03 81.330.69
9Vertical21.502.8881.500.62
 Horizontal18.76 85.730.56
10Vertical25.142.7984.340.76
 Horizontal24.98 79.420.76
11Vertical21.132.8079.320.63
 Horizontal122.89 98.350.68
 Horizontal219.21 84.430.71
12Vertical22.172.7888.000.73
 Horizontal122.83 79.900.75
 Horizontal223.08 85.860.75
 Horizontal326.15 96.000.77
13Vertical17.082.7648.220.69
 Horizontal16.14 43.730.66
14Vertical18.072.7780.090.64
 Horizontal19.49 80.900.71

The hydraulic anisotropy ratio (rk) was assessed in undisturbed tailings samples using a constant-head, flexible-wall permeameter. This device consists of pressure and volume change sensors, as well as a signal conditioning and data acquisition system, as shown in Figures 1(a) and 1(b).

Figure 1.
An image depicts a laboratory setup with various valves and gauges, a computer, and an exploded view of a mechanical assembly labelled with parts including drains and O rings.The image depicts a laboratory setup divided into 3 parts. The first part features a metallic panel equipped with several valves and gauges, clearly labelled with numbers for identification, including 2 vertical cylinders and valves for controlling fluid flow. The second part includes a computer and monitor setup where a software interface related to the experiment is visible. The third section depicts an exploded view of a mechanical assembly showing components such as a Cover, Top drain, O-ring, Bolster, Porous stones, Tailings, Pedestal, Lower drain, and Base. Each component is marked with arrows to indicate its arrangement and relationships clearly.

(a) Constant head permeameter: 1 – lower back pressure regulator, 2 – upper back pressure regulator, 3 – confinement regulator, 4 – bottom flow burette, 5 – top flow burette, 6 – air–water interface, 7 – air–water interface, 8 – lower volumetric change sensor, 9 – upper volumetric change sensor, 10 – test chamber, 11 – lower back pressure sensor, 12 – upper back pressure sensor, 13 – confining pressure sensor; (b) data acquisition system; and (c) mounting of tailings specimens in the chamber

Figure 1.
An image depicts a laboratory setup with various valves and gauges, a computer, and an exploded view of a mechanical assembly labelled with parts including drains and O rings.The image depicts a laboratory setup divided into 3 parts. The first part features a metallic panel equipped with several valves and gauges, clearly labelled with numbers for identification, including 2 vertical cylinders and valves for controlling fluid flow. The second part includes a computer and monitor setup where a software interface related to the experiment is visible. The third section depicts an exploded view of a mechanical assembly showing components such as a Cover, Top drain, O-ring, Bolster, Porous stones, Tailings, Pedestal, Lower drain, and Base. Each component is marked with arrows to indicate its arrangement and relationships clearly.

(a) Constant head permeameter: 1 – lower back pressure regulator, 2 – upper back pressure regulator, 3 – confinement regulator, 4 – bottom flow burette, 5 – top flow burette, 6 – air–water interface, 7 – air–water interface, 8 – lower volumetric change sensor, 9 – upper volumetric change sensor, 10 – test chamber, 11 – lower back pressure sensor, 12 – upper back pressure sensor, 13 – confining pressure sensor; (b) data acquisition system; and (c) mounting of tailings specimens in the chamber

Close modal

Results of the hydraulic anisotropy ratio of soils and rocks using a variety of laboratory techniques have been reported in the literature (Adams et al., 2016; Chan and Kenney, 1973; Lucero, 2022), as well as of field techniques (Chapuis and Gill, 1989).

To evaluate the hydraulic anisotropy ratio in undisturbed tailings samples, the experimental method of trimming specimens in different directions presented by Chapuis and Gill (1989) has been used; this methodology allows the determination of rk values in different materials such as clays (Lucero, 2022) and rocks, but it is not known that it has been applied to find rk in tailings. The procedure consists of first determining the vertical permeability, followed by the horizontal permeability. Once permeability has been measured in both directions, the hydraulic anisotropy ratio (rk) is calculated. It is important to note that the specimens tested for kv and kh were taken from the same Shelby tube but not from the same horizon, since the specimen dimensions do not allow this; however, they were tested under identical stress conditions and hydraulic gradients. Figure 2 shows two specimens of different tailings, trimmed in both directions, to illustrate the highly stratified nature of the material.

Figure 2.
Four laboratory images display two cylindrical samples under pressure testing, one sample is a brown material and the other a grey material.The image consists of 4 laboratory photographs depicting 2 cylindrical samples being tested under pressure. The left side depicts the process of compressing a brown material sample, with 2 views, one depicts the sample within a cylindrical testing apparatus, and the other depicts the sample after removal, showcasing its texture. The right side depicts a grey material sample undergoing a similar compression test, with the first image depicting the apparatus and the second depicting the cylindrical sample post testing. Each set of images illustrates the physical properties of the respective materials being tested, with clear focus on the structure and texture of the samples.

Cutting of undisturbed tailings samples: (a) evaluation of kv and (b) evaluation of kh

Figure 2.
Four laboratory images display two cylindrical samples under pressure testing, one sample is a brown material and the other a grey material.The image consists of 4 laboratory photographs depicting 2 cylindrical samples being tested under pressure. The left side depicts the process of compressing a brown material sample, with 2 views, one depicts the sample within a cylindrical testing apparatus, and the other depicts the sample after removal, showcasing its texture. The right side depicts a grey material sample undergoing a similar compression test, with the first image depicting the apparatus and the second depicting the cylindrical sample post testing. Each set of images illustrates the physical properties of the respective materials being tested, with clear focus on the structure and texture of the samples.

Cutting of undisturbed tailings samples: (a) evaluation of kv and (b) evaluation of kh

Close modal

The experimental programme to obtain the hydraulic anisotropy ratio consisted of 31 tests, of which 14 were tested in the vertical direction and 17 in the horizontal direction.

The experimental procedure consists of six stages. The first stage is the formation and mounting of samples, where the specimen is worked in the desired direction, in the direction of the field extraction for vertical specimens and in the transverse direction of the field extraction for horizontal specimens. The samples are mounted in the constant head permeability equipment, confined by a latex membrane (Figure 1(c)).

Most of the tailing samples tested had saturation levels above 70%, so to facilitate saturation, water was circulated and saturated with back pressure. To indirectly estimate sample saturation, Skempton’s β was measured, such that the sample is considered saturated when this parameter is equal to or greater than 0.95 [ASTM D4767-11 (ASTM, 2020)].

Because the test is non-destructive during the consolidation phase, all samples were tested at six effective stress levels: σ′ = 19.6, 39.2, 78.5, 156.9, 245.2, and 313 kPa. These stresses were assigned to observe the variation in the void ratio, compressibility parameters, and permeability with stress, as well as to compare the samples, ensuring that all were evaluated under the same experimental conditions. One of the limitations of the equipment is that it is not possible to measure the displacements sustained by the specimen when the effective stress is applied, but a record of the output volumes is kept as established by the standard ASTM D5084-24 (ASTM, 2024); therefore, it is assumed that the total volumetric changes affect only the cross-sectional area of the specimen, leaving the height of the specimen constant.

In the testing stage, the vertical and horizontal permeabilities are determined for each of the effective stresses applied. The test is carried out by applying an upward flow through the specimens, measuring the permeability under four hydraulic heads: h = 150, 217, 284, and 350 cm of water column; the flow rate is calculated from the slope of the time against volumetric change line, as shown in Figure 3.

Figure 3.
Two graphs depict volumetric change over time for different heights, each with linear regression equations. The graphs display 4 distinct lines, each representing a different height.The image features two graphs labelled A and B, illustrating the relationship between volumetric change in cubic centimetres and time in seconds for varying heights. The left graph A indicates a volumetric change percentage of 99.4 percent, while the right graph B indicates 99.6 percent. Each graph contains 4 lines, corresponding to heights of 150 centimetres, 217 centimetres, 284 centimetres, and 350 centimetres, represented by different colours. Each line is accompanied by a linear regression equation in the form of y equals m x plus b, where m is the slope and b is the y intercept. The X-axis and Y-axis are clearly marked, and a legend provides clarification on the colour coding for the heights. The data is organised into 2 separate but similar graphs, arranged left to right for comparison.

Applied hydraulic heads (sample 3): (a) evaluation of kv and (b) evaluation of kh

Figure 3.
Two graphs depict volumetric change over time for different heights, each with linear regression equations. The graphs display 4 distinct lines, each representing a different height.The image features two graphs labelled A and B, illustrating the relationship between volumetric change in cubic centimetres and time in seconds for varying heights. The left graph A indicates a volumetric change percentage of 99.4 percent, while the right graph B indicates 99.6 percent. Each graph contains 4 lines, corresponding to heights of 150 centimetres, 217 centimetres, 284 centimetres, and 350 centimetres, represented by different colours. Each line is accompanied by a linear regression equation in the form of y equals m x plus b, where m is the slope and b is the y intercept. The X-axis and Y-axis are clearly marked, and a legend provides clarification on the colour coding for the heights. The data is organised into 2 separate but similar graphs, arranged left to right for comparison.

Applied hydraulic heads (sample 3): (a) evaluation of kv and (b) evaluation of kh

Close modal

The permeability as a function of the effective stress applied is obtained by plotting the linear relationship between the discharge velocity U and the hydraulic gradient i; the slope of the line becomes the permeability, as shown in Figure 4. For each effective stress, the hydraulic anisotropy ratio rk is calculated from Equation 1. After each test, the sample is retrieved and dried for 24 h, and it is then analysed by wet sieving [ASTM D1140-17 (ASTM, 2017b)] to determine its fines content.

Figure 4.
Two graphs depict discharge velocity against hydraulic gradient, with plotted lines representing different stress values. Each line is labelled with a linear equation.The image contains two graphs, labelled A and B, depicting the relationship between discharge velocity, presented in units of metres per second, and hydraulic gradient, expressed as a numeric value. Each graph features multiple lines, colour-coded to represent different stress values, ranging from 19.6 kilopascals to 313.8 kilopascals. Each line is accompanied by its corresponding linear equation written in scientific notation. Both graphs have the same axis scales, with the hydraulic gradient shown on the X-axis labelled from 10 to 90, and discharge velocity on the Y-axis labelled from 0 to 6 times 10 to the power of negative 7, with incremental markings.

Determination of the permeability for each applied stress (sample 3): (a) evaluation of kv and (b) evaluation of kh

Figure 4.
Two graphs depict discharge velocity against hydraulic gradient, with plotted lines representing different stress values. Each line is labelled with a linear equation.The image contains two graphs, labelled A and B, depicting the relationship between discharge velocity, presented in units of metres per second, and hydraulic gradient, expressed as a numeric value. Each graph features multiple lines, colour-coded to represent different stress values, ranging from 19.6 kilopascals to 313.8 kilopascals. Each line is accompanied by its corresponding linear equation written in scientific notation. Both graphs have the same axis scales, with the hydraulic gradient shown on the X-axis labelled from 10 to 90, and discharge velocity on the Y-axis labelled from 0 to 6 times 10 to the power of negative 7, with incremental markings.

Determination of the permeability for each applied stress (sample 3): (a) evaluation of kv and (b) evaluation of kh

Close modal

Compressibility is a measure of the volume change of tailings in response to an increase in stress. Understanding compressibility parameters (such as the compression index, Cc, or the recompression index, Cs) is essential for the design and operation of tailings dams, as it allows predicting the magnitude and rate of settlements and estimating the dissipation of pore pressures. These factors are crucial for ensuring the static stability of dams, optimising drainage systems, planning growth stages, and preventing failures associated with excessive deformations or critical outlet gradients, both during operation and during dam closure (Bussiere, 2004; Hu et al., 2017).

Hu et al. (2017) analysed the geotechnical properties of four tailings mixtures from the Yuhezhai and Bahuerachi mines, which were tested at stresses of 12.5, 25, 50, 100, 200, 400, 800, and 1600 kPa. For the tailings from the Yuhezhai mine, they obtained Cc values of 0.046 and 0.260, while for the samples from the Bahuerachi mine, they obtained Cc values of 0.025 and 0.085. In another study conducted by Fan et al. (2022), the influence of fines content on the compressibility of tailings mixtures with variable fines content, from 0% to 100%, was determined. The samples were tested at stresses of 90, 180, 390, and 610 kPa, for tailings with between 0 and 50% fines, Cc values between 0.051 and 0.085 were obtained, and for tailings with fines greater than 50%, Cc values between 0.091 and 0.15 were obtained. Figure 5 illustrates the effect of fines content on the compression index (Cc) for the samples tested in this study. For relatively dense samples (e ≤ 1), Cc shows a clear increasing trend with fines content: the sample with the lowest fines content (Fc = 13%) exhibited a Cc = 0.056, whereas for samples with Fc = 100%, Cc ranged from 0.065 to 0.102. In contrast, for looser tailings samples with void ratios greater than one, Cc values ranged from 0.194 to 0.440, indicating significantly higher compressibility. This behaviour is consistent with the findings of Hu et al. (2017) and Fan et al. (2022), who also reported that the compression index increases with fines content. This trend reflects the greater compressibility of fine-grained tailings.

Figure 5.
A scatter plot depicts the relationship between fines content percent on the X-axis and compression ratio on the Y-axis, with different sample types indicated by colours and shapes.This scatter plot illustrates the relationship between the percentage of fines content, labelled as F c, on the X-axis, which ranges from 0 to 120 percent, and the compression ratio, denoted as C c, on the Y-axis, spanning from 0 to a maximum of 0.45. The data points are represented by different shapes and colours, blue circles for Vertical samples, red circles for Horizontal samples, and black diamonds indicating specific studies, Fan et al., 2022 and Hu et al., 2017. A dashed red horizontal line at C c equals 0.15 serves as a reference, separating the areas where e is greater than 1 and less than 1, indicated by vertical arrows pointing to the respective regions. The distribution of the data points shows varying compression ratios across different fines content levels, with some samples clustering below the reference line and others above it.

Fines content against compression index for consolidated samples from 19.6 to 313.8 kPa

Figure 5.
A scatter plot depicts the relationship between fines content percent on the X-axis and compression ratio on the Y-axis, with different sample types indicated by colours and shapes.This scatter plot illustrates the relationship between the percentage of fines content, labelled as F c, on the X-axis, which ranges from 0 to 120 percent, and the compression ratio, denoted as C c, on the Y-axis, spanning from 0 to a maximum of 0.45. The data points are represented by different shapes and colours, blue circles for Vertical samples, red circles for Horizontal samples, and black diamonds indicating specific studies, Fan et al., 2022 and Hu et al., 2017. A dashed red horizontal line at C c equals 0.15 serves as a reference, separating the areas where e is greater than 1 and less than 1, indicated by vertical arrows pointing to the respective regions. The distribution of the data points shows varying compression ratios across different fines content levels, with some samples clustering below the reference line and others above it.

Fines content against compression index for consolidated samples from 19.6 to 313.8 kPa

Close modal

Tailings consolidation results from particle deformation and the expulsion of water or air from voids and is governed by permeability and compressibility (Mosquera, 2013; Vermeulen, 2001; Vick, 1990).

According to Terzaghi, the consolidation process can be divided into primary and secondary phases. Primary consolidation controls the rate of dissipation of excess pore pressures under constant loading, while secondary consolidation is additional deformation under constant loading, even after pore pressure has been completely dissipated, and is associated with grain-to-grain sliding and particle rearrangement (Terzaghi et al., 1996). Primary and secondary consolidation occur simultaneously; however, most of the effects of secondary consolidation in tailings are very small and are often considered insignificant (Vick, 1990; Xu, 2019).

According to Saad (2008), sand tailings (coarse fraction) dissipate excess pore pressure quickly since they usually have high permeability and lower compressibility (according to the results presented in the previous section); consequently, the values of the consolidation coefficient Cv are relatively high. On the other hand, in saturated fine tailings, when subjected to an increase in load, the excess pore pressure generated dissipates much more slowly due to their low permeability and high compressibility, which causes the Cv to decrease. This can be seen in the results shown in Figure 6(a) for the vertical (V-S) and horizontal (H-S) samples; for tailings with fines contents less than 50%, Cv values of 2.6 × 10−2 to 5 × 10−5 m2/s were obtained; while for tailings with fines contents greater than 50%, Cv values of the range of 4.5 × 10−4 to 4.4 × 10−7 m2/s were obtained.

Figure 6.
Scatter plots depict the relationship between the coefficient of consolidation and 2 variables, fines content and void ratio, with varying symbols for different data sets.The image consists of two scatter plots labelled a and b illustrating the relationship between the coefficient of consolidation, c v, and two variables, fines content, F c, on the left and void ratio, e, on the right. The Y-axis in both plots represents the coefficient of consolidation measured in square metres per second, ranging from 1 times 10 to the power of negative 7 to 1 times 10 to the power of 1. The X-axis for plot a indicates fines content, which ranges from 0 to 120 percent, while plot b shows the void ratio ranging from 0 to 2. Each data point is presented using different symbols and colours, as described in the legend below the plots, distinguishing various data sets such as V S 1, H S 1, V S 2, and H S 3, among others.

(a) Effect of fines content on the coefficient of consolidation and (b) effect of void ratio on the coefficient of consolidation

Figure 6.
Scatter plots depict the relationship between the coefficient of consolidation and 2 variables, fines content and void ratio, with varying symbols for different data sets.The image consists of two scatter plots labelled a and b illustrating the relationship between the coefficient of consolidation, c v, and two variables, fines content, F c, on the left and void ratio, e, on the right. The Y-axis in both plots represents the coefficient of consolidation measured in square metres per second, ranging from 1 times 10 to the power of negative 7 to 1 times 10 to the power of 1. The X-axis for plot a indicates fines content, which ranges from 0 to 120 percent, while plot b shows the void ratio ranging from 0 to 2. Each data point is presented using different symbols and colours, as described in the legend below the plots, distinguishing various data sets such as V S 1, H S 1, V S 2, and H S 3, among others.

(a) Effect of fines content on the coefficient of consolidation and (b) effect of void ratio on the coefficient of consolidation

Close modal

During consolidation, the void ratio decreases, and permeability and volumetric compressibility typically decrease. The opposing interaction of these two functions determines whether Cv remains constant, increases, or decreases with the void ratio (Fan et al., 2022; Wang et al., 2015). Figure 6(b) shows the influence of the void ratio on Cv for the tailings tested.

The works by Bussiere (2004), Mosquera (2013), Vermeulen (2001), Vick (1990), and Qiu and Sego (2001) report values for the coefficients Cc and Cv for different types of tailings, both for undisturbed and reconstituted samples.

Cubrinovski and Ishihara (2002) used data from more than 300 natural sands (clean sands, sands with fines, and sands with small amounts of clay-sized particles) to assess the influence of fines, grain size composition, and particle shape on the void ratio (emaxemin). They found that both emax and emin for mixtures of sand with fines decrease as the fines content increases from 0% to about 20%, as small particles fill the gaps between the larger particles. Within the range from 20% to 40% of fines, emax and emin show a pattern change indicating a transition from the void-filling process to the solid particle replacement process. In mixtures of sands with fines content greater than 40%, emax and emin constantly increase until they reach the highest values corresponding to a fines content of 100%.

Figure 7(a) shows undisturbed tailings samples with fines contents ranging from 13% to 100%, with their initial void ratio (corresponding to σ′ = 19.6 kPa) and final void ratio (σ′ = 313.8 kPa); it can be observed that the tailings show a behaviour similar to that described by Cubrinovski and Ishihara (2002). The void ratios (ei and ef) tend to decrease as Fc increases in tailings samples with very low fines content (between 0% and 30%); then, there is a transition zone in the tailings samples with Fc between 40% and 50%, in which ei and ef present their lowest value and begin to increase, until they reach the highest values for tailings with Fc equal to 100%.

Figure 7.
A two-panel plot depicts the relationship between void ratio and fines content in the first panel, and depth in the second panel. Data points represent various measurements with distinct shapes and symbols.The image features a two panel scatter plot. The left panel A charts void ratio on the Y-axis against fines content, expressed as a percent, on the X-axis. Various data points are represented by different shapes and colours, corresponding to distinct sample designations listed below the plot. The curves, one for a pressure of 19.6 kilopascals and another for 313.8 kilopascals, depict how void ratio varies with fines content, indicating two relationships. The right panel B depicts void ratio against depth in metres, with similar distinct data point representations. This arrangement facilitates comparison across different conditions and sample types.

(a) Variation of ei and ef in terms of fines content for undisturbed tailings samples and (b) variation of ei and ef for undisturbed tailings samples with fines content ranging from 80% to 100%

Figure 7.
A two-panel plot depicts the relationship between void ratio and fines content in the first panel, and depth in the second panel. Data points represent various measurements with distinct shapes and symbols.The image features a two panel scatter plot. The left panel A charts void ratio on the Y-axis against fines content, expressed as a percent, on the X-axis. Various data points are represented by different shapes and colours, corresponding to distinct sample designations listed below the plot. The curves, one for a pressure of 19.6 kilopascals and another for 313.8 kilopascals, depict how void ratio varies with fines content, indicating two relationships. The right panel B depicts void ratio against depth in metres, with similar distinct data point representations. This arrangement facilitates comparison across different conditions and sample types.

(a) Variation of ei and ef in terms of fines content for undisturbed tailings samples and (b) variation of ei and ef for undisturbed tailings samples with fines content ranging from 80% to 100%

Close modal

The highest values of ei and ef are related to soils with 100% fines, as can be seen in Figure 7(a). According to Vick (1990), lower void ratios generally correspond to greater depths within a tailings dam, and higher void ratios are generally associated with near-surface materials. Figure 7(b) shows the values of ei and ef for tailings with fine contents between 80% and 100%; samples 1, 2, and 3 have the highest void ratios corresponding to the shallowest samples.

As previously mentioned, higher permeability favours elevated Cv values and, consequently, faster consolidation, whereas low permeabilities result in reduced Cv and prolonged dissipation times (Qiu and Sego, 2001). This relationship between Cv and k has direct implications for the stability of tailings deposits; tailings with high permeability and low compressibility promote rapid dissipation of pore pressures, allowing a more stable response under loading. In contrast, tailings with low permeability and reduced Cv cause excess pore pressure to persist for long periods, thereby decreasing effective stresses and increasing the likelihood of structural instability and even static liquefaction phenomena (Mosquera, 2013; Qiu and Sego, 2001; Vermeulen, 2001). Therefore, jointly analysing Cv and k is essential to design drainage strategies, define safe loading rates, and evaluate the short- and long-term stability of tailings dams (Vermeulen, 2001). Figure 8 illustrates this relationship for the tested samples; in tailings with high permeabilities (between 1 × 10−4 and 1 × 10−6 m/s; Fell, 2005), Cv values range from 9.4 × 10−3 to 4.05 × 10−3 m2/s, while in low-permeability tailings Cv values decrease to a range of 4.0 × 10−4 to 5.9 × 10−7 m2/s.

Figure 8.
A scatter plot depicts the relationship between permeability and the coefficient of consolidation, with various symbols representing different data points.The image displays a scatter plot illustrating the relationship between permeability, measured in metres per second, on the Y-axis, and the coefficient of consolidation, measured in square metres per second, on the X-axis. The plot features various symbols, including squares, diamonds, and triangles, each representing different data sets. The symbols are colour coded to indicate specific groups, with a legend at the bottom that identifies each group by its respective symbol and colour. The axes are marked with exponential scales, noting ranges from 1 times 10 to the power of negative 10 to 1 times 10 to the power of 0 on the Y-axis, and from 1 times 10 to the power of negative 7 to 1 times 10 to the power of 0 on the X-axis. The plot has no gridlines, emphasising the distribution of data points across the graph.

Effect of permeability on the coefficient of consolidation

Figure 8.
A scatter plot depicts the relationship between permeability and the coefficient of consolidation, with various symbols representing different data points.The image displays a scatter plot illustrating the relationship between permeability, measured in metres per second, on the Y-axis, and the coefficient of consolidation, measured in square metres per second, on the X-axis. The plot features various symbols, including squares, diamonds, and triangles, each representing different data sets. The symbols are colour coded to indicate specific groups, with a legend at the bottom that identifies each group by its respective symbol and colour. The axes are marked with exponential scales, noting ranges from 1 times 10 to the power of negative 10 to 1 times 10 to the power of 0 on the Y-axis, and from 1 times 10 to the power of negative 7 to 1 times 10 to the power of 0 on the X-axis. The plot has no gridlines, emphasising the distribution of data points across the graph.

Effect of permeability on the coefficient of consolidation

Close modal

Tailings permeability is difficult to generalise; on the average, it varies from five or more orders of magnitude, that is from 10−4 m/s for clean coarse sands to 10−9 m/s for consolidated sludge (Vick, 1990). There are important aspects that require consideration, such as the influence of the void ratio, grain size distribution, plasticity, fines content, deposition method, and dam thickness (Mittal, 1975; Vick, 1990). From a geotechnical point of view, tailings are considered to have low to very low hydraulic conductivity (Lambe, 1958).

Hydraulic deposition affects the properties of the tailings, so that the parameters change with distance to the discharge point; material near the discharge point has a void ratio ranging from 0.6 to 1.0, as opposed to the furthest points where the void ratio varies from 1.0 to 1.6 or more (Espinoza, 2005; Jantzer et al., 2008). The tailings generally evidence a loose state (Blight, 2010).

Tailings are more compressible compared with natural soils of equivalent grain size distribution, due to the loose state of deposition and their high angularity; therefore, tailings require measurements of permeability at various stress levels to estimate reductions in permeability as a function of the void ratio (Vick, 1990).

Effective stress and void ratio are two factors that are associated to and directly influence permeability. As the confinement stress in the sample increases, the specimen becomes densified, and its volume of voids decreases. As the volume of voids decreases, the permeability in the sample decreases too, since the area through which the water flows decreases. There is a linear relationship between the void ratio and the effective stress applied to the tailing samples (Qiu and Sego, 2001), as well as between permeability and the effective stress in tailings (Ma et al., 2023). The vertical and horizontal samples tested show these linear relationships. As it can be seen in Figures 9(a) and 9(b), the tailings specimens tested presented void ratios of 1.5–0.56 in their loosest condition (σ′ = 19.6 kPa) and void ratios between 1.23 and 0.48 for the densest condition (σ ′ = 313.8 kPa). For both conditions, permeabilities ranging from 1 × 10−4 to 1 × 10−9 m/s were obtained.

Figure 9.
Graphs depict void ratio and permeability as functions of effective stress, with different markers representing various samples. The data is presented in 2 separate plots.The image features two graphs that depict the relationship between void ratio and permeability against effective stress. The left graph, labelled a, plots void ratio on the Y-axis and effective stress on the X-axis, with values ranging from 10 to 1000 kilopascals. Various symbols in different colours represent distinct samples, including filled and unfilled circles, triangles, and diamonds. The right graph, labelled b, plots permeability on the Y-axis, with a scale from 1 times 10 to the power of negative 9 to 1 times 10 to the power of negative 4 metres per second, against effective stress on the X-axis. Both graphs include a legend explaining the marker symbols for each sample, with labels such as V S 1, H S 1, and H S 12 2. The layout allows comparison across the two datasets using consistent axis arrangement and presentation style.

(a) Effect of void ratio on permeability in undisturbed tailings samples and (b) effect of effective stress on permeability in undisturbed tailings samples

Figure 9.
Graphs depict void ratio and permeability as functions of effective stress, with different markers representing various samples. The data is presented in 2 separate plots.The image features two graphs that depict the relationship between void ratio and permeability against effective stress. The left graph, labelled a, plots void ratio on the Y-axis and effective stress on the X-axis, with values ranging from 10 to 1000 kilopascals. Various symbols in different colours represent distinct samples, including filled and unfilled circles, triangles, and diamonds. The right graph, labelled b, plots permeability on the Y-axis, with a scale from 1 times 10 to the power of negative 9 to 1 times 10 to the power of negative 4 metres per second, against effective stress on the X-axis. Both graphs include a legend explaining the marker symbols for each sample, with labels such as V S 1, H S 1, and H S 12 2. The layout allows comparison across the two datasets using consistent axis arrangement and presentation style.

(a) Effect of void ratio on permeability in undisturbed tailings samples and (b) effect of effective stress on permeability in undisturbed tailings samples

Close modal

With improvements in beneficiation technology, ores are finely ground, and the fines content of the tailings is increasing. Fine-grained tailings have low permeability, a long consolidation time, and a slow increase in strength, which implies high safety risks to the dams (Ma et al., 2023; Zhang et al., 2020).

Valenzuela (2015) presented results on saturated copper tailings from different mines in Chile and Peru. The most remarkable differences observed were in tailings with fines contents between 30% and 40%, probably because of a change in the soil fabric, since fines can occupy the spaces between the larger particles, thus decreasing the size of the voids as the size of the fines is reduced.

Figure 10 shows the kv and kh permeabilities obtained for tailings samples with fines contents between 12% and 100%, as well as the results and boundaries proposed by Valenzuela (2015). It is important to mention that these limits are exclusively for kv.

Figure 10.
A scatter plot depicts permeability against void ratio, with distinct data points categorised by shape and colour, and labelled lines indicating ranges of fines.The image displays a scatter plot where the X-axis represents void ratio, labelled as Void ratio, e, and the Y-axis represents permeability, labelled as Permeability, k, metres per second. The data points vary in shape and colour, including circles, diamonds, squares, and triangles, representing different datasets. A grey background indicates points from Valenzuela, 2015, while other points are identified by codes such as H S 1 through H S 14 and V S 1 through V S 14. The plot includes two black lines that demarcate regions for 15-30 percent fines and 40-100 percent fines, illustrating groupings based on permeability and void ratio relationships. The permeability values range from approximately 1 times 10 to the power of negative 10 to 1 times 10 to the power of negative 4 along the Y-axis.

Permeability as a function of the fines content for undisturbed tailings samples

Figure 10.
A scatter plot depicts permeability against void ratio, with distinct data points categorised by shape and colour, and labelled lines indicating ranges of fines.The image displays a scatter plot where the X-axis represents void ratio, labelled as Void ratio, e, and the Y-axis represents permeability, labelled as Permeability, k, metres per second. The data points vary in shape and colour, including circles, diamonds, squares, and triangles, representing different datasets. A grey background indicates points from Valenzuela, 2015, while other points are identified by codes such as H S 1 through H S 14 and V S 1 through V S 14. The plot includes two black lines that demarcate regions for 15-30 percent fines and 40-100 percent fines, illustrating groupings based on permeability and void ratio relationships. The permeability values range from approximately 1 times 10 to the power of negative 10 to 1 times 10 to the power of negative 4 along the Y-axis.

Permeability as a function of the fines content for undisturbed tailings samples

Close modal

Increasing the fines content in a tailings dam reduces permeability and jeopardises the stability of the dam (Singh et al., 2021; Saad, 2008). The effect of the fine content Fc on the permeability (kv and kh) for all test conditions, that is for the six stress increments (19.6, 39.2, 78.5, 156.9, 245.2, and 313.8 kPa), is shown in Figure 11, permeability in tailings samples tends to decrease as fines increase, regardless of the orientation of the test. The permeability in the vertical sample with the lowest fines content (Fc = 12.2%) is of about 1 × 10−5 m/s, while for the samples with the highest fines content (Fc = 100%), permeabilities of the order of 1 × 10−9 m/s were obtained. As for the horizontal samples, they showed a permeability of 1 × 10−5 m/s in the sample with a fines content of 18.3% and from 1 × 10−8 to 1 × 10−9 cm/s for the specimen with fines content of 100%.

Figure 11.
A graph depicts the relationship between fines content percent and permeability, showing two types of samples, Vertical and Horizontal, plotted to compare their permeability behaviour.The graph depicts the relationship between fines content, indicated on the X-axis as percent from 0 to 100, and permeability measured in kilopascals per metre per second on the Y-axis, ranging from 1 times 10 to the power of negative 10 to 1 times 10 to the power of negative 4. Two sets of data points represent Vertical samples in blue and Horizontal samples in red, with dotted lines indicating the trend for each sample type. Bands of light shading accompany the lines to denote uncertainty or variability in the data. The layout moves left to right for fines content and bottom to top for permeability.

Influence of the fine contents in kv and kh in undisturbed tailings samples

Figure 11.
A graph depicts the relationship between fines content percent and permeability, showing two types of samples, Vertical and Horizontal, plotted to compare their permeability behaviour.The graph depicts the relationship between fines content, indicated on the X-axis as percent from 0 to 100, and permeability measured in kilopascals per metre per second on the Y-axis, ranging from 1 times 10 to the power of negative 10 to 1 times 10 to the power of negative 4. Two sets of data points represent Vertical samples in blue and Horizontal samples in red, with dotted lines indicating the trend for each sample type. Bands of light shading accompany the lines to denote uncertainty or variability in the data. The layout moves left to right for fines content and bottom to top for permeability.

Influence of the fine contents in kv and kh in undisturbed tailings samples

Close modal

Due to the stratified nature of the tailings, differences were found in the fines content among samples trimmed either vertically or horizontally, extracted from the same Shelby tube. Because the amount of fines in adjacent samples of the same tube is not always equal, the term ‘difference in fines content’ (ΔFc) was defined as follows:

2

where Fchs is the fines content in the horizontal sample, and Fcvs is the content of fines in the vertical sample.

Figure 12 shows the influence of the difference in fines content (ΔFc) on the hydraulic anisotropy ratio (rk). The hydraulic anisotropy ratio tends to increase when the fines content in the horizontal sample is lower than the fines content in the vertical sample, and the rk value can be above 100. The opposite is the case when the fines content in the horizontal sample is greater than the fines content in the vertical sample; rk decreases, and values of rk less than unity may occur.

Figure 12.
A scatter plot depicts the relationship between permeability anisotropy and the difference in fines content for various samples, labelled with different shapes and colour.The image depicts a scatter plot illustrating the relationship between permeability anisotropy on the Y-axis and the difference in fines content on the X-axis. The Y-axis uses a logarithmic scale ranging from 1 to 1000, while the X-axis ranges from negative 30 to positive 30. Different samples are represented using varied shapes and colours, with a legend at the bottom identifying each sample number and its corresponding marker. Directional arrows indicate the conditions F c s h greater than F c s v and F c s h less than F c s v, along with annotations identifying specific sample groupings. The data points are distributed across the plot, showing varying permeability anisotropy values relative to differences in fines content.

Influence of the difference in the fine contents on the hydraulic anisotropy ratio

Figure 12.
A scatter plot depicts the relationship between permeability anisotropy and the difference in fines content for various samples, labelled with different shapes and colour.The image depicts a scatter plot illustrating the relationship between permeability anisotropy on the Y-axis and the difference in fines content on the X-axis. The Y-axis uses a logarithmic scale ranging from 1 to 1000, while the X-axis ranges from negative 30 to positive 30. Different samples are represented using varied shapes and colours, with a legend at the bottom identifying each sample number and its corresponding marker. Directional arrows indicate the conditions F c s h greater than F c s v and F c s h less than F c s v, along with annotations identifying specific sample groupings. The data points are distributed across the plot, showing varying permeability anisotropy values relative to differences in fines content.

Influence of the difference in the fine contents on the hydraulic anisotropy ratio

Close modal

As already mentioned, in the literature, hydraulic anisotropy ratio values are found in tailings varying from 2 to 100. In the present study, for most of the samples tested, hydraulic anisotropy ratios were found in the intervals of 1 to 10; as shown in Figure 13, hydraulic anisotropy ratios above this range are due to a significant difference in fines content (ΔFc > 10%) of the samples tested.

Figure 13.
A scatter plot showing the relationship between horizontal and vertical permeability, with experimental data points and shaded regions indicating different permeability ranges.The scatter plot illustrates the relationship between horizontal permeability, labelled k h in meters per second, on the horizontal axis and vertical permeability, labelled k v in meters per second, on the vertical axis. The axes are logarithmic, with values ranging from one times ten to the power of negative nine to one times ten to the power of negative three for both axes. The experimental data points are represented by black dots. There are several shaded regions indicating varying permeability categories: yellow for k less than one, light blue for one less than or equal to k less than or equal to ten, pink for ten less than or equal to k less than or equal to one hundred, and green for k greater than one hundred. The regions are separated by dashed lines that outline the different ranges. The plot provides a clear visual representation of how permeability varies in relation to the two types of measurements.

Hydraulic anisotropy ratio in tailings

Figure 13.
A scatter plot showing the relationship between horizontal and vertical permeability, with experimental data points and shaded regions indicating different permeability ranges.The scatter plot illustrates the relationship between horizontal permeability, labelled k h in meters per second, on the horizontal axis and vertical permeability, labelled k v in meters per second, on the vertical axis. The axes are logarithmic, with values ranging from one times ten to the power of negative nine to one times ten to the power of negative three for both axes. The experimental data points are represented by black dots. There are several shaded regions indicating varying permeability categories: yellow for k less than one, light blue for one less than or equal to k less than or equal to ten, pink for ten less than or equal to k less than or equal to one hundred, and green for k greater than one hundred. The regions are separated by dashed lines that outline the different ranges. The plot provides a clear visual representation of how permeability varies in relation to the two types of measurements.

Hydraulic anisotropy ratio in tailings

Close modal

In this work, the results of 31 permeability tests (14 in a vertical direction and 17 in a horizontal direction) performed to evaluate the hydraulic anisotropy ratio rk in undisturbed tailings samples were presented. The main findings of the research are summarised as follows:

  • In the literature, it is possible to find various values of rk for clays, sands, and rocks, but the results available for tailings are scarce. In this paper, the experimental procedure to obtain rk in undisturbed tailings specimens subjected to permeability tests carried out on a constant head permeameter using the methodology of sample cutting in different directions has been presented. The samples were tested in a saturated state, because the test is non-destructive, six levels of effective stress equal to 19.6, 39.2, 78.5, 156.9, 245.2, and 313.8 kPa were applied to each specimen. For each stress, the vertical or horizontal permeability was obtained by applying hydraulic heads (h) of 150, 217, 284, and 350 cm of water column.

  • The experimental method of cutting samples in different directions presented by Chapuis and Gill (1989) has been used to determine rk in materials such as clays and rocks, but it is not known that this methodology has been used to determine rk in tailings, which was implemented in this research.

  • From the consolidation stage, the values of the Cc and Cv coefficients were obtained for the different samples tested. The influence of the fines content on Cc was evaluated; for samples with e ≤ 1, the Cc increases with the fines content from 0.056 to 0.102; for tailings samples with void ratios greater than unity, the Cc varies from 0.194 to 0.44.

  • The effect of the void ratio and the fines content on the Cv coefficient was also evaluated; for tailings with fines contents less than 50%, Cv values of 2.6 × 10−2 to 5 × 10−5 m2/s were obtained; while for tailings with fine content greater than 50%, Cv values of the order of 4.5 × 10−4 to 4.4 × 10−7 m2/s were obtained.

  • In the experimental results presented in this work, it was shown that the fines content is a factor that has a significant influence on rk. Therefore, because the tailings have a high stratification, it is recommended to obtain the Fc value of each of the specimens after determining kv and kh.

  • Tailings were found to have a behavior similar to that described by Cubrinovski and Ishihara (2002) for mixtures of sand with fines. It was determined that the relationship of initial and final void ratios (ei and ef) tends to decrease as Fc increases in tailings samples with very low fines contents (between 0% and 30%); then, there is a transition zone in tailings with Fc ranging from 40% to 50%, in which ei and ef present their lowest value and begin to increase until reaching the highest values of void ratio for tailings with an Fc equal to 100%.

  • The tailings studied with fines contents between 80% and 100% showed a significant variation in ei and ef (from 1.4 to 0.5); the void ratios exceeding unity belong to the most superficial samples tested. This is like what was found by Vick (1990), who pointed out that the highest void ratios are associated with tailing samples close to the surface of the dam.

  • According to the literature, higher permeability favours high Cv values, while low permeabilities result in reduced Cv; the results obtained in this study show that for tailings with permeabilities between 1 × 10−4 and 1 × 10−6 m/s, Cv values vary from 9.4 × 10−3 to 4.05 × 10−3 m2/s, while in low permeability tailings, Cv values decrease to a range of 4.0 × 10−4 to 5.9 × 10−7 m2/s.

  • The linear relationship between the void ratio and the effective stress applied in the tailings samples was verified, as well as that between permeability and effective stress; the tailings samples tested presented void ratios from 1.5 to 0.56 in their loosest condition (σ′ = 19.6 kPa) and void ratios between 1.23 and 0.48 for the densest condition (313.8 kPa). For both conditions, permeability values of the order of 1 × 10−4 to 1 × 10−9 m/s were determined.

  • The permeability results as a function of the fines content were located within the boundaries established by Valenzuela (2015) for tailings with 15% and 30% fines content, and for tailings having from 40% to 100% of fines content.

  • From the test results, it was determined that the permeability in the tailings samples tends to decrease as the fines content increases, regardless of the orientation of the test specimens. The permeability in the vertical specimen with the lowest fines content (Fc = 12.2%) was in the order of 1 × 10−5 m/s, while for the samples with the highest fines content (Fc = 100%) permeabilities obtained were in the order of 1 × 10−9 m/s. As for the horizontal samples, they presented a permeability of 1 × 10−5 m/s in the specimen with fines content of 18.3% and from 1 × 10−8 to 1 × 10−9 cm/s for the sample with fines content equal to 100%.

  • The vertical and horizontal samples tested to evaluate rk are from the same Shelby tube, but do not correspond to the same horizon; due to this and to the highly stratified nature of the tailings, in most of the tests, a difference in fines content was obtained between a vertical and a horizontal specimen. Therefore, the term difference in fines contentFc) was defined, representing the difference in the fines content of the horizontal sample and the fines content of the vertical specimen.

  • It was observed that the hydraulic anisotropy ratio tends to increase when the fines content in the horizontal sample is smaller than the fines content in the vertical specimen, and the rk value can exceed 100. The opposite case arises when the fines content in the horizontal sample is larger than the fines content in the vertical specimen; rk decreases, and values of rk less than unity may be found.

  • Values of hydraulic anisotropy ratios, rk, for tailings from 2 to 100 are found in the literature. In this study, for most of the samples tested, hydraulic anisotropy ratios were found to be in the range from 1 to 10; hydraulic anisotropy ratios above this range are due to a difference in fines content between samples (ΔFc > 10%).

  • Nowadays, the accuracy required in water flow models has been increasing, so it is necessary to determine realistic values of the hydraulic anisotropy ratio in tailings, such as those found in the present work, since the hydraulic anisotropy ratio within the tailings dam has an influence on aspects that can even affect the general stability of the dam.

Abadjiev
CB
(
1976
)
Seepage through mill tailings dams
.
Commission Internationale Des Grands Barrages. México
1
(R42)
:
381
393
.
Adams
AL
,
Nordquist
TJ
,
Germaine
JT
and
Flemings
PB
(
2016
)
Permeability anisotropy and resistivity anisotropy of mechanically compressed mudrocks
.
Canadian Geotechnical Journal
53
(9)
:
1474
1482
.
ASTM
(
2017
b) ASTM D1140-17 Standard test methods for determining the amount of material finer than 75-μm (no. 200) sieve in soils by washing.
ASTM
,
West Conshohocken, PA
.
ASTM
(
2019
) ASTM D2216-19 Standard Test Method for Laboratory Determination of Water (Moisture) Content of Soil and Rock by Mass.
ASTM
,
West Conshohocken, PA
.
ASTM
(
2020
) ASTM D4767-11 Standard Test Method for Consolidated Undrained Triaxial Compression Test for Cohesive Soils.
ASTM
,
West Conshohocken, PA
.
ASTM
(
2023
) ASTM D854-23 Standard test methods for specific gravity of soil solids by water pycnometer
ASTM
,
West Conshohocken, PA
.
ASTM
(
2024
) ASTM D5084 − 24 Standard Test Methods for Measurement of Hydraulic Conductivity of Saturated Porous Materials Using a Flexible Wall Permeameter.
ASTM
,
West Conshohocken, PA
.
Blight
GE
(
2010
)
Geotechnical Engineering for Mine Waste Storage Facilities
.
CRC Press
.
Bolton
AJ
,
Maltman
AJ
and
Fisher
Q
(
2000
)
Anisotropic permeability and bimodal pore-size distributions of fine-grained marine sediments
.
Marine and Petroleum Geology
17
(6)
:
657
672
, .
Bussiere
B
(
2004
)
Colloquium 2004: Hydrogeotechnical properties of hard rock tailings from metal mines and emerging geoenvironmental disposal approaches
.
Canadian Geotechnical Journal
.
Candelaria
J
,
Promotor
A
,
Medina
S
et al.
(
2022
)
Aspectos Relevantes Para el Análisis de Estabilidad de Un Depósito De Jales Mediante El Método De Equilibrio Límite
.
XXXI Reunión Nacional de Ingeniería Geotécnica
.
Chan
HT
and
Kenney
TC
(
1973
)
Laboratory investigation of permeability ratio of new liskeard varved soil
.
Canadian Geotechnical Journal
10
(3)
:
453
472
, .
Chapuis
RP
and
Gill
DE
(
1989
)
Hydraulic anisotropy of homogeneous soils and rocks: influence of the densification process
.
Bulletin of the International Association of Engineering Geology
39
(1)
:
75
86
, .
Clavaud
J
,
Maineult
A
,
Zamora
M
,
Rasolofosaon
P
and
Schlitter
C
(
2008
)
Permeability anisotropy and its relations with porous medium structure
.
Journal of Geophysical Research: Solid Earth
113
(B1)
, .
Cubrinovski
M
and
Ishihara
K
(
2002
)
Maximum and minimum void ratio characteristics of sands
.
Soils and Foundations
42
(6)
:
65
78
, .
Espinoza
I
(
2005
)
Análisis Del Comportamiento de Residuos Mineros a Partir de Estudios Experimentales
.
UNAM
.
Fan
J
,
Rowe
RK
and
Brachman
RWI
(
2022
)
Compressibility and permeability of sand–silt tailings mixtures
.
Canadian Geotechnical Journal
59
(8)
:
1348
1357
, .
Fell
R
(
2005
)
Geotechnical Engineering of Dams
.
CRC press
.
Gu
X
,
Liang
X
and
Hu
J
(
2022
)
Quantifying fabric anisotropy of granular materials using wave velocity anisotropy: a numerical investigation
.
Géotechnique
74
(12)
:
1263
1275
, .
Hu
L
,
Wu
H
,
Zhang
L
,
Zhang
P
and
Wen
Q
(
2017
)
Geotechnical properties of mine tailings
.
Journal of Materials in Civil Engineering
29
(2)
:
04016220
, .
ICOLD
(
2001
)
Tailings Dams – Risk of Dangerous Occurrences: Lessons Learnt from Practical Experiences (Bulletin 121)
.
International Commission on Large Dams
.
Jamshidi Chenari
R
and
Mahigir
A
(
2014
)
The effect of spatial variability and anisotropy of soils on bearing capacity of shallow foundations
.
Civil Engineering Infrastructures Journal
47
(2)
:
199
213
.
Jamshidi Chenari
R
and
Behfar
B
(
2017
)
Stochastic analysis of seepage through natural alluvial deposits considering mechanical anisotropy
.
Civil Engineering Infrastructures Journal
50
(2)
:
233
253
.
Jantzer
I
,
Bjelkevik
A
and
Pousette
K
(
2008
)
Material Properties of Tailings from Swedish Mines
.
Academia
.
Kossoff
D
,
Dubbin
WE
,
Alfredsson
M
et al.
(
2014
)
Mine tailings dams: Characteristics, failure, environmental impacts, and remediation
.
Applied Geochemistry
51
:
229
245
.
Lambe
TW
(
1958
)
The structure of compacted clay
.
Jnl. of the Soil Mech. and Foundn. Div ASCE
84
(SM2)
:
10
34
.
López-Acosta
NP
and
Mendoza-Promotor
JA
(
2016
)
Flujo de Agua En Suelos Parcialmente Saturados y su Aplicación a la Ingeniería Geotécnica
.
Serie Investigación y Desarrollo del Instituto de Ingeniería, UNAM
.
López-Acosta
NP
(
2017
)
Manual de Diseño de Obras Civiles de la Comisión Federal de Electricidad (CFE), Sección B. Geotecnia. Tema 2. Mecánica de Suelos. Capítulo 9. “Flujo De Agua En Suelos”
,
México
.
Lucero
R
(
2022
)
Permeabilidad y su Relación de Anisotropía Obtenida Experimentalmente En Suelos Finos
.
Tesis de maestría, Universidad Nacional Autónoma de México
,
México D.F., México
.
Ma
C
,
Zhang
C
,
Guo
X
,
Zhou
H
and
Gan
S
(
2023
)
High-stress permeability and consolidation characteristics of flocculated fine tailings
.
Proceedings of the Institution of Civil Engineers – Geotechnical Engineering
176
(3)
:
284
294
, .
Mittal
HK
(
1975
)
Effect of fines on permeability and consolidation characteristics of sandy soils
.
Indian Geotechnical Journal
5
(2)
:
135
150
.
Mosquera
J
(
2013
)
Static and Pseudo-Static Stability Analysis of Tailings Storage Facilities Using Deterministic and Probabilistic Methods
.
Department of Mining and Materials Engineering McGill University
,
Montreal
.
Nelson
J
,
Shepherd
T
and
Charlie
W
(
1977
)
Parameters Affecting Stability of Tailings Dams
. Proceedings of Geotechnical Practice for Disposal of Solid Waste Materials,
ASCE, University of Michigan
, pp.
444
460
.
Nordquist
TJ
(
2015
)
Permeability Anisotropy of Resedimented Mudrocks
.
Thesis of Master of Science in Civil and Environmental Engineering
.
Massachusetts Institute Of Technology
.
Oliva-González
AO
,
Butlanska
J
,
Fernández-Merodo
JA
and
Rodriguez-Pacheco
R
(
2024
)
Flow Failure of “La Luciana”
,
Spain
.
Orozco
R
(
2010
)
La ingeniería geotécnica en las presas de jales mexicanas
,
Memorias de la XXV Reunión Nacional de Mecánica de Suelos e Ingeniería Geotécnica
.
Acapulco Guerrero, México
.
Perlea
V
and
Botea
E
(
1984
)
Stability problems of tailings dams
, in
Proceedings 11th International. Conference of the ISSMFE
,
San Francisco
, pp.
1275
1280
.
Pettibone
HC
and
Kealy
CD
(
1971
)
Permeability and consolidation characteristics of sand – fines mixtures
.
Journal of the Soil Mechanics and Foundations Division, ASCE
97
(SM8)
:
1111
1130
.
Qiu
YY
and
Sego
DC
(
2001
)
Laboratory properties of mine tailings
.
Canadian Geotechnical Journal
38
(1)
:
183
190
, .
Saad
B
(
2008
)
Transient coupled analysis of upstream tailings disposal facilites construction
.
Thesis
.
Faculty of Graduate Studies and Research
.
Scholes
ON
,
Clayton
SA
,
Hoadley
AFA
and
Tiu
C
(
2007
)
Permeability anisotropy due to consolidation of compressible porous media
.
Transport in Porous Media
68
(3)
:
365
387
, .
Singh
S
,
Zurakowski
Z
,
Dai
S
and
Zhang
Y
(
2021
)
Effect of grain crushing on the hydraulic conductivity of tailings sand
.
Journal of Geotechnical and Geoenvironmental Engineering
147
(12)
, .
Terzaghi
K
,
Peck
RB
and
Mesri
G
(
1996
)
Soil Mechanics in Engineering Practice
.
John Wiley & Sons
.
Valenzuela
L
(
2015
)
Design, construction, operation and the effect of fines content and permeability on the seismic performance of tailings sand dams in Chile
.
Obras y Proyectos
(19)
:
6
22
, .
Vermeulen
NJ
(
2001
)
The composition and state of gold tailings
. A thesis submitted in partial fulfilment of the requirements for the degree Philosophiae Doctor (Engineering).
Vick
SG
(
1990
)
Planning, Design and Analysis of Tailings Dams
.
BiTech Publishers Ltd
, .
Villavicencio
G
,
Espinace
R
,
Palma
J
,
Fourie
A
and
Valenzuela
P
(
2014
)
Failures of sand tailings dams in a highly seismic country
.
Canadian Geotechnical Journal
51
(4)
:
449
464
.
Wang
D
,
Randolph
MF
and
Gourvenec
S
(
2015
)
Coefficient of consolidation for soil-that elusive quantity. In COUPLED VI
: Proceedings of the VI International Conference on Computational Methods for Coupled Problems in Science and Engineering (pp.
1218
1231
).
Cimne
.
Witt
KJ
and
Brauns
J
(
1983
)
Permeability anisotropy due to particle shape
.
Journal of Geotechnical Engineering
109
(9)
:
1181
1187
, .
Witt
KJ
and
Schönhardt
M
(
2004
)
Sustainable Improvement in Safety of Tailings Facilities
.
United Nations Economic Commission For Europe
.
Xu
C
(
2019
)
Long-Term Seepage Assessment Using Numerical Modeling for Upstream-Type Tailings Dams
,
Tesis for the Degree of Master of Applied Science
.
Concordia University
,
Montreal, Quebec, Canada
.
Xu
BH
,
He
N
,
Jiang
YB
,
Zhou
YZ
and
Zhan
XJ
(
2020
)
Experimental study on the clogging effect of dredged fill surrounding the PVD under vacuum preloading
.
Geotextiles and Geomembranes
48
(5)
:
614
624
, .
Zhang
C
,
Ma
C
,
Chen
Q
et al.
(
2020
)
Influence of rock percentage on strength and permeability of tailing-waste rock mixtures
.
Bulletin of Engineering Geology and the Environment
80
(1)
:
399
411
, .
Published by Emerald Publishing Limited. This article is published under the Creative Commons Attribution (CC BY 4.0) licence. Anyone may reproduce, distribute, translate and create derivative works of this article (for both commercial and non-commercial purposes), subject to full attribution to the original publication and authors. The full terms of this licence may be seen at Link to the terms of the CC BY 4.0 licenceLink to the terms of the CC BY 4.0 licence.

or Create an Account

Close Modal
Close Modal