2016年12月30日，国家发改委、能源局对外公布的《煤炭工业发展“十三五”规划》中明确要求生态文明绿色矿区的建设要取得积极进展，要最大程度减轻煤炭生产开发对环境的影响，提升资源综合利用水平，其中矿井水的利用率要达到80%^[1]。由于我国煤炭资源主要采用地下开采方式，为避免突水灾害的发生，大量的矿井水必须排出井外，据统计我国矿井水每年排放量高达74亿m³，随着矿井水的大量疏排，地下含水系统的平衡与稳定受到影响，由于矿井水水质特殊、成分复杂，易对矿区生态环境造成影响，而水质较好的矿井水在多数情况下也未能得到有效合理的利用，造成地下水及洁净矿井水资源的浪费^[2]，因此，开展煤矿区地下水水化学特征及成因机制的研究，可为矿井水排水处理、资源化利用提供科学依据。20世纪90年代，为查明煤矿区水文地质条件，主要通过水文地球化学研究判断地下水补径排关系和含水层水力联系^[3-4]，后来逐渐广泛应用于矿井突水水源判别和水害防治方面^[5-7]。近些年，随着对矿区生态环境的重视以及对煤矿区地下水资源的保护和矿井水资源化利用的要求，水文地球化学的研究逐渐侧重于煤矿区地下水和矿井水成因机制方面^[8-9]。目前，国内外对地下水水文地球化学特征的研究方法可分为野外调查研究、水化学数据分析、模拟实验以及计算机数值模拟^[10]，其中水化学数据分析方法研究和应用最为广泛：Piper图解、Gibbs图解等水化学类型划分法可以直观地判断地下水水化学类型、组分及水质状况，判断地下水水源、含水层间水力联系等^[11-14]；多元数理统计法利用统计学手段从大量水化学数据中获取重要信息，快捷、高效地分析地下含水层水化学组分间的内在规律^[15]；离子比例系数法通过分析部分离子比值呈规律性变化的特点，判断地下水的形成因素^[16-17]；同位素分析方法用于判断地下水的补给来源、成因以及检测地下水污染^[18-19]等。

华北型石炭-二叠纪煤田作为我国主要的原煤产地^[20]，在面临严重水害问题的同时，也面临地下水水质恶化、水量衰减、矿井水污染等多方面问题，为适应国家对煤矿区环境保护的要求，深部地下水水化学特征、富水性特征及成因机制是目前研究和探索的目标^[21]。山东巨野煤田新巨龙煤矿位于山东省菏泽市巨野县境内，地下水TDS在3~10 g/L，硫酸盐浓度高，受地下水影响，矿井水TDS偏高，水处理难度大，随着煤矿开采，地下含水层还将受开采影响，对其水化学特征及成因机制的分析和研究十分有必要。

笔者以新巨龙煤矿地下水为研究对象，对场区地下水水化学历史数据进行统计、取样进行分析测试，综合氢氧同位素和水化学分析测试结果，通过变异系数分析、Piper图解、同位素分析、相关性分析、Gibbs图解、离子比值关系和矿物饱和指数等分析方法，探究煤矿区场地地下水水化学特征及成因机制，以期为该区地下水、矿井水演化和特征污染物运移扩散的研究奠定基础，同时为该煤矿高TDS、高硫酸盐矿井水处理、资源化利用以及矿井水水源判别提供理论依据，可指导煤矿防治水工作。

1   新巨龙煤矿区场地概况

1.1   自然地理

新巨龙煤矿位于菏泽市巨野县龙堌镇境内，属黄河冲洪积平原，地形平坦，地势略呈西北高东南低，地面高程为+40.01~+46.14 m，平均+43.26 m，自然地形坡度为2‰；井口高程+44.8 m。地表水系比较发育，河流沟渠纵横成网，且多系人工开掘的季节性河流，主要有洙水河等，并与区内各沟渠贯通，旱季可引水灌溉，雨季可防洪排涝。

1.2   水文地质

矿井为全隐蔽的华北型石炭-二叠纪煤田，煤系以奥陶系石灰岩为基底，沉积石炭系本溪组、石炭-二叠系太原组，二叠系山西组和石盒子群，上覆地层为新近系、第四系；主要含煤地层为太原组和山西组，主采二叠系山西组3煤，开采深度800~1 300 m^[22]。区内地层大致呈走向南北、倾向东的单斜构造，发育有次一级宽缓褶曲并伴有一定数量的断层，构造复杂程度中等。

