重庆大学学报  2023, Vol. 46 Issue (2): 43-56  DOI: 10.11835/j.issn.1000-582X.2022.263 RIS(文献管理工具)
0

引用本文 

秦冰, 叶栋林, 石占山, 孙维吉, 李刚, 张晋京. 双重上保护层叠加开采应力分布规律[J]. 重庆大学学报, 2023, 46(2): 43-56. DOI: 10.11835/j.issn.1000-582X.2022.263.
QIN Bing, YE Donglin, SHI Zhanshan, SUN Weiji, LI Gang, ZHANG Jinjing. Pressure relief law of double upper protective layer superposition mining[J]. Journal of Chongqing University, 2023, 46(2): 43-56. DOI: 10.11835/j.issn.1000-582X.2022.263.

基金项目

国家自然科学基金青年基金资助项目(52004118);辽宁省教育厅青年科技人才"育苗"资助项目(LJ2020QNL009);辽宁工程技术大学学科创新团队项目(LNTU20TD-11)

通信作者

石占山, 男, 博士, 主要从事瓦斯灾害防治方面的研究, (E-mail)shizhanshan@163.com

作者简介

秦冰(1987—), 博士, 主要从事渗流力学及煤与瓦斯共采方面的研究, (E-mail)qinbing20071111@163.com

文章历史

收稿日期: 2022-02-19
双重上保护层叠加开采应力分布规律
秦冰 1a, 叶栋林 1b, 石占山 1b,1c, 孙维吉 1a, 李刚 1b, 张晋京 2     
1a. 辽宁工程技术大学 力学与工程学院, 辽宁 阜新 123000;
1b. 辽宁工程技术大学 矿业学院, 辽宁 阜新 123000;
1c. 辽宁工程技术大学 矿产资源开发利用技术及装备研究院, 辽宁 阜新 123000;
4. 中国平煤神马能源化工集团有限责任公司, 河南 平顶山 467000
摘要: 为获取双重上保护层重叠采动作用下的煤层卸压规律及保护层间的相互影响,以平煤八矿一采区为原型,采用FLAC3D软件模拟了丁、戊组煤层多工作面重叠开采过程。研究结果表明,仅丁组煤层开采时,采区边界煤柱对应范围出现应力集中现象,最大应力值达到19 MPa,影响范围达到下方80 m,不利于被保护层卸压。工作面间区段煤柱最大应力值达39 MPa,但向底板传递范围较小。丁组单独开展卸压区域能够影响至己组,己组煤层应力卸压值约为1 MPa;丁组、戊组煤层重叠开采,当戊组工作面位于丁组煤层区段煤柱下方,同时丁组工作面位于戊组工作面区段煤柱上方,使丁组煤层工作面间区段煤柱应力集中减弱利于卸压,丁组煤层区段煤柱应力值由39 MPa卸载至7.5~10 MPa之间;同时当己组煤层位于戊组单独保护范围时,垂直应力卸压值为2~3 MPa,当己组煤层位于丁戊共同保护范围时,垂直应力卸压值为4~6.5 MPa。
关键词: 双重保护层    采动应力    叠加开采    数值模拟    煤柱集中应力    
Pressure relief law of double upper protective layer superposition mining
QIN Bing 1a, YE Donglin 1b, SHI Zhanshan 1b,1c, SUN Weiji 1a, LI Gang 1b, ZHANG Jinjing 2     
1a. School of Mechanics & Engineering, Liaoning Technical University, Fuxin, Liaoning 123000, P. R. China;
1b. School of Mines, Liaoning Technical University, Fuxin, Liaoning 123000, P. R. China;
1c. Liaoning Academy of Mineral Resources Development and Utilization Technical and Equipment Research Institute, Liaoning Technical University, Fuxin, Liaoning 123000, P. R. China;
4. China Pingmei Shenma Group, Pingdingshan, Henan 467000, P. R. China
Abstract: To obtain the law of stress redistribution of coal seam under the action of overlapping mining in multi seam working face, taking Pingmei No.8 Coal Mine No.1 mining area as the prototype, the process of overlapping mining in group D and group E coal seams was simulated using FLAC3D software. The results show that the stress concentration occurs in the corresponding range of the coal pillar at the boundary of the mining area only when the coal seam in group D is mined. The maximum stress value of the coal pillar reaches 19 MPa, and the influence range is 80 m below it. The maximum stress value of the coal pillar in the section between working faces reaches 39 MPa. The pressure relief affects J group, and the stress relief value of J Group coal seam is about 1 MPa. When the D group and E group coal seams are overlapped, and when the E group working face is located under the coal pillar of the D group coal seam section, and the D group working face is located above the coal pillar of the E group working face section, the stress concentration of the coal pillar of the D group coal seam section is weakened, and the stress value of the coal pillar of the D group coal seam section is unloaded from 39 MPa to the range of 7.5 MPa to 10 MPa. At the same time, when the coal seam of group I is located in the independent protection range of group E, the vertical stress relief value is 2 MPa to 3 MPa. When the coal seam of group I is located in the joint protection range of group E, the vertical stress relief value is from 4 MPa to 6.5 MPa.
Keywords: double upper protective layer    mining stress    superposition mining    numerical simulation    concentrated stress of coal pillar    

