采动覆岩卸压瓦斯运储通道属性参数拓扑关系研究

2024-04-08赵鹏翔卓日升李树刚林海飞常泽晨贾永勇刘元嘉

煤炭科学技术 2024年2期

赵鹏翔，卓日升，李树刚,2,3 ，林海飞,2,3 ，常泽晨，贾永勇，金权，刘元嘉

（1.西安科技大学安全科学与工程学院, 陕西西安 710054；2.教育部西部矿井开采及灾害防治重点实验室, 陕西西安 710054；3.西部矿井瓦斯智能抽采工程研究中心, 陕西西安 710054；4.新疆工程学院新疆煤矿灾害智能防控与应急重点实验室，新疆乌鲁木齐 830023；5.新疆工程学院新疆煤炭资源绿色开采教育部重点实验室, 新疆乌鲁木齐 830023；6.新疆维吾尔自治区煤炭科学研究所，新疆乌鲁木齐 830091；7.新疆工程学院，新疆乌鲁木齐 830023）

0 引言

采动卸压瓦斯运储的主要方式表现为瓦斯在煤岩体裂隙中进行解吸—渗流—升浮—扩散—储集[1-2]。煤岩体渗透性受内部微细观裂隙影响较大[3-4]，微细观裂隙在开采扰动作用下逐渐扩展为宏观裂隙，同时为瓦斯在煤岩体中运储提供了通道与场所[5-6]。

早在19 世纪初期，国内外学者针对煤岩体裂隙扩展发育特征做了大量研究，同时通过开展煤岩体力学测试试验研究了裂隙的微观结构特征。首先，周世宁院士探讨了微细观裂隙的产生及发育机理，得到裂隙系统是构成瓦斯流动的主要通道[7]。TIEN[8]、JHA 等[9]、YASITLI 等[10-11]、XIE 等[12]、SMYTH 等[13]采用数值模拟方法，分析了煤岩体变形后应力变化规律及内部裂隙发育特征。国内学者齐庆新等[14-15]采用测窗法和钻孔窥视相结合的方法，得到了覆岩受采动影响下的内部裂隙分布规律及演化特征。受采动影响，覆岩裂隙发生扩展、连通构成了卸压瓦斯运储的通道，使得煤岩层的渗透性得到改善[16]。近年来，国内学者对于卸压瓦斯运储通道的定义及识别方法的研究也取得一定的成果，姜振学等[17]运用物理实验和现场测量相结合的方法得到不同地质条件下油气的运移方向，并定义了级差、分隔、流向、流压和断面优势综合作用下的油气优势运移通道。牛彦良等[18]、魏建平等[19]、李宏义等[20]、JIA 等[21]对于油气在岩体中运移特征进行了系统性总结，并开展了不同影响因素条件下油气运移特性的研究。张勇等[22-23]将采动影响下的煤岩体中瓦斯流通通道划分为宏观通道和细观通道，同时结合岩体断裂力学，对所划分的区域开展进一步的机理研究[24-25]，发现瓦斯运移行为主要发生在宏观通道中。李立[26]从空间尺度上出发，将瓦斯通道划分为3 个区域，自上而下分别是瓦斯紊流区、瓦斯过渡流通道区以及瓦斯渗流通道区，并通过UDEC 数值模拟方法探究了瓦斯在各区域的流通特性。肖康等[27]在微观尺度上提出用微粒运移理论来解释卸压瓦斯运储通道的形成机理，并结合压力损耗规律建立动态孔隙网络模型。基于李树刚教授提出的采动裂隙椭抛带[28]，刘洪永等[29]采用数值模拟方法，开展不同推进速度条件的卸压瓦斯运储通道的演化规律实验，发现覆岩中的垂向裂隙为瓦斯运移提供主要通道。笔者[30-31]运用物理模拟和数值模拟相结合的方法探究了开采高度、推进速度等不同因素对采空区空间形态的影响，深入分析了多因素影响条件下覆岩裂隙网络分布特性，并采用分形理论将覆岩划分为瓦斯活跃区、瓦斯流通区以及瓦斯富集区[32]。

