Slurry Preparation Effects on the Cemented Phosphogypsum Backfill through an Orthogonal Experiment

Abstract

The cemented phosphogypsum (PG) backfill technique provides a new method for massive consumption of PG, and therefore alleviating the environmental pollution of PG. This study considered the effects of slurry preparation on the performance of cemented PG backfill. A L₁₆(4⁴) orthogonal experiment was designed to analyze four factors, namely the solid content, phosphogypsum-to-binder ratio (PG/B ratio), stirring time and stirring speed, with each factor having four levels. According to the range analysis, the solid content played the dominant role in controlling the bleeding rate, while the setting times strongly depended on the PG/B ratio. In terms of strength development of the backfill, the PG/B ratio was shown to be the most significant factor determining the unconfined compressive strength (UCS), followed by the solid content, stirring time and stirring speed. Furthermore, the results showed that the slurry preparation affected the environmental behavior of impurities that originated in PG. By analyzing the concentrations of impurities in the bleeding water of the slurry as well as the leachates of the tank leaching test, the results showed that the release of F⁻ and SO₄²⁻ was aggravated clearly with the increase in the PG/B ratio, while the release of PO₄³⁻ always remained at relatively low levels.

1. Introduction

In the phosphate industry, phosphoric acid is usually extracted from phosphate ore by using concentrated sulfuric acid, leaving a by-product of phosphogypsum (PG) which is mainly composed of CaSO₄·2H₂O. Approximately 4~5 t of PG is generated when 1 t of phosphoric acid is produced, resulting in the global production of PG being about 100 to 280 Mt every year [1,2,3]. Besides CaSO₄·2H₂O (>90%), PG also contains some impurities such as phosphate (PO₄³⁻), fluoride (F⁻), organic matters, heavy metals and radioactive components [4,5,6,7]. These minor compounds have not only caused serious environmental pollution at storage sites [8,9], but also hindered the reuse of PG, and thus only 15% of PG is recycled worldwide [10,11]. The local government of Guizhou, China recently launched a new policy stipulating that no new surface land would be allocated for PG storage. Therefore, it is extremely urgent to find an effective way to consume such a large amount of phosphogypsum. Alternatively, it has been estimated that about 60% of PG generated could be consumed by using PG as the aggregates in the backfill process [12]. In the cemented PG backfill process, the PG together with the hydraulic binder, water and some additives were mixed homogeneously to form a slurry on the land surface, and then the prepared slurry was transported to the underground stopes by gravity or pumping. In the stopes, the strength was developed by the hydration reaction of the binder, ensuring the stability of underground mine stopes and realizing the maximum recovery of ores. In addition, hardened backfill can well stabilize/solidify the minor compounds in PG, which is also a significant advantage for environment protection.

In the cemented PG backfill process, the first step is to prepare a qualified backfill slurry which could develop a desired strength after the placement. The effects of slurry preparation on the strength of cemented backfill have been investigated by some studies. For instance, Fall and Benzaazoua [13] demonstrated that an increase in the binder dosage could effectively improve the strength of cemented backfill samples. Ercikdi et al. [14] found that the solid content used for backfill process should be determined by the balance between the maximization of the strength and the effectiveness of pumping. Cao et al. studied the relations between strength performance of cemented tailings backfill and solid content, cemented-to-tailings ratio and curing time [15]. Furthermore, the stirring is also a crucial factor in the backfill preparation. Inhomogeneous stirring would lead to inconsistent distribution of solid and water, and thus reduce the strength of the backfill [16]. Although research on cemented PG backfill has been carried out for several years, the majority of the studies considered only a single factor for the optimization of backfill properties. Indeed, the slurry preparation contains multi-factors (such as solid content, aggregates-to-binder ratio and stirring), so it is important to determine which one is the most significant factor determining the slurry proprieties and the strength of the hardened backfill.

Furthermore, previous studies focused mainly on fluidity and/or strength of the backfill process, while relatively few studies paid attention to the environmental pollution related to the cemented backfill process. The potential pollution of backfill mainly stems from two aspects. First, when the slurry is placed into the stopes, the excess water used for improving the slurry fluidity will be secreted. Second, the backfill would suffer from the percolation of groundwater after forming the hardened backfill. Therefore, impurities in the bleeding water and the leachates would transfer into the groundwater, causing environmental pollution. Li et al. [12] showed that impurities in PG could be well solidified by comparing the concentrations of phosphate, sulfate and fluoride before and after the cementation. However, the question of whether the slurry preparation could affect the environmental behavior of cemented PG backfill has not been studied yet. In this case, this study will consider the influence of different preparing conditions on the bleeding water and the leachates in the backfill process.

