宁正煤田侏罗系煤层顶板富水性及涌（突）水风险评价

李奇峰 崔恒 肖乐乐 刘凯 许鑫

摘 要：为研究宁正煤田侏罗系煤层顶板含水层的富水性，运用主客观因素综合赋权的方法，选取含水层厚度、粗砂岩含量、冲洗液消耗量、岩芯采取率、砂泥岩互层数、泥岩含量作为富水性评价的指标，对8煤层的顶板富水性进行分析。结果表明：利用模糊层次分析法-熵权法综合赋权得到各项指标的综合权重，结合逼近理想解排序法的基本原理建立了煤层顶板富水性评价模型，南部地区涌水量约为0.2L/S，中部和北部涌水量约为0.27L/S，西北部涌水量约为0.36L/S，西南部涌水量约为0.43L/S，在煤田内涌水量最大。并通过已知涌水点的分布，对所建立的评价模型进行检验，验证了该模型的准确性。宁正煤田8煤层顶板富水性在西部和东部部分地区N703、N109钻孔周边的富水性最强，在南部补406、补601周边的富水性最弱。

关键词：风险评价；模糊层次分析法；熵权法；宁正煤田

中图分类号：TD 745文献标志码：A

文章编号：1672-9315（2023）03-0523-07

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

Evaluation of water-abundance and water gushing risk of Jurassic coal seam roof in Ningzheng coal field

LI Qifeng¹，CUI Heng²，XIAO Lele³，LIU Kai³，XU Xin³

（1.146 Team，Gansu Coal Geology Bureau，Pingliang 744000，China；2.Binxian Shuiliandong Coal Co.，Ltd.，Xianyang 712000，China；3.College of Geology and Environment，Xian University of Science and Technology，Xian 710054，China）

Abstract：This research attempts to investigate the water richness of the roof aquifer of Jurassic coal seam in Ningzheng coal field and evaluate the water richness of the roof aquifer accurately.By combining subjective and objective factors with comprehensive weighting method，aquifer thickness，coarse-grained sandstone content，flushing fluid consumption，drilling core rate，sand and mudstone interbedding number，mudstone content  are selected as indicators to evaluate the water richness of No.8 coal roof.The results show that：The comprehensive weight of each index is obtained by using fuzzy analytic hierarchy process and entropy weight method.Based on the fundamented principle of approximate ideal solution sorting method，a water-rich evaluation model of coal seam roof is established.The water inflow in southern area is about 0.2L/S，the water inflow in central and northern areas is about 0.27L/S，and the water inflow in northwestern area is about 0.36L/S.The water inflow in the southwest is about 0.43L/S，which is the largest in the coal field.Through the distribution of known inrush points，the evaluation model is the tested and the accuracy of the model is verified.The water-rich area of No.8 coal seam roof in Ningzheng coalfield is the strongest around N703 and N109 boreholes in western and eastern parts.In the south，the water richness is weakest around drilling holes 406 and 601.

Key words：risk assessment；fuzzy analytic hierarchy process；entropy weight method；Ningzheng coal field

0 引 言