近年来，国内外诸多学者从多个角度，针对采动覆岩裂隙特征进行了较为有效的表征[33-34]，系统地阐述了覆岩裂隙网络通常是由主干裂隙延伸、扩展、联通次级裂隙构成，不同地层岩性和构造条件下的裂隙发育密度、隙宽差异很大，总结出通常煤岩体裂隙网络具有非均质性、连续性、各向异性和定向性。以上研究从不同角度，针对不同裂隙网络进行了较为有效地表征，但都未考虑对裂隙网络性能具有重要影响的开采条件及网络联接状态。覆岩裂隙网络的联接状态多用具有拓扑特征来表示[35-36]，裂隙网络可以看作是一系列相互连接的裂隙的集合，用节点和线来表示，每个节点表示一个裂隙，线表示裂隙之间的连接关系。裂隙的属性包括长度、宽度、厚度等。岩石裂隙网络的连接拓扑结构具有局部聚类的特征，在不同开采条件影响下，裂隙网络的外部几何参数及内部裂隙属性参数（破断裂隙密度、离层量、贯通度以及分形维数等）会伴生对卸压瓦斯运储通道影响程度的演变，进一步促进卸压瓦斯运储通道的形成。

因此，笔者采用二维物理相似模拟试验，结合拓扑学方法，从卸压瓦斯运储通道属性参数拓扑关系角度，分析各因素对优势通道属性参数作用，建立卸压瓦斯运储通道网络干系结构模型，表征不同开采条件下卸压瓦斯运储通道的拓扑演化结构，分析其几何特征和裂隙网络发育参数，实现了覆岩裂隙网络干系影响因素的量化描述，为研究覆岩裂隙网络形成、扩展机制提供新的思路，并为优化卸压瓦斯抽采系统布置提供了理论支撑，助力实现卸压瓦斯的精准抽采。

1 试验设计

1.1 试验工作面概况

以山西和顺天池能源有限责任公司15 号煤层302 工作面作为试验原型。通过查阅前期地质勘探相关资料，发现该矿井地质构造简单，地层总体呈走向北东，倾向北西的单斜构造，地层倾角平缓，一般在15°左右，褶皱构造发育程度较低，多为宽缓褶曲，工作面内仅有1 条不甚发育的断层，不影响工作面的正常布置。

根据矿井岩层钻探数据分析，15 号煤层距离K2灰岩（关键层）平均约20 m，煤层直接顶为泥岩，直接底为铝质泥岩，基本底为泥岩，煤岩物理力学特性如图1 所示。该工作面平均煤层厚度为5.1 m，平均煤层倾角为8°，属于近水平煤层，工作面日推进度平均为3 m/d。工作面相对瓦斯涌出量7.8 m3/t，瓦斯压力0.45 MPa，煤层透气性系数为0.87 m2/(MPa2·d)，瓦斯含量实测值最大16 m3/t，属于高瓦斯煤层。

图1 工作面煤岩层物理力学性质Fig.1 Physical and mechanical properties of prototype coal formations

1.2 试验设计

为探究不同开采条件下卸压瓦斯运储通道属性参数拓扑关系，采用自主研发智能多角度平面物理相似模拟综合实验台，该实验台主要由动态底座系统、高强度框架搭建系统、智能变角度系统、数据监测采集系统等组成（图2）。

图2 智能多角度平面物理相似模拟综合实验台Fig.2 Intelligent multi angle plane physics similarity simulation comprehensive experimental platform

根据已有的相似准则及302 主采工作面的实际情况，结合拓扑学理论进一步探究卸压瓦斯运储通道几何、裂隙参数的拓扑关系特征。该工作面的平均推进速度为5 m/d，煤层平均厚度为6.15 m，平均采高为4 m。由于煤层赋存的非均质性、开采条件的复杂性，为更好地模拟现场开采情况，针对采高、推速、K2层位和倾角4 个因素，分别结合矿上实际条件每个因素设计3 个梯度，试验设计见表1。

