煤矿采空区积水水源的识别研究

任海峰 吴青海 李西民 严由吉

摘 要：為了能够准确判别采空区积水水源，以陕西黄陵二号煤矿为例，在各含水层水质分析的基础上，运用模糊综合评判法和水化学特征分析法准确判别出采空区积水的水源。研究认为，在每个含水层采集了3个水样，共计18个，并对各含水层水样的K+、Na+、Ca2+、Mg2+、Cl-、SO2-4、HCO-3、CO2-3、Fe3+、Fe2+各种阴阳离子进行了水质分析化验，在此基础上，选取各水样阴阳例子浓度指标作为分类因素进行模糊综合评判，依据最大隶属度原则来确定积水水源，其结果为采空区积水来自侏罗系延安组含水层。对各含水层水样进行水化学特征分析，并绘制水化学背景Piper三线图和水源判别Piper三线图，发现采空区积水水源以侏罗系直罗组和延安组含水层为主。再者，顶板导水裂隙带是采空区充水的主要途径，导水裂隙带绝大部分在延安组发育，少数区段发育至直罗组底部。立足模糊综合评判法和水化学特征分析法，结合导水裂隙带发育高度，得出采空区积水主要水源是侏罗系延安组含水层，次要水源是侏罗系直罗组含水层底部。关键词：采空区积水;水化学分析;模糊综合评判;含水层;水源 中图分类号：TD 741

文献标志码：A

文章编号：1672-9315（2021）02-0246-07

DOI：10.13800/j.cnki.xakjdxxb.2021.0208开放科学（资源服务）标识码（OSID）：

Idenfification of water source in coal mine goaf

REN Haifeng1，WU Qinghai1，LI Ximin2，YAN Youji1

（1.College of Safety Science and Engineering，Xian University of Science and Technology，Xian 710054，China;

2.Shaanxi Huangling No.2 Coal Mine Co.，Ltd.，Huangling 721307，China）

Abstract：In order to accurately identify the water source of goaf ponding，taking Huangling No.2 Coal Mine as an example，the fuzzy comprehensive evaluation method and hydrochemical characteristic analysis method are used to discriminate the water source of goaf ponding，based on the water quality analysis of each aquifer.According to the research，three water samples were collected in each aquifer，a total of 18 samples，and the water quality of K+，Na+，Ca2+，Mg2+，Cl-，SO2-4，HCO-3，CO2-3，Fe3+，Fe2+ ions in each aquifer was analyzed.On this basis，the concentration indexes of Yin and Yang samples of each water sample were selected as classification factors for fuzzy comprehensive evaluation，and the product was determined according to the principle of maximum membership degree The result shows that the water in goaf comes from the aquifer of Yanan formation of Jurassic system.Then，the hydrochemical characteristics of water samples of each aquifer are analyzed，and the piper threeline map of hydrochemical background and the piper threeline map of water source discrimination are drawn.It is found that the main water sources in the goaf are aquifers of Jurassic Zhiluo Formation and Yanan formation.Moreover，the roof water diversion fracture zone is the main way to fill water in the goaf.Most of the water diversion fracture zones are developed in Yanan Formation，and a few sections are developed to the bottom of Zhiluo formation.Finally，the main source of water in goaf is the aquifer of Yanan formation of Jurassic system，and the secondary source is the bottom of the aquifer of Zhiluo Formation of Jurassic system.Key words：goaf water;hydrochemical analysis;fuzzy comprehensive evaluation;aquifer;water source

0 引 言随着煤田开采逐渐由浅部转向深部，采掘活动的影响范围逐渐扩大，采空区积水对矿井安全生产的影响日益凸显。采空区积水是威胁煤矿安全生产的严重水害因素，如何准确判别采空区积水的水源是预防矿井水害的关键性因素之一。水化学特征分析法、模糊综合评判法是判别采空区积水的主要方法[1-3]。水化学特征分析法，可以判断地下水的来源，揭示不同含水层间水力联系，有助于认识区域地下水水化学特征[4-6]。

