Research on change law of fire smoke characteristics during dynamic sealing process of mine tunnels
-
摘要:
为研究煤矿井下巷道火灾动态密闭过程中烟气特性变化规律，通过FDS软件对矿井巷道火灾进行全尺寸数值模拟；利用缩尺实验建立了巷道火灾CO体积分数衰减速率模型，研究2、3、4 MW火源功率下，巷道在50%至100%的密闭过程中横向、纵向温度分布及烟气蔓延规律；并建立巷道完全密闭后CO体积分数衰减模型，预测火灾熄灭时间。结果表明：密闭比例由50%至100%的变化过程中，顶棚中心温度先急剧下降后缓慢升高；密闭后火源上方温度呈上高下低的分布状态，且距离火源越近温度衰减速率越快；顶棚温度纵向分布出现整体升高现象，最高温度点向密闭侧偏移1~3 m；烟气在密闭侧沉积现象明显，分层现象消失；单侧完全密闭后，CO体积分数与时间呈线性关系，通过CO体积分数衰减速率模型预测2、3、4 MW火源功率下的火源分别于103.16、94.40、86.21 s后熄灭。
Abstract:In order to study the change law of smoke characteristics in the dynamic sealing process of coal mine roadway fire, FDS software was used to conduct full-scale numerical simulation of mine roadway fire, and the CO volume fraction decay rate model of roadway fire was established by scale experiment. The transverse and longitudinal temperature distribution and smoke spread law in the sealing process of 50% to 100% of the roadway were established, and the CO volume fraction attenuation model was established to predict the fire extinguishing time after the roadway was completely sealed. The results show that during the change of sealing ratio from 50% to 100%, the temperature at the center of the ceiling first sharply decreases and then slowly increases; after being sealed, the temperature above the fire source shows a distribution pattern of high temperature above and low temperature below, and the closer it is to the fire source, the faster the temperature decay rate; the vertical distribution of ceiling temperature shows an overall increase, the maximum temperature point shifts 1-3 m to the sealed side; the deposition phenomenon of smoke on the sealed side is obvious, and the stratification phenomenon disappears; the relationship between CO volume fraction and time is linear when one side is completely sealed. According to the CO volume fraction decay rate model, it is predicted that the ignition source of 2, 3 and 4 MW is extinguished after 103.16, 94.40 and 86.21 s, respectively.
-
Keywords:
- roadway fire /
- dynamic sealing process /
- temperature distribution /
- rapid sealing /
- flue gas
-
在煤矿井下，煤层瓦斯压力测定（比如胶囊-黏液封孔等）和钻孔瓦斯涌出初速度测定均需要向钻孔内送入封孔器。当封孔器在钻孔内移动时，封孔器端头与钻孔残留煤渣相互作用；一方面，钻孔残留煤渣会阻碍封孔器移动，另一方面，受封孔器端头力学作用，煤渣出现堆积现象；如果堆积严重，封孔器难以快速到达钻孔预定位置，参数测定可能失败（比如钻孔瓦斯涌出初速度测定，要求2 min内将封孔器送入预定位置并密封钻孔）。
为减轻钻孔残留煤渣对封孔器移动的影响，学者们从排渣和减阻2个方面进行了研究。在排渣方面：刘涛[1]提出了环缝式引射器输送钻孔煤渣技术；王永龙等[2-3]提出了双动力低螺旋钻杆排渣技术和刻槽钻杆排渣技术；杜长胜等[4]提出了脉动通风与钻杆分风联合作用排渣技术，李国强等[5]提出了旋风分离除渣技术；李栋等[6]提出了射流排水排渣技术；张辉等[7]提出了泵吸反循环钻进排渣技术；潘竞涛等[8]提出了下行定向钻孔氮气泡沫幂律多相流携渣解堵技术；汪义龙等[9]提出了煤层定向钻孔气举排渣技术；李成成[10]提出了“逐级升压切缝、低压旋转排渣”的排渣技术；张宏杰等[11]提出了“高螺旋机械排渣为主、液动排渣为辅”的复合排渣技术；童碧等[12]提出了基于整体式螺旋钻杆的复合排渣工艺技术；何将福等[13]提出了三通道反循环排渣技术；张宏钧等[14]提出了螺旋钻杆与压风复合排渣技术。在减阻方面：齐黎明等[15]提出了在封孔器前端安设锥面端头的减阻技术；周福宝等[16]提出了让高压气流从钻孔孔壁与抽采管之间流出，从而减轻煤渣阻碍抽采管移动的技术方案。
如果排渣技术过度使用，会扰动钻孔孔壁并产生新的煤渣甚至于塌孔；由此可见，排渣技术虽然可降低钻孔残留煤渣数量，但不能将残留煤渣清除；因此，仅采用排渣技术是不够的，需补充应用减阻技术。在封孔器前端安设锥面端头可以有效降低残留煤渣阻力[15]，是一种简单、实用的减阻技术；封孔器锥面端头（简称锥面端头）与钻孔残留煤渣的力学作用及其减阻机理是该技术推广应用中亟待解决的基础理论问题，并且锥面端头设计也需要进一步优化，在公开的文献资料中，未见相关研究成果。为此，从理论上分析锥面端头与钻孔残留煤渣的力学作用，探索锥面端头减阻机理，并进行工程应用。
1. 锥面端头前方钻孔残留煤渣堆积状况分类
锥面端头前方残留煤渣堆积状况如图1所示。
由图1可以看出，随着锥面端头移动，锥面端头前方堆积的残留煤渣（简称堆积煤渣）越来越厚。根据堆积煤渣是否超过锥面端头锥顶，锥面端头前方钻孔残留煤渣堆积状况可分为2类：堆积煤渣未超过锥顶和堆积煤渣超过锥顶。
2. 锥面端头移动最小推力理论分析
2.1 堆积煤渣未超过锥顶时的最小推力
假定锥面端头在钻孔内呈水平状态，以锥面端头锥底最低点O为原点，锥面端头移动方向（锥面端头轴线方向）为x轴正方向，垂直向上为z轴正方向，y轴正方向垂直于平面xOz，并指向右侧，建立直角坐标系，锥面端头受力分析如图2所示。
在图2中，锥面CAOB为锥面端头与煤渣的接触面（简称接触面），O1点为锥底圆心，2θ为锥面角度；C点为锥面上接触面与非接触面分界线的顶点，位于x轴正上方；h为堆积煤渣厚度。
锥面端头受力如下：封孔器推力F1整体作用于x轴正方向；重力G整体作用于z轴负方向；煤渣压力F2（挤压应力σ的总称）垂直于接触面并指向锥面端头内部；煤渣摩擦阻力F3与接触面相切并垂直于z轴（锥面端头以近水平形式移动，并不断破坏堆积煤渣的稳定性，与滑块在固定斜面的移动有着本质区别）。
为便于受力分析，需对煤渣压力与煤渣摩擦阻力沿x轴、y轴、z轴3个方向进行分解，并求出各分力计算方法。
y轴、z轴方向分力微元挤压应力分析如图3所示。由于图2中锥面CAO和锥面CBO沿平面CDO对称分布，因此，取锥面CBO进行分析即可。先作2个平行于yOz平面的截面（间距为dx），再从圆锥的中心轴线作2个截面（夹角为dβ），在锥面CBO上取微元（面积为ds），并进行受力分析。
![]()
图 3 y轴、z轴方向分力微元挤压应力分析Figure 3. Extrusion stress analysis of differential element for y-axis and z-axis component of force根据图3（a），微元面积可表示为：
$$ \mathrm{d}s=\frac{\mathrm{d}x}{\mathrm{c}\mathrm{o}\mathrm{s}\;\varphi }（H-x）\frac{r}{H}\mathrm{d}\beta $$ (1) 式中：s为锥面CBO的面积，m2；x为x轴坐标，m；φ为沿x轴方向的锥面角度，（°）；H为圆锥高度，m；r为锥面端头锥底半径，m；β为夹角，（°）。
根据图3（b），微元承受的挤压应力可分解为平行于x轴的分应力dF2x和沿着截面指向圆心的分应力dF2yz，dF2yz可表示为：
$$ {\mathrm{d}}{F}_{2{yz}}=\sigma （H-x）\frac{r}{H}{\mathrm{d}}x{\mathrm{d}}\beta $$ (2) 式中：F2yz为锥面CBO的煤渣压力沿着截面指向圆心分力，N；σ为接触面挤压应力，Pa。
对F2yz进行分解，力分解示意图如图4所示。
根据图4，F2yz平行于y轴的分力F2y可表示为：
$$ F_{2{{{y}}}}=\int_0^{\frac{h}{\mathrm{t}\mathrm{a}\mathrm{n}\; \theta}}\; \sigma（H-x）\frac{r}{H}（1-\mathrm{cos}\; \alpha）\mathrm{d}x $$ (3) 式中：F2y为F2yz平行于y轴的分力（锥面CBO的煤渣压力沿y轴方向分力），N；h为堆积煤渣厚度，m；θ为锥面端头锥面角度的1/2值，（°）；α为弧线对应夹角，（°）。
根据三角函数关系，有：
$$ （H-x）\frac{r}{H}（1-\mathrm{c}\mathrm{o}\mathrm{s}\;\alpha ）=h-x\mathrm{t}\mathrm{a}\mathrm{n}\;\theta $$ (4) 将式（4）代入式（3）并积分，有：
$$ {F}_{2{y}}=\sigma \frac{{h}^{2}}{2\mathrm{t}\mathrm{a}\mathrm{n}\;\theta }=\sigma {S}_{\vartriangle CDO} $$ (5) 式中：${S}_{\vartriangle CDO} $为$\vartriangle $CDO的面积，m2。
$\vartriangle $CDO为锥面CBO沿y轴的投影，因此，因此，锥面CBO的煤渣压力沿y轴方向分力等于挤压应力σ与投影面积SCDO之积。
同理，F2yz平行于z轴的分力F2z可表示为：
$$ {F}_{2{z}}={\int }_{0}^{\frac{h}{\mathrm{t}\mathrm{a}\mathrm{n}\;\theta }}\sigma （H-x）\frac{r}{H}\mathrm{s}\mathrm{i}\mathrm{n}\;\mathrm{\alpha }\mathrm{d}x $$ (6) 式中：F2z为F2yz平行于z轴的分力（锥面CBO的煤渣压力沿z轴方向分力），N。
根据三角函数关系，有：
$$ （H-x）\frac{r}{H}\mathrm{s}\mathrm{i}\mathrm{n}\;\mathrm{\alpha }=\sqrt{2Rh-{h}^{2}}（1-\frac{\mathrm{t}\mathrm{a}\mathrm{n}\;\theta }{h}x） $$ (7) 将式（7）代入式（6）并积分，有：
$$ {F}_{2{z}}=\sigma \frac{h}{2\mathrm{t}\mathrm{a}\mathrm{n}\;\theta }\;\sqrt{2Rh-{h}^{2}}=\sigma {S}_{\vartriangle CBD} $$ (8) 式中：${S}_{\vartriangle CBD} $为三角形CBD的面积，m2。
$\vartriangle $CBD为锥面CBO沿z轴的投影，因此，锥面CBO的煤渣压力沿z轴方向分力等于挤压应力σ与投影面积SCBD之积。
x轴方向分力微元挤压应力分析图如图5所示。图3中平行于x轴的微元分应力dF2x是φ的函数，由于φ非固定值，不便求解；为此，先作2个平行于xOy平面的截面（间距为dz），再作2个平行于xOz平面的截面（间距为dy），在锥面CBO上取微元（面积为ds），并进行受力分析。
![]()
图 5 x轴方向分力微元挤压应力分析图Figure 5. Extrusion stress analysis of differential element for x-axis component of force
同理，图5中微元面积可表示为：
$$ \mathrm{d}s=\frac{\mathrm{d}y}{\mathrm{s}\mathrm{i}\mathrm{n}\;\varphi }\mathrm{d}z $$ (9) 式中：y为y轴坐标，m；z为z轴坐标，m。
图5中微元承受的挤压应力平行于x轴的分应力dF2x可表示为：
$$ \mathrm{d}{F}_{2{x}}=\sigma \mathrm{d}y\mathrm{d}z $$ (10) 式中：F2x为锥面CBO的煤渣压力沿x轴方向分力，N。
对式（10）进行积分，有：
$$ {F}_{2{x}}=\frac{\sigma }{2}（{r}^{2}\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{c}\mathrm{o}\mathrm{s}\frac{r-h}{r}-\sqrt{2rh-{h}^{2}}（r-h））=\sigma {S }_{BDO} $$ (11) 式中：SBDO为割圆BDO的面积，m2。
割圆BDO为锥面CBO沿x轴的投影，因此，锥面CBO的煤渣压力沿x轴方向分力等于挤压应力σ与投影面积SBDO之积。
同理，锥面CBO的煤渣摩擦阻力在x轴、y轴2个方向的分力F3x、F3y等于挤压应力σ、摩擦阻力系数f和投影面积（投影平面与分力方向平行且包含z轴，投影面积分别为SCDO、SBDO）之积，可表示为：
$$ {F}_{3{x}}=\sigma f{S}_{CDO}=\mathrm{\sigma }f\frac{{h}^{2}}{2\mathrm{tan}\;\theta } $$ (12) $$ {F}_{3{y}}=\frac{\sigma }{2}f（{r}^{2}\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{c}\mathrm{o}\mathrm{s}\frac{r-h}{r}-\sqrt{2rh-{h}^{2}}（r-h）） $$ (13) 式中：F3x为锥面CBO的煤渣摩擦阻力沿x轴方向分力，N；F3y为锥面CBO的煤渣摩擦阻力沿y轴方向分力，N；f为锥面端头与煤渣的最大静摩擦系数。
当封孔器推力F1达到一定数值（最小推力F）时，锥面端头才能突破煤渣阻力向前移动，x轴、y轴、z轴3个方向的受力均处于平衡状态。y轴方向受力呈对称分布，自然平衡；x轴、z轴方向受力平衡方程可表示为：
$$ F=2{F}_{2{x}}+2{F}_{3{x}} $$ (14) $$ G=2{F}_{2\textit{z}} $$ (15) 式中：F为锥面端头移动最小推力，N；G为锥面端头重力，N。
将式（8）、式（11）、式（12）代入式（14）和式（15），并解方程，有：
$$\begin{split} &\begin{gathered} F = G\left(\mathrm{t}\mathrm{a}\mathrm{n}\; \theta\left(\frac{r^2\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{c}\mathrm{o}\mathrm{s}\; \frac{r-h}{r}}{h\sqrt{2rh-h^2}} - \frac{r-h}{h}\right) + f\sqrt{\frac{h}{2r-h}}\right) \\ 0 < h \leqslant r \end{gathered}\\[-16pt]& \end{split} $$ (16) 根据式（16）可知，锥面端头移动最小推力F与端头重力G、端头锥面角度的一半值θ、端头锥面与煤渣的最大静摩擦系数f属于正相关，与堆积煤渣厚度h和端头锥底半径r的相关性，无法直接判定。
假设端头重力G为50 N，端头锥面角度的一半值θ为30°，最大静摩擦系数f为0.1，通过数值计算并绘制得到的锥面端头移动最小推力F随堆积煤渣厚度h和端头锥底半径r的变化曲线如图6所示。
![]()
图 6 锥面端头移动最小推力与堆积煤渣厚度和端头锥底半径的关系曲线Figure 6. Relation curves between minimum thrust of cone end movement and height of accumulated coal dregs and cone base radius of end根据图6可知：锥面端头移动最小推力F与堆积煤渣厚度h和端头锥底半径r也属于正相关。
2.2 堆积煤渣超过锥顶时的受力分析
当堆积煤渣超过锥顶时，根据是否超过锥面端头锥底最高点，可将残留煤渣堆积情况进一步划分，本研究仅分析未超过锥面端头锥底最高点的情形。
锥面端头与煤渣发生力学作用的面积相当于整个锥面与非接触面的面积之差，并忽略锥面端头锥面上方煤渣重力。在式（16）中，令h=r，求得的F值乘以2，即为整个锥面与煤渣接触条件下的锥面端头移动最小推力；再令h=2r−h，求得的F值，即为非接触面与煤渣接触条件下的锥面端头移动最小推力；上述两者之差即为堆积煤渣超过锥顶时的锥面端头移动最小推力F，如式（17）：
$$\begin{split} &\begin{gathered} F=G\mathrm{t}\mathrm{a}\mathrm{n}\;\theta \left(\pi +\frac{f}{\mathrm{t}\mathrm{a}\mathrm{n}\;\theta}\left(1-\sqrt{\frac{2r}{h}-1}\right)-\right.\\ \left.\left. \frac{{r}^{2}\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{c}\mathrm{o}\mathrm{s}\;\frac{h-r}{r}}{\left(2r-h\right)\sqrt{2rh-{h}^{2}}}+\frac{h-r}{2r-h}\right)\right. \qquad r < h \leqslant 2r \end{gathered}\\[-18pt]& \end{split} $$ (17) 3. 锥面端头移动最小推力试验研究
锥面端头移动最小推力F的影响因素较多，因篇幅受限，仅试验研究锥面端头移动最小推力随端头锥面角度的变化规律，并验证理论分析结果的可靠性。
试验方案如下：采用弹簧秤牵拉带有锥面端头的封孔器并让其在含煤渣的管道内匀速移动，弹簧秤显示读数即为锥面端头移动最小推力；锥面角度分别是60°、45°、30°、20°、18°、15°，共6组试验。试验装置设计图如图7所示，锥面端头移动最小推力与锥面角度的关系曲线如图8所示。
![]()
图 8 锥面端头移动最小推力与锥面角度的关系曲线Figure 8. Relation curves between minimum thrust of cone end movement and conical surface angle由图8可知：锥面角度增加，锥面端头移动最小推力升高，并且上升速率加快；移动距离越远，上升速率加快趋势越明显。根据式（16），锥面端头移动最小推力F与tan θ成正比，tan θ为正切函数，随θ增大而升高并且上升速度逐渐加快。由此可见，试验研究结果与理论分析结果基本一致。
4. 锥面端头减阻机理
与平面端头相比，锥面端头可以通过接触面对堆积煤渣施加平行于y轴的侧向力；因此，可通过分析侧向力来探索锥面端头减阻机理。因篇幅受限，以堆积煤渣未超过锥顶的情形为例，进行分析。
堆积煤渣承受的侧向力Fy是锥面端头y轴方向受力的反作用力，方向由接触面指向堆积煤渣，并平行于y轴，有：
$$ {F}_{{y}}={F}_{2{y}}-{F}_{3{y}} $$ (18) 将式（5）、式（13）、式（15）代入式（18），有式（19）：
$$\begin{split} &\begin{gathered} {F}_{{y}}=G\frac{h}{2\sqrt{2rh-{h}^{2}}}-\frac{G\mathrm{t}\mathrm{a}\mathrm{n}\;\theta }{h\sqrt{2rh-{h}^{2}}}f（{r}^{2}\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{c}\mathrm{o}\mathrm{s}\;\frac{r-h}{r}-\\ \sqrt{2rh-{h}^{2}}（r-h）） \end{gathered}\\[-19pt]& \end{split} $$ (19) 式中：Fy为堆积煤渣承受的侧向力，N。
根据式（19）可知，θ越小，Fy越大，即锥面端头越尖，Fy越大。
与锥面端头相接触的煤渣，在侧向力Fy的推动下，向钻孔孔壁侧移动（或被压缩）；导致锥面端头通过后，锥面端头后方残留煤渣增加，锥面端头前方堆积煤渣减少，堆积煤渣厚度h降低；锥面端头越尖，侧向力Fy越大，堆积煤渣厚度h越低。
锥面端头移动的突破点通常位于接触面顶端（堆积煤渣未超过锥顶时，在C点；堆积煤渣超过锥顶，在锥面端头锥顶），如果接触面顶端承受的应力超过该点的堆积煤渣剪切应力，则堆积煤渣遭到破坏、锥面端头可向前移动。
由于堆积煤渣相邻颗粒之间具有一定的咬合作用，图9（b）的某片煤渣薄板受侧向力挤压后，受压效应向前传递（距离越远，效应越弱），导致该煤渣薄板至钻孔孔壁一侧应力升高，而至钻孔中心一侧应力下降（相当于1根受压弹簧，从某个位置施加新力后，受压一侧弹簧绷紧程度加强，而另一侧弹簧绷紧程度减弱）；虽然煤渣属于散体，单片煤渣薄板微观上的受压效应传递能力较低，但是，当众多煤渣薄板的微观受压效应传递累积起来时，宏观上可以在一定程度上影响应力分布。
在全体煤渣薄板受压效应传递综合作用下，从C点到钻孔孔壁（E1、E2）的堆积煤渣应力变化如下：靠近C点一侧应力下降，靠近钻孔孔壁（E1、E2）一侧应力上升，C点（位于钻孔中心，并且距侧向力最近）应力下降大。根据散体动力学理论[17]，摩尔库伦破坏准则适用于散体，堆积煤渣剪切应力与正应力成正比。因此，受侧向力挤压效应传递的影响，堆积煤渣C点及其邻近区域的剪切应力降低，C点下降值最大。
综上所述，锥面端头减阻机理如下：锥面端头越尖（θ越小），侧向力Fy越大，堆积煤渣厚度h越低、接触面顶端及其临近区域堆积煤渣剪切应力越小，堆积煤渣对锥面端头移动的阻碍能力越弱。
注：C点及其邻近区域堆积煤渣受力破坏后，堆积煤渣整体失去稳定性，靠近钻孔孔壁（E1、E2）一侧堆积煤渣应力自然得到释放。
5. 封孔器端头优化改进与现场试验
根据锥面端头移动最小推力与影响因素的关系，可从以下几个方面对锥面端头开展优化设计：①采用低密度材料制作锥面端头，降低端头重力；②适当减小端头锥面角度；③采用喷漆等措施提高端头锥面的光滑性，降低端头锥面与煤渣的最大静摩擦系数；④采用更细的封孔器，从而减小端头锥底半径。
在文献[15]中，将1个空心外罩（锥面角度约20°）安设于钻孔瓦斯涌出初速度测定装置的封孔器端头，端头由台阶式转变为锥面式；在此基础上，根据锥面端头优化设计思路，对封孔器端头作进一步改进。
根据图8，当锥面角度较低时，降低锥面角度改变锥面端头移动最小推力的作用不大；由于端头锥底半径与封孔器主体部分一致，而主体部分需要考虑膨胀性能、耐压性能和膨胀液输送、内置管道气体流动阻力等因素，截面尺寸（目前为38 mm）很难进一步缩小。因此，从如下2个方面对封孔器端头进行改进：①将空心外罩的制作材料由不锈钢改合金，减轻端头质量；②对封孔器锥面端头打蜡，降低端头锥面与煤渣的最大静摩擦系数。
在开滦矿区钱家营矿1355工作面进行了钻孔瓦斯涌出初速度测定现场试验，现场试验数据见表1。
表 1 现场试验数据Table 1. Experimental data in field test封孔器型号 最高送入孔深/m 孔深/m 时间/min 原始型 16 8 1.2 10 1.3 12 1.5 14 1.7 16 2.0 改进型 >20 8 0.9 10 1.0 12 1.2 14 1.3 16 1.5 18 1.7 20 1.9 由表1可知：改进型封孔器送入阻力小、推进速度快、送入钻孔深度大。
6. 结 语
1）锥面端头与煤渣接触面所承受的煤渣压力沿x轴、y轴、z轴的分力均等于挤压应力与投影面积之积；锥面端头与煤渣接触面所承受的摩擦阻力沿x轴、y轴的分力也均等于挤压应力、摩擦阻力系数和投影面积之积。
2）从理论上推导了锥面端头移动最小推力计算公式，它与端头质量、端头锥面角度、端头锥面与煤渣的最大静摩擦系数、堆积煤渣厚度和端头锥底半径均具有正相关性，为锥面端头优化设计指明了方向。
3）锥面端头移动最小推力与锥面角度呈正切函数关系，随锥面角度增大而升高并且上升速度逐渐加快。
4）通过侧向力分析，揭示了锥面端头减阻机理；锥面端头越尖，侧向力越大，堆积煤渣厚度越低、接触面顶端及其临近区域堆积煤渣剪切应力越小，堆积煤渣对锥面端头移动的阻碍能力越弱。
5）根据锥面端头优化设计思路改进后的封孔器端头，可以降低送入阻力、提高送入速度与孔深。
-
![]()
图 3 不同火源功率下顶棚中心点温度变化图
Figure 3. Temperature variation diagram of ceiling center point under different fire power conditions
![]()
图 7 巷道纵向温度分布对比曲线
Figure 7. Comparison curves of longitudinal temperature distribution in tunnels
表 2 CO衰减速率模拟合结果
Table 2 Simulation results of CO attenuation rate
序号 火源功率/MW 拟合公式 斜率 相关系数R2 1 2 C1=−0.18t1+18.57 0.18 0.915 2 3 C2=−0.25t2+23.60 0.25 0.879 3 4 C3=−0.33t3+28.45 0.33 0.975 注：C1、C2、C3为火源功率2、3、4 MW时的CO体积分数；t1、t2、t3为火源功率2、3、4 MW时的反应时间。 
下载: 导出CSV
-
[1] 薛彦平. 巷道分岔角度对火灾烟气蔓延的影响数值模拟[J]. 煤矿安全,2020,51(6):179−183. XUE Yanping. Numerical simulation of effect of tunnel diverging angles on fire smoke flow[J]. Safety in Coal Mines, 2020, 51(6): 179−183.
[2] 董炳燕,吴晋湘,黄有波. 封堵比例对隧道火灾温度影响的数值模拟研究[J]. 消防科学与技术,2018,37(8):1051−1054. doi: 10.3969/j.issn.1009-0029.2018.08.012 DONG Bingyan, WU Jinxiang, HUANG Youbo. Numerical investigation on the effect of tunnel portal sealing ratio on gas temperature during tunnel fire[J]. Fire Science and Technology, 2018, 37(8): 1051−1054. doi: 10.3969/j.issn.1009-0029.2018.08.012
[3] 马砺,李超华,张鹏宇,等. 封堵比例对隧道火灾顶棚温度分布的影响[J]. 中国安全科学学报,2020,30(2):79−85. MA Li, LI Chaohua, ZHANG Pengyu, et al. Influence of sealing ratio on ceiling temperature distribution of tunnel fire[J]. China Safety Science Journal, 2020, 30(2): 79−85.
[4] CHEN C K, XIAO H, WANG N N, et al. Experimental investigation of pool fire behavior to different tunnel-end ventilation opening areas by sealing[J]. Tunnelling and Underground Space Technology, 2017, 63: 106−117. doi: 10.1016/j.tust.2017.01.001
[5] YAO Y Z, CHENG X D, SHI L, et al. Experimental study on the effects of initial sealing time on fire behaviors in channel fires[J]. International Journal of Thermal Sciences, 2018, 125: 273−282. doi: 10.1016/j.ijthermalsci.2017.11.031
[6] 彭敏,吴振坤,蒋顺,等. 纵向通风对隧道内地铁列车车厢烟气温度特性的影响[J]. 消防科学与技术,2023,42(12):1656−1662. doi: 10.3969/j.issn.1009-0029.2023.12.011 PENG Min, WU Zhenkun, JIANG Shun, et al. Study on the effect of longitudinal ventilation on the smoke temperature characteristics of subway train fires inside tunnel[J]. Fire Science and Technology, 2023, 42(12): 1656−1662. doi: 10.3969/j.issn.1009-0029.2023.12.011
[7] 李晴,康建宏,周福宝,等. 全尺寸巷/隧道火灾风烟流温度预测模型与验证[J]. 中国安全生产科学技术,2022,18(8):5−12. LI Qing, KANG Jianhong, ZHOU Fubao, et al. Prediction model and verification of smoke flow temperature in full-scale roadway/tunnel fires[J]. Journal of Safety Science and Technology, 2022, 18(8): 5−12.
[8] HUANG Y B, LI Y F, DONG B Y, et al. Numerical investigation on the maximum ceiling temperature and longitudinal decay in a sealing tunnel fire[J]. Tunnelling and Underground Space Technology, 2018, 72: 120−130. doi: 10.1016/j.tust.2017.11.021
[9] JIANG X P, ZHANG H X, JING A. Effect of blockage ratio on critical velocity in tunnel model fire tests[J]. Tunnelling and Underground Space Technology, 2018, 82: 584−591. doi: 10.1016/j.tust.2018.09.001
[10] TANG Z Q, GAO K, SHAN Y F, et al. Study of the fire behavior of multilayer cables in a mine tunnel[J]. Energies, 2022, 15(6): 2059. doi: 10.3390/en15062059
[11] LI J X, LI Y F, LI J M, et al. Numerical investigation on the smoke behaviour and longitudinal temperature decay in tilted tunnel fire with portal sealing[J]. Indoor and Built Environment, 2023, 32(1): 133−148. doi: 10.1177/1420326X211034894
[12] 马砺,张鹏宇,郭睿智,等. 巷道火灾密闭过程烟气温度预测的GA-SVM模型[J]. 中国矿业大学学报,2021,50(4):641−648. MA Li, ZHANG Pengyu, GUO Ruizhi, et al. GA-SVM model for prediction flue gas temperature of roadway fire under sealing process[J]. Journal of China University of Mining & Technology, 2021, 50(4): 641−648.
[13] SHI J K, ZUO C, XIONG Y Y, et al. Experimental study of different sealing ratios on the self-extinction of tunnel fires[J]. Tunnelling and Underground Space Technology, 2021, 112: 103894. doi: 10.1016/j.tust.2021.103894
[14] 李智胜,蒋浩锴,高云骥,等. 分岔隧道火灾烟气流动特性研究[J]. 中国安全生产科学技术,2021,17(5):118−122. LI Zhisheng, JIANG Haokai, GAO Yunji, et al. Study on flow characteristics of fire smoke in bifurcated tunnel[J]. Journal of Safety Science and Technology, 2021, 17(5): 118−122.
[15] 沈云鸽,王德明. 基于FDS的矿井巷道火灾烟气致灾的数值模拟[J]. 煤矿安全,2020,51(2):183−187. SHEN Yunge, WANG Deming. Numerical simulation of smoke disaster caused by mine roadway fire based on FDS[J]. Safety in Coal Mines, 2020, 51(2): 183−187.
[16] 李超华. 隧道火灾封闭过程中烟气温度与流动特性试验研究[D]. 西安:西安科技大学,2020. [17] HU L H, TANG F, YANG D, et al. Longitudinal distributions of CO concentration and difference with temperature field in a tunnel fire smoke flow[J]. International Journal of Heat and Mass Transfer, 2010, 53(13/14): 2844−2855.
[18] LI T, CHEN L F, ZHANG Y C, et al. Study on the smoke movement characteristics in large scale interchange tunnel fire[J]. International Journal of Ventilation, 2020, 19(3): 224−232. doi: 10.1080/14733315.2019.1693174
计量
- 文章访问数: 36
- HTML全文浏览量: 0
- PDF下载量: 2