煤矿水害是我国煤矿主要五种自然灾害之一，破坏性巨大，给煤矿安全生产带来了很大隐患和损失^［1］。曾一凡等从影响含水层富水性的岩性差异、水力特性、构造因素、地球物理参数4个因素出发，构建了能够较为全面，真实地刻画富水性分布规律的四大类，16小类的主控因素体系^[2]；周全超等提出利用CRITIC法与AHP法分别得到评价指标权重，运用拉格朗日乘子法进行耦合，最后绘制了研究区含水层富水性分区图^[3]；汪子涛等结合熵权法和TOPSIS分区，建立了某盆地的涌水量预测模型，对研究区地下水空间富水性进行了分区^[4]；唐李斌等建立了AHP-熵权法耦合的分界砂岩含水层富水性评价模型，将研究区分界砂岩含水层划分为强富水区、较强富水区、中等富水区、较弱富水区和弱富水区５个等级^[5]；毕尧山等通过主、客观有机结合得到综合权重，构建了基于层次分析法与独立性权系数法综合确权的含水层富水性指数法评价模型，利用ArcGIS绘制了研究区顶板砂岩含水层富水性分区图^[6]；武强等采用GIS技术，提出将渗透系数，砂岩厚度，冲洗液消耗量，岩心采取率以及脆塑性岩厚度比作为顶板上的涌（突）水的主控因素，构建了含水层富水性分区，同时结合多源地学信息复合叠加原理，提出9煤层顶板冒落涌（突）水综合分区的划分方案^[7-8]；李哲等采用BP人工神经网络方法，选取了砂岩含水层厚度、岩芯采取率、脆性岩厚度比和风化影响指数４个指标，引入具有较好自主学习、非线性映射能力的BP人工神经网络，使用训练好的神经网络对研究区含水层富水性进行预测^[9]；张红梅等采取岩性结构指数、断裂构造分维值和褶皱变形系数３个定量指标，结合模糊聚类综合评判方法分析了含水层富水性，基于岩芯采取率和冲洗液消耗量等水文地质信息，验证了模糊聚类综合评判结果的可靠性^[10]；韩承豪等利用集对分析-可变模糊集耦合的评价方法，对直罗组砂岩含水层富水性进行预测，建立顶板砂岩的富水性评价模型^[11]；肖樂乐等以富水结构指数为基础，选取含水层等效厚度、岩心采取率、砂泥比以及砂泥岩交互层数4个影响因素作为独立评价指标，提出了侏罗系煤田顶板直罗组砂岩含水层富水评价方法^[12]；李世峰等通过研究浓缩因子的分区原则，结合矿区富集度，构建了矿区地下水富水性进行分区^[13]；卜庆林等结合3#煤顶板中砂岩厚度、冲洗液最大消耗量、岩心采取率、平均孔隙度和断层密度5个因素的影响作用，利用信息拟合方法，构建了煤层顶板富水性分区^[14]；石守桥等采用断层分维值对断层构造特征进行定量描述，构建了岩性结构指标模型，对煤层顶板富水进行分区^[15]；根据模糊层次、模糊聚类^[16-17]等数据方法的分析，对煤层顶板的富水性进行分区，提出了各影响因子权重的确定问题。

笔者采用模糊层次分析法-熵权法综合赋权-逼近理想解排序分区的方法，对宁正煤田侏罗系8号煤层的顶板富水性进行分区，以期达到对煤田内顶板富水性强弱有具体的理论认知的目的。

1 地质概况

宁正煤田包括5个含水层和3个隔水层，其中第四系全新统砂砾层孔隙潜水含水层为井田第一含水层，下白垩统志丹群碎屑岩类孔隙裂隙潜水含水层为井田第二含水层，下白垩统志丹群孔隙、裂隙承压含水层为井田第三含水层，中侏罗统直罗组、延安组上、中部（8煤层顶板以上）砂岩复合承压含水层为井田第四含水层，中侏罗统延安组下部8煤层底板以下—三叠系砂岩、砾岩孔隙、裂隙复合承压含水层为井田第五含水层（图1）。

中侏罗统直罗组、延安组上、中部砂岩复合承压含水层该层属煤层顶板含水层，主要为细砂岩、中砂岩、粗砂岩，一般直罗组下部及延安组每层煤的间接顶板都为中—粗粒砂岩，含水层厚0.7～140 m，平均60.11 m。

