赵固二矿非对称底鼓破坏规律研究.pdf
收稿日期2020 06 03 作者简介毕 鹏1981 - ꎬ男ꎬ山西高平人ꎬ工程师ꎬ从事生产技术及管理工作ꎮ doi10. 3969/ j. issn. 1005 -2798. 2020. 10. 003 赵固二矿非对称底鼓破坏规律研究 毕 鹏1ꎬ魏文胜2 1. 中煤第五建设有限公司第三工程处ꎬ江苏 徐州 221002ꎻ 2. 中国矿业大学北京共伴生能源精准开采北京市重点实验室ꎬ北京 100083 摘 要针对赵固二矿采动巷道底板破坏严重的问题ꎬ以 11060 工作面轨道巷为研究对象ꎬ采用离散元模 拟软件 3DECꎬ结合现场实测、理论分析等方法ꎬ研究 11060 轨道巷在巷道掘进影响稳定阶段和采动阶段下 底鼓破坏形态、应力场及位移场的差异ꎮ 结果表明受采动影响ꎬ赵固二矿深部采动巷道在高地应力作用 下极易发生非对称底鼓ꎬ应力及位移场分布向采空区一侧偏移ꎮ 关键词3DEC 模拟软件ꎻ采动影响ꎻ高应力ꎻ非对称底鼓 中图分类号TD353 文献标识码A 文章编号1005 2798202010 0008 04 Failure Rule of Asymmetric Floor Heave in Zhaogu No. 2 Coal Mine BI Peng1ꎬWEI Wen sheng2 1. No. 3 Engineering Office of China Coal Fifth Construction CompanyꎬXuzhou 221002ꎬChinaꎻ 2. Beijing Key Laboratory of Companion Energy of China Mining University BeijingꎬBeijing 100083ꎬChina Abstract In view of the serious damage to the floor of the mining roadway in Zhaogu No. 2 Coal Mineꎬtaking the track roadway of 11060 working face as the research objectꎬthe discrete element simulation software 3DEC is adoptedꎬand the methods of field measure ̄ ment and theoretical analysis are combinedꎬthe paper studies the difference of floor heave failure formꎬstress field and displacement field of 11060 track roadway in the influence of roadway excavation on stability stage and mining stage. The results show that the deep mining roadway in Zhaogu No. 2 Coal Mine is prone to asymmetric floor heave and the distribution of stress and displacement field devi ̄ ates to the side of mined - out area under high ground stress. Key words3DEC simulation softwareꎻmining influenceꎻhigh stressꎻasymmetrical floor heave 我国深部矿井资源赋存丰富ꎬ尽管近些年煤矿 智能化无人化不断发展ꎬ但巷道仍是工人与设备在 井下重要的工作环境ꎬ因此研究深部巷道的稳定性 对保障矿井安全生产至关重要[1 -3]ꎮ 当载荷超过深 部巷道底板破坏强度时底板就会发生破坏ꎬ加之一 次巷道掘进使巷道周围应力发生变化ꎬ岩石处于破 坏阶段ꎬ在高应力的挤压作用下极易发生底鼓[4]ꎮ 武精科等[5]认为煤岩层赋存环境差、岩性差、强采动 高应力是发生非对称破坏的原因ꎬ提出了围岩结构 分级分区耦合支护关键技术ꎻ马念杰等[6]研究了围 岩偏应力场和塑性区分布规律ꎬ提出了“蝶”形非对 称破坏理论ꎻ陈上元等[7]认为巷道产生非对称变形 主要是采动影响使巷道围岩应力场大小和方向发生 改变导致的ꎮ 虽然众多学者对深部巷道底鼓破坏特 征研究较多ꎬ但多采用连续元 Flac 进行研究ꎬ本文 采用最新离散元 3DEC 数值模拟软件及地质雷达监 测设备ꎬ研究了赵固二矿底板非破坏特征ꎬ确定了底 板破坏深度ꎬ对赵固二矿的安全生产具有重要意义ꎮ 1 工程概况 赵固二矿为焦作矿区的新建矿井ꎬ底板富水性 强ꎬ水压大ꎬ隔水层薄ꎬ威胁着工作面安全生产ꎮ 11060 工作面轨道巷位于 11041 工作面西侧ꎬ11041 工作面主采二1煤层ꎬ煤层平均厚度约 6. 5 mꎬ平均 倾角 3ꎬ埋深 700 m 左右ꎮ 11060 工作面轨道巷受 11041 工作面采动影响ꎬ经常发生巷道大面积变形ꎬ 底鼓严重ꎬ影响矿井的安全生产和工人的生命安全ꎮ 11041 工作面采掘平面如图 1 所示ꎬ工作面顶底板 情况见表 1ꎮ 8 试试验验研研究究 总第 254 期 图 1 11041 工作面采掘平面 表 1 11041 工作面顶底板情况 顶底板名称岩性特征平均厚度/ m岩石名称 基本底 灰色ꎬ成分为石英ꎬ 含少量暗色矿物ꎬ分 选性 较 好ꎬ 泥 质 胶 结ꎮ 7.46中、细粒砂岩 直接底 灰黑色ꎬ含植物化石 及菱铁质结核ꎮ 7.54砂质泥岩 直接顶 灰黑色ꎬ含植物化石 及大量菱铁质结核ꎮ 6.01砂质泥岩、泥岩 伪顶 黑色ꎬ局部炭砂质ꎬ 水平层理ꎬ植物化石 较多ꎮ 0.5泥岩 基本顶 深灰色ꎬ分选磨圆度 中等ꎬ含大量泥质条 带ꎬ泥硅质胶结ꎮ 5.76大占砂岩 2 底板非对称底鼓破坏分析 2. 1 非对称底鼓特征 图 2 为赵固二矿 11060 轨道巷采用防爆相机拍 摄的底鼓图ꎬ图 3 是经过专业处理后的地质雷达扫 描特征图ꎬ图 3 颜色的差异代表底板不同破坏特征ꎮ 对图 2、图 3 的分析表明ꎬ11060 轨道巷在 11041 工 作面回采期间底板围岩变形破坏严重ꎬ需要大面积 挖底翻修作业ꎬ并且肉眼可见靠近采空区一侧的巷 道底板底鼓程度及最大破坏范围明显大于靠近实体 煤的一侧ꎬ赵固二矿深部 11060 采动巷道底鼓破坏 具有非对称性ꎬ底板最大裂隙发育深度为 20. 5 mꎮ 2. 2 顶板对底板破坏影响 深部开采具有强扰动特性ꎬ工作面的强扰动性 必然会通过煤柱对底板产生一定影响[8]ꎬ主要表现 为工作面顶板的破坏触矸以力或者能量的方式传递 到底板ꎬ加大底板破坏深度及破坏宽度ꎬ底板受强扰 动破坏将会形成塑性屈服区、卸荷破坏区、卸荷扩展 区、触矸破坏区及峰前损伤区见图 4 [9]ꎮ 深部开 采时各区的范围远大于浅部开采ꎬ深部开采强扰动 性产生的触矸也会进一步加大底板破坏范围ꎮ 极易 产生非对称破坏ꎮ 根据砌体梁理论、半无限体理论、滑移线理论可 以建立底板破坏力学模型如图 5 所示ꎬ可得出底板 最大破坏深度 Hmax [10 -11]计算公式为 Hmax= Lsinφ 2cos π 4 + φ 2 e π 4 + φ 2tanθ 1 式中φ 为内摩擦角ꎬꎻL 为极限平衡区宽 度ꎬmꎮ 依据现场实际地质情况及矿方资料知赵固二矿 11060 工作面轨道巷底板岩体内摩擦角为 36ꎬ极限 平衡区宽度 18 mꎬ带入公式计算得底板最大破坏深 度为 21 mꎮ 图 2 赵固二矿 11060 轨道巷非对称底鼓特征 图 3 采动阶段下巷道地质雷达扫描特征 图 4 深部底板强扰动破坏分区 图 5 底板破坏力学模型 9 2020 年 10 月 毕 鹏等赵固二矿非对称底鼓破坏规律研究 第 29 卷第 10 期 3 3DEC 数值模型建立及模拟方案 根据 11060 轨道巷顶底板地质情况及电脑硬件 限制建立 3DEC 数值模拟模型ꎬ模型尺寸为长 宽 高 =334. 8 m 2 m 150 mꎬ模型岩层本构采用 Mohr - Coulomb 模型cons = 1ꎬ节理裂隙采用库伦 滑移模型jcons = 2ꎬ上覆岩层载荷 16. 4 MPaꎮ 模 型的边界条件如图 6ꎬ模型的岩层和力学参数如表 2 所示ꎮ 数值模拟方案模型建立平衡后先开挖 11060 轨道巷ꎬ平衡后保存掘进影响稳定阶段 sav 文件ꎬ再 开挖 11041 工作面ꎬ充分开采后保存采动阶段巷道 底板 sav 文件ꎬ进行对比分析研究ꎮ 4 3DEC 数值模拟分析 4. 1 不同阶段下应力场对比分析 运用离散元数值模拟软件 3DEC 导出巷道在掘 进影响稳定阶段和采动影响阶段下垂直应力和偏应 力云图ꎬ如图 7 所示ꎮ 图 6 模型边界条件 表 2 煤岩体力学参数 岩性厚度/ m密度/ kgm-3 体积模量/ GPa剪切模量/ GPa粘聚力/ MPa抗拉强度/ MPa内摩擦角/ 砂质泥岩8.52 50010. 444.545.362.636 泥岩8.62 3508.825.055.241.532 砂质泥岩6.72 50010. 444.545.362.636 细粒砂岩2.02 6007.606.50134.535 中粒砂岩6.82 7008.675.678.501.735 粗粒砂岩2.62 8009.977.978.201.635 细粒砂岩3.42 6007.606.50134.535 泥岩5.52 3508.825.055.241.532 细粒砂岩5.02 6007.606.50134.535 泥岩11.02 3508.825.055.241.532 香炭砂岩5.82 7005.675.678.501.735 砂质泥岩8.12 50010. 444.545.362.636 图 7 不同阶段巷道应力场对比 由图 7 可知ꎬ深部巷道在掘进影响稳定阶段和 采动阶段围岩应力场存在巨大差异ꎬ掘进影响稳定 阶段的巷道垂直应力、偏应力云图均呈现明显的以 巷道中心为中心的对称性ꎬ采动影响后ꎬ各应力场普 遍向煤柱一侧发生偏移ꎮ 巷道底鼓破坏严重ꎬ产生 非对称底鼓ꎮ 由图 7a可知ꎬ掘进影响稳定阶段的巷道开挖 稳定后垂直应力场呈现规则的对称“V”型分布ꎬ最 大压应力和拉应力分别为 3. 50 107Pa 和 1. 72 106Paꎬ低拉应力主要分布在巷道顶底板处ꎬ主要表 现为拉应力ꎬ以巷道中心向顶底板逐渐递增ꎬ并逐渐 表现为由拉伸破坏向剪切破坏变化ꎬ发生剪切破坏 高压应力主要集中在巷道两帮ꎬ并以巷道中心逐渐 向两帮侧递减ꎻ由图 7b可知ꎬ采动影响的巷道垂 直应力场呈现不规则的偏转“V”形分布ꎬ顶板应力 场向煤柱侧偏移ꎬ底板应力场向实体煤帮偏移ꎬ应力 场值整体增大ꎬ最大压应力和拉应力分别为 9. 40 107和 2. 41 106Paꎬ压应力显著增大ꎬ增幅分别为 168. 6%和 40. 1%ꎮ 压应力和拉应力区分布位置基 本不变ꎮ 岩土塑性力学认为[12 -13]ꎬ岩体的破坏是由偏应 力控制的ꎬ巷道塑性区的发展与偏应力场的演化密 切相关ꎬ由图 7c、d可知ꎬ偏应力场分布规律与 垂直应力场基本一致ꎬ受采动影响深部巷道偏应力 值及分布区域显著增大ꎬ偏应力场向煤柱方“倾倒” 01 2020 年 10 月 毕 鹏等赵固二矿非对称底鼓破坏规律研究 第 29 卷第 10 期 严重ꎬ偏转角度达到 40ꎬ可明显看出采动对巷道煤 柱侧底板挤压作用ꎮ 4. 2 不同阶段下位移场对比分析 运用离散元数值模拟软件 3DEC 导出巷道在掘 进影响稳定阶段和采动影响阶段下底鼓效果、垂直 位移云图ꎬ如图 8 所示ꎮ 图 8 不同阶段巷道位移场对比 由图 8 可知ꎬ深部巷道在掘进影响稳定阶段和 采动阶段围岩位移场存在巨大差异ꎬ掘进影响稳定 阶段的巷道垂直位移、水平位移云图均呈现明显的 以巷道中心为中心的对称性ꎬ采动影响后ꎬ巷道底鼓 效果明显ꎬ周围垂直值增大ꎬ位移场均匀分布被破 坏ꎮ 巷道非对称底鼓破坏严重ꎮ 由图 8c可知ꎬ掘进影响稳定阶段的巷道开挖 稳定后垂直位移场呈均匀上下对称“钟”型分布ꎬ顶 板岩层主要以负向位移即垮落剪切为主ꎬ巷道两帮 及底板主要以正位移即拉伸为主ꎬ并以巷道中心向 顶板和底板两侧变形量都逐渐减小ꎬ最大底鼓位于 巷道底板中心ꎬ底鼓量 9. 46 10 -2 mꎬ最大垮落位 置位于顶板中心ꎬ最大垮落量为 4. 78 10 -2 mꎻ由 图 8d可知受采动影响后ꎬ巷道位移场呈不均匀 “波浪”形ꎬ巷道周围垂直位移值普遍增大ꎬ尤其是 煤柱下底板及巷道靠近煤柱部分底板ꎬ煤柱侧底鼓 范围及底鼓量远大于实体煤一侧ꎬ最大底鼓处位于 距煤柱侧底板1. 8 m 处ꎬ底鼓量达到5. 09 10 -1 mꎬ 增幅达到 438%ꎮ 4. 3 不同阶段下塑性区对比分析 运用离散元数值模拟软件 3DEC 导出巷道在掘 进影响稳定阶段和采动影响阶段下塑性区图ꎬ如图 9 所示ꎮ 由图 9 可知ꎬ深部巷道在掘进影响稳定阶 段和采动阶段塑性区分布存在巨大差异①采动阶 段巷道顶底板及两帮的塑性区破坏范围远大于掘进 影响稳定阶段ꎬ巷道左下角煤柱侧塑性区破坏范 围和深度远大于巷道右下角实体煤侧ꎻ②巷道底 板煤柱侧塑性破坏严重ꎬ破坏深度大ꎬ远大于实体煤 一侧ꎬ呈现非对称破坏ꎮ 图 9 不同阶段塑性区对比 5 结 语 1 通过防爆相机和地质雷达实测发现赵固 二矿 11060 轨道巷在 11041 工作面回采期间底板围 岩变形破坏严重ꎬ巷道呈现非对称底鼓ꎮ 2 深部开采具有强扰动特性ꎬ深部开采时底 板塑性屈服区、卸荷破坏区、卸荷扩展区、触矸破坏 区及峰前损伤区的范围远大于浅部开采ꎬ深部开采 强扰动性产生的触矸也会进一步加大底板破坏范 围ꎬ易产生非对称破坏ꎮ 3 3DEC 数值模拟表明ꎬ深部巷道在掘进影 响稳定阶段和采动阶段围岩应力场、位移场、塑性区 存在巨大差异ꎬ掘进影响稳定阶段的巷道垂直应力、 偏应力云图、垂直位移云图均呈现明显的以巷道中 心为中心的对称性ꎬ采动影响后ꎬ各应力及位移场普 遍向煤柱一侧发生偏移ꎬ巷道底鼓破坏更严重ꎬ产生 非对称底鼓ꎮ 参考文献 [1] 谢和平ꎬ彭苏萍ꎬ何满潮. 深部开采基础理论与工程实 践[M]. 北京科学出版社ꎬ2006. [2] 谢和平ꎬ高 峰ꎬ鞠 杨ꎬ等. 深部开采的定量界定与 分析[J]. 煤炭学报ꎬ2015ꎬ4011 -10. [3] 申李华ꎬ宋作文. 王庄煤矿 91 采区合理煤柱尺寸研究 [J]. 煤ꎬ2019ꎬ28126 -29. [4] 钱鸣高ꎬ石平五. 矿山压力与岩层控制[M]. 徐州中 国矿业大学出版社ꎬ200358 -60. [5] 武精科ꎬ阚甲广ꎬ谢福星ꎬ等. 深井沿空留巷顶板变形 破坏特征与控制对策研究[J]. 采矿与安全工程学报ꎬ 2017116 -23. [6] 马念杰ꎬ赵希栋ꎬ赵志强ꎬ等. 深部采动巷道顶板稳定 性分析与控制[J]. 煤炭学报ꎬ2015ꎬ40102 287 - 2 295. [7] 陈上元ꎬ常 胜ꎬ郭志飚ꎬ等. 深部动压巷道非对称变 形力学机制及控制对策[J]. 煤炭学报ꎬ2016ꎬ411 246 -254. 下转第 14 页 11 2020 年 10 月 毕 鹏等赵固二矿非对称底鼓破坏规律研究 第 29 卷第 10 期 但是增加的幅度逐渐减小ꎮ 出现这种现象的原因 是ꎬ随着吸附时间的增加ꎬ煤体孔隙内瓦斯的吸附量 逐渐增加ꎬ同时裂隙内被一定量的瓦斯所填充ꎬ造成 煤体孔隙、裂隙膨胀ꎬ因此膨胀变形的量逐渐增加ꎮ 膨胀变形增加的幅度逐渐变小的原因是ꎬ随着吸附 的进行ꎬ煤体对瓦斯的吸附能力逐渐减弱ꎬ因此膨胀 变形增加的幅度逐渐变小ꎮ 吸附平衡压力越大ꎬ试 样的体应变越大ꎬ因为吸附平衡压力越大ꎬ煤体孔隙 内吸附的瓦斯越多ꎬ因此其体应变越大ꎮ 煤体的膨 胀变形ꎬ即体应变与时间满足幂函数的函数关系ꎬ采 用幂函数可以对散点进行较好的拟合ꎬ拟合公式如 式2所示ꎬ拟合参数如表 2 所示ꎮ ε = ctd2 式中ε 为体应变ꎻc 和 d 为拟合参数ꎻt 为抽真 空时间ꎮ 图 3 抽真空期间的煤样的收缩变形随时间的变化规律 图 4 吸附期间煤样的膨胀变形随时间的变化规律 表 2 式2中的拟合参数 煤样 编号 吸附平衡 压力/ MPa cdR2 1 号0.712. 047 51 10 -4 0. 195 510.987 47 2 号1.131. 207 73 10 -4 0.321 10.998 65 3 号1.522. 475 85 10 -4 0. 241 080.998 92 3 结 语 1 对煤体进行抽真空时ꎬ煤体会发生收缩变 形ꎬ煤体的体应变随着时间的增加而逐渐增加ꎬ煤样 的体应变与时间满足 Langmuir 的函数关系ꎻ 2 煤体进行吸附瓦斯时ꎬ煤体会发生膨胀变 形ꎬ煤体的体应变随着时间的增加而逐渐增加ꎬ煤体 的体应变与时间满足幂函数的函数关系ꎻ 3 随着吸附平衡压力的增加ꎬ煤体吸附的瓦 斯量逐渐增多ꎬ煤体的体应变逐渐增加ꎮ 参考文献 [1] 罗 维. 双重孔隙结构煤体瓦斯解吸流动规律研究 [D]. 北京中国矿业大学ꎬ2013. [2] 唐志远ꎬ杨延强ꎬ刘国昌. 基于双重孔隙结构的煤层气 储层震波速度预测研究[J]. 煤炭科学技术ꎬ2018ꎬ46 12187 -191. [3] 杨华平. 煤体甲烷吸附解吸机理研究[D]. 西安西安 科技大学ꎬ2014. [4] 张遵国ꎬ赵 丹ꎬ陈 毅. 不同含水率条件下软煤等温 吸附特性及膨胀变形特性[J]. 煤炭学报ꎬ2020. [5] 栗 婧ꎬ汪 振ꎬ戚瑞康ꎬ等. 不同温度下型煤吸附解吸 变形规律的研究[J]. 煤炭技术ꎬ2020ꎬ39183 -88. [6] 翟盛锐. 不同粒度型煤煤样瓦斯吸附 - 解吸变形特征 实验研究[J]. 中国安全生产科学技术ꎬ2018ꎬ146 84 -89. [7] 刑俊旺. 超临界 CO2与 CH4吸附解吸引起煤体变形特 性的对比研究[D]. 太原太原理工大学ꎬ2018. [责任编辑常丽芳] 上接第 11 页 [8] 李春元ꎬ张 勇ꎬ张国军ꎬ等. 深部开采动力扰动下底板 应力演化及裂隙扩展机制[J]. 岩土工程学报ꎬ2018ꎬ40 112 031 -2 040. [9] 李春元ꎬ张 勇ꎬ左建平ꎬ等. 深部开采砌体梁失稳扰动 底板破坏力学行为及分区特征[J]. 煤炭学报ꎬ2019ꎬ44 51 508 -1 520. [10] 魏文胜ꎬ李春元ꎬ司艳龙ꎬ等. 下伏煤层空巷对上覆煤 层开采的影响[J]. 煤矿安全ꎬ2019ꎬ506258 -262. [11] 李春元ꎬ张 勇. 深埋薄基岩顶板来压与底板破坏深 度关系[J]. 煤炭科学技术ꎬ2016ꎬ44874 -79. [12] 陈希哲. 土力学地基基础[M]. 北京清华大学出版 社ꎬ2004. [13] 王 仁. 塑性力学引论[M]. 北京北京大学出版社ꎬ 1982. [责任编辑常丽芳] 41 2020 年 10 月 李会云煤体的吸附膨胀变形特性研究 第 29 卷第 10 期