研究区位于安徽省宿州市埇桥区朱仙庄镇东侧采煤沉陷区内。朱仙庄矿位于市东南13 km处，朱仙庄镇，北起植物园南边缘(界洪河)，穿303省道向东南延伸，地表为平原，宿灵公路贯穿矿区，专用运煤铁路与京沪铁路连接，交通运输便利。朱仙庄煤矿于1983年建成投产，煤种为三分之一焦煤和气煤，产能为120万t/年，矿区井田面积约21.56 km²，是淮北煤田的重要组成部分，位于宿东向斜北翼，轴部为二叠系，两翼由石炭系和二叠系组成，总体呈SN向。区内地层有奥陶系、石炭系、二叠系、侏罗系、新近系和第四系。朱仙庄煤矿含水层由新生界松散含水层、侏罗系含水层、二叠纪砂岩含水层、太灰含水层以及奥灰含水层组成，由于第四纪含水层的特殊位置，造成其是矿井水的主要来源，也是塌陷区域水来源。目前，宿州朱仙庄煤矿塌陷塘面积约为23.36 km。研究区内居民的经济来源主要以种植业为主，并在塌陷塘内有水产养殖活动进行。

2023年12月在宿州市埇桥区东部采煤塌陷塘塘水及附近村庄与农田地下水进行系统水样采集工作，共采集水样36件，采样点分布见 图1 。在采集地下水水样时，先用预采集水样对事先用纯水清洗过的样品瓶润洗3次后装瓶至溢出，贴好标签密封保存。使用GPS工具箱(中国小狼信息技术有限公司)对采样点经纬度信息进行现场记录。利用便携式测试仪，对pH、TDS (溶解性总固体)、ORP (氧化还原电位)进行现场测定、记录。其余水化学指标在安徽省矿井水资源化利用重点实验室测试。采集样品经0.45 μm滤膜抽滤后低温保存，随后进行主离子测试。采用美国赛默飞公司ICS-600型和ICS-900型离子色谱仪分别对阳离子(Na⁺、K⁺、Ca²⁺、Mg²⁺、Li⁺)、阴离子( ${\text{SO}}_4^{2-}$、Cl⁻、F⁻、 ${\text{NO}}_3^{-}$)检测分析， ${\text{HCO}}_3^{-}$和 ${\text{CO}}_3^{2-}$采用滴定法测定 [7] 。

对研究区采集的36件水样水化学指标进行统计分析结果见 表1 ，研究区地下水中阳离子主要以Na⁺为主，其离子质量浓度占阳离子总量的58.09%，平均浓度为131.34 mg/L，阳离子平均质量浓度呈现Na⁺ > Ca²⁺ > Mg²⁺ > K⁺的关系；阴离子主要以 ${\text{HCO}}_3^{-}$为主，其离子质量浓度占阴离子总量的60.66%，平均浓度为346.63 mg/L，阴离子呈现 ${\text{HCO}}_3^{-}$ > ${\text{SO}}_4^{2-}$ > Cl⁻ > ${\text{NO}}_3^{-}$ > F⁻的关系。地下水pH值范围为7.32~8.67，平均值为8.02，整体呈弱碱性。36件样品中有4件样品TDS值超出仪器测量范围，剩余32件样品的TDS含量介于451.00~978.00 mg/L之间，平均值为529.64 mg/L，属于淡水类型 [8] 。氧化还原电位(ORP)为94.00~187.00 mV，平均值为136.81 mV，适于硝化反应的进行 [9] ，有助于净化水质，减少水中的有害物质。

Piper三线图是一种常用于反映水化学主离子的相对丰度和分布特征，分析水化学成分演化规律的水化学类型图示方法 [10] 。基于地下水测试结果，绘制地下水Piper图，如 图2 所示。研究区地下水阳离子三角图中，水样点沿Na⁺-K⁺边线分布，Na⁺为优势阳离子；阴离子主要沿阴离子三角图 ${\text{HCO}}_3^{-}$ + ${\text{CO}}_3^{-}$-Cl⁻边线分布，表明阴离子以 ${\text{HCO}}_3^{-}$和 ${\text{SO}}_4^{2-}$为主， ${\text{HCO}}_3^{-}$占主导优势。根据舒卡列夫分类，研究区地下水样品水化学类型大部分为Na-HCO₃型，共21件，占总水样的58.33%，部分为Ca-HCO₃型，共8件，占总水样的22.22%，少量为Na-SO₄型，共3件，极少量为Na-Cl型和Mg-HCO₃型，皆为2件。

