高温热源表面点燃甲烷和煤尘实验研究

王毅, 王海燕, 张雷, 胡浪

王毅,王海燕,张雷,等.高温热源表面点燃甲烷和煤尘实验研究[J].工矿自动化,2021,47(9):85-90.. DOI: 10.13272/j.issn.1671-251x.17823
引用本文: 王毅,王海燕,张雷,等.高温热源表面点燃甲烷和煤尘实验研究[J].工矿自动化,2021,47(9):85-90.. DOI: 10.13272/j.issn.1671-251x.17823
WANG Yi, WANG Haiyan, ZHANG Lei, HU Lang. Experimental study on ignition of methane and coal dust on the surface of high temperature heat source[J]. Journal of Mine Automation, 2021, 47(9): 85-90. DOI: 10.13272/j.issn.1671-251x.17823
Citation: WANG Yi, WANG Haiyan, ZHANG Lei, HU Lang. Experimental study on ignition of methane and coal dust on the surface of high temperature heat source[J]. Journal of Mine Automation, 2021, 47(9): 85-90. DOI: 10.13272/j.issn.1671-251x.17823

高温热源表面点燃甲烷和煤尘实验研究

基金项目: 

国家重点研发计划资助项目(2018YFC0808100)

国家自然科学基金项目(51874313)

详细信息
  • 中图分类号: TD714.5