An orthogonal experiment is usually designed to study multi-factor conditions for obtaining effective degree of each single factor on the results [17,18,19,20]. Therefore, the current study aims to optimize the slurry preparation by considering both the properties and environmental behaviors of cemented PG backfill. After discussions with the operators in mines, four factors related to the slurry preparation are considered in this study, namely phosphogypsum/binder (PG/B) ratio, solid content, stirring speed and stirring time. According to the orthogonal experimental design, the cemented PG backfill samples were reconstituted in the laboratory. The slurry properties were measured, including bleeding rate and setting times. The strength of hardened backfill was analyzed by measuring the unconfined compressive strength (UCS). The concentrations of phosphate (PO₄³⁻), fluoride (F⁻) and sulfate (SO₄²⁻) in the bleeding water and leachates of tank leaching test (TLT) were also determined. Then, the range analysis was applied to calculate the effective degree of each single factor, and therefore figuring out the effects of preparation conditions of the process of cemented PG backfill.

2. Materials and Methods

2.1. Raw Materials

Phosphogypsum and a composite binder (the mix proportion of binder is yellow phosphorous slag: fly ash: cement clinker = 4:1:1, and 16–20% lime of the yellow phosphorus slag mass ratio is added) [21] were collected from Guizhou Kailin (Group) Co., Ltd., Guiyang, China. The particle sizes of the PG and binder were measured by a Malvern Mastersizer 2000 particle size analyzer, as shown in Figure 1 and

Table 1. Chemical compositions and physical characteristics of phosphogypsum and binder.

Chemical Composition	PG	Binder
	%	%
Na2O	0.22	0.39
MgO	0.41	1.70
Al2O3	0.47	4.97
SiO2	1.76	23.31
P2O5	2.61	1.58
SO3	55.00	5.40
K2O	0.05	0.97
CaO	37.69	51.44
TiO2	0.03	0.46
Fe2O3	0.25	2.60
Physical Characteristic
D10 (μm)	19.40	6.08
D30 (μm)	54.78	14.18
D60 (μm)	102.41	33.66
Cu = D60/D10	5.28	5.54
Cc = D302/(D60 × D10)	1.51	0.98

. The main chemical compositions and physical characteristics of PG and binder used in this study are shown in Table 1. Coefficient of uniformity (Cu) and coefficient of curvature (Cc) are used to reflect the distribution of particle size.

2.2. Orthogonal Experiment

In this study, the cemented PG backfill slurry was prepared based on an orthogonal array [L₁₆(4⁴) matrix]. The following four factors were studied: solid content (factor A), phosphogypsum-to-binder ratio (PG/B ratio, factor B), stirring time (factor C) and stirring speed (factor D). The solid content refers to the mass percentage of PG and binder in slurry. A total of 16 formulations of cemented PG backfill were prepared in this study. The designed levels and factors are listed in

Table 2. Factors and levels in the orthogonal experiment.

Factors		Levels
		1	2	3	4
A	Solid content (%)	45	50	55	60
B	PG/B ratio	2:1	3:1	4:1	5:1
C	Stirring time (min)	5	30	60	120
D	Stirring speed (rpm)	300	400	500	600

.

In order to figure out which factor is the most significant one, range analysis is essential. Two parameters k_ij and R_j were used for evaluation. k_ij is defined as the sum of the evaluation indexes of all levels (j, j = 1, 2, 3, 4) in each factor (i, i = A, B, C, D), and ${\overline{\mathrm{k}}}_{\mathrm{i}\mathrm{j}}$ (mean value of k_ij) is used to determine the optimal level and combination of factors. When ${\overline{\mathrm{k}}}_{\mathrm{i}\mathrm{j}}$ is the largest, it can be considered that this level is most optimal for this factor. R_j is defined as the range between the maximum and minimum value of ${\overline{\mathrm{k}}}_{\mathrm{i}\mathrm{j}}$, and it is used for evaluating the importance of the factors to each evaluation index. The larger R_j value, the greater importance of this factor [22,23]. M is described as the evaluation index. For this L₁₆(4⁴) matrix, the relevant calculations are as follows (factor B, for example):

