★ 科技与工程 ★

大采高综放小煤柱沿空掘巷围岩稳定性数值模拟研究

彭林军1 周光华2 宋振骐1 王金平3 何维胜3

(1.大连大学院士创业园，辽宁省大连市，116622；2.国能宁夏煤业集团有限责任公司，宁夏自治区银川市，750011；3.国能宁夏煤业集团有限责任公司羊场湾煤矿，宁夏自治区银川市，750011)

摘 要 以羊场湾煤矿130205超大采高综放工作面为背景，以实用矿山压力理论为指导，运用FLAC3D数值模拟方法对该工作面回风巷道6 m小煤柱沿空掘巷围岩应力分布进行研究。研究表明：小煤柱风巷岩体塑性区范围为0～10 m，高应力区范围为10～14 m；工作面回采期间巷道两帮最大移近量为600 mm、顶底板最大移近量为430 mm。回采期间6m小煤柱中心处约有1.5 m稳定区。通过围岩应力分布规律研究，证明巷道变形量在可控范围满足生产要求，为大采高综放小煤柱沿空掘巷围岩稳定性提供了科学依据。

关键词 大采高综放开采 小煤柱 沿空掘巷 数值模拟 支承压力分析

中图分类号 TD353

文献标识码 A

基金项目：国家自然科学基金资助项目(51674058)，辽宁省自然科学基金资助项目(201602026)

引用格式：彭林军，周光华，宋振骐等. 大采高综放小煤柱沿空掘巷围岩稳定性数值模拟研究[J]. 中国煤炭，2019,45(12)：94-99.

Peng Linjun, Zhou Guanghua, Song Zhenqi, et al. Numerical simulation research on surrounding rock stability of gob-side entry driving in small pillar in fully mechanized caving face with large mining height[J]. China Coal，2019,45(12)：94-99.

Numerical simulation research on surrounding rock stability of gob-side entry driving in small pillar in fully mechanized caving face with large mining height

Peng Linjun1, Zhou Guanghua2, Song Zhenqi1, Wang Jinping3, He Weisheng3

(1. Academician Pioneering Park, Dalian University, Dalian, Liaoning 116622, China; 2. China Energy Ningxia Coal Industry (Group) Co., Ltd., Yinchuan, Ningxia 750011, China; 3. Yangchangwan Coal Mine, China Energy Ningxia Coal Industry (Group) Co., Ltd., Yinchuan, Ningxia 750011, China)

Abstract This paper took 130205 fully mechanized caving face with extremely large mining height in Yangchangwan Coal Mine as the background and the practical mine pressure theory as the guidance, studied the stress distribution of the surrounding rock of gob-side entry driving in 6 m width small coal pillar in the return air roadway of the working face by FLAC 3D numerical simulation method. The results showed that the plastic zone in small pillar in air roadway was 0～10 m, and the high stress zone was 10～14 m; the maximum displacement of two sides of roadway was 600 mm, and the maximum displacement of roof and floor was 430 mm. There was about 1.5 m stable area in the center of 6 m width small coal pillar during mining. Through the study of the stress distribution of the surrounding rock, it was proved that the deformation of the roadway could meet the production requirements within the controllable range, which provided a scientific basis for the stability of the surrounding rock of the gob-side entry driving in the small coal pillar in fully mechanized caving face with large mining height.

Key words fully mechanized caving mining with large mining height, small coal pillar, gob-side entry driving, numerical simulation, support pressure analysis

近年来,无(小)煤柱巷道在我国取得了较大的发展[1-3]。随着煤矿开采强度和深度的增加，沿空掘巷小煤柱巷道受二次采动动压影响更为剧烈。因此，研究在不同采动条件下，掌握以“岩层运动和支承压力分布”为核心[4-11]，建立其动态结构力学模型,实现由大煤柱护巷的传统开采系统向无(小)煤柱护巷开采体系的变革,是煤矿开采围岩破裂灾变机理及其控制的关键。

