0.   引　　言

高强度、大规模的煤炭开采活动破坏矿区原始岩层结构，扰动地下含水层，产生大量矿井水^[1-3]，据不完全统计，仅2014—2019年，我国由煤炭开采所产生的矿井水总量就高达288.09亿t^[4]，大量未经处理的矿井水被抽排到地面，造成严重水资源浪费和环境污染，同时伴随地面塌陷和采空区积水等一系列问题，大幅增加了矿山修复与治理的成本和难度，严重制约煤矿的可持续发展^[5-6]。

矿井水资源化利用技术手段的提出，为解决煤矿水资源浪费和污染问题开辟了新的途径，该技术切实贯彻了煤矿开采“安全、高效、绿色”发展理念^[7-9]。我国关于矿井水资源化利用方面的研究，成果丰硕，武强^[10]对我国矿井水资源化利用的存在的问题、研究进展、发展趋势做了详细的分析；何绪文^[11-13]一直致力于探究矿井水资源化综合利用和矿井水处理技术创新；顾大钊^[14-15]提出了煤矿地下水库的理论，为解决西部煤矿开采水资源短缺问题提供了有效的创新性技术手段；孙亚军^[5,16]根据神东矿区水文地质特征，分析了如何针对不同矿井水类型采取相应的处理手段，并提出具有针对性的保水采煤关键技术，为干旱半干旱生态环境脆弱区的煤矿开采和水资源保护提供宝贵的实践经验。

基于以上研究成果，以淮南煤田张集煤矿为研究区，首先利用Ca²⁺、Mg²⁺、Na⁺+K⁺、Cl⁻、SO₄ ²⁻、CO₃ ²⁻、HCO₃ ⁻七大离子含量和pH、溶解性总固体(TDS)、电导率(EC)、总硬度(TH)11个指标对矿区矿井水进行水化学特征分析，然后依据我国《地下水水质标准》，进行矿井水水质分级分类判别和评价，最后根据钠吸附比(SAR)和EC关系，并结合离子来源，探求张集煤矿矿井水资源化利用途径。张集煤矿矿井水水量大，地表积水面积广，探求张集煤矿矿井水资源化利用的途径，对矿区的可持续发展意义重大，另外，到目前为止，该矿区还未有矿井水资源化利用的相关研究，较以往研究，本研究不仅结合水化学特征分析，和遵循矿井水资源化利用“分级分类”原则，而且进一步判定了研究区矿井水实际可用途径，为实现张集煤矿区矿井水资源化利用提供一定的理论依据。

1.   研究区概况

张集煤矿位于中国安徽省淮南市淮南煤田潘谢矿区的西南部(图1)，矿区东西走向长约12 km，南北倾斜宽约9 km，面积约71 km²。研究区地表水系发达，西淝河从西到东贯穿整个矿区，内涝区从西到东横跨了研究区，由于采矿活动矿井水抽排，积水区面积也不断扩大，更加重了矿区水域承载负担，截止2018年预计矿区总积水面积可达53.17万m²，积水量预计将达到25.44万t。地下水埋藏情况变化较大，张集煤矿区由3个主要充水含水层组成，分别为新生界松散含水层、煤系砂岩裂隙含水层、灰岩岩溶含水层，地下水在新生界松散层中以水平运动为主，类型上属于孔隙承压水，是本区的水源补给。

图  1  研究区示意

Figure  1.  Sketch map of study area

下载: 全尺寸图片 幻灯片

2.   取样和方法

根据现有资料整理分析，分别在张集煤矿的西区、中央区和北区3个分区选取33个取样点（S1～S33），如图2所示，取样时间从2017-01-18—2018-07-25，其中同一个取样点不同时间取样多次，一共取得82个矿井水水样，并按采样的先后顺序将样品从ZJ135～ZJ227进行编号，水样来源有蓄水池、采空区、排水孔、水源井、放水孔、塌陷区等，水样来源几乎包含矿区所有矿井水来源。

