Study on the overburden failure features and microseismic measurements in non-pillar gob-side entry retaining by roof cutting
-
摘要: 为进一步研究无煤柱切顶留巷技术开采后的覆岩破坏规律,以柠条塔煤矿S1201−Ⅱ工作面为工程背景,采用物理相似模拟与数值模拟的研究手段,结合现场微震监测技术建立了微震波形数据库,研究了随工作面持续开采,无煤柱切顶留巷不同阶段的覆岩采动裂隙演化及应力空间展布特征,得出了工作面覆岩周期性破断规律。研究结果表明:工作面发生初次来压时的覆岩裂隙发育高度为57.6 m,切顶前中部裂隙带发育高度为95.5~96.1 m,裂采比为23.8~24.0,边缘侧裂隙发育高度为105.9~106.4 m,裂采比为26.4~26.6。切顶后工作面两侧裂隙带最终发育高度为104.3~105.2 m,裂采比为26.1~26.3,工作面中部裂隙带由于上覆岩层的不断压实弥合,最终发育高度为94.3~95.2 m,裂采比为23.6~23.8。当巷道分别处于掘进、切缝阶段,顶板位移基本没有产生改变;当其进入顶板下沉、切顶成巷阶段,顶板位移不断增大。切顶卸压完成后,巷道侧支承压力峰值增大,表明切缝之后的工作面跨度进一步增大,倾向支承压力不断增大;工作面顶板卸压效果显著,顶板产生大范围应力释放现象。在该工作面布置了微震监测系统,发现微震事件的周期性产生与工作面周期来压有强关联性,其发展过程可划分为萌芽期—发展期—高潮期,进一步综合得出覆岩的周期性破断演化规律。Abstract: In order to further study the failure law of overburden after the mining of non-pillar gob-side entry retaining by roof cutting technology, taking the S1201-II working face of Ningtiaota Coal Mine as the engineering background, physical similarity simulation and numerical simulation research methods are used. Combined with on-site microseismic monitoring technology, a microseismic waveform data library is established. With continuous mining of the working face, the evolution of overburden mining induced cracks and stress spatial distribution features at different stages of non-pillar gob-side entry retaining by roof cutting are studied. The periodic crack law of the overburden in the working face has been obtained. The research results show that the height of the overburden cracks during the initial pressure on the working face is about 57.6 m, the height of the middle crack zone before cutting is 95.5-96.1 m, the crack mining ratio is 23.8-24.0, the height of the edge side cracks is 105.9-106.4 m, and the crack mining ratio is 26.4-26.6 m. After the roof cutting, the final development height of the crack zone on both sides of the working face is about 104.3-105.2 m, with a crack mining ratio of 26.1-26.3. Due to the continuous compaction and closure of the overburden layer, the final development height of the crack zone in the middle of the working face is 94.3-95.2 m, with a crack mining ratio of 23.6-23.8. When the roadways are in the excavation and cutting stages respectively, there is basically no change in the displacement of the roof. When it enters the sinking and roadway formation stage, the displacement value of the roof continuously increases. After the completion of roof cutting and pressure relief, the peak support pressure on the side of the roadway increases, indicating that the span of the working face further increases after the cutting seam, and the inclined support pressure continues to increase. The pressure relief effect of the working face roof is significant, and the roof produces a large-scale stress release phenomenon. A microseismic monitoring system is installed in the working face, and it is found that there is a strong correlation between the periodic occurrence of microseismic events and the periodic pressure of the working face. The development process can be divided into the budding stage, development stage, and climax stage. Further comprehensive analysis can be conducted to obtain the periodic crack evolution law of the overburden.
-
0. 引言
煤矿巷道的形变和破坏是煤矿生产实践中需要着重考虑解决的问题。煤矿巷道形变是指巷道岩层在受到外力因素作用下产生的形状和尺寸的变化,巷道形变会造成井下交通中断、生产设备损坏、人员伤亡等事故,严重威胁煤矿正常生产。因此,在煤矿开采过程中,及时、科学、准确地监测巷道形变,掌握巷道状态及形变趋势,对保障煤矿安全生产具有十分重要的意义[1]。
随着计算机技术的快速发展,多种数值模拟技术被广泛应用于巷道形变的研究中,如有限元法、离散元法、边界元法、有限差分法、十字观测法等。王伟等[2]运用理论计算与数值模拟分析,分析了动压巷道的应力分布规律,为巷道的布设位置选取与支护提供了指导。姜耀东等[3]对开滦矿区赵各庄矿、唐山矿深部开采过程中巷道形变、破坏特征进行了研究,采用FLAC2D数值模拟软件研究地应力与巷道岩层之间的密切关系。余伟健等[4]利用离散裂隙网络在FLAC3D软件中对矿山现场调查得到的裂隙进行了重构还原,研究分析了裂隙发育岩体巷道及典型支护结构体的变形破坏特征。杨景贺[5]采用离散元数值模拟方法,建立了不同支护条件下的高应力软岩巷道模型,研究了高应力软岩巷道变形破坏和支护加固的机理。
目前,针对巷道形变破坏的研究以现场实测为主,通过实测可掌握巷道形变破坏情况,直观了解导致巷道形变破坏的因素及作用效果。叶美图等[6]针对现有矿井巷道表面位移十字交叉测量方法存在的测量不便、耗时长、工作强度大、误差不可控等问题,设计了一种矿井巷道表面位移激光测量装置。该装置基于激光测距原理,利用激光测距仪完成十字定心和变形测量,通过单个激光测距仪可分别测得巷道中每个测试站顶板、底板、左帮和右帮4个方向的表面位移量。井下试验结果表明,该装置实现了无绳测量,测量效率高,结果准确。宋康磊等[7]采用地质雷达测定巷道松动圈,研究了极薄煤层跨采条件下的松动圈大小及巷道形变规律。荣耀等[8]提出了综采工作面三维激光扫描建模总体方案,对煤壁与顶板交线提取、标靶球检测、点云拼接及坐标转换等关键技术进行了研究,实现了三维地质坐标系下煤壁与顶板交线信息的近实时获取,该信息可直接发送给采煤机滚筒,为采煤机下一刀截割提供数据参考。杜江丽等[9]提出了多边形分割法,将模型的重点关注区域划分为多个多边形,对矿区边坡形变监测数据进行研究。亓玉浩等[10]提出了一种综采工作面实时三维建图方法,可满足综采工作面监控的精度需求。刘晓阳等[11]针对岩性变化和构造的分析监测越来越困难的问题,研究了当前巷道顶板稳定性监测方法。王峰[12]提出了以煤层透明化、综采装备透明化、决策及控制透明化为核心的基于透明工作面的智能化开采概念。邹筱瑜等[13]面向移动机器人在煤矿巷道的建图需求,提出了一种基于集成式因子图优化的煤矿巷道移动机器人三维地图构建方法。高晓进等[14]提出了软岩巷道超前爆破预裂顶板防控技术,以解决软岩巷道结构的形变问题。
上述方法为巷道形变监测提供了可行的解决方案,但还存在如下问题:① 单一的形变量指标难以准确反映巷道形变演化规律及围岩状态与发展趋势。② 巷道形变具有区域性、连续性、突发性的特点,现有方法存在监测站点单一、测点有限的问题,难以实现大范围、多站点、区域性同步测量。因此需要选取多个指标建立巷道形变监测指标体系,且需要研究一种具有区域性、连续性、实时性特点的巷道形变监测方法,及时、准确、全面地监测巷道形变以确保煤矿开采安全。
三维激光扫描技术可快速获取被测物体表面三维坐标,点云数据精度高、密度大,包含反射强度等信息[15],可实现对巷道形变连续监测。因此,针对煤矿巷道形变监测面临的问题,本文提出了一种基于三维激光扫描的煤矿巷道形变监测方法。首先,利用深度学习模型VoxelNet对三维激光扫描数据进行检测去噪。然后,利用alphashape算法对断面离散点进行拟合,并采用基于差值的方法,对数据进行多维差值计算,得到巷道形变的具体数据,实现矿区巷道形变监测的全覆盖。最后,利用中煤大同能源有限责任公司塔山煤矿为期1个月的实际监测数据,对巷道形变进行截面分析和三维整体分析。
1. 煤矿巷道形变监测方法
1.1 三维激光扫描技术
三维激光扫描技术是一种先进的全自动、高精度立体扫描技术,具有操作简单、速度快、灵活性高、精度高等特点[16],数字化、自动化、实时性强,解决了目前空间信息技术发展的瓶颈。
三维激光扫描技术借助水平方位偏转控制器及对应的高度角偏转控制器科学控制反射棱镜的实际转动,使激光测距系统发射的激光能够顺着不同的坐标轴实时移动并扫描测量。扫描仪自身定义的坐标系统中, X轴在横向扫描面内,Y轴在横向扫描面内与X轴垂直,Z轴与横向扫描面垂直。借助测得的三维扫描仪中心到目标点的斜距S、激光束水平方向偏转角a和竖直方向偏转角b来校准并计算激光脚点的实际三维坐标P(X,Y,Z)[17]。
$$ X = S \cos b \cos a $$ (1) $$ Y = S \cos b \sin a $$ (2) $$ Z = S \sin b $$ (3) 三维激光点坐标计算原理如图1所示,通过测量大量物体表面点位的三维坐标,读取并集成现场物体的反射率和纹理等信息,生成被测物点云模型,获取高精度、高分辨率的三维巷道点云数据。
1.2 点云去噪
对三维巷道点云数据去噪是煤矿巷道形变监测得以有效应用的重要预处理步骤[18]。将噪点设为巷道内的人物和内壁附着的螺丝钉和管道等,通过VoxelNet网络检测出噪点后,再对数据进行分离。
首先使用点云标注工具Semantic Segmentation Editor制造数据集,标注人物和管道等信息后送入VoxelNet训练。然后使用VoxelNet中的Voxel将三维空间中的一定空间划分为一个格子,并使用PointNet网络对格子的数据进行特征提取。最后用提取出来的特征代表这个格子,并放回到三维空间。将无序的点云数据变成高维特征数据,并使用三维卷积抽取有序的三维Voxel数据。
1.3 点云旋转
给定2个来自不同坐标系的三维数据点集,找到2个三维数据点集的空间变换关系,使2个三维数据点集处于同一坐标系中,即可实现三维点云配准[19]。配准目标是在全局坐标框架中找到单独获取视图的相对位置和方向,使2个三维数据点集之间的相交区域完全重叠。通过求点云M和点云N 之间的旋转矩阵和平移向量,使点云M和点云N的坐标处于同一坐标系下。通过欧拉角生成旋转矩阵,将点云数据进行平移和旋转。设新的空间直角坐标轴分别为X',Y',Z',绕X',Y',Z' 3个轴旋转的角度分别为α,β,γ,绕X',Y',Z'旋转分别得到的旋转矩阵为
$$ {{\boldsymbol{R}}_{X'}}\left( \alpha \right) = \left[ \begin{gathered} 1\quad \quad 0\quad \quad 0 \\ 0\quad \cos \alpha \quad - \sin \alpha \\ 0\quad \sin \alpha \quad \cos \alpha \\ \end{gathered} \right] $$ (4) $$ {{\boldsymbol{R}}_{Y'}}\left( \beta \right) = \left[ \begin{gathered} \cos \beta \quad \;0\quad \sin \beta \\ \quad 0\quad \quad 1\quad \;\;0 \\ - \sin \beta \;\;\;0\quad \cos \beta \\ \end{gathered} \right] $$ (5) $$ {{\boldsymbol{R}}_{Z'}}\left( \gamma \right) = \left[ \begin{gathered} \cos \gamma \quad - \sin \gamma \;\;\;0 \\ \sin \gamma \quad \;\cos \gamma \quad \;0 \\ 0\;\;\;\;\;\;\;\;\;\;\;\;0\quad \,\,\cos \gamma \\ \end{gathered} \right] $$ (6) 对3个旋转矩阵的点云坐标进行逐元素相乘,得到最终的旋转矩阵R。
$$ {\boldsymbol{R}} = {{\boldsymbol{R}}_{Z'}}\left( \gamma \right) {{\boldsymbol{R}}_{Y'}}\left( \beta \right) {{\boldsymbol{R}}_{X'}}\left( \alpha \right) $$ (7) 1.4 点云连续断面
将巷道切分为连续的断面,根据断面对巷道进行形变分析,分别向点云的X'OY'和Y'OZ'面进行数据投影,计算当前断面的中心点$ {P^*} \left( {{x^ * },{y^ * },{{\textit{z}}^ * }} \right) $,其中,$ {x^ * },{y^ * },{{\textit{z}}^ * } $为中心点$ {P^*} $在X',Y',Z' 3个坐标轴上对应的数值。
随后,构造与Z'轴平行的平面,平面与Z'轴平行且与Y'轴相切,则该平面的法向量u=[0 1 0],该平面的方程为$ y - {y^ * } = 0 $。由于落于平面的点比较稀疏,所以赋予平面一定的厚度,将点到平面一定距离d内的点都算作该平面的点。
1.5 形变分析
在实际扫描作业中,由于障碍物遮挡或受巷道空间扫描角度的限制,会出现部分点云空洞,难以保证断面上处处都有激光斑点,所以需要对投影后的断面离散点集进行拟合。alphashape算法可控制点集边界拟合以创建非凸区域,采用该算法拟合投影后的断面离散点集,其计算圆心的原理如图2所示,计算步骤如下:
1) 选取任意一点$ p \left( {x,y} \right) $,以半径为r的圆进行滚动,在断面点集中搜索所有距离$ p $点2r以内的点的集合,形成点集Q。
2) 选取点集Q中的任意点$ {p_1} \left( {{x_1},{y_1}} \right) $,根据$ p $和$ {p_1} $2个点的坐标和半径r,计算出经过$ p $和$ {p_1} $2个点且半径为r的圆心${p_2} \left( {{x_2},{y_2}} \right) $,${p_3} \left( {{x_3},{y_3}} \right) $。
$$ \left\{ \begin{array}{l} {x_2} = x + 1/2({x_1} - x) - H ({y_1} - y) \\ {y_2} = y + 1/2({y_1} - y) - H (x - {x_1}) \\ {x_3} = x + 1/2({x_1} - x) + H ({y_1} - y) \\ {y_3} = y + 1/2({y_1} - y) + H (x - {x_1}) \\ \end{array} \right. $$ (8) 式中H为参数。
3) 计算点集Q中(除$ {p_1} $点)剩余点分别到$ {p_2},{p_3} $点的距离。若所有点到$ {p_2} $或$ {p_3} $的距离均大于半径r,则证明$ p $点为边界点。
4) 若剩余的点到$ {p_2} $或$ {p_3} $点的距离不全都大于半径r,则遍历点集Q内的所有点,将其轮换作为$ {p_1} $点,重复上述过程。若存在某点满足2)和3),则该点为边界点,终止该点的判断,继续判断下一点。若点集Q中所有临近点中均不存在$ {p_1} $这样的点,则$ p $点为非边界点。
根据alphashape算法拟合巷道断面轮廓图,构建不同监测时期同一位置的巷道断面,将其叠加后,采用对比不同时期同一位置巷道断面的高度和宽度等信息的方法,可以更直观地分析巷道形变趋势。具体分析过程:首先,由alphashapes算法计算第一期巷道断面的轮廓点,分别计算相邻2个点$p'\left( {{x_0},{y_0},{{\textit{z}}_0}} \right), q\left( {{x_1},{y_1},{{\textit{z}}_1}} \right)$所构成的线段$ {l_1} $的方向向量,默认$p',q$处于同一平面内,故只保留X轴和Z轴值,记为$p'\left( {{x_0},{{\textit{z}}_0}} \right), q\left( {{x_1},{{\textit{z}}_1}} \right)$。然后,计算过$p '$点与$ {l_1} $垂直的直线$ {l_2} $,$ {l_2} $即法向直线。最后,计算第二期与第一期相同法向直线的点的坐标,计算2个点间的距离,将其记为形变值。
2. 工程实践
2.1 工程概况
中煤大同能源有限责任公司塔山煤矿的开采工作面为30507工作面,工作面采用单巷布置方式,即在工作面分别布置胶带巷和回风巷。工作面西侧布置回风巷,东侧布置胶带巷,胶带巷东侧为实煤区,回风巷西侧为30509工作面采空区(2018年6月停采),30507工作面与30509工作面之间区段煤柱为8 m,因此将30507工作面回风巷作为监测地点,进行为期1个月的实际形变监测。
通过对30507工作面进行实地勘查,其具有以下特点:
1) 靠近工作面的回风巷区域形变量大,已有控制点坐标经过一定时间会发生变化,无法直接利用。
2) 对于井下空间地理坐标系,除用全站仪进行测点定位外,没有更好的方式将点云数据和矿井坐标系进行联系。
3) 矿区井下环境相对复杂,煤尘、空气湿度、振动等工况极易影响工作效率和数据精度。
2.2 实践方案
塔山煤矿正在回采的工作面为30507工作面,开始监测时采煤机所处位置在里程1 000 m处,监测范围为里程600~900 m。在测绘区域起始位置即里程600 m处,前后上下交错架设多个360°棱镜。根据30507回风巷实际情况,考虑到工作面附近棱镜位置会产生变化,因此以风门处于终采线以内的基准点为基准,测量各个360°棱镜的三维坐标,使棱镜与井下基准点处于同一坐标系下。采用后方交会的方式,在每次进行巷道形变监测数据采集时,以相同基准点作为初始测点的定向点,将矿井坐标系赋予点云数据,以保证每次进行形变监测工作所采集的点云数据都处于同一个空间坐标系下。
受巷道结构影响,只能选取接近于直线的布站路线进行三维激光扫描。为避免由于线性布站及井下煤尘、湿气等因素造成的拼接误差,相邻测站之间采用棱镜定位的方式进行拼接定位,选择应用可识别棱镜并记录棱镜信息的三维激光扫描仪。相邻测站间需保证有2个及以上共同的棱镜作为拼接基准。
将周期性采集的点云数据置于同一坐标系下,对多次的点云数据进行对比,得出巷道顶底板和煤柱侧、回采侧两侧壁的形变量。
2.3 数据采集
1) 控制测量。由于三维激光扫描仪具有独立的坐标系统,为了将点云数据坐标系转换为实际地理坐标系,需结合全站仪进行控制测量。在里程600 m处做3个能够固定在锚杆上的棱镜放置装置,并在其上安置3个360°棱镜。以风门处7号、9号、11号基准点为基准,采用全站仪测出3个棱镜在矿区地测坐标系中的坐标。
2) 在里程600~900 m段,按30 m的间距,在巷道两帮粉刷并做标记,作为扫描的架站位置;在里程600 m处开始架设三维激光扫描仪,并进行对中整平等操作。
3) 定位棱镜测量。扫描仪架设完毕后,使用棱镜扫描功能采集设置的360°棱镜点位数据,将定位棱镜信息同步至点云模型中。
4) 根据项目要求选择合适的扫描参数,建立测站文件夹,使用棱镜的扫描功能,识别测站间用于点云拼接所架设的棱镜,识别完成后开始三维激光扫描,获取该测站全景点云数据。
5) 检查数据质量,若数据质量合格,即可将扫描仪迁移至下一预设测站,并以同样步骤进行该测站的三维激光扫描。依次循环,直至待测区域全部扫描完成,即数据采集工作完成。
2.4 数据处理
数据处理过程只需将相邻测站间相同的棱镜目标进行匹配,即可达到点云拼接目的,过程如下:
1) 控制点坐标录入。在数据处理软件中,使用起始测站中采集的定位棱镜点创建目标,将全站仪测量的定位棱镜坐标赋予对应目标,控制点坐标录入,如图3所示。
2) 后方交会拼接。起始测站点云利用其中的定位棱镜点与已赋予现场坐标系坐标的目标点位匹配。第二测站点云以起始测站为基准,利用与起始测站共视的2个或2个以上棱镜进行后方交会,将其与起始测站拼接[20]。后续测站依次进行,得到定位于现场坐标系中的整体测区点云模型。
3) 数据去噪。为采集完整的巷道顶板、底板及煤柱侧、回采侧两侧壁数据,采用360°全景扫描仪采集现场所有可视数据。由于噪点数据[21](图4)会对实施效率及结果精度造成影响,所以需对噪点数据进行去噪,去噪后数据如图5所示。
4) 数据输出。对数据进行拼接、去噪后,在确保点云定位信息固定的前提下,导出相应格式的点云模型,如图6和图7所示。
2.5 数据分析
在数据处理过程中,将所有三维点云模型置于同一坐标系下,因此在将同一测区的2期三维点云数据导入对比分析软件后,2期三维点云模型会自动放置在同一位置。
1) 三维整体分析。将置于同一位置的2期点云数据沿其法线方向进行差值分析,得到整体形变量,以色谱的形式将差值进行显示,不同颜色的分布表示不同区域2期点云模型偏差的分布,颜色的深浅表示该区域偏差大小。通过对模型上颜色的分布和深浅分析,得到该测段巷道形变的整体趋势,如图8所示。如果巷道受到向内的挤压,会变成红色,红色越深,表明挤压程度越大;如果巷道趋近于无变化,则会变成绿色;如果巷道向外扩张,会变成蓝色,蓝色越深,表明向外扩张程度越大。
2) 截面分析。以三维整体分析为基础,在整体巷道点云模型的竖直方向上按一定距离布置截面,截取2期点云在截面上的偏差,以色谱的形式将差值进行显示,截面分析结果如图9所示。可看出离工作面距离越近,附色模型的颜色越偏向红色和蓝色,且颜色越深,代表巷道形变越大。
两期巷道截面对比如图10所示,通过2个巷道轮廓的相对偏差,可看出该区域主要存在的形变:如果上轮廓向内偏离,则顶板发生塌陷;如果下轮廓向外偏离,则底板发生鼓起。
3) 预警记录。点选预警记录可查看预警位置、里程位置和最大形变值等信息,数据精确,可随时观察巷道形变情况,以便采取相应维护措施,巷道形变预警记录如图11所示。
3. 结论
1) 将三维激光扫描技术应用于塔山煤矿30507工作面的形变监测,将人工智能与传统方法相结合,提出了基于三维激光扫描的煤矿巷道形变监测方法。利用深度学习VoxelNet模型对点云数据进行去噪,利用alphashape算法对断面离散点进行拟合,并采用基于差值的方法,将数据进行多维差值计算以得到巷道形变数据。
2) 利用塔山煤矿的实际监测数据,对巷道三维点云数据进行了截面分析和三维整体分析。截面分析主要基于巷道三维点云按固定间距选取截面,比对2期点云截面数据,通过截面轮廓差异判断巷道形变情况;三维整体分析利用2期巷道三维点云数据进行整体比对,通过色谱显示比对结果,分析巷道形变区域和变化趋势。随着测点距工作面的距离越近,附色模型越偏向红色和蓝色,且颜色越深,表明巷道形变越大。
-
表 1 模拟实验相似材料配比方案
Table 1 Simulation experiment for similar material matching scheme
岩层 厚度/cm 材料质量/kg 河砂 石膏 大白粉 粉煤灰 细粒砂岩 6 8.40 0.48 0.72 — 砂质泥岩 2 8.64 0.29 0.67 — 细粒砂岩 8 8.40 0.48 0.72 — 粗粒砂岩 2 8.64 0.29 0.67 — 粉砂岩 16 8.53 0.21 0.86 — 中粒砂岩 14 8.64 0.29 0.67 — 细粒砂岩 6 8.40 0.48 0.72 — 1−2上煤 3 1.96 0.10 0.49 1.96 粉砂岩 11 8.53 0.21 0.86 — 细粒砂岩 9 8.40 0.48 0.72 — 中粒砂岩 19 8.64 0.29 0.67 — 2−2煤 4 1.96 0.10 0.49 1.96 粉砂岩 16 8.53 0.21 0.86 — 表 2 煤岩体物理力学参数
Table 2 Mechanical parameters of coal and rock
岩石 密度/
(kg·m-3)体积模
量/MPa剪切模
量/MPa抗拉强
度/MPa黏聚
力/MPa内摩擦
角/(°)细砂岩 2 420 4 167 2 869 1.8 3.5 37 中粒砂岩 2 550 3 435 3 876 1.9 3.9 37.5 煤层 1 350 2 381 1 163 0.6 1.3 32.9 粉砂岩 2 530 3 372 3 816 1.7 4.5 20.4 细粒砂岩 2 640 9 302 9 137 2.1 4.2 28 表 3 周期来压步距和周期来压强度统计
Table 3 Statistical of the periodic weighting step and pressure strength
来压时简( 年−月−日) 工作面上端头侧 工作面中部 工作面下端头侧 来压强度/MPa 来压步距/m 来压强度/MPa 来压步距/m 来压强度/MPa 来压步距/m 2018−12−12 29.24 15.0 38.6 17.3 39.2 15.6 2018−12−15 31.6 16.3 36.9 14.8 36.4 17.4 2018−12−17 28.6 15.8 38.6 18.9 37.4 21.3 2018−12−19 29.3 20.3 37.7 26.0 36.5 22.9 2018−12−22 32.6 22.7 37.4 25.0 37.5 23.0 2018−12−24 35.7 24.0 37.5 24.7 38.9 24.7 2018−12−26 39.8 24.8 38.2 23.1 34.9 15.7 2018−12−28 34.7 19.6 38.1 17.0 35.6 21.4 2019−01−02 35.0 15.6 39.5 18.2 36.8 16.9 2019−01−04 37.6 17.0 37.9 18.9 37.6 15.4 2019−01−06 36.4 15.0 36.8 17.3 34.8 15.6 -
[1] 何满潮,宋振骐,王安,等. 长壁开采切顶短壁梁理论及其110工法——第三次矿业科学技术变革[J]. 煤炭科技,2017(1):1-9,13. HE Manchao,SONG Zhenqi,WANG An,et al. Theory of longwall mining by using roof cuting shortwall team and 110 method-the third mining science and technology reform[J]. Coal Science & Technology Magazine,2017(1):1-9,13.
[2] 张国锋,何满潮,俞学平,等. 白皎矿保护层沿空切顶成巷无煤柱开采技术研究[J]. 采矿与安全工程学报,2011,28(4):511-516. ZHANG Guofeng,HE Manchao,YU Xueping,et al. Research on the technique of no-pillar mining with gob-side entry formed by advanced roof caving in the protective seam in Baijiao Coal Mine[J]. Journal of Mining & Safety Engineering,2011,28(4):511-516.
[3] 赵萌烨,黄庆享,朱磊,等. 无煤柱切顶沿空留巷顶板“短砌体—铰接”结构及支撑阻力研究[J]. 煤炭学报,2021,46(增刊1):84-93. ZHAO Mengye,HUANG Qingxiang,ZHU Lei,et al. Research on roof structure and support resistance of gob-side entry retaining with roof cutting and non-pillar mining[J]. Journal of China Coal Society,2021,46(S1):84-93.
[4] 华心祝. 我国沿空留巷支护技术发展现状及改进建议[J]. 煤炭科学技术,2006,34(12):78-81. HUA Xinzhu. Development status and improved proposals on gob-side entry retaining support technology in China[J]. Coal Science and Technology,2006,34(12):78-81.
[5] 王亚军,何满潮,张科学,等. 无煤柱自成巷开采巷道矿压显现特征及控制对策[J]. 采矿与安全工程学报,2018,35(4):677-685. WANG Yajun,HE Manchao,ZHANG Kexue,et al. Strata behavior characteristics and control countermeasures for the gateroad surroundings in innovative non-pillar mining method with gateroad formed automatically[J]. Journal of Mining & Safety Engineering,2018,35(4):677-685.
[6] 朱珍,何满潮,王琦,等. 柠条塔煤矿自动成巷无煤柱开采新方法[J]. 中国矿业大学学报,2019,48(1):46-53. ZHU Zhen,HE Manchao,WANG Qi,et al. An innovative non-pillar mining method for gateroad formation automatically and its application in Ningtiaota Coal Mine[J]. Journal of China University of Mining & Technology,2019,48(1):46-53.
[7] 何满潮,马资敏,郭志飚,等. 深部中厚煤层切顶留巷关键技术参数研究[J]. 中国矿业大学学报,2018,47(3):468-477. HE Manchao,MA Zimin,GUO Zhibiao,et al. Key parameters of the gob-side entry retaining formed by roof cutting and pressure release in deep medium-thickness coal seams[J]. Journal of China University of Mining & Technology,2018,47(3):468-477.
[8] 迟宝锁,周开放,何满潮,等. 大采高工作面切顶留巷支护参数优化研究[J]. 煤炭科学技术,2017,45(8):128-133. CHI Baosuo,ZHOU Kaifang,HE Manchao,et al. Optimization research on supporting parameters of roof cutting entry retaining with large mining height face[J]. Coal Science and Technology,2017,45(8):128-133.
[9] WANG Qi,HE Manchao,YANG Jun,et al. Study of a no-pillar mining technique with automatically formed gob-side entry retaining for longwall mining in coal mines[J]. International Journal of Rock Mechanics and Mining Sciences,2018,110:1-8. DOI: 10.1016/j.ijrmms.2018.07.005
[10] HE Manchao,GONG Weili,WANG Jiong,et al. Development of a novel energy-absorbing bolt with extraordinarily large elongation and constant resistance[J]. International Journal of Rock Mechanics and Mining Sciences. 2014,67:29-42.
[11] 唐文胜,张国锋,刘小强,等. 异形断面短臂梁卸压效应模拟分析[J]. 中国高新技术企业,2016(9):162-163. TANG Wensheng,ZHANG Guofeng,LIU Xiaoqiang,et al. Simulation analysis of pressure relief effect of special short arm beam[J]. China High-Tech Enterprises,2016(9):162-163.
[12] 金龙. 厚煤层综放工作面安全高效开采技术及应用[D]. 徐州:中国矿业大学,2020. JIN Long. Technology and application of safe and efficient mining in fully mechanized caving face of thick coal seam[D]. Xuzhou:China University of Mining and Technology,2020.
[13] 顾有富,王慧芹. 切顶卸压沿空留巷技术在土城矿的应用[J]. 能源与节能,2013(1):6-7. GU Youfu,WANG Huiqin. Application of top cutting pressure-relief gob-side entry retaining technology in Tucheng Mine[J]. Energy and Energy Conservation,2013(1):6-7.
[14] 刘衍利,黎卫兵,黄星源. 切顶卸压爆破技术在沿空留巷中的应用[J]. 煤矿安全,2014,45(6):132-135. LIU Yanli,LI Weibing,HUANG Xingyuan. The application of roof cutting and pressure relief blasting technology in gob-side entry retaining[J]. Safety in Coal Mines,2014,45(6):132-135.
[15] 宋润权,谢家鹏. 切顶卸压技术在工作面及沿空巷道维护中的应用[J]. 煤炭科技,2012(3):52-54. SONG Runquan,XIE Jiapeng. Application of roof cutting and pressure unloading technology in the maintenance of working face and along the empty roadway[J]. Coal Science & Technology Magazine,2012(3):52-54.
[16] 蔡洪林,尹贤坤,汤朝均,等. 切顶卸压沿空留巷无煤柱开采技术研究与应用[J]. 矿业安全与环保,2012,39(5):15-18. CAI Honglin,YIN Xiankun,TANG Chaojun,et al. Research and application of coal pillar mining technology for cutting pressure discharge along the empty lane[J]. Mining Safety & Environmental Protection,2012,39(5):15-18.
[17] 吴海波,赵凯. 微震监测技术在复杂地质煤矿开采优化及顶板控制中的应用[J]. 现代矿业,2023,39(9):14-17,22. WU Haibo,ZHAO Kai. Application of microseismic monitoring technology in mining optimization and roof control of complex geological coal mine[J]. Modern Mining,2023,39(9):14-17,22.
[18] 李浩荡,蓝航,杜涛涛,等. 宽沟煤矿坚硬厚层顶板下冲击地压危险时期的微震特征及解危措施[J]. 煤炭学报,2013,38(增刊1):6-11. LI Haodang,LAN Hang,DU Taotao,et al. Micro-seismic characteristic and danger-relief method in rock-burst danger period of mining face under hard and thick roof of Kuangou Mine[J]. Journal of China Coal Society,2013,38(S1):6-11.
[19] 孔令海,李峰,欧阳振华,等. 采动覆岩裂隙分布特征的微震监测研究[J]. 煤炭科学技术,2016,44(1):109-113,143. KONG Linghai,LI Feng,OUYANG Zhenhua,et al. Study on microseismic monitoring and measuring of fracture distribution features in mining overburden strata[J]. Coal Science and Technology,2016,44(1):109-113,143.
[20] 孔令海,王永仁,李少刚. 房柱采空区下回采工作面覆岩运动规律研究[J]. 煤炭科学技术,2015,43(5):26-29. KONG Linghai,WANG Yongren,LI Shaogang. Analysis on overburden strata movement law of coal mining face under goaf of room and pillar mining face[J]. Coal Science and Technology,2015,43(5):26-29.
[21] 李鸿昌. 压力的相似模拟实验[M]. 徐州:中国矿业大学出版社,1988. LI Hongchang. A similar simulation experiment for the pressure[M]. Xuzhou:China University of Mining and Technology Press,1988.