我国学者针对采动巷道围岩控制问题，进行了大量的相关研究工作[12-15]，丰富了采动巷道围岩控制理论,解决了大量的小煤柱巷道围岩控制问题；但对于大采高综放开采煤柱尺寸对采动巷道的影响程度以及小煤柱异型巷道产生非均匀大变形的力学本质仍需进一步研究，要为小煤柱巷道找到更加合理的支护方案[16-18]。本文针对羊场湾煤矿小煤柱巷道非均匀大变形控制难题,研究沿空留设6 m小煤柱时，回风巷道围岩塑性区分布形态与应力分布规律以及采场覆岩结构运动特征，揭示了小煤柱巷道非均匀变形的机理，提出具有针对性的控制对策并进行井下试验，对支护效果进行评价。

1 矿井工作面概况

130205工作面位于羊场湾一号井田东部，地形低缓平坦，起伏不大，地表为沙丘覆盖，开采2#煤层，工作面宽度350 m，埋深650 m，煤厚8.2～10.7 m，平均煤厚8.4 m。2#煤层伪顶为炭质泥岩，直接顶为粉砂岩、细砂岩，其上部为中砂岩、细砂岩。

130205风巷开口位于13采区回风下山，机巷开口位于13采区胶运下山，按设计方位349°施工。130205工作面上邻130203工作面(截至2015年1月19日，工作面已回采完毕)，130205风巷与130203机巷原留设35 m保安煤柱；工作面北以F201正断层为界，南侧以13采区井筒保安煤柱为界，130205工作面下部为原始煤层未进行采动，无采掘活动影响工作面掘进。距离130205工作面最近钻孔柱状图为1916#钻孔，如图1所示；1916#钻孔柱状岩石力学参数如表1所示。