表1 试验设计方案Table 1 Experimental design scheme

依据302 工作面实际开采工艺，结合矿井地质因素设计试验，试验流程如下：

1）模型搭建。以石膏和淀粉为凝胶剂，河砂为骨架。参照相应的配比号，模拟不同条件下岩体特性。每层模拟厚度保证为1 cm，并且分层材料选用效果较好的云母片，模拟岩层的层理结构，在确保材料相似性的同时也提高了试验测量的准确性。自然风干15～20 d，使得模型充分晾干，并根据工作面埋深对覆岩顶部均匀布置载荷。

2）开采设计。为了去除边界效应的影响，两端分别预留10 cm 的煤柱，开切眼距离为8 cm。为了避免来压步距的奇偶性，采取2、3 cm 交替开采的设计，推进速度按照时间相似比（1:10）进行计算。

3）数据采集。在试验过程中，分别对底板应力变化及各下沉量的测点位移进行记录，同时测量覆岩采动裂隙起裂位置、长度及张开度，并采用高清摄像机Alpha 9 记录采动覆岩裂隙的发育过程。为深入分析不同开采影响条件下覆岩裂隙网络的演化规律提供基础

1.3 试验结果

通过探究不同影响条件下（采高、推进速度、煤层倾角、关键层位）覆岩裂隙网络分布形态，分析不同影响因素下瓦斯运储通道分布区域位置（图3），得到上覆岩层采动裂隙“纵两带”的演化高度及形态。

图3 多因素条件下采动覆岩裂隙“纵两带”分布特征Fig.3 Distribution characteristics of the “two zones” of the overburden strata

采动覆岩“纵两带”在采高、推速、K2层位以及倾角的影响下发育范围不同，其中裂隙带发育高度受采高影响较大，6 m 采高条件下裂隙带的发育高度是2 m 采高的1.39 倍，同时裂隙带的变化趋势明显大于冒落带，如图3a 所示。

加快推进速度，裂隙带发育高度逐渐减小。工作面快速的推进导致裂隙发育时间较短，裂隙发育不充分，如图3b 所示。在3 m/d 的推进速度下，裂隙带高度发育距离煤层底板59 m，分别是7、5 m/d 的1.51、1.15 倍，且裂隙带所占的区域比例由60%减小至38%。

K2位于不同层位时，裂隙的发育及岩层移动并不完全同步，随着工作面的推进，覆岩离层及岩层垮落向上发展趋势依赖于关键层层位分布。随着K2距煤层顶板距离的增加，裂隙带及冒落带发育高度也随之增大，如图3c 所示。

不同的煤层倾角主要影响采动裂隙“纵两带”的发育形态，在倾斜煤层回采过程中由于覆岩不是垂直向下垮落，因此覆岩的垮落形态呈现不对称变化，随着煤层倾角的增大，这种不对称形态变化的愈加明显，如图3d，裂隙带的发育范围也随着煤层倾角的增大而呈现增大的变化趋势。

2 卸压瓦斯运储通道裂隙网络全域拓扑构架

煤岩体在开采扰动影响下，原始应力重新分布，导致煤岩体发生弯曲、破裂、垮落、压实等力学行为，煤岩体经历了形变、损伤、破裂的全过程[37]，其中裂隙的发育呈现产生、发育、扩展、贯通、张开及压实的空间演化过程，通过拓扑学的思想将各阶段抽象为节点与线段的连接方式，便形成了煤岩体裂隙拓扑结构。煤、岩中的裂缝网络发育是复杂、无序的，其拓扑结构亦多样变化，不同的拓扑结构会使裂隙网络系统的连续性、连通性等性能产生较大差异，进而对瓦斯运储产生较大影响。要想比较充分地表征岩石裂缝网络，既要考虑裂隙的几何、裂隙参数，还要掌握裂隙网络的拓扑结构。实现采动覆岩卸压瓦斯运移通道的简化表征，可充分描述煤岩体内部裂纹、裂隙的位置、形状及发育情况，能够为覆岩裂隙场中的卸压瓦斯运储通道的孪生映射提供基础，大幅提升瓦斯抽采系统布置的透明度。需要注意的是，实际煤岩体的拓扑网络结构非常复杂，通常需要通过煤岩体性质的测试、煤岩心观察及覆岩钻孔窥视等手段来获得更详细的拓扑结构信息[38-39]。本文通过开展不同因素下二维物理模拟实验，获得采动覆岩裂隙的拓扑结构分布特征。将煤岩体拓扑结构细分为原生裂纹结构、次生裂隙结构、重新压实结构、离层裂隙结构、破断裂隙结构和网状结构[40-41]：