保护层开采后,下伏煤岩体应力降低、采动裂隙发育、煤岩层透气性增大,使煤层瓦斯解吸、运移[1-2],利于通过瓦斯抽采降低煤层瓦斯含量,在国内已有较成熟的应用。目前,针对保护层开采应力重分布规律已有大量研究。石必明等[3-4]研究了保护层开采上覆岩层的应力变化以及膨胀变形规律;施龙青等[5]和朱第植等[6]从工程力学角度出发,对岩石损伤及断裂特征进行了大量的实验,并结合工程实践将釆场下伏煤岩体分为“四带”;沈明荣等[7]分析了工作面回釆后采场底板及下覆煤岩体向釆场无应力区弹性恢复的过程;李树清等[8]和袁志刚等[9]应用岩土力学数值计算软件FLAC模拟了保护层的开采过程并给出卸压角及保护效果影响参数;黄光利[10]运用COMSOL Multiphysics软件的固气耦合模型计算了煤矿上保护层开采的保护范围;朱志洁等[11]开展了远距离重叠开采条件下,煤柱集中应力影响下的强矿压显现机制,分析了煤柱对强矿压显现的影响机理;Gao等[12]通过数值计算分析了应力分布特征,并考虑了上覆煤柱的存在,同时揭示了采空区应力释放与上覆煤柱应力集中对最终应力分布的共同影响;李杨等[13]开展了多煤层开采中,中间岩层对覆岩移动的影响机理研究,提出了上覆岩层与中间岩层厚度之比并不是控制覆岩移动的决定性因素;Yin等[14]研究考察了重庆地区煤样在开采过程中应力路径下的应力应变渗透率关系。康钦容等[15]利用自行研制的“多场耦合煤矿动力灾害大型模拟试验系统”,进行了三维采动应力条件下的三维模拟开采实验,对多煤层卸压采动后的底板岩层破坏规律进行了研究。

然而当前研究中大多考虑单一上保护开采工况[16-21],且以单一工作面为研究对象,忽略了上保护层多工作面布置时区段煤柱产生的应力集中对下被保护层卸压的影响,及多重上保护层开采的卸压增强效应。同时在保护层效果判定中,《AQ1050—2008保护层开采技术规范》[22]中给出,上保护层开采时对于缓倾斜煤层保护最大有效垂距为50 m,规范给出的保护有效距离是基于单一保护层开采提出的,对于多重上保护层叠加开采保护有效距离的评价未见相关规范指南,多重上保护层开采的卸压效果需要具体分析。

双重上保护层开采应力重分布不同于单一煤层开采。对于双重上保护层开采工作面,不仅能够对下方煤层进行双重卸压,而且在工作布置时能够对双重保护层内区段煤柱进行相互掩护,双重卸压规律及相互掩护对卸压的影响有待进一步研究。因此笔者充分考虑了现场实际工况,以平煤八矿丁5-6煤层、戊9-10煤层、己15煤层组成的煤层群开采为工程背景,依据现场多煤层多工作面的实际采掘布置,建立了双重上保护层多工作面开采数值计算模型,对双重上保护层开采采动应力分布规律开展研究,进一步依据应力重分布规律给出被保护煤层的瓦斯治理措施及采掘布置计划,达到瓦斯精准治理[23]的目标。

1 FLAC3D模型的建立及参数确定 1.1 模型的建立

平煤八矿为高瓦斯矿井,主要可采煤层有丁5-6、戊9-10、己15和己16-17煤层,本次研究重点区域为丁5-6、戊9-10、己15煤层组成的煤层群,其中丁5-6煤层和戊9-10煤层间距是70~90 m,戊9-10煤层和己15煤层间距平均约170 m,煤层倾角平均9°。研究区域丁5-6煤层布置有4个工作面,戊9-10煤层布置有2个工作面,且均已开采完毕。现阶段拟开采己15煤层21030工作面。采区内丁5-6、戊9-10煤层的工作面走向相同,长短不一,在倾向上交错布置。己15煤层21030工作面走向与丁5-6、戊9-10煤层工作面走向呈17°交角。拟分析丁5-6、戊9-10煤层双重上保护层开采对己15煤层工作面应力分布的影响。各工作面空间位置如图 1所示。

图 1 工作面空间分布 Fig. 1 Spatial distribution of working faces