QUANG等应用水化学分析，揭示了德国鲁尔地区煤矿地下水水文化学演变过程[7]。杨建通过对葫芦素煤矿各含水层的水化学特征分析，判别了该矿矿井水来源[8]。郭瑞等在大量水化学资料分析的基础上，将水化学分析与Fisher分析法结合判别水源[9]。VERMEULEN等在南非的一个煤矿采用水化学分析法分析采动是否破坏浅部含水层和较深含水层之间的水文地质屏障，从而确定了2个含水层之间潜在的导水区域[10]。ZHANG等采用水化学分析、Fisher判别分析等方法，对安徽古北煤矿突水水源进行了识别[11]。

DELIA等应用水化学分析法分析了地表水、地下水和矿井水之间的关系[12]，魏大勇等运用分形理论与模糊综合评判方法对恒源矿底板突水的危险区域进行分区与预测[13]。GUAN等应用模糊综合评价法和聚类分析法，对内蒙古民东一矿的突水水源进行了识别[14]。姬亚东基于聚类分析的模糊综合评判法，对顶板涌水区进行了划分[15]。李凤莲等为了提高水源判别准确率，提出一种改进的模糊综合评判模型[16]。张红梅等针对卧龙湖煤矿砂岩裂隙含水层水文地质条件，采用模糊聚类综合评判方法对该含水层富水性进行了综合评判[17]。代革联等采用水质分析对比和系统聚类分析相结合的方法，判别了象山煤矿280排矸石门突水的水源[18]。水化学特征分析法、模糊综合评判法具有各自的优点，但是，水化学特征分析法受水文地质条件的干扰较大，模糊综合评判法中样本选取、权重系数确定对判定结果影响较大。单一分析法难以顾全各种影响因素，为了更加准确判别采空区积水水源，本研究立足于将模糊综合评判方法和水化学特征分析方法有机结合起来，以陕西黄陵二号煤矿为例，综合判别采空区积水的水源。

1 水文地质特征及水样特征

1.1 水文地质特征黄陵二号煤矿地处陕北黄土高原南部，属于黄陇侏罗系煤田，位于鄂尔多斯盆地南缘。侏罗系延安组为本区含煤地层，主采煤层为2#煤层，平均厚度2.32 m，层位稳定且全区分布。开采方式为单水平、分盘区开采。地层总体倾向北西，地层倾角小于5°。按岩性组合及含水层水力性质埋藏条件等，该区含水层自上而下主要有：地表水、第四系、白垩系、侏罗系。第四系松散层含水层岩性为细砂，属弱-极弱富水的含水层。白垩系洛河组含水层岩性以砂岩为主，富水性属中等-强。侏罗系自上而下分为侏罗系直罗组、侏罗系延安组，富水性为弱-极弱。

1.2 水样采集与特征分析前期采集水样18个，并对各水样的水质特征进行了化验。第1～15号水样分别采自于地表水、第四系松散层、白垩系洛河组、侏罗系直罗组、侏罗系延安组;第16～18号水样采自于采空区积水。主要检测指标包括：K+、Na+、Ca2+、Mg2+、Cl-、SO2-4、HCO-3。各水样水质特征见表1。

2 模糊綜合评判法判别水源模糊综合评判法是一种基于模糊数学的综合评判方法，该方法根据模糊数学的隶属度理论把定性评判转化为定量评判，即用模糊数学对受到多种因素制约的事物或对象做出一个总体的评判[19-20]。

选取采集的各个水样阴阳例子浓度指标作为分类因素进行综合评判，并建立影响因素集合U={Ca2+，Mg2+，Na++K+，SO2-4，HCO-3，Cl-};所以选取地表水、第四系水以及各含水层水样阴阳例子浓度平均值作为对应水源的背景值。因侏罗系延安组与侏罗系直罗组含水层距离采空区相对更近，因此建立评判集合B={洛河组水Ⅰ，直罗组水Ⅱ，延安组水Ⅲ}，V={第四系水Ⅰ，直罗组水Ⅱ，延安组水Ⅲ}，F={地表水Ⅰ，直罗组水Ⅱ，延安组水Ⅲ}。因计算过程相同，就以评判集合B为例进行计算，并选取16、17、18号采空区水样作为因素指标实测值（表2）。