含水层自上而下有新生界砂砾层孔隙含水层、二叠系石盒子群砂岩、山西组3煤层顶、底板砂岩(简称3砂)裂隙含水层、石炭-二叠系太原组石灰岩(三灰、十_下灰)及奥陶系石灰岩(简称奥灰)岩溶裂隙含水层(图 1)，主采3煤主要充水水源有底板奥灰水、太原组三灰水、顶板山西组3煤顶底板砂岩水和新近系砂岩水。根据区域规律，未开发前浅层地下水顺地势径流，深层地下水随着补给与排泄区的分布变化而变化，地下水循环慢，径流微弱。受矿区排水影响，各含水层水位将发生分异。各含水层静止水位差别不大，水力联系程度差，多数无直接联系。

图  1  新巨龙煤矿水文地质柱状图

Figure  1.  Hydrogeological histogram of Xinjulong Coal Mine

下载: 全尺寸图片 幻灯片

2   材料与方法

2.1   水质资料来源及水样采集

为分析研究区地下含水层水化学特征及水质变化情况，本文整理了该煤矿场区勘探阶段至建井(1999—2008)，以及正式投产至今(2008—2018)2个阶段的所有水质分析数据，并于2020年10月在新巨龙煤矿目前已开挖的巷道及工作面进行取样。历史和取样测试取得的水质数据共计76个，其中新近系水样15个、山西组3砂水样29个、太原组三灰水样21个、奥灰水样11个。

2.2   样品测试

δ²H和δ¹⁸O同位素测试采用气体稳定同位素比质谱仪MAT253，测定结果以相对于维也纳标准平均海水(VSMOW)的千分差(‰)表示，测定精度分别为±2.0‰和±0.1‰；水质常规阴离子Cl^-、SO₄^2-、NO₃^-、F^-、Br^-、PO₄^3-采用离子色谱(MagIC Net)测定；HCO₃^-、CO₃^-根据MT/T 204—2011《煤矿水碱度的测定方法》使用877 Titino plus滴定仪进行测试；Na、K、Ca、Mg、Fe、Sr等30种金属阳离子使用ICP-700ES电感耦合等离子光谱仪测试。

3   结果与讨论

3.1   煤矿区场地地下水水化学特征

3.1.1   水化学特征

为研究该煤矿区场地地下含水层水化学特征及水质演化情况，根据整理的研究区历史水质资料，结合本次研究的取样测试结果，重点分析对目前正在开采的山西组3煤有影响的新近系底部孔隙(新近系)含水层、3煤顶底板砂岩裂隙(3砂)含水层、太原组第三层灰岩岩溶裂隙(三灰)含水层及奥陶系灰岩(奥灰)岩溶裂隙含水层的水质特征。

通过SPSS软件对水化学各个指标进行数理统计分析，统计结果见表 1；运用AqQA水化学分析软件，分析各含水层中的常规水化学成分，得到地下水的Durov图(图 2)。

表  1  含水层水化学参数统计特征值

Table  1.  Hydrochemical analysis data in groundwater of each aquifer

含水层	指标	各离子质量浓度ρ/(mg·L-1)						TDS/(mg·L-1)	pH值
		Na++K+	Ca2+	Mg2+	Cl-	SO42-	HCO3-
新近系含水层	Min	362.43	288.12	23.53	50.88	2 087.64	5.61	3 190.75	6.50
	Max	1 980.78	499.90	257.56	106.32	5 729.88	196.55	8 429.66	9.90
	Avg	1 723.48	416.81	159.68	85.86	4 991.42	119.24	7 446.65	8.21
	SD	382.97	71.38	57.97	13.54	907.76	48.19	1 315.89	0.71
	CV	0.22	0.17	0.36	0.16	0.18	0.40	0.18	0.09
山西组3砂含水层	Min	1 249.77	22.47	7.34	50.81	2 141.11	141.01	4 045.74	7.30
	Max	2 030.28	506.21	224.43	673.91	5 637.86	1 119.98	8 467.84	8.40
	Avg	1 603.62	216.18	87.09	121.51	3 777.36	325.33	5 972.35	7.82
	SD	223.63	126.95	56.25	114.39	998.27	208.18	1 282.29	0.33
	CV	0.14	0.59	0.65	0.94	0.26	0.64	0.21	0.04
