• (010)84657853 84658665 84657900
                              • mt@www.larryacraig.com
                              • 北京市朝阳区芍药居35号中煤信息大厦(100029)

                              ★ 煤矿安全 ★

                              宁东煤田鸳鸯湖矿区煤层顶板水害防控技术与应用

                              赵宝峰1,2,黄忠正3,宗伟琴3

                              (1.中煤科工集团西安研究院有限公司,陕西省西安市,710054;2.陕西省煤矿水害防治技术重点实验室,陕西省西安市,710177;3.国家能源集团宁夏煤业有限责任公司,宁夏回族自治区银川市,750011)

                              摘 要 为了避免宁东煤田鸳鸯湖矿区矿井发生顶板水害事故,对浅部煤层顶板含水层进行了富水性评价、涌水量预测和顶板水疏放效果评价和系统研究。提出了基于地质钻孔水文信息的含水层富水性分区方法,选取影响直罗组下段含水层富水性的6个评价指标,分析各评价指标与钻孔单位涌水量之间的相关关系,借助语气算子比较法确定各指标的权重,利用三角形隶属度函数构建了模糊评判矩阵,采用模糊综合评判法对直罗组下段含水层富水性进行了分区。充分考虑了顶板砂岩含水层的水文地质特征,提出了静储量结合动态补给量的工作面涌水量预测方法;根据矿井涌水量组成和变化特点,提出了基于时空分区的矿井涌水量动态预测方法;利用顶板水疏放的水量观测值,结合工作面静储量和动态补给量预测值,提出了工作面疏放水效果“双指标四因素”定量化评价方法。结果表明:基于地质钻孔水文信息的含水层富水性分区方法、静储量结合动态补给量的工作面涌水量预测方法,基于时空分区的矿井涌水量预测方法和顶板水疏放效果定量化评价方法可以对顶板水害起到有效的防控作用。

                              关键词 顶板水害;富水性分区;涌水量预测;疏放水效果评价;钻孔单位涌水量

                              煤炭在我国一次能源结构中占60%左右,在未来几十年内仍将是我国的主体能源。随着我国东部地区煤炭资源逐渐枯竭,煤炭生产重点逐步向西转移,西部侏罗纪煤炭资源开发已经成为我国煤炭工业可持续发展和国民经济建设的重要支撑[1-2]。侏罗纪煤田主要面临煤层顶板水害威胁,许多学者围绕顶板水害防控开展了大量研究,在含水层富水性分区方面主要采用物探[3-4]、富水性指数[5-6]、综合探查法[7-8]等,涌水量预测方面主要采用比拟法[9]、解析法[10]和数值法[11]等,以上研究成果对于顶板水害防控发挥了重要作用。

                              宁东煤田是我国批准建设的14个亿吨级煤炭生产基地之一,同时也是宁东煤化工基地的供煤产地,包括鸳鸯湖、马家滩、碎石井、积家井等矿区,其探明储量达331亿t,占宁夏煤炭总储量的85%[12]。宁东煤田大部分矿井主采侏罗系延安组煤层,其上部普遍发育直罗组砂岩含水层,在大规模机械化开采条件下,工作面开采后对顶板覆岩扰动强度大,形成的导水通道容易沟通上覆含水层,造成工作面顶板水害事故。鸳鸯湖矿区内各矿井均面临不同程度的顶板水害威胁,需要解决的防治水难题包括顶板含水层科学合理的富水性分区、工作面与矿井涌水量动态预测和顶板水疏放效果定量化评价等,为保障矿井的安全生产,需要对顶板水害防控方法与技术进行系统研究。

                              1 鸳鸯湖矿区及其水文地质条件

                              1.1 鸳鸯湖矿区概况

                              1.1.1 基本条件

                              鸳鸯湖矿区位于宁夏回族自治区宁东煤田中部,总体呈南北条带状展布,鸳鸯湖矿区及各矿井平面位置如图1所示,由北向南分别为清水营、梅花井、石槽村、红柳和麦垛山煤矿,规划产能44 Mt/a。

                              1.1.2 地层与构造

                              根据钻孔揭露资料,研究区地层由上向下分别为第四系、古近系、白垩系、侏罗系和三叠系,其中侏罗系延安组为矿区的含煤地层,2、6、10、12和18号煤层全区发育,全部或大部分可采。研究区位于鄂尔多斯盆地西缘冲断、推覆构造带中北部,北部以褶皱为主,向南逐渐断层发育,破坏了褶皱的完整性。

                              1.2 鸳鸯湖矿区水文地质条件

                              1.2.1 含水层

                              研究区含水层可划分为孔隙潜水含水层和孔隙裂隙承压含水层,后者又可划分为白垩系、侏罗系和三叠系含水层,其中侏罗系含水层根据主要可采煤层可进一步划分为直罗组含水层、2~6号煤层、6~18号煤层和18号煤层以下含水层。

                              近年来,研究区主采2号煤层,其主要充水含水层为直罗组下段含水层,其平均厚度、钻孔单位涌水量和渗透系数的平均值由北向南逐渐增大,其水文地质条件也随之逐渐复杂。鸳鸯湖矿区及各矿位置如图1所示。鸳鸯湖矿区各矿井直罗组下段含水层概况见表1。

                              图1 鸳鸯湖矿区及各矿井平面位置

                              表1 鸳鸯湖矿区各矿井直罗组下段含水层概况

                              煤矿平均厚度/m单位涌水量/(L·s-1·m-1)渗透系数/(m·d-1)清水营47.780.009 90.009 6梅花井46.140.033 20.047 5石槽村67.780.058 00.059 4红柳68.990.060 50.077 4麦垛山100.700.148 70.301 6

                              1.2.2 隔水层

                              研究区隔水层包括古近系清水营组、侏罗系安定组、直罗组、延安组煤层顶底板隔水层,岩性主要为泥岩、粉砂岩和砂质泥岩。

                              1.2.3 地下水补径排

                              潜水含水层地下水主要接受大气降水补给,沿地势由高向低径流,在沟谷洼地以泉的形式排出,还有一部分越流补给承压含水层;承压含水层地下水主要接受露头区大气降水补给、潜水含水层的越流补给,径流方向大致与潜水含水层保持一致,最终以矿井排水的形式排出。

                              2 鸳鸯湖矿区顶板水害防控技术与应用

                              2.1 顶板含水层富水性分区

                              2.1.1 含水层富水性分区存在的问题

                              钻孔单位涌水量是含水层富水性评价的唯一直接指标,主要是根据抽水试验得到的。研究区各矿井在水文补勘阶段均施工了水文地质钻孔,在抽水试验后留设为地下水位长观孔,为了避免受到采掘活动的影响和破坏,水文地质钻孔一般位于采区边界外的保护煤柱内,但是对矿井防治水工作具有指导意义的却是采区内煤层顶板含水层富水性条件,怎样利用采区外水文地质钻孔的信息评价采区内的煤层顶板含水层富水性是需要解决的重要问题。通常井田范围内地质勘探钻孔较多,而水文地质钻孔较少,怎样有效利用大量地质勘探钻孔的水文信息对含水层富水性进行评价也是防治水面临的主要问题。

                              2.1.2 含水层富水性分区思路

                              建立采区外水文地质钻孔与采区内地质勘探钻孔的水文信息相关关系,进而利用地质勘探钻孔的水文信息对采区内含水层进行富水性分区成为首选思路。基于以上思路,提出了基于地质勘探钻孔水文信息的含水层富水性分区技术,即选取与含水层富水性相关性较强的评价指标,然后利用灰色关联度、语气算子比较法、隶属度函数和模糊综合评判法对煤层顶板含水层富水性进行分区[13]。基于地质勘探钻孔水文信息的含水层富水性分区流程如图2所示。

                              图2 基于地质勘探钻孔水文信息
                              的含水层富水性分区流程

                              2.1.3 含水层富水性分区步骤

                              (1)各评价指标与钻孔单位涌水量的关联度计算。灰色关联分析的具体计算步骤如下:设参考数列为:x0={x0(k)|k=1,2,…j},各水文孔的单位涌水量为参考数列;设比较数列为:x1={x1(k)|k=1,2,…j},x2={x2(k)|k=1,2,…j},…,xi={xi(k)|k=1,2,…j},含水层砂地比、含水层厚度、粗砂岩厚度、粗砂岩层数、断层分维值和褶皱分维值为比较数列。

                              关联系数表达式如下:

                              (1)

                              式中:ρ——分辨系数,通常取0.5。

                              关联度r0i作为比较数列与参考数列关联程度的数量表示,其计算公式如下:

                              (2)

                              根据式(2)计算的关联度大小排序,如果r0n<r0m,说明参考数列m比参考数列n与比较数列更为接近,对比较数列的影响更大。

                              (2)各评价指标权重的确定。在确定各评价指标的权重向量时,采用语气算子比较法,具体算法如下:设有6项含水层富水性评价指标组成的指标集:D=(d1,d2,…,d6),首先研究指标集D对重要性的二元比较定性排序。指标集D中的元素dkdl就“重要性”作二元比较,若dkdl重要,记定性标度ekl=1,elk=0;dkdl同样重要,记ekl=0.5,elk=0.5;若dldk重要,记ekl=0,elk=1。其中,k=1,2,…,6;l=1,2,…,6。以下矩阵为指标集D对重要性作二元比较的定性排序标度矩阵E

                              (3)

                              在二元比较过程中要求判断思维不出现矛盾,即要求逻辑判断的一致性,其一致性检验条件为:若ehk>ehl,有ekl=0;若ehk<ehl,有ekl=1;若ehk=ehl=0.5,有ekl=0.5。其中,h=1,2,…,6。

                              为了在二元定量对比中更便于语言习惯,大连理工大学的陈守煜教授给出了语气算子与重要性定量标度之间的对应关系:语气算子为同样、稍稍、略为、较为、明显、显著、十分、非常、及其、极端、无可比拟时,对应的相对隶属度分别为1.000、0.818、0.667、0.538、0.429、0.333、0.250、0.176、0.111、0.053、0。

                              (3)建立模糊评判矩阵。选取研究区内地质勘探钻孔的砂地比、砂岩厚度、粗砂岩厚度、砂岩层数、断层分维值和褶皱分维值作为含水层富水性评价的指标,在此基础上建立模糊评判矩阵R[10]

                              (4)

                              式中:rij——第i个因素对第j个评语的隶属度。

                              式中R通常称为模糊关系矩阵,由各单因素评判结果得到,隶属函数的建立是用来刻画模糊集合的,即用来计算rij。对模糊对象只有给出切合实际的隶属函数,才能应用模糊数学方法进行计算。三角形隶属度函数是最常见最简单的一种模糊隶属度函数,在此Ai均选取三角形隶属度函数。

                              (4)含水层富水性模糊综合评判。模糊综合评判数学模型的基本形式为:

                              (5)

                              式中:B——评判对象的综合评判结果;

                              A——n个参与模糊综合评判因素的权重;

                              R——模糊综合评判矩阵。

                              由于加权平均型算子对所有因素权重大小均衡兼顾,因此,在这里采用了加权平均算子。

                              2.1.4 现场应用

                              利用模糊综合评判结果,绘制红柳煤矿首采区直罗组下段含水层富水性分区,如图3所示。

                              图3 红柳煤矿首采区直罗组下段含水层富水性分区

                              根据对研究区内各工作面涌水量的观测台账,位于五级和四级富水性的工作面涌水量明显大于位于三级和二级富水区的工作面,基于地质勘探钻孔水文信息的含水层富水性分区结果与实际井下工作面涌水量较为一致。

                              2.2 工作面与矿井涌水量动态预测

                              2.2.1 工作面与矿井涌水量预测存在的问题

                              以往针对鸳鸯湖矿区工作面涌水量预测主要采用大井法或集水廊道法,其预测值通常比实际值小,通过分析工作面顶板充水含水层水文地质特征,认为工作面顶板直罗组下段砂岩含水层以孔隙为主,裂隙发育,地下水主要为古封存水,加之鸳鸯湖矿区构造发育,含水层平面上被断层切割,连通性较差,以静储量为主。采用大井法或集水廊道法计算的涌水量仅是工作面的动态补给量,而忽视了顶板含水层的静储量。

                              以往预测矿井涌水量时,通常将整个预测区域作为大井或集水廊道,所得到的涌水量预测值为固定值,不能实时指导各采区或矿井排水系统的建立。而矿井涌水量实际上是一个随着采掘空间和时间变化的动态值。

                              2.2.2 静储量结合动态补给量的工作面涌水量预测

                              基于以上分析,提出了静储量叠加动态补给量的“动静结合”工作面涌水量预测方法,工作面涌水量包含垮落带周边的侧向动态补给水量和垮落带内含水层的静储量,使工作面涌水量预测结果更加符合实际。

                              (1)静储量预测。静态储存量主要是指当煤层回采产生的导水裂缝带影响范围内的含水层由于重力作用而释放出来的水量,计算式如下:

                              Q=LB1

                              (6)

                              式中:Q——含水层的静储量,m3/h;

                              L——工作面走向长度,m;

                              B1——垮落区宽度,m;

                              M——含水层厚度,m;

                              μ——含水层的给水度。

                              (2)动态补给量预测。对于工作面动态补给量的预测可以采用“大井法”,其计算式如下:

                              Q

                              (7)

                              式中:Q——工作面动态补给量,m3/h;

                              K——含水层的渗透系数,m/d;

                              H——水头高度,m;

                              R0——引用影响半径,m;

                              r0——引用半径,m。

                              也可以采用集水廊道法,其计算式如下:

                              Q

                              (8)

                              式中:c——巷道底板周长,m;

                              R1——影响半径,m。

                              (3)工作面涌水量预测。工作面涌水量由顶板含水层的静储量和动态补给量组成,计算式如下:

                              Q=Q+Q

                              (9)

                              2.2.3 基于时空分区的矿井涌水量预测

                              为了使矿井涌水量能够实时指导防治水工作和排水系统的设置,提出了基于时空分区的矿井涌水量预测方法,即根据矿井采掘接续计划将预测期划分为若干时段,再分别预测各时段内各分区(工作面或巷道)的涌水量,最终得到一个随时间变化的矿井涌水量动态预测值[14],计算式如下:

                              (10)

                              i=1,2,3…mj=1,2,3…n

                              式中:Qi——第i时段矿井涌水量,m3/h;

                              qij——第i时段j处涌水量,m3/h。

                              基于时空分区的矿井涌水量也可以用于采区或不同水平涌水量预测。

                              2.2.4 现场应用

                              麦垛山煤矿11采区是2号煤层首采区,为了给防治水工作和排水系统建立提供依据,需要对11采区进行涌水量预测。根据11采区2021-2028年采掘接续方案,将预测期按照年度划分为8个时间段,分别对各时间段内不同工作面涌水量进行预测,最终得到11采区涌水量动态预测值,如图4所示。

                              图4 麦垛山煤矿11采区涌水量预测曲线

                              2.3 工作面顶板水害威胁程度定量化评价

                              2.3.1 工作面顶板含水层疏放水效果评价存在的问题

                              以往对工作面受底板灰岩水害威胁程度通常采用突水系数法,底板受构造破坏的地段突水系数一般不得大于0.06 MPa/m,隔水层完整无断裂构造破坏的地段不得大于0.10 MPa/m。由于突水系数简单易用,是评价工作面底板水害威胁程度的常用定量化方法。

                              受顶板砂岩水害威胁的工作面常采用顶板水疏放的防治水方法,但是由于缺少顶板水疏放效果定量化评价方法,一方面在未达到疏放水效果的情况下盲目开采会导致工作面涌水量较大,甚至发生水害事故,另一方面在达到疏放水效果时继续对顶板含水层进行疏放,不仅破坏了地下水资源,同时还使矿井涌水量增大,对矿区周边生态环境造成不良影响。

                              2.3.2 工作面顶板含水层疏放水效果定量化评价

                              基于受顶板水害威胁工作面的涌水量组成分析,影响和威胁工作面安全回采的主要为顶板含水层的静储量,这部分水量如果未经疏放,在工作面回采时会经过导水裂缝带进入采空区,导致涌水量增大。工作面回采时顶板含水层的动态补给量在短时间内是无法疏干的,配备满足要求的排水系统即可。

                              建立“双因素四指标”疏放水效果评价方法:工作面疏放总水量大于等于静态储存量,并且疏放水残余总水量小于等于动态补给量时,可以判断工作面顶板含水层的静态储存量得到了充分疏放,利用工作面疏放水观测值和涌水量预测值双因素中的静态储存量、动态补给量、疏放总水量和残余总水量4个指标进行比对,实现了对顶板水疏放效果的定量化评价[15],评价流程如图5所示。

                              图5 工作面顶板水疏放效果定量化评价流程

                              2.3.3 现场应用

                              利用“双因素四指标”法对鸳鸯湖矿区典型工作面疏放水效果进行了定量化评价,见表2。

                              表2 典型工作面顶板水疏放效果定量化评价

                              煤矿工作面工作面涌水量预测值静储量/104m3动态补给量/(m3·h-1)工作面疏放水观测值疏放总水量/104m3残余水量/(m3·h-1)疏放水效果评价梅花井23220445.1635.7660.8530.00好石槽村11220235.84196.0072.00183.00好石槽村11220623.21179.3346.79202.00较好红柳02020214.54161.40124.60143.80好麦垛山110207385.52344.85503.78296.40好

                              由表2可以看出 ,各工作面疏放总水量均大于静储量预测值,疏放水钻孔的残余水量与动态补给量预测值接近,说明工作面顶板含水层的静储量得到了有效疏放,钻孔残余水量为动态补给量,各工作面在回采过程中无集中涌水现象发生,均实现了安全回采。

                              3 结语

                              (1)通过建立水文地质钻孔与水文信息之间的相关关系,在此基础上利用地质勘探钻孔的水文信息对含水层进行富水性分区,有效地利用了井田范围内的地质勘探钻孔,为顶板水害防控提供了依据。

                              (2)考虑了静储量和动态补给量工作面涌水量预测与实际涌水量更为接近,而基于时空分区的矿井涌水量预测提供的动态预测值能够更好地指导矿井防治水工作和排水系统建立。

                              (3)“双因素四指标”评价方法实现了对顶板水疏放效果的定量化评价,避免了顶板水疏放不足造成水害事故,同时也能够显著减少顶板水的过度疏放。

                              参考文献:

                              [1] 董书宁,姬亚东,王皓,等.鄂尔多斯盆地侏罗纪煤田典型顶板水害防控技术与应用[J].煤炭学报,2020,45(7):2367-2375.

                              [2] 董书宁.煤矿安全高效生产地质保障的新技术新装备[J].中国煤炭,2020,46(9):15-23.

                              [3] 侯恩科,樊江伟,高利军,等.地面核磁共振技术在隐伏火烧区富水性探测中的应用[J].煤田地质与勘探,2021,49(5):230-237.

                              [4] 方刚,高波.榆横北区地面瞬变电磁法探查含水层特征应用[J].中国煤炭,2018,44(11):27-32,43.

                              [5] 武强,樊振丽,刘守强,等.基于GIS的信息融合型含水层富水性评价方法——富水性指数法[J].煤炭学报,2011,36(7):1124-1128.

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

                              [7] 梁向阳.基于综合探查手段的深埋矿井工作面顶板富水性分区[J].煤矿安全,2020,51(7):161-165.

                              [8] 郭小铭,王皓,周麟晟.煤层顶板巨厚基岩含水层空间富水性评价[J].煤炭科学技术,2021,49(9):167-175.

                              [9] 段东伟,丁湘,蒲治国,等.基于充水条件比拟的工作面涌水量预计方法[J].煤矿安全,2019,50(5):209-213.

                              [10] 尹尚先,张祥维,徐慧,等.“大井法”中渗透系数及含水层厚度的优化[J].煤田地质与勘探,2015,43(5):53-56.

                              [11] 李明泽.基于导水裂缝带的地下水动态模拟及预测[J].中国煤炭,2017,43(12):55-62.

                              [12] 李莉,李海霞,马兰.宁东煤田矿井水资源及其利用现状分析[J].干旱区资源与环境,2021,35(8):108-113.

                              [13] 赵宝峰.基于含水层沉积和构造特征的富水性分区[J].中国煤炭地质,2015,27(4):30-34.

                              [14] 赵宝峰,牛光亮,丁湘,等.基于时空分区的矿井涌水量预测方法研究[J].煤炭技术,2016,35(12):155-157.

                              [15] 赵宝峰.基于AHP-物元可拓模型的顶板水疏放程度评价[J].煤田地质与勘探,2017,45(3):82-86.

                              Prevention technology and its application of water hazard in coal seam roof in Yuanyanghu mining area of Ningdong coalfield

                              ZHAO Baofeng1,2, HUANG Zhongzheng3, ZONG Weiqin3

                              (1. Xi'an Research Institute of China Coal Technology & Engineering Group, Xi'an, Shaanxi 710054, China;2. Shaanxi Key Laboratory of Coal Mine Water Hazard Prevention and Control Technology, Xi'an, Shaanxi 710077, China;3. China Energy Ningxia Coal Industry Co., Ltd., Yinchuan, Ningxia 750011, China)

                              Abstract In order to avoid water hazard in roof in Yuanyanghu mining area of Ningdong coalfield, the water-abundance evaluation, water inflow prediction and water drainage effect evaluation of roof aquifer in shallow coal seam were systematically researched. Based on the hydrologic information of geological boreholes, a method of aquifer water-abundance zoning was proposed. By selecting six evaluation indexes affecting the water-abundance of the lower section of Zhiluo formation, the correlation between each evaluation index and the unit water inflow of boreholes was analyzed. The weight of each index was determined by mood operator comparison method. A fuzzy evaluation matrix was constructed by using triangle subordinating degree function, and the aquifer-abundance zones of the lower section of Zhiluo formaiton were divided by fuzzy comprehensive evaluation method. Considering the hydrogeological characteristics of roof sandstone aquifer, a prediction method of water inflow of working face based on static reserves and dynamic recharge was put forward; according to the composition and variation characteristics of mine water inflow, a dynamic prediction method of mine water inflow based on spatio-temporal zoning was proposed, based on the observed water volume of roof water drainage, combined with the predicted value of static reserves and dynamic recharge of the working face, a quantitative evaluation method of two indexes and four factors for the drainage effect of the working face was proposed. The results showed that the aquifer water-abundance zoning method, which was based on the hydrological information of geological boreholes, the water inflow prediction method of working face with static reserves and dynamic recharge, the mine water inflow prediction method based on spatio-temporal zoning and the quantitative evaluation method of roof water drainage effect, could effectively prevent and control the roof water hazard.

                              Key words roof water hazard; water-abundance zoning; water inflow prediction; water drainage effect evaluation; unit water inflow of borehole

                              移动扫码阅读

                              引用格式:赵宝峰,黄忠正,宗伟琴.宁东煤田鸳鸯湖矿区煤层顶板水害防控技术与应用[J].中国煤炭,2022,48(3)∶23-29. doi:10.19880/j.cnki.ccm.2022.03.004

                              ZHAO Baofeng, HUANG Zhongzheng, ZONG Weiqin.Prevention technology and its application of water hazard in coal seam roof in Yuanyanghu mining area of Ningdong coalfield[J].China Coal,2022,48(3):23-29. doi:10.19880/j.cnki.ccm.2022.03.004

                              基金项目:国家重点研发计划(2017YFC0804102),国家自然科学基金(41807221),陕西省自然科学基础研究计划(2018JQ5150),宁夏人才小高地项目

                              作者简介:赵宝峰(1981-),男,河北涉县人,博士,研究员,主要从事矿井防治水工作。E-mail:sunman1220@163.com

                              中图分类号 TD742.1

                              文献标志码 A

                              (责任编辑 张艳华)

                              九州彩票