基于环境稳定同位素法分析煤矿开采对岩溶水的影响量

来源:用户上传      作者:张建友 张鑫 高旭波 李成城

　　摘要：针对煤矿开采对岩溶水资源造成的影响，以太原市西山某煤矿为例，综合采用硫、锶和氢氧同位素特征分析法，研究了岩溶地下水中同位素的主要来源和同位素演化的主要水岩相互作用，分析计算了矿坑水和矿区裂隙水中岩溶水的混合比例，定量评估了煤矿开采对岩溶地下水的影响量。结果表明：西山某煤矿矿坑水中岩溶水（含地表水）的比例为51％，裂隙水中岩溶水（含地表水）的比例为60％，据此计算得出煤矿开采影响岩溶地下水量为129.5万 m3／a。
　　关键词：环境同位素法；岩溶水影响量；煤矿开采；山西省
　　中图分类号：P641.134；TV211.1＋2 文献标志码：A doi：10.3969／j.issn.1000－1379.2022.10.019
　　引用格式：张建友，张鑫，高旭波，等.基于环境稳定同位素法分析煤矿开采对岩溶水的影响量[J].人民黄河，2022，44（10）：101－104.
　　Analysis of the Influence of Coal Mining to the Karst Water Volume Based on Environmental Stabilization Isotope Method
　　ZHANG Jianyou 1，ZHANG Xin2，GAO Xubo2，LI Chengcheng2
　　（1.Shanxi Water Resources Research Institute Co.，Ltd.，Taiyuan 030001，China； 2.School of Environmental Studies ，China University of Geosciences ，Wuhan 430074，China）
　　Abstract：In view of the serious impact of coal mining on karst water resources ，this paper took a coal mine of Xishan in Taiyuan as an ex?ample and adopted the characteristic analysis method of sulfur strontium and hydrogen and oxygen isotopes to study the main sources of iso?topes in karst groundwater calculate the mine drainage and the mixing ratio of karst water in the fissure water in the mining area and quanti?tatively evaluate the impact on karst groundwater.The results show that in the coal mine in Xishan，the proportion of karst water （ including surface water） in the mine drainage is 51％，and the proportion of karst water （ including surface water） in the fissure water is 60％，and theamount of karst water affected by coal mining is calculated as 1.295 million m3／a.
　　Key words：environmental isotopic method；karst water volume；coal mining；Shanxi Province
　　山西作为国家新型综合能源与制造业基地，1949年以来全省范围内对煤炭资源进行了大规模、高强度开采，造成岩溶大泉衰减，甚至干枯断流。我国目前存在煤矿开采破坏岩溶泉域系统方面的研究空白，无法定量评估煤矿开采对岩溶水的影响量，管理部门在缺乏依据的情况下，无法要求煤矿企业开展岩溶泉域系统生态修复工作。环境同位素法作为研究地下水系统水文地球化学过程的有效手段，不仅可以探究地下水的起源、形成、埋藏条件、水质变化等，而且可以判定地下水的补给来源、补给强度、地下水与地表水之间的联系[1－4]。笔者以太原西山某煤矿为例，采用硫、锶和氢氧同位素分析法，揭示影响地下水中这些同位素演化的主要水－岩相互作用和矿坑水对岩溶地下水的影响，定量评估煤矿开采对岩溶水的影响量。
　　1 研究区概况
　　研究区位于山西省太原市古交市、吕梁山东翼，属中低山区，区内地形切割强烈。全区地势西南高、东北低，最高处海拔为1400 m ，最低处海拔为970 m 。研究区属暖温带半干旱大陆性季风气候区，多年平均降水量为419.4 mm ，降水年内分布极不均匀，主要集中在7―9月。多年平均蒸发能力1871.8 mm 。多年平均气温8.1℃，极端最高气温39.4℃，极端最低气温－25.5℃。
　　某煤矿位于古交市，1988年原国家计划委员会批准原煤生产能力为400万 t／a，煤矿于1997年年底建成并开始试生产，2002年10月正式投产。2017年经山西省煤炭工业厅重新核定，原煤生产能力为450万 t／a，井田面积为73.3 km2，开采标高为540～950 m。1998―2017年累计生产原煤6070.6万 t，累计矿井排水量为3154.4万 m3，煤富水系数为0.52 m3／t 。2017年原煤产量为306.7万 t，矿井排水量为214.1万 m3，其中矿坑排水量为124.0万 m3、裂隙水排水量为59.4万 m3，井下水位观测井漏失岩溶地下水量为30.7万 m3，煤富水系数为0.698 m3／t。