太原组三灰含水层	Min	703.96	19.77	9.33	92.29	845.19	146.92	2 720.44	7.15
	Max	1 489.24	376.34	113.27	294.19	3 397.29	1 239.70	5 383.68	8.40
	Avg	1 184.29	157.18	60.20	198.28	2 500.10	344.41	4 274.74	7.80
	SD	192.63	91.58	26.59	54.79	488.44	244.67	526.90	0.38
	CV	0.16	0.58	0.44	0.28	0.20	0.71	0.12	0.05
奥陶系含水层	Min	422.28	1.60	1.01	132.44	310.35	138.64	1 326.22	7.20
	Max	1 630.91	578.01	154.92	394.17	4 516.08	300.23	6 918.55	11.10
	Avg	762.85	299.52	74.75	313.34	1 980.66	215.71	3 539.01	8.07
	SD	335.26	154.07	46.19	65.64	983.35	46.98	1 304.01	1.31
	CV	0.44	0.51	0.62	0.21	0.50	0.22	0.37	0.16
注：Min为最小值；Max为最大值；Avg为平均值；SD为标准差；CV为变异系数，无量纲。

下载: 导出CSV 

| 显示表格

图  2  新巨龙煤矿地下水Durov图

Figure  2.  Durov diagram of groundwater in Xinjulong Coal Mine

下载: 全尺寸图片 幻灯片

根据SPSS对水化学参数的梳理统计结果，结合各含水层的Durov图(图 2)，从常规水化学指标来看，pH值在6.5~11.1，呈中性略偏碱性，其中奥灰含水层pH值相差较大，最小值为7.2，最大值为11.1，说明其赋存环境可能存在不均的情况。

据表 1可知，矿区各地下含水层水体TDS普遍偏高，均属咸水范畴，含水层由浅至深TDS呈减小趋势。分析和取样测试结果表明：深部奥灰含水层在高pH值的还原环境下，水体中有明显的刺激性气味，SO₄^2-含量少，TDS低；低pH值的氧化环境下，水体中SO₄^2-含量高，TDS高。

除新近系外，各含水层离子浓度变化较大，但总体上，水样中阳离子质量浓度由大到小为Na⁺+K⁺、Ca²⁺、Mg²⁺，阴离子质量浓度表现为${\rm{SO}}_{\rm{4}}^{{\rm{2}} - }$远大于${\rm{HCO}}_{\rm{3}}^ - $和Cl^-，${\rm{HCO}}_{\rm{3}}^ - $大于Cl^-。常规离子变异系数大于0.5的有3砂含水层中的Ca²⁺、Mg²⁺、Cl^-和${\rm{HCO}}_{\rm{3}}^ - $，三灰含水层中的Ca²⁺、${\rm{HCO}}_{\rm{3}}^ - $，奥灰含水层中的Ca²⁺、Mg²⁺和${\rm{SO}}_{\rm{4}}^{{\rm{2}} - }$。说明深部含水层水化学组分受开采扰动影响较大，水化学组分中Na⁺+K⁺的变化不明显，Ca²⁺、Mg²⁺和${\rm{SO}}_{\rm{4}}^{{\rm{2}} - }$的变化较为明显，存在时空分布不均的情况，同时易受水文地质条件改变及水化学环境变化的影响。

3.1.2   水化学类型

为分析该矿区地下水水质的变化情况，以煤矿建井时间为界，根据水质资料绘制各含水层建井前后的地下水水化学Piper三线图，如图 3所示。

图  3  新巨龙煤矿地下水Piper三线图

Figure  3.  Piper diagram of groundwater in Xinjulong Coal Mine

下载: 全尺寸图片 幻灯片