根据平煤八矿确定的研究区域工作面布置情况,采用FLAC3D数值计算软件,建立了等比例模型,模型尺寸为800 m(x)×1 400 m(y)×500 m(z)生成网格884 800个,节点913 680。丁5-6、戊9-10煤层间距选取80 m,戊9-10、己15煤层间距选取170 m。模型上表面距离地表约400 m,覆岩平均容重取20 kN/m3,上表面施加载荷为8 MPa。模拟煤层倾角为9°,采用摩尔库仑本构模型,模型前后、左右4个面约束其法向自由度,底面约束xyz 3个方向自由度。建立的模型如图 2所示。

图 2 数值模型建立 Fig. 2 Numerical modeling
1.2 计算参数及数值模拟设计方案 1.2.1 计算参数

依据现场取样在实验室进行的参数测试结果,并结合现有论文中关于平煤股份八矿煤岩物理力学参数测试结果[24-25],考虑岩层与岩石试件在应力换算中存在的尺度效应,依据数值模拟计算结果反演模拟计算采用的岩体力学参数如表 1所示。

表 1 数值计算岩体力学参数 Table 1 Rock mass mechanical parameters of numerical calculation
1.2.2 数值模拟设计方案

开挖过程如下:①对丁5-6煤层工作面全部开采;②开挖后计算平衡,再进行戊9-10煤层工作面全部开采;③戊9-10煤层工作面开采后计算平衡,对己15煤层21030工作面进行分次开采模拟回采过程,工作面每次推进距离为100 m,每次推进后计算平衡,共计开采12次。

为获取实际工作面布置条件下双重上保护层开采应力分布规律,设计了5种方案进行对比。①丁5-6煤层开挖后,通过选取截面A-A分析丁5-6煤层单一工作面保护时,己15煤层应力分布规律。通过选取截面B-B分析丁5-6煤层多工作面共同保护时,己15煤层应力分布规律;②丁5-6、戊9-10煤层均开挖后,通过选取截面C-C分析戊9-10煤层单一工作面保护时,己15煤层应力分布规律。通过选取截面D-D分析丁5-6、戊9-10煤层多工作面交错掩护布置共同保护时,己15煤层应力分布规律。通过选取截面E-E分析丁5-6、戊9-10煤层工作面未相互掩护共同保护时,己15煤层应力分布规律。截面布置如图 1(b)所示,通过多方案对比阐明双重上保护层开采对被保护层应力分布的影响,如表 2所示。

表 2 模拟方案 Table 2 Simulation schemes
2 双重保护层叠加开采应力分布规律研究 2.1 单一保护层开采被保护层倾向应力分布规律

模型初始应力平衡后,选取了y=700 m截面即模型中部,给出各煤层工作面位置的初始应力分布状态,如图 3所示。其中应力值为负值表示该位置受压,反之受拉,卸压值为卸压后应力值与初始应力值之差。图 3中水平应力值方向对应图 1(a)x轴所示的方向。图 3中标记了截面对应的工作面分布位置,图中数字为各工作面名称。

图 3 倾向y=700 m截面初始应力分布状态 Fig. 3 The initial stress distribution state of tendency y=700 m section
2.1.1 丁5-6煤层单独保护己15煤层应力分布规律

为对比双重保护层开采与单一煤层保护层开采应力分布的差异,开展了丁5-6煤层工作面单独保护己15煤层时的应力分布规律研究。在丁5-6煤层工作面开挖后戊9-10煤层工作面未开挖前,丁5-6煤层对己15煤层存在两种保护形式,分别为丁5-6单一工作面开采保护己15煤层(y=200 m截面);丁5-6多工作面开采保护己15煤层(y=500 m截面)。

1) 丁5-6单一工作面保护己15煤层。如图 4所示,该截面为表 2中方案1。丁5-6单一工作面开挖,戊9-10、己15煤层未开挖,丁5-6开挖工作面两侧煤柱发生应力集中现象,应力分布状态如表 3所示。单一工作面开挖后两侧为采区边界煤柱,煤柱集中应力最大值为-19 MPa,应力集中传递至底板下方80 m。丁5-6保护下,己15煤层卸压区域垂直应力值由-14~-16 MPa下降至-13~-15 MPa,下降1 MPa。由于单一保护层边界煤柱集中应力传递范围较大的影响,对己15煤层21030工作面布置位置卸压作用不明显。

图 45-6单一工作面保护己15工作面应力分布云图(y=200 m截面) Fig. 4 The stress distribution nephogram of the single working face of D5-6 protecting J15 working face (y=200 m section)
表 35-6单一工作面保护己15工作面应力分布 Table 3 The stress distribution of the single working face of D 5-6 protecting J 15 working face  

2) 丁5-6多工作面保护己15煤层。如图 5所示,该截面为表 2中方案2。丁5-6多工作面开挖,戊9-10、己15煤层未开挖,应力分布状态如表 4所示。丁5-6煤层工作面间的区段煤柱应力集中值增加为-39 MPa,应力值较大,但其向底板传递的距离较短,约为26 m。下方己15煤层工作面所在位置最大卸压值2 MPa。由于丁5-6煤层多工作面开采,且己15煤层工作面不在边界煤柱下方,己15煤层21030工作面位置卸压值有所增加。

