1. 引言

水是地球上珍贵的自然资源之一，它承载着生物的生命活动，同时也是各种产业不可或缺的要素。不同类型的水资源，如工业用水、灌溉用水、水产养殖用水等都在各自的领域中发挥着重要的作用。然而，近年来随着社会的快速发展，人类对水资源的大规模开采和利用产生了双重影响。一方面，水资源为居民生活和社会发展提供了有力保障，是确保健康生活的核心要素之一；另一方面，它也带来了一些不利因素，成为了影响安全生产的隐患。

随着水资源利用目标的不同，许多领域的研究都在不断深入。这些领域包括水资源管理、水污染防治、水文地球化学、水质监测等。在这些研究中，水污染问题逐渐被揭露出来，它与农药、化肥的使用以及工业废水排放等因素密切相关。这些问题不仅对水资源本身造成污染，还会对整个生态环境产生深远的影响。地下水的水化学特性及其组成是地下水与其周围环境长期相互作用的结果。这些特性可以反映特定区域的演化历史，而水化学信息也为研究者提供了线索，帮助其了解地下水的起源和形成过程。通过使用Piper图、Stiff图、平均当量变化雷达图、描述性统计、重金属污染指数等方法，我们可以有效地解释地下水的起源、演化和污染程度等，这些方法有助于相关管理和决策部门更好地管理和保护地下水资源 [1] [2] 。

煤炭是世界上储量最丰富的化石能源，全球约有80多个国家拥有丰富的煤炭资源。在中国，矿井突水是六大矿难(水灾、火灾、瓦斯爆炸、粉尘危害、地热灾害和顶板垮塌)中最为严重的灾害之一。矿井突水是指掘进或采矿过程中，巷道揭穿导水断裂、富水溶洞、积水老空，大量地下水突然涌入矿山井巷的现象。矿井突水是煤矿生产过程最具威胁的灾害之一，人员伤亡大，经济损失列于煤矿三大事故的榜首。一旦矿井发生突水，如何及时准确地判断突水成因，查找突水水源，是解决和进一步预防突水灾害的关键问题。通过分析离子浓度，我们可以更好地了解含水层内隐藏的水化学信息。因此，将丰富的水化学数据与水文地球化学分析方法相结合，可帮助我们更有效地识别突水含水层水化学特征及突水水源，精准的水源判别分析结果将为实施有效的突水防治措施提供重要的参考依据。

两淮煤田煤炭资源丰富，开采出的煤炭资源是华东地区经济持续发展的重要能源保障。青东煤矿是淮北煤田的井工煤矿之一，随着煤炭资源的开采，矿井地质条件水文地质条件发生了一定的变化，受采掘作业的扰动下，其含隔水层组性能及水流场、水化学场等多场环境也发生了明显变化。含水层(组)流场及水化学场变化给矿山的安全生产带来了挑战，因此深入研究和理解矿井含水层(组)的特性对于制定有效的突水控制措施至关重要 [3] 。近年来，研究区及周边相继开展了土壤、塌陷塘、河流、地下水等相关研究。孙浩运用地累积指数法、内梅罗综合指数法和潜在生态风险评价法等多种方法评价了淮北矿区土壤的生态环境质量，揭示了微量元素在固废堆积区周边土壤中的分布规律，探讨了固废类型和堆积时间对土壤中微量元素分布的影响 [4] ；李涛通过对淮北地区塌陷区水域渔业开发及利用现状的分析，探讨了渔业可持续发展的对策 [5] ；李兴旺等人采用有限单元法、面状井系法、水均衡法、系统分析等方法，研究地表水与地下水联合运用，建立了亳州涡河以北地下水位预测模型 [6] ；潘强通过对不同土壤侵蚀类型区的各项水土流失影响因子统计、归纳，来研究调查水土流失现状、水土流失规律、水土保持措施及其生态建设模式，并提出平原地区水土保持生态建设模式 [7] 。陈家玉等人利用16SrRNA基因V4区测序方法开展深层地下水微生物群落多样性测试与分析，结果显示煤系水、太灰水是淮北煤田煤矿开采的主要充水水源，而且具有不同的微生物群落特征 [8] ；张杰等人采用水化学及同位素测试并收集相关成果进行数据补充研究华北型煤田奥陶系岩溶水水文地球化学特征及其对地热的指示意义，结果表明淮北煤田岩溶地热水主要来源于大气降水补给，北部闸河矿区以直接入渗补给为主，水化学表现为HCO₃-Ca型，水–岩作用以溶滤溶解为主 [9] 。煤矿砂岩含水层作为矿井充水水源之一，虽有相关矿区主要充水含水层的相关报到，但针对青东矿砂岩含水层的研究较少。基于此，本文以青东矿主采煤层的直接充水水源–砂岩含水层常规离子浓度作为研究对象，讨论了该含水层的水化学类型和水岩相互作用，以期为矿井突水水源的快速识别提供给水文地球化学判别依据。