$$\begin{array}{l}{\mathrm{k}}_{\mathrm{B}1}={\mathrm{M}}_1+{\mathrm{M}}_5+{\mathrm{M}}_9+{\mathrm{M}}_{13};\\ {\mathrm{k}}_{\mathrm{B}2}={\mathrm{M}}_2+{\mathrm{M}}_6+{\mathrm{M}}_{10}+{\mathrm{M}}_{14};\\ {\mathrm{k}}_{\mathrm{B}3}={\mathrm{M}}_3+{\mathrm{M}}_7+{\mathrm{M}}_{11}+{\mathrm{M}}_{15};\\ {\mathrm{k}}_{\mathrm{B}4}={\mathrm{M}}_4+{\mathrm{M}}_8+{\mathrm{M}}_{12}+{\mathrm{M}}_{16};\\ \overline{{\mathrm{k}}_{\mathrm{B}1}}=\frac{{\mathrm{k}}_{\mathrm{B}1}}{4}\\ \overline{{\mathrm{k}}_{\mathrm{B}2}}=\frac{{\mathrm{k}}_{\mathrm{B}2}}{4}\\ \overline{{\mathrm{k}}_{\mathrm{B}3}}=\frac{{\mathrm{k}}_{\mathrm{B}3}}{4}\\ \overline{{\mathrm{k}}_{\mathrm{B}4}}=\frac{{\mathrm{k}}_{\mathrm{B}4}}{4}\\ {\mathrm{R}}_{\mathrm{B}}=\max \left(\overline{{\mathrm{k}}_{\mathrm{Bj}}}\right)-\min \left(\overline{{\mathrm{k}}_{\mathrm{Bj}}}\right)\end{array}$$

Then, the range analyses applied the above method, and M is replaced by the index that needs to be evaluated, such as the bleeding rate, setting times, strength and concentrations of impurities.

2.3. Sample Preparation

The fresh slurry was prepared by mixing PG, binder and deionized water. After mixing for certain times according to the orthogonal design, the slurry was collected for the analysis of bleeding rate, initial setting time and final setting time.

Then, the prepared slurry was poured into a plastic mold with internal dimensions of 40 mm × 40 mm × 40 mm. There was a 2 mm hole at the bottom of each mold to allow the drainage of excess water. The bleeding water draining from the mold was collected and filtered through a 0.45 μm filter for further analysis to evaluate the environmental behavior. After the initial solidification, the samples were demolded and cured in a chamber at a constant temperature of 20 ± 2 °C and a humidity of 90 ± 5%.

2.4. Tank Leaching Test

The tank leaching test (TLT) is commonly used to study the release properties of minor compounds and to predict the number of released ions in the hardened backfill [24,25,26]. The potential hazards related to the leachable impurities from hardened backfill samples were measured by conducting compliance leaching test according to EA NEN 7375. The backfill samples cured for 28 days were placed in 500 mL plastic bottles and immersed in the bottle with plastic thin threads. The deionized water was used as leachate, and the liquid/solid ratio was 5 cm³ of solution per cm² of exposed solid. The leachate was changed with the same volume of deionized water after a cumulative leaching time of 0.25, 1, 2.25, 4, 9, 16, 36 and 64 days. Then the leachate at each period was collected and filtered (0.45 μm) for the subsequent analysis. The measured leaching per fraction was calculated separately for each factor using the Equation (1), and the cumulative leaching quantity was calculated by Equation (2), as specified in the EA NEN 7375:2004 standard:

$$E_i=\frac{c_i\times V}{f\times A}$$

(1)

$$T_n=\frac{{\displaystyle \sum _{i=1}^n{E}_i•A}}{f}\begin{array}{cc}& n=1\sim N\end{array}$$

(2)

where $E_i$ is the measured leaching of a component in fraction i in mg/m²; $C_i$ is the concentration of the component in fraction i in ug/L; V is the volume of the leachate in l; A is the surface area of the test sample in m²; f is a conversion factor in 1000 μg/mg; T_n is the cumulative leaching quantity of component in fraction i in g.

2.5. Analytical Methods

2.5.1. Bleeding Rate and Setting Times

The bleeding rate was measured in accordance with the Chinese standard GB/T 50080-2016. The slurry was poured into a container with a lid, which was then shaken for 20 s on a shaking table. The bleeding water was drawn with a syringe at intervals until no more water was secreted. The initial slurry and slurry after bleeding in containers were weighed, and then the bleeding rate was calculated.

The initial setting time (IST) and final setting time (FST) were determined according to the Chinese standard GB/T1346-2001 using a Vicat apparatus. The slurry was poured into the Vicat round mold, and bleeding water was secreted over a period of time. The IST was considered as the time when initial Vicat needle penetrated the sample to 5 ± 1 mm from the bottom of the mold. The FST was defined as no visible mark being left when the final Vicat needle was on the surface of the sample [27]. Tests were carried out in triplicate and the reported values were the averages of these three tests.

2.5.2. Unconfined Compressive Strength of Cemented PG Backfill Samples