nlc202212201110

　　煤矿井田位于晋祠泉域古交市岩溶水径流区，分布有背斜、向斜、断层等地质构造，位于西部的鲜则沟背斜和元家山向斜都为压性结构面，导水性弱，水力联系差；位于北部的风坪岭断层岩溶透水性较差，地下水流动缓慢；位于东部的古交断层两侧水位相差20 m，岩溶水由北东向西南缓慢流动。按照成因、岩性和地下水赋存特征可把井田内地下水划分为松散岩类孔隙含水层组、碎屑岩类裂隙含水层组和碳酸盐岩类裂隙岩溶含水层组。松散岩类孔隙水分布于原平河河谷及其一级阶地，岩性主要为冲、洪积作用堆积的砂砾石、中粗砂等，厚度20～30 m；碎屑岩类裂隙水主要分布于西山山区及太原断陷盆地，二叠系的砂、泥岩互层至上石炭统山西组地层含风化裂隙潜水和层间裂隙承压水；碳酸盐岩类裂隙岩溶水在整个研究区均有分布，含水层主要为中奥陶统石灰岩，其次为泥灰岩夹石膏层。其中碳酸盐岩类裂隙岩溶地下水流场总体呈现向北突起的弧形，岩溶水位为925～950 m 。井田内8号煤层为现状主采煤层，底板标高为540～860 m，低于岩溶水位90～385 m，均为岩溶水带压开采区。
　　2 样品采集和测试
　　在对研究区水文地质l件详细调查的基础上，按照补给区、径流区和排泄区顺序采集地下水（含岩溶水、裂隙水和矿坑水）和地表水。共采集10组地下水样品和1组地表水样品（ S1），其中：在补给―径流区（矿区上游）采集4组岩溶水样品（ Y1、Y2、Y3、Y4），在矿区采集1组矿坑水样品（ K1）和1组裂隙水样品（ L1），在径流―排泄区（矿区下游）采集4组裂隙水样品（ L2、L3、L4、L5，即西山裂隙水），采样点分布见图1。地下水样品采集时，所有水样均在现场用0.45μm 的滤膜过滤后装入聚乙烯瓶。滤膜提前在实验室用3％的稀硝酸浸泡24 h，再用去离子水多次冲洗，备用。水位、水量、水温、pH 值、碱度、电导率、浑浊度、色、臭和味、肉眼可见物等指标现场完成监测。野外样品采集后在4℃条件下冷藏，在一周内完成主要离子和同位素检测工作。主要阴离子采用离子色谱仪测定；水样氢氧同位素（δ D 和δ18 O）采用水同位素分析仪测定，分析误差为δ D＜0.5‰，δ18 O＜0.1‰；水样硫同位素（δ34 S）采用元素分析仪－稳定同位素比值质谱仪（ EA－ IRMS）测定；锶同位素（87 Sr／86 Sr）采用多接收电感耦合等离子体质谱仪（ MC－ICP－MS）测定。
　　3 矿坑水稳定同位素特征分析
　　3.1 氢氧同位素特征
　　蒸发过程中，轻的水分子 H 2（1）6 O 比包含有重同位素的分子（即H2 H 16 O 或 H 2（1）8 O）更为活跃。因此，从海面蒸发的水蒸气比海水中的δ18 O 少12‰～15‰、比δ D 少80‰～120‰。受同位素分馏作用影响，大气降水中δ18 O 、δ D 含量由降水形成时的温度控制。通过对地下水中δ18 O 和δ D 的关系与当地大气降水线进行对比，可以分析得出地下水的补给来源[5－6]。
　　为了研究矿区裂隙水和矿坑水的稳定同位素特征，对矿区水（矿坑水和裂隙水）和矿区上游岩溶水氢氧同位素进行分析。全球大气降水线方程为δ D＝7.9δ18 O＋8.2[7]，太原地区大气降水线方程为δ D ＝6.42δ18 O－4.66[8]。由图2可以看出，所采水样基本落在太原地区大气降水线附近或下方，说明研究区含水层有统一的补给来源，都来自大气降水补给，而δ D 和δ18 O 值不同说明受蒸发作用影响不同。同时可以看出，矿区裂隙水δ D 和δ18 O 值与部分岩溶水的值接近，说明矿区裂隙水和岩溶水均受到较强的水岩作用或外源地表水的入渗补给。矿坑水落在了太原地区大气降水线附近，表明其补给来源主要为大气降水。
　　3.2 硫同位素特征
　　在矿区地下水中，SO4（2）－是地下水中重要的阴离子。自然界δ34 S 值具有变化范围广（－65‰～120‰）、分馏作用明显的特点：现代海洋水δ34 S 值为21.0‰；沉积岩中硫化物、石油和煤中的δ34 S 多为负值。硫同位素分馏过程记录了各种环境信息，可以揭示水循环系统中硫的来源、迁移和转化，有助于认识自然界中各种含硫物质的演变历史，在研究硫酸盐的来源时具有独特的优势。我国北方煤总硫的δ34 S 值平均为3.68‰[9] 。山西中奥陶统石膏层δ34 S 值为23.8‰～31.4‰，上覆地层中石炭－二叠系煤系地层黄铁矿中δ34 S 值为－9.1‰－4.8‰，矿坑水中δ34 S 值为－13.6‰～8.0‰，新近系黄土中δ34 S 值为5.5‰～9.5‰，汾河及其支流与引黄水中δ34 S 值为7.9‰～10.7‰[10]。这些数据说明在山西中奥陶统与上覆含硫较高的地层中，δ34 S 同位素存在较为明显的分馏作用，对北方岩溶水中 SO4（2）－的两个主要物源黄铁矿氧化、石膏溶解具有良好的指示作用，可以利用δ34 S 值来判断和解析岩溶水中 SO4（2）－的来源[2，11] 。
　　为了研究矿区地下水与矿坑水δ34 S 同位素的特征，对部分水样的δ34 S 同位素进行分析，δ34 S 与 SO4（2）－的关系见图3。矿区矿坑水和裂隙水具有较低的δ34 S 值和较大的 SO4（2）－质量浓度，δ34 S 值分别为－1.0‰和－2.9‰，SO4（2）－质量浓度均高于450 mg／L，说明矿区矿坑水 SO4（2）－主要来源于黄铁矿氧化。矿区裂隙水除受上游岩溶水补给外，矿坑水的混入是其δ34 S 值偏负的主要原因。
　　3.3 锶同位素特征分析
　　Sr 同位素不会因物理、化学风化和生物过程而发生分馏。也就是说，水体中溶解的 Sr 同位素在地下水演化过程中，不会因衰变而发生变化，也不会因 Sr 从水体中分离（矿物沉淀或离子交换）而发生87 Sr 或86 Sr 的变化。不同来源水体的化学及同位素特征不同，其演化主要受下列条件影响：水体在岩层中滞留和迁移时间（速度）、初始水成分、水流路径的曲折程度、水体间的混合作用、所流经区域矿物成分特征及其与水体反应情况。因此 Sr 同位素常被用于示踪地下水的水岩相互作用程度等。硅酸盐风化来源的87 Sr／86 Sr值为0.714～0.715，碳酸盐风化来源的87 Sr／86 Sr 值为0.7102～0.7108。同时研究发现，奥陶统纯灰岩的87 Sr／86 Sr 值为0.70840～0.70994，黄铁矿的87 Sr／86 Sr值为0.71400～0.72120[1，12]。

