The samples used in this study from No. 11 mine include mixed water from the sump (denoted C1, C2, and C3), drilling into Cambrian limestone (H1, H2, and H3), roof sandstone fissures water (D1) and goaf water (K1). The positions of these sampling points are shown (Figure 1 and Table 1); All samples were analyzed at the Henan Geological Environment Monitoring Institute. Simple analyses were performed on samples C2 and C3 and the remainder were subjected to full analysis. The simple analysis includes the common ions, and the full analysis includes all ions, etc.

Test results of conventional hydrochemical components are shown Table 1. A piper diagram was constructed based on the six major ions was shown as Figure 2 [3] . Thus, four water sample groups within this analysis were characterized as HCO₃-Na (D1, K1, C1, and C3), three were HCO₃-Ca (H1, H3, and C2), and one was HCO₃-Mg (H2). The data presented in Table 1 shows that because the pH of all samples is between 7.1 and 8.0, this is neutral or weak alkaline water. In addition, the turbidity of mixed water in the central water sump (C1) is 20, while the remainder of samples are less than 1; this is because the central water sump contains a large volume of coal dust, which leads to an increase in turbidity. The TDS values for samples D1 and K1 were 2777.83 mg/L and 1156.65 mg/L, respectively, and so these are classified as brackish water (1000 mg/L < TDS < 3000 mg/L), while the rest of these samples have TDS values less than 1000 mg/L and so are classified as freshwater (TDS < 1000 mg/L) [18] [19] . The reason is that the D1 and K1 samples were taken from coal roof and goaf, respectively. The total hardness of these samples falls between 18 mg/L and 542.5 mg/L (average: 309 mg/L), within the range of hard water.

The data presented in Table 2 shows the test results of trace components. Obviously, Trace component content of NO 3 − and F⁻ in Ordovician limestone water samples are far lower than that of other water samples (D1, K1 and C1). In addition, content of NH 4 + in roof crack water samples is 70 times as other samples. Generally speaking, content of F⁻ in goaf (K1) and NH 4 + and F⁻ in roof crack, all exceed class V water standards as defined by Chinese Groundwater Quality Standards (GB/T 14848-93).

A correlation analysis was performed to measure the relationship between two groundwater variables and to demonstrate source consistency and variability. According to mathematical statistics theory, when r > 0.9, it indicates that there is a significant correlation between the two variables; when 0.6 < r ≤ 0.9, it indicates moderatet correlation between the two variables; When 0.4 < r ≤ 0.6, it indicates low level correlation between the two variables, similar, When r ≤ 0.46, it indicates little correlation between the two variables [20] [21] .

Pearson correlation coefficients for Cambrian limestone groundwater were calculated using the software SPSS (Table 3). These data show that the correlation coefficients of TDS and total hardness and K⁺ + Na⁺ are 0.985 and 0.988, respectively. The correlation coefficients of TDS and Ca²⁺ reached 1. And, the correlation coefficients of TDS and total hardness and HCO 3 − are 0.998 and 0.992, respectively. In addition, the correlation coefficients of TDS and total hardness and SO 4 2 − are 0.999 and 1, respectively. This indicates that TDS and total hardness, K⁺ + Na⁺, Ca²⁺, HCO 3 − , SO 4 2 − are significantly correlated. In other words, The concentration of TDS and total hardness in Ordovician limestone groundwater is mainly determined by the content of K⁺ + Na⁺, Ca²⁺, HCO 3 − , SO 4 2 − . This is mainly because the several ions, K⁺ + Na⁺, Ca²⁺, HCO 3 − , SO 4 2 − , control the nature of groundwater.

The data presented in Figure 1 show that mine water flows into a sump in the Ji 2 and Ji 4 mining areas. This flow then raises mine water into the −593 m west tunnel through the pump house, flows to the central water sump, and then into the ground. Results from the software Aq・QA used to study the water chemical ions in samples are shown in Table 4.