宁正煤田内含水层自上而下主要有第四系孔隙潜水含水层、第四系中更新统黄土孔隙-裂隙潜水含水层、新近系砂砾岩含水层段、白垩系下统华池组砂岩裂隙含水层段、白垩系下统洛河组砂岩孔隙-裂隙承压含水层段、白垩系下统宜君组砾岩裂隙含水层段、侏罗系中统直罗组砂岩裂隙承压含水层、侏罗系中统延安组煤层及顶板砂岩承压含水层和三叠系上统胡家村组砂岩裂隙含水层。其中侏罗系中统直罗组砂岩裂隙承压含水层、侏罗系中统延安组煤层及其顶板砂岩承压含水层是开采8煤层的直接充水含水层。

2 顶板富水性权重确定

2.1 主观权重

模糊层次分析法是在层次分析法的基础上添加模糊数学理论，对目标指标进行评价以减少个人主观影响^[18]。层次分析法是一种多因素分析、多指标分级的评价方法^[19]。层次分析法的关键在于构造各层次的判断矩阵，将模糊数学思想和方法引入^[20-21]，得到综合考虑人为判断的模糊性的定性与定量相结合的系统分析方法，即模糊层次分析法。

2.2 客观权重

熵权法能反应出信息无序化的程度，能最大化地降低权重计算中人为因素的干扰^[22]。熵权法是根据各主控因素之间变化幅度的大小来确定客观权重，熵是系统无序程度的一个度量^[23]，能够较好地反映各主控因素之间的内部客观逻辑，信息熵值越小，指标的离散程度越小。

2.3 综合权重

由于主客观权重的赋权依据不同，为兼顾主客观权重的优点，需要对主客观因素进行综合赋权。当评价指标个数较多，权数分配比较均匀时往往采用乘法合成归一化法将主客观权重有机结合^[24]，计算主控因素的综合权重。

3 顶板富水性风险因素

3.1 含水層厚度

含水层厚度直接影响含水层的储水量，在一定程度上能够代表含水层储水空间的大小，是确定含水层富水性的重要指标之一。含水层较厚的地方主要集中在井田的西北部（图2（a））。在N701、N703、N602附近含水层厚度为189～360 m；在N005、N008、补402附近含水层厚度在118～189 m；在N801、N603、N801、N111附近含水层厚度在65～118 m；在N506、N406、N304、N205、N306、补406附近含水层厚度为0～65 m。

3.2 粗砂岩含量

粗砂岩含量主要指岩石粗砂岩的含量，粗砂岩含量越多，则涌入采动空间的水量越大，富水性越强。粗砂岩厚度较大的区域主要在井田西南部（图2（b））。在N703的附近138～240 m；在N701附近，粗砂岩含量较大，粗砂岩含量在88～138 m；在N601、N506、N406、N304、N306、N205、补101、补402附近，粗砂岩为0～40 m。

3.3 砂泥岩互层数

砂泥岩互层数是指煤层顶板含水层与隔水层的交互层数，含水层、隔水层互层层数越多，富水性越弱。砂泥岩互层数较多的区域集中在西北部和东部（图2（c））。在N802、N106、N102的附近，砂泥互层数为42～75 m；在N701、N602、N703、补402的附近，砂泥岩互层数为29～42层；在N601、N506、N406、N304、N112的附近，砂泥岩互层数为20～29层；在N801、N306、N205、补406附近，砂泥岩互层数为10～20层。

3.4 岩芯采取率

岩芯采取率可以用来表示岩石的完整程度，岩石越破碎，取芯率越低，连通性越好，富水性越好，井田西南部的取芯率较高（图2（d））。在N506、N406、补406、N010的附近，岩芯采取率为80%～92%；在N304、N306、N201、N112的附近，岩芯采取率为71%～80%；N701、N603、N701的附近，岩芯采取率为38%～63%。

3.5 冲洗液消耗量

冲洗液消耗量能反应岩层透水性能，冲洗液消耗量大通常渗透系数大。冲洗液消耗量大的区域主要在井田北部及中部区域（图2（e））。在补403的附近，冲洗液消耗量为6.5～15 L·min^-1；在N304、N703、N108的附近，冲洗液消耗量为3～6.5  L·min^-1；在N111、N110、N801、N802的附近，冲洗液消耗量为0～1.4  L·min^-1。

