综放沿空煤柱内的自燃特征.docxVIP

综放沿空煤柱内的自燃特征 据统计，中国大部分天然采矿自然发芽率几乎来自道路，其中道路围岩的自然发芽率难以预测和控制。在采矿准备期间，靠近百叶窗周围的采矿道路周围的煤体和炉柱将被压缩，从而产生严重的扎架现象。这在很大程度上取决于采矿准备期间总暴露的煤体的氧气强度和自然发芽率。了解完井煤场周围总暴露的氧气和自然发芽率，有助于评估火灾风险带的深度和可能发生的位置，并使用适当的技术手动干预和控制。例如，将阻滞剂注入煤炭柱中，以防止煤炭氧化和自发释放。在这里，使用数值法模拟分析了疏松煤体内的漏风运动、氧浓度分布和加热过程。 1 边界条件及物理模型的确定 取单位长度的巷道段为研究对象,煤层顶、底板作为不透气边界(图1),模型几何边界的右侧边为采空区煤壁,并将压酥的煤柱体壁近似为一直线边坡,左侧边为采放的煤层,以估算的压酥范围为界,另一边界是巷道的3壁.松散煤体中漏风风流移动为低压差、低速度的稳定渗流,视其为不可压缩气体,根据质量守恒定律和达西定律,以压力p表示的漏风渗流连续性微分方程及边界条件为 式中,p为采空区风压,Pa;k为松散煤体的渗透系数,m2/(Pa·h);Ω为计算区域;Γ为计算边界,Γ=Γ1∪Γ2,Γ1,Γ2分别代表第1,2类边界;n为漏风边界法方向;pL为巷道边界的风压,Pa;qn为沿巷道边界上的漏风量,m3/(min·m). 在同一流场中,假设氧浓度与速度场可独立求解,氧扩散-消耗变化过程的定解模型为mdcdτ+divmdcdτ+div 煤的实际耗氧速度除与氧气浓度有直接关系外,随着煤体氧化升温还会使煤的氧化能力成倍增加,为充分考虑这一规律,这里近似取煤耗氧量计算式为 当耗氧为W(O2)时的总放热量为Qs=η2(b1+βb2)1+2βW(O2)Qs=η2(b1+βb2)1+2βW(Ο2) 边界条件的给定有两种情况:煤柱两侧压差或漏风量,由现场实测确定,可用间接方法测定沿空侧煤柱的漏风强度.两种边界条件一个用于计算,另一个作为模拟结果的检验.计算方案将区域按变形程度不同,采用变渗透系数处理方法.计算时物理模型进一步简化:煤柱被整体压酥呈半松散体状态,可近似为均质各向同性体;未考虑煤体温度升高所产生热风压的影响. 上述建立的数学模型用有限元数值解法求解,计算分两步:① 根据漏风流场求速度场;② 按时间循环,每个循环中由煤耗氧生热计算氧气浓度分布,再计算温度场.该对流扩散方程在数学描述上存在一个非对称、非线性的对流项,求解计算时可能出现对流效应超强于扩散,容易造成数值解振荡失真,这里采用迎风有限元解法.用专门编制的分析此问题的程序计算,计算区域自动剖分,有限元数值结果由计算机以图形方式显示出来.在输出图中,风压等值线差距为5Pa,流函数线的流量差距为0.025m3/min. 2 煤体氧传导系数模拟 算例煤层厚度11.5m,倾角11°,沿空煤柱宽度为14m.由于大采放高度使矿压显现剧烈,煤柱整体被压酥,漏风严重,通过观测巷道沿程风量损失量推定每米巷道平均漏风量为0.171m3/(min·m).利用向煤柱打钻逐段测温并抽取气样得到煤体内部的氧浓度,反求拟合得γ0=0.037mol/(m3·h),在钻孔钻透后观测得煤柱两端压差为134.5Pa(漏风指向采空区),进而得煤柱内k=0.11m2/(Pa·h),右侧煤体内k=0.009m2/(Pa·h),并逐渐过渡,m=0.093.取氧沿速度纵向扩散系数0.070m2/h,横向扩散系数0.007m2/h,模拟时设采空区侧的风压为0,巷道壁侧边界上的氧浓度始终维持在c0,原始煤体和新风的温度为17 ℃.传热参数λs=2.13W/(m·℃),Cs=0.207 J/(kg·℃), λg=0.026 5W/(m·℃), Cg=1 207J/(kg·℃),b1=221.08kJ/mol,b2=393.12kJ/mol,η=0.8.考虑煤巷片帮和顶板离层,计算时巷道形状尺寸略加改变(图2(a)). 图2为巷道周围煤体漏风流场流态的数值计算结果,风压分布等值线与流函数线正交,煤柱内部的漏风风压的分布及其流态基本呈均匀的态势,巷道左侧和上部接顶处相对势弱,沿空侧煤柱内漏风严重,流速较高.整个煤柱体内氧的浓度分布不是均匀的,氧气浓度在7mol/m3以上,煤体处于较强的氧化作用之中,巷道左侧煤体氧浓度迅速衰减.如图3所示,首先进入自燃高温区的是沿空侧半包围区,在煤动态过程上,从假设的初始漏风氧化开始,模拟大约第32天后温度分布趋于自燃状态,煤柱内的自燃高温区在3.5～6.5m处(图3(b)),煤柱内最高温度88.7 ℃.因此,巷道沿空侧3.5～6.5m范围是煤自然发火的危险区域.煤体内的升温过程是加速的(图4). 3 自然燃放期与耗氧速度回归分析 进一步模拟煤体内高温点的升温过程与漏风量、煤氧化速度常数的关系