图 55-6多工作面保护己15工作面应力分布云图(y=500 m截面) Fig. 5 The stress distribution nephogram of the multiple working faces of D5-6 protecting J15 working face (y=500 m section)
表 45-6多工作面保护己15工作面应力分布 Table 4 The stress distribution of the multiple working faces of D5-6 protecting J15 working faced  
2.1.2 戊9-10煤层单独保护己15煤层应力分布规律

图 6所示,该截面为表 2中方案3。丁5-6工作面开挖后对戊9-10工作面进行了开挖。由于戊9-10煤层工作面走向长度不同,该截面戊9-10煤层对应单一工作面,同时戊9-10煤层单一工作面左半部分上方有丁5-6工作面,丁5-6煤层工作面右侧煤柱位于戊9-10煤层工作面采空区上方,戊9-10工作面左侧煤柱位于丁5-6工作面采空区下方,使得区段煤柱应力集中现象减弱。戊9-10工作面右半部分形成单一煤层保护己15煤层工作面的条件,该保护区域内己15被保护层垂直应力由-14~-16 MPa卸载至-7.5~-12 MPa之间,如表 5所示,右侧单独保护区域卸压最大值为4 MPa。

图 69-10单一工作面保护己15工作面应力分布云图(y=400 m截面) Fig. 6 The stress distribution nephogram of the single working face of E9-10 protecting J15 working face (y=400 m section)
表 59-10单一工作面保护己15工作面应力分布 Table 5 The stress distribution of the single working face of E9-10 protecting J15 working face  
2.2 双重保护层开采被保护层倾向应力分布规律 2.2.1 丁5-6、戊9-10工作面区段煤柱相互掩护

图 7所示,该截面为表 2中方案4。此时丁5-6、戊9-10工作面均开采完毕,己15煤层21030工作面位于双重保护层下方。由表 6可知,当己15工作面位于丁5-6、戊9-10煤层工作面叠加采动卸压区域下方时,卸压范围进一步向下延伸,己15煤层卸压值进一步增加,由-14~-16 MPa卸压至-7.5~-10 MPa,卸压值最大达6.5 MPa。同时,对应的丁5-6、戊9-10煤层工作面间的区段煤柱的集中应力降低。丁5-6工作面右侧集中应力是由于丁5-6右侧存在一个较短工作面,该截面恰好位于该工作面走向边界导致。

图 75-6、戊9-10工作面交错布置保护己15工作面应力分布云图(y=580 m截面) Fig. 7 The stress distribution nephogram of D5-6 and E9-10 working faces staggered arrangement to protect Ji15 working face (y=580 m section)
表 65-6、戊9-10工作面交错布置保护己15工作面应力分布 Table 6 The stress distribution of D5-6 and E9-10 working faces staggered arrangement to protect J15 working face  
2.2.2 丁5-6、戊9-10煤层工作面区段煤柱未相互掩护

图 8所示,该截面为表 2中方案5。此时丁5-6、戊9-10工作面均开采完毕,丁5-6、戊9-10煤层工作面区段煤柱未相互掩护。由表 7可知,己15煤层由-14~-16 MPa卸压至-11~-13 MPa,卸压值最大达3 MPa。同时,对应的丁5-6、戊9-10煤层工作面间区段煤柱集中应力增加,最大集中应力达到34 MPa。

图 85-6、戊9-10工作面未交错布置保护己15工作面应力分布云图(y=260 m截面) Fig. 8 The stress distribution nephogram of D5-6 and E9-10 working faces without staggered arrangement to protect J15 working face (y=260 m section)
表 75-6、戊9-10工作面未交错布置保护己15工作面应力分布 Table 7 The stress distribution of D5-6 and E9-10 working faces without staggered arrangement to protect J15 working face  
2.3 双重保护层开采被保护层走向应力分布规律

前述给出了丁5-6、戊9-10煤层工作面开采不同保护类型时,己15煤层工作面倾向应力分布规律。为进一步确定己15煤层工作面走向应力分布情况,取丁5-6、戊9-10煤层开采后,己15煤层21030工作面的走向切面应力进行分析,如图 9所示。

图 95-6、戊9-10工作面开采后己15-21030工作面走向切面 Fig. 9 The trend section of J15-21030 working face after the mining of D5-6 and E9-10

图 9(b)(c)可知,己15煤层工作面在走向方向上卸压角不同。当己15煤层工作面走向边界位于戊9-10煤层工作面单独保护区域时,卸压角较大为65°;位于丁5-6煤层工作面单独保护区域时,对应的卸压角较小为50°,如图 9(b)中虚线所示,且水平应力变化较小。