Experimental study on ignition of methane and coal dust on the surface of high temperature heat source

  • 摘要: 煤矿采空区煤自燃高温区域可能导致甲烷/煤尘爆炸,有必要开展煤自燃诱发瓦斯/煤尘爆炸的研究。目前基于电火花诱发方式对瓦斯/煤尘爆炸特性的研究不能充分反映采空区内煤自燃诱发瓦斯/煤尘爆炸的特性,针对该问题,利用自主研制的40 L气体爆炸设备研究在高温热源表面点火方式下甲烷/煤尘爆炸特性。结果表明:① 通过高温热源表面点火方式得到的甲烷/煤尘爆炸过程与甲烷/空气爆炸过程一致,均可分为3个阶段。② 甲烷/煤尘的3个爆炸参数(爆炸压力、压力上升速率和爆炸温度)随甲烷浓度的变化规律一致,即先随着甲烷浓度的增加而升高,达到峰值后随甲烷浓度的增加而降低;当甲烷浓度小于转折浓度时,甲烷/煤尘的3个爆炸参数分别大于甲烷/空气相对应的爆炸参数;当甲烷浓度大于转折浓度时,甲烷/煤尘的3个爆炸参数分别小于甲烷/空气相对应的爆炸参数;由于煤尘的加入,甲烷/空气的3个爆炸参数的转折体积分数由115%分别降至95%,105%和95%,峰值分别增加了08%,69%和08%,且规律性下降。③ 由于点火方式不同,爆炸温度的变化规律不同:电火花引爆甲烷/煤尘后爆炸温度在下降阶段(Ⅲ阶段)再次出现升高的现象;高温热源表面引爆甲烷/煤尘后爆炸温度变化规律与未加入煤尘时保持一致。④ 对煤尘点燃温度预测函数进行了修正,修正后的拟合优度提高了2005%;结合实验数据得到了煤尘在高温热源表面点火方式下的点燃温度预测函数。
    Abstract: The high temperature area of coal spontaneous combustion in coal mine goaf may cause methane/coal dust explosion. Therefore, it is necessary to carry out research on coal spontaneous combustion induced methane/coal dust explosion. The current research on the characteristics of methane/coal dust explosion based on the electric spark induced method cannot fully reflect the characteristics of methane/coal dust explosion induced by coal spontaneous combustion in the goaf. In order to solve this problem, the self-developed 40 L gas explosion equipment is used to study methane/coal dust explosion characteristics under the high temperature surface ignition of heat source. The results are showed as follows. ① The methane/coal dust explosion process obtained by the surface ignition method of the high-temperature heat source is consistent with that of methane/air explosion, which can be divided into 3 stages. ② The three explosion parameters of methane/coal dust (explosion pressure, pressure rise rate and explosion temperature) are consistent with the change law of methane concentration, that is, the parameters first increase with the increase of methane concentration, reach a peak and then decrease with the increase of methane concentration. When the methane concentration is less than the turning concentration, the three explosion parameters of methane/coal dust are greater than the corresponding explosion parameters of methane/air respectively. When the methane concentration is greater than the turning concentration, the three explosion parameters of methane/coal dust are less than the corresponding explosion parameters of methane/air respectively. Due to the addition of coal dust, the turning concentrations of the three explosion parameters of methane/air decrease from 11.5% to 9.5%, 10.5% and 9.5% respectively, and the peak values increase by 0.8%, 6.9% and 0.8% respectively, and decrease regularly. ③ Due to the different ignition methods, the explosion temperature change law is different. The explosion temperature rises again in the falling stage (stage III) after the methane/coal dust is detonated by the electric spark. The explosion temperature change law after detonation of methane/coal dust on the surface of high-temperature heat source remains the same as that before the addition of coal dust. ④ The coal dust ignition temperature prediction function is corrected, and the goodness of fit is improved by 20.05% after correction. Combined with the experimental data, the ignition temperature prediction function of coal dust in the ignition mode of the high temperature heat source surface is obtained.
  • 我国西部矿区具有埋藏浅、煤层厚、地质条件简单等特点。为提高工作面开采效率,缓解接续紧张局面,回采巷道多采用双巷布置方式,双巷间留设一定宽度的区段煤柱,用于维护回采巷道的稳定[1-3]。煤柱宽度过小则不能有效支撑顶板,在采动影响产生的高集中应力作用下破坏失稳[4-5];煤柱宽度过大,不仅造成煤炭资源浪费,还有可能引发冲击地压等动力灾害[6-7]

    近年来,学者们针对区段煤柱合理宽度设计开展了大量研究,主要运用极限平衡理论计算煤柱塑性区宽度,在考虑一定安全系数的情况下,预留一定宽度的弹性核区,将两者宽度之和作为煤柱合理宽度[8-10]。张念超等[11]在统一强度理论的基础上,提出了区段煤柱合理宽度理论计算方法。王琦等[12]提出了最大临界尺寸的概念,认为区段煤柱合理宽度应小于最大临界尺寸。姚强岭等[13]考虑采空区积水弱化作用,提出了煤矿地下水库煤柱坝体塑性区宽度计算公式。周光华等[14]针对特定的近距离倾斜煤层群采掘地质条件,研究了煤矿地下水库煤柱坝体的合理尺寸及其稳定性问题。研究表明,区段煤柱合理宽度的确定不仅受到工作面采动的影响[15-16],而且与区段煤柱一侧采空区积水弱化作用有关[17-19]。采空区积水与煤岩作用会弱化区段煤柱强度而引起煤柱逐渐破坏和失效,水岩作用是区段煤柱合理宽度设计中必须考虑的关键因素。

    本文考虑内蒙古鄂尔多斯新街矿区某矿31采区采空区积水对煤体强度的弱化作用,通过单轴压缩实验获得了煤体强度弱化系数,分析了水岩作用下区段煤柱塑性区宽度特征;利用FLAC3D数值计算模型实现煤体强度的动态弱化,并在此基础上开展了不同宽度煤柱稳定性数值模拟研究,可为相似工程条件下煤柱宽度设计提供参考。

    新街矿区某矿31301工作面位于31采区(主采3−1煤层),为31采区的首采工作面,走向长度为2 503.74 m,宽度为300.58 m,面积为752 573 m2。煤层平均厚度为6 m,倾角为1~3°。直接顶为厚5.96 m的砂质泥岩,基本顶为厚18.7 m的中粒砂岩,直接底为厚1.23 m的泥岩,基本底为厚9.2 m的砂质泥岩。

    31采区煤层底板总体呈东北部高、西南部低的趋势。31301采空区积水高度超过3 m,部分区域达7 m以上,采空区积水对3−1煤层始终存在强度弱化作用。目前,根据矿井的接续安排,需要在31301采空区一侧布置33采区33301工作面,如图1所示。

    图  1  工作面布置
    Figure  1.  Layout of working face

    在31采区同一位置选择大小适中(300 mm×300 mm×300 mm)、形状规则、层理分布均匀且表面无明显裂纹的煤块,将现场选择的煤块运回实验室后加工处理成标准试样,用于单轴压缩实验,结果如图2所示。3−1煤层试样在干燥、饱和条件下的抗压强度分别为18.45,7.89 MPa,弹性模量分别为2.1,1.1 GPa。根据实验结果,确定煤层饱和条件下强度弱化系数为0.43,弹性模量弱化系数为0.52。

    图  2  煤样应力−应变曲线及破坏形态
    Figure  2.  Stress-strain curves and failure mode of coal sample

    区段煤柱服务期间,先后经历掘进影响、31301工作面一次采动影响和33301工作面二次采动影响阶段,其中,31301工作面回采结束后,由于侧向顶板破断,导致煤柱侧产生宽度为x0的塑性破坏区,在33301工作面回采期间,31301采空区积水对煤柱强度进一步弱化,产生新的塑性区(宽度为x1),根据极限平衡理论,可得[13]

    $$ {x_1} = \frac{{{M}}}{{2\lambda f}}\ln \left[ {\frac{{{{{K}}_{\text{1}}}\gamma H}}{{\lambda p + {\sigma _{\rm{r}}}}} - \frac{{{V_{\rm{m}}}{S_{\rm{g}}}\left( {1 - {\exp\left({\dfrac{{2\lambda f}}{{{S_{\rm{m}}}{S_{\rm{g}}}}}({\sigma _{\rm{c}}} - {\sigma _{\rm{r}}})}\right)}} \right)}}{{2\lambda f\left( {\lambda p + {\sigma _{\rm{r}}}} \right)}}} \right] $$ (1)

    式中:M为煤层开采厚度;λ为应力系数,λ=(1+sin φ)/(1−sin φ),φ为内摩擦角; f为煤层与顶底板的摩擦因数;K1为煤柱积水侧应力集中系数;γ为覆岩容重;H为煤层埋深;p为积水侧静水压力;σr为煤体残余强度;Vm为煤体软化模量;Sg为塑性区煤体应变梯度;σc为弹性阶段煤体的单轴抗压强度。

    由式(1)可知,区段煤柱积水侧塑性区宽度随煤体强度弱化程度的增加而扩大。31301工作面具体参数:M=6 m,λ=2.5,f=0.56,K1=3,γ=25 kN/m3H=380 m,p=0.06 MPa,Vm=800 MPa,Sg=0.08,σr=1.81 MPa,σc=7.89 MPa。将以上参数代入式(1),可得x1=6.42 m。

    由于33301工作面为33采区首采工作面,区段煤柱需要长期保持稳定,以防止31301采空区积水溃入,区段煤柱保持稳定的基本条件:区段煤柱两侧产生塑性变形后,在煤柱中央存在一定宽度的弹性核区。保持煤柱稳定的最小宽度为

    $$ B = 2{x_0} + kM + {x_1} $$ (2)

    式中k为安全系数,取3。

    根据31301工作面现场测试结果,x0取14.6 m。通过理论计算,确定31301工作面和33301工作面间的区段煤柱宽度为53.62 m。

    根据实际条件建立不同宽度煤柱FLAC3D数值计算模型,模型大小为880 m×400 m×190 m(长×宽×高),在模型的顶面施加6.8 MPa载荷来代替上方未建立岩层的重力作用,模型单元体材料采用Mohr-Coulomb强度准则,煤岩体物理力学参数见表1

    表  1  数值计算模型煤岩体物理力学参数
    Table  1.  Physical and mechanical parameters of coal and rock mass in numerical calculation model
    岩层厚度/m弹性模
    量/GPa
    泊松比黏聚
    力/MPa
    内摩擦
    角/(°)
    抗拉强
    度/MPa
    粉砂岩343.510.175.1634.581.500
    细粒砂岩81.610.174.3731.830.510
    粉砂岩345.860.155.6636.000.890
    中粒砂岩63.450.094.2431.340.480
    砂质泥岩265.350.196.2637.441.080
    泥岩26.300.194.9233.830.670
    60.990.213.0526.090.610
    煤(浸水)0.490.210.429.530.004
    泥岩66.750.185.8031.770.680
    粉砂岩264.690.176.4633.580.860
    下载: 导出CSV 
    | 显示表格

    分别取煤柱宽度为50,56,62,68,74 m进行模拟,通过编写Fish语言,对煤柱塑性区单元强度进行弱化,并在数值计算过程中进行迭代,实现模拟水岩作用的动态过程,即随着33301工作面推进,逐次弱化位于31301采空区积水侧的煤柱,直至33301工作面回采结束。 31301及33301工作面模拟走向长度均为300 m,推进长度为100 m。

    31301工作面回采后煤柱塑性区分布如图3所示。可看出当31301工作面回采80 m左右时,煤柱塑性区范围达到最大,随着工作面继续推进,煤柱塑性区范围变化不大,当工作面推进100 m时,塑性区宽度约为12 m。

    图  3  31301工作面回采后煤柱塑性区分布
    Figure  3.  Distribution of plastic zone after 31301 working face mining

    煤柱宽度为50,56 m时塑性区分布如图4所示。由图4(a)可知,当33301工作面回采70 m时,50 m宽煤柱距离开切眼25 m处塑性区发生贯通,同时,该处直接顶也产生了塑性区的贯通,直接底具有塑性区贯通的危险,33301工作面回采结束后,整个煤柱已大面积发生破坏,失去承载能力。由图4(b)可知,当33301工作面回采90 m时,56 m宽煤柱距离开切眼65 m处塑性区发生贯通,同时,该处直接顶也产生了塑性区的贯通,直接底具有塑性区贯通的危险,33301工作面回采结束后,煤柱大部分发生破坏,承载能力极低。

    图  4  煤柱宽度为50,56 m时塑性区分布
    Figure  4.  Distribution of plastic zone when the width of coal pillar is 50 m and 56 m

    煤柱宽度为62 m时塑性区分布如图5所示。可看出当33301工作面回采结束时,若不考虑采空区积水弱化作用,煤柱未产生塑性区贯通,中央保留约14 m宽的弹性核区;若考虑采空区积水弱化作用,煤柱中部距开切眼45 m处塑性区产生贯通,同时,该处直接顶和直接底具有塑性区贯通的危险性较大,煤柱仅保留较弱的承载能力。

    图  5  煤柱宽度为62 m时塑性区分布
    Figure  5.  Distribution of plastic zone when the width of coal pillar is 62 m

    煤柱宽度为68 m时塑性区分布如图6所示。可看出当33301工作面回采结束时,68 m宽煤柱采空区积水侧塑性区宽度为24 m,33301工作面回采全过程中煤柱及直接顶、直接底均无贯通危险,完整弹性核区宽度为32 m,33301工作面回采结束后,煤柱虽具备一定的承载能力,但采空区积水弱化作用对煤柱强度仍有一定影响,煤柱中部塑性区有凸出发育的趋势。

    图  6  煤柱宽度为68 m时塑性区分布
    Figure  6.  Distribution of plastic zone when the width of coal pillar is 68 m

    煤柱宽度为74 m时塑性区分布如图7所示。可看出33301工作面回采结束后,74 m宽煤柱采空区积水侧距开切眼50 m处塑性区宽度为24 m,弹性核区宽度为40 m,采空区积水弱化作用对煤柱稳定性的影响较弱,煤柱整体稳定性强。

    图  7  煤柱宽度为74 m时塑性区分布
    Figure  7.  Distribution of plastic zone when the width of coal pillar is 74 m

    综上所述,在弹性核区处于相对较高的集中应力承载状态时,采空区积水弱化作用对煤柱的稳定性影响显著,易造成煤柱及直接顶、直接底的贯通型渗漏,并大幅降低了整体承载性能。煤柱宽度在68 m以下时,煤柱均会发生不同程度的贯通破坏,且直接顶、直接底也存在贯通破坏的危险,采空区积水渗漏途径不会仅局限于煤壁,加之煤柱的承载能力很弱,稳定性得不到可靠保障;68 m宽煤柱虽然弹性核区尚可,但煤柱中部塑性区仍有凸出发育的趋势,采空区积水长期作用下存在一定的贯通渗漏隐患;74 m宽煤柱稳定性更好,采空区积水弱化作用不显著。

    不同宽度煤柱的应力分布云图、垂直应力集中系数分布曲线分别如图8图9所示。

    图  8  煤柱应力分布云图
    Figure  8.  Nephogram of coal pillar stress distribution
    图  9  垂直应力集中系数分布曲线
    Figure  9.  Vertical stress concentration coefficient distribution curve

    不同宽度煤柱在33301工作面回风巷侧垂直应力集中系数远高于采空区积水侧,最高达2.89;随着煤柱宽度增加,垂直应力集中系数不断下降,平均降幅为4.2%;50~68 m宽煤柱最高垂直应力集中位置距离33301工作面回风巷右帮8 m左右;煤柱宽度为74 m时,最高垂直应力集中系数大幅下降至2.08,降幅为15.5%;煤柱采空区积水侧除62 m宽煤柱垂直应力集中程度较低外(由距开切眼45 m处小范围贯通破坏影响所致),其他宽度煤柱垂直应力集中系数稳定在1.9左右,且随着煤柱宽度的增加,垂直应力集中位置逐渐远离33301工作面回风巷侧。

    随着煤柱宽度的增大,其承载能力增强,两侧形成的应力集中差异性减小,弹性核区较低应力集中区域范围逐渐增大。68 m宽煤柱双侧最高垂直应力集中系数相差0.45,弹性核区最低垂直应力集中系数为1.63,垂直应力分布均匀;74 m宽煤柱双侧最高垂直应力集中系数相差0.15,弹性核区最低垂直应力集中系数为1.47,垂直应力分布均匀。煤柱宽度较小时,采空区积水弱化作用对较高应力集中的弹性核区具有更强的破坏能力;随着煤柱宽度的增大,弹性核区应力集中程度降低,采空区积水侧垂直应力低于原岩垂直应力的区域范围则有所增大,煤柱两侧应力集中分布趋于均匀,采空区积水弱化作用对弹性核区的影响不再显著。综合以上数值模拟结果,确定煤柱宽度应不小于68 m。

    根据理论计算和数值模拟结果,确定31301采空区与33301工作面之间的煤柱宽度为70 m。在33301回风巷内布置测站,对二次采动影响条件下矿压显现规律进行监测,结果如图10所示。可看出70 m留设煤柱可以有效承载顶板压力,33301回风巷顶板最大变形量为103 mm,两帮最大移近量为39 mm,巷道围岩变形小,锚索受力稳定,为矿井安全高效生产提供了保障。

    图  10  33301回风巷矿压显现规律
    Figure  10.  Mine pressure behavior law of 33301 return airway

    (1) 水岩作用对煤体强度弱化产生显著影响,区段煤柱积水侧塑性区宽度随煤体强度弱化程度的增加而扩大。基于区段煤柱保持稳定的基本条件,通过理论计算确定了区段煤柱的合理宽度为53.62 m。

    (2) 数值模拟结果表明:煤柱宽度小于62 m时,采空区积水弱化作用对较高应力集中的弹性核区具有更强的破坏能力;随着煤柱宽度的增大,弹性核区应力集中程度降低,采空区积水侧垂直应力低于原岩垂直应力的区域范围则有所增大,煤柱两侧应力集中分布趋于均匀,采空区积水弱化作用对弹性核区的影响不再显著。

    (3) 综合理论计算与数值模拟结果,确定了区段煤柱宽度为70 m。工程应用结果表明,70 m宽留设煤柱可以有效承载顶板压力,巷道围岩变形小,锚索受力稳定,保障了矿井安全生产。

  • 期刊类型引用(6)

    1. 石龙飞. 工作面停采煤柱宽度的确定与实践. 煤炭与化工. 2024(04): 14-17+24 . 百度学术
    2. 刘华锋. 干湿循环处理后煤在差速循环荷载作用下的力学特征. 工矿自动化. 2024(11): 152-160 . 本站查看
    3. 屈世甲,杨欢. 巷道内瓦斯爆炸状态下人工坝体的力学响应研究. 工矿自动化. 2023(09): 132-139 . 本站查看
    4. 贾尚伟,徐壮飞,王海洋,贺富领. 特厚煤层大采高综放开采区段煤柱留设宽度研究. 河南科学. 2023(10): 1446-1453 . 百度学术
    5. 高振东,郑文翔. 煤柱尺寸对沿空巷道稳定性影响研究. 内蒙古煤炭经济. 2023(19): 7-9 . 百度学术
    6. 胡平,孙爱东,温国惠,朱建平,连平小,段会强,唐共龙,戴永帅,史振,杨超,常超. 东胜矿区浅埋煤层区段煤柱宽度优化设计. 煤炭技术. 2023(12): 56-60 . 百度学术

    其他类型引用(2)

计量
  • 文章访问数:  159
  • HTML全文浏览量:  18
  • PDF下载量:  14
  • 被引次数: 8
出版历程
  • 发布日期:  2021-09-19
  • 刊出日期:  2021-12-28

目录

/

返回文章
返回