氟是最轻的卤族元素之一,水溶液中常常以氟化物(F-)的形式存在。地下水中F-的高低对人类健康都有重要的影响,摄入不足或过量均会导致健康问题[1]。摄入低质量浓度的F-(0.05~0.07 mg/L)对牙齿有好处,可以帮助硬化牙釉质,有效防止蛀牙的发生[2-3];摄入过量的F-(>4.00 mg/L)会引发神经系统损伤、生育率下降,以及氟骨症等健康问题[4-6]。饮用一直被认为是人体摄入F-最主要的途径[6],为此,世界卫生组织(WHO)建议饮用水中F-质量浓度不得高于1.50 mg/L,而中国《生活饮用水卫生标准》(GB5749—2006)规定,饮用水中F-的质量浓度限值为1.00 mg/L[4,7-8]。为此,通常将水体中F-质量浓度超过1.00 mg/L的统称为高氟水,据不完全统计,中国约有4 000万人长期饮用高氟地下水,主要分布在内蒙古河套平原[9-11],甘肃张掖[1],山西运城[12],山西大同[13],陕西榆林[14],河北沧州[15],安徽皖北平原[16]等干旱半干旱地区。
F-在地下水的质量浓度高低主要受含氟矿物的溶解度、温度、pH、阳离子交换强度、蒸发强度、地下水停留时间、盐度等因素的影响[7-8]。地下水中F-的来源一般分为天然源和人为源2种。天然源主要来自萤石、氟磷灰石、黄玉、云母和闪长石等含氟矿物在水中的溶解[5,17],除了受溶解度控制之外,蒸发浓缩,离子交换和竞争吸附等地质作用也是控制F-在地下水中质量浓度的重要因素[10,18]。人为源有很多种,主要包括密集的农业活动、无计划的城市开发、高强度的工业生产和采矿活动等人为释放的F-[3,12,19-20]。
神东矿区煤炭资源十分丰富,但淡水资源相对匮乏,生态环境脆弱,人均水资源量非常有限。由于地表水资源缺失,采矿产生的矿井水会全部被储存于地下水水库中,作为当地生活饮用水循环利用[21-22]。然而,随着煤矿的高强度大规模开采,高质量浓度的F-在矿井水中被检出,其中布尔台矿区矿井水中F-质量浓度高达7.40 mg/L[23],超过了《生活饮用水卫生标准》7.4倍,高氟矿井水已成为威胁矿区居民饮用水安全和制约资源化循环利用的关键问题。然而其分布特征,来源和形成机制还未进行过系统的分析和研究。为此,笔者以神东矿区矿井水中F-为主要研究对象,系统分析高氟矿井水的含量特征和空间分布,探讨了其来源和形成机制,为实现矿区水资源循环利用,保护矿区居民饮用水水质安全,防止地下水F-污染,维护区域生态环境健康发展提供技术指导和科学支撑。
神东矿区位于陕北高原北缘及鄂尔多斯高原东南部,地处黄土高原北缘与毛乌素沙漠过渡地带东段,处于东经110.07°~110.18°,北纬39.33°~39.50°,总面积约3 481 km2。神东矿区西北部为库布其沙漠,中部为群湖高平原,西南部为毛乌素沙漠,地形总体呈西北高,东南低,海拔平均 800~1 385 m。神东矿区属温带大陆性干旱半干旱气候,多年平均降水量约437.20 mm,多集中在7—9月,多年平均蒸发量2 065.10 mm,多年平均气温9.91 ℃。水系属于黄河水系,主要河流是乌兰木伦河,年平均水流量4.52~12.20 m3/s。
神东矿区主采煤层为侏罗系中下统延安组,煤层赋存稳定,呈近水平,可采储量约为78.38亿t,目前主采煤层为2-2煤,3-1煤和5-2煤,煤层平均厚度4~6 m。神东矿区位于鄂尔多斯盆地之次级构造单元陕北斜坡,矿区构造总体为向西倾斜的单斜地层。主要含水层组自上到下依次为萨拉乌苏组孔隙含水层(Q3s),直罗组(含安定组)裂隙含水层(J2z)和延安组裂隙含水层(J2y)等3组,水文地质剖面如图1所示。因富水性较强,水质较好,Q3s和J2y含水层一直是矿区居民最重要的地下饮用含水层,也是矿井水的主要充水水源。
图1 水体采样点分布和水文地质剖面[24]
Fig.1 Locating of the sampling sites in the Shendong Coal-mining district and Hydrogeological profile[24]
2020年7—8月,在神东矿区共采集了不同类型的水样品62组,其中矿井水45组,延安组(延安组裂隙含水层)6组,萨拉乌苏组(萨拉乌苏组孔隙含水层)4组,地表水(乌兰木伦河水)4组和直罗组(直罗组(含安定组)裂隙含水层)3组,采样点位置如图1所示,图1中C(F-)为F-质量浓度。为获取新鲜的代表性水样品,采样前必须保证水是流动的,若静置,需排水5~10 min后再采。采样前,先用蒸馏水润洗取样瓶(聚乙烯棕色瓶)2~3遍,再用待取水样润洗2~3遍后取样。每个采样点,分别采集3瓶(每瓶500 mL)的水样,其中1瓶用于测定阴离子,1瓶用于测定阳离子,1瓶用于测定氢氧同位素。为了提高F-质量浓度的检测精度,每瓶阳离子待测水样中加入5 mL F-标准溶液(1 mg/L)用于加标回收。采样过程中,现场使用便携式pH仪(OHAUS ST20)和便携式TDS仪(OHAUS ST20T-B)测量矿井水的pH和总溶解固体(TDS)值。
测试前,采集的所有水样均需通过0.45 μm的水膜过滤。阴离子和阳离子(Ca2+,Mg2+,Na+,K+)选用离子色谱法(Dionex Integrion IC,Thermo Fisher,USA)进行测定,测量精度为0.10 通过酸碱滴定法测定。水体的氢氧同位素δD和δ18O送往中国地质大学(武汉)生物地质环境地质国家重点实验室,选用同位素比质谱仪(Isotope Ratio Mass Spectrometer,IRMS)进行测试,测量精度分为±1‰和±0.10‰。
为了方便研究,按照F-质量浓度的高低将神东矿区矿井水分为4类:第1类F-质量浓度小于1 mg/L,称为低氟矿井水;第2类F-质量浓度在1~4 mg/L;第3类F-质量浓度在4~10 mg/L;第4类F-质量浓度大于10 mg/L,后3类统称为高氟矿井水。
每个待测水样均采用3次重复测量取平均值的模式,标准偏差控制在10%以下。F-标准溶液的加标回收率控制在95%~105%。最后,所有离子的结果需经过离子平衡误差计算,要求误差不得高于5%。
神东矿区不同水体水化学特征分析结果见表1。
表1 神东矿区不同水体水化学数据
Table 1 Geochemistry data in different water of Shendong mine water samples
类别指标质量浓度/(mg·L-1)Na+K+Ca2+Mg2+Cl-SO2-4HCO-3F-NO-3TDSpHδD/‰δ18O/‰低氟矿井水(n′=10)最大值989.007.5484.9031.60483.00904.001 031.000.958.152 967.008.24-68.00-8.30最小值8.940.6314.847.285.676.57159.000.160.21186.006.20-72.00-9.20平均值211.063.8256.0714.2267.93229.94393.650.412.99789.207.79-70.00-8.75标准差295.142.8819.677.36146.90270.85263.860.272.79827.460.612.830.64高氟矿井水(n′=35)最大值956.0012.0078.2027.13369.001 114.001 524.0012.757.752 648.008.50-65.00-8.40最小值159.000.110.780.6116.7021.20203.320.160480.006.25-86.00-10.60平均值530.794.4616.647.13199.06284.52798.215.011.191 455.407.97-73.00-9.26标准差225.373.2720.247.91131.82252.34360.633.761.65557.990.545.770.55延安组(n′=6)最大值1 245.005.3653.3023.701 069.00988.00497.907.981.333 408.009.00-59.00-7.70最小值100.800.911.090.485.523.65131.000.410195.007.12-75.00-9.40平均值681.443.0225.686.62266.27208.14232.485.190.761 059.678.23-67.50-8.53标准差452.812.0222.779.32408.07384.03139.715.120.551 175.830.675.890.59萨拉乌苏组(n′=4)最大值439.002.6443.6020.07302.00296.25529.001.3322.541 017.008.34-58.00-7.30最小值22.000.834.780.526.964.58196.000.430209.007.00-62.00-7.70平均值196.751.9122.579.33111.78120.28352.200.955.92640.007.79-60.50-7.43标准差179.540.8616.579.01135.08135.86137.100.4111.08334.360.621.730.19直罗组(n′=3)最大值71.481.9369.6024.2024.2228.40418.340.641.80638.5010.01-61.00-7.20最小值13.871.531.290.106.6219.7618.000.480.20113.008.00-61.00-7.30平均值47.051.7638.3011.8215.9523.49207.020.551.07345.028.73-61.00-7.25标准差29.790.2134.5112.068.844.44201.100.080.81268.091.180.000.07地表水(n′=4)最大值326.0019.1067.3032.40244.00137.00445.009.942.33756.009.07-67.00-8.10最小值131.004.663.360.32124.0038.90254.000.830.35598.008.20-69.00-8.40平均值227.758.4730.6115.70162.0084.13306.255.031.31687.508.47-68.00-8.25标准差101.957.1029.8017.0155.3450.7092.634.750.8474.430.412.830.16
注:n′为各类水的样品个数;低于检出限的元素浓度按“0”标记。
由表1可知,神东矿区97.14%的高氟矿井水和90.00%低氟矿井水样品呈现弱碱性或碱性,pH平均值分别为7.97±0.54和7.79±0.61。低氟矿井水TDS值为186.00~2 967.00 mg/L;高氟矿井水的TDS范围较窄,为480.00~2 648.00 mg/L。根据《生活饮用水卫生标准》中TDS限值(1 000 mg/L),高氟和低氟矿井水分别有77.14%和10.00%的样品超标。高氟和低氟矿井水中,Na++K+和Ca2+均是最主要的阳离子,其次为Mg2+。低氟矿井水中阴离子质量浓度值大小顺序为(229.94±270.85 mg/L)>Cl-(67.93±146.90 mg/L)(2.99±2.79 mg/L)。虽然高氟矿井水阴离子质量浓度普遍偏高于低氟地下水,但大小顺序未发生变化,依然为(284.52±252.34 mg/L)>Cl-(199.06±131.82 mg/L)
Piper三线图(图2)中,矿井水样品的阳离子以Na+和K+为主,阳离子质量分数大于40%的占77.78%,阴离子变化范围较大,总体上以和Cl-为主要阴离子,质量分数分别为3.20%~43.88%和30.15%~92.07%,水化学类型主要为Na-HCO3-Cl型水(37.78%),Na-HCO3-SO4型水(20.00%),Ca-Mg-HCO3型(15.56%)等为主(括号中数字为该类型水占所有类型水的比例)。
图2 不同水体Piper三线图
Fig.2 Piper plots of different water of the Shendong mine samples
高氟和低氟矿井水水化学类型差异性较大,低氟矿井水的水化学类型,主要为Ca-Mg-HCO3(60.00%)和Na-HCO3型(40.00%);高氟水的水化学类型相对复杂,主要为Na-HCO3-Cl型水(54.29%),Na-HCO3-SO4型(22.86%)和Na-SO4-HCO3(17.15%)等,结合高氟矿井水较高的TDS值(表1),表明高氟矿井水具有更为强烈的水-岩相互作用过程[18,20]。
神东矿区矿井水中F-质量浓度为0.16~12.75 mg/L,平均值为5.01 mg/L,最高值偏高于苏双清2020年7.40 mg/L的调查结果[23],根据《生活饮用水卫生标准》,有77.78%的样品超过了1.00 mg/L的标准。高氟矿井水F-平均质量浓度(5.01±3.76 mg/L)约为低氟矿井水F-平均质量浓度(0.41±0.27 mg/L)的10.66倍。矿井水中C(F-)<1 mg/L,1 mg/L<C(F-)<4 mg/L,4 mg/L<C(F-)<10 mg/L和C(F-)>10 mg/L的样品数分别占总样品数的24.44%,20.00%,44.44% 和11.12%,表明神东矿区矿井水中F-质量浓度主要分布在4~10 mg/L,具有典型高氟水的特征,与LI等[25]研究运城盆地高氟地下水F-质量浓度分布特征相似。
从空间分布看(图1),神东矿区矿井水F-质量浓度呈现西北高,东南低的态势。高氟矿井水在空间上表现出高度的异质性,主要分布在西北部的布尔台矿和乌兰木伦矿等矿区。而从垂向上看(图3),矿井水中F-质量浓度也具有明显的垂向分带性,由浅至深呈逐渐升高趋势,高氟矿井水主要集中分布在埋深150~300 m,据调查发现,主要为3-1煤和5-2煤开采后产生的矿井水的埋深范围,可能与煤层沉积环境有关。
图3 不同埋深的矿井水氟化物分布
Fig.3 Depth distribution of F- in Shendong mine water samples
不同水体F-丰度水平也存有巨大的差异,地表水、延安组、萨拉乌苏组和直罗组地下水中F-的质量浓度分别为0.83~9.94,0.41~7.98,0.43~1.33和0.48~0.64 mg/L,依照《生活饮用水卫生标准》III类标准(1.0 mg/L),各类水中超标个数分别占该类水总数的66.67%,50.00%,25.00%和0。相比矿井水,萨拉乌苏组和直罗组中F-质量浓度水平偏低,作为高氟矿井水的F-源的可能性低,而延安组和地表水中F-质量浓度却与矿井水接近,则表明延安组和地表水可能为高氟矿井水中F-的来源。
神东矿区高氟和低氟矿井水氢氧同位素δD值分别为-86.00‰~-65.00‰和-72.00‰~-68.00‰,δ18O值分别为-10.60‰~-8.40‰和-9.20‰~-8.30‰。97.14%高氟矿井水和100%低氟矿井水样品点均分布于全球降雨线(δD=8δ18O+10.00)和当地降雨线(δD=8.36δ18O+14.50)下方,呈现出明显的蒸发效应[26]。图4中高氟和低氟矿井水的氢氧同位素被分成了2块区域,其中低氟矿井水样品与萨拉乌苏组,直罗组和16.67%延安组的氢氧同位素样品聚集在图4小圈内,表明低氟矿井水可能受其共同补给影响,可能与采煤引发的塌陷导致顶板含水层串通有关[27];而高氟矿井水样品与83.37%延安组,25%地表水样品聚集在图4大圈内,说明高氟矿井水可能受延安组和地表水双重补给影响。经勘查发现,该地区地表水为神东矿区最低排泄点,其标高(历史最高1 039 m)普遍低于现有矿井水(1 060~1 120 m),故地表水作为矿井水补给源的可能性极低。地表水氢氧同位素特征与矿井水相似,可能受矿井水排泄于地表水有关[24]。随着矿井水中F-质量浓度的增加,氢氧同位素从萨拉乌苏组为主导的端元过渡到延安组为主导的端元,可证实延安组是高氟矿井水最主要的补给源。
图4 不同水体氢氧同位素分布
Fig.4 Distribution of hydrogen and oxygen isotopes in different water samples
矿井水中F-与其它水化学离子相互关系如图5所示。图5中矿井水和与F-均呈正相关关系,其中TDS,Na+和与F-呈现出中等正相关关系,表明矿井水中高TDS,Na+和的质量浓度对F-富集有积极作用,而F-与Ca2+,K+和Mg2+呈负相关,其中F-仅与Ca2+呈现中等负相关关系,说明矿井水中高Ca2+对F-的富集也有抑制影响。结合表1可知,偏碱性,高低Ca2+是形成高氟矿井水的主要水化学环境特征,该结果与前人研究高氟地下水的水环境特征基本保持一致[1,9,14,18]。
图5 矿井水中F-与其它水化学离子相互关系
Fig.5 Correlation coefficients of the different hydrogeochemical parameters with F- in Shendong mine water samples
神东矿区高氟和低氟矿井水中F-质量浓度均与Ca2+呈负相关关系(图6),表明F-与Ca2+溶解平衡关系是控制矿井水中F-质量浓度的重要因素。88.89%的高氟矿井水和100%的低氟样品均位于萤石溶解平衡线(Ksp=10-10.6)下方,表明萤石溶解未饱和,而11.12%的高氟样品位于萤石溶解平衡线上方,可能与其它含氟矿物的溶解有关。前人研究发现神东矿区地层中广泛存有含氟矿物,主要有萤石(CaF2),氟磷灰石(Ca5(PO4)3F)和闪长石(NaCa2(Mg,Fe,Al)5(Al,Si)8O22F2)等,总含氟量平均值为803.5 mg/kg,可溶性氟含量占总氟含量的4.10%~9.45%[28-29]。这些含氟矿物的溶解过程主要如式(1)~(4)表述:
(1)
CaF2+H2O+CO2(g)→CaCO3+2F-+2H+
(2)
NaCa2(Mg,Fe,Al)5 (Al,Si)8O22F2+2OH-→
NaCa2(Mg,Fe,Al)5(Al,Si)8O22(OH)2+2F-
(3)
Ca5(PO4)3F+OH-→Ca5(PO4)3(OH)+F-
(4)
图6 矿井水中Ca2+和F-相互关系
Fig.6 Relationship between Ca2+ and F- in mine water samples
无论含氟矿物进行哪种水-岩反应,式(1)~(4)均导致了矿井水中F-质量浓度升高,和Ca2+质量浓度的降低(沉淀),因此矿井水中F-的来源受萤石,氟磷灰石和闪长石等含氟矿物的溶解控制。
为一步验证含氟矿物的溶解与沉淀作用,运用水化学PHREEQC软件模拟了矿井水中方解石,白云石,石膏,岩盐和萤石等矿物的溶解过程。图7(a)中,低氟水和高氟水中绝大多数的方解石和白云石饱和系数SI值>0,呈过饱和状态,表明方解石和白云石过饱和是促进矿井水中F-质量浓度升高的重要因素[7,9,12]。这是因为方解石和白云石的饱和降低了矿井水中Ca2+和Mg2+的质量浓度,加速了萤石,氟磷灰石和闪长石等含氟矿物的溶解,导致式(1)~(4)将向右移动,造成矿井水中F-质量浓度的升高。图7(b)在方解石过饱和区高氟矿井水F-的活度值随F-质量浓度的升高逐渐增大,进一步说明了可以通过改变方解石的溶解度来控制矿井水形成较高质量浓度的F-[30]。图7中石膏和岩盐的SI值均<0,呈未饱和状态,暗示矿井水溶解了更多的和Na+,Cl-等,而矿井水中Na+质量浓度的升高,也有助于F-质量浓度的进一步提升。
图7 F-与多种矿物饱和系数SI关系
Fig.7 Relationship between F- and saturation coefficient SI of various minerals
萤石的饱和系数SI在不同F-质量浓度的高氟矿井水中存有明显差异,当1 mg/L<C(F-)<4 mg/L时,所有的高氟矿井水样品SI值均小于0,呈未饱和状态;当4 mg/L<C(F-)<10 mg/L时,13.64%样品SI值大于0,部分样品开始出现过饱和状态;当C(F-)>10 mg/L,40.00%样品SI值大于0,越来越多的样品出现过饱和状态。同时图7(c)中,F-质量浓度与萤石的饱和系数SI呈正相关关系,矿井水中F-质量浓度越高,萤石的饱和系数SI值就会越大,并逐渐向过饱和过渡,进一步说明了含氟矿物(如萤石)的溶解对矿井水中氟浓度的贡献较大[10-11]。60%的C(F-)>10 mg/L的高氟矿井水样品未出现萤石过饱和状态,其SI值为-0.38~-0.13,暗示除了萤石溶解外,可能还存有其他含氟矿物的溶解。
Gibbs图常常用来反映地下水的蒸发效应[1,3]。由图8(a)可知,神东矿区高氟和低氟矿井水样品均位于含氟矿物溶解区域和蒸发效应区间,表明其水化学成分主要受含氟矿物溶解和蒸发效应共同控制。低氟矿井水n(Na+)/[n(Na+)+n(Ca2+)]值范围较宽(0.21~0.82),其中,n为物质的量。90%样品的TDS 值小于1 000 mg/L,表明含氟矿物溶解作用是控制低氟矿井水水化学成分的主要影响因素。随着矿井水中F-质量浓度的增加,1 mg/L<C(F-)<4 mg/L,4 mg/L<C(F-)<10 mg/L和C(F-)>10 mg/L样品点沿着含氟矿物溶解区向右平移,n(Na+)/[n(Na+)+n(Ca2+)]值变大并趋于1.0,TDS值也有升高,表明高氟矿井水具有更剧烈的水-岩作用[18,20]。所有高氟矿井水样品位于含氟矿物溶解区和蒸发效应区之间,而濒临含氟矿物溶解区,表明蒸发浓缩作用不是控制高氟矿井水产生的重要因素。
图8 矿井水中Gibbs图和F-与C(F-)/C(Cl-)对比关系
Fig.8 Plots of Gibbs and F- concentration vs C(F-)/ C(Cl-) of the Shendong mine water samples
由图8(b)可知,所有高氟矿井水的样品点均落在含氟矿物溶解和F-地质成因区(F-质量浓度和C(F-)/C(Cl-)值均较高),而远离蒸发作用区(F-质量浓度较高,而C(F-)/C(Cl-)值较低),进一步验证了蒸发浓缩作用不是控制高氟矿井水产生的重要因素[12,20]。所有1 mg/L<C(F-)<4 mg/L样品点位于图8(b)的左侧(含氟矿物溶解区),且随着C(F-)/C(Cl-)值的升高,F-质量浓度值没有明显变化,呈现主要受含氟矿物溶解控制;而4 mg/L<C(F-)<10 mg/L和C(F-)>10 mg/L样品点却落在含氟矿物溶解和F-地质成因之间,随着F-质量浓度的增加,逐渐远离海水线(C(F-)/C(Cl-)=0.02),C(F-)/C(Cl-)值也随之升高,说明4 mg/L<C(F-)<10 mg/L和C(F-)>10 mg/L的高氟矿井水主要受含氟矿物溶解和F-地质成因双重控制。F-地质成因可能与蒸发岩和岩盐的溶解有关,因为高氟水中石膏和岩盐的饱和溶解系数SI均小于0(图7(a)),同时水化学离子质量浓度均值也明显高于低氟水(表1),这两组数据均证实了高氟矿井水中确实存在蒸发岩和岩盐溶解过程,而蒸发岩和岩盐的溶解往往也会提高高氟矿井水中的TDS值[17,19,31],造成了图8(a)中所有高氟矿井水的TDS值均偏高于含氟矿物溶解区。
阳离子交换作用可以降低地下水中Ca2+/Na+比值,引起方解石饱和沉淀,促进萤石,氟磷灰石和闪长石等含氟矿物的进一步溶解,释放出更多的F-到地下水中。阳离子交换的强度常常使用Scholler方法中的CAI1和CAI2指标来表达[1,25],其计算方法分别为
(4)
(5)
当CAI1和CAI2的值均为正数时,表明地下水中的Na+,K+被固相中的Ca2+交换;相反,当CAI1和CAI2的值均为负数时,表明地下水的Ca2+被固相中的Na+,K+阳离子交换,CAI1和CAI2的值为0时,表明阳离子交换作用未发生;CAI1和CAI2的绝对值越大,阳离子交换作用强度越强[1,25]。图9(a)中,所有的高氟矿井水样品CAI 1和 CAI2值均为负值,表明高氟矿井水中的Ca2+均被固相中的Na+,K+进行了阳离子交换,导致Na+,K+质量浓度的增加。随着矿井水中F-质量浓度的升高,CAI1和CAI2绝对值逐渐越大,表明阳离子交换作用越来越强,其中C(F-)>10 mg/L样品点CAI1和CAI2绝对值最大,分别为63.69和2.50,离子交换作用也最激烈,表明离子交换作用也是控制高氟矿井水产生的重要因素。
图9 矿井水中CAI1与CAI2和TDS与n(Na+)/ [n(Na+)+n(Cl-)]相互关系
Fig.9 Plots of CAI1 vs CAI2 and TDS vs n(Na+)/[n(Na+)+ n(Cl-)] of water samples
阳离子交换作用也可以用n(Na+)/[n(Na+)+n(Cl-)]与TDS的关系来阐述,如果地下水中Na+全部来源于海水,n(Na+)/[n(Na+)+n(Cl-)]应等于或接近0.50。若地下水中n(Na+)/[n(Na+)+n(Cl-)]大于0.50,则发生了阳离子交换作用,此时地下水的Ca2+会被固相中的Na+阳离子交换;n(Na+)/[n((Na+)+n(Cl-)]越接近1.00,阳离子交换作用就越强[32-33]。图9(b)中高氟矿井水样品n(Na+)/[n(Na+)+n(Cl-)]均大于0.50,表明高氟矿井水的Ca2+均被固相中的Na+阳离子交换。C(F-)>10 mg/L样品点n(Na+)/[n(Na+)+n(Cl-)]最高(0.86~0.98),表明离子交换作用也最激烈,进一步验证了离子交换作用也是控制高氟水产生的重要因素。
F-为阴离子,离子半径为0.133 nm,与OH-离子半径相当,在地下水中极易与同为阴离子的OH-和发生竞争吸附作用[34]。图10(a)中高氟矿井水的F-质量浓度与质量浓度呈正相关关系(R=0.64),表明高氟矿井水中质量浓度增高有利于F-质量浓度的提升,发生了明显的竞争吸附作用,主要是因为高氟矿井水与空气充分接触,溶解了大量的CO2,矿井水中质量浓度随之升高,而高质量浓度的能够置换出黏土矿物(煤)表面吸附态的F-,导致F-被释放到矿井水中,有利于高氟矿井水中F-的富集[10-12]。
图10(b)中高氟矿井水的F-质量浓度与pH关系较弱(R=0.16),表明高氟矿井水的F-质量浓度与pH并未发生明显的竞争吸附作用。研究区高氟矿井水的pH值在6.25~8.50(表1),虽然高氟矿井水呈弱碱性或碱性,但OH-并未抑制F-在黏土矿物表面的吸附。
和Cl-常常用于分析地下水人为污染的程度[12,35]。一般认为地下水质量浓度超过5 mg/L时,存在人为生活污染[12,36]。由图11(a)可知,95.56%的高氟矿井水样品质量浓度未超过5 mg/L,且与Cl-相关性很弱,表明绝大多数高氟矿井水并未受到人为污染,而高的Cl-质量浓度可能与岩盐溶解有关,图7(a)中岩盐饱和系数SI值小于0可验证。图11(b)中高氟矿井水中质量浓度并未随F-质量浓度的增加而升高,表明矿井水虽极易受人类活动影响,但人为活动并不是形成高氟矿井水的主要因素。
图10 矿井水中和pH分别与F-相互关系
Fig.10 Plots of vs F- and pH vs F- of the Shendong mine water samples
图11 矿井水中F-和Cl-分别与相互关系
Fig.11 Plots of vs Cl-and vs F- of the Shendong mine water samples
(1)神东矿区矿井水中F-质量浓度在0.16~12.75 mg/L,平均值为5.01 mg/L,有77.78%的样品超过了《生活饮用水卫生标准》(GB5749—2006)。高氟矿井水F-质量浓度主要分布在4~10 mg/L,占总样品数的44.44%。从空间分布看,神东矿区矿井水中F-质量浓度呈现西北高,东南低的态势,高氟矿井水主要分布在布尔台矿和乌兰木伦矿等矿区。垂向上看,高氟地下水主要集中分布在埋深150~300 m,主要为3-1煤和5-2煤开采产生的矿井水的埋深范围。
(2)偏碱性,高和低Ca2+是形成高氟矿井水的水化学环境。低氟矿井水与萨拉乌苏组,直罗组和16.67%延安组的氢氧同位素特征相似,表明低氟矿井水可能受其共同补给影响;而高氟矿井水与83.37%延安组和25%地表水样品的氢氧同位素特征相似,说明高氟矿井水受延安组和地表水双重补给影响,而地表水为该地区的最低排泄点,故延安组是高氟矿井水最主要的补给源。
(3)高氟矿井水的形成主要受含氟矿物的溶解控制,方解石和白云岩的溶解饱和影响了萤石等含氟矿物的溶解平衡,造成矿井水中F-质量浓度的升高;阳离子交换作用降低了矿井水中Ca2+的质量浓度,的竞争吸附置换出了黏土矿物表面吸附态的F-,促进了F-在矿井水中富集,而蒸发浓缩和人类活动对高氟矿井水形成贡献较小。
[1] HE Jin,AN Yonghui,ZHANG Fucun.Geochemical characteristics and fluoride distribution in the groundwater of the Zhangye Basin in Northwestern China[J].Journal of Geochemical Exploration,2013,135:22-30.
[2] AGHAPOUR S,BINA B,TARRAHI M J,et al.Distribution and health risk assessment of natural fluoride of drinking groundwater resources of Isfahan,Iran,using GIS[J].Environmental Monitoring and Assessment,2018,190(3):137.
[3] LAXMANKUMAR D,SATYANARAYANA E,DHAKATE R.Hydrogeochemical characteristics with respect to fluoride contamination in groundwater of Maheshwarm mandal,RR district,Telangana state,India[J].Groundwater for Sustainable Development,2019,8:474-483.
[4] PANKAJ K,SINGH C K,SARASW C,et al.Evaluation of aqueous geochemistry of fluoride enriched groundwater:A case study of the Patan district,Gujarat,Western India[J].Water Science,2019,31(2):215-229.
[5] KNAPPETT P S K,LI Yanmei,HERNANDEZ H,et al.Changing recharge pathways within an intensively pumped aquifer with high fluoride concentrations in Central Mexico[J].Science of the Total Environment,2018,622-623:1029-1045.
[6] MEENAKSHI,GARG V K,KAVITA,et al.Groundwater quality in some villages of Haryana,India:Focus on fluoride and fluorosis[J].Journal of Hazardous Materials,2004,106(1):55-60.
[7] MONDAL D,GUPTA S,REDDY D V,et al.Geochemical controls on fluoride concentrations in groundwater from alluvial aquifers of the Birbhum district,West Bengal,India[J].Journal of Geochemical Exploration,2014,145:190-206.
[8] KALPANA L,BRINDHA K,FIMAR L E.A new fluoride index to mitigate geogenic contamination by managed aquifer recharge[J].Chemosphere,2019,220:381-390.
[9] HE Xin,MA Tenghua,WANG Yanxin,et al.Hydrogeochemistry of high fluoride groundwater in shallow aquifers,Hangjinhouqi,Hetao Plain[J].Journal of Geochemical Exploration,2013,135:63-70.
[10] 毛若愚,郭华明,贾永锋,等.内蒙古河套盆地含氟地下水分布特点及成因[J].地学前缘,2016,23(2):260-268.
MAO Ruoyu,GUO Huaming,JIA Yongfeng,et al.Distribution characteristics and genesis of flouoride groundwater in the Hetao basin,Inner Mongoli[J].Earth Science Frontiers (in Chinese),2016,23(2):260-268.
[11] GUO Huaming,ZHANG Yang,XING Lina,et al.Spatial variation in arsenic and fluoride concentrations of shallow groundwater from the town of Shahai in the Hetao Basin,Inner Mongolia[J].Applied Geochemistry,2012,27(11):2187-2196.
[12] LI Danni,GAO Xubo,WANG Yanxin,et al.Diverse mechanisms drive fluoride enrichment in groundwater in two neighboring sites in northern China[J].Environmental Pollution,2018,237:430-441.
[13] PI Kunfu,WANG Yanxin,XIE Xianjun,et al.Hydrogeochemistry of co-occurring geogenic arsenic,fluoride and iodine in groundwater at Datong Basin,Northern China[J].Journal of Hazardous Materials,2015,300:652-661.
[14] XIAO Jun,JIN Zhangdong,ZHANG Fei.Geochemical controls on fluoride concentrations in natural waters from the middle Loess Plateau,China[J].Journal of Geochemical Exploration,2015,159:252-261.
[15] 邢丽娜,郭华明,魏亮,等.华北平原浅层含氟地下水演化特点及成因[J].地球科学与环境学报,2012,34(4):57-67.
XING Lina,GUO Huaming,WEI Liang,et al.Evolution feature and gensis of fluoride groundwater in shallow aquifer from North China Plain[J].Journal of Earth Sciences and Environment (in Chinese),2012,34(4):57-67.
[16] GAO Hongjian,JIN Youqian,WEI Junling.Health risk assessment of fluoride in drinking water from Anhui Province in China[J].Environmental Monitoring and Assessment,2013,185(5):3687-3695.
[17] KUMAR S,VENKATESH A S,SINGH R,et al.Geochemical signatures and isotopic systematics constraining dynamics of fluoride contamination in groundwater across Jamui district,Indo-Gangetic alluvial plains,India[J].Chemosphere,2018,205:493-505.
[18] OLAKA L A,WILKE F D H,OLAGO D O,et al.Groundwater fluoride enrichment in an active rift setting:Central Kenya Rift case study[J].Science of the Total Environment,2016,545-546:641-653.
[19] RASHID A,GUAN DONGXING,FAROOQI A,et al.Fluoride prevalence in groundwater around a fluorite mining area in the flood plain of the River Swat,Pakistan[J].Science of the Total Environment,2018,635:203-215.
[20] DEHBANDI R,MOORE F,KESHAVARZI B.Geochemical sources,hydrogeochemical behavior,and health risk assessment of fluoride in an endemic fluorosis area,central Iran[J].Chemosphere,2018,193:763-776.
[21] 赵焰,陆梦楠,孙斌,等.典型煤矿含氟矿井水工业化除氟技术研究[J/OL].煤炭科学技术:1-8[2021-06-17].http://kns.cnki.net/kcms/detail/11.2402.TD.20200619.1508.006.html.
ZHAO Yan,LU Mengnan,SUN Bing,et al.Research on practical fluoride removal technology for typical fluorides-contaminated mine drainage[J/OL].Coal Science and Technology:1-8[2021-06-17]. http://kns.cnki.net/kcms/detail/11.2402.TD.20200619.1508.006.html.
[22] 孙亚军,陈歌,徐智敏,等.我国煤矿区水环境现状及矿井水处理利用研究进展[J].煤炭学报,2020,45(1):304-316.
SUN Yajun,CHEN Ge,XU Zhimin,et al.Research progress of water environment,treatment and utilization in coal mining areas of China[J].Journal of China Coal Society,2020,45(1):304-316.
[23] 苏双青,赵焰,徐志清,等.我国煤矿矿井水氟污染现状及除氟技术研究[J].能源与环保,2020,42(11):9-14.
SU Shuangqing,ZHAO Yan,XU Zhiqing,et al.Status quo of fluoride pollution of coal mine water in China and research on fluoride removal technology[J].China Energy and Environmental Protection,2020,42(11):9-14.
[24] 付昌昌,李向全,马剑飞,等.窟野河流域中游煤矿区地下水质量及补给来源研究[J].水文,2018,38(6):42-47.
FU Changchang,LI Xiangquan,MA Jianfei,et al.Groundwater quality assessment and recharge sources of groundwater in coal mining area of Kuye River Basin[J].Journal of China hydrology (in Chinese),2018,38(6):42-47.
[25] LI Chengcheng,GAO Xubo,WANG Yanxin.Hydrogeochemistry of high-fluoride groundwater at Yuncheng Basin,northern China[J].Science of the Total Environment,2015,508:155-165.
[26] AYADI Y,MOKADEM N,BESSER H,et al.Hydrochemistry and stable isotopes (δ18O and δ2H) tools applied to the study of karst aquifers in southern mediterranean basin (Teboursouk area,NW Tunisia)[J].Journal of African Earth Sciences,2018,137:208-217.
[27] 徐祝贺,李全生,李晓斌,等.浅埋高强度开采覆岩结构演化及地表损伤研究[J].煤炭学报,2020,45(8):2728-2739.
XU Zhuhe,LI Quansheng,LI Xiaobin,et al.Structural evolution of overburden and surface damage caused by high-intensity mining with shallow depth[J].Journal of China Coal Society,2020,45(8):2728-2739.
[28] 段磊,王晓轩,孙亚乔,等.煤矸石中氟的赋存特性及生态风险评价[J].煤炭转化,2021,44(2):1-10.
DUAN Lei,WANG Xiaoxuan,SUN Yaqiao,et al.Occurrence characteristics and ecological risk assessment of fluorine in coal gangue[J].Coal Conversion,2021,44(2):1-10.
[29] 赵波,魏宁.陕西煤中氟的含量分布及赋存形态研究[J].煤质技术,2020,35(6):34-39.
ZHAO Bo,WEI Ning.Study on content distribution and occurrence form of fluorine in Shaanxi coal[J].Coal Quality Technology (in Chinese),2020,35(6):34-39.
[30] BIDDAU R,CIDU R,LORRAI M,et al.Assessing background values of chloride,sulfate and fluoride in groundwater:A geochemical-statistical approach at a regional scale[J].Journal of Geochemical Exploration,2017,181:243-255.
[31] RAFIQUE T,NASEEM S,OZSYATH D,et al.Geochemical controls of high fluoride groundwater in Umarkot Sub-District,Thar Desert,Pakistan[J].Science of the Total Environment,2015,530-531:271-278.
[32] ALI S,SHEKHAR S,BHATTACHARYA P,et al.Elevated fluoride in groundwater of Siwani Block,Western Haryana,India:A potential concern for sustainable water supplies for drinking and irrigation[J].Groundwater for Sustainable Development,2018,7:410-420.
[33] HAO Chunming,ZHANG Wei,GUI Herong.Hydrogeochemistry characteristic contrasts between low-and high-antimony in shallow drinkable groundwater at the largest antimony mine in Hunan Province,China[J].Applied Geochemistry,2020,117:104584.
[34] 涂成龙,何令令,崔丽峰,等.氟的环境地球化学行为及其对生态环境的影响[J].应用生态技术,2019,30(1):21-29.
TU Chenglong,HE Lingling,CUI Lifeng,et al.Environmental and geochemical behaviors of fluorine and its impacts on ecological environment[J].Chinese Journal of Applied Ecology,2019,30(1):21-29.
[35] ALI W,ASLAM M W,JUNAID M,et al.Elucidating various geochemical mechanisms drive fluoride contamination in unconfined aquifers along the major rivers in Sindh and Punjab,Pakistan[J].Environmental Pollution,2019,249:535-549.
[36] HEATON T H E.Isotopic studies of nitrogen pollution in the hydrosphere and atmosphere:A review[J].Chemical Geology:Isotope Geoscience Section,1986,59:87-102.
HAO Chunming,ZHANG Wei,HE Ruimin,et al.Formation mechanisms for elevated fluoride in the mine water in Shendong coal-mining district[J].Journal of China Coal Society,2021,46(6):1966-1977.