图 10为丁5-6、戊9-10煤层开采前后己15煤层21030工作面的走向应力分布对比。由图 10可知,最大卸压值位于工作面走向300~700 m范围内,此处为双重保护层叠加保护区域,卸压影响范围为工作面走向100 m到1 150 m范围内。

图 105-6、戊9-10煤层开采前后己15-21030工作面走向切面应力值对比 Fig. 10 Comparison of the stress values of trend section of J15-21030 working face before and after the mining of D5-6 and E9-10

图 11为丁5-6、戊9-10煤层开采后,己15煤层全区的应力分布情况。由图 11可知,在双重保护层保护区域的卸压值较大,且呈椭圆形。当丁5-6煤层工作面或戊9-10煤层工作面单独保护时,卸压区域不规则延伸,且卸压值减小。

图 115-6、戊9-10开采后己组煤层全区应力分布 Fig. 11 Stress distribution in the whole area of the Group J coal seam after the mining of D5-6 and E9-10
2.4 双重保护层开采后己15煤层工作面回采应力演化规律

为指导己15煤层工作面采掘计划设计并优化瓦斯抽采钻孔布置,进一步分析了己15煤层21030工作面回采过程中工作面前方支承压力应力集中系数的演化规律。依据图 12所示应力云图,统计应力集中系数如表 8所示。

图 1215煤层回采应力演化规律 Fig. 12 Law of evolution of mining stress in coal seam of J15
表 8 工作面推进应力集中系数 Table 8 Stress concentration factor during working face advance

推进过程中距开切眼300~1 000 m范围内应力集中系数相对较小。最大应力集中系数为3.40,最小应力集中系数为1.98,在回采中应注意对应力集中系数较大区域提前给出安全防治措施。

3 应力分布规律现场验证 3.1 测试方法及原理

钻屑量是一种综合反映煤体应力、瓦斯压力和煤的物理力学性质的指标。钻孔形成后,在其周围将形成破碎区、塑性软化区和弹性区,不同围岩应力环境下钻孔三区分布范围不同,导致钻屑量不同,钻屑量测试能间接反映钻孔围岩应力分布情况,围岩应力值越大,钻屑量越大。

3.2 测试方案

针对丁5-6、戊9-10煤层对己15煤层的保护类型,共设计了7个考察区域布置了13个施工位置开展钻孔钻屑量考察,如表 9所示。

表 9 分区情况说明 Table 9 Partition description
3.3 测试过程及结果

静水压力下巷道围岩处于弹性变形阶段,切向及径向集中应力影响范围为3~5倍巷道半径,巷道宽度为4.8 m,估算应力集中影响范围为14.4~24 m。剔除巷道周围应力重分布影响,从钻进至24 m处开始取屑,直至钻进至79 m终止,测试位置如图 13所示,部分测试结果如图 14所示。

注:图中红色线为丁5-6煤层工作面;蓝色线为戊9-10煤层工作面;绿色线为己15煤层工作面 图 1315-21030工作面钻屑量测试位置 Fig. 13 Drilling cuttings weight test position at J15-21030 working face
图 14 不同区域钻屑量测试结果 Fig. 14 Test results of drilling cuttings weight in different areas
3.4 测试结果分析

数值计算表明己15煤层原始区域应力值为14~16 MPa,区域五、区域六钻屑量测试对应该方案,平均钻屑量16.5 kg/m,如图 14(a)所示;区域四钻屑量测试对应方案2,最大卸压值为2 MPa,平均钻屑量14.8 kg/m;区域二钻屑量测试对应方案4,叠加采空区下方己15煤层卸压情况,卸压区域最大卸压值6.5 MPa,平均钻屑量10.8 kg/m;区域一钻屑量测试对应方案4,掩护煤柱下方己15煤层卸压情况,卸压值为4 MPa,平均钻屑量为11.3 kg/m;区域三钻屑量测试对应方案5,丁5-6、戊9-10工作面煤柱未相互掩护,区域最大集中应力值可达到34 MPa,最大卸压值3 MPa,平均钻屑量12.3 kg/m;区域七对应该方案3,戊9-10单一工作面保护己15煤层,卸压最大值为4 MPa,平均钻屑量11.7 kg/m,如图 14(b)所示。图 14中给出了区域六、七的钻屑量测试结果,其他区域钻屑量值测试可参考文献[26],各区域测试结果均能反映卸压值变化情况。

4 现场应用

现场实测表明,己15-21030工作面瓦斯含量及压力较大,因此需要科学制定采掘计划及抽采钻孔布置方案。依据数值计算结果,卸压区范围为工作面走向100 m到1 150 m,最大卸压值位于工作面走向300 m到700 m范围内,应力集中系数在300~1 000 m范围内较小。当前正常情况下巷道日掘进进尺为8 m/d,根据卸压区分布情况,将卸压区外掘进时的进尺调整为4 m/d,工作面走向300~700 m范围内的掘进进尺按原计划8 m/d进行,其余区域掘进进尺调整为6 m/d。原设计本煤层预抽钻孔间距为2.5 m,在卸压区外将钻孔间距调整为1.5 m,走向300~700 m范围内按原计划2.5 m间距设计,其余区域调整为2.0 m间距。回采过程中在300~1 000 m范围内保持正常回采速度,在该范围外降低工作面推进速度。调整设计方案后,依据后续掘进反馈数据,未发生瓦斯突出及超限事故。

