APP下载

震动载荷下煤体瓦斯瞬时放散特征实验研究

2022-03-14孔祥国邓岱雨王恩元季鹏飞周雨璇

西安科技大学学报 2022年1期
关键词:煤体轴向震动

孔祥国,邓岱雨,王恩元,季鹏飞,王 旭,周雨璇

(1.西安科技大学 安全科学与工程学院,陕西 西安 710054;2.西安科技大学 教育部西部矿井开采及灾害防治重点实验室,陕西 西安 710054;3.中国矿业大学 安全工程学院,江苏 徐州 221116)

0 引 言

随着煤矿开采深度的逐渐增加,地应力、瓦斯压力和瓦斯含量不断增大,采场结构越来越复杂,应力和瓦斯复合型动力灾害日趋严重及复杂[1-4]。含瓦斯煤层开采过程中,顶底板断裂、断层活化和煤体破坏,均会产生应力波扰动(动载),进一步引起瓦斯异常放散,甚至煤与瓦斯突出[5-6]。

煤体瓦斯放散是煤基质瓦斯解吸、裂隙瓦斯渗流、基质裂隙瓦斯交互作用的结果[7]。林海飞等研究了煤的吸附孔结构对瓦斯放散的影响,平均孔径越大,瓦斯扩散阻力越小,瓦斯放散初速度越大[8];刘彦伟等采用压汞法研究了软、硬煤粒孔隙结构特征的差异,分析了粒径对软硬煤瓦斯扩散速度的影响[9];张小东等探讨了粒度、表面张力对原生结构煤和构造煤的瓦斯放散初速度的影响,发现吸附孔隙增多,煤基质的表面张力增大,吸附能力增强,不利于气体的解吸[10];秦跃平等根据煤粒瓦斯解吸放散的特点,建立了煤粒瓦斯放散方程,研究了瓦斯解吸速度和瓦斯累积解吸量变化规律,与试验结果一致[11-12];基于达西定律和菲克定律,王健建立了圆柱形煤粒瓦斯放散数学模型,证实了煤粒瓦斯放散服从达西定律[13];为了量化煤层瓦斯运移能力,刘永茜等探究了煤层瓦斯运移的关键参数,发现瓦斯运移速度和衰减系数取决于煤的吸附能力和煤质组分,煤体瓦斯质量增量与运移时间成负指数关系,且随煤阶升高而增大[14]。

煤体瓦斯由基质解吸到裂隙渗流再到脱离煤体自由放散,受到应力场、瓦斯场和外在条件等多种因素影响。赵洪宝等研究了固定应力环境下瓦斯放散过程中煤样变形规律,轴向应变和横向应变均表现出明显的非线性特征,并与瓦斯放散速度呈关联性[15];张冰等开展了煤粒瓦斯非常压放散特征试验,发现放散环境压力越大,瓦斯放散速度越小[16];周宏伟等研究了采动应力下煤体渗透率演化,并认为强扰动因素是确定深部煤体渗透率的关键[17];任伟光等研究发现上保护层开采扰动力学行为使得被保护层煤体状态从弹性态变到塑性态,进而导致其渗透率显著变化[18];程远平等认为深部煤层地应力主导有效应力的变化,渗透率增加的根本原因是应力的降低[19];李铁等通过现场案例和试验证实了冲击地压会诱发煤与瓦斯突出,形成瞬时瓦斯涌出现象[20];谢广祥等研究表明深部含瓦斯煤采动扩容行为致使瓦斯压力突变是诱发动力灾害的本质特征[21]。

前述研究为完善深部煤体瓦斯运移机制奠定了基础,解决了矿井生产中一系列安全问题。煤炭开采过程中,采动扰动行为有时并未引起煤体宏观破坏,但会出现瓦斯异常放散的现象。因此,震动载荷下三维含瓦斯煤瓦斯瞬间放散特征研究显得尤为必要。

1 试验系统及过程

1.1 试验系统