2. 材料与方法

2.1. 研究区概况

青东煤矿位于淮北煤田(图1)，气候属于暖温带半湿润季风气候区。该区域气候温和，年气温约为14.5℃，降水量达到254.4毫米。研究区位于东经116˚25′44″~116˚34′44″，北纬33˚36′34″~33˚40′29″的范围内。该煤矿地层形成的顺序从古老的奥陶纪开始，历经石炭纪、二叠系，直至第四纪 [10] 。其中，二叠系下石盒子组5号煤层和山西组6号煤层均形成于二叠系，是该矿区最主要的可开采煤层。该层位的煤层不仅具有重要的工业价值，也为当地经济发展提供了有利保障。

青东煤矿矿井突水的四种主要类型，即1) 第四系松散含水层：该含水层松散且易于渗透，含水量丰富，在浅部煤层开采时因覆岩破坏或岩层岩性变化导致突水；2) 5号煤顶板砂岩含水层：该含水层位于5号煤层的顶部，因其形成于陆源环境，其含水量以静储量为主；3) 6号煤底板灰岩含水层：该含水层位于6号煤层的底部，因沉积环境为浅海环境，故其含水量以动储量为主；4) 采空区积水：当工作面回采完成后，开采后形成的空间易形成高位积水区，当老空水未及时疏放时，可对接替采区或工作面形成安全隐患。

2.2. 采样与分析

本次主要研究煤系砂岩含水层，共收集23个样品，样品采集均来井下出露点。现场采集的水用聚乙烯瓶进行密封，随后运至指定的测试中心(具备CMA实验室资质)进行常规离子测试工作，各常规离子的测试方法参照相关规范进行 [11] 。离子平衡误差精度控制在±5%以内。

2.3. 描述性统计

本文拟采用多元统计分析软件SPSS，进行矿区砂岩水的描述性统计、主成分分析和相关性分析。描述性统计可用于讨论常规离子含量的变化；主成分分析用于提取出控制含水层水化学特征的影响因素；相关性分析用来检验变量之间的相互关联性，以皮尔逊系数来反映变量之间的关系程度 [12] 。

2.4. 水化学分析

采用Origin软件，绘制Piper三线图，对砂岩含水层的水化学类型进行分析；利用离子含量图对砂岩含水层的常规离子含量进行对比分析。

3. 结果与讨论

3.1. 水化学类型分析

Piper三线图是研究地下水化学类型常用的方法。如图2所示，阳离子三角图显示Na⁺的含量均高于其他离子，阴离子三角图显示 $S{O}_4^{2-}$ 、HCO₃含量最大。水化学类型中，20个样品为Na-HCO₃型、3个样品为Na-SO₄型 [13] 。从水化学类型分析可得，研究区水岩相互作用主要以碳酸盐岩、硫酸盐岩及硅酸盐岩矿物的风化溶解为主。

由图3研究区离子浓度含量可知，采样点离子含量排序： $S{O}_4^{2-}$ > Na⁺ + K⁺ > $HC{O}_3^{-}$ > Cl⁻ > Ca²⁺ > $C{O}_3^{2-}$ > Mg²⁺。阳离子Na⁺ + K⁺含量明显高于其他阳离子含量，阴离子以 $S{O}_4^{2-}$ ， $HC{O}_3^{-}$ 以及Cl⁻为主要阴离子。在阳离子中，Na⁺ + K⁺的含量占比大，占总体阳离子含量的92%，主要来源于盐岩矿物溶解或钠长石等硅酸岩矿物风化；Ca²⁺和Mg²⁺的含量占比相同，相对于总体阳离子含量占比较低，主要来源于碳酸岩矿物风化；在阴离子中， $S{O}_4^{2-}$ 和 $HC{O}_3^{-}$ 的含量相对于Cl⁻占比大，主要来源于硫酸盐矿物和盐岩矿物溶解； $C{O}_3^{2-}$ 的含量占比低于其他阴离子含量占比，主要来源于碳酸岩矿物或硅酸岩矿物风化溶解。