3.6 泥岩含量

泥岩作为隔水层可以阻止含水层的水进入矿井，泥岩的厚度直接影响着其隔水能力，泥岩含量较大区域主要在井田的西南部（图2（f））。在N802、N701、N602、N111的附近，泥岩含量为93～170 m；在N005、N801、N108、N004的附近，泥岩含量为57～93 m；在N304、N506、N406、N019、N112的附近，泥岩含量为27～57 m；在N011、N006、补406钻孔的附近，泥岩含量为0～27 m。

4 富水性评价体系

4.1 评价模型

为建立正确合理的顶板富水性评价模型，将顶板富水性评价模型分为3个层次结构，目标层为煤层顶板富水性影响因素，准则层包括有富水介质类型、储水空间大小和水力联系强度，决策层对准则层每一步条件进行细化，充分考虑每一个影响因素，包括有含水层厚度、粗砂岩含量、冲洗液消耗量、岩芯采取率、泥岩含量和砂泥互层数（图3）。

4.2 主控因素指标权重的确定

4.2.1 模糊层次分析法确定主观权重

根据煤层顶板富水主控指标建立模糊一致判断矩阵。利用0.1～0.9模糊标度法，结合宁正煤田自身地质条件，构建富水预测的模糊互补判断矩阵。从富水介质类型来看水力联系强度比储水空间大小对富水性的影响更大，从储水空间大小角度来看水力联系强度比富水介质类型对富水性的影响更大，从水力联系强度来看储水空间大小比富水介质类型对富水性影响更大（表1～表4）。

4.2.2 熵权法确定客观权重

熵权法不参杂人为主观因素，熵值越大，富水性影响因素的占比越重，越重要；熵值越小，富水性因素占比越轻，越不重要。熵权法所得到的客观权重冲洗液消耗量的占比最大，泥岩含量占比的最小（表5）。

4.2.3 模糊层次分析法-熵权法确定综合权重

模糊层次分析法是含水层富水性分析的一种方法，熵权法是不考虑专家意见而进行的客观打分，采用综合权重可以有效的规避单独使用二者之一的缺陷。冲洗液消耗量的综合权重最大，最重要；砂泥岩互层数的综合权重最小，对富水性的影响最小。岩芯采取率值最小，岩芯采取率对富水性影响最弱，冲洗液消耗量值最大，对富水性影响最大（表5）。

4.3 宁正煤田顶板富水性风险分区

宁正煤田已知有56个钻孔数据，结合宁正煤田含水层富水性影响因素权重（表5），利用逼近理想解排序法对宁正煤田顶板富水性进行风险区，得出各钻孔涌水风险系数为0.2～0.6，涌水风险系数越大，富水性越强；涌水风险系数越小，富水性越弱。宁正煤田内8煤层顶板富水风险主要在西北部和东部地区，南部地区富水风险极低，北部地区风险中等，钻孔涌水量大小与富水性分区一致（图4）。

5 结 论

1）提出宁正煤田8煤层顶板富水的影响因素，主要有含水层厚度、粗砂岩含量、冲洗液消耗量、岩芯采取率、砂泥岩互层数、泥岩含量。

2）利用模糊综合评价-熵权法综合赋权，逼近理想解排序法分区的方法对宁正煤田进行顶板富水性评价，构建了煤层顶板富水性风险评价模型。

3）对宁正煤田的8煤层顶板富水性进行分区，南部地区涌水量最小，中部和北部涌水量較小，西南部在矿井内涌水量最大。

参考文献（References）：

[1] 曹主军.随钻测量定向钻进技术在矿井顶板水患治理中的应用研究[J].探矿工程，2015，42（2）：23-27.

CAO Zhujun.Application research on directional drilling with MWD in roof water control[J].Exploration Engineering，2015，42（2）：23-27.