2.1 建立权重模糊矩阵

E={e1，e2，e3，e4，e5，e6}，其中ei=（xi/ai）/

∑6i=1

xiai

（∑6i=1ei=1，

ei>0）

（1）

式中 ei为权重向量;xi为采空区积水样中第i个因素指标实测值;ai为对应的第i个因素指标在3个已知含水层中的平均值。依据以上公式并以16号采空区水样为例得出E1={0.10，0.05，0.22，006，0.09，0.48}。

2.2 确定隶属度函数依据已知信息，采用专家推理法，确定隶属度函数为梯形模糊分布函数。以16号水样中的Ca2+为例进行计算，其隶属度函数r采用降半阶梯型分布函数。

r=

1x≤50.17

61.08-x61.08-50.1750.17

0x≥61.08

（2）

式中 x为采空区积水源因素指标实测值;r为Ca2+对Ⅰ、Ⅱ、Ⅲ水源隶属度。计算结果以矩阵形式表达。以采空区16号水样为例，计算出Rca=（0，0.970 0，0.030 0）。同理，分别计算出Mg2+、Na++K+、SO2-4、HCO-3、Cl-作为影响因素的评判结果，进而求出模糊评价矩阵R1。

R1

=

00.870 90.129 1

00.970 00.030 0

001

00.868 20.131 8

0.382 50.617 50

001

2.3 模糊矩阵复合运算通过模糊矩阵复合运算得到16号水样的综合评判结果为

B1=E1R1=（0.03，0.24，0.72）同理分别求出17、18号采空区积水水样的权重模糊矩阵E2和E3以及模糊评价矩阵R2和R3，得出综合评判结果B2和B3。16、17、18号采空区积水水样的评判结果见表3。

2.4 评判结果评判结果B1中0.72>0.24>0.03，其中003、0.24、0.72分别代表白垩系洛河组、侏罗系直罗组、侏罗系延安组水样对于采空区水样的隶属度，依据最大隶属度原则分析得出：16号水样来源为侏罗系延安组。同理分析B2和B3可得17号、18号水样来源均为侏罗系延安组，即采空区积水来自侏罗系延安组。各计算结果表明侏罗系延安组水样的隶属度最大，而侏罗系直罗组次之。但侏罗系直罗组水是否作为采空区积水来源之一，模糊综合评判法因模糊性的干扰无法进行有效的判别，因此下文将进行进一步的验证，并结合导水通道研究综合判别水源。

3 水化学特征分析法判别水源含水层的岩性特征、水文地质条件以及赋存环境等对地下水化学成分的形成产生重要影响[21-23]。不同含水层因围岩性质、循环条件等因素的作用，具有不同的水化学特征，通过分析各含水层的水化学特征并与采空区积水进行比较，进而判别水源[24-25]。

3.1 各含水层水化学特征分析

3.1.1 地表水水化学特征地表水属于弱碱性水，矿化度较低。水质类型为HCO3-Mg·Ca、HCO3·SO4-Mg·Ca型水，与大气降水具有相同的水质。

3.1.2 第四系含水层水水化学特征第四系含水层水为弱碱性，低矿化度。水化学类型为HCO3-Mg·Ca、HCO3-Na·Ca，水化学类型与地表水接近。

3.1.3 白垩系洛河组含水层水化学特征白垩系洛河组水质总体表现为低矿化度、弱碱性、重碳酸钙型水特征。水质类型以HCO3·SO4-Na·Ca（Mg）、HCO3-Mg·Ca为主。该含水层与地表水、第四系松散层水力联系非常密切，水循环较快，可以将地表水、第四系和白垩系含水层视为一个整体的浅部含水层。

3.1.4 侏罗系直罗组含水层水化学特征

侏罗系直罗组含水层水质呈弱碱性，矿化度较高，含水层水化学类型以SO4-Na、Cl-Na为主，其水化学特征与第四系和白垩系含水层存在显著区别，说明本地区浅部含水层与深部含水层在地下水运移和水文地球化学演化等方面存在较大差异。