5 结论

1) 丁5-6煤层单一工作面保护己15煤层时,煤柱为边界煤柱,煤柱应力集中值达到19 MPa,影响范围达到下方80 m。卸压影响至己15煤层,己15煤层垂直应力最大卸压值为1 MPa。戊9-10单一工作面保护己15煤层时,垂直应力卸压值最大达4 MPa。丁5-6煤层多工作面保护己15煤层时,丁5-6工作面区段煤柱应力集中值达39 MPa,向下影响范围较小为26 m,对己15煤层最大卸压值为2 MPa,水平应力变化不明显。边界煤柱集中应力较区段煤柱小,但向下影响范围大。

2) 丁5-6、戊9-10煤层双重保护时,工作面区段煤柱相互掩护时,垂直应力卸压值最大达6.5 MPa。丁5-6、戊9-10煤层工作面相互掩护使得工作面区段煤柱集中应力由39 MPa卸载至10 MPa。水平应力卸压值变化不明显。工作面区段煤柱未相互掩护时,垂直应力卸压值最大为3 MPa,煤柱集中应力制约了被保护层的卸压。丁、戊煤层双重叠加掩护卸压增强了卸压效应降低了区段煤柱的应力集中,较单一煤层开采明显提高了被保护层的卸压值。

3) 开展了现场工作面钻屑量测试,验证了双重保护层开采应力分布规律。得到己15煤层21030工作面卸压影响范围为工作面100 m到1 150 m,最大卸压值位于工作面走向300 m到700 m范围内,应力集中系数在300~1 000 m范围内较小。依据上述规律制定了采掘计划及抽采钻孔设计方案。