1）原生裂纹结构：主要为天然裂纹，原生裂纹结构对瓦斯的渗流有重要影响。

2）次生裂隙结构：裂纹机构遭受拉伸或伸缩应力时，裂纹会以煤岩体应力集中点或弱点开始发育、扩展，形成次生裂隙结构。

3）重新压实结构：当垮落岩体在扰动应力的作用下重新压实，次生裂隙网络结构逐渐闭合，煤岩体形成紧密堆积的结构，这种结构增加了垮落岩体的密度和强度，瓦斯渗流作用减弱。

4）离层裂隙结构：由于上下相邻岩层具有不同的弯曲刚度，在采动应力作用下，相邻岩层出现变形、弯曲、破断及下沉，岩层间出现离层裂隙，这种结构为瓦斯储集提供场所。

5）破断裂隙结构：岩层在张拉破坏或受压屈服后产生破断、回转，破断岩块与相邻岩块产生纵向破断裂隙，这种结构对瓦斯运移起着重要作用。

6）网状结构：煤岩体在开采扰动作用下，离层裂隙与破断裂隙会形成交叉连接的网状结构，这种结构是瓦斯运储通道的重要组成结构。

文中的拓扑结构即指卸压瓦斯运储通道发育网络M（Mesh）通过映射所表现的抽象结构M′。这种简单的拓扑结构的拓扑性质可以类比卸压瓦斯运储通道网络发育的拓扑性质，影响卸压瓦斯运储通道拓扑结构的主要因素包括开采条件、岩层性质等。随着回采作业的进行，裂隙网络拓扑结构也逐渐发育，煤岩体内的原生裂纹结构发育态势较缓，在采动应力作用下扩展而成的次生裂隙结构逐渐转化为离层裂隙结构和破断裂隙结构，采空区内部垮落岩体的裂隙结构在扰动应力作用下形成重新压实结构，且分布范围随采动作业的进行而不断扩大。同时，离层裂隙和破断裂隙结构逐渐相互贯通成为网状结构，卸压瓦斯在网状结构内发生大量的运移、储集行为，是卸压瓦斯运储的主要场所，如图4 所示。基于复杂网络演化理论可以简化现实干系的拓扑变换[42]，通过对拓扑网络的可视化探明卸压瓦斯运储网络的演化因素，同时结合复杂网络特征参数实现运储通道干系结构和因素节点的演化特征分析。

图4 卸压瓦斯运储通道网络拓扑结构Fig.4 Gas storage and migration channel network topology

3 卸压瓦斯运储通道属性参数拓扑结构特征

为进一步明晰卸压瓦斯运储通道发育特征在不同影响因素下的演化规律，选取运储通道的高度和宽度作为表征其发育范围趋势的几何参数，选取破断裂隙密度、离层量、贯通度以及分形维数的变化特征来定量描述瓦斯运储能力的大小，见表2。统计、分析不同影响因素下运储通道几何参数和裂隙参数的变化趋势，为进一步分析运储通道干系结构和因素节点的演化特征提供了基础试验数据。

表2 不同影响因素下卸压瓦斯运储通道属性参数的选取Table 2 Selection of attribute parameters of channels under different influencing factors

3.1 卸压瓦斯运储通道的几何参数与各因素的拓扑关系特征

3.1.1 卸压瓦斯运储通道高度演化的多因素拓扑架构