The minimum strength of backfill into the underground must meet the strength requirements in order to ensure the stability of underground mine stopes. Therefore, the strength development is one of the most important mechanical properties in the backfill system design [28]. According to the Chinese Standard JGJ/T 70-2009 [29], cemented backfill cured for 28 days was used for testing the UCS. UCS tests were conducted using a servo-hydraulic machine of 200 kN loading capability at a constant displacement rate of 0.1 mm/min. In order to avoid the randomness and contingency of the test data, three samples were prepared for each test, and the average values of UCS were calculated.

2.5.3. Microstructural Analysis

After UCS tests, the broken samples were immediately soaked in an ethanol solution to stop the hydration reaction. The samples were then dried at 50 °C in an oven. SEM analysis with the HELIOS NamoLab 600i (FEI, Lake Oswego, OR, USA). was used to find the micromorphologic development inside the cemented backfill samples [30]. Owing to the poor conducting performance of backfill, it was necessary to cover a layer of gold coating on samples to enhance the conductive properties of the samples in order to meet the inspection requirements.

2.5.4. Chemical Measurements

The pH of bleeding waters and leachates after TLT tests were measured by Ohaus TSARTR 300 pH meter (Ohaus, Parsippany, NJ, USA). The concentrations of F⁻ were measured with the Leici PF-1-01 fluorine ion-selective electrode (Leici, Shanghai, China). The concentrations of SO₄²⁻ and PO₄³⁻ in bleeding waters and leachates of TLT tests were measured using baryta yellow spectrophotometric method and ammonium molybdate tetrahydrate spectrophotometric method, respectively, via a Shimadzu UV1800 spectrophotometric instrument (Shimadzu, Kyoto, Japan).

3. Results and Discussion

3.1. Differences of the Slurry Properties and Strength of Cemented PG Backfill

According to the L₁₆(4⁴) matrix, 16 batches were carried out, and the results are shown in

Table 3. Test data of evaluation indexes under different experimental conditions.

Batch No.	Factors				Evaluation Indexes
	A Solid Content	B PG/B Ratio	C Stirring Time	D Stirring Speed	Backfill		Slurry Properties			Bleeding Waters			Leachates of TLT
					UCS (MPa)	Wet Unit Weight (kN/m3)	IST (h)	FST (h)	Bleeding Rate (%)	F− (mg/L)	SO42− (mg/L)	PO43− (mg/L)	pH	F− (mg)	SO42− (mg)	PO43− (mg)
1	45%	2:1	5	300	0.83	13.22	87	105	40.10	257	2149	2.55	10.62	28	876	0.25
2	45%	3:1	30	400	0.81	13.42	98	112	35.17	267	2433	2.50	10.75	38	1482	0.26
3	45%	4:1	60	500	0.74	12.89	103	122	46.40	282	2300	5.00	10.46	43	1597	0.28
4	45%	5:1	120	600	0.77	12.65	99	124	42.54	454	2888	1.80	10.13	50	1751	0.32
5	50%	2:1	30	500	1.57	13.26	80	95	38.70	277	2168	2.05	11.18	27	1080	0.33
6	50%	3:1	5	600	0.92	13.44	99	115	35.73	287	2266	2.25	11.26	37	1601	0.37
7	50%	4:1	120	300	0.87	13.02	98	106	43.78	356	2518	3.15	10.71	39	1341	0.36
8	50%	5:1	60	400	0.80	12.58	92	108	42.29	533	2870	1.75	10.08	60	1848	0.41
9	55%	2:1	60	600	2.06	14.22	76	88	37.49	305	2489	3.10	11.40	25	939	0.34
10	55%	3:1	120	500	1.17	13.78	88	96	40.42	307	2168	2.85	11.02	44	1651	0.38
11	55%	4:1	5	400	0.96	13.04	92	112	32.29	376	2540	2.35	10.51	47	1732	0.40
12	55%	5:1	30	300	0.93	12.54	97	125	32.25	545	3167	2.75	10.35	48	1890	0.40
13	60%	2:1	120	400	2.26	14.22	72	85	27.60	306	2395	2.75	11.52	27	1060	0.37
14	60%	3:1	60	300	1.41	13.86	85	100	26.34	335	2462	2.40	11.09	41	1812	0.39
15	60%	4:1	30	600	1.37	13.22	83	94	30.64	415	2529	2.65	10.71	49	2036	0.40
16	60%	5:1	5	500	1.02	12.98	98	116	26.45	615	3483	4.10	10.28	51	2054	0.41

. It can be seen that the bleeding rate ranged from 26.34~46.40%. The initial setting time (IST) and final setting time (FST) were in the ranges from 72 h to103 h and 85 h to 125 h, respectively. The unconfined compressive strength (UCS) of samples cured for 28 d ranged from 0.74 MPa to 2.26 MPa. These results indicate that the slurry preparation conditions have significant effects on the cemented PG backfill, including both properties of slurry and the strength of hardened backfill. Different mean values indexes ${\overline{\mathrm{k}}}_{\mathrm{i}\mathrm{j}}$ at different levels of each factor are listed in