[2]曾一凡，武强，杜鑫，等.再论含水层富水性评价的“富水性指数法”[J].煤炭学报，2020，45（7）：2423-2431.

ZENG Yifan，WU Qiang，DU Xin，et al.Further research on “water-richness index method”for evaluation of aquifer water abundance[J].Journal of China Coal Society，2020，45（7）：2423-2431.

[3]周全超，顾雷雨，殷裁云，等.基于AHP-CRITIC-TOPSIS法的含水层富水性评价模型[J].人民黄河，2022，44（S2）：116-118，128.

ZHOU Quanchao，GU Leiyu，YIN Caiyun，et al.Evaluation model of aquifer water abundance based on AHP-CRITIC-TOPSIS method[J].Peoples Yellow River，2022，44（S2）：116-118，128.

[4]汪子涛，王建萍，韩光.基于熵权-TOPSIS法的柴达木盆地地下水富水性评价[J].盐湖研究，2022，30（2）：42-51.

WANG Zitao，WANG Jianping，HAN Guang.Evaluation of groundwater water abundance in Qaidam Basin based on entropy weight-TOPSIS method[J].Salt Lake Study，2022，30（2）：42-51.

[5]唐李斌，吴基文，毕尧山，等.基于AHP-熵权法耦合的含水层富水性评价研究[J].中国矿业，2020，29（12）：147-152.

TANG Li Bin，WU Jiwen，BI Yaoshan，et al.Aquifer water abundance evaluation based on AHP-entropy weight method coupling[J].China Mining Industry，2020，29（12）：147-152.

[6]毕尧山，吴基文，翟晓荣，等.基于AHP与独立性权系数综合确权的煤矿含水层富水性评价[J].水文，2020，40（4）：40-45.

BI Yaoshan，WU Jiwen，Zhai Xiaorong，et al.Water abundance evaluation of coal mine aquifer based on AHP and independent weight coefficient[J].Hydrology， 2020，40（4）：40-45.

[7]武强，许珂，张维.再论煤层顶板涌（突）水危险性预测评价的“三图一双预测法”[J].煤炭学报，2016，41（6）：1341-1347.

WU Qiang，XU Ke，ZHANG Wei.Further research on “three maps-two predictions” method for prediction on coal seam roof water bursting risk[J].Journal of China Coal Society，2016，41（6）：1341-1347.

[8]武强，黄晓玲，董东林.评价煤层顶板涌（突）水条件的“三图-双预测法”[J].煤炭学报，2000，25（1）：60-65.

WU Qiang，HUANG Xiaoling，DONG Donglin.“Three maps-two predictions”method to evaluate water bursting conditions on roof coal[J].Journal of China Coal Society，2000，25（1）：60-65.

[9]李哲，曾一凡，刘守强，等.BP人工神经网络在富水性评价中的应用[J].煤炭工程，2018，50（8）：114-118.

LI Zhe，ZENG Yifan，LIU Shouqiang，et a1.Application of back propagation artificial neural network in water abundance evaluation[J].Coal Engineering，2018，50（8）：114-118.

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

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

[11]韩承豪，魏久传，谢道雷，等.基于集对分析-可变模糊集藕合法的砂岩含水层富水性评价——以宁东矿区金家渠井田侏罗系直罗组含水层为例[J].煤炭学报，2020，45（7）：2432-2443.

HAN Chenghao，WEI Jiuchuan，XIE Daolei，et al.Water-richness evaluation of sandstone aquifer based on set pair analysis-variable fuzzy set coupling method A case from Jurassic Zhiluo Formation of Jinjiaqu coal mine in Ningdong mining area[J].Journal of China Coal Society，2020，45（7）：2432-2443.

[12]肖乐乐，牛超，代革联，等.基于富水性结构指数法的直罗组地层富水性评价[J].煤炭科学技术，2018，46（11）：207-213.

XIAO Lele，NIU Chao，DAI Celian，et al.Evaluation of water abundance in Zhiluo formation based on watery structure index method[J].Coal Science and Technology，2018，46（11）：207-213.

