采空区充水过程水位变化特征及影响因素研究.pdf

1804factors during water filling in goafJ.Safety in Coal Mines,2023,54(9):180-186.移动扫码阅读GAO Xin,DENG Cunbao,XING Yuzhong.Research on water level variation characteristics and influencing54(9):180-186.高鑫，邓存水过程水位变化特征及系研分全.2023SafetyinCoal MinesSep.20232023年9 月煤防发全No.9Vol.54第9 期第5 4 卷DOI:10.13347/ki.mkaq.2023.09.024采空区充水过程水位变化特征及影响因素研究高鑫1 2，邓存宝，邢玉忠1.21.太原理工大学安全与应急管理工程学院，山西太原0 3 0 0 2 4 2.山西省煤矿安全研究生教育创新中心，山西太原0 3 0 0 2摘要：影响矿井突水水位的重要因素有突水强度、矿井巷道和矿井采空区；为研究矿井突水后采空区充水过程中的内部水位变化特征及影响因素，搭建了实验平台，进行采空区充水实验研究。结果表明：在研石对水流的阻滞作用下，采空区充水过程分为初期水流蔓延阶段和后续充水阶段，整个过程中近工作面处水位较高而深部水位较低，水面出现坡度且随时间逐渐减小；采空区充水率随时间增长且近似二次函数变化，突水量和矿井空间结构是影响充水过程的重要因素。关键词：矿井突水；水位；相似模拟；变化特征；影响因素中图分类号：TD745文献标志码：A文章编号：1 0 0 3-4 9 6 X（2 0 2 3）0 9-0 1 8 0-0 7Research on water level variation characteristics and influencing factors during water filling in goafGAO Xin2,DENG Cunbao,XING Yuzhong1,2(1.College of Safety and Emergency Management Engineering,Taiyuan University of Technology,Taiyuan 030024,China;2.Center of Shanxi Mine Safety for Graduate Education Innovation,Taiyuan 030024,China)Abstract:The important factors that affect the water inrush level in mines are the intensity of water inrush,mine roadways and minegoafs.In order to study the characteristics and influencing factors of internal water level changes in the process of water filling ingoafs after mine water inrush,an experimental platform was built to conduct experimental research on water filling in goafs.The res-ults show that:under the blocking effect of gangue on water flow,the water filling process of the goaf is divided into the initial wa-ter flow spreading stage and the subsequent water filling stage.In the whole process,the water level near the working face is higherwhile the deep water level is lower,and the water surface slope appears and gradually decreases with time;the water filling rate ofgoaf increases with time and approximates a quadratic function change.Water inrush and mine spatial structure are important factorsaffecting the water filling process.Key words:mine water inrush;water level;similar simulation;variation features;influencing factors随着煤矿开采深度的增加，井下水文地质条件日益复杂，突水事故危害日益严重，一旦发生突水事故，往往会造成巨额财产损失和大量人员伤亡；突水后避难与救援等相关工作的开展很大程度上取决于矿井水位发展情况，有效的矿井突水水位预测及水灾路径搜索手段具有重要意义1-3 。目前，以基于图论及网络理论的研究为主；马恒等 4 以巷道交叉点的标高为搜索权重提出了水流下向、上向蔓延路径搜索算法，并应用于未来时间段的突水范围和水位标高点预测；蔡明杰等 5 采用无向图和邻接表对矿井巷道网络进行描述和存储，使用优化SPFA算法进行了单源路线搜索；于丹等 6 结合随时间变化的水位高度对逃生路径选择的影响，建立了以巷道初始当量长度与危险性程度为衡量最优路径的权值时变数学模型，并结合改进的Dijkstra算法进行求解，优化后的模型收稿日期：2 0 2 2-0 9-0 9责任编辑：庚馨杰基金项目：国家自然科学基金-山西煤基低碳联合基金重点支持资助项目（U1710258)作者简介：高鑫（1 9 9 8），男，山西大同人，硕士研究生，研究方向为矿井灾害防治。E-mail：7 9 3 8 7 7 5 1 6 q q.c o m181SafetyinCoal MinesSep.20232023年9 月煤砺发全No.9Vo1.54第9 期第5 4 卷更能反映突水时巷道内的实际情景。突水强度、矿井巷道和矿井采空区是影响矿井突水水位发展的重要因素，采空区具有不同于巷道空间的孔隙结构，突水涌入采空区后的流动过程将不同于其在巷道内的流动情况，且将对矿井水位的上升过程起到重要影响 7-8 。综上，为研究突水后采空区充水过程中的内部水位变化特征及其影响因素，以高河能源E1305工作面为工程背景，以单个采空区为切人点设计了采空区充水实验；研究成果对矿井突水水位变化发展研究具有一定的借鉴意义。1采空区充水实验研究1.1采空区两带发育及储水结构特征采空区两带结构示意图如图1。曲下沉带裂隙带0000000000000000跨落000000000000000000500000图1 买采空区两带结构示意图Fig.1Schematic diagram of two-zone structure of goaf煤层开采过程中，工作面的持续推进使得顶板悬顶跨度逐渐增大，挠曲变形逐渐加剧，直接顶岩层下沉并发生断裂、垮落、充填；随着工作面的持续推进，上覆岩层裂隙逐渐发育扩展，最终形成垮落带、裂隙带和弯曲下沉带；裂隙带位于垮落带之上，关键层切分成块且整体跨落，跨落带内的岩体破碎成块度不一的散体，并发生堆积，体积较原体积有所增大，使其内部具有大量的孔隙结构，从而发挥主要的储水作用 9-1 。中硬覆岩下较厚煤层综放开采形成的垮落带及导水裂隙带高度可由式(1)和式(2)计算 1 2 :100MHk=4.71(1)0.49M+19.12100MHi=11.49(2)0.26M+6.88式中：H为垮落带高度，m；M为煤层厚度，m；H 为裂隙带高度，m。跨落带未完全充水及完全充水状态示意图如图2。顶板落带媒层（a）垮落带未完全充水状态顶板跨落带媒层（b）垮落带完全充水状态图2垮落带未完全充水及完全充水状态示意图Fig.2Schematic diagrams of caving zone not fully filledwithwaterand fullyfilledwith water矿井发生突水时，水位上升至煤层底板标高后水流涌人工作面采空区，并主要对垮落带内的孔隙结构进行充填，此过程即为采空区的充水过程；由于内部研石的阻滞作用，水流在采空区的流动情况将与其在巷道等无阻碍空间的流动情况有着显著不同，这就导致采空区出现外部水位快速上升，而内部水位低于外部水位，且自由水面出现坡度的情况；可以预知的是，在充水完成后，水面坡度减小，但该过程的具体水流蔓延情况、水位变化特征、突水强度以及空间结构对该过程的影响尚不明确。1.2研究背景选取璐安集团高河能源E1305工作面为研究背景，该工作面采用走向长壁、后退式综放一次采全高采煤法，主采煤层为二叠系山西组3 煤层，煤层回采厚度稳定且构造简单；一般倾角为315，平均倾角为5，可视作近水平煤层。3 煤层顶板、底板岩层柱状图如图3。工作面落带最高高度为3 3.5 m，破碎的岩石散体体积相较原本岩体明显变大的现象即为岩体的碎胀特性，常用破碎后体积与原体积之比，即碎胀系数表示 1 3 。岩体的碎胀系数满足式(3)：VRIkp=VRo(3)式中：kp为碎胀系数；Vr为破碎后体积，m；Vro为原体积,m。182SafetyinCoal MinesSep.20232023年9 月煤防发全No.9Vol.54第9 期第5 4 卷层序均厚/m柱状岩性713.3粉砂岩：局部夹砂质泥岩67.2砂质泥岩：夹细砂条带51.9泥岩：水平纹理，含碳屑46.43#煤层：主采煤层32.9粉砂岩28.9细砂岩117.1砂质泥岩图3 3*煤层顶板、底板岩层柱状图Fig.3Column diagram of roof and floor of No.3 coal seam岩体的碎胀受到埋深、岩性、厚度等条件的多方面影响，对于由下及上为泥岩、砂质泥岩、砂岩这类分布情况的顶板，覆岩垮落后下部岩石散体的块度相对较小，而碎胀系数则相对较大 1 4 ；结合E1305工作面顶板岩性，泥岩、砂质泥岩和粉砂岩是垮落带的主要堆积岩石。1.3实验装置为高效模拟采空区的充水过程，简化实验复杂程度，将所研究采空区以外的部分简化为1 个整体结构，使用底部尺寸为2 0 0 mm200mm的亚克力材质水箱作为替代物。水箱一侧外壁由多个可活动挡板构成，能够通过增减挡板数量实现高度调节，以对不同的巷道空间结构高度进行模拟；其余组成包括采空区模型、供水装置、数据采集及分析仪器等；采空区模型使用亚克力板材搭建并进行防漏水处理，采空区与水箱由单侧巷道相连接，整体几何相似比为1：2 0 0；结合E1305工作面走向长度实际情况，主要研究采空区走向方向的水位变化情况，对工作面倾向长度即模型宽度进行一定了取舍；采空区尺寸（长宽高）为2 0 0 0 mm220mm150mm，工作面尺寸（长宽高）为2 2 0 mm39mmx18mm，进（回）风巷的尺寸（长宽高）为2 0 0mm25mm18mm；根据上文垮落带高度计算结果，并参照顶板岩性及模型参数进行岩石铺设，从下到上分别为泥岩层、砂质泥岩层和粉砂岩层，共计1 5 0 mm。模型的岩层铺设情况见表1。实验过程中，使用水泵以某流量值向水箱供水使水位上升、水流通过巷道流人工作面采空区，多余水流将在到达水箱预设高度后溢出以维持该表1模型的岩层铺设情况Table 1 Laying of model rock formations层序铺设岩石粒径/mm铺设厚度/mm3粉砂岩2030852砂质泥岩1020501泥岩51015水位；为达到多测点数据监测的目的，选用多个液位传感器及配套装置进行数据采集及记录，液位传感器测量范围为2 0 0 mm，外径为5 mm，测量精度为0.5%FS；依据模型大小及研究需要，将9 个传感器居中布置并设置测点间距为2 5 0 mm，近工作面处为1#测点，并以此类推，每秒记录1次水位数据。2实验结果2.1采空区内部水位变化特征通过分析各测点的水位上升曲线研究水位变化特征，将水量设置为1 7 0 0 L/h，水箱高度设置为1 5 0 mm。不同测点的水位随时间变化曲线如图4。t-26st-141 s1601II120测点1测点280测点3测点440测点5测点6测点7050100150时间/s图4 不同测点的水位随时间变化曲线Fig.4Curves of water level changing with time atdifferent measuring points由图4 可知：从传感器探测水流到水位上升至垮落带高度，不同测点处水位随时间的变化具有相似的趋势；但在充水初期，不同位置水位上升的起始时间有先后差别；近工作面端1 测点水位最先上升后，其余测点依次开始变化，出现这种现象的原因在于，在采空区充水过程中，石对水流向远端的流动起到阻滞作用，使得水流到达深部不同位置的时间有着明显的先后差异；7#183SafetyinCoalMines2023年9 月Sep.2023煤防岁全Vol.54No.9第9 期第5 4 卷测点处距离工作面更远，因此其水位最后上升，同时，从曲线的变化趋势可以看出，近工作面端水位初期上升快速而后期速度有所下降，远端水位虽起始时间较晚，但曲线更为快速平稳；在充水末期，两端水位达到一致高度即完成充水过程。结合各测点的水位上升曲线，可将采空区充水过程分为2 个不同的阶段：水流向深部蔓延阶段：此阶段研石将阻滞水流向远端的流动，工作面端水位持续升高的同时水流依次流过1#7 测点处；后续充水：完成水流向深部的蔓延后，近工作面端水位明显高于远端水位，在此阶段中，近工作面端水位上升速度放缓，而远端水位快速上升，直至完成充水过程。为研究充水过程的具体发展情况，取1#7 测点的水位测量值均值为某时刻的平均水位高度，通过平均水位高度计算可得充水完成率，计算公式如下：Volk=(4)VGO式中：k为充水率；Vai为已充水的采空区体积，m；V c o 为采空区总体积，m。随着充水过程的进行，采空区内部逐渐被水淹没，直观体现在充水率随时间的变化而逐渐增长，并绘制出变化曲线。充水率随时间变化曲线如图5。1.00,J(-3.5011710-5)x2+0.012.2x-0.014 35R2=0.998 90.5充水率拟合曲线0255075100125150时间图5 3充水率随时间变化曲线Fig.5Curveofwaterfillingratewithtime由图5 可知：在充水初期，完成率曲线上升趋势明显，随着充水过程的继续，上升速度有所放缓，并在后期呈现缓慢上涨趋势；原因在于随着充水过程的进行，采空区内部水位逐渐升高，导致进水压力相对降低，充水速率相对下降，充水完成率的变化曲线呈近似二次函数形式。采空区内水面的坡度将随时间发生变化，在初期蔓延阶段近工作面端水位持续上升，水流在向采空区深部流动的过程中，水面与水平地面形成的夹角随时间变化；在后续充水阶段中，由于不同位置的水位上升速率不同，水面夹角仍将继续变化；水位夹角正切值的大小能够有效反映水面坡度大小，水面夹角正切值随时间的变化趋势则能够有效反映该过程中的水位坡度变化情况夹角正切值随时间变化曲线如图6。0.08y=（5.6 1 0 3 3 1 0-）x 2 0.0 0 1 4 1 x+0.0 8 9 5R2=0.9950.060.04夹角正切值拟合曲线0.020255075100125150时间/s图6 夹角正切值随时间变化曲线Fig.6Curve of tangent value of included angle changingwith time由图6 可知：充水初期充水速率较快，水面坡度变化较为强烈，正切值随时间快速降低；随着充水过程后期充水速率的降低，水位夹角正切值的减小速度也相应放缓，并在充水完成时到达最小值；比较夹角正切值的变化趋势与充水率的变化趋势，可以发现两者具有良好的相关性。结合上述变化过程可以认为：在整个充水过程中，水面总是呈现近工作面端水位较高而远工作面端水位较低的坡度情况，随着充水过程的逐步完成，自由水面坡度由陡峭而逐渐平缓。2.2影响因素井下空间的复杂性不仅体现在突水点的不确定性上，也体现在突水后井下巷道空间结构对于矿井水位上升的影响；对于单个采空区而言，突水后与该采空区相连巷道的空间结构将决定突水后采空区人口水位的高低，这是影响充水过程的重要因素；在此次实验中，通过改变水箱高度来实现变量控制。设置水量为2 0 0 0 L/h，不同水箱高度下采空区完成充水用时见表2。由表2 可知：随着水箱高度的增加，完成采空区充水所用的时长减少；以水箱高度1 5 0 mm时184SafetyinCoalMinesSep.20232023年9 月No.9Vol.54煤砺发全第9 期第5 4 卷表2不同水箱高度下采空区完成充水用时Table 2Water filling time at different heights ofwater tank水量/(Lh)水箱高度/mm用时/s降幅/%1501280200017510914.82009823.4所用时长1 2 8 s为基准，提升至1 7 5 mm后，用时缩短了1 4.8%；继续提高至2 0 0 mm后，则缩短了23.4%。在重力作用下，更高的入口水位将使进水压力相对更大，充水的速度相对更快。为研究不同水箱高度条件对于充水过程的影响，绘制各条件下的充水率随时间变化曲线，在充水时间相同时，水箱高度2 0 0 mm的实验组充水率最高，1 7 5 mm条件下次之，1 5 0 mm条件下则最低。不同水箱高度条件下的充水率变化曲线如图7。1.00.5一2 0 0 mm-175 mm-150mm0306090120150时间/s图7不同水箱高度条件下的充水率变化曲线Fig.7Variation curves of water filling rate fordifferent tank heights由图7 可知：较早阶段3 条曲线的一致率较高，这主要是因为在充水初期，外部水位的高度尚未达到预设条件而处于上升过程，充水速度相近导致的；在充水时间较短时，充水率极为相近，而在后期则出现水箱高度越高充水率越高的现象。对充水过程后期不同充水率的所用时长进行定量分析，以水箱高度1 5 0 mm实验组为基准，研究在不同条件下完成同样充水率时所用时长的降幅特征。不同充水率所需时长见表3，依据表中数据得到以低水箱高度条件为基准，不同水箱高度所需时长的降幅曲线如图8。表3不同充水率所需时长Table3Time required for different filling rates时长/s水箱高度/mm50%60%70%80%90%100%1504458718499128175425365768810920041516373839825.0-200 mm20.0-175 mm%/刺斗15.010.05.02.5.0.40.50.60.70.80.91.0充水率图：不同水箱高度所需时长的降幅曲线Fig.8Curves of decreasing time required fordifferent tank heights以水箱高度1 5 0 mm实验组为参照，相比之下，水箱高度1 7 5 mm的实验组所用时长的降幅呈现增长趋势；充水率为0.5 时，降幅仅为4.5%；到充水完成后，降幅达到最大值1 4.8%；水箱高度2 0 0mm实验组具有同样的变化趋势，充水率为0.5 时，降幅为6.8%；到充水完成后，降幅达到最大值23.4%，且2 0 0 mm实验组的降幅始终大于1 7 5 mm实验组；期间虽然出现个别异常点，但整体上升趋势明显。从整个变化过程可以看出，水箱高度越高，完成同样充水效果所需的时间越短；前期相同充水率所用时长的降幅很小，随着充水过程的进行，降幅逐渐变大，降幅曲线呈上升趋势。矿井突水水量的大小是影响矿井水位发展过程的重要因素。将影响采空区的充水过程直观反映在充水用时的长短，对不同水量下充水所用时长进行定量分析，以此来研究突水量对于充水过程的影响。将水箱高度设置为1 5 0 mm，不同水量下采空区完成充水用时见表4，不同水量条件的充水率变化曲线如图9。由表4 可知：随着水量的增加，采空区充水所用的时长有所减少，这意味着大水量将使采空185Safety in Coal MinesSep.20232023年9 月No.9Vo1.54煤砺发全第9 期第5 4 卷表4不同水量下采空区完成充水用时Table4Water filling time at different water volumes水箱高度/mm水量/(Lh)用时/s降幅/%1400153015017001417.8200012816.31.00.52 000 L/h1 700 L/h-1 400 L/h04080120160时间/s图9不同水量条件的充水率变化曲线Fig.99Variation curves of water filling rate underdifferent water volumes区充水速度有所加快。以小水量用时1 5 3 s为基准，当提升至中水量后，用时缩短了7.8%；提升至大水量后，用时则缩短了1 6.3%。由图9 可知：水量越大，相同充水时间后的充水完成率越高，且在整个过程中，2 0 0 0 L/h条件下的充水率均大于1 7 0 0 L/h条件下的充水率，1400L/h条件下的充水率则始终最低。不同水量条件下完成相同充水效果所需的时间差异较大，对不同充水率的所用时长进行定量分析，并以水量1 4 0 0 L/h实验组为基准，研究所用时长的降幅特征。不同充水率所需时长见表5，不同水量所需时长的降幅曲线如图1 0。表5不同充水率所需时长Table5Time required for different filling rates水量/时长/s(L:h)10%20%30%40%50%60%70%80%90%100%14001323324255698398116153170012202837486276 90107 1412000101826334458718499128由图1 0 可知：不同水量条件下的时长降幅较为稳定，曲线呈现非增长趋势；以水量1 4 0 0 L/h302000L/h1700 L/h%/刻斗201000.20.40.60.81.0充水率图1 0不同水量所需时长的降幅曲线Fig.10Curves of decreasing time required fordifferent watervolumes实验组为参照，对于水量1 7 0 0 L/h的实验组，充水率为0.1 时所用时长降幅为7.7%，充水完成时的降幅也仅为7.8%；水量2 0 0 0 L/h实验组具有同样的变化趋势，充水率为0.1 时所用时长降幅为23%，充水完成时降幅为1 6.3%，且在整个过程中，大水量条件下的降幅始终大于中等水量实验组；从整个变化过程可以看出，突水量越大，采空区的充水速率越快，完成同样充水率所需的时长越短，且不同充水率所用时长的降幅较为一致，时长降幅曲线呈平稳趋势。3结语1)研石对水流的阻滞作用将使充水过程包含初期水流蔓延阶段和后续充水阶段，该过程中充水完成率随时间逐渐增长并近似二次函数变化；该过程中近工作面端处水位较高而远端处水位较低，水面出现坡度且随时间逐渐平缓，水面夹角正切值呈现非线性变化。2)矿井空间结构和突水量是影响采空区充水速率的重要因素。通过改变水箱高度研究矿井空间结构的影响，发现提高水箱高度将缩短充水所用时长，以水箱高度1 5 0 mm实验组为基准，1 7 5、200mm实验组用时分别缩短了1 4.8%和2 3.4%；大突水量也将加快充水速率，以水量1 4 0 0 L/h为基准，1 7 0 0、2 0 0 0 L/h实验组的用时则分别缩短了 7.8%和 1 6.3%。参考文献(References):1武强.我国矿井水防控与资源化利用的研究进展、问186SafetyinCoalMinesSep.20232023年9 月煤防发全No.9Vol.54第9 期第5 4 卷题和展望 J.煤炭学报，2 0 1 4,3 9(5)：7 9 5-8 0 5.WU Qiang.Progress,problems and prospects of preven-tion and control technology of mine water and reutiliza-tion in ChinaJ.Journal of China Coal Society,2014,39(5):795-805.2李文平，乔伟，李小琴，等.深部矿井水害特征、评价方法与治水勘探方向 J.煤炭学报，2 0 1 9，4 4(8）：2437-2448.LI Wenping,QIAO Wei,LI Xiaoqin,et al.Characterist-ics of water disaster,evaluation methods and explora-tion direction for controlling groundwater in deep min-ingJ.Journal of China Coal Society,2019,44(8):2437-2448.3张培森，朱慧聪，李复兴，等.2 0 0 8 2 0 1 9 年我国煤矿水害事故统计及演变趋势分析 J.煤矿安全，2 0 2 1，52(8):194-200.ZHANG Peisen,ZHU Huicong,LI Fuxing,et al.Evolu-tion trend and statistical analysis of coal mine water dis-aster accidents in China from 2008 to 2019J.Safety inCoal Mines,2021,52(8):194-200.4马恒，刘超.矿井突水路径搜索算法与应用 J.中国安全生产科学技术，2 0 1 5,1 1(1 1)：3 5-4 0.MA Heng,LIU Chao.Searching algorithm of mine wa-ter inrush paths and its applicationJ.Journal of SafetyScience and Technology,2015,11(11):35-40.5蔡明杰，朱希安，刘德民，等.基于优化SPFA算法的矿井突水救援模型 J.煤田地质与勘探，2 0 1 9,4 7(6)：78-83.CAI Mingjie,ZHU Xian,LIU Demin,et al.Mine wa-ter inrush rescue model based on optimized SPFA al-gorithmJ.Coal Geology&Exploration,2019,47(6):78-83.6于丹，颜伟，李劭昱.基于权值时变模型的矿井突水最优逃生路径的动态选择 J.科学技术与工程，2 0 2 2,22(12):4762-4771.YU Dan,YAN Wei,LI Shaoyu.Dynamic selection ofoptimal escape path of mine water inrush based onweight time-varying modelJ.Science Technology andEngineering,2022,22(12):4762-4771.7周建军，虎维岳，张壮路.废弃矿井地下水回灌模拟分析研究 J.西安科技大学学报，2 0 0 8,2 8（3)：4 3 4-4 3 8.ZHOU Jianjun,HU Weiyue,ZHANG Zhuanglu.Simu-lation of grounder water rebound in abandoned mineJ.Journal of Xian University of Science and Technology,2008,28(3):434-438.8周建军，虎维岳，侯大勇.废弃矿井地下水淹没过程的水流与水位数值模拟 J.煤田地质与勘探，2 0 1 1,39(4):28-31.ZHOU Jianjun,HU Weiyue,HOU Dayong.Numericalsimulation of groundwater rebound process and watertable value in abandoned minesJ.Coal Geology&Exploration,2011,39(4):28-31.9李连崇，唐春安，梁正召.考虑岩体碎胀效应的采场覆岩冒落规律分析 J.岩土力学，2 0 1 0，3 1(1 1)：3537-3541.LI Lianchong,TANG Chunan,LIANG Zhengzhao.In-vestigation on overburden strata collapse around coalface considering effect of broken expansion of rockJ.Rock and Soil Mechanics,2010,31(11):3537-3541.1 0 张海涛，许光泉.采空区冒落结构与充水性特征 .煤田地质与勘探，2 0 1 8,4 6(4)：9 9-1 0 2.ZHANG Haitao,XU Guangquan.Research on cavingstructure and water filling characteristic of goafJ.Coal Geology&Exploration,2018,46(4):99-102.1 1 梁冰，汪北方，姜利国，等.浅埋采空区垮落岩体碎胀特性研究 J.中国矿业大学学报，2 0 1 6，4 5(3)：475-482.LIANG Bing,WANG Beifang,JIANG Liguo,et al.Broken expand properties of caving rock in shallow bur-ied goafJ.Journal of China University of Mining&Technology,2016,45(3):475-482.【1 2 许延春，李俊成，刘世奇，等.综放开采覆岩“两带”高度的计算公式及适用性分析 J.煤矿开采，2 0 1 1,16(2):4-7.XU Yanchun,LI Juncheng,LIU Shiqi,et al.Calcula-tion formula of“two-zone”height of overlying strataand its adaptability analysisJ.Coal Mining Techno-logy,2011,16(2):4-7.【1 3 苏承东，顾明，唐旭，等.煤层顶板破碎岩石压实特征的试验研究 J.岩石力学与工程学报，2 0 1 2,3 1(1)：18-26.SU Chengdong,GU Ming,TANG Xu,et al.Experi-ment study of compaction characteristics of crushedstones from coal seam roofJ.Chinese Journal of RockMechanics and Engineering,2012,31(1):18-26.1 4 马资敏，郭志飚，刘晓宇，等.切顶留巷冒落研石碎胀规律与力学特性研究 J.煤炭工程，2 0 1 7,4 9(1 1)：94-98.MA Zimin,GUO Zhibiao,LIU Xiaoyu,et al.Study onbulking rule and mechanical characteristics of fallengangue in roof-cutting entry retainingJ.Coal Engin-eering,2017,49(11):94-98.