Table 4. Range analysis data of mechanical and physical properties.

Evaluation Index	Value Name	A Solid Content (%)	B PG/B Ratio	C Stirring Time (h)	D Stirring Speed (rpm)
UCS (MPa)	k1	0.79	1.68	0.93	1.01
	k2	1.04	1.08	1.17	1.21
	k3	1.28	0.99	1.25	1.13
	k4	1.52	0.88	1.27	1.28
	Rj	0.73	0.80	0.33	0.27
Ranking	B > A > C > D
IST (h)	k1	97	79	94	92
	k2	92	93	90	89
	k3	88	94	89	92
	k4	85	97	89	89
	Rj	12	18	5	3
Ranking	B > A > C > D
FST (h)	k1	116	93	112	109
	k2	106	106	107	104
	k3	105	109	104	107
	k4	99	118	103	105
	Rj	17	25	9	5
Ranking	B > A > C > D
Bleeding rate (%)	k1	41.05	35.97	33.64	35.62
	k2	40.13	34.41	34.19	34.34
	k3	35.61	38.28	38.13	37.99
	k4	27.76	35.88	38.59	36.60
	Rj	13.29	3.87	4.95	3.65
Ranking	A > C > B > D

. By comparing R_j of different factors of IST, FST, and UCS, the order of factor significance is as follows: the PG/B ratio > solid content > stirring time > stirring speed. However, the order of significant factors for bleeding rate is the solid content> stirring time>PG/B ratio > stirring speed.

3.1.1. Effects of Slurry Preparation on Slurry Properties

The setting times (IST and FST),which are related to the time periods for slurry transportation, are important slurry properties to evaluate the performance of backfill [31]. As shown in Figure 2a,b, PG/B ratio was the most significant factor, and the IST and FST increased with the increase in PG/B ratio. It is well known that residual acid and other minor compounds in PG could retard the cementation process [32], leading to the extension of IST and FST with the high PG/B ratio. Furthermore, the IST and FST both reduced evidently with the increase in solid content. The increase in solid content meant a relatively higher binder to water ratio in the slurry, which would accelerate the hydration reactions, resulting in a reduction in IST and FST. In previous studies using cement as the binder, the setting times of the slurry were usually less than 1 day [27,33]. However, the setting times in this study ranged from 2 days to 5 days, which could be attributed to the retarding effect of PG. Some studies have found that sulfate significantly affected the setting time of cement [34]. As such, the sulfate in phosphogypsum might be one of the reasons for the increase of setting times. In this case, the smaller the IST and FST of cemented PG backfill, the better the filling efficiency. Hence, according to the above analysis, the optimal combination of slurry preparation conditions for shortening IST and FST is respectively, a solid content of 60% and a PG/B ratio of 2:1.

As shown in Figure 2c, the bleeding rate is also affected by the slurry preparation. The bleeding rate decreased sharply with the increase in solid content, illustrating that the solid content was the key to determining how much water would be secreted into stopes after the placement of the slurry. With the increase in the solid content, more water would be consumed for the hydration reactions, and would thus lower the bleeding rate. Compared to the solid content as the most significant factor, variations of bleeding rate related to other factors (PG/B ratio, stirring time and stirring speed) were not significant. In the actual filling process, the excess water usually needs to be discharged from the mine by pumping machine, which would increase the cost. Meanwhile, impurities in the bleeding water would transfer into the groundwater. Therefore, considering from these aspects, a relatively low bleeding rate is recommended.

3.1.2. Effects of Slurry Preparation on Unconfined Compressive Strength

The strength of the backfill is one of the most important parameters that directly affect the safety performance of backfill [35]. According to the range analysis as shown in Table 4, four factors were all affecting the strength of cemented PG backfill samples. For the PG/B ratio, the UCS decreased from 1.68 MPa to 0.88 MPa (decreased by 48%) with the PG/B ratio increasing from 2:1 to 5:1. However, for the solid content increasing from 45% to 60%, the UCS increased from 0.79 MPa to 1.52 MPa (increased by 92%). For the stirring time increasing from 5 min to 120 min, the UCS increased from 0.93 MPa to 1.27 MPa (increased by 37%). For the stirring speed, the UCS increased from 1.01 MPa to 1.28 MPa (increased by 27%) with the stirring speed increasing from 300 rpm to 600 rpm.