卸压瓦斯运储通道的高度直接决定瓦斯抽采钻孔终孔位置的确定，对比不同因素影响下卸压瓦斯运储通道几何参数的变化规律，发现随着采高、K2层位及煤层倾角的增加，卸压瓦斯运储通道的高度发育范围也不断增加，如图5 所示。

图5 运储通道高度与影响因素的关系拓扑Fig.5 Topology of relationship between channel height and influencing factors

采出空间大，岩层拥有充分的自由垮落空间，卸压瓦斯运储通道不断向上发育，直到发育到K2关键层的位置，表现为采高6 m 时的瓦斯运储通道高度是2 m 的1.17 倍，在K2距煤层顶板30 m 条件下，通道高度仅发育至45 m。

工作面快速推进，覆岩裂隙带内过早的出现悬臂梁[43]，从而更早的形成岩层压力拱，因此裂隙网络发育不充分，发育范围增长减缓，3 m/d 推速下的卸压瓦斯运储通道发育高度是7 m/d 的1.73 倍。

通过计算绘制不同因素条件下瓦斯运储通道发育高度雷达图，如图6 所示。发现采高和推速对高度发育影响最大，高度分别发育至距煤层顶板59.5、57.5 m。而采高2 m 和推速7 m/d 时的发育高度仅为30.2 m 和37.9 m，高度变化率分别为1.97%和1.52%，而不同倾角和K2层位下变化率仅为1.15%、1.21%。采高、倾角、K2层位与瓦斯运储通道高度变化规律成正比，和推速呈反比。

图6 瓦斯运储通道发育高度因素雷达图Fig.6 Radar map of channel development height factor

结合不同因素对与运储通道高度的影响规律，用通道高度和影响因素作为节点，节点间的距离表示之间的关系，其中线段长度表示关系的紧密程度，实线表示为正相关变化趋势，虚线表示负相关的变化趋势，如图4 所示。采高对通道的高度影响最大，K2层位对其影响次之，推速和倾角的变化对于通道高度影响最小，其中随着推速加快，卸压瓦斯运储通道发育高度呈现降低的趋势。

3.1.2 卸压瓦斯运储通道宽度演化的多因素拓扑架构

抽采钻孔的偏角主要依据卸压瓦斯运储通道宽度的变化特征来进行设计，不同因素下瓦斯运储通道的宽度变化趋势不同。对比不同因素下卸压瓦斯运储通道宽度变化规律，发现回采前期卸压瓦斯运储通道宽度变化率都较小，快速的发育主要集中在2～3 次周期来压后，如图7 所示。

图7 通道宽度与开采条件关系拓扑图Fig.7 Topology of channel width versus mining conditions

通过绘制瓦斯运储通道宽度变化因素雷达图，如图8 所示。发现K2关键层距离煤层越远，覆岩裂隙向上发育受关键层影响越晚，通道内的裂隙有着充分的时间进行横向发育，最宽达到35.4 m。分析瓦斯运储通道在采空区两侧的几何参数的差异，发现倾角30°时，内外边界最大可达12 m，这一差异为卸压瓦斯提供了优质的储集场所，且煤层倾角越大，这种差异越明显。

图8 通道宽度变化因素雷达图Fig.8 Radar map of variation factors of channel width

通过对比不同因素下卸压瓦斯运储通道的宽度变化趋势，得到卸压瓦斯运储通道宽度的影响因素影响拓扑关系（图7），图中“～”为连接折点，用来表示关系微弱的两个节点。可以看出，通道主要受倾角和K2层位来共同控制，受推速和采高的影响较小。

3.2 卸压瓦斯运储通道的裂隙参数与各因素的拓扑关系特征

3.2.1 卸压瓦斯运储通道破断裂隙密度变化规律及拓扑构架

为了定量描述破断裂隙发育的程度，以裂隙密度 ρf（条/m）表示破断裂隙随工作面推进过程中的变化特征[44]。依据试验统计数据，在模型走向方向，按10 m 为1 个计数单元统计该范围内裂隙条数。通过计算绘制不同影响因素下的破断裂隙的密度，发现破断裂隙密度均发生4 次突变（图9）。