自主搭建了含瓦斯煤霍普金森压杆试验系统,实现了含瓦斯煤受载过程中瓦斯放散的监测。该系统主要由轴向静载子系统、围压加载子系统、瓦斯充放子系统和动载荷加载子系统、超动态应变采集子系统(采样率40 Msps,采集开始时间-0.019 97 s,采集时间0.08 s,触发电平-0.094 V)和红外测速子系统构成,如图1所示。轴向静载主要通过固定入射杆,由缓冲装置一端加压推动透射杆前移,将试样夹在入射杆和透射杆的同心轴线上;围压通过独立的液压站向试样腔体内注入液压油加压,液压油和试样通过围压胶套进行隔离;瓦斯放散量通过接在放气管路上的流量计进行实时采集;动载荷通过压缩空气发射撞击杆向入射杆撞击以施加动载荷;超动态应变采集系统通过贴于入射杆和透射杆上的应变片采集应力波信号;红外测速子系统主要用来测试撞击杆的速度,以衡量动载荷的大小。

图1 含瓦斯煤霍普金森压杆实验系统

1.2 试验试样

所选大块煤样根据煤岩动力学测试要求,制成φ100 mm×50 mm的试样,经过打磨将试样两端平行度控制在0.02 mm公差范围内。为了减少试样的离散性,选择密度相近的试样进行试验,试验时在试样两端涂抹黄油以减少摩擦效应。

1.3 试验方案

文中主要研究不同轴向静载、围压、瓦斯压力和震动载荷组合加载下含瓦斯煤受载过程中瓦斯放散特征,分析各因素对含瓦斯煤瓦斯最大放散量、瓦斯放散持续时间的影响,试验方案见表1。

表1 试验方案

2 试验结果及分析

2.1 瓦斯放散随时间变化规律

不同轴向静载、围压、瓦斯压力和动载荷组合加载下含瓦斯煤受载过程中瓦斯放散如图2(a)~(d)所示。图2(a)呈现了瓦斯放散随着轴向静载的变化规律,随着轴向静载的增加,最大瓦斯放散量逐渐减小,说明动载荷不足以引起煤体破裂时,轴向静载有助于提高煤体完整性使其不易破坏,轴向静载越大,原始裂隙越致密,动载荷不易引起煤体裂隙扩展而导致瓦斯放散;不同轴向静载下瓦斯放散持续时间大约在2 s之内,瓦斯放散持续时间与轴向静载关系不太明显,但是可以反映瓦斯放散的瞬时性,解释了震动诱导含瓦斯煤瓦斯异常放散的瞬时特性。

图2(b)展现了瓦斯放散随围压的关系,最大瓦斯放散量随着围压的增加呈减小的趋势,围压的作用使得煤体更加致密,而且围压限制煤体变形,围压越大,煤体新的裂隙越难形成,瓦斯不易放散,但是在动载荷作用瞬间仍然有瓦斯瞬时增大时刻,说明动载荷会引起煤体微损伤,但是在围压的作用下,新生裂隙又闭合,瓦斯放散持续时间在0~4 s之内,围压越小,瓦斯放散持续时间相对越长。

图2(c)反映了含瓦斯煤瓦斯放散随煤体内部瓦斯压力的关系,瓦斯压力越大,冲击瞬间瓦斯放散量越大,主要是瓦斯压力越大,放散时瓦斯压力梯度越大,促进瓦斯流动;瓦斯放散持续时间控制在2 s以内,放散持续时间随着瓦斯压力增大而相对延长,说明高瓦斯煤层受动载荷影响将更加危险。

图2(d)刻画了瓦斯放散与动载荷的关系,动载荷越大,最大瓦斯放散量越大,说明动载荷瞬间冲击作用更强,对煤体微损伤影响较为明显,不同动载荷下含瓦斯煤瓦斯放散持续时间大约在0~4 s之内,仍然反映了震动载荷作用下瓦斯放散的瞬时性。

图2 震动载荷下含瓦斯煤瓦斯放散规律

2.2 最大瓦斯放散速度与各因素关系