图  2  研究区采样点

Figure  2.  Sampling points in the study area

下载: 全尺寸图片 幻灯片

研究区矿井水水样的采集和测定依照《煤矿水水质分析的一般规定》及《地下水质量标准》。采取水样前需现场测定水温、pH、EC值并记录，pH采用玻璃电极法测定，水温用水温计测量，EC值用电导率仪可直接测定；现场测定结束后采集水样，采样前需用去离子水冲洗采样容器(按标准采用1 L聚乙烯瓶)3次，然后用待采矿井水冲洗容器3次，最后采集水样并装样密封。将水样带回实验室进行离子含量测定，TDS值在105 ℃下烘干称重得到，TH值采用EDTA容量法获得，采用酸碱滴定法测定HCO₃ ⁻和CO₃ ²⁻；离子色谱法应用于Cl⁻和SO₄ ²⁻的测定；Ca²⁺和Mg²⁺采用EDTA滴定法；而Na⁺和K⁺则用火焰原子吸附分光光度法测定。

3.   水化学特征

3.1   一般水文化学

82个矿井水样品的常规水化学指标检测结果，由Ca²⁺、Mg²⁺、Na⁺+K⁺、Cl⁻、SO₄ ²⁻、CO₃ ²⁻、HCO₃ ⁻七种离子含量和pH、TDS、EC、TH共11个指标组成，然后对检测结果做一般性的水化学特征分析，统计分析结果见表1。由表1可看出矿井水的pH于7.16～12.24，均值为8.33，中值为8.25，矿区矿井水均为碱性水，pH最大的水样检测含有大量OH⁻，无HCO₃ ⁻，为张集煤矿锅炉房排水，另外pH的变异系数仅为0.08，说明研究区矿井水pH变化很小，基本保持在8左右。TDS含量为255.47～7 291.74 mg/L，平均值为2 346.60 mg/L，变异系数为0.51，多数水样的TDS与中位数2 584.22 mg/L接近。

表  1  张集煤矿矿井水一般水化学分析结果

Table  1.  General hydrochemical analysis results of mine water in Zhangji Coal Mine

参数	最大值	最小值	平均值	中位数	标准差	变异系数
Ca2+质量浓度/(mg·L−1)	530.66	1.60	45.76	33.48	62.18	1.36
Mg2+质量浓度/(mg·L−1)	145.86	0.00	16.73	16.53	17.77	1.06
K++Na+质量浓度/(mg·L−1)	2 964.18	1.00	887.54	1 009.40	514.69	0.58
Cl−质量浓度/(mg·L−1)	3 247.22	7.09	752.71	831.31	479.48	0.64
SO4 2−质量浓度/(mg·L−1)	1 594.60	0.00	269.69	186.36	295.25	1.09
CO3 2−质量浓度/(mg·L−1)	420.07	0.00	29.05	0.00	72.49	2.50
HCO3 −质量浓度/(mg·L−1)	3 368.67	0.00	690.03	395.11	698.65	1.01
TDS质量浓度/(mg·L−1)	7 291.74	255.47	2 346.60	2 584.22	1 199.76	0.51
pH	12.24	7.16	8.33	8.25	0.66	0.08
EC/(μS·cm−1)	10 560.00	466.00	3 309.66	4 310.00	1 797.41	0.47
TH/(mg·L−1)	1 321.06	8.01	183.96	176.14	183.80	1.00

下载: 导出CSV 

| 显示表格

所有水样中阳离子均以K⁺+Na⁺占绝对优势，Ca²⁺次之，Mg²⁺最少，而阴离子以Cl⁻、HCO₃ ⁻为主，两者相差不多，其次是SO₄ ²⁻，CO₃ ²⁻含量最少，但研究区所有离子含量的变异系数皆大于0.5，说明离子含量极不稳定，极易随时间和空间的变化而发生改变，其中CO₃ ²⁻含量少，但变异系数达到2.5，在整个采样范围内随时空变化特征最为明显。

