Energy is one of the important factors supporting the rapid and stable development of China’s economy, and coal resources occupy an important dominant position in China’s energy structure . The continuous loss of resources accelerates the process of deep mining of coal resources, while deep mining of coal faced with the complex environment of “three highs and one disturbance”  , complex environment makes the deep mining of coal in China is directly threatened by the karst water damage in the lower part. In recent years, China’s karst water disasters frequently occur, causing huge loss of life and property. Therefore, it is particularly important to correctly understand the hydrogeological conditions of the mining area. Hydrogeochemical method is one of the most applications and effective methods to analyze the hydrogeological conditions of research areas.
Group A in Huainan Coalfield is rich in coal resources and good in coal quality. At present, the mining in Panxie mining area is mainly concentrated in the mining areas of Gubei, Zhangji and Xieqiao, but it is vulnerable to karst water damage due to the direct contact between the bottom plate and the lower karst aquifer, and relatively difficult to mines. According to statistics, there have been 19 floor karst water inrush accidents in Huainan coalfield in the past 10 years . Previous studies on hydrogeochemistry in Huainan coalfield mostly focused on water source identification, water chemistry, water flow infiltration and hydrodynamic characteristics  -  , and the study on the genesis of aquifers and the effects on environmental formation less.
Therefore, this paper collects 199 Taiyuan karst water samples from the three mines of Xieqiao-Zhangji-Gubei in the past 10 years, using multivariate statistical analysis, water chemistry (piper diagram, ion ratio) and hydrogeochemical simulation method to analyze its chemical characteristics and to explore the causes in panxie mining area coal seam group A bottom karst water. The research results can provide a theoretical basis for the prevention and control of karst water damage in coal seam of Group A in Huainan coalfield.
2. Overview of the Study Area
At present, the coal mining of Group A in Huainan coalfield is mainly concentrated in the three mines of Xieqiao-Zhangji-Gubei in Panxie mining area. The overall east-west and south-north boundary faults of the three mines are changfeng fault, chenqiao fault, fufeng fault, shangtang-minglongshan fault. There are unevenly distributed secondary wide and slow folds and faults in the area, and the overall trend is EW and NNE. The strata in the area belongs to the typical Huainan sedimentary strata, and the bottom to the top is Cambrian, Ordovician, Carboniferous, Permian and Quaternary loose layers. The total thickness of the Taiyuan Formation is 99.99 - 129.11 m, with an average of 113.99 m; limestone, mudstone, siltstone and thin coal are the main components, including 10 - 13 layers of thin layer of limestone. The Initial water level by xie qiao, zhang ji and gu bei during detailed and detailed exploration is all higher than the surface. The karst water level is 23.08 - 27.06 m; the unit water is 0.0000405 - 1.764 L/s·m, the water-rich water is weak to medium, and the water temperature is 29˚C - 36.5˚C. It does not touch with the upper cover, and the replenishment conditions are poor.
The aquifers affecting coal mining in group A are mainly coal-series sandstone fissure aquifer, taiyuan inter-layer limestone karst fissure confined aquifer, Ordovician carbonate karst fissure confined aquifer and Cambrian carbon rock karst-fissure confined aquifer. Due to the direct contact between the limestone karst-fissure aquifer in taiyuan formation and the lower coal floor of group A, there is a great potential flood hazard.
There are currently eight pairs of mines in the Panxie mining area that have been put into production, and accumulated more geological and hydrogeological data. There are 4 karst collapse columns measured and 6 suspected, and group A of coal in this area has a large mining prospect. Therefore, it is of great significance to study the chemical composition and genesis of karst fissure water in taiyuan formation at the bottom of group A coal in this area for the prevention and control of karst water damage.
3. Sampling and Analysis
A total of 199 karst water samples from the Taiyuan Formation of the lower coal of Group A were collected, and their distribution is shown in Figure 1. The pH, total dissolved solids (TDS), conductivity and water temperature were measured in site; deionized before sampling wash the 1 L polyethylene bottle 3 times in water, rinse it three times with the sample to be taken, seal the sample, and bring the collected water sample back to the laboratory for routine ion measurement. and by acid-base titration, Cl− and by ion chromatography, Ca2+ and Mg2+ using EDTA titration, K+ + Na+ was measured by flame atomic absorption spectrophotometry. After testing, the relative error of the anion and cation balance of all water samples is less than 5%, which meets the standard.
Figure 1. Sampling point distribution map in Xieqiao-Zhangji-Gubei Mine Water of Huainan Coalfield.
4. Hydrochemical Characteristics of Karst Water
4.1. Statistical Characteristics of Different Components
Through the statistical analysis of the 199 Taiyuan water samples of the Xieqiao-Zhangji-Gubei mine, the statistical analysis results are shown in Table 1. The total dissolved solids (TDS) content ranges from 708.37 mg/L to 4061.35 mg/L, the average value is 2544.43 mg/L, the cations are dominated by Na+ + K+, the Ca2+ ion content is the second, the Mg2+ content is the least; the anion content was mainly Cl−, followed by , and , the average pH value is 8.73, is alkaline water.
The coefficient of variation is the ratio of the standard deviation to the average. The ions with a coefficient of variation greater than 0.5 are Ca2+, Mg2+, , , , indicating that there is a certain difference in the spatial distribution of each ion content. Among them, the coefficient of variation of Ca2+ and is much higher than 1, indicating that the spatial distribution is more dispersed and the ion concentration varies greatly in space.
The water quality types of karst water in the research area are shown in Figure 2, in which Na+ + K+ and Cl− are dominant. The water quality types are mainly distributed in the lower right side of the diamond. The water quality type is mainly Na+ + K+, Cl·HCO3-Na + K, Cl·SO4-Na + K, respectively accounted for 21.61%, 44.22%, and 23.31%, and other water quality types accounted about for 10%.
Figure 2. Piper map of research area.
Table 1. Chemical statistic results of karst water in Taiyuan Formation of study area.
4.2. Spatial Distribution Characteristics of Different Components
The TDS content in the study area was between 708.37 - 4061.35 mg/L, and the centralized distribution was between 2000 - 3000 mg/L. With the deepening, the TDS content increased, as shown in Figure 3(a), and the TDS content (c) with depth (h) the growth relationship is shown in the following Formula (1), which shows that with the increase of depth, the TDS content is continuously accumulated, and the water quality changes further. With the change of depth, Cl− shows a linear decrease trend, as shown in Figure 3(b), and the ion content is mostly distributed between 750 - 1250 mg/L. The relationship between the ion content and the depth change between is shown in Equation (2), which indicates that the aquifer may have replenishment; the change trend of ion content of and with depth is just the opposite and their distribution is shown in Figure 3(c) and Figure 3(d), with the increase of depth, the content of showed a slow growth trend, and showed a slow decreasing trend, indicating that the hydrochemical environment of aquifers is constantly changing with the increase of depth.
According to the above analysis, with the increase of vertical depth, TDS and show an increasing trend, Cl− and show a decreasing trend, which is related to the hydrochemical environment where groundwater runoff is, the condition of diameter and discharge, lithology, water and rock action and other factors. When groundwater flows through different areas, water-rock reactions will occur with different minerals, and the content of ions will change correspondingly. The total amount of ions continues to accumulate, but the Cl− content shows a decreasing trend with the depth, which is related to the situation that the aquifer is dominated by interlayer runoff, but some areas have vertical recharge through cracks. The increase of content and the decrease of content are related to the water chemical environment of the aquifer,
Figure 3. Relationship between ion content and depth.
the deep mining of the coal seam makes the originally closed reducing environment become a partially open oxidizing environment, which provides corresponding environmental conditions with the chemical reactions of related minerals. Therefore, as the depth increases, the flow rate slows down, the water-rock effect is further sufficient, and the water chemical environment is further changed, so that the contents of TDS, Cl−, and in the aquifer are spatially Make a certain change.
5. Genetic Analysis
5.1. Genetic Effect
According to the ratio of , the specific conditions of dissolution of carbonate rock and salt rock in aquifer can be judged . As shown in Figure 4(a), the ratio of in the water sample of the study area falls below the 1:1 line, and a few are above the line, indicating that the aquifer is mainly dissolved by salt rock and contains a small amount of carbonate rock to dissolve. The taiyuan formation in the lower part of group A coal seam is mainly composed of salt rock and calcite, with less dolomite content. Therefore, the dissolution of this area is mainly caused by the dissolution of salt rock.
Figure 4. Main ion ratio map of Karst water in study area.
According to the ratio of , the source of Ca2+, Mg2+ and ions in the aquifer can be determined . As shown in Figure 4(b), some water samples fall near the line 1:1, indicating that the dissolution of sulfate in the aquifer is not the only source of Ca2+, Mg2+, and may be derived from the dissolution of carbonate, the ratio of some water samples is much less than 1, indicating that there are other sources of in the aquifer. The Taiyuan Formation in Huainan Mining Area has less sulfate and more pyrite, so the sulfate is mostly derived from pyrite. The oxidation of the equation is:
The aquifer is entirely in a reducing environment as a whole, but due to deep coal seam mining, the aquifer becomes a partially open oxidizing environment, allowing gases such as O2 and CO2 enter the aquifer, thus enabling pyrite to generate under the action of oxygen and water.
According to the relationship between the ratio of and the line 1:1, the sources of Ca2+, Mg2+ and ions in aquifers can be inferred  , it can be seen from Figure 4(c) that the water sample points in the aquifer are distributed on both sides of the 1:1 line, indicating that the dissolution of carbonate rocks is not the only source of Ca2+, Mg2+ and , according to the spatial distribution of the water chemical components in the previous section, the ion content in the aquifer is concentrated in the position of −500 - −700 m, and the whole is in a closed reduction environment, making ions produce under the action of desulfuric bacteria.
5.1.3. Cation Exchange and Adsorption
Different cations in groundwater will alternately adsorb under certain conditions, according to is close to −1 to determine the strength of its alternate adsorption . As shown in Figure 4(d), the slope of the asymptote of the water sample in the study area is −0.48, and the negative correlation is weak, indicating that there is alternating adsorption of cations in the aquifer, but its effect is not obvious.
5.2. Simulation Analysis
The reverse simulation in the reaction path can verify the chemical reaction and composition of the material in the aquifer from the reverse   . Aquifer in the study area do bedding movement from north to south, Therefore, according to the run-off path of the study area (Figure 5), E121 is selected as the starting point and E157 is used as the termination point. The specific data is shown in Table 2, and the simulation software is used to calculate the saturation index, the results are shown in Table 3, it can be seen that the saturation index of vermiculite, calcite and dolomite in the aquifer is greater than 0, indicating that it is saturated in the aquifer; anhydrite, CO2, rock salt and O2 are all less than 0 and are in a dissolved state. On the run-off path, the saturation index of calcite, vermiculite, dolomite, anhydrite and gypsum is continuously decreasing, indicating that the water-rock interaction along the way is continuing.
Calculate the saturation index and the geological data of the mining area, and determine the “possible mineral phase” required for the simulation as calcite, dolomite, pyrite, salt rock, O2, CO2, H2S, Fe(OH)3 and a cation exchanger (NaX, CaX2).
Figure 5. Map of reverse simulation path in study area.
From this, the inversion simulation results are shown in Table 4. From the table, it can be seen that the CaCO3 value is less than 0, and the CaMg (CO3)2 is greater than 0, indicating that calcite precipitation and dolomite dissolution occur in the aquifer of the study area, the phenomenon of dolomite dissolution is different from the calculation of the saturation index, which is related to the phenomenon of “non-conformity dissolution” in this area . In this area, rock salt dissolution occurred mainly along the runoff path, and the contents of Cl− and Na+ increased, but cation exchange, dissolved CaX2, precipitated NaX, decreased Na+, and increased Ca2+ content. The sum of the cation exchange and the precipitation of calcite is still less than the dissolved amount of rock salt. Therefore, Na+ is still growing, but the content is less than Cl− content; FeS2 and O2 are dissolved, Fe(OH)3 is in a precipitated state, without the production of H2S gas, indicating that the oxidation of pyrite occurs in the aquifer, and the content of increases, but the desulfurization is not very obvious. From this, it can be inferred that the leaching action mainly occurs in the aquifer, and there are alternating cation adsorption and desulfurization, which is consistent with the analysis results above.
Table 2. Hydrochemiscal data of the simulation points in different zones from different aquifers.
Table 3. Mineral saturation index SI of the water sample points in Taiyuan formation.
Table 4. Summary of simulation results of karst water simulation path in Taiyuan formation (unit: 10−3 mmol/L).
Note: Positive values represent the amount of dissolution, negative values indicate the amount of precipitation, indicating that there is no such value in the simulation calculation.
Through the analysis of the chemical characteristics and genesis of the karst water in the coal of group A in the study area, the following conclusions are drawn:
1) The karst water cations of Taiyuan Formation are dominated by Na+ + K+, followed by Ca2+, Mg2+, and the anions of Cl− and . The ion content is ;
2) According to the relationship of ion ratio, it is known that the leaching effect mainly occurs in the karst aquifer of Taiyuan Formation, the dissolution of rock salt and the oxidation of pyrite are dominant, and desulfurization and cation alternate adsorption are present;
3) The water chemistry occurring in the aquifer was reversed by the runoff path in the simulation software, and the results were consistent with the previous genetic analysis, which also confirmed the rationality of the above analysis from the side.
 Li, T., Li, W.P., Gao, Y., et al. (2010) Characteristics of Karst—Fissure Water Bodies Deeply Seated in the Floor of No. 6 Coal Seam in Yangzhuang Coal Mine. Journal of Mining & Safety Engineering, 27, 98-103.
 Ma, L., Qian, J.Z., Zhao, W.D., et al. (2013) Multivariate Statistical Analysis of Chemical Characteristics of Groundwater in Gubei Coal Mine. Journal of Hefei University of Technology, No. 12, 1495-1498.
 Chu, T.T., Ge, T., Liu, G.J., et al. (2014) Coal Group A Floor Limestone Component Structure and Hydrochemical Characteristics in Xieqiao Coalmine, Huainan Mining Area. Coal Geology of China, No. 3, 27-31.
 Yu, Y.L., Song, X.F., Guo, J., et al. (2017) Characteristics and Controlling Factors of Water Chemistry in Surface Waters in the Yellow River Delta. Journal of Arid Land Resources and Environment, No. 10, 58-64.
 Chen, L.W., Xu, D.Q., Yin, X.X., et al. (2017) Analysis on Hydrochemistry and Its Control Factors in the Concealed Coal Mining Area in North China: A Case Study of Dominant Inrush Aquifers in Suxian Mining Area. Journal of China Coal Society, 42, 996-1004.
 Luan, F.J., Zhou, J.L., Jia, R.L., et al. (2017) Hydrochemical Characteristics and Formation Mechanism of Groundwater in Plain Areas of Barkol Yiwu Basin, Xinjiang. Environmental Chemistry, 36, 380-389.
 Matthew, M.U. and John, M.S. (2001) Tracing Regional Flov Paths to Major Springs in Trans-Pecos Texas Using Geochemical Data and Geochemical Models. Chemical Geology, 179, 53-72.
 Meng, S.X. and Maynard, J.B. (2001) Use of Statistical Analysis to for Mulate Conceptual Models of Geochemical Behavior: Water Chemical Data from the Botucatu Aquifer in Sao Paulo State. Brazilian Journal of Hydrology, 250, 78-97.
 Kenoyer, G.J. and Bowler, C.J. (1992) Groundwater Chemical Evolution in a Sandy Silicate Aquifer in Northern Wisconsin: 2. Reaction Modeling. Water Resources Research, 28, 591-600.