3.1.5 侏罗系延安组含水层水化学特征侏罗系延安组水质呈弱碱性，水质类型为SO4-Na、Cl-Na型。延安组水的矿化度在矿区含水层中最高。长期的滞流型地下水运移和水文地球化学演化过程中，形成了高矿化度、弱碱性、硫酸钠（氯化钠）型水。

3.1.6 采空区积水水化学特征采空区积水的水化学类型以SO4·Cl-Na、Cl·SO4-Na为主。水化学特征表明，采空区积水主要受深部含水层岩性、地下水运移、水岩作用综合影响，以及顶板水进入采空区后，与岩石垮落物发生的水岩作用，导致采空区储水形成了高矿化度、硫酸钠（氯化钠）型水，与浅部含水层（包括第四系和白垩系）存在较大差异，水质与侏罗系直罗组、侏罗系延安组接近。

3.2 采空区积水水源分析依据表1所有水样数据，利用Aquachem软件绘制水化学背景Piper三线图和水源判别Piper三线图，如图1和图2所示。

从图1可分析出地表水主要处于菱形的左侧，白垩系洛河组水和第四系水处于左中部，水质类型主要呈现出HCO3-Ca（Mg）的特点，可与地表水视为同一浅部含水层。而直罗组和延安组含水层水样处于菱形的右侧，水质类型呈现出SO4-Na型和Cl-Na型的特点。图2直观显示出采空区积水样位于菱形右端，水质类型为SO4-Na型和Cl-Na型，与侏罗系直罗组和侏罗系延安组含水层水样相似。从两图分析可以得出：采空区积水水源来自侏罗系直罗组和延安组。采空区积水水化学特征与浅部含水层（地表水、第四系、白垩系）存在非常大的差异，与侏罗系直罗组和侏罗系延安组水化学特征相近，表明采空区积水水源以侏罗系直罗组和延安组含水层水为主。

4 导水通道判断矿井充水分别受大气降水、地表水和地下水等因素的影响，通过以上模糊综合评判法和水化学特征分析得出，采空区积水水源与侏罗系的直罗组和延安组含水层有关，与浅部的地表水、第四系水、白垩系差异很大。这就排除了大气降水和地表水作为矿井充水的可能性。调研地质资料可知，该矿未发现落差较大的断层、岩溶陷落柱以及封闭质量不良的钻孔等。由此推断，该矿区采空区积水的主要导水通道是采空区上部的导水裂隙带。根据《矿区水文地质工程地质勘探规范》（GB12719-91）中的冒落带、导水裂隙带最大高度经验公式计算可得，导水裂隙带最大高度为7.35～105.80 m，一般72 m。根据该矿的钻孔资料，绝大部分钻孔导水裂隙带位于侏罗系延安组，仅西南部的4个钻孔导水裂隙带发育至侏罗系直罗组下段（表4），突破高度分别为6.30 m、6.18 m、517 m、6.63 m。

5 結 论1）模糊综合评判法得出采空区积水来自侏罗系延安组含水层。2）水化学特征分析结果发现采空区积水水化学特征与侏罗系直罗组和延安组非常接近。3）导水裂隙带是采空区积水的主要导水通道，导水裂隙带主要在侏罗系延安组发育，少数区段发育至侏罗系直罗组底部。4）根据模糊综合评判法以及水化学特征的分析结果，结合导水裂隙带发育高度，能准确地判断出采空区积水水源主要来自侏罗系延安组含水层，少部分来自直罗组含水层。

参考文献（References）：

[1] 王心义，赵伟，刘小满，等.基于熵权-模糊可变集理论的煤矿井突水水源识别[J].煤炭学报，2017，42（9）：2433-2439.

WANG Xinyi，ZHAO Wei，LIU Xiaoman，et al.Identification of water inrush source from coalfield based on entropy weightfuzzy variable set theory[J].Journal of China Coal Society，2017，42（9）：2433-2439.

[2]任鄧君，孙亚岳，李建阳.高家堡煤矿煤层顶板水水化学特征及其水害防治技术[J].煤田地质与勘探，2019，47（S1）：26-31.