瓦斯放散是一个时间过程,动载荷作用于煤体之后,煤体内部发生微损伤,当动载荷较大时,裂隙扩展和贯通导致煤体宏观破裂,因此瓦斯流动亦随着时间而不断演化。在震动载荷下,含瓦斯煤最大瓦斯放散速度随轴向静载、围压、瓦斯压力和震动载荷冲击速度变化如图3(a)~(d)所示。瓦斯最大放散速度随轴向静载呈线性递减,拟合度高达94.7%,震动载荷不足以使煤体宏观破裂,轴向静载提高了其动态抗压强度,瓦斯不易放散,相应放散速度较小。由图3(b)可得,含瓦斯煤最大瓦斯放散速度随围压呈指数衰减,拟合度高达92.84%,围压不利于瓦斯放散,瓦斯放散瞬间最大瓦斯放散速度相应较小。由图3(c)和(d)可得,含瓦斯煤最大瓦斯放散速度随瓦斯压力和动载荷冲击速度的增加均呈现“线性增加”的趋势,瓦斯压力梯度有助于瓦斯流动,震动载荷冲击速度对煤体形成微损伤,增加了裂隙瓦斯流动通道。

图3 震动载荷下含瓦斯煤瓦斯放散速度与各因素关系

2.3 瓦斯放散率与各因素关系

文中类比渗透率稳态法,测定含瓦斯煤瞬时放散率,控制煤体两侧瓦斯压力不变,来监测通过煤体的瓦斯流量,通过式(1)计算煤体瓦斯瞬时放散率,由于本实验震动载荷作用时间较短,因此,只考虑震动载荷冲击瞬间的最大放散率,称为瞬时放散率。

(1)

式中k′为瞬时放散率,m2;q为瞬时瓦斯流量,m3/s;p0为大气压力,MPa;μ为瓦斯动力粘度,1.84×10-11MPa·s;L0为试样长度,m;A0为试样横截面积,m2;p2为进气端瓦斯压力,MPa;p1为试样出气端瓦斯压力,MPa。由于试样出气端为自由端与大气联通,因此这里p1=p0。

不同轴向静载、围压、瓦斯压力和震动载荷冲击速度组合加载下,含瓦斯煤瓦斯放散率如图4(a)~(d)所示。类似于瓦斯放散速度,轴向静载的作用使得煤体不易损伤,瓦斯放散率随着轴向静载的增加而减小;围压使煤体更加致密并限制变形造成瓦斯放散率降低;虽然瓦斯压力有助于体高瓦斯放散速度,但瓦斯压力平方差大于瓦斯放散量增加的程度,造成瓦斯放散率随瓦斯压力增加而非线性减小;震动载荷冲击速度越大,动载荷携带能量越多,对煤体损伤越严重,因此瓦斯放散率将随着动载荷的增大而增加。

图4 震动载荷下含瓦斯煤瓦斯放散率

3 结 论

1)含瓦斯煤在震动载荷作用瞬间瓦斯放散表现出瞬时性特征,瓦斯放散持续时间在0~4 s内;震动载荷会诱导煤体微裂隙增加而引起瓦斯异常放散,但围压作用很快又使煤体微裂隙闭合。

2)当震动载荷不足以使煤体破裂时,轴向静载有助于提高煤体动态抗压强度,裂隙难以形成,最大瓦斯放散速度将随轴向静载增加而减小;围压的限制变形不利于瓦斯流动造成最大瓦斯放散速度随围压的增加而减小;瓦斯压力梯度促进瓦斯流动,震动载荷作用过程中最大瓦斯放散速度随着瓦斯压力增大而增加;震动载荷冲击速度越大,携带能量越多,煤体微损伤越严重,最大瓦斯放散速度将随震动载荷冲击速度增加而增大。

3)含瓦斯煤瞬时放散率随着轴向静载、围压的增大而线性减小,降低幅度分别达到44.80%和89.90%;随着瓦斯压力的增加呈非线性减小,主要是由于瓦斯压力的平方差增加程度大于瓦斯放散量的增加程度;随着震动载荷冲击速度增加基本呈增加的趋势,增加幅度达90.89%。

猜你喜欢

煤体轴向震动
逆断层区域煤体应力分布规律
松散煤体声发射特征与损伤本构模型
新型轴向磁通转子错角斜极SRM研究
受载突出煤体的力学状态演变及破坏倾向性
千分尺轴向窜动和径向摆动检定装置的研制
基于串联刚度模型的涡轮泵轴向力计算方法
画与理
确定性地震动空间差异对重力坝地震响应影响研究
双楔式闸阀阀杆轴向力的计算
双楔式闸阀阀杆轴向力的计算