新近系水化学特征较为简单，含水层水样落点较为集中，水化学类型为SO₄-Na型，随着建井开采，水化学类型未发生明显变化；山西组3砂含水层水样落点也相对集中，水化学类型也为SO₄-Na型，随着建井开采，含水层中发生溶滤及离子交换作用，阳离子组分发生变化，阴离子变化则不太明显，水化学类型尚未发生大的变化，仍以SO₄-Na型为主；太原组三灰和奥灰含水层建井后的水化学类型则较建井前发生比较大的变化，建井前水化学类型均为SO₄-Na型；建井后部分三灰含水层水样的HCO₃^-浓度有明显升高，水化学类型变为SO₄-Na和SO₄·HCO₃-Na型，部分奥灰含水层水样中Ca²⁺和Mg²⁺浓度变化较明显，主要的水化学类型为SO₄-Na和SO₄-Ca·Mg型。

总体来看，除新近系外，随着矿井建设和煤矿开采，3砂、三灰和奥灰含水层的水化学类型发生一定程度的变化，且深部含水层的变化更为明显，说明受开采影响，地下含水层的水文地球化学条件发生了改变，含水层水体受到扰动，水化学组分随之发生变化。

3.1.3   同位素特征

根据氢氧稳定同位素特征可以更好地了解地下水的水文地球化学作用，通过分析其组成特征可以研究深层地下水水循环及动态特征。研究区取样测得的氢氧同位素数据见表 2，根据各水样的δ²H、δ¹⁸O值绘制δ²H-δ¹⁸O关系图(图 4)。从图 4可以看出，研究区所取的3个地下水水样：3砂、三灰和奥灰水，在图中的落点位置相对集中，位于大气降水线的右下方，存在¹⁸O正向漂移现象，研究区所取水样的含水层温度较高，在43~54℃，地下水在高温情况下长时间与岩石矿物接触发生¹⁸O交换，其中砂岩含水层中水体主要与硅酸盐矿物发生反应，灰岩含水层中水体主要与碳酸盐岩发生反应，从而导致水体中的¹⁸O含量增加^[23]。

表  2  氢氧同位素数据

Table  2.  Hydrogen and oxygen isotope data

水样	δ2H/‰	δ18O/‰	d=δ2H-8δ18O
雨水	-28.1	-5.09	12.62
01地表排水口	-65.4	-9.35	9.40
02洙水河交汇处	-62.6	-8.64	6.52
03洙水河上游	-54.1	-6.91	1.18
04洙水河下游	-62.8	-8.36	4.08
05水仓	-64.2	-8.60	4.60
06岩巷淋水	-60.1	-8.07	4.46
07煤巷淋水	-58.5	-7.83	4.14
08老空水	-56.2	-7.12	0.76
09奥灰水	-55.4	-6.48	-3.56
10 3砂水	-53.4	-6.40	-2.20
11三灰水	-55.2	-6.32	-4.64
注：d为氘过量参数。

下载: 导出CSV 

| 显示表格

图  4  新巨龙煤矿区地下水δ²H与δ¹⁸O的相关关系

Figure  4.  Relations between δ²H and δ¹⁸O of groundwater in Xinjulong Coal Mine

下载: 全尺寸图片 幻灯片

地下水中，氘过量参数d可以总体反映水岩反应中¹⁸O的交换程度，地下水中¹⁸O的含量取决于岩石的含氧化合物组分、含水层的温度和地下水在含水层内滞留时间的长短。地热水在热储中滞留时间越长，径流速度越慢，其水岩相互作用越明显，d值越小^[24]。表 2中对各水样的d值进行了计算，结果显示只有3砂、三灰和奥灰地下水的d值为负，为-2.20‰~-4.64‰；地表水和矿井水的d均为正值，说明区域地下含水层中，水体有明显的热交换趋势，水循环缓慢，水交替作用弱，随着水体和围岩间相互作用时间越久，与含氧矿物中的¹⁸O交换越多。

3.2   煤矿区场地地下水水化学成因

3.2.1   水化学指标相关性

各含水层的TDS浓度存在明显不同，同时TDS值与各离子浓度的相关性也不尽相同，通过分析各离子与TDS之间的相关性可以解释离子的来源和TDS变化的原因^[25]。利用SPSS软件中双变量相关性分析方法对各含水层地下水的TDS和常规离子进行相关性分析，得到各含水层水化学成分相关系数矩阵，见表 3。

表  3  含水层水化学指标相关系数矩阵

Table  3.  Correlation coefficient matrix of aquifer water chemical composition