Gibbs图解法是揭示研究区水体来源机制分析与水化学形成作用的重要方法之一 [11] 。Gibbs [12] 通过对世界多数地下水的TDS与Na⁺/(Na⁺ + Ca²⁺)及Cl⁻/(Cl⁻ + ${\text{HCO}}_3^{-}$)关系研究，总结出影响地下水化学组成的三大控制因素：蒸发浓缩作用控制型、岩石风化作用控制型、降水作用控制型。其中，溶解性物质含量较低且Na⁺/(Na⁺ + Ca²⁺)的比值或Cl⁻/(Cl⁻ + ${\text{HCO}}_3^{-}$)的比值接近1时，具有该特征的点分布于Gibbs模型的右下角，表明此种地下水离子来源以降水控制为主；溶解性物质含量中等且Na ⁺/(Na⁺ + Ca²⁺)的比值或Cl⁻/(Cl⁻ + ${\text{HCO}}_3^{-}$)的比值小于0.5或位于0.5附近时，该地下水的水化学组分以岩石风化作用控制为主导，此种特征点分布于Gibbs模型中部靠左；溶解性物质含量较高且Na⁺/(Na⁺ + Ca²⁺)的比值或Cl⁻/(Cl⁻ + ${\text{HCO}}_3^{-}$)的比值接近1时，该地下水的特征点位于模型的右上角，反映该地下水离子来源受蒸发–结晶控制为主。研究区溶解性物质质量浓度中等，Na⁺/(Na⁺ + Ca²⁺)的比值在0.93~0.22，Cl⁻/(Cl⁻ + ${\text{HCO}}_3^{-}$)的比值在0.02~0.47。如 图3 所示，研究区地下水化学组成主要受岩石风化作用控制。

地下水水化学成分主要来源于岩石的风化溶解，Ca²⁺、Mg²⁺来源以碳酸盐岩、硅酸盐岩和蒸发岩的溶解为主，Na⁺、K⁺来源以蒸发岩和硅酸盐岩的溶解为主， ${\text{SO}}_4^{2-}$、Cl⁻来源以蒸发岩的溶解为主， ${\text{HCO}}_3^{-}$来源以碳酸盐岩和硅酸盐岩溶解为主，通常使用混合图揭示地下水水域中，化学风化作用产生的离子来源 [13] 。由 图4 研究区水样中Ca²⁺/Na⁺与 ${\text{HCO}}_3^{-}$/Na⁺、Mg²⁺/Na⁺的离子浓度比值主要分布在硅酸盐岩附近，表明其离子来源受硅酸盐岩、碳酸盐岩、硫酸盐岩的风化溶解作用影响，以硅酸盐岩的风化溶解作用影响为主。

综合水质指数(WQI)被广泛运用于地表水、地下水水质评价工作中，根据世界卫生组织(WHO)水质参数标准 [14] ，见 表2 ，采用加权算术指数法计算每个样品的WQI，本研究选取Na⁺、K⁺、Mg²⁺、Ca²⁺、Cl⁻、 ${\text{SO}}_4^{2-}$、 ${\text{HCO}}_3^{-}$、pH、TDS、 ${\text{NO}}_3^{-}$，10个指标进行地下水WQI饮用评价，其计算过程见文献 [15] 。基于WQI可将水质情况分为5类：WQI < 25，表明水质极好；26 < WQI < 50，表明水质好；51 < WQI < 75，表明水质差；76 < WQI < 100，表明水质极差；WQI > 100，表明水质不适宜饮用。由 图5 可知，研究区地下水样本中有68.75%分布于26 < WQI < 50区间内，表明水质较好；剩余水样WQI值均大于51，最大值达到164.83，说明样品水质差或极差，不适合饮用。

灌溉水中化学组分较高的水体进行灌溉时会导致土壤盐碱化，进而影响植物生长 [3] 。对此，选取EC (电导率)、Na% (钠百分比)作为研究区灌溉水质评价参数，用Wilcox图解法绘制灌溉水分类图( 图6 )对研究区地下水进行灌溉水质评价。结果显示， 图6 中水样位于水质过滤类和水质良好类两个区域，表明研究区地下水可直接作为灌溉水使用，且合理利用不会使土壤盐碱化。