REN Dengjun，SUN Yayue，LI Jianyang.Hydrochemical characteristics and control of water hazard from coal seam roof in Gaojiabao coal mine[J].Coal Geology & Exploration，2019，47（S1）：26-31.

[3]张淑莹，胡友彪，琚棋定.基于水化学特征分析判别朱集矿矿井突水水源[J].矿业安全与环保，2018，45（6）：53-56.

ZHANG Shuying，HU Youbiao，JU Qiding.Distinguishing the sources of water inrush in Zhuji coal mine by hydrochemical characteristics[J].Mining Safety & Environmental Protection，2018，45（6）：53-56.

[4]杨婷婷，许光泉，余世滔，等.煤层下部太原组岩溶水化学组分特征及其成因分析[J].水文地质工程地质，2019，46（2）：100-108.

YANG Tingting，XU Guangquan，YU Shitao，et al.An analysis of the chemical composition characteristics and formationof the karst groundwater in the Taiyuan group in the lower part of a coal seam[J].Hydrogeology & Engineering Geology，2019，46（2）：100-108.

[5]张雁.水化学分析法判定矿井突水水源及示踪法验证[J].煤炭科学技术，2014，42（10）：98-100，105.

ZHANG Yan.Determination of mine water inrush source by hydrochemical analysis and verification by tracer method[J].Coal Science and Technology，2014，42（10）：98-100，105.

[6]GOLCHIN I，MOGHADDAM M A.Hydrogeochemical characteristics and groundwater quality assessment in Iranshahr plain aquifer，Iran[J].Environmental Earth Sciences，2016，75（4）：1-14.

[7]QUANG T T，ANDRE B，FRANK W K，et al.Mine water hydrogeochemistry of abandoned coal mines in the outcropped Carboniferous formations，Ruhr Area，Germany[J].Environmental Earth Sciences，2020，79（4）.

[8]杨建.呼吉尔特矿区葫芦素煤矿水文地球化学特征研究[J].煤矿安全，2016，47（12）：203-206.

YANG Jian.Research on hydrogeochemical characteristics in Hulusu coal mine of Hujierte coal field[J].Safety in Coal Mines，2016，47（12）：203-206.

[9]郭瑞，王永申，许刚.不连沟煤矿充水水源水化学特征分析及突水水源判别模型研究及应用[J].煤炭科技，2018（3）：128-131.

GUO Rui，WANG Yongshen，XU Gang.Chemicalcharacteristics analysis of waterfilled source water in Buliangou coal mine and research and application of water source water source discrimination model[J].Coal Science & Technology Magazine，2018（3）：128-131.

[10]VERMEULEN P D，BURGER M，VAN WYK A.Potential hydrological interaction of a gold mine and a coal mine in south africa[J].Mine Water and the Environment，2014，33：3-14.

[11]ZHANG H T，XU G Q，WEI J，et al.Hydrogeochemical characteristics and groundwater inrush source identification for a multiaquifer system in a coalmine[J].Acta Geologica Sinica English Edition，2019，93（6）：1922-1932.

[12]DELIA C P，COCIUBA I，

ALEXANDRA C，et al.Origin and geochemistry of mine water and its impact on the groundwater and surface running water in postmining environments：Zlatna Gold Mining Area（Romania）[J].Aquatic Geochemistry，2017，23：247-270.

[13]魏大勇，王飞，许进鹏，等.基于分形与模糊综合评价法的矿井突水危险性评价[J].煤矿安全，2013，44（8）：184-186.

WEI Dayong，WANG Fei，XU Jinpeng，et al.Risk assessment of mine water inrush based on Fractal and fuzzy comprehensive evaluation[J].Safety in Coal Mines，2013，44（8）：184-186.

[14]GUAN Z L，JIA Z F，YOU Q Y，et al.Identification of inrush water recharge sources using hydrochemistry and stable isotopes：A case study of Mindong No.1 coal mine in northeast inner Mongolia，China[J].Journal of Earth System Science，2019，128（7）.

[15]姬亚东.基于聚类分析与模糊综合评判的煤层顶板涌水危险性评价[J].矿业安全与环保，2019，46（4）：68-72.