含水层	指标	Na++K+	Ca2+	Mg2+	Cl-	SO42-	HCO3-	TDS
新近系含水层	Na++K+	1
	Ca2+	0.011	1
	Mg2+	0.097	0.808**	1
	Cl-	-0.198	0.149	-0.005	1
	SO42-	0.907**	0.395	0.487	-0.100	1
	HCO3-	0.121	0.539*	0.759**	-0.169	0.395	1
	TDS	0.919**	0.375	0.463	-0.109	0.999**	0.384	1
山西组3砂含水层	Na++K+	1
	Ca2+	0.726**	1
	Mg2+	0.610**	0.856**	1
	Cl-	-0.379*	-0.318	0.090	1
	SO42-	0.894**	0.927**	0.748**	-0.485**	1
	HCO3-	-0.213	-0.365	-0.022	0.709**	-0.453*	1
	TDS	0.919**	0.927**	0.824**	-0.325	0.980**	-0.285	1
太原组三灰含水层	Na++K+	1
	Ca2+	-0.604**	1
	Mg2+	-0.496*	0.850**	1
	Cl-	-0.669**	0.443*	0.330	1
	SO42-	0.510*	0.306	0.428	-0.319	1
	HCO3-	0.051	-0.642**	-0.722**	-0.176	-0.778**	1
	TDS	0.652**	0.176	0.280	-0.384	0.976**	-0.639**	1
奥灰含水层	Na++K+	1
	Ca2+	-0.086	1
	Mg2+	0.468	0.828**	1
	Cl-	-0.952**	-0.206	-.690*	1
	SO42-	0.833**	0.475	0.873**	-0.951**	1
	HCO3-	0.361	-0.176	-0.015	-0.300	0.194	1
	TDS	0.845**	0.458	0.864**	-0.957**	1.000**	0.211	1
注：*表示相关性显著水平在0.05(双尾)；**表示相关性显著水平在0.01(双尾)。

下载: 导出CSV 

| 显示表格

从表 3可以看出，新近系水TDS与K⁺+Na⁺和${\rm{SO}}_{\rm{4}}^{{\rm{2}} - }$正相关，3砂水TDS与K⁺+Na⁺、Ca²⁺、Mg²⁺和${\rm{SO}}_{\rm{4}}^{{\rm{2}} - }$正相关，三灰水TDS与K⁺+Na⁺和${\rm{SO}}_{\rm{4}}^{{\rm{2}} - }$正相关，奥灰水TDS与K⁺+Na⁺、Mg²⁺和${\rm{SO}}_{\rm{4}}^{{\rm{2}} - }$正相关；且研究区各含水层TDS和${\rm{SO}}_{\rm{4}}^{{\rm{2}} - }$的相关性接近1，K⁺+Na⁺和${\rm{SO}}_{\rm{4}}^{{\rm{2}} - }$的相关性除三灰外在0.9左右；TDS与${\rm{HCO}}_{\rm{3}}^ - $和Cl^–表现为明显的负相关或非常小的相关性，说明该区域地下水TDS主要受K⁺+Na⁺和${\rm{SO}}_{\rm{4}}^{{\rm{2}} - }$影响，Ca²⁺和Mg²⁺的影响程度次之；含水层存在一定程度的硫酸盐、碳酸盐和岩盐溶滤作用。同时，由于地下水中Na⁺与TDS的相关性明显大于Ca²⁺和Mg²⁺，Na⁺与Ca²⁺和Mg²⁺呈负相关或低相关性，说明还存在阳离子交换作用导致的Ca²⁺、Mg²⁺组分含量减少。Na⁺、${\rm{SO}}_{\rm{4}}^{{\rm{2}} - }$组分含量逐渐升高，成为地下水主要成分，故形成SO₄-Na型水。

3.2.2   水动力因素

Gibbs图可以用来判断地下水化学类型及演化机制的主要影响因素，通过分析ρ(Na⁺)/ρ (Na⁺+Ca²⁺)、ρ(Cl^-)/ρ(Cl^-+HCO₃^-)与TDS的比值关系，可以判断水体水化学成分的变化具体受大气降水作用、岩石风化作用还是蒸发浓缩作用的影响。对于深层地下含水层来说，随着深度增加，蒸发作用的影响逐渐降低，水动力因素成为水化学成分变化的主导因素，据此，H. K. 依格纳托维奇提出3个具有代表性的“水交替带”，即：水积极交替作用带、水迟缓交替作用带和水消极交替作用带^[26]。