图9 充分采动时覆岩裂隙破断裂隙密度分布Fig.9 Distribution of overburden fracture breaking fracture density

随着距煤柱距离增加，切眼侧与煤壁侧裂隙密度不断增加，采动裂隙随采空区面积的增加逐步向高位发育。工作面从35 m 推进至70 m 时，此区域内的裂隙密度普遍较低，说明在采空区中部受上覆载荷的作用，采动裂隙发生压实、闭合，破断裂隙密度骤降。可以看出，裂隙富集区域主要集中在煤壁前后20 m 范围，煤壁侧的裂隙数量大于切眼侧，裂隙密度变化规律呈类“马鞍”状。

结合图9 破断裂隙分布特征，得到卸压瓦斯运储通道内破断裂隙密度的影响因素影响拓扑构架（图10），发现卸压瓦斯运储通道内破断裂隙密度的变化主要受煤层倾角影响，受推速的影响较小，而不同采高下运储通道的破断裂隙密度的变化不显著。破断裂隙密度高的区域瓦斯更易进行运移，根据多因素条件下破断裂隙密度的变化趋势，进一步优化瓦斯抽采钻孔终孔位置，使得终孔位置始终处于高密度区域（破断裂隙密度>4.7 条/m），有效地起到瓦斯截留的作用。

图10 破断裂隙密度影响因素拓扑干系图Fig.10 Topological stem diagram of factors influencing fracture density

3.2.2 卸压瓦斯运储通道离层裂隙的变化规律拓扑特征

以离层量（S）[31]（相邻岩层间的离层高度，m）来定量描述采动过程中卸压瓦斯运储通道中离层裂隙的变化特征。根据工作面各岩层间的下沉量，可得各岩层间的离层量分布，不同开采条件下上覆岩层离层量随距煤柱距离的变化规律如图11 所示。

图11 充分采动时多因素下离层量的变化曲线Fig.11 Variation curve of overburden detachment volume under multiple factors

通过计算充分采动时多因素下瓦斯运储通道内离层量的均值（S），发现随采高、倾角以及关键层层位的增大，S逐渐升高，如图10 所示。覆岩在开采扰动下发生岩层间出现离层现象，在岩体压实的区域覆岩离层量较小，而在岩体交接、错落的地方离层量较高，岩层间裂隙张开度也较大（0.5 m≤S≤2.0 m），为气体的储集提供储集场所。对比分析可得，采高的增加对S的影响最大，6 m 采高时相较于4 m 采高时增幅达到146.9%，K2层位的变化对S影响较小，增幅分别为3.1%，12.3%。同时，随着推速加快，覆岩发育时间缩小，S呈现降低的变化趋势，7 m/d 推速下S仅为1.34 m。

图12 为卸压瓦斯运储通道内离层量的影响因素影响拓扑干系图，发现离层量主要受采高和煤层倾角影响，推速和K2层位下离层量变化不明显，如图12 所示。高离层量的区域靠近煤壁及开切眼侧，此处收应力影响显著，瓦斯运储通道结构变化频繁，更易实现瓦斯抽采孔（巷）“延伸–破断”自主调节的目的，进而避免孔（巷）伸入采空区过长或过短对瓦斯抽采效果所造成的影响。

图12 离层量影响因素拓扑干系图Fig.12 Topological stem diagram of factors influencing Variation curve

3.2.3 卸压瓦斯运储通道中覆岩裂隙贯通度变化规律及拓扑构架

煤层回采过程中，上覆岩层破断后在采空区弯曲下沉，影响卸压瓦斯运储效果。因此引入破断裂隙贯通度（Di）[45]来描述岩层的破断程度及破断裂隙内卸压瓦斯运储能力的强弱，岩层的裂隙贯通度主要受下位岩层碎胀的填充程度、岩层厚度、岩体力学性质以及岩层的周期来压距共同影响。不同采高、推速、煤层倾角以及关键层层位条件下，下位岩层碎胀的填充程度和周期破断距不同，导致不同开采条件下卸压瓦斯运储通道的贯通度呈不同的变化规律，如图13 所示。

