赵艳秋, 王昌建, 张 旭, 吴 越, 於 星
(合肥工业大学 土木与水利工程学院,安徽 合肥 230009)
地形条件是影响森林地表火蔓延的重要因素,文献[1]的研究表明,世界各地由森林火灾引起的致命事故大都发生在山谷或者山脊中,而山谷和山脊地形在实际中很常见,因此研究山地火蔓延对预防和扑救森林火灾具有重要意义;文献[2]研究了地形条件对火蔓延特性的影响,结果表明地形的坡度是影响火灾行为非常重要的因素。目前,许多学者已经开展了坡度条件下火蔓延的实验研究和模型分析。文献[3-4]改变地形坡度进行了地表火蔓延实验,测量了火焰长度、火焰倾斜角、火焰辐射和温度分布等特性,发现火焰在斜坡上蔓延时,辐射传热在火蔓延预热阶段占主导地位;文献[5]认为外界风和坡度可以改变辐射和对流传热效率,从而改变火蔓延速率(rate of fire spread,ROS),将实验得到的有坡度与无坡度时的ROS之比定义为坡度附加系数,并将坡度附加系数引入地表火灾蔓延数学模型中,使该模型可适用于有坡度的地形;文献[6-7]研究了风和坡度对蔓延行为的作用,认为坡度对ROS的影响与环境风的影响相似;文献[8-10]通过线形火点燃的上坡火蔓延实验,分析了坡度对ROS、火焰滞留时间、火焰温度和几何形状等火焰特性的影响,研究了火焰前锋对前方未燃烧区域的辐射和对流作用;文献[11-12]比较分析了点火源引燃和线形火点燃的火焰前锋的变化规律,结果发现辐射传热和对流传热在上坡火蔓延阶段共同作用,导致ROS不稳定,较大坡度时,对流加热效应的增强会引发爆发火;文献[13-14]在实验室和野外进行了实验,测量了上坡火蔓延过程中火焰前锋不同位置的热流,并观察了坡度对火焰动力学和火焰前方传热的影响。实验发现随着坡度增加,火焰动力和传热均发生显著变化。测得的总热流曲线与辐射热流曲线差距很小,即在该实验中,辐射传热对未燃烧的燃料的预热起主要作用。文献[15]在实验中将粒子图像测速(particle image velocity,PIV)技术测量系统与视频相结合,研究上坡火蔓延中火焰的动力学特性和燃烧机制;文献[16-18]通过在燃料床两侧布置相同高度的挡板抑制火焰的侧向卷吸,研究了沟槽地形的上坡火蔓延,实验结果表明,在较高的坡度下火焰会向燃料床表面附着,导致火蔓延加速,火焰前方产生热气流,增强对流加热作用从而引发爆发火;文献[19]研究了坡度对沟槽地形火灾行为的影响,结果表明坡度越大,爆发火喷射得越剧烈,达到最高温度的速度越快;文献[20-22]开展了一系列的实验,研究不同坡度下的下坡火灾行为,结果表明下坡ROS较低,随着坡度的增加,ROS在20°达到最小值,随后增加,其变化趋势呈抛物线形。前人关于山地火灾的实验研究主要集中在单一的地形,对复杂地形的火蔓延的实验和理论研究还不够完善。为了研究山地火灾的蔓延规律,本文利用一个山脊和山谷相连的模型进行实验,通过分析实验数据研究坡度对火灾行为和ROS的影响,并分析火焰前方的传热机理。
本文的所有实验均在一个长为30.2 m、宽为8.4 m、高为6.5 m的大厅内进行。实验时,大厅内的门窗处于关闭状态,排除气流对火蔓延过程的影响,模拟无风条件下的山地火蔓延。山地火蔓延实验在1个山脊和1个山谷相连的模型上进行,实验台是用3面1.0 m×2.0 m的相同尺寸金属板通过折页连接而成,山地的几何结构可以用中心坡度角α和侧面坡度角δ来描述,其中:中心坡度角α为山脊和山谷线与水平地面之间的夹角;侧面坡度角δ为山脊和山谷倾斜底边与水平地面之间的夹角。
在实验台上铺设防火板,并在防火板上放置细钢丝网,以防止燃料滑落,使燃料铺设更均匀。实验中使用原产地四川的干松针作为燃料,燃料铺设厚度约5 cm。实验前使用温湿度计以及含水率测定仪测得环境温度、空气相对湿度和燃料的含水率,实验参数见表1所列。每次实验时,同时开展实验台尺寸为1.0 m×1.0 m的平坡下的火蔓延实验,得到的火蔓延速率在0.27 ~0.30 cm/s范围内,重复性较好。因此,认为可以排除环境条件对实验的影响。
表1 实验参数
用浸有酒精的棉球点火,点火位置位于山脊中心线上距底部20 cm。改变中心坡度角α为0°、10°、20°、30°、40°,侧面坡度角δ为0°、10°、20°、30°、40°,进行10种工况的重复实验。
实验仪器布置如图1所示,采用K型热电偶(直径为1 mm,探头长度为300 mm,测温范围为-200~1 300 ℃,测温精度为±1.5 ℃)测量不同位置处的火焰前沿温度。热电偶从燃料床下穿至燃料床表面上方3 cm处,沿山脊线和山谷线以20 cm等距离分布。使用Captec公司生产的2个HS34-CB总热流计,量程为±200 kW/m2,灵敏度为2.5 μV/(W/m2))和2个TS-34C辐射热流计量程为±200 kW/m2,灵敏度为2 μV/(W/m2))测量火焰前方燃料表面接收到的总热量和辐射热量。在山脊和山谷延长线上,距离燃料床末端40 cm处分别安装1个总热流计和1个辐射热流计,测量窗口朝上。 使用NI公司生产的NI9213热电偶输入模块(测量灵敏度为0.02 ℃,电压测量范围为±78 mV)以及C-7018数据采集模块(电压范围为±15 mV,采样频率为10 次/s,精确度为±0.1%)采集热电偶和热流计的信号。在燃料床侧面和前面架设DV数字摄像机来记录火焰的燃烧过程,视频记录格式为25 帧/s。
图1 实验仪器布置图
坡度为0°时的火蔓延过程如图2a所示,平坡时各方向的ROS相同,因此会形成以点火位置为圆心,半径不断扩大的半圆形火线。坡度条件下,山脊火蔓延过程如图2b所示,上坡ROS较快,火线呈现V形,火焰附着在燃料床上,自点火位置向其他方向的火蔓延速率较慢,火线呈现半圆形,因此在山脊地形上的火线会形成水滴形,2条火线之间的夹角会先逐渐增大直至稳定不变。山谷地形火蔓延过程如图2c、图2d所示,山脊侧坡的火线继续蔓延至山谷线后,会迅速沿着山谷侧坡的上坡方向以及山谷线蔓延,在山谷的侧坡上形成V形火线。随着α的增加,侧坡上的火头会向上坡推移,快速蔓延至燃料床末端。
图2 山脊和山谷火线轮廓随时间的变化
本文利用山脊线和山谷线上的热电偶阵来确定火焰前锋到达某点和离开该点的时间,将热电偶的温度初次达到燃料着火温度(350 ℃)对应的时刻作为火焰前锋到达这个热电偶的时刻,然后根据热电偶的位置坐标以及所对应的时刻,对这些数据点进行一次线性拟合,拟合曲线的斜率为ROS。 不同中心坡度下山脊、山谷火焰前锋位置与点火时刻的线性拟合曲线以及拟合曲线的斜率如图3所示。图3中,R为ROS值,下同。
图3 山脊、山谷火焰前锋位置与点火时刻的拟合曲线
因为线性回归的相关系数均大于0.9,所以在这段距离内火焰前锋在山脊和山谷线上的蔓延速率是稳定的。对比图3a、图3b可以看出,山脊火蔓延得比较缓慢、稳定,山谷ROS远高于山脊的。
ROS随α和δ的变化如图4所示。从图4可以看出,山脊和山谷ROS均会随着α的增加而增加,随着δ的增加而减小。山脊ROS随α和δ的变化较小,山谷ROS的变化较大。文献[14]的研究表明,火焰在上坡蔓延时,随着坡度增加,倾斜角减小,火焰向上坡方向的燃料床倾斜,ROS增加。本文中,火焰在山脊和山谷上蔓延时,当α不变,δ不断增加,甚至大于α时,山脊和山谷火焰由向山脊线和山谷线倾斜转而向侧坡倾斜,火蔓延方向随之发生改变,导致山脊和山谷线上的ROS随着δ的增加而减小。
图4 ROS随中心坡度和侧面坡度的变化
将V形火线之间的夹角定义为火线夹角。火焰蔓延至山谷线后向山谷侧坡上坡方向蔓延,侧坡上火线与山谷线之间形成的夹角也是V形,为了方便计算,将这个夹角作为山谷火焰的火线夹角的1/2。通过对火线演化的研究发现,α和δ是影响火线夹角的重要因素。从实验视频的图像中提取火蔓延达到稳定时的火线夹角如图5所示,山脊和山谷火线夹角随坡度的变化如图6所示。从图5、图6可以看出,随着α的增加,山脊与山谷火线夹角明显减小;随着δ的增加,山脊火线夹角缓慢减小,而山谷火焰的火线夹角快速增加,这是因为侧面坡度低,中心坡度高,山脊侧坡的火线蔓延过山谷线后快速向上坡蔓延,而其他方向上的ROS较低,所以山谷火线夹角会随着δ的增加而增加。
图5 不同坡度下的火线夹角
图6 火线夹角随坡度的变化
将火头位置与火焰底部中心之间的长度定义为火焰长度。文献[23]关于火焰长度的研究表明,火焰长度与火焰的强度呈正比关系,文献[24]给出火焰强度的经验公式,即
IF=QfmfR
(1)
其中:IF为火焰强度;mf为单位面积消耗的燃料质量;R为ROS;Qf为燃料的热值,对于大部分燃料,燃料热值的差别不大。
由式(1)可知,火焰长度会随着ROS的增加而增加。本文提取了山脊和山谷火蔓延达到稳定阶段的火焰长度,山脊和山谷火线长度随坡度的变化如图7所示。
图7 山脊和山谷火线长度随坡度的变化
从图7可以看出,在坡度条件下,山谷火焰较长,其长度明显高于山脊火焰长度。随着α和δ的增加,山脊火焰长度缓慢增加,山谷火线长度迅速增加。
平坡时,山脊线和山谷线上的火焰基本保持竖直,ROS较小,火焰长度也小,随着α的增加,火焰向燃料床倾斜,ROS加快,火焰长度明显增加,与式(1)得到的结论一致。δ增加时,火焰向燃料床侧坡的上坡方向倾斜,虽然ROS缓慢减小,但火焰长度仍然增加。
在地表火蔓延的预热阶段,火焰前锋对于前方未燃烧的燃料的预热作用是影响火焰特性的重要因素,因此研究火蔓延特性主要分析预热阶段的传热机理。文献[25]将燃料层看作热薄型,认为在燃料床的厚度方向上没有温差,燃料床内部的热传导作用可以忽略。本文主要考虑来自火焰表面的辐射热、燃料表面与气体之间对流传热。
利用安装在燃料床末端40 cm处的总热流计和辐射热流计来模拟火焰前方未燃烧区域的燃料表面接收到的总热流和辐射热流。实验测得的热流曲线存在剧烈的波动,为了减小误差,数据分析时对热流曲线进行平滑处理。不同坡度下山脊和山谷的热流随时间的变化如图8、图9所示。
图8 不同坡度下山脊火热流随时间的变化
图9 不同坡度下山谷火热流曲线随时间的变化
对比图8、图9中可知,辐射热流峰值总是高于总热流峰值,这表明山脊与山谷火蔓延预热阶段的对流传热表现为对流冷却作用。由图8可知,随着α和δ增加,对流冷却不断增加,然而辐射传热量远高于对流传热量,因此在山脊火蔓延预热阶段辐射传热占主导作用。由图9可知:α和δ低于30°时,对流传热量很小,辐射传热在山谷火蔓延预热阶段占主导;α和δ高于30°时,对流冷却作用明显增强;在α=40°,δ=30°时,对流传热量与辐射热量相当;因此当α和δ高于30°时,在山谷火蔓延预热阶段辐射传热和对流传热共同作用。
山脊火和山谷火的辐射热流峰值随中心坡度和侧面坡度的变化如图10所示。由图10可知,随着α和δ的增高,火焰辐射热流峰值增大,当α和δ大于20°,山脊火辐射热流峰值缓慢增长,山谷火的辐射热流峰值快速增长,且山谷火焰辐射热量远高于山脊火焰。
图10 辐射热量峰值随中心坡度和侧面坡度的变化
本文改变中心坡度和侧面坡度,开展了较高坡度下的山脊和山谷火蔓延特性实验研究,得出如下结论:
1) 山脊火在坡度条件下会形成水滴形火线,上坡V形火线夹角会逐渐增大直至保持稳定,山脊侧坡火线蔓延至山谷线后会快速沿着山谷侧坡上坡和山谷线蔓延,在山谷侧坡上形成V形火线。
2) 中心坡度角α和侧面坡度角δ都对火蔓延行为特性产生影响。随着α的增加,山脊ROS缓慢增加,山谷ROS快速增加,山脊和山谷的火线夹角快速减小、火线长度明显增加。随着δ的增加,山脊和山谷ROS减小,火焰长度增加;山脊火线夹角减小,而山谷火焰前锋向侧坡快速蔓延,导致火线夹角增加。
3) 辐射传热在山脊火蔓延的预热阶段占主导地位,当α和δ低于30°时,辐射传热在山谷火蔓延起主要作用;高于30°时,对流传热量明显增加,辐射和对流传热共同作用。当α和δ高于20°时,山谷火焰对燃料的辐射作用明显增强,远高于山脊火。因此山谷ROS、火焰长度和火焰辐射热量远高于山脊火。