nlc202212201110

　　为了研究矿区地下水与矿坑水的 Sr 同位素特征，对部分水样的 Sr 同位素进行分析，得到87 Sr／86 Sr 值与 Sr 的关系见图4。从图4可以看出，矿区附近岩溶水 Sr 质量浓度、87 Sr／86 Sr 值分别为0.2～1.1 mg／L、0.71092～0.71179。西山裂隙水87 Sr／86 Sr 比值为0.71315～0.71532，明显高于岩溶水的。矿区矿坑水87 Sr／86 Sr 值与裂隙水的接近且高于矿区裂隙水的，同时接近黄铁矿中的值，可能受裂隙水补给和黄铁矿溶解影响。矿区裂隙水的87 Sr／86 Sr值低于西山裂隙水的，接近岩溶水的，同时 Sr 质量浓度较大，初步判断其受到岩溶水补给的同时，还受地表水等其他来源水补给。
　　4 应用锶同位素分析岩溶水影响量
　　在分析计算岩溶水影响量时，利用质量守恒定律和同位素守恒定律计算岩溶水在矿区裂隙水和矿坑水中的混合比例。假定地下水中某一同位素组分的浓度、质量分别为 C、N，根据质量守恒定律和同位素守恒定律可得
　　式中：Ca 、Na 和Cb、Nb 分别为 a 端元和 b 端元水中同位素组分的浓度、质量，其中Cb＞Ca 。
　　根据式（1）、式（2），在图5中做投影变换，将矿区矿坑水样点垂直投影到黄铁矿与岩溶水（含地表水）混合线（图中直线）上，按照式（3）计算岩溶水（含地表水）对矿坑水的影响量。
　　式中：P 为岩溶水的混合流量；My 为矿区上游岩溶水锶同位素组成；Mh为黄铁矿的锶同位素组成；M 为可能受到岩溶水影响的锶同位素组成。
　　矿区矿坑水和裂隙水中岩溶水的混合比例计算结果为，矿坑水中岩溶水（含地表水）混合比例约为51％，裂隙水中岩溶水（含地表水）的混合比例约为60％。矿坑排水量为124万 m3／a，裂隙水排水量为59.4万 m3／a，据此可以计算出矿坑水中岩溶地下水量为63.2万 m3／a，裂隙水中岩溶地下水量为35.6万 m3／a，加上井下观测井漏失的岩溶地下水量30.7万 m3／a，岩溶地下水影响量为129.5万 m3／a。
　　5 结语
　　通过稳定同位素法分析计算，太原市西山某煤矿矿坑水中岩溶水（含地表水）的比例为51％，裂隙水中岩溶水（含地表水）的比例为60％，计算煤矿开采影响岩溶地下水量为129.5万 m3／a 。值得注意的是，利用 Sr 同位素分析矿坑水与岩溶水的混合比例，可能受混合比例中岩溶水包含地表水等其他因素影响。为此，建议作为半定量的评估手段，分析研究煤矿开采对岩溶地下水的影响量，为生态修复与岩溶泉域保护提供科技支撑。
　　参考文献：
　　[1] 高旭波，向绚丽，侯保俊，等.水化学―稳定同位素技术在岩溶水文地质研究中的应用[J].中国岩溶，2020，39（5）：629－636.
　　[2] POROWSKI Adam，POROWSKA Dorota，HALAS Stanislaw.Identification of Sulfate Sources and Biogeochemical Processes in an Aquifer Affected by Peatland：Insights from Monitoring the Isotopic Composition of Groundwater Sulfate in Kampinos National Park ，Poland [J].Water，2019，11（7）：1388.
　　[3] GAT J R.Comments on the Stable Isotope Method in RegionalGroundwater Investigations [J].Water Resources Research，1971，7（4）：980－993.
　　[4] 臧红飞，连泽俭.柳林泉域地表水与岩溶地下水的混合作用研究[J].人民黄河，2020，42（1）：53－56，62.
　　[5] 束龙仓，宫荣，栾佳文，等.地下水与地表水水量交换识别及交换量计算：以新汴河宿州段为例[J].水科学进展，2022，33（1）：57－67.
　　[6] 侯国华，高茂生，叶思源，等.黄河三角洲浅层地下水盐分来源及咸化过程研究[J].地学前缘，2022，29（3）：145－154.
　　[7] CRAIG H.Isotopic Variations in Meteoric Waters[J].Science，1961，133（3465）：1702－1703.
　　[8] 贾振兴，臧红飞，郑秀清，等.太原地区大气降雨的氢氧同位素特征研究[J].水资源与水工程学报，2015，26（2）：22－25.
　　[9] 洪业汤，张鸿斌，朱煊，等.中国煤的硫同位素组成特征及燃煤过程硫同位素分馏[J].中国科学（ B 辑化学生命科学地学），1992，22（8）：868－873.
　　[10] 张江华，梁永平，王维泰，等.硫同位素技术在北方岩溶水资源调查中的应用实例[J].中国岩溶，2009，28（3）：235－241.
　　[11] LI C C，GAO X B ，WANG W Z ，et al.Hydro?BiogeochemicalProcesses of Surface Water Leakage into Groundwater in Large Scale Karst Water System：A Case Study at Jinci， Northern China [J].Journal of Hydrology ，2020，596：125691.
　　[12] WU X，GAO X，TAN T，et al.Sources and Pollution PathIdentification of PAHs in Karst Aquifers：An Example fromLiulin Karst Water System，Northern China[J].Journal ofContaminant Hydrology 2021241：103810.
　　【责任编辑 吕艳梅】

nlc202212201110

转载注明来源:https://www.xzbu.com/1/view-15443525.htm