图13 多因素条件下运储通道贯通度变化规律Fig.13 Variation pattern of storage and migration channel penetration under different mining conditions

采高直接影响着下位岩层碎胀的填充程度，采高越大，垮落的自由空间越大，岩层的回转角和弯曲曲率越大，岩层破断裂隙的贯通度越大。采高6 m时，距离煤层顶板35 m 处，贯通度出现下降的拐点，而采高2 m 在距离煤层顶板15 m 时贯通度出现下降拐点。K2层位越高，关键层的下位岩层拥有充足的填充空间，自由空间阈值较高，贯通度也随之增大。推速越快，岩体未得到充分的时间发育，卸压瓦斯运储通道发育高度较低，7 m/d 的覆岩自由空间阈值远小于3 m/d，可以看出贯通度与采高及K2层位呈正相关变化趋势，而与推速呈负相关变化趋势。

岩层破断的会转角受岩层周期来压距和煤层倾角共同影响，随着倾角增大，周期破断距逐渐减小[30]，岩层的破断角度逐渐增大，贯通度逐渐增大，水平煤层的同一高度卸压瓦斯运储通道贯通度仅为倾斜煤层时的60%。

图14 为卸压瓦斯运储通道内贯通度的影响因素影响拓扑干系图，发现贯通度受采高影响最大，煤层倾角、K2层位次之，推速对于贯通度影响最小。结合图13 贯通度变化规律可知，卸压瓦斯运储通道内低贯通度的位置瓦斯难以大量储集。因此，贯通度大的区域岩层破断程度高，同时此处采动裂隙导气能力强，瓦斯运移更为活跃，将孔（巷）的终孔位置动态调整布置在高贯通度区域，有效地提高“截留”瓦斯的浓度，减少了瓦斯抽采效率的波动性。

图14 贯通度影响因素拓扑干系图Fig.14 Topological stem diagram of factors influencing penetration

3.2.4 卸压瓦斯运储通道分形维数变化规律及拓扑构架

受开采扰动，卸压瓦斯运储通道在自重力和水平应力的作用下随覆岩的垮落、滑移、下沉及断裂而逐渐发育，卸压瓦斯运储通道内的裂隙网络逐步向采空区上方扩展。在周期来压后，岩体出现大面积的弯曲下沉及破断，裂隙网络充分发育，分形维数急剧增大，计算得到不同因素下卸压瓦斯运储通道的分形维数变化规律[46]，如图15 所示。

图15 运储通道的分形维数变化曲线Fig.15 Variation curve of fractal dimension of channels

结果表明，随着采高的增加，运储通道分维值有增大的趋势。裂隙网络发育复杂程度为：6.0 m>4.0 m>2.0 m。K2层位越高，卸压瓦斯运储通道有着较多的时间去缓慢发育，分形维数也随着增大，与采高影响下的分形维数变化规律一致。

随着开采的不断前进，裂隙网络在来压后发育充分，也逐步完善，因此推进速度减缓，裂隙的发育速度逐渐却加快，此时推速3 m/d 的分形维数逐渐超过7 m/d，并且增长趋势最大。可以看出推速越缓慢，卸压瓦斯运储通道发育越完善。

通过分析运储通道的分形维数变化趋势以及与影响因素拓扑特征，如图16 所示。发现随着采高增加，卸压瓦斯运储通道的分形维数逐渐升高，关键层的位置对于卸压瓦斯运储通道分形维数的影响最小，表明在大采高的条件下，裂隙网络发育的更为复杂。分形维数的变化直观表征着裂隙网络的演化趋势，高维区域的裂隙横纵交错复杂，为瓦斯的运储提供良好的场所，明确不同因素下高低维分布区域，可增大孔（巷）的有效抽采距离，避免将孔（巷）在抽采作业中出现坍塌等现象，提高瓦斯的抽采效率。