图1 130205工作面(1916#)钻孔柱状图

表1 130205工作面(1916#)钻孔岩石力学参数

编号岩石名称厚度/m深度/m抗压强度/MPa抗拉强度/MPa弹性模量/MPa泊松比粘聚力/MPa容重/(kN·m-3)M10粉砂岩14.68474.90------M9中砂岩22.40497.30------M8粉砂岩35.29532.59------M7中砂岩3.70536.29------M6粉砂岩7.59543.88------M5粗砂岩39.27583.15------M41#煤层1.15584.30------M3粉砂岩5.24589.5445.24.5214.60.234.6623.5M2中砂岩12.00601.5438.63.8614.40.234.8323.5M1粉砂岩8.56610.1035.73.5712.50.224.2722.0M02#煤层8.70618.8013.01.301.860.250.7515.0F1粉砂岩1.45620.45------

2 小煤柱沿空掘巷计算模型建立

数值模型中，各层岩体力学参数依据表1取值。模型的前后、左右边界施加水平约束，模型的底边界施加水平及垂直约束，根据模型的几何尺寸划分计算网格，考虑到模型单元总数的限制和对小煤柱沿空巷道围岩应力分析进行区域集中性研究，建立数值模拟模型如图2所示。工作面与风巷数值计算模型如图3所示。

图2 三维数值模型网格图

图3 工作面与风巷数值计算模型

3 风巷围岩稳定性数值模拟分析

3.1 风巷掘进期间稳定性数值模拟分析

由地层变形云图和岩体位移矢量图可以看出，风巷掘进期间，风巷周围岩体变形逐渐增大，尤其是拱顶附近变形值最大，风巷掘进至距开切眼350 m时，风巷附近变形最大值为114.4 mm。拱顶往上延伸，地层变形收敛幅度很快，向岩体内延伸5 m 后，围岩变形基本在10 mm以内。底部最大隆起量约45 mm，由于断面异形，两侧边墙呈现明显的变形不对称分布。掘进期间风巷小煤柱帮最大变形为45 mm，实体煤帮最大变形为158 mm；顶板下沉量为114 mm，底鼓量为75 mm。岩体塑性区分布图见图4。

图4 风巷掘进350 m岩体塑性区分布图

风巷掘进期间岩体垂直应力SZZ(表示岩体最大垂直应力)和水平应力SXX(表示岩体最大水平应力)分布如图5所示。

由图5可知，岩体垂直应力SZZ分布：风巷两帮应力大于顶部和底部，高应力区最大应力33.13 MPa，高应力区边缘距离风巷7 m左右，高应力区核心距离风巷8.5 m左右；岩体水平应力SXX分布， 6 m煤柱在130203工作面帮的顶部SXX应力较大，最大值为19.76 MPa。在现有支护条件下，围岩未出现明显较大的拉应力，风巷整体稳定性较好。

3.2 回采期间围岩稳定性数值模拟分析

130205工作面回采期间三维数值示意图见图6。

3.2.1 回采期间工作面超前支承压力分析

对130205工作面回采期间工作面超前支承压力分析如图7所示。

图5 风巷掘进期间岩体应力分布图

图6 工作面数值模型开挖图

回采期间工作面超前煤壁0～4 m为塑性区、低应力区，4～10 m为塑性区、高应力区，10～14 m为应力峰值区，14～50 m为应力高值区、缓降区。回采期间工作面超前支承压力分布如图8所示。

3.2.2 回采期间工作面处围岩应力分布图

回采期间工作面处岩体应力分布见图9。数值计算结果表明，6 m小煤柱中心处SZZ最大，应力值达8.05 MPa；岩体应力SXX主要在6 m小煤柱上部，尤其是靠近130203开采区上部，最大值为23.55 MPa。另外，应力计算表明，现有支护条件下，风巷受二次动压影响，围岩并未出现较大的拉应力区(压应力为负值，拉应力为正值)，6 m小煤帮侧主要受压，只是在小煤柱煤帮出现较小范围的拉应力区，其他部位拉应力仍小于1.0 MPa，这说明现有支护条件下，回采期间工作面处6 m小煤柱能够保持整体稳定性。

图7 回采工作面超前支承压力及塑性区分布图

图8 回采期间工作面超前支承压力分布曲线图

3.2.3 回采期间应力峰值处风巷围岩塑性区

回采期间应力峰值处围岩塑性区分布见图10。由图10可知，靠近巷道实体煤帮回采期间塑性区并未明显扩展，片帮深度约0.3 m。在小煤柱侧帮的剪切和拉伸塑性区深度约0.8 m，回采期间应力峰值处风巷变形计算结果也表明小煤柱侧帮变形较大，最大值为560 mm，因此，回采期间应力峰值处小煤柱侧帮鼓帮深度约0.6 m。

回采期间应力峰值围岩变形趋势见图11。从图11可以看出，风巷不仅仅在实体煤帮和6 m煤柱之间存在相对移近变形，在风巷顶部和6 m煤柱之间也发生了相对变形，如图11(a)箭头所示，风巷顶部和6 m煤柱间夹角会变小，如图11(b)蓝色的圈表示，以夹角为中心，越远离夹角位置的岩体变形趋势会越大。

图9 风巷回采期间岩体应力分布图

回采期间应力峰值处岩体变形位移矢量图见图12。由图12可以看出，随着工作面的推进，回采期间应力峰值处风巷周围岩体变形逐渐增大，6 m小煤柱帮变形明显增大，在靠近小煤柱帮深度2.5 m范围以内，风巷顶底板移近量为730 mm，其中顶板下沉量为460 mm，底鼓量为270 mm。两帮移近量为870 mm，其中小煤柱帮为560 mm，实体煤帮为310 mm。回采期间6 m小煤柱支护结构受力会增大，需要重点对6 m小煤柱进行加固。

图10 回采期间应力峰值处岩体塑性区分布图

图11 回采期间应力峰值围岩变形趋势

3.2.4 回采期间应力峰值处风巷稳定性分析

回采期间工作面处煤层开采岩体塑性区分布图与裂隙发展区域示意图见图12，大采高工作面推进过程中，工作面煤壁在支承压力作用下，沿推进方向在一定范围内产生破坏。

因此，根据煤层开采岩体最大应力云图的分布特点以及相对应的煤层塑性区分布图，计算得到的煤壁片帮情况如下：

(1)130205开采面侧6 m煤壁(图中I区)深度0～2.0 m为拉伸和剪切塑性区，为典型的鼓出形式破坏，在支承压力作用下，煤体自身强度的降低导致工作面煤壁外溢，可能造成煤壁的片帮；

(2)130203开采面6 m小煤柱(图中II区)深度0～1.5 m为拉伸和剪切塑性区，在支承压力作用下，煤体自身强度降低导致工作面煤壁外溢，可能造成煤壁的片帮；

(3)风巷实体煤壁(图中III区)深度0～2.5 m和顶部深度范围内，岩体破碎较严重，存在剥落趋势；

(4)6 m小煤柱中心2.0～4.5 m范围内尽管部分煤层进入塑性区，但主要为剪切塑性区，最小主应力为压应力，并不存在片帮，岩体破碎程度不剧烈。

综合分析表明，6 m小煤柱中心处约有1.5 m是非破碎区。

图12 岩体塑性区分布图与裂隙重点发展区域

4 结论

(1)通过数值模拟分析得出，小煤柱风巷岩体塑性区范围为0～10 m，高应力区范围为10～14 m；14～50 m为应力高值区、缓降区。130205沿空巷道区段煤柱留设6 m宽度在塑形区范围内，尺寸合理。

(2)风巷掘进期间顶底板移近量189 mm，两帮移近量203 mm。回采期间工作面峰值处，风巷顶底板移近量730 mm，两帮移近量870 mm。6 m小煤柱中心处约有1.5 m是非破碎区。

(3)区段煤柱应在满足安全生产的前提下，尽量减少煤柱宽度，从而将回采巷道完全布置在应力降低区内，更有利于巷道的维护、煤柱的稳定性，达到安全高效生产。

参考文献：

[1] 何满潮,王炯,孙晓明等. 负泊松比效应锚索的学特性及其在冲击地压防治中的应用研究[J]. 煤炭学报,2014,39(2):214-221.

[2] 侯朝炯, 马念杰. 煤层巷道两帮煤体应力和极限平衡区的探讨[J]. 煤炭学报, 1989, 14(4): 21-29.

[3] 许家林,钱鸣高.关键层运动对覆岩及地表移动影响的研究[J]. 煤炭学报,2000,25(2):122-126.

[4] 宋振骐,蒋宇静,彭林军等. 煤矿重大事故预测和控制的动力信息基础的研究[M]. 北京:煤炭工业出版社,2003:120-146.

[5] Wang Jinhua. Key technology for fully-mechanized top coal caving with large mining height in extra-thick coal seam[ J]. Journal of China Coal Society, 2013,38(12):2089-2098.

[6] 康红普,姜鹏飞,蔡嘉芳. 锚杆支护应力场测试与分析[J]. 煤炭学报，2014，39(8)：1521-1529.

[7] 秦永洋，许少东，杨张杰.深井沿空掘巷煤柱合理宽度确定及支护参数优化[J].煤炭科学技术，2010, 38(2)：20-23.

[8] 李化敏. 沿空留巷顶板岩层控制设计[J]. 岩石力学与工程学报,2000,19(5):651-654.

[9] 姜福兴,王建超,孙广京等. 深部开采沿空巷道冲击危险性的工程判据[J]. 煤炭学报,2015,40(8): 1729-1736.

[10] 何满潮，郭志飚. 恒阻大变形锚杆力学特性及其工程应用[J]. 岩石力学与工程学报，2014，33(7)：1297-1308.

[11] 张东升,马立强,缪协兴等. 综放沿空留巷围岩变形影响因素的分析[J]. 中国矿业大学学报,2006,35(1):1-6.

[12] 冯国瑞,王鲜霞,康立勋. 采场上覆岩层面接触块体结构的力学机理分析[J]. 煤炭学报, 2008,33 (1): 33-37.

[13] 彭林军,张东峰,郭志飚等.特厚煤层小煤柱沿空掘巷数值分析及应用[J].岩土力学,2013，34(12)：3609-3616.

[14] 史红，姜福兴. 综放采场上覆岩厚层坚硬岩层破断规律的分析及应用[J]. 岩土工程学报，2006，28(4)：525-528.

[15] 鞠金峰,许家林,朱卫兵等. 近距离煤层工作面出倾向煤柱动载矿压机理研究[J]. 煤炭学报,2010,35(1):15-20.

[16] 李迎富, 华心祝. 沿空留巷上覆岩层关键块稳定性力学分析及巷旁充填体宽度确定[J]. 岩土力学, 2012,33(4): 1134-1140.

[17] 成云海, 姜福兴, 庞继禄. 特厚煤层综放开采采空区侧向矿压特征及应用[J]. 煤炭学报, 2012, 37(7): 1088-1093.

[18] 肖同强, 柏建彪, 王襄禹等. 深部大断面厚顶煤巷道围岩稳定原理及控制[J]. 岩土力学, 2011, 32(6): 1874-1880.

作者简介：彭林军(1973-)，男，河南灵宝人，博士，副教授，主要研究方向为矿山压力与岩层控制。E-mail:penglinjun969@sina.com。

(责任编辑 郭东芝)