Main cation variation characteristics. Figure 1. shows the changes in the main cation content in oil shale-water reaction solutions at different reaction temperatures. The figure shows that under the reaction conditions of 20 ℃ and 50 ℃, the contents of the four cations in the aqueous solution have little difference. At reaction temperatures of 20°C and 50°C, the content of Ca2+ is the highest, followed by K+ and Na+, and the content of Mg2+ is the lowest. At reaction temperatures of 80°C and 100°C, the content of Na+ in the oil shale-water solution changes strongly; the Na+ content is the highest, followed by K+ and Ca2+, and the Mg2+ content is the lowest, possibly because Ca2+ and Mg2+ in aqueous solution come mainly from the dissolution of dolomite and calcite. With increasing temperature, the solubility of dolomite and calcite decreases, and the reaction equations are as follows:
Moreover, K+ changes little at 20°C, 50°C, and 80°C but increases rapidly at 100°C, indicating that a temperature of approximately 100°C has the greatest influence on the K+ content. The cation content in the aqueous solution tends to be stable at approximately 15 days.
As shown in Fig. 2, compared with the cation content in the oil shale-aqueous solution, the cation content in the oil shale ash-aqueous solution shows a different trend with increasing reaction temperature. Under different reaction temperatures, the K+ content in oil shale ash-aqueous solution always remains the highest, and the K+ content at 50 ℃, 80 ℃, and 100 ℃ is twice the K+ content at 20 ℃. At 20°C, the average content of Na+ (5.79 mg/L) was lower than the average content of Ca2+ (7.10 mg/L), but with increasing reaction temperature, the Na+ content gradually increased and remained higher than the Ca2+ content. Ca2+ has a declining trend with increasing temperature. Under a reaction temperature of 100°C, the Na+ content in the oil shale ash-aqueous solution increases rapidly with increasing reaction time. After 30 days of reaction, the Na+ content reached 16.10 mg/L, which was close to the K+ content, indicating that high temperature is favorable for the dissolution of Na and K minerals in oil shale ash.
The variation in the content of the main cations in the surrounding rock-aqueous solution with different reaction temperatures is shown in Fig. 3. With increasing reaction time, the content of cations showed an upwards trend. Under the reaction conditions of 20°C and 50°C, the content of Ca2+ in the aqueous solution was the highest, followed by Na+ and K+, and Mg2+ was the lowest. The content of Na+ is the highest at 80°C and 100°C, followed by Ca2+, K+, and Mg2+. With increasing reaction time, the cationic content in aqueous solution increased steadily under the reaction condition of 80°C but fluctuated greatly under the reaction condition of 100°C.
Main anion variation characteristics.The variation trend of the major anion content in the oil shale-aqueous solution under different temperature reaction conditions is shown in Fig. 4. The content of HCO3− is always the highest, followed by Cl−, and the content of HCO3− is more than 15 mg/L in the aqueous solutions at different reaction temperatures. When the reaction temperature was 20 ℃, 50 ℃ and 100 ℃, the content of Cl− increased significantly with increasing reaction time, especially when the reaction temperature was 100 ℃. After 30 days of reaction, the content of Cl− reached 15.4 mg/L. The content change of SO42− is similar to the content change of Cl−, and the change is small at reaction temperatures of 20°C, 50°C, and 80°C. The content of SO42− in aqueous solution is less than 3 mg/L, but the content of SO42− increases sharply at a reaction temperature of 100°C. The content of F− and NO3− in the aqueous solution are always low at different temperature gradients.
The variation in the main anion content in oil shale ash slag-aqueous solution with reaction temperature is shown in Fig. 5. HCO3− is the main anion in aqueous solution, the content of HCO3− is greater than the content of oil shale aqueous solution, and the content is more than 20 mg/L. The content of HCO3− decreases with increasing reaction time in aqueous solutions with higher reaction temperatures (80 ℃ and 100 ℃). Different from oil shale aqueous solution, the content of SO42− in aqueous solution increases significantly, and SO42− becomes the second-largest anion, surpassing the content of Cl−. With increasing reaction time, the content of SO42− also increases, especially after 30 days of reaction at 80°C, and the content of SO42− in aqueous solution reaches 36.97 mg/L. SO42− in aqueous solution generally comes from the dissolution of minerals containing gypsum or other sulfates. In addition, oil shale and its surrounding rock contain a large amount of pyrite. Oxidation of pyrite results in the presence of water-insoluble sulfur in water as SO42−. At 100°C water temperature, the growth rule of Cl− is similar to that of SO42−. The F− and NO3− content is still low, and the NO3− content is almost zero.
Figure 6. describes the changes in the main anion content in the surrounding rock-aqueous solution at different reaction temperatures. Under the reaction conditions of 20 ℃, 50 ℃, and 80 ℃, the variation trend of the main anions in the oil shale surrounding rock water solution is similar to the variation trend of the main anions in the oil shale water solution. HCO3− is the main anion, followed by Cl−, and the content of HCO3− and Cl− gradually increases with time. However, at the reaction temperature of 100°C, HCO3− and Cl− content changes differently from other temperature conditions, and both of these ions increase first and then decrease with the reaction time. The content of SO42−, F- and NO3− in the aqueous solution did not change significantly and remained at a low level.
Analysis of migration mechanism. The correlation coefficient can accurately describe the degree of correlation between variables in numerical form[30]. The correlation coefficient matrix can reflect the correlation between each research parameter and characterize the coexistence of the whole dataset or the relationship between any two indices[31]. The Pearson correlation coefficient matrix of each hydrogeochemical index was calculated by SPSS software, as shown in Table 1.
Table 1
The Pearson correlation coefficient matrix
|
K
|
Na
|
Mg
|
Ca
|
HCO3−
|
F
|
Cl
|
NO3
|
SO4
|
pH
|
conductivity
|
reaction time
|
reaction temperature
|
lithology
|
K
|
1.000
|
0.731**
|
-0.204*
|
0.290**
|
0.821**
|
0.029
|
0.451**
|
-0.272**
|
0.857**
|
0.355**
|
0.765**
|
0.078
|
0.310**
|
0.830**
|
Na
|
|
1.000
|
-0.341**
|
0.240*
|
0.674**
|
0.254**
|
0.796**
|
-0.406**
|
0.697**
|
0.159
|
0.859**
|
0.420**
|
0.576**
|
0.495**
|
Mg
|
|
|
1.000
|
0.472**
|
0.068
|
-0.248**
|
-0.218*
|
0.076
|
-0.129
|
-0.011
|
-0.165
|
0.129
|
-0.639**
|
-0.053
|
Ca
|
|
|
|
1.000
|
0.688**
|
-0.193*
|
0.149
|
-0.210*
|
0.239*
|
0.130
|
0.340**
|
0.204*
|
-0.186
|
0.269**
|
HCO3-
|
|
|
|
|
1.000
|
-0.046
|
0.342**
|
-0.359**
|
0.612**
|
0.329**
|
0.678**
|
0.131
|
0.186
|
0.648**
|
F
|
|
|
|
|
|
1.000
|
0.291**
|
0.087
|
-0.027
|
0.116
|
0.188
|
0.143
|
0.338**
|
-0.135
|
Cl
|
|
|
|
|
|
|
1.000
|
-0.334**
|
0.509**
|
-0.067
|
0.672**
|
0.506**
|
0.491**
|
0.273**
|
NO3
|
|
|
|
|
|
|
|
1.000
|
− .277**
|
0.126
|
− .368**
|
-0.103
|
-0.335**
|
0.007
|
SO4
|
|
|
|
|
|
|
|
|
1.000
|
0.200*
|
0.790**
|
0.206*
|
0.227*
|
0.758**
|
pH
|
|
|
|
|
|
|
|
|
|
1.000
|
0.240*
|
-0.079
|
0.017
|
0.285**
|
conductivity
|
|
|
|
|
|
|
|
|
|
|
1.000
|
0.287**
|
0.424**
|
0.601**
|
reaction time
|
|
|
|
|
|
|
|
|
|
|
|
1.000
|
0.000
|
0.000
|
reaction temperature
|
|
|
|
|
|
|
|
|
|
|
|
|
1.000
|
0.000
|
lithology
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1.000
|
*represents a significant correlation at the 0.01 level (bilateral). |
** represents a significant correlation at the 0.05 level (bilateral). |
According to the correlation coefficient matrix, lithology has a great influence on K+, HCO3-, and SO42- content and electrical conductivity in the water-rock interaction, showing a strong positive correlation.
The reaction time was positively correlated with Cl- content. In the aqueous solution, Cl- comes mainly from the dissolution of rock salt (NaCl) or other chlorides (MgCl2, CaCl2). Chlorine salts are highly soluble and do not easily precipitate from water. The reaction temperature has a great influence on the content of Na+ and Mg2+. With increasing reaction temperature, the Na+ content increases, while the Mg2+ content shows the opposite trend. Ca2+ and Mg2+ in aqueous solution come mainly from the dissolution of dolomite and calcite. The solubility of dolomite and calcite decreases with increasing temperature. Although the reaction temperature is negatively correlated with Ca2+, the correlation is not strong, indicating that there are other factors affecting the Ca2+ content in the aqueous solution.
There is a strong positive correlation between Ca2+ and HCO3- in the aqueous solution because carbon dioxide in an aqueous solution will increase the solubility of CaCO3. The content of HCO3- depends on the partial pressure of CO2 in aqueous solution. Figure 7. depicts the ionic ratio relationship between Ca2+ and HCO3-. The results showed that only one water sample point was between 1:1 and 1:2, and the mole ratios of Ca2+ and HCO3- in other reaction solutions were all below the 1:2 line, indicating that calcite in aqueous solution may dissolve faster than dolomite, and calcite reaches supersaturation earlier than dolomite and precipitates calcite. The dissolution of dolomite also reacts with CO2 in an aqueous solution to produce HCO3-, making the mole ratios of Ca2+, Mg2+, and HCO3- between 1:1 and 1:2. + in aqueous solution.
The ratio relationship between Ca2+, Mg2+, and HCO3- is shown in Fig. 8. In the water-rock interaction solution, only part of the samples are between 1:1 and 1:2 under the reaction condition of 20 ℃, and most of the sample points are below the ratio of 1:2, suggesting that solutions of dolomite and calcite are only part of the source of HCO3-. In addition, the solubility of Na2CO3 is very high. When the concentration of Na+ in an aqueous solution is high, the content of HCO3- will exceed the upper limit related to Ca2+.
K+ was positively correlated with Na+, HCO3- and SO42-, Na+ was positively correlated with Cl-, and HCO3- was positively correlated with SO42-, suggesting that they have the same material source or formation. The K+ in aqueous solution comes mainly from the dissolved K+ of illite:
Figure 9. shows the scattering of γNa/γCl in the aqueous solution under different reaction conditions. The results show that most of the points fall below the 1:1 line. The content of Na+ in the aqueous solution is higher than the content of Cl-, especially in the aqueous solution with a higher reaction temperature, indicating that rock salt is not the only source of sodium ions in aqueous solution, and sodium ions also come from the dissolution of plagioclase and sodium montmorillonite. An increase in reaction temperature will promote the dissolution of these minerals:
SO42- in solution usually comes from the dissolution of gypsum or other sulfates. In addition, the oxidation of pyrite in oil shale and surrounding rock leads to the existence of water-insoluble sulfur in water in the form of SO42-.
Fluorine comes mainly from the hydrolysis of fluorite (CaF2) and apatite (Ca5(Cl,F,OH)(PO4)3). F- in aqueous solution can react with Ca2+ to form CaF2. The solubility product of CaF2 is 10-10.4. Therefore, the Ca2+ concentration is an important factor in controlling the F- content in aqueous solutions.
Analysis of the characteristics of water chemical evolution.The Piper diagram can reflect the chemical composition of water-rock interactions in aqueous solution, as shown by the Piper diagrams of different types of aqueous solutions (Fig. 10). With the progress of the reaction, the hydrochemical type of oil shale water melts gradually and changes from HCO3-Ca type to HCO3·Cl-Ca·Na type at 20 ℃ and 50 ℃ and from HCO3-Ca·Na type to HCO3-Na·Ca type at 80 ℃ and at 100 ℃, HCO3-Na·Ca type to HCO3·Cl-Na type. The hydrochemical types of oil shale ash aqueous solution transition from HCO3·SO4-Ca type to HCO3-Ca·Na type at 20 ℃, from HCO3-Ca·Na + K type to HCO3-Na + K·Ca type at 50 ℃, and from HCO3-Na + K type to SO4·HCO3-Na + K type at 80 ℃ and 100 ℃. At 20 ℃, the hydrochemical type of oil shale surrounding rock solution transitions from HCO3-Ca type to HCO3·Cl-Ca·Na type and remains HCO3-Ca·Na type at 50 ℃. At 80 ℃ and 100 ℃, the solution transitions from the HCO3·Cl-Ca·Na type and HCO3-Na·Ca type to the HCO3·Cl-Na·Ca type, respectively.
With increasing temperature and the progress of the reaction, the content of Ca2+, Mg2+ and HCO3- in aqueous solution gradually decreases, while the content of Na+ and K+ increases. The content of SO42- in different aqueous solutions is obviously different. The content of SO42- in the oil shale water solution is the highest, followed by the oil shale water solution, and the content of SO42- in the oil shale surrounding rock water solution is the lowest.
Analysis of mineral saturation index change. The chemical composition of an aqueous solution is affected mainly by mineral dissolution and precipitation in the process of water-rock interaction. The saturation index (SI) is a commonly used index to judge the saturation of mineral components in solution. Figure 11. shows the variation in the saturation indices of calcite, stone salt, dolomite, gypsum, and fluorite with total dissolved solids (TDS) in water-rock interaction aqueous solutions.
As the figure shows, the saturation index of calcite, salt, dolomite, gypsum and fluorite in each aqueous solution is less than zero, indicating that these minerals are unsaturated in the aqueous reaction solution and can continue to dissolve. In the aqueous solutions of different reaction substances, except for dolomite, the saturation index of other minerals shows little change, presenting a rising trend with the increase in TDS of aqueous solutions and is closer to 0.