3.2. 主成分分析

经过旋转后，选取特征值大于1的主成分，得到3个主成分组成(表1)。结果表明：3个主成分方差累计为80.757%，主成分1的方差贡献率为28.35%，总计为2.268，其中主成分1中Ca²⁺和Cl⁻均有较强的正相关，Cl⁻含量与岩盐矿物溶解有关，而Ca²⁺含量较高则说明阳离子交替吸附作用也较为明显 [14] 。主成分2的方差贡献率为27.45%，总计为2.196，其中Ca²⁺、Mg²⁺荷载值较高，其次为 $S{O}_4^{2-}$ ，说明主成分2的水化学作用主要为硫酸盐岩矿物溶解 [15] 。主成分3的方差贡献率为24.957%，总计为1.997，其中K⁺ + Na⁺、 $S{O}_4^{2-}$ 有强正相关，与Cl⁻有中相关，代表了硫酸岩矿物、岩盐等矿物的溶解 [16] 。一般情况下，碳酸盐岩矿物的溶解也会在水岩相互作用中发生一定程度的作用，但从主成分分析未发现明显的作用，说明该含水层中的碳酸盐岩矿物风化溶解作用强度较低。

3.3. 相关性分析

相关性分析可指示出两个参数之间的相关性，从而判定两个参数是否属于同源物质，从而更好地了解水质的影响因素，是一种研究事物之间相互依存关系的数理统计方法。由表2分析可知，pH和 $C{O}_3^{2-}$ 在α = 0.01相关性显著，相关系数为0.713，说明它们有极高的相关性，表明研究砂岩含水层中 $C{O}_3^{2-}$ 对pH起主要影响作用 [17] 。在显著水平0.01上，K⁺ + Na⁺和 $S{O}_4^{2-}$ (0.774)、Ca²⁺和Mg²⁺ (0.812)之间的正相关性较强，主要源于硅酸盐矿物的风化和硫酸盐岩矿物的溶解 [18] 。K⁺ + Na⁺和Cl⁻ (0.441)呈中等正相关，说明这些元素可能来自氯酸盐和岩盐。在显著水平0.05上，pH和Ca²⁺ (−0.465)呈负相关，说明上述离子可能与 $C{O}_3^{2-}$ 发生了交换吸附 [19] 。 $S{O}_4^{2-}$ 和 $HC{O}_3^{-}$ (−0.500)呈负相关 [20] 。

注：^*在0.05级别(双尾)，相关性显著，^**在0.01级别(双尾)，相关性显著。

4. 结论

本文采用多元统计分析(包括描述性统计分析、相关性分析和主成分分析)和水化学类型分析，讨论了矿区砂岩含水层的水文地球化学特征及水岩相互作用类型。研究得出结论如下：

1) 水化学类型分析中，20个样品为Na-HCO₃型、3个样品为Na-SO₄型。

2) 砂岩含水层中的Na⁺ + K⁺来自钠长石等硅酸盐矿物的风化， $S{O}_4^{2-}$ ， $HC{O}_3^{-}$ 以及Cl⁻为主要阴离子，来自石膏、盐岩等硫酸盐矿物以及碳酸盐岩矿物的溶解。

3) 主成分分析提取出三个主成分，主成分1的方差贡献率为28.35%，总计为2.268；主成分2的方差贡献率为27.45%，总计为2.196；主成分3的方差贡献率为24.957%，总计为1.997。成分1与岩盐矿物溶解和阳离子交替吸附作用有关；主成分2与硫酸盐岩矿物溶解有关，主成分3与硫酸岩矿物、岩盐等矿物的溶解有关。

4) 相关性分析中， $C{O}_3^{2-}$ 对pH起主要影响作用，K⁺ + Na⁺与 $S{O}_4^{2-}$ 、Ca²⁺，以及Mg²⁺之间的正相关性较强，主要源于硅酸盐矿物的风化和硫酸盐岩矿物的溶解；在显著水平0.05上，pH和Ca²⁺ (−0.465)呈负相关，说明上述离子可能与 $C{O}_3^{2-}$ 发生了交换吸附。

基金项目

项目资助：宿州学院大学生创新创业训练计划项目，YBXM23-022；宿州学院课程思政建设研究项目，szxy2021ksjy031；宿州学院新建专业质量提升项目(szxy2023xjzy02)。

参考文献

NOTES

^*通讯作者。

参考文献