动态推进U型采场煤自燃多场耦合数值研究.pdf
《动态推进U型采场煤自燃多场耦合数值研究.pdf》由会员分享,可在线阅读,更多相关《动态推进U型采场煤自燃多场耦合数值研究.pdf(8页珍藏版)》请在文库网上搜索。
1、92eMeouscombustioninnamicSafetvinCoalMin54(10):92-99shaadvancel移动扫码阅读QIN Jianyun,SONGShuangWANGYongletal.Multi-fieldnumericalstudyyofcoal spontan-2023,54(10):92-99.秦剑云,宋双林,王永敬推进灯型米场煤数研煤矿安全Safetyin Coal MinesOct.20232023年10 月No.10Vo1.54煤矿岁全第54 卷第10 期DOI:10.13347/ki.mkaq.2023.10.013动态推进U型采场煤自燃多场耦合数值研究
2、秦剑云,宋双林2 3,王永敬2 3,李小超,党龙2 3,朱鹏,李世豪2.3(1.库车县榆树岭煤矿有限责任公司,新疆库车8 4 2 0 9 9;2.中煤科工集团沈阳研究院有限公司,辽宁抚顺11312 2;3.煤矿安全技术国家重点实验室,辽宁抚顺11312 2)摘要:针对动态推进过程中U型采场的煤自然发火问题,建立包含采场进回风巷道、工作面与采空区的多场耦合数学模型,并通过现场实测数据验证了模型的合理性。模拟分析了采空区内温度场、氧气体积分数场随工作面动态推进的变化规律,并进一步研究了围岩温度、风流温度、工作面推进速度对采空区煤氧化自热的影响。结果表明:工作面推进过程中,采空区的进风侧温度场形成“
3、机翼断面”形状的高温区域拖尾现象,并且高温区域在采空区的堆积压实和围岩散热作用下温度会逐渐降低;停采后采空区内继续升温但高温区域向工作面迁移;不同通风温度对采空区最高温度的影响,在开采的初期几乎无影响,而在开采的后期,通风温度越高,采空区内最高温度越高;工作面推进速度越快,采空区内最高温度越低。关键词:U型采场;煤自燃;动态推进;浓度场;温度场中图分类号:TD75+2.2文献标志码:A文章编号:10 0 3-4 9 6 X(2023)10-0092-08Iulti-field coupling numerical study of coal spontaneous combustion in
4、U-shaped stope dynamic advancQIN Jianyun,SONG Shuanglin?3,WANG Yongjing,LI Xiaochao,DANG Long3,ZHU Peng,LI Shihao?2,3(1.Kuqa County Yushuling Coal Mine Co.,Ltd.,Kuqa 842099,China;2.China Coal Technology&Engineering Group ShenyangResearch Institute,Fushun 113122,China;3.State Key Laboratory of Coal M
5、ine Safety Technology,Fushun 113122,China)Abstract:In order to solve the problem of spontaneous combustion of coal in U-shaped stope during dynamic advance,a multiphys-ical field coupled mathematical model including the inlet and return airway,working face and goaf is established,and the reasonable-
6、ness of the model is verified by field measurement data.The simulation analyzed the law of temperature field and oxygen concentra-tion field in the goaf with the dynamic advance of working face,and further studied the influence of surrounding rock temperature,air flow temperature and working face ad
7、vance speed on the self-heating of coal oxidation in the mining area.The results indicatethat the temperature field on the inlet side of the goaf forms the phenomenon of high temperature region tailing in the shape of“wingsection in the process of advancing the working face,and the temperature of th
8、e high-temperature area gradually decreases underthe effect of compaction of the goaf and heat dissipation of the surrounding rock;the temperature in the goaf continues to rise afterthe stoppage of mining,but the high-temperature area migrates to the working face;in the early stage of mining,differe
9、nt ventilationtemperature has almost no effect on the maximum temperature in the mining area,while in the later stage of mining,the higher theventilation temperature,the higher the maximum temperature in the mining area;the faster the working face advances,the lower themaximum temperature intheminin
10、g area.Key words:U-shaped stope;spontaneous combustion of coal;dynamic advance;concentration field;temperature field收稿日期:2 0 2 2-0 9-2 3责任编辑:陈洋作者简介:秦剑云(19 8 3),男,山西长治人,工程师,硕士,从事瓦斯综合治理与瓦斯利用方面的工作。E-mail:154 4 9 4 0 9 9 q q.c o m93SafetyinCoalMines2023年10 月Oct.2023煤砺全第54 卷Vol.54第10 期No.10煤炭作为我国主体能源,其安全
11、供给直接关系到国家能源安全!。煤矿火灾是威胁煤矿安全生产的主要原因之一,已严重威胁到煤矿企业的安全生产,其中煤炭自燃是引发煤矿火灾的最主要原因 2-3,采动影响下的垮落岩体和遗煤在地下形成的多孔松散采空区是瓦斯与煤自燃致灾的重灾区,井下采空区遗煤自燃占井下火灾事故的60%4针对采空区的自然发火问题,国内外学者开展了一系列的研究,包括防灭火措施、自燃规律的数值模拟及实验分析等。在数值模拟方面,前期的研究以研究风流在采场的稳态分布规律为主 5-7 ,根据求解的漏风速度对自燃区域进行大致区分。目前的研究多是针对采空区自然发火的多场耦合模拟研究 8-,将流场、温度场、浓度场、应力场等进行耦合,从矿井尺
12、度上研究煤自燃高温区域的发生及发展的规律10-12 。由于采空区的物理尺寸和物理特性都会随着工作面的推进而改变,一些学者通过引人移动坐标系来处理物理边界的变化 8.13,或者利用CFD软件的动网格技术 14 ,得到孔隙度在时间和空间上的非均质分布,体现高温区域在小范围推进距离下的动态移动规律,模拟结果更符合实际。文献 15-17 通过COMSOLMultiphysics与MATLAB联用,在生成新的采空区网格后,调用上一时间步的场信息作为当前时间步的初值,实现模拟的物理量在时间和空间上的连续性,较大程度呈现了采空区煤自燃升温过程,但这些研究未将工作面与采空区的流动、传热及传质问题联合求解,从而
13、也缺少围岩温度、通风温度等因素对采空区煤自燃升温过程影响的研究。为此,考虑采空区动态演化过程,以巷道出人口作为边界,将工作面与采空区的物理场联合求解,模拟工作面动态推进下采空区煤岩裂隙场中煤自燃灾害的发生、发展及演化规律;通过对采空区遗煤氧化升温情况的预测预报,可确定采空区煤自燃的危险性和高温热源的位置,进而为采空区这一隐蔽火区煤自燃的治理提供及时有效的手段。1模型构建U型采场可分为几个区域,即煤体固相区、采空多孔介质区和采煤作业区。煤体固相区包括面间煤柱、边界煤柱、待回采实体煤、顶板岩层和底板岩层;采煤作业区包括进风巷、工作面和回风巷。采空区二维示意图如图1。围岩回风巷工实体媒面开一(x,y
14、)切眼围岩XXo图1采空区二维示意图Fig.1Two-dimensional schematic diagram ofthe mining area1.1基本假设1)因采空区长度和宽度远大于煤层开挖引起的上覆煤岩层垮落高度,将采空区视为高度上平均的二维非均质多孔介质。2)假设气体为理想气体,并忽略水蒸气对煤自热过程的影响1.2采空区孔隙率及渗透率的动态演化随着工作面向前推进,自由垮落带的砰石、遗煤在重力和矿压作用下会被逐渐压实,垮落碎胀系数在压实过程中不断减小,其数值满足式(1)18 ;Kp(x,y)=Kp.min+(Kp,max-Kp,min)exp(-xd,1-exp(-$i,d,)(1)
15、式中:Kp(x,y)为采空区跨落碎胀系数;Kp.max为初始垮落碎胀系数;Kp,min为垮落压实后的碎胀系数;x、分别为采空区碎胀系数倾向和走向方向的衰减率,1/m;d x、d,分别为采空区任意点(,)到采空区某一固定围岩边界和工作面边界的距离,m;$为控制“o”型圈模型分布形态的调整系数,si1。如果工作面推进速度为u,则随着时间推进,采空区内任意一点到采空区上下边界和工作面的距离分别可表示为:dx=x-Xo+utd,=0.5L-ly-0.5Ll(2)式中:u为工作面推进速度,m/s;t 为工作面回采时间,d;L为工作面长度,m;x o 为工作面初始位置,m。垮落采空区的高度H和孔隙率的表达
16、式分别为:94Safetyin Coal Mines2023年10 月Oct.2023第54 卷第10 期No.10Vo1.54煤矿发全MK,(x,y)H=(3)Kp.max-118b=1(4)K,(x,y)式中:M为工作面采高,m;8 为采空区孔隙率。根据Ergun的实验研究 19 ,采空区渗透率K可以根据孔隙率、体积平均粒径计算得到,计算公式为:K=(5)150(1-8b)2式中:d为采空区多孔介质的平均颗粒直径,m。1.3采空区流动及传热方程采空区内渗透率具有明显的非均质特征,固壁处的渗透率和采空区深处的孔隙率相差1个数量级,研究表明在靠近工作面的采空区内具有较大的空气流速,空气处于湍流
17、状态,而采空区内部的空气流速较低,空气处于层流或者蠕动流状态,因此使用Darcy-Brinkman方程描述采空区内气体的流动 2 0 。对于采空区内空气与煤岩体之间的换热过程,采用局部热平衡假设,即假定局部流体温度与煤岩体温度一致。采空区内的流动、传热及传质方程 2 1-2 3 为:0(pf8b)V.(prpV)+=0tV1CFPPf=-Vp+兰2 Vbt8bKO.5Ho(pc)V.VT=V.(kmVT)(Qb-Qo)moteff,fHHo(ci)+V.Vc;-V(Dp,i+De,)Vci=R;otH(6)式中:Pr为空气密度,kg/m;V为气流速度,m/s;IVI=Vu+u 2;u、v 分别
18、为走向和倾向方向上1.75的速度,m/s;t 为时间;CF=Q.为造煤150%与顶底板接触损失的热量;T为温度,K;(o c)m为煤体有效热容,J/(mK);k m为煤体有效导热系数,W/(mK);H。为采空区遗煤厚度,m;D p.为单煤颗粒间气体弥散扩散系数,m/s;D e 为单煤颗粒间气体有效扩散系数,m/s;p 为压力,Pa;为动力黏度,Pas;(p c p)e r t r 为煤体表面与参与空气对流导热的有效热容,J/(mK);H为垮落采空区的高度,m;Q.为遗煤产生的热量,J;c,为组分i的浓度,mol/m;R,为遗煤向顶底板扩散的组分i的浓度看,mol/m。氧气消耗量与氧气浓度和煤自
19、燃参数有关,煤氧反应热与氧气消耗量和反应热有关,分别计算如下:()Ro,=-Aexpco.Sv(7)Qa=Ro,4Q(8)式中:Ro,为氧气消耗量,mol;Q a 为煤氧仅应热,J;c o,为氧气浓度,mol/m;A 为指前因子,s-l;E为活化能,J/mol;R为气体常数,J/(mol-K);n为反应级数;Q为煤氧化学反应过程中生热量,J/mol;s v 为形状因子1.4工作面及巷道控制方程煤矿井下采场工作面和巷道内区域气体一般都处于瑞流状态,因此需要使用瑞流模型进行模拟。选用L-VEL湍流模型 2 4 ,该模型基于普朗特混合长度理论,根据局部流速和与最近壁面的距离来计算流黏度,稳定且计算效
20、率高,适合用于计算内部流动。流动、传热及传质方程如下:p(V.V)V=V.-pl+(+r)(VV+(VV)T)pV.V=0TPfV.VT-V.(aVT)=Q(9)otciV.(D,Vc)+V.Vc;=0ot式中:r为有效湍流黏度,kgm/s;I为压力在空向上的分布;p为流体密度,kg/m;C,为空气比热容,J/(kgK);2 为空气导热系数,W/(mK);D,为组分i的扩散系数,m/s。1.5围岩(煤)传热控制方程围岩(煤)固相传热以热传导为主,可用下式表示:8T(oC)+V.(a,VT)=0(10)t式中:(pC)为围岩(煤)热容,J/(mK);为导热系数,J/(mK)1.6多场耦合关系建立
21、的模型具有3个方面的耦合关系。首先,在采空区内部,采空区能量方程和气体组分方程相互耦合,并且均与采空区流动方程所求解的速度场相关联。其次,工作面和巷道气体传热方程和气体组分方程均与其内部速度场相关联。最后,温度、风速、压力、气体体积分数在采空区和工作面区域的界面上数值相等,即具有连续性。95Safety in Coal MinesOct.20232023年10 月Vol.54煤砺发全No.10第10 期第54 卷1.7工作面推进过程中计算网格的处理随着工作面向前推进,采空区面积不断增大,计算区域物理尺寸不断变化,因此计算网格也应随之调整。利用COMSOLMultiphysics软件的变形几何模
22、型,通过设定计算边界的移动速度,实现网格的自由变形,从而实现对采空区动态增长过程的模拟。1.8初始条件及边界条件通过模型耦合关系的分析,明确了采空区、面间煤柱、边界煤柱、底板岩层、顶板岩层、待回采实体煤、进风巷、回风巷等地点的初始条件和边界条件,如下:1)初始条件。对于采空区,各点温度T、氧气浓度Co,和一氧化碳浓度Cco设定为:t=0时,T=Tw,Co.=Co(C。为进风流O,浓度),Cco=0;对于面间煤柱、边界煤柱、底板岩层、顶板岩层和待回采实体煤:t=0时,T=Tw;对于进风巷、进风巷和工作面:t=0时,T=Tin,Co,=Co,Cc o=O。2)边界条件。对于进风巷及回风巷面间煤柱壁
23、面、始采线边界煤柱固壁和终采线边界煤柱固壁:n(a,VT)=O(n为与边界垂直的向量);对于进风巷人口处:u=Uin(uin为人口风速度),Co,=Co,Cc o=O,T=T i n(T n 为空气人口温度);对于回风巷出口处:p=O,n(a,VT)=O,n D;ci=O。进风巷及回风巷面间煤柱固壁、工作面固壁、始采线边界煤柱固壁和终采线边界煤柱固壁满足无滑移边界条件。2模型验证利用在开滦集团崔家寨煤矿E12604工作面获得的实测数据与建立的数学模型计算结果进行对比分析。崔家寨煤矿E12604工作面推进距离57 0m,倾向长度135m,煤层厚度3.0 3.2 m,属易自燃煤层,最短自然发火期9
24、 0 d,采煤方法为综采倾斜长壁、一次采全高,采高为2.8 3.0 m。现场数据监测点距离E12604工作面开切眼12 0 m,距回风巷底板2.5m。数值模拟参数为:入口风速vin=1.01m/s;空气人口温度Tm=295.85K;活化能E=21.6kJ/mol;进风流O,浓度C=9.375mol/m;指前因子A=0.85s;氧化反应热4 0 0.11kJ/mol;围岩初始温度2 9 5.8 5K;煤岩密度130 0 kg/m;煤岩比热容10 0 3J/(kgK);煤岩导热系数=0.2 W/(mK);初始碎胀系数Kp.mx=1.5;压实碎胀系数Kp.min=1.15;碎胀系数衰减率x=0.26
25、8ml;碎胀系数衰减率,=0.0368m;工作面采高3m;工作面推进速度7.5m/d;煤研石直径d,=0.04m;工作面长度2 38 m。崔家寨煤矿E12604工作面采空区气体温度和O,体积分数的模拟与实测结果如图2。50实测数据模拟数据40/30201004080120160200推进距离/m(a)监测点温度随推进距离的变化25实测数据模拟数据20%/1510504080120160200推进距离/m(b)监测点O2体积分数随推进距离的变化图2模拟与实测数据对比Fig.2Comparison of simulated and measured data从图2 可以看出:在工作面回采初期采空区
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 动态 推进 型采场煤 自燃 耦合 数值 研究