JI Yadong.The risk assessment of roofwater inrush based on cluster analysis andfuzzy comprehensive evaluation[J].MiningSafety & Environmental Protection，2019，46（4）：68-72.

[16]李凤莲，冯琳，张雪英，等.模糊综合评判法的改进及在水源判别中的应用[J].太原理工大学学报，2015，46（4）：444-447.

LI Fenglian，FENG Lin，ZHANG Xueying，et al.Improvement of the fuzzy comprehensive evaluation and its application in waterbursting source discrimination[J].Journal of Taiyuan University of Technology，2015，46（4）：444-447.

[17]張红梅，吴基文，翟晓荣，等.煤系砂岩含水层富水性模糊综合预测与评价——以卧龙湖煤矿七含为例[J].矿业安全与环保，2018，45（4）：70-74，79.

ZHANG Hongmei，WU Jiwen，ZHAI Xiaorong，et al.Fuzzy comprehensive prediction and evaluation of the water abundance of coalmeasures sandstone aquifer：Taking Qihan of Wolonghu coal mine for example[J].Mining Safety & Environmental Protection，2018，45（4）：70-74，79.

[18]代革联，薛小渊，牛超.基于水化学特征分析的象山矿井突水水源判别[J].西安科技大学学报，2017，37（2）：213-218.

DAI Gelian，XUE Xiaoyuan，NIU Chao.Discrimination of water inrush source in Xiangshan coal mine based on chemical characteristics of groundwater[J].Journal of Xian University of Science and Technology，2017，37（2）：213-218.

[19]马占新，斯琴.基于模糊综合评判的改进策略及优化方法[J].系统工程与电子技术，2018，40（9）：2016-2025.

MA Zhanxin，SI Qin.Optimization methodand improving strategy based on fuzzycomprehensive evaluation[J].Systems Engineering and Electronics，2018，40（9）：2016-2025.

[20]刘秀娥，韩科明，樊振丽.基于熵权模糊综合评价的矿井突水水源判别方法[J].煤矿开采，2017，22（6）：82-84，9.

LIU Xiue，HAN Keming，FAN Zhenli.Discriminated method of mine water inrush source based on entropy weight fuzzy comprehensive analysis[J].Coal Mining Technology，2017，22（6）：82-84，9.

[21]刘剑民，王继仁，刘银朋，等.基于水化学分析的煤矿矿井突水水源判别[J].安全与环境学报，2015，15（1）：31-35.

LIU Jianmin，WANG Jiren，LIU Yinpeng，et al.Coal mine based on hydrochemical analysis discrimination of mine water inrush source[J].Journal of Safety and Environment，2015，15（1）：31-35.

[22]贾亚琳，边凯，李冲，等.基于GIS和相关性分析峰峰矿区岩溶水水化学特征[J].中国煤炭，2019，45（6）：95-101.

JIA Yalin，BIAN Kai，LI Chong.Based on GIS and correlation analysis of karst water chemical characteristics in Fengfeng mining area[J].China Coal，2019，45（6）：95-101.

[23]陈凯，孙林华.任楼煤矿地下水化学组成及其控制因素分析[J].煤炭科学技术，2019，47（10）：240-244.

CHEN Kai，SUN Linhua.Analysis of chemical composition and control factors of groundwater in Renlou coal mine[J].Coal Science and Technology，2019，47（10）：240-244.

[24]劉洋，方刚，杨建.准格尔矿区酸刺沟煤矿水文地球化学特征研究[J].煤矿安全，2017，48（7）：204-207.

LIU Yang，FANG Gang，YANG Jian.Hydrogeoc hemical characteristics of Suancigou coal mine in Zhungeer mining area[J].Safety in Coal Minesy，2017，48（7）：204-207.

[25]杨志斌，边凯.基于特定水化学特征的矿井涌水量构成比例分析[J].西安科技大学学报，2015，35（6）：727-731.

YANG Zhibin，BIAN Kai.Analysis of mine water inflow composition proportion based on the specific water chemical characteristics[J].Journal of Xian Universityof Science and Technology，2015，35（6）：727-731.