3.2   主要离子成分

piper三线图可用来指示研究区矿井水水样的主要离子组成和化学成分，并可作为判别水样类型的方法^[17]。在研究区33个取样点共取得33组水样得到82个矿井水样品，图3为33组水样的piper三线图，由图3可知，87.8%的水样碱金属离子大于碱土金属离子，即Na⁺+K⁺>Ca²⁺、Mg²⁺，33组水样(S1～S33)阳离子值绝大多数落在Na+K轴，比例在80%～98%，而Ca²⁺值变化范围较大，阴离子含量则以Cl⁻和HCO₃ ⁻+CO₃ ²⁻为主要成分，有超过30组水样点分布在Cl轴80%和HCO₃+CO₃轴80%交界处。由图3可明显看出，除了第27组水样2个样品类型不同(分别是Cl·HCO₃-Na+K型和Cl·HCO₃-Ca型)，其余同一取样点所采水样类型基本一致，研究区矿井水水样水化学类型主要有Cl·HCO₃-Na+K型、HCO₃·SO₄-Na+K型、Cl·HCO₃-Ca型和SO₄-Na+K型，所占比例分别为59.8%、17.6%、13.4%、9.2%。

图  3  研究区矿井水主要离子Piper图

Figure  3.  Piper diagram of main ions in mine water of study area

下载: 全尺寸图片 幻灯片

3.3   水化学特征的动态变化

通过比对于不同采样时间、同一采样点采取的矿井水指标含量，探讨研究区矿井水指标含量的动态变化特征，推断矿井水指标变化规律及受控因素。选取了1412A工作面灰岩增补孔(图4a)、北区水仓蓄水池(图4b)、西二定向长钻孔(图4c)、西二疏水巷排水孔(图4d)、西三疏水巷定向长钻孔(图4e)、西二封闭墙排水孔(图4f)6个采样点，共采取33个矿井水水样作为分析对象，绘制图4。

图  4  化学指标-时间变化

Figure  4.  Hydrochemical indicators-time variation chart

下载: 全尺寸图片 幻灯片

由图4中4a、4c、4d、4e、4f五张柱状图可看出同一采样点的各项离子浓度随时间的起伏变化不明显，范围基本在30 mg/L以内，说明同一采样点的水源补给较单一、稳定，而图4b蓄水池中矿井水离子指标含量从2017年6月至2017年12月出现显著增加而后又明显减少，据分析其原因可能是与研究区气候变化有关，6月到12月研究区自然降水逐渐减少，采样点位于蓄水池，受外界影响显著，池中矿井水经过6个月的蒸发作用，离子浓度逐步升高，而后突然降低可能是受采矿活动矿井水排放增大，水量增多使离子浓度在短时间内降低。

3.4   相关性分析

对变量进行相关性分析的目的，在于表征两个变量之间的相互关系和变化趋势^[18]，另外能够反映变化趋势的方向和程度，并在一定程度上对试验样本做来源预测分析^[19]。试验对82个水样的10组指标数据依次采取双变量进行相关分析，由于Pearson和Spearman相关系数只能应用于连续变量，而水样数据为无序离散变量，因此，选用Kendall公式计算水化学指标参数的相关系数^[20-21]，然后依据所得相关系数值对样品水化学指标数据之间的关联性进行分析，推断各离子指标之间的相互关系度，每种水化学指标的影响作用范围，水质检测标准TDS最主要的控制因素是哪项指标，并根据实际取样环境进一步推断离子来源。

由计算结果(表2)可得，矿井水水样品中对TDS值贡献度最高的是钠离子和钾离子，相关系数值高达0.821，结果极显著，且TDS总体上与各离子指标含量相关性最大，阴离子中与TDS相关系数值比较高的是Cl⁻、SO₄ ²⁻和HCO₃ ⁻，反映这3种离子在TDS中所占比重较高；而水质检测标准TH与钙、镁离子密切相关，这也与实际测量结果一致，同时，钙、镁离子之间呈正相关，两者的来源应该多是碳酸盐岩溶解，而二者与CO₃ ²⁻、HCO₃ ⁻的含量呈负相关关系也证实了这一点，前者与后者会生成沉淀，因此CO₃ ²⁻、HCO₃ ⁻的含量越多，Ca²⁺和Mg²⁺含量越少，而碳酸盐岩溶解度也与游离离子含量成反比。