将该煤矿区地下水水化学数据投到Gibbs图上(图 5)，地下水中，TDS分布在3 000~9 000 mg/L，ρ(Na⁺)/ρ(Na⁺+Ca²⁺)主要在0.4~1.0，ρ(Cl^-)/ρ(Cl^-+HCO₃^-)主要在0.1~1.0。从图 5a中可以明显看到大部分地下水水样点位于水消极交替作用带范围内，蒸发浓缩作用机制占主导，部分水样点落在分区之外，主要表现为Na⁺浓度偏大，说明对于深层地下水而言，Na⁺还有其他来源，结合前文分析，Na⁺主要受蒸发浓缩作用和岩盐溶解作用影响。图 5b中地下水水样点除奥灰外大部分位于水消极交替作用带范围外，主要是因为除奥灰含水层外，其余含水层水体中Cl^-含量偏低，进一步说明含水层中岩盐溶解Na⁺的贡献程度一般，还受水交替消极作用导致的高温蒸发浓缩作用影响。

图  5  新巨龙煤矿地下水Gibbs图

Figure  5.  Gibbs diagram of groundwater in Xinjulong Coal Mine

下载: 全尺寸图片 幻灯片

结合3.1.3节中同位素特征分析的结果可以看出，该煤矿区地下水动力条件差，蒸发浓缩作用较强，地下含水层中水交替作用消极，长时间的高温水岩及蒸发浓缩作用导致Na⁺、SO₄^2-组分增加，地下水TDS升高^[27]。

3.2.3   水化学组分来源

不同水体中离子浓度的毫克当量γ比值在数值上会存在明显的差异，可以用来判断地下水的形成原因，进一步了解地下水水化学常量组分的来源。γ(Na⁺)/γ(Cl^–)称为地下水的成因系数，表征地下水中Na⁺的富集程度^[28]，自然条件下，岩盐溶解会释放出等量的Na⁺和Cl^–进入溶液^[29]，比值接近1时，表明盐岩矿物的溶滤在地下水中起主导作用，图 6a中γ(Na⁺)/γ(Cl^–)>1，说明Na⁺不仅源于岩盐溶解，还受蒸发浓缩作用影响。γ(Ca²⁺+Mg²⁺)/γ(${\rm{HCO}}_{\rm{3}}^ - $)可判断Ca²⁺、Mg²⁺的主要来源，若比值大于1，说明Ca²⁺、Mg²⁺不仅来源于碳酸岩盐溶解还有其他来源；图 6b中除奥灰含水层外，地下水中Ca²⁺、Mg²⁺除碳酸盐外还有其他来源。γ(Ca²⁺+Mg²⁺)/γ(${\rm{SO}}_{\rm{4}}^{{\rm{2}} - }$+${\rm{HCO}}_{\rm{3}}^ - $)可以进一步判断Ca²⁺、Mg²⁺的主要来源，若比值接近1，说明Ca²⁺、Mg²⁺主要来自碳酸盐和硅酸盐溶解^[30]；图 6c中，γ(Ca²⁺+Mg²⁺)/γ(${\rm{SO}}_{\rm{4}}^{{\rm{2}} - }$+${\rm{HCO}}_{\rm{3}}^ - $) < 1，说明Ca²⁺、Mg²⁺还来源于硫酸盐溶解，${\rm{SO}}_{\rm{4}}^{{\rm{2}} - }$和${\rm{HCO}}_{\rm{3}}^ - $则还受氧化还原作用和CO₂溶解影响。

图  6  新巨龙煤矿地下水离子比值关系

Figure  6.  Diagram of ion raties of groundwater in Xinjulong Coal Mine

下载: 全尺寸图片 幻灯片

地下水成因除溶滤作用外，还有离子的交换作用。确定含水层中是否发生阳离子交换作用可以通过计算氯碱指数CAI-1和CAI-2进行分析。若CAI-1和CAI-2均为正值，进行正向阳离子交换，发生如式(1)的反应；若CAI-1和CAI-2均为负值，则进行反向阳离子交换作用，发生式(2)的反应^[31]。

地下水中阳离子交换作用的程度可以用γ(Na⁺-Cl^-)和γ(Ca²⁺+Mg²⁺-SO₄^2--HCO₃^-)比值关系进行判断，如果反向阳离子交换作用显著，2个参数的比值应接近-1^[32]。