图16 分形维数与开采条件的关系拓扑图Fig.16 Topological stem diagram of factors influencing.

4 多因素下卸压瓦斯运储通道网络全域影响干系拓扑因素模型

结合拓扑学方法，从拓扑演化的角度出发[47-48]，分析多因素条件下卸压瓦斯运储通道几何参数、裂隙参数的变化规律及拓扑干系结构，实现了覆岩裂隙网络干系影响因素的定量描述，建立了卸压瓦斯运储通道网络干系结构矩阵及模型，如图17 所示。煤层开采条件–卸压瓦斯运储通道参数拓扑关系存在多个层次，分别为：①单开采条件下通道不同特征呈现不同的变化关系；②单个通道参数与不同开采条件下的影响规律不同；③分属不同开采条件下，其开采因素的组合也不同；④通道内的不同参数（如：几何参数、裂隙参数等）存在深度并发关系、串联关系、补偿互馈关系等。

图17 卸压瓦斯运储通道干系网络矩阵Fig.17 Gas storage and migration channel stem network matrix

层次①②中的拓扑关系可以用数字的大小表示干系紧密程度，数字“1”表示两者相互关系紧密级别第一，“2”表示两者相互关系紧密级别第二，数字的大小直接表示了关系间的紧密程度，其中“－∞”表示开采条件对通道=参数的影响微乎其微。不同开采条件下瓦斯运储通道属性参数间的拓扑关系，可以结合第三章的拓扑图中节点之间的距离进行计算，在同一条件下，数字的序号从小到大，则表示了此参数与开采条件影响关系逐渐减弱，当节点间拓扑关系出现“～”，表示相互关系极其微弱，在干系网络矩阵中用“－∞”表示。

①②层次处于交互状态，为了明确不同开采条件与卸压瓦斯运储通道的参数特征之间的拓扑关系，需要对干系网络进行统一编码，得到卸压瓦斯运储通道发育干系网络矩阵，如图17 所示。通过对干系网络矩阵进行遍历，即可得到卸压瓦斯运储通道不同开采条件下的连接拓扑关系。干系网络关系的统一编码确定了卸压瓦斯运储通道特征参数演化行为因素影响规律，瓦斯抽采系统的布置根据通道参数与开采条件的拓扑关系以及其他规则，进行动态优化调整，从而提高瓦斯抽采效率。为了更为直观地表示二者之间的拓扑关系，同时便于实际工程上的应用，将卸压瓦斯运储通道的干系网络矩阵转化为卸压瓦斯运储通道干系因素结构图，如图18 所示。卸压瓦斯运储通道的裂隙发育网络可抽象成有节点和线组成的拓扑结构，以表示运储通道网络发育干系因素影响规律。其中，节点表示主控因素，线可以象征影响因素与运储通道的干系结构。干系紧密程度统一编码确定了开采条件对瓦斯运储通道参数的影响次序，瓦斯抽采工程根据其与开采条件的拓扑关系以及其他规则，进行布置层位、施工工艺的逐时段调节计算，从而优化瓦斯抽采工程对应开采条件下的抽采效果。

图18 卸压瓦斯运储通道干系因素拓扑模型Fig.18 Structure of gas storage and migration channel stem factors

根据上述方法，确定卸压瓦斯运储通道不同开采条件下的连接拓扑关系，由此构建卸压瓦斯运储通道发育干系网络模型。可以看出卸压瓦斯运储通道的形态由采高、推速、K2层位以及煤层倾角共同决定。

结合上文中试验所分析的数据，发现采高对于高度、离层量以及贯通度影响较大，推速主要影响运储通道的高度、破断裂隙密度，关键层层位主要控制卸压瓦斯运储通道的高度、宽度及贯通度，同时和煤层倾角共同影响运储通道的宽度发育范围，而不同搭配的开采条件综合影响着卸压瓦斯运储通道分形维数的变化趋势，控制着运储通道网络发育的复杂程度，依此为不同开采条件下卸压瓦斯抽采钻孔的布置及优化提供了一定的理论指导作用。