表  2  水化学指标相关系数矩阵

Table  2.  Correlation coefficient matrix of hydrochemical indexes

离子	Ca2+	Mg2+	Na++K+	Cl−	SO4 2−	CO3 2−	HCO3 −	EC	TH	PH
TDS	−0.229**	−0.215**	0.821**	0.456**	0.465**	0.225**	0.435**	0.786**	−0.196**	0.100
Ca2+		0.509**	−0.403**	0.063	0.080	−0.357**	−0.523**	−0.210**	0.847**	−0.268**
Mg2+			−0.307**	0.148	−0.071	−0.219*	−0.466**	−0.200**	0.663**	−0.133
Na++K+				0.341**	0.307**	0.301**	0.555**	0.681**	−0.368**	0.160*
Cl−					0.413**	−0.061	0.009	0.461**	0.129	−0.046
SO42-						−0.193*	0.060	0.357**	0.027	−0.221
CO32-							0.275**	0.207*	−0.316**	0.653**
HCO3−								0.428**	−0.515**	0.133
EC									−0.189*	0.094
TH										−0.231**
注：**相关性在0.01水平上显著，*相关性在0.05水平上显著。

下载: 导出CSV 

| 显示表格

PH对水质是否合格至关重要，而根据计算结果发现，CO₃ ²⁻离子是决定本次水样PH的最主要因素，这也为后期矿井水水处理提供理论依据，如果想要改变矿井水的酸碱度，可通过去除矿井水中的CO₃ ²⁻来达到目的；EC一点程度上可反应TDS，同样可作为水质检验的标准，由相关系数值也可看出，EC与TDS的值大致相同，说明二者作用意义基本一致；SO₄ ²⁻跟所有指标的相关性皆小于0.5，即与其他指标的关联性都不大，说明其缘由可能是农业灌溉、生活污水等人类活动而产生的废水汇入矿井水，属于外来物质，与矿区矿井水本身来源关系不显著。

3.5   主成分分析

完成了以上样品指标的水化学特征分析，我们还需要总体把握研究区矿井水水样的水化学特征，由于本次试验选取的样品指标较多，为减少变量个数，降低工作强度，同时尽可能完整保留原始数据信息，达到综合研究全部样品数据的目的，我们采用了主成分分析方法。将水样重新进行组分分析，得到2个新的综合变量F₁、F₂，主成分F₁、F₂的方差贡献率分别为35.121%、25.503%，然后对Ca²⁺、Mg²⁺、Na⁺+K⁺、Cl⁻、SO₄ ²⁻、CO₃ ²⁻、HCO₃ ⁻七个离子指标的主成分综合得分值再分析，并绘制分析变量荷载分布图(图5)。

图  5  矿井水分析变量荷载分布

Figure  5.  Distribution of variable load of mine water analysis

下载: 全尺寸图片 幻灯片

由图5可知，Na⁺+K⁺、Cl⁻在2个主成分上都有较高荷载值，即Na⁺+K⁺、Cl⁻是研究区矿井水的主要离子组分，反映张集矿矿井水出现高度“咸化”，水质逐渐恶化，结合采样点水文地质资料获悉，Na⁺+K⁺离子含量较高由于煤炭开采、钻井勘探等原因，导通各含水层，富含Mg²⁺、Ca²⁺离子的地下水流经岩土层将Na⁺+K⁺置换出来，不断地进行阳离子交换作用，使得矿井水含盐量逐渐增高;综合来看，HCO₃> ⁻、CO₃ ²⁻荷载得分值与SO₄ ²⁻、Mg²⁺、Ca²⁺荷载得分值在2条主成分轴上均相反，反映在水样中HCO₃ ⁻、CO₃ ²⁻离子的含量与SO₄ ²⁻、Mg²⁺和Ca²⁺离子的含量之间存在明显的相互作用，原始的封闭还原环境中，SO₄ ²⁻发生脱硫酸作用被还原成H₂S，使得HCO₃ ⁻含量增高，而由于采矿活动，使还原环境变为氧化环境，水中携带了大量的O₂和CO₂，导致碳酸盐岩被溶解产生Ca²⁺，而HCO₃ ⁻与CO₃ ²⁻是相辅相成的存在。