Simulation calculation value of a mixture of various sources of water (H2 + C3 + K1 + H3 + C2 + D1 + H1), based on Aq・QA software and measured value are shown in Table 4 and Figure 3. Obviously, except for SO 4 2 − and HCO 3 − , the difference between the calculated value and the measured value of most ions is small. In addition, the trend of the connection between the simulated value and the measured value is consistent. This shows that the software simulation results of mine water confluence can be used to estimate the content of ion in mine water.

From Table 4, we know that values of TDS changed from 667.16 mg/L (H2) to 795.3 mg/L when Cambrian groundwater (H2) was mixed with water from the Ji 4 water sump (C3) and the goaf (K1). Thus, TDS increased from 683.3 mg/L (H3) to 1375 mg/L while water quality changed from freshwater to brackish when the Cambrian limestone groundwater (H3) was mixed with that from the Ji 2 water sump (C2) as well as 16,150 roof water from the working face (D1). In addition, TDS fell to 1071 mg/L and total hardness increased to 305.9 mg/L when this was mixed with water from H2 + C3 + K1 and Cambrian limestone groundwater (H1) as the latter was combined with sake. The above results comprehensively show that the mine water with different ion concentrations will change its ion concentration due to the confluence and affect the overall quality of the water, ultimately.

A set of evaluation factors were created on the basis of the chemical composition of test results from sampling points. Thus, on this basis, those with greater impact on the water quality of chemical indicators were selected, as follows: U = {Cl⁻, SO 4 2 − , NH 4 + , NO 3 − , F⁻, Ba, TDS, Total hardness, Turbidity}. Chinese Groundwater Quality Standards (GB/T 14848-93) have established a set of five evaluation factors V = {I, II, III, IV, V} [22] ; Table 5 shows selected standard values for groundwater chemical grade classifications.

Take water sample H1 is as an example to state evaluation process. According to the comparison of the evaluation factor matrices and the evaluation rating matrices, the linear membership function is used to calculate the membership of each pollution factor [23] , the applied formula as follows:

μ i ( x ) = { 1 − ( m i + 1 − x ) / ( m i + 1 − m i )               m i − 1 < x < m i ( m i + 1 − x ) / ( m i + 1 − m i )                       m i < x < m i + 1 0                                                                         x ≤ m m − 1 , x ≥ m i + 1

where x is the measured value; i is evaluation rank vaule, μ is the membership of the evaluation factor.

Thus, the corresponding Fuzzy relationship matrix is calculated as follows:

T 1 = [ t i j ] = [ 0.883 0.117 0.000 0.000 0.000 0.000 0.300 0.700 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 1.000 1.000 0.000 0.000 0.000 0.000 0.000 1.000 0.000 0.000 0.000 0.000 0.034 0.966 0.000 0.000 0.000 0.000 0.075 0.925 0.000 1.000 0.000 0.000 0.000 0.000 ]

And the membership matrix of other water samples are calculated as in sample H1.

The AHP method [24] was initially used to compare each element and to establish a comparison matrix, as follows:

P = [ p i j ] = [ 0 − 1 0 − 1 − 1 1 − 1 − 1 0 1 0 1 − 1 − 1 1 − 1 − 1 0 0 − 1 0 − 1 − 1 1 − 1 − 1 0 1 1 1 0 − 1 1 0 1 1 1 1 1 1 0 1 0 1 1 − 1 − 1 − 1 − 1 − 1 0 − 1 − 1 − 1 1 1 1 0 0 1 0 1 1 1 1 1 − 1 − 1 1 − 1 0 1 0 0 0 − 1 − 1 1 − 1 − 1 0 ]

Thus, in this matrix p i j = { 1               p i   is more important than   p j 0               p i   and   p j   are just as important − 1             p i   is not important for   p j

The following formula is then used to build an optimal transfer matrix, q i k , q j k are the element in the comparison matrix Q.

q i j = 1 n ∑ k = 1 n ( p i k − p j k ) = 1 n ∑ k = 1 n ( p i k + p j k )

Q = [ 0.000 − 0.333 0.000 − 1.000 − 1.222 0.444 − 1.111 − 0.667 − 0.111 0.333 0.000 0.333 − 0.667 − 0.889 0.778 − 0.778 − 0.333 0.222 0.000 − 0.333 0.000 − 1.000 − 1.222 0.444 − 1.111 − 0.667 − 0.111 1.000 0.667 1.000 0.000 − 0.222 1.444 − 0.111 0.333 0.889 1.222 0.889 1.222 0.222 0.000 1.667 0.111 0.556 1.111 − 0.444 − 0.778 − 0.444 − 1.444 − 1.667 0.000 − 1.556 − 1.111 − 0.556 1.111 0.778 1.111 0.111 − 0.111 1.556 0.000 0.444 1.000 0.667 0.333 0.667 − 0.333 − 0.556 1.111 − 0.444 0.000 0.556 0.111 − 0.222 0.111 − 0.889 − 1.111 0.556 − 1.000 − 0.556 0.000 ]

Next, after an optimal transfer matrix is established, a judgment matrix is calculated using the formula r i j = exp ( q i j ) , as follows:

R = [ 1.000 0.717 1.000 0.368 0.295 1.560 0.329 0.513 0.895 1.396 1.000 1.396 0.513 0.411 2.177 0.450 0.717 1.250 1.000 0.717 1.000 0.368 0.295 1.560 0.329 0.513 0.895 2.718 1.948 2.718 1.000 0.801 4.240 0.895 1.396 2.432 3.395 2.432 3.395 1.249 1.000 5.295 1.18 1.743 3.038 0.641 0.459 0.641 0.236 0.189 1.000 0.211 0.329 0.574 3.038 2.177 3.038 1.118 0.895 4.738 1.000 1.560 2.718 1.948 1.396 1.945 0.717 0.574 3.038 0.641 1.000 1.743 1.118 0.801 1.118 0.411 0.329 1.743 0.368 0.574 1.000 ]

Finally, in order to determine the weight of water chemical indicators, the following formula is applied [25] :

w i = ∏ k = 1 n r i k n ∑ k = 1 n ∏ k = 1 n r i k n

w i = ( 0.062 , 0.086 , 0.062 , 0.167 , 0.209 , 0.040 , 0.187 , 0.120 , 0.069 )

Based on B = T i × w i , water quality can be evaluated as B = { b 1 , b 2 , ⋯ , b j } , where b_j is the corresponding grade corresponding to a maximum B value. T_i is the membership matrix of water sample, w_i is the weight of each evaluation factor [26] . On this basis, water quality evaluation of Cambrian limestone groundwater (H1) is B 1 = ( 0.394 , 0.079 , 0.250 , 0.111 , 0.167 ) and can be classified as type I. All other results are shown in Table 6.

The water samples are all tested by Henan Provincial Geological Environment Monitoring Institute with Grade A qualification, and the test results are reliable. All water quality assessments are graded according to China’s Groundwater Quality Standard (GB/T 14848-93), the criteria are consistent. The weights of the evaluation factors are determined using the AHP method, avoiding man-made influence. Based on the above three reasons, the evaluation results are reliable and can be used for the treatment and utilization of mine water.

Research indicates Cambrian limestone groundwater from NO. 11 mine is Class I, and the water quality is superior. This is mainly due to the excellent circulation of the Cambrian limestone groundwater. Then, roof water is ClassⅤ, inferiorin addition water from goaf and central water sump are Class III. This is because Sandstone fissure water has recharge and poor runoff conditions.

Data show that the large difference in water quality depends on gushing source. This means that roof water from the working face and Cambrian limestone groundwater should be stored separately before being distinctively treated and utilized according to their respective levels of quality. This approach will guarantee the direct use of Cambrian limestone water and will also reduce mine drainage treatment costs.