综放工作面区段煤柱合理宽度的数值模拟研究.pdf
第 4 2卷第 8期能 源 与 环 保 V o l 4 2 N o 8 2 0 2 0年8月 C h i n aE n e r g ya n dE n v i r o n m e n t a l P r o t e c t i o nA u g . 2 0 2 0 收稿日期 2 0 2 0- 0 5- 1 8 ; 责任编辑 郭海霞 D O I 1 0 . 1 9 3 8 9 / j . c n k i . 1 0 0 3- 0 5 0 6 . 2 0 2 0 . 0 8 . 0 5 6 作者简介 王雨豪( 1 9 9 1 ) , 男, 山西长治人, 助理工程师, 2 0 1 5年毕业于中国矿业大学, 主要从事煤矿生产管理工作。 引用格式 王雨豪, 居昌波, 王亮亮. 综放工作面区段煤柱合理宽度的数值模拟研究[ J ] . 能源与环保, 2 0 2 0 , 4 2 ( 8 ) 2 5 3 2 5 8 . Wa n g Y u h a o , J uC h a n g b o , Wa n g L i a n g l i a n g . N u m e r i c a l s i m u l a t i o ns t u d y o nr e a s o n a b l e w i d t ho f c o a l p i l l a r i ns e c t i o no f f u l l y m e c h a n i z e dc a v i n g f a c e [ J ] . C h i n aE n e r g ya n dE n v i r o n m e n t a l P r o t e c t i o n , 2 0 2 0 , 4 2 ( 8 ) 2 5 3 2 5 8 . 综放工作面区段煤柱合理宽度的数值模拟研究 王雨豪1, 居昌波2, 王亮亮3 ( 1 . 山西潞安集团司马煤业有限公司, 山西 长治 0 4 6 0 0 0 ; 2 . 山东李楼煤业有限公司, 山东 菏泽 2 7 4 7 0 0 ; 3 . 华北科技学院安全工程学院( 中心) , 河北 三河 0 6 5 2 0 1 ) 摘要 为了有效测定厚煤层工作面的煤柱合理宽度, 针对1 2 0 8 工作面生产地质条件, 采用理论计算和 数值模拟相结合的方法, 对合理的煤柱宽度进行了研究。通过理论计算, 初步确定了区段煤柱留设 2 5 m , 运用数值模拟分析了区段煤柱内应力分布规律, 在煤柱宽度在 1 5~ 2 2m时, 应力值处于“ 谷底” , 且处于弹性区。现场应用表明, 巷道顶底板移近量最大为1 6 9m m , 两帮移近量最大为3 8 3m m , 巷道维 护效果好, 煤柱稳定。研究成果可为该矿后续的工作面区段煤柱尺寸留设提供借鉴。 关键词 综放工作面; 区段煤柱宽度; 厚煤层; 数值模拟 中图分类号 T D 3 5 3 ; T D 8 2 2 文献标志码 A 文章编号 1 0 0 3- 0 5 0 6 ( 2 0 2 0 ) 0 8- 0 2 5 3- 0 6 N u me r i c a l s i mu l a t i o ns t u d yo nr e a s o n a b l ew i d t ho f c o a l p i l l a ri ns e c t i o no f f u l l y me c h a n i z e dc a v i n gf a c e Wa n gY u h a o 1, J uC h a n g b o2, Wa n gL i a n g l i a n g3 ( 1 . S h a n x i L u ′ a nG r o u pS i m aC o a l I n d u s t r yC o . , L t d . , C h a n g z h i 0 4 6 0 0 0 , C h i n a ; 2 . S h a n d o n gL i l o uC o a l I n d u s t r yC o . , L t d . , H e z e 2 7 4 7 0 0 , C h i n a ; 3 . S c h o o l o f S a f e t yE n g i n e e r i n g , N o r t hC h i n aU n i v e r s i t yo f S c i e n c e a n dT e c h n o l o g y ( C e n t e r ) , S a n h e 0 6 5 2 0 1 , C h i n a ) A b s t r a c t I no r d e r t or e a s o n a b l yd e t e r m i n et h e r e a s o n a b l e w i d t ho f c o a l p i l l a r i ns e c t i o no f t h i c kc o a l s e a mw o r k i n g f a c e , t h e r e a s o n a b l e w i d t ho f c o a l p i l l a r i ns e c t i o nw a s a n a l y z e db yt h e o r e t i c a l c a l c u l a t i o na n dn u m e r i c a l s i m u l a t i o na c c o r d i n gt ot h ep r o d u c t i o ng e o l o g i c a l c o n d i t i o n s o f 1 2 0 8w o r k i n g f a c e . T h r o u g ht h e o r e t i c a l c a l c u l a t i o n , 2 5mw a s p r e l i m i n a r i l y d e t e r m i n e df o r t h e c o a l p i l l a r i nt h e s e c t i o n , a n d t h es t r e s s d i s t r i b u t i o nl a wi nt h e c o a l p i l l a r w a s a n a l y z e db y n u m e r i c a l s i m u l a t i o n . T h e s t r e s s v a l u e i nt h e r a n g e o f 1 5~ 2 2mi nt h e c o a l p i l l a r i s a t t h e " b o t t o m "a n di nt h e e l a s t i c z o n e . T h e f i e l da p p l i c a t i o ns h o w e dt h a t , t h e m a x i m u md i s p l a c e m e n t o f r o a d w a y r o o f a n df l o o r i s 1 6 9m m , a n dt h em a x i m u md i s p l a c e m e n t o nb o t hs i d e s i s 3 8 3m m , a n dt h e r o a d w a y m a i n t e n a n c e e f f e c t i s g o o da n dt h e c o a l p i l l a r i s s t a b l e . T h er e s e a r c hr e s u l t s c a np r o v i d ear e a s o n a b l er e f e r e n c ef o r t h ec o a l p i l l a r s i z ed e s i g no f t h es u b s e q u e n t w o r k i n gf a c es e c t i o no f t h em i n e . K e y w o r d s f u l l y m e c h a n i z e dc a v i n gf a c e ; s e c t i o n a l p i l l a r w i d t h ; t h i c kc o a l s e a m ; n u m e r i c a l s i m u l a t i o n 0 引言 区段煤柱是保留在 2个相邻长壁工作面之间的 保护性煤柱, 其主要功能是隔离采空区并保持巷道 的稳定性[ 1 5 ]。煤柱宽度过小, 工作面回采过程中, 不能有效支撑上区段侧向应力及下区段前支撑压 力, 导致煤柱塑性变形区大于煤柱宽度, 煤柱被压裂 破坏[ 6 ]; 煤柱宽度过大, 处于上区段应力升高区范 围, 难以维护, 并造成资源浪费。近年来, 国内外学 者通过理论计算分析、 数值模拟、 现场实测等方法确 定煤柱宽度, 取得了许多研究成果[ 7 1 0 ]。运用计算 机数值模拟和理论计算, 王爱文等[ 1 1 ]分析了煤柱的 稳定性和支护应力分布, 确定了三软倾斜厚煤层综 放工作面煤柱宽度。王建辉等[ 1 2 ]建立了计算煤柱 两侧塑性破坏带的理论公式, 结合现场实测资料, 提 出了不同的煤柱宽度, 并进行了数值模拟分析, 确定 了高瓦斯综放工作面煤柱的合理宽度。蒋昊良 等[ 1 3 1 8 ]采用三维有限差分法模拟分析了工作面推 352 2 0 2 0年第 8期 能 源 与 环 保第 4 2卷 进 1 5 0m时工作面后不同距离开采对侧向围岩的 破坏特征和应力分布, 研究了综采顶煤煤巷煤柱宽 度。 司马煤业3 号煤层平均采高 6 3m , 回采过程中 矿压显现强烈, 以往的设计中工作面的区段煤柱一般 设置在3 0 ~ 3 5m , 煤柱的承载条件和合理的设计宽 度, 需要有效的数据支撑。因此, 研究煤柱的合理宽 度对厚煤层放顶煤工作面的优化具有重要意义[ 1 9 ]。 1 工程概况 1 2 0 8工作面平均开采高度为 6 3m , 开采下二 叠统山西组下部的 3号煤层。工作面走向长度为 2 0 0m , 工作面东侧倾斜长度为12 4 7 5m , 西侧倾斜 长度为 12 4 7 0m , 可采区面积为 2 4 94 5 0m 2。1 2 0 8 工作面的东部为已完成回采的 1 2 0 7工作面, 西部为 未采区, 区间煤柱留设为 3 5m , 在矿区的南部, 南部 有保护煤矿边界的煤柱, 而在北部, 第 2矿区有 3条 主要巷道。1 2 0 8工作面布置如图 1所示, 地层柱状 如图 2所示。 图 1 1 2 0 8工作面布置 F i g 1 L a y o u t o f w o r k i n gf a c e 1 2 0 8 2 区段煤柱合理宽度理论分析 煤柱宽度计算理论模型如图 3所示。煤柱宽度 计算公式为 B= X 1+ X2+ X3 ( 1 ) 式中, X 1为上区段工作面开采后在采空区侧煤体中 产生的塑性区宽度, 根据式( 2 ) 计算。 X 1= M A 2 t a nφ 0 l n k γ H+ C 0 t a nφ 0 C 0 t a nφ 0+ P z A ( 2 ) 式中, m为巷道高度, 取 3m ; A为侧压系数, 取 图 2 煤层顶底板地层综合柱状 F i g 2 C o mp r e h e n s i v ec o l u mnd i a g r a mo f c o a l s e a mr o o f a n df l o o rs t r a t a 图 3 理论煤柱宽度计算模型 F i g 3 T h e o r e t i c a l p i l l a rw i d t hc a l c u l a t i o nmo d e l 0 8 ; φ 0为界面的内摩擦角, 取 3 0 5 ; C0为黏聚力, 取 1 5M P a ; k 为应力集中系数, 取 3 ; γ为平均容重, 取 2 5k N / m 3 ; H为埋深, 4 1 2 8 5 1~ 4 2 4 2 7 8m , 计算中取 4 1 5m ; P z为阻力, 取 0 0 8M P a ; X2为帮锚杆有效长 度, 现有帮锚杆长 2m ; X 3为考虑到由于煤层厚度大 而煤柱宽度增加, 取 X 1+ X2值的1 5 % ~ 3 5 %计算。 将数据代入公式, 可得到理论煤柱宽度值, 计算 出煤柱宽度值为 2 4 2 5~ 2 8 4 7m 。从资源回收角 度考虑, 初步考虑区段煤柱留设 2 5m 。 3 数值模拟分析 3 1 模型的建立 根据司马矿的生产地质条件, 建立了 F L A C 3 D数 值计算模型。该模型包括 2 8层 3号煤及其上下地 层, 总厚度为1 1 0m , 宽2 4 5m ( 结合现场及所建模型 452 2 0 2 0年第 8期王雨豪, 等 综放工作面区段煤柱合理宽度的数值模拟研究 第 4 2卷 的复杂性, 工作面长度取实际长度一半, 则 1 2 0 8和 1 2 0 7工作面的长为 1 0 0m , 1 2 0 7运输巷和 1 2 0 8回 风巷的截面尺寸为 5m 3m , 煤柱宽度为 2 5m ) 。 其中, 1 2 0 8工作面推进长度 1 2 0m , 开挖应逐步进 行, 每前进 2 0m应计算10 0 0步, 当采矿达到1 0 0m ( 即模型尺寸为 2 4 5m 1 2 0m 1 1 0m ) 时, 计算应 保持平衡。为了方便巷道应力应变的数值计算和分 析条件, 将均匀宽度网格布置在水平方向上, 煤柱巷 道附近的网格较密于工作面, 均匀分布的网格也分 布在垂直方向。3号煤层网格划分相对较密, 上、 下 地层网格划分相对较稀疏。 将水平位移约束加在模型两侧, 垂直位移约束 施加在前后边界, 在上边界上施加了均匀载荷。根 据上覆地层的重力计算均匀荷载的大小。司马煤矿 上覆层地层密度为 25 0 0k g / m 3, 3号煤层埋深为 3 9 4m 。根据模型层的分布, 施加在模型上部的均匀 分布载荷为 8 6 5M P a 。模拟了上覆地层的自重, 侧 压系数为 1 。根据司马煤矿的地质资料和煤参数实 测结果, 模型的煤岩体力学参数见表 1 。 表 1 模型的煤岩力学参数 T a b 1 C o a l a n dr o c kp a r a me t e r s o f t h emo d e l 岩层名称 抗拉强度/ M P a 黏聚力/ M P a 内摩擦角/ ( ゜) 体积模量 K / G P a 剪切模量 G / G P a 砂质泥岩1 5 52 . 0 03 04 1 71 9 2 粉砂岩5 4 63 . 0 03 23 0 0 31 3 8 6 细粒砂岩3 5 99 7 83 52 0 0 21 5 0 1 煤0 7 21 . 0 02 80 8 30 3 8 数值计算模型与初始应力平衡如图 4所示。 图 4 F L A C 3 D数值模拟模型 F i g 4 E s t a b l i s h me n t o f F L A C 3 Dn u me r i c a l s i mu l a t i o nmo d e l 1 2 0 8回风巷道支护情况如图 5所示。 图 5 1 2 0 8巷道支护情况 F i g 5 S u p p o r t i n gc o n d i t i o n s o f 1 2 0 8r o a d w a y 3 2 模拟方案及结果分析 对1 2 0 8工作面分别按推进0 、 2 0 、 4 0 、 6 0 、 8 0 、 1 0 0 m煤柱内部应力分布进行了模拟计算。1 2 0 8工作 面回采后煤柱垂直应力分云图如图 6所示, 1 2 0 8工 作面推进中垂直应力分布曲线如图 7所示。 图 6 1 2 0 8工作面回采后煤柱垂直应力分云图 F i g 6 V e r t i c a l s t r e s s c l o u dd i a g r a mo f p i l l a ra f t e rmi n i n go n1 2 0 8w o r k i n gf a c e 552 2 0 2 0年第 8期 能 源 与 环 保第 4 2卷 图 7 1 2 0 8工作面推进中垂直应力分布曲线 F i g 7 V e r t i c a l s t r e s s d i s t r i b u t i o nc u r v eo f 1 2 0 8 w o r k i n gf a c ed u r i n gp r o p u l s i o n 由垂直应力分布曲线可得 煤柱宽度为 2 5m 时, 应力分布仍呈现“ 双峰一谷” , “ 双峰” 位置分别 出现在面间煤柱内 4~ 2 4m处; “ 谷” 在 1 5~ 2 2m 分布。此宽度范围内, 煤柱的支承压力处于平缓水 平, 并且是 2个开采扰动的叠加, 可以认为是煤柱内 的一定宽度的弹性区, 具有较强的承载力[ 2 0 ]; 在煤 柱的其他区域, 煤体具有完全的塑性屈服力, 仅具有 一定的承载能力。依据煤柱留设原理, 当弹性区的 宽度大于或等于煤层开采高度的 2倍时, 煤柱的稳 定性较好。 综合以上分析, 得出煤柱中间有一个弹性岩心, 宽度为 7m , 仍然有一定的剩余量, 所以合理的留设 煤柱宽度为 2 5m 。 4 现场应用效果分析 司马煤矿在 1 2 0 8工作面留设了宽 2 5m的区段 煤柱, 监测分析巷道顶底板和两帮表面移近量可以 判定巷道围岩的稳定性和锚杆支护效果。表面位移 测量站设置在距工作面 1 0 0m的巷道断面, 如图 8 所示。 图 8 表面位移测量站布置 F i g 8 L a y o u t o f s u r f a c ed i s p l a c e me n t me a s u r e me n t s t a t i o n 顶板和两帮位移的监测, 在 1 2 0 8工作面回风巷 棚架号为 1 0 6 5处的断面进行( 图 8 ) 。得到实测数 据后, 根据结果绘制巷道顶底板及两帮变形曲线, 如 图 9 图 1 2所示。 图 9 顶底板变形量曲线 F i g 9 C u r v eo f r o o f a n df l o o rd e f o r ma t i o n 图 1 0 顶底板变形速度曲线 F i g 1 0 D e f o r ma t i o ns p e e do f r o o f a n df l o o r 图 1 1 两帮的变形量曲线 F i g 1 1 D e f o r ma t i o nc u r v e s o f t h et w os i d e s 652 2 0 2 0年第 8期王雨豪, 等 综放工作面区段煤柱合理宽度的数值模拟研究 第 4 2卷 图 1 2 两帮的变形速度曲线 F i g 1 2 D e f o r ma t i o nv e l o c i t yc u r v e s o f t h et w os i d e s 由图 9 图 1 2可知, 监测 2 1d , 顶底板最大累 计移近量为 1 6 9m m , 两测最大累计移近量为 3 8 3 m m ; 监测 1 6d 内, 巷道围岩变形趋于基本稳定。顶 板和底板的最大移动速率为 1 1 1m m/ d , 两帮的最大 移动速率达到 2 6 6m m/ d , 均发生在测量站布置后第 2 1天。初期, 两侧变形速率为 3 9m m/ d , 顶底板的 变形速率为 1 7m m/ d , 后期顶底板和两帮移动的速 度较大。监测站围岩整体位移不大。分析可得, 两 帮和顶底板移近量均在巷道允许变形范围内, 快速 变形速率时间短, 巷道留设效果较好。 5 结论 ( 1 ) 通过理论分析, 得出合理的煤柱留设宽度 为 2 5m 。 ( 2 ) 针对 2 5m的区段煤柱宽度进行了数值模 拟验证, 数值计算表明, 在煤柱 1 5~ 2 2m应力值处 于“ 谷底” , 且处于弹性区, 证明留设宽 2 5m的煤柱 是合理的。 ( 3 ) 通过现场工业性应用, 巷道两帮最大移近 量为 3 8 3m m , 顶底板最大移近量为 1 6 9m m , 巷道维 护效果好, 煤柱稳定。 参考文献( R e f e r e n c e s ) [ 1 ] 柏建彪. 沿空掘巷围岩控制[ M] . 徐州 中国矿业大学出版社, 2 0 0 6 . [ 2 ] 马立强, 张东升, 王红胜, 等. 厚煤层巷内预置充填带无煤柱开 采技术[ J ] . 岩石力学与工程学报, 2 0 1 0 , 2 9 ( 4 ) 6 7 4 6 8 0 . M aL i q i a n g , Z h a n gD o n g s h e n g , Wa n gH o n g s h e n g , e t a l . N o n p i l l a r m i n i n g t e c h n o l o g y o f p r e f i l l e df i l l i n g b e l t i nt h i c kc o a l s e a mr o a d w a y [ J ] . C h i n e s eJ o u r n a lo fR o c kM e c h a n i c sa n dE n g i n e e r i n g , 2 0 1 0 , 2 9 ( 4 ) 6 7 4 6 8 0 . [ 3 ] 董文敏. 高瓦斯矿大采高工作面区段煤柱尺寸合理确定 [ J ] . 煤炭科学技术, 2 0 0 7 , 3 5 ( 1 2 ) 6 7 7 0 . D o n g We n m i n . P i l l a r s i z eo f l a r g em i n i n gf a c es e c t i o ni nh i g hg a s m i n ei s r e a s o n a b l yd e t e r m i n e d [ J ] . C o a l S c i e n c ea n dT e c h n o l o g y , 2 0 0 7 , 3 5 ( 1 2 ) 6 7 7 0 . [ 4 ] 秦永洋, 许少东, 杨张杰. 深井沿空掘巷煤柱合理宽度确定 及 支护参数优化[ J ] . 煤炭科学技术, 2 0 1 0 , 3 8 ( 2 ) 1 5 1 8 . Q i nY o n g y a n g , X uS h a o d o n g , Y a n gZ h a n g j i e . D e t e r m i n a t i o no f r e a s o n a b l e w i d t ho f c o a l p i l l a r a n do p t i m i z a t i o no f s u p p o r t p a r a m e t e r s a l o n gt h er o a d w a yi nd e e pw e l l [ J ] . C o a l S c i e n c e a n dT e c h n o l o g y , 2 0 1 0 , 3 8 ( 2 ) 1 5 1 8 . [ 5 ] 申志平. 沿空掘巷煤柱宽度与支护参数的研究与应用[ J ] . 煤 炭科学技术, 2 0 1 0 , 3 8 ( 7 ) 3 1 3 4 . S h e nZ h i p i n g . R e s e a r c ha n da p p l i c a t i o no f c o a l p i l l a rw i d t ha n d s u p p o r t p a r a m e t e r s i nr o a d w a y e x c a v a t i o na l o n gt h ev o i d [ J ] . C o a l S c i e n c ea n dT e c h n o l o g y , 2 0 1 0 , 3 8 ( 7 ) 3 1 3 4 . [ 6 ] 陈科. 沿空掘巷小煤柱破坏规律及合理宽度的确定[ J ] . 煤矿 安全, 2 0 0 9 ( 2 ) 1 0 0 1 0 2 . C h e nK e . D e t e r m i n a t i o no ft h ed e s t r u c t i o nl a w a n dr e a s o n a b l e w i d t ho f s m a l l c o a l p i l l a r sd u ga l o n gt h er o a d w a y [ J ] . S a f e t yi n C o a l M i n e s , 2 0 0 9 ( 2 ) 1 0 0 1 0 2 . [ 7 ] 张兴润. 综放工作面采空区侧煤柱合理宽度留设研究[ J ] . 能 源技术与管理, 2 0 1 9 , 4 4 ( 3 ) 9 6 9 8 . Z h a n gX i n g r u n . S t u d yo nr e a s o n a b l ew i d t ho f c o a l p i l l a ri ng o a f s i d e o f f u l l y m e c h a n i z e dc a v i n gf a c e [ J ] . E n e r g yT e c h n o l o g ya n d M a n a g e m e n t , 2 0 1 9 , 4 4 ( 3 ) 9 6 9 8 . [ 8 ] 梁翱翔. 司马煤业二采区综放工作面区段煤柱尺寸优化研究 [ D ] . 徐州 中国矿业大学, 2 0 1 9 . [ 9 ] 吕晓波. 采动应力影响条件下护巷煤柱留设尺寸研究[ J ] . 煤 矿安全, 2 0 1 9 , 5 0 ( 1 2 ) 1 8 9 1 9 1 , 1 9 6 . L X i a o b o . S t u d yo nt h es i z eo f r e t a i n i n gp i l l a r i nr o a d w a yp r o t e c t i o nu n d e r m i n i n g s t r e s s [ J ] . S a f e t yi nC o a l M i n e s , 2 0 1 9 , 5 0 ( 1 2 ) 1 8 9 1 9 1 , 1 9 6 . [ 1 0 ] 宋虎森, 李金华, 李晓鹏, 等. 司马煤矿 1 2 0 8风巷回采期间矿 压监测研究[ J ] . 煤炭与化工, 2 0 1 8 , 4 1 ( 4 ) 3 0 3 2 , 3 6 . S o n gH u s e n , L iJ i n h u a , L iX i a o p e n g , e ta l . M o n i t o r i n go fm i n e p r e s s u r e d u r i n gs t o p i n g i n 1 2 0 8w i n dl a n e o f S i m a C o a l M i n e [ J ] . C o a l a n dC h e m i c a l I n d u s t r y , 2 0 1 8 , 4 1 ( 4 ) 3 0 3 2 , 3 6 . [ 1 1 ] 王爱文, 潘一山, 李忠华. 深部孤岛工作面区段煤柱冲击地压 的试验研究[ J ] . 安全与环境学报, 2 0 1 7 , 1 7 ( 3 ) 9 0 7 9 1 1 . Wa n gA i w e n , P a nY i s h a n , L iZ h o n g h u a . E x p e r i m e n t a ls t u d yo n c o a l p i l l a r r o c kb u r s t i nd e e pi s l a n dw o r k i n g f a c e s e c t i o n [ J ] . J o u r n a l o f S a f e t ya n dE n v i r o n m e n t , 2 0 1 7 , 1 7 ( 3 ) 9 0 7 9 1 1 . [ 1 2 ] 王建辉, 武谋达. 文家坡煤矿区段煤柱合理留设宽度研究[ J ] . 煤矿开采, 2 0 1 8 , 2 3 ( 4 ) 2 9 3 2 . Wa n g J i a n h u i , WuM o u d a . S t u d y o f r a t i o n a l l a y o u t w i d t h o f s e c t i o n a l c o a l p i l l a r o f We n j i a p o C o a l M i n e [ J ] . C o a l M i n i n g T e c h n o l o g y , 2 0 1 8 , 2 3 ( 4 ) 2 9 3 2 . [ 1 3 ] 蒋昊良. 沿空掘进小煤柱合理宽度留设及支护分析[ J ] . 内蒙 古煤炭经济, 2 0 1 9 ( 1 0 ) 1 4 4 1 4 5 . J i a n g H a o l i a n g . A n a l y s i s o nt h e r e a s o n a b l e w i d t ho f s m a l l c o a l p i l 752 2 0 2 0年第 8期 能 源 与 环 保第 4 2卷 l a r s f o r r e t a i n i n ga n ds u p p o r t i n ga l o n gt h ee x c a v a t i o n [ J ] . I n n e r M o n g o l i aC o a l E c o n o m y , 2 0 1 9 ( 1 0 ) 1 4 4 1 4 5 . [ 1 4 ] 王方田, 屠世浩, 李召鑫, 等. 浅埋煤层房式开采遗留煤柱突变 失稳机理研究[ J ] . 采矿与安全工程学报, 2 0 1 2 , 2 9 ( 6 ) 7 7 0 7 7 5 . Wa n g F a n g t i a n , T uS h i h a o , L i Z h a o x i n , e t a l . S t u d yo nm e c h a n i s m o f a b r u p t i n s t a b i l i t yo f c o a l p i l l a r l e f t o v e r f r o mr o o m t y p em i n i n g i ns h a l l o ws e a m[ J ] . J o u r n a lo fM i n i n g& S a f e t yE n g i n e e r i n g , 2 0 1 2 , 2 9 ( 6 ) 7 7 0 7 7 5 . [ 1 5 ] 赵维生, 梁维, 金志远, 等. 遗留煤柱下伏近距煤层巷道布置及 围岩稳定性[ J ] . 煤矿安