为消除变量之间的量纲影响，使其具有可比性，对Ca²⁺、Mg²⁺、Na⁺+K⁺、Cl⁻、SO₄ ²⁻、CO₃ ²⁻、HCO₃ ⁻指标含量数据作标准化处理，将标准化处理后的数据代入主成分计算公式(式(1)和式(2))，求得共82组主成分得分数据，并得到研究区离子指标主成分得分散点图(图6)。

图  6  主成分计算结果荷载投影图

Figure  6.  Load projection of principal component calculation results

下载: 全尺寸图片 幻灯片

主成分计算表达式如下：

式中：$ F_1 $、$ F_2 $为两个主成分；$ {x}_{i} $（$ i=\mathrm{1,2},\cdots, 7 $）为标准化处理后的7个离子指标含量。

图6显示，水样整体沿主成分F₁轴展布，说明水样来源含水层中存在脱硫酸作用和阳离子交换作用，导致矿井水盐度增高；另有多数水样点集中在F₂轴上半部分，结合图5可知，矿井水流经碳酸盐岩层，伴随石灰岩、白云岩溶解，水质发生“硬化”，综上所述，张集煤矿矿井水存在“咸化”、“硬化”现象，水质情况不容乐观。

4.   结果与分析

4.1   矿井水水质评价及离子来源分析

为遵循矿井水“分级应用”的原则，本研究参照以上研究结果及现行《地下水质量标准》(表3)，分别对研究区矿井水水样中的TH、TDS、SO₄ ²⁻、Cl⁻、Na⁺指标含量进行判别分类，然后根据单个指标分类，得到82个样品的综合水质分级结果。

表  3  地下水质量一般化学指标及限值(部分)

Table  3.  General chemical indexes and limits of groundwater quality (part)

类别	Ⅰ类	Ⅱ类	Ⅲ类	Ⅳ类	Ⅴ类
TH/(mg·L−1)	≤150	≤300	≤450	≤650	＞650
TDS/(mg·L−1)	≤300	≤500	≤1 000	≤2 000	＞2 000
SO4 2−/(mg·L−1)	≤50	≤150	≤250	≤350	＞350
Cl−/(mg·L−1)	≤50	≤150	≤250	≤350	＞350
Na+/(mg·L−1)	≤100	≤150	≤200	≤400	＞400

下载: 导出CSV 

| 显示表格

图7显示，SO₄ ²⁻、Na⁺、Cl⁻和TDS指标的分类结果中第Ⅴ类所占例较大，Ⅴ类的样品个数分别为27、65、63和57，水质较差，而TH指标等级较低，共有36和40个Ⅰ类和Ⅱ类，说明研究区TH指标良好。每组样品水质等级划分如图7b所示，张集煤矿矿井水水样均不符合Ⅰ类水质标准，而Ⅱ、Ⅲ类占13%，Ⅴ类占比高达81%，Ⅳ类为6%。

图  7  水化学指标类别堆积柱形图及样品水质分类饼图

Figure  7.  Stacked column chart of hydrochemical index categories and pie chart of sample water quality classification

下载: 全尺寸图片 幻灯片

由水质评价结果，可知研究区离子组分含量较高，为判断研究区矿井水中主要组分的来源，了解研究区的水环境条件，以便为矿井水水质改善和资源化利用提供理论依据。