It is clear that the PG/B ratio is the most significant factor affecting the strength of the backfill, as indicated by the R_j in Table 4 and the trend charts in Figure 2d. Previous studies showed that the strength development of the backfill was mainly ascribed to the formed cementitious products (such as C-S-H gel and ettringite) by hydration reactions [36,37]. A high PG/B ratio meant a low binder proportion in the mixture, resulting in less hydration products and therefore low strength. In this study, two backfill samples with the highest PG/B ratio (13# at a ratio of 2:1) and the lowest PG/B ratio (16# at a ratio of 5:1) were taken for SEM-EDS analysis. As shown in Figure 3a,b, the SEM images showed that the ettringite and C-S-H gel in sample 13#were much more than those in sample 16#, corresponding to the UCS of 2.26 MPa and 1.02 MPa, respectively. EDS analyses were done to identify CaSO₄, ettringite and C-S-H gel in the backfill samples, as shown in Figure 3c–e. According to the EDS analysis, in sample 13# with PG/B ratio of 2:1, great amount of hydration products (C-S-H gel and ettringite) were wrapping around the PG, providing strength development for cemented PG backfill. In addition, it is notable that the PG/B ratio could also affect the physical properties of the backfill samples, which in turn influenced its mechanical properties. In this study, the wet unit weights of the backfill samples were measured, as shown in Table 3. Orthogonal range analysis and the trend charts of wet unit weight (provided in the Annex File) show that the PG/B ratio has the most significant effect on the wet unit. Kitazume et al. [38] has demonstrated that strength ratio increased with the wet unit weight ratio in soil and Portland cement mixture. Therefore, it seems that the PG/B ratio influences the mechanical performance of the backfill samples by affecting both the physical properties and the cementitious products. More binder could effectively improve the strength, which could benefit the stability of underground mine stopes. However, it is suggested that the binder accounts for over 75% of the total cost of backfill process [31], so high binder dosage might burden the backfill process [39]. Thus, considering both cost-effectiveness and strength requirements, a reasonable PG/B ratio must be decided according to the actual conditions of mines.

3.2. Impurities in Bleeding Water

The concentrations of impurities of the bleeding water were measured, as listed in Table 3. The concentrations of F⁻, SO₄²⁻ and PO₄³⁻ were in the ranges of 257–615 mg/L, 2149–3483 mg/L and 1.75–5.00 mg/L, respectively. Compared to the impurities in PG leachate (the concentrations of F⁻, SO₄²⁻ and PO₄³⁻ of PG were1355 mg/L, 23,842 mg/L and 315 mg/L, respectively), the majority of these ions were solidified/stabilized when the PG was cemented in the backfill. The average solidification rates were 73% for F⁻, 89% for SO₄²⁻ and 99% for PO₄³⁻. These results indicate that cementation was an effective way to control PG pollution.

To understand the degrees of effects of slurry preparation on the quality of bleeding water, ${\overline{\mathrm{k}}}_{\mathrm{i}\mathrm{j}}$ and R_j of four factors are listed in

Table 5. Range analysis data of conductivity and impurities concentration in bleeding water.

Evaluation Index	Value Name	A	B	C	D
F− (mg/L)	k1	315	286	384	373
	k2	363	299	376	371
	k3	383	357	364	370
	k4	418	537	356	365
	Rj	103	251	28	8
Ranking	B > A > C > D
SO42− (mg/L)	k1	2443	2300	2610	2574
	k2	2456	2332	2574	2559
	k3	2591	2472	2530	2530
	k4	2717	3102	2492	2543
	Rj	274	802	118	44
Ranking	B > A > C > D
PO43− (mg/L)	k1	2.96	2.61	2.81	2.71
	k2	2.30	2.50	2.49	2.34
	k3	2.76	3.29	3.06	3.50
	k4	2.98	2.60	2.64	2.45
	Rj	0.68	0.79	0.58	1.16
Ranking	D > B > A > C

. As shown in Table 5 and Figure 4a,b, the concentrations of SO₄²⁻ and F⁻ in the bleeding water rose clearly with the increase in solid content and the PG/B ratio, but dropped slightly with the increase in stirring time and stirring speed. To be specific, for solid content, higher solid content meant more solids in the slurry system, leading to more dissolved ions in the bleeding water. The PG/B ratio was the most significant factor influencing the quality of the bleeding water, as shown in Figure 4a,b. The higher PG/B ratio, the higher concentrations of F⁻ and SO₄²⁻ that were detected in the bleeding water. The reason for this might be that the F⁻ and SO₄²⁻ were mainly sourced from PG, so higher PG dosage led to higher concentrations of F⁻ and SO₄²⁻ in the bleeding water. For stirring time, the impurities concentrations decreased within a certain small range with the increase in stirring time. This can be explained by the fact that the longer the stirring time, the longer the reaction time for hydration, resulting in more ions consolidated in the backfill and thus less free ions dissolved in the bleeding water [40]. Similarly, with the increase in the stirring speed of slurry, the contact frequency between the impurities and cementitious materials increased, leading to the enhancement of the cementation reactions, which would reduce the concentrations of F⁻ and SO₄²⁻ in the bleeding water.

As opposed to the SO₄²⁻ and F⁻ varying with slurry preparation conditions, the concentrations of PO₄³⁻ always kept at very low levels at all experimental batches, which was likely due to the different solidification mechanisms of these ions. It is reported in previous studies that the PO₄³⁻ reacts rapidly with the large amounts of calcium in the alkaline environment, leading to the precipitation of dissolved PO₄³⁻ and thus a low concentration of PO₄³⁻ in the bleeding water [41,42]. At the same time, the EDS analysis revealed that the peaks of Ca, O and P were found, inferring the precipitation of calcium phosphate. Although cemented PG backfill technique has a strong capacity to consolidate impurities in backfill, there is still a certain amount of F⁻ and SO₄²⁻, which might transfer into the groundwater. Therefore, in bleeding water, the concentrations of impurities should be kept as low as possible. Thus, the slurry preparation conditions with 45% of the solid content, 2:1 of PG/B ratio, 120 min of stirring time and 600 rpm of stirring speed was thought to be the optimal combination by considering impurities in the bleeding water.

3.3. Impurities in the Leachates of Tank Leaching Test (TLT)

After the initial bleeding water secretion, the backfill would harden and suffer from the underground water passing through, so it is significant to understand the leaching behavior of impurities in the long-time period. The tank leaching test (TLT) is commonly used to study the dynamic and static properties of impurities in the cement materials [43]. Therefore, in this study, the TLT test with eight leaching periods was adopted to study if the slurry preparation could affect the impurities performance in the long-term water immersion. The concentrations of impurities after TLT were measured as listed in Table 3.

3.3.1. pH Variation

pH values during total eight leaching periods are listed in

Table 6. pH values in eight leaching periods.

Trial No.	pH
	Period 1 (0.25 d)	Period 2 (1 d)	Period 3 (2.25 d)	Period 4 (4 d)	Period 5 (9 d)	Period 6 (16 d)	Period 7 (36 d)	Period 8 (64 d)
1	8.81	9.87	10.28	10.19	10.80	11.60	11.27	10.62
2	9.20	10.17	10.27	10.61	10.89	11.74	11.40	10.75
3	9.44	9.97	10.57	10.42	10.63	11.43	11.12	10.46
4	9.49	9.44	9.88	9.87	10.48	11.08	10.75	10.13
5	10.47	10.06	10.54	10.45	10.85	11.82	11.64	11.18
6	8.98	9.27	10.63	10.60	10.64	11.64	11.59	11.26
7	8.77	8.91	9.95	8.82	10.29	11.27	10.92	10.71
8	8.65	7.84	9.05	8.61	9.58	10.17	10.18	10.08
9	9.27	9.63	10.33	10.03	10.68	11.68	11.60	11.40
10	9.02	9.93	10.19	10.22	10.64	11.70	11.42	11.02
11	8.87	8.66	9.21	9.26	10.33	11.27	10.86	10.51
12	8.40	9.53	9.06	8.71	10.03	10.76	10.46	10.35
13	9.88	9.32	9.83	10.37	10.60	11.75	11.47	11.52
14	8.08	9.46	9.55	9.28	10.50	11.37	11.13	11.09
15	8.01	8.29	9.25	9.28	10.37	11.01	10.87	10.71
16	8.94	8.48	9.10	8.85	9.77	10.42	10.34	10.28

. It can be seen that with the replacement of leachate, pH increased from the first leaching period to the sixth leaching period, which was likely due to the constant release of the hydroxide ions generated in the hydration process. The pH values reached the peak at the sixth leaching on 16 d, and then decreased from the sixth leaching period to the last leaching period. This indicates the finish of the hydration process on 16 d, and less hydroxide ions could be released into the leachates [44].

To investigate if the slurry preparation conditions could affect the pH of the leachates in the long run, the pH values at the last leaching period were used to approximate the pH of the underground water after long-term immersion. Different evaluation indexes of four factors are shown in

Table 7. Range analysis of pH and cumulative leaching quantity of TLT.

Evaluation Index	Value Name	A	B	C	D
pH	k1	10.49	11.18	10.67	10.69
	k2	10.81	11.03	10.75	10.72
	k3	10.82	10.60	10.76	10.74
	k4	10.90	10.21	10.85	10.88
	Rj	0.41	0.97	0.18	0.18
Ranking	B > A > C > D
F− (mg)	k1	39	27	41	39
	k2	41	40	40	43
	k3	41	44	42	41
	k4	42	52	40	40
	Rj	3	25	2	4
Ranking	B > D > A > C
SO42− (mg)	k1	1427	989	1566	1480
	k2	1468	1637	1622	1530
	k3	1553	1677	1549	1596
	k4	1741	1885	1451	1582
	Rj	314	896	171	116
Ranking	B > A > C > D
PO43− (mg)	k1	0.28	0.32	0.36	0.35
	k2	0.37	0.35	0.35	0.36
	k3	0.38	0.36	0.36	0.34
	k4	0.39	0.39	0.36	0.36
	Rj	0.11	0.07	0.01	0.02
Ranking	A > B > D > C

and Figure 5a. The pH decreased sharply with the increase in PG/B ratio, and gradually increased with the increase in solid content, stirring time and stirring speed. Apparently, the PG/B ratio was the most significant factor influencing the pH of leachate. It is well known that the PG is usually acidic, while the binder is alkaline [45]. Therefore, a high PG/B ratio meant a high dosage of acid and a low dosage of alkali, leading to the low pH values of leachates. At the same time, the solid content is also a significant factor influencing the pH value of TLT. As shown in Figure 5a, the pH raised with the increase in solid content, which might be explained by more hydroxide ions being produced due to the higher solid content.

3.3.2. Cumulative Effects of Impurities on the Environment

In order to evaluate the long-term environmental behavior of impurities in the backfill samples, the cumulative leaching quantities of impurities in the PG and those in the cemented PG backfill (taking sample 13# as an example) was used for comparison, as shown in Figure 6. It is clear that the total quantities of impurities in leachates of backfill samples were much less than those in the leachates of PG. The cumulative leaching quantity of SO₄²⁻ in PG is 5 times as much as that in backfill, the cumulative leaching quantity of F⁻ in PG is 81 times in backfill, and the cumulative leaching quantity of PO₄³⁻ in PG is 1678 times in backfill. The reason for this is that, when PG formed a dense structure, the impurities were precipitated and/or incorporated in the hydration products, and thus impurities were less likely to escape from the cemented PG backfill [46].

3.3.3. Effects of Slurry Preparation on Leaching Behavior of Impurities

According to the tank leaching test, the ranges of cumulative leaching quantities of F⁻, SO₄²⁻ and PO₄³⁻ were 25~60 mg, 876~2054 mg and 0.25~0.41 mg, respectively. The variation trends of four factors are shown in Figure 5. It is clear that the PG/B ratio was the most influential factor, and the cumulative leaching quantities of F⁻ and SO₄²⁻ increased with the increase in PG/B ratio. It is likely due to that the low strength of backfill related to the high PG/B ratio could not well solidify the impurities in PG [47], resulting in the increase in the cumulative leaching quantities of F⁻ and SO₄²⁻. However, the cumulative quantity of PO₄³⁻ was much less compared to those of F⁻ and SO₄²⁻ in the leachates, as shown in Figure 5. The hydration reactions would provide an alkaline environment and a great amount of calcium ions for the precipitation of F⁻, SO₄²⁻ and PO₄³⁻ [12]. However, the solubility product constant of calcium phosphate (2.0 × 10⁻²⁹) was several orders of magnitude lower than those of calcium fluoride (5.3 × 10⁻⁹) and calcium sulfate (9.1 × 10⁻⁶), resulting in a lower leaching quantity of PO₄³⁻ in the leachate under long-term immersion (as shown in Figure 5d).

4. Conclusions

In this study, the orthogonal experiment was designed to understand the effect of slurry preparation conditions on the performance of cemented PG backfill. Four factors were examined: solid content, PG/B ratio, stirring time and stirring speed. The following conclusions can be made. Firstly, according to the range analysis, the most significant factor affecting the setting times and strength was the PG/B ratio, followed by the solid content. The solid content has the most significant effect on the bleeding rate. Considering both the slurry properties and strength development, it is determined the optimal condition is the slurry concentration of 60% and PG/B ratio of 2:1. Secondly, the lower the concentrations of impurities in bleeding water and TLT leachate, the friendlier the conditions are for the ground water. Considering environmental behaviors, it is determined that optimal combination is a solid concentration of 45%, PG/B ratio of 2:1, stirring time of 120 min and stirring speed of 600 rpm. Results showed that the differently optimal combinations exist for considering the slurry properties, mechanical strength and environmental behaviors of cemented PG backfill. Therefore, it is recommended that mines choose different optimization conditions according to their actual demands.