参考文献
[1]
耿朝光, 马建宏, 侯江涛. 中泰矿业下保护层开采数值模拟研究[J]. 能源技术与管理, 2018, 43(1): 4-6, 31.
Gen C G, Ma J H, Hou J T. Research on Numerical simulation of mining protective layer under Zhongtai mining[J]. Energy Technology and Management, 2018, 43(1): 4-6, 31. (in Chinese)
[2]
Zhang H T, Wen Z H, Yao B H, et al. Numerical simulation on stress evolution and deformation of overlying coal seam in lower protective layer mining[J]. Alexandria Engineering Journal, 2020, 59(5): 3623-3633. DOI:10.1016/j.aej.2020.06.011
[3]
石必明, 俞启香. 远距离保护层开采煤岩移动变形特性的试验研究[J]. 煤炭科学技术, 2005, 33(2): 39-41, 45.
Shi B M, Yu Q X. Test and research on coal and rock movement deformation characteristics in coal mining with long distance protection layer[J]. Coal Science and Technology, 2005, 33(2): 39-41, 45. (in Chinese)
[4]
石必明, 俞启香, 周世宁. 保护层开采远距离煤岩破裂变形数值模拟[J]. 中国矿业大学学报, 2004, 33(3): 259-263.
Shi B M, Yu Q X, Zhou S N. Numerical simulation of far-distance rock strata failure and deformation caused by mining protecting stratum[J]. Journal of China University of Mining & Technology, 2004, 33(3): 259-263. (in Chinese)
[5]
施龙青, 韩进. 底板突水机理及预测预报[M]. 徐州: 中国矿业大学出版社, 2004.
Shi L Q, Han J. Floor Water-inrush Mechanism and Prediction[M]. Xuzhou, China: China University of Mining & Technology Press, 2004. (in Chinese)
[6]
朱第植, 王成绪. 原位应力测试在底板突水预测中的应用[J]. 煤炭学报, 1998, 23(3): 295-299.
Zhu D Z, Wang C X. Application of in situ stress measurement to forecast of water inrush from floor[J]. Journal of China Coal Society, 1998, 23(3): 295-299. (in Chinese)
[7]
沈明荣, 陈建峰. 岩体力学[M]. 上海: 同济大学出版社, 2006.
Shen M R, Chen J F. Rock mechanics[M]. Shanghai: Tongji University Press, 2006. (in Chinese)
[8]
李树清, 龙祖根, 罗卫东, 等. 煤层群下保护层开采保护范围的数值模拟[J]. 中国安全科学学报, 2012, 22(6): 34-40.
Li S Q, Long Z G, Luo W D, et al. Numerical simulation of protection scope when lower-protective layer mined in coal seams[J]. China Safety Science Journal, 2012, 22(6): 34-40. (in Chinese)
[9]
袁志刚, 王宏图, 胡国忠, 等. 急倾斜多煤层上保护层保护范围的数值模拟[J]. 煤炭学报, 2009, 34(5): 594-598.
Yuan Z G, Wang H T, Hu G Z, et al. Numerical simulation for protection scope of upper protective seam in steeply inclined multi-coal seam[J]. Journal of China Coal Society, 2009, 34(5): 594-598. (in Chinese)
[10]
黄光利. 砚石台煤矿急倾斜俯伪斜上保护层开采保护范围研究[D]. 重庆: 重庆大学, 2014.
Huang G L. Study on protection scope of steep-inclined upper-protective layer with pseudo-incline technique for yanshitai coal mine[D]. Chongqing: Chongqing University, 2014. (in Chinese)
[11]
朱志洁, 张宏伟, 陈蓥, 等. 远距离重叠煤柱作用下的强矿压显现作用机制: 以大同矿区多煤层开采为例[J]. 中国安全科学学报, 2017, 27(3): 83-88.
Zhu Z J, Zhang H W, Chen Y, et al. Strong strata behavior mechanism under long distance overlapping coal pillar: taking multiple coal seams mining of Datong mining area as an example[J]. China Safety Science Journal, 2017, 27(3): 83-88. (in Chinese)
[12]
Gao R, Yu B, Meng X B. Stress distribution and surrounding rock control of mining near to the overlying coal pillar in the working face[J]. International Journal of Mining Science and Technology, 2019, 29(6): 881-887.
[13]
李杨, 王建鹏, 陈一鼎, 等. 多煤层开采中间岩层对覆岩移动的影响研究[J]. 煤炭科学技术, 2020, 48(4): 246-255.
Li Y, Wang J P, Chen Y D, et al. Study on effect of interburden on movement of overburden in multiple coal seams[J]. Coal Science and Technology, 2020, 48(4): 246-255. (in Chinese)
[14]
Yin G Z, Li M H, Wang J G, et al. Mechanical behavior and permeability evolution of gas infiltrated coals during protective layer mining[J]. International Journal of Rock Mechanics and Mining Sciences, 2015, 80: 292-301.
[15]
康钦容, 张卫中, 韩佩博, 等. 三维应力条件下煤层底板采动破坏规律实验研究[J]. 地下空间与工程学报, 2017, 13(3): 612-618.
Kang Q R, Zhang W Z, Han P B, et al. Experimental study of mine-induced underlying strata failure regularity in three dimensional stress condition[J]. Chinese Journal of Underground Space and Engineering, 2017, 13(3): 612-618. (in Chinese)
[16]
Yin W, Miao X X, Zhang J X, et al. Mechanical analysis of effective pressure relief protection range of upper protective seam mining[J]. International Journal of Mining Science and Technology, 2017, 27(3): 537-543.
[17]
Wu X Q, Dou L M, Lv C G, et al. Research on pressure-relief effort of mining upper-protective seam on protected seam[J]. Procedia Engineering, 2011, 26: 1089-1096.
[18]
Yao B H, Ma Q Q, Wei J P, et al. Effect of protective coal seam mining and gas extraction on gas transport in a coal seam[J]. International Journal of Mining Science and Technology, 2016, 26(4): 637-643.
[19]
Liu Z, Yang H, Cheng W M, et al. Stress distribution characteristic analysis and control of coal and gas outburst disaster in a pressure-relief boundary area in protective layer mining[J]. Arabian Journal of Geosciences, 2017, 10(16): 1-15.
[20]
Gao R, Yu B, Xia H C, et al. Reduction of stress acting on a thick, deep coal seam by protective-seam mining[J]. Energies, 2017, 10(8): 1209.
[21]
Song Z M, Gao X C, Tian K Y. Determination of protection range of mining upper protective layers and its numerical simulation[J]. Journal of Coal Science and Engineering (China), 2012, 18(4): 368-373.
[22]
国家安全生产监督管理总局. 保护层开采技术规范: AQ 1050—2008[S]. 北京: 煤炭工业出版社, 2009.
State Administration of Quality and Technical Supervision of the People's Republic of China. Technical criterion of protective coal seam exploitation: AQ 1050—2008[S]. Beijing: China Coal Industry Publishing House, 2009. (in Chinese)
[23]
刘军, 周东平, 蒋和财, 等. 小甘沟煤矿工作面瓦斯涌出分析与精准治理研究[C]//2019年西南五省(市、区)煤炭学术年会(重庆部分)论文集. 重庆: 重庆市煤炭学会, 2019: 98-104.
Liu J, Zhou D P, Jiang H C, et al. Study on gas emission analysis and precision control in Xiaogangou Coal Mine working face[C]//Proceedings of the 2019 Southwest Five Provinces (City, District) Coal Academic Annual Conference (Chongqing Part). Chongqing: Chongqing Coal Society, 2019: 98-104. (in Chinese)
[24]
李建伟, 刘艳辉. 平煤八矿远距离保护层开采可行性分析[J]. 矿业安全与环保, 2011, 38(2): 80-81, 84.
Li J W, Liu Y H. Feasibility analysis of long-distance protective seam mining in Pingmei No. 8 Mine[J]. Mining Safety & Environmental Protection, 2011, 38(2): 80-81, 84. (in Chinese)
[25]
翟新献, 涂兴子, 李如波, 等. 平煤八矿硐室围岩控制联合支护技术研究[J]. 煤炭技术, 2018, 37(5): 1-3.
Zhai X X, Tu X Z, Li R B, et al. Study on combined support technology of surrounding rocks deformation control of chamber in Pingdingshan No.8 coal mine[J]. Coal Technology, 2018, 37(5): 1-3. (in Chinese)
[26]
周洋, 梁冰, 石占山, 等. 考虑煤层倾角和钻孔倾角的钻屑量理论研究[J]. 矿业安全与环保, 2020, 47(1): 26-30, 35.
Zhou Y, Liang B, Shi Z S, et al. Theoretical study on the drilling cuttings considering coal seam angle and borehole angle[J]. Mining Safety & Environmental Protection, 2020, 47(1): 26-30, 35. (in Chinese)
图 1 工作面空间分布 Fig. 1 Spatial distribution of working faces
图 2 数值模型建立 Fig. 2 Numerical modeling
表 1 数值计算岩体力学参数 Table 1 Rock mass mechanical parameters of numerical calculation
表 2 模拟方案 Table 2 Simulation schemes
图 3 倾向y=700 m截面初始应力分布状态 Fig. 3 The initial stress distribution state of tendency y=700 m section
图 45-6单一工作面保护己15工作面应力分布云图(y=200 m截面) Fig. 4 The stress distribution nephogram of the single working face of D5-6 protecting J15 working face (y=200 m section)
表 35-6单一工作面保护己15工作面应力分布 Table 3 The stress distribution of the single working face of D 5-6 protecting J 15 working face  
图 55-6多工作面保护己15工作面应力分布云图(y=500 m截面) Fig. 5 The stress distribution nephogram of the multiple working faces of D5-6 protecting J15 working face (y=500 m section)
表 45-6多工作面保护己15工作面应力分布 Table 4 The stress distribution of the multiple working faces of D5-6 protecting J15 working faced  
图 69-10单一工作面保护己15工作面应力分布云图(y=400 m截面) Fig. 6 The stress distribution nephogram of the single working face of E9-10 protecting J15 working face (y=400 m section)
表 59-10单一工作面保护己15工作面应力分布 Table 5 The stress distribution of the single working face of E9-10 protecting J15 working face  
图 75-6、戊9-10工作面交错布置保护己15工作面应力分布云图(y=580 m截面) Fig. 7 The stress distribution nephogram of D5-6 and E9-10 working faces staggered arrangement to protect Ji15 working face (y=580 m section)
表 65-6、戊9-10工作面交错布置保护己15工作面应力分布 Table 6 The stress distribution of D5-6 and E9-10 working faces staggered arrangement to protect J15 working face  
图 85-6、戊9-10工作面未交错布置保护己15工作面应力分布云图(y=260 m截面) Fig. 8 The stress distribution nephogram of D5-6 and E9-10 working faces without staggered arrangement to protect J15 working face (y=260 m section)
表 75-6、戊9-10工作面未交错布置保护己15工作面应力分布 Table 7 The stress distribution of D5-6 and E9-10 working faces without staggered arrangement to protect J15 working face  
图 95-6、戊9-10工作面开采后己15-21030工作面走向切面 Fig. 9 The trend section of J15-21030 working face after the mining of D5-6 and E9-10
图 105-6、戊9-10煤层开采前后己15-21030工作面走向切面应力值对比 Fig. 10 Comparison of the stress values of trend section of J15-21030 working face before and after the mining of D5-6 and E9-10
图 115-6、戊9-10开采后己组煤层全区应力分布 Fig. 11 Stress distribution in the whole area of the Group J coal seam after the mining of D5-6 and E9-10
图 1215煤层回采应力演化规律 Fig. 12 Law of evolution of mining stress in coal seam of J15
表 8 工作面推进应力集中系数 Table 8 Stress concentration factor during working face advance
表 9 分区情况说明 Table 9 Partition description
注:图中红色线为丁5-6煤层工作面;蓝色线为戊9-10煤层工作面;绿色线为己15煤层工作面 图 1315-21030工作面钻屑量测试位置 Fig. 13 Drilling cuttings weight test position at J15-21030 working face
图 14 不同区域钻屑量测试结果 Fig. 14 Test results of drilling cuttings weight in different areas
双重上保护层叠加开采应力分布规律
秦冰 , 叶栋林 , 石占山 , 孙维吉 , 李刚 , 张晋京