根据氯碱指数图(图 7)中易看出各含水层氯碱指数均为负值，说明地下含水层中存在反向离子交换作用，图 8显示，煤矿区地下水的γ(Ca²⁺+Mg²⁺- SO₄^2--HCO₃^-)和γ(Na⁺-Cl^-)呈负相关性，比值在-0.5左右，结合氯碱指数法得出的结论，说明地下水中存在一定程度的反向阳离子交换作用，但作用程度并不显著，溶滤作用是研究区地下水水化学组分特征的主要成因。

图  7  氯碱指数

Figure  7.  Chloro-Alkaline indices

下载: 全尺寸图片 幻灯片

图  8  新巨龙煤矿γ(Ca²⁺+Mg²⁺-SO₄^2--HCO₃^-)和γ(Na⁺-Cl^-)比值关系

Figure  8.  Diagram of γ(Ca²⁺+Mg²⁺-SO₄^2--HCO₃^-)/γ(Na⁺-Cl^-) of groundwater in Xinjulong Coal Mine

下载: 全尺寸图片 幻灯片

3.2.4   矿物溶解平衡

矿物饱和指数可以用于评估矿物与地下水之间的平衡和反应^[33]，根据矿物相对地下水饱和指数的计算结果，可以确定地下水系统中的反应矿物。水溶液相对矿物的饱和指数SI计算公式为：

式中：IAP为矿物所含组分在水溶液中的离子活度积；K为矿物溶解反应时的平衡常数。

当SI < 0时，表示该矿物相对于水溶液未达到饱和，为溶解状态；

当SI=0时，表示水溶液与矿物正好处于平衡状态；

当SI>0时，表示该矿物相对水溶液处于过饱和状态。

利用PHREEQC软件计算出该区各含水层方解石、白云石、石膏和盐岩的饱和指数，如图 9所示。

图  9  新巨龙煤矿地下水SI与TDS关系

Figure  9.  Correlation between SI and TDS of groundwater in Xinjulong Coal Mine

下载: 全尺寸图片 幻灯片

图 9显示该区方解石和白云石的饱和指数均大于0，处于饱和状态；石膏的饱和指数小于0，除三灰外，饱和指数接近0，虽处于未饱和状态，但有达到饱和的趋势；盐岩的饱和指数在-5.0左右，处于未饱和状态。从石膏和盐岩处于未饱和状态可以说明，地下水中石膏和盐岩的溶解占离子作用的主导地位，也是导致该地区SO₄-Na型水呈现出高TDS的原因。

饱和指数法的分析结果和前文离子相关性和离子比值分析的结果一致，该煤矿区地下水最主要组分的Na⁺和SO₄^2-主要来源于岩盐和石膏矿物的溶解，方解石和白云岩处于饱和状态呈现沉淀趋势，结合前文分析，钙、镁类碳酸盐岩虽存在一定量的溶解，但同时随着蒸发浓缩作用和反向阳离子交换作用，水体中Na⁺不断增加，Ca²⁺和Mg²⁺相应减少，地下含水层中TDS升高。

4   结论

a. 新巨龙煤矿区地下水均为弱碱性水，TDS普遍偏高，深部含水层地下水TDS小于浅部，水化学类型以SO₄-Na型为主，随着煤矿开采，水化学环境发生变化，部分组分离子变异系数较大，新近系和3砂水水化学类型未发生明显变化，三灰和奥灰水的水化学类型由单一的SO₄-Na型增加了SO₄·HCO₃-Na和SO₄-Ca·Mg型，同时由于深层地下水温度偏高，含水层水体中存在¹⁸O的正向漂移。

b. 煤矿区地下水水动力条件差，含水层温度高，长时间高温水岩作用、蒸发浓缩作用、岩盐和硫酸盐矿物溶滤作用以及一定程度的阳离子交换作用，是该区高硫酸盐、高TDS的SO₄-Na型地下水的主要成因。

c. 对煤矿区地下水水化学特征及成因机制的研究可为该煤矿高TDS、高硫酸盐矿井水处理、资源化利用以及污染防控提供重要支撑；同时可为后续矿井水水源判别提供理论依据，指导煤矿防治水工作的进行。