研究采用Gibbs半对数图探究离子成因机制，由张集煤矿矿井水水样Gibbs图(图8)显示，水样点主要分布在图的上半部分，多数落在蒸发结晶作用带，少数在岩石风化带，说明研究区矿井水受蒸发结晶和岩石风化控制作用明显，这也与矿井水的实际来源相符，研究区矿井水水样多数从地下抽排至地表，离子来自含水层围岩和第四纪松散层土壤，到达地表后被蒸发结晶，另外，取样时间几乎都在研究区降水较少时节，受大气降水作用影响较小，这也解释了检测结果中矿井水TDS含量偏高的原因。

图  8  矿井水Gibbs半对数图

Figure  8.  Gibbs semilogarithmic diagram of mine water

下载: 全尺寸图片 幻灯片

4.2   矿井水资源化利用

根据水质分级结果，研究区矿井水样品等级有Ⅱ、Ⅲ、Ⅳ和Ⅴ类4个等级，但由于指标不全，不能保证研究区矿井水其他指标达标，因此作为Ⅱ、Ⅲ类用水功能有待进一步研究，而可以确定的是，水样全部符合《地下水质量标准》Ⅴ类标准，即不宜做生活饮用水，水域功能可作为一般工业用水区、人体非直接接触娱乐用水区、农业用水或者一般景观水域等。

矿井水在地面积聚以后受蒸发作用明显，常规离子含量过高，植物生长需要适宜的水化学环境，因此需进一步判断是否可做农业灌溉用水。从相关性分析可知，组成TDS的主要离子是Na⁺+K⁺，由主成分分析可知研究区矿井水咸化严重，钠离子或钾离子含量过高对土壤和植物存在毒害作用，从而限制矿井水资源化利用，因此，采用钠吸附比(SAR)来进行水质检测，钠吸附比常应用于地表及地下水灌溉性水质评价^[22]，计算公式如下：

式中：$ {\rm{SAR}} $为钠离子吸附比；$ c $为离子浓度，mg/L。

研究利用SAR与EC之间的关系，进行张集矿矿井水盐渍化和钠离子毒害效应评价，并作为判断矿井水资源化利用可行性的依据。图9显示，整体来看研究区矿井水SAR和EC值基本偏高，虽然水样中有钠毒害较低风险的存在，但EC值皆在E1-E2界线之上，钠毒害风险较低但会受到盐渍化威胁，另外，79.3%的水样点分布在S3-E4区和S4-E4区，说明钠毒害及盐渍化风险均较高。即研究区矿井水不适宜用做农业灌溉用水，可作为工业生产、城市绿化、水体景观等用途。

图  9  SAR-EC风险判别分区

E1、E2、E3、E4—盐渍化低、中、较高、高风险;S1、S2、S3、S4—钠毒害低、中、较高、高风险

Figure  9.  Risk discrimination division of SAR-EC

下载: 全尺寸图片 幻灯片

5.   结　　论

1）张集煤矿矿井水TDS主要贡献阳离子为Na⁺和K⁺，阴离子主要组成是Cl⁻和HCO₃ ⁻，水化学类型中Cl·HCO₃-Na+K型占总数的59.8%，水中离子主要来源为蒸发结晶和岩石风化作用。

2）研究区矿井水仅达到Ⅴ类水质标准，且水质“咸化”“硬化”现象严重，存在较高的钠离子毒害与盐渍化风险，不适宜用做农业灌溉用水，参考历来煤矿矿井水资源化利用的成功案例和研究区自身特点，张集矿矿井水实际可用作道路洒水、矿区绿化、冷却降温、洗煤，向周边煤化工厂或工业园区提供工业用水等。

3）位于地面下的采样点，受外界环境影响较小，矿井水来源稳定，在治理过程中，严格把控水源，针对TDS及Na⁺、K⁺、HCO₃ ⁻含量较高的特点，选取合适的物理、化学或生物方法进行净化处理，提高水质，使矿井水得到更广泛的应用。