[13]李世峰，陈志斌.利用浓缩因子区划矿井富水性[J].煤炭工程，2005（1）：43-45.

LI Shifeng，CHEN Zhibin.Concentrated factors applied to identify mine rich in water[J].Coal Engineering，2005（1）：43-45.

[14]卜庆林，陈成星，杨成超，等.煤层顶板岩层富水性分区指标及其涌水量预侧[J].山东科技大学学报（自然科学版），2005，24（3）：28-31.

BU Qinglin，CHEN Chengxing，YANG Chengchao，et al.Zoning index of rich water bearing roof strata and water in rush forecasting in coal seams[J].Journal of Shandong University of Science and Technology（Natural Science），2005，24（3）：28-31.

[15]石守桥，魏久传，尹会永，等.济三煤矿煤层顶板砂岩含水层富水性预测[J].煤田地质与勘探，2017，45（5）：100-104.

SHI Shouqiao，WEI Jiuchuan，YIN Huiyong，et al.Forecasting of water abundance of coal roof sandstone aquifer in Jining No.3 mine[J].Coal Geology & Exploration，2017，45（5）：100-104.

[16]张吉军.模糊层次分析法（FAHP）[J].模糊系统与数学，2000，14（2）：80-88.

ZHANG Jijun.Fuzzy analytical hierarchy process[J].Fuzzy Systems and Mathematics，2000，14（2）：80-88.

[17]韩超，伴晓华，李国梁，等.基于GIS多源信息集成的含水层富水性模糊层次分析法[J].水文地质工程地质，2012，39（4）：19-25.

HAN Chao，BAN Xiaohua，LI Cuoliang，et al.The fuzzy analytic hierarchy process of water abundance of an aquifer based on CIS and multiource information fusion techniques[J].Hydrogeology & Engineering Geology，2012，39（4）：19-25.

[18]郭金玉，张忠彬，孙庆云.层次分析法的研究与应用[J].中国安全科学学报，2008，18（5）：148-153.

GUO Jinyu，ZHANG Zhongbin，SUN Qingyun.Study and applications of analytic hierarchy process[J].China Safety Science Journal，2008，18（5）：148-153.

[19]张桂欣，孙柏涛.基于模糊层次分析的建筑物单体震害预测方法研究[J].工程力学，2018，35（12）：185-193，202.

ZHANG Guixin，SUN Baitao.Seismic damage prediction for a single building based on a fuzzy analytical hierarchy approach[J].Engineering Mechanics，2018，35（12）：185-193，202.

[20]李博.GRA-FAHP模型的煤層底板突水危险性评价[J].地质论评，2015，61（5）：1128-1134.

LI Bo.Risk assessment model of coal floor water-irruption based on GREA-FAHP[J].Geological Review，2015，61（5）：1128-1134.

[21]姬东朝，宋笔锋，喻天翔.模糊层次分析法及其在设计方案选优中的应用[J].系统工程与电子技术，2006（11）：1692-1694，1755.

JI Dongchao，SONG Bifeng，YU Tianxiang.FAHP and its application in the selection of design scheme[J].Systems Engineering and Electronics，2006（11）：1692-1694，1755.

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

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

[23]王汉斌.基于ArcGIS的首旺煤矿突水水害危险性评价[D].北京：中国地质大学，2017.

WANG Hanbin.Risk assessment of water inrush from coal mine of the Shouwang colliery based on ArcGIS[D].Beijing：China University of Mining and Technology，2017.

[24]刘守强，武强，曾一凡，等.基于GIS的改进AHP型脆弱性指数法[J].地球科学，2017，42（4）：625-633.

LIU Shuqiang，WU Qiang，ZENG Yifan，et al.The improved AHP vulnerable index method based on GIS[J].Earth Science，2017，42（4）：625-633.

（责任编辑：李克永）