-
森林火灾是一种突发性强、破坏性大的自然灾害，能烧毁森林资源、破坏全球生态平衡、危及人类生命财产安全，还需投入大量人力、物力和财力进行防范和扑救，其造成的环境污染问题和生态系统危害需要几十年甚至上百年才能恢复[1]。目前世界各国都在探索有效的森林防火防范方法，其中建设科学合理的林道网并发挥其防火隔离带作用是一种便捷、高效防治森林火灾的方法。国外林道网建设研究从18世纪开始，工业革命的兴起和木材需求使奥地利、德国、澳大利亚等发达国家在林区修建林道，到20世纪进入了飞速发展阶段，20世纪60年代汽车的出现推动林区道路大规模建设。亚洲地区许多国家的林区道路建设起步较晚，韩国林道建设始于20世纪中期，末期进入全面加速发展阶段。近年来，欧洲林业发达国家的林道密度已达到40 m/hm2，甚至100 m/hm2。当前林道研究主要聚焦于路网规划[2-3]，林道综合费用[4-5]，林道密度[6-7]，自动选线[8]，林道损伤[9]，林道生态[10]，林道防火和救援功能[11]，以及林道宽度计算[12]等。地中海地区许多森林火灾频发的国家已发展到利用高密度的生物防火林带、林区道路和防火线将森林网格化以减小森林区域面积和降低火险达到防火目的，再配合较强的空中、地面灭火力量实现高水平防火和灭火。同时林区道路还具有阻隔森林火灾蔓延的隔离带功能，为林火巡护监测提供交通条件，以及发生火灾时快速运送扑救人员和装备，缩短到达火场时间，实现及时扑救森林火灾。
我国林道网研究始于20世纪60年代初，数量、质量与发达国家相比差距都很大，初期主要满足木材采运，研究内容聚焦于平均集材距离和林道密度。当前我国林道网处于建设阶段，水平相对落后，存在林道网总量不足，林道等级和密度低，缺乏林道网络优化配置，忽略林道宽度差异设置，忽视林道阻火功能和生态效益等诸多问题[13]，仍停留在经济效益阶段。罗康[14]分析北京林道密度并提出须增建道路和林区道路。王国华[15]分析广州市林区道路密度，发现该市没有达到标准林道密度。邱荣祖等[16]分析林道网阻隔林火蔓延的效益和森林火灾损失费用的关系。随着国家政策、法律和法规的颁布和完善，在地方各级政府及林业主管部门的共同努力下，我国林道网络防火阻隔系统的建设取得较大进展。截止到2005年底，全国建有林区防火道路8.95 × 105 km，其中林区等级公路4.62 × 105 km，工程阻隔带5.04 × 105 km，生物防火林带1.503 × 106 km，其他隔离带1.05 × 105 km[17]。林区防火道路、工程阻隔带和生物防火林带的发展，为森林防火工作稳步发展奠定了良好基础，森林火灾隐患显著减少，一定程度上提高了我国森林火灾的综合防控能力。
目前应用火灾动态仿真模拟软件（fire dynamic simulation，FDS）的扩展软件Thunderhead Engineering PyroSim模拟森林火灾的研究较少。Sun 等[18]比较FDS和CCWM（Clark coupled wildfire model）两个模型的模拟效果，表明FDS在较粗网格水平仍能与火灾真实数据吻合，CCWM则对表面垂直网格大小以及能量释放更为敏感。国内学者尚超等[19]提出树木建模方法。关祥毅等[20]应用FDS模拟樟子松（Pinus sylvestris var. mongolica）树冠火燃烧的速度矢量场和温度场。辛喆等[21]基于FDS建立草原火灾模型并分析影响草原火灾的因素及其蔓延规律。李世友[22]应用FDS研究华山松(P. armandii）林、云南松（P. yunnanensis）林细小可燃物层燃烧初始蔓延速度。FDS扩展软件WFDS广泛应用于森林和草原火灾研究[23-27]。在我国林道网络建设较为薄弱的背景下，依据林区公路设计标准，科学指导林道规划和配置，提高林道密度和等级，合理设置林道宽度，使林道宽度及其网络既能满足生产、经济、交通、旅游等方面要求，兼具防火隔离带功能，阻隔森林火灾蔓延和提供应急救援保障，提升森林火灾的防控能力。
云南松是云南省主要植被树种之一，富含油脂，极易点燃蔓延。2006年云南省安宁市“3·29”重大森林火灾的火场植被以云南松、地盘松（P. yunnanensis
var. pygmaea ）、栎类 （Quercus spp.）灌木为主，平均坡度约25°，强风作用下频发飞火，其距离最长达1 000 m。2011年云南省大理县剑川县“3·2”森林火灾的林相主要为云南松中幼林及杂草，火场平均坡度约30°。云南省森林覆盖率54.64%，其中云南松、华山松、草灌等易燃植被广布，火灾多发生在坡度25°以上山区，高山峡谷内风速多变，在风力作用下易发生极端火行为。据云南省减灾年鉴及中国森林火灾案例统计，云南松林火灾中出现飞火的概率较高，但飞火的距离从十多米到上百米，甚至上千米，用林区道路防范极端火行为不符合实际。目前我国林区道路的防火阻隔功能研究大多在火灾发生后，依据实地调查或经验，以及依据林区公路路线设计规范设置统一林道宽度，忽视不同坡度、不同风速、树种差异对林道防火阻隔功能的需求。选择云南省云南松林火灾区域的地形数据和气象数据，应用火灾动态仿真模拟软件PyroSim构建不同坡度、不同风速下的微观森林模型，通过仿真模拟燃烧确定实现林道阻火隔离功能需求的宽度。 -
PyroSim是美国国家标准与技术研究院（National Institute of Standards and Technology，NIST）研发的专用于FDS的软件，该软件以流体动力学为理论依据，将研究空间分成多个小的三维矩形控制体或计算单元，每个计算单元的气体密度、速度、温度、压力和组分浓度用质量守恒、动量守恒、能量守恒和组分守恒的偏微分方程进行近似有限差分，通过对同一网格使用有限体积技术计算热辐射、烟气流动、热传递过程，追踪预测火灾气体的产生和移动，并结合燃料特性计算火灾的增长和蔓延。模拟求解后可获得设置测点的温度、CO浓度、CO2浓度、能见度、烟气运动等数据，并提供火灾动态模拟图形界面。
-
（1） t2火模型
国内外常用稳态和非稳态两类模型描述火灾增长阶段热释放速率（heat release rate，HRR）的变化规律。稳态增长模型指在火灾发展过程中，热释放速率始终保持一个恒定值，将整个火灾过程理想化。非稳增长模型是火灾发展过程中热释放速率随时间发生变化，主要有t2火模型，MRFC（multi-room-fire-code）模型，FFB（forschungsstelle fuer brandschutztechnik）指数模型[28]，完整发展模型CFAST（consolidate fire and smoke transport）[29−30]等。火灾实际增长过程中热释放速率是发展变化的，可采用t2火模型描述火灾增长，该模型热释放速率变化按时间的平方规律增长。通常火灾研究不考虑火灾酝酿期，从火灾出现即算有效燃烧，公式为：
$$ Q = \alpha {t^2} $$ (1) 式中：Q为热释放速率（kW）；α为火灾增长系数（kW/s2）；t为着火后的时间（s）。
不同可燃物的燃烧热值、热释放速率和增长系数存在较大差异，火灾增长类型按增长速率分为慢速、中速、快速、超快速4种，其增长系数值以及典型可燃材料见下表1[31]。实际应用中依据可燃物材料特性选择对应参数值，森林火灾的材料类别属于木质类，选择中速增长模型，设置运行时间为300 s。
表 1 4种标准t2火
Table 1. Four kinds of t2 standard fire
火灾增长类型
Fire growth type增长系数
Growth coefficient/(kW·s−2)达到1 MW需要时间
Time of reaching 1 MW/s典型可燃材料
Typical combustible material超快速 Ultra-fast speed 0.187 800 75 油池火、易燃的装饰家具、轻的窗帘
Pool fire, flammable decorative furniture, light curtains快速 Fast speed 0.046 890 150 装满东西的邮袋、塑料泡沫、叠放的木架
A mail bag full of things, plastic foam, stacked wooden shelves中速 Medium speed 0.011 270 300 棉与聚酯纤维弹簧床垫、木制办公桌
Cotton and polyester spring mattress, wooden desk慢速 Slow speed 0.002 931 600 厚重木制品
Heavy woodwork（2） 燃烧速率
燃烧速率是指单位时间内可燃物燃烧的质量，决定火灾的热释放速率，公式为：
$$ G={m}/{t} $$ (2) 式中：G为燃烧速率（kg/s）；m为可燃物燃烧的质量（kg）；t为燃烧时间（s）。
-
树木一般由树冠、树干和树根组成，森林火灾最常见的类型为地表火和树冠火，地下火较为少见。森林火灾主要燃烧地表有效可燃物，即树木地上的树干和树冠两部分。树干模型简化为圆柱体。树冠是树木可燃成分的主要组成部分，树冠火是森林燃烧危害最严重的类型。不同树种的树冠存在一定差别，按照前人研究将其形态分为锥形、纺锤形和半球形三类[21]。树种间树高、胸径、冠幅等形态和指标存在差异进而影响森林燃烧性，文中研究云南松林区道路的防火阻隔功能，忽略树木差异，选择相同参数构建树木模型。以我国云南省分布较广的云南松中龄林为建模对象，选择锥形树冠，取树高H = 10 m，树冠高h = 2 m，冠幅L = 2 m（树木模型必须与0.5网格贴合，否则物体若不规则，可能会导致数据无法读取或不精确），胸径D = 0.43 m，云南松可燃物相关参数见表2[22]。
表 2 云南松可燃物参数
Table 2. Forest fuel parameters of P. yunnanensis
参数
Parameter燃烧热
Combustion heat/
(kJ·g−1)参考温度
Reference
temperature/℃密实度
Density/
(kg·m−3)反应热
Reaction heat/
(kJ·g−1)热传导率
Heat conductivity/
(W·km−1)比热容
Specific heat capacity/
(J·kg−1)树干 Trunk 8 260 900 3 0.28 3.0 树叶 Leaf 13.28 260 15 1.5 0.05 3.2 -
影响森林火灾发生、发展的因素众多，以PyroSim软件建立火场全真模型模拟测算具体燃烧数值计算量大、时间消耗长。其中模型大小、网格精度、设置时长等决定仿真模拟燃烧运行时间，通常选择最能代表研究区域的一个地段或区域为对象构建简化模型。通过调查现行《LYJ 113—1992林区公路路线设计规范》[32]，研究区森林火灾蔓延速率和风速等基础数据，以及云南松森林火灾案例，建立符合现实情况的林区道路微观梯形山体简化模型。如图1所示，四周网格设置全部开启以模拟野外森林实况的自然环境，设置树木4排5列的小片森林，横间距为2.5 m，前3排纵间距为2.5 m，3排与4排中间设置不同宽度的道路模拟林道阻隔森林火灾蔓延，风向设置为顺风，即山体模型最下方面向山顶方向的进风口。依据林区公路路线设计规范，道路宽度选取低于标准2 m、最低标准4.5 m、加宽6 m三个类别，坡度以云南松分布区域、火灾多发点的地形数据为参考，选择坡度25°、35°、45° 三种类别，以及西南林区低风速 1 m/s、中风速 2.5 m/s、高风速 4.5 m/s[33]共构建27种工况。
图 1 模型整体外观示意图
Figure 1. Overall appearance diagram of model
-
沿微观山地模型中轴线，分别在上坡位树最左端、林道两侧、下坡位每排树与树的中央共设置纵向5列测点，每列间隔1 m共设置16个测点，总计最多90个测点（随山体坡度增加，测点被遮挡逐渐减少）探测模型中各个位置的温度、CO浓度、CO2浓度等指标，如图2所示。并在整个模型的纵向中轴、每排树木和道路的两侧设置温度、可见度和粒子移动情况共计7片切片，如图3所示。切片的功能是在轴对齐平面上测量压力、速度、温度等数据。
图 2 监测点设置示意图
Figure 2. Schematic diagram of monitoring point setting
图 3 切片监测设置示意图
Figure 3. Schematic diagram of slice monitoring setting
-
森林模型设置为4排5列的微观梯度小片森林，即使全部成片燃烧也难以达到真实森林火灾燃烧的林火强度，以致林道网作为防火隔离带的作用难以完全体现，且用PyroSim仿真模拟森林火灾的研究极少，无法查得发生云南松森林火灾时单位面积的热释放速率。基于上述模型和相关基础数据进行模型调试，以期选择合适的火源功率模拟中等强度的林火蔓延。
调用 FDS/Smokeview 运行模型，模型越大、越复杂、网格越细、设置火灾模拟时间越长，运算时间越长。一般情况下复杂模型先用粗网格短时间试运行，依据粗略模拟结果预判模型是否符合预期。模型试运行参数见表3。
表 3 模型试运行参数
Table 3. Model testing running parameters
工况
Working condition网格
Mesh/m3火源面积
Fire source area/m2热释放速率
Heat release rate/(kW·m−2)蔓延速度
Spreading rate/(m·min−1)1 0.5 1 4.464 × 104 0.2 2 0.5 12 2.0 × 104 1.0 3 0.5 12 4.6 × 104 8.0 试运行工况1：设置点燃树木的火源功率，取α = 0.496，t = 300 s，Q = 4.464 × 104 kW，并设置火源1 m2，运行后发现树木燃烧时间过长，热释放速率较小，蔓延速度慢，无法模拟森林火灾全盛时期的蔓延情况，模拟效果差。
试运行工况2：设置长12 m、宽1 m、HRR为2 × 104 kW/m2的带状火源模拟燃烧充分的成片火源，布置于模型最右侧下坡位置。燃烧5 min后，火势开始蔓延，但蔓延速度较低，模拟效果差。
试运行工况3：调高HRR至4.6 × 104 kW/m2，其他指标不变，并设置树冠材料HRR至4.6 × 104 kW/m2，火势在1 min左右向上坡方向蔓延8 m，符合实际中等强度的森林火灾蔓延规律。试运行的火蔓延速率结果见表3。
经过上述数十次模型试运行并结合参考文献，设置长12 m、宽1 m、热释放速率为4.6 × 104 kW/m2的带状火源时，蔓延速率与学者王秋华[34]研究云南松林火灾蔓延速率10 ~ 80 m/min的结果相一致，以此火源模拟中等强度林火分析不同林道宽度在不同坡度、不同风速下的防火阻隔功能。在确保模型运行结果精度和节约时间的前提下，模型网格划分为0.5 m × 0.5 m × 0.5 m，工况设置见表4。
表 4 模拟工况表
Table 4. Simulated working conditions
工况号
Working condition No.类别
Class工况号
Working condition No.类别
Class工况号
Working condition No.类别
Class工况号
Working condition No.类别
Class1 25°-2.0 m-1.0 m/s 8 25°-6.0 m-2.5 m/s 15 35°-4.5 m-4.5 m/s 22 45°-4.5 m-1.0 m/s 2 25°-2.0 m-2.5 m/s 9 25°-6.0 m-4.5 m/s 16 35°-6.0 m-1.0 m/s 23 45°-4.5 m-2.5 m/s 3 25°-2.0 m-4.5 m/s 10 35°-2.0 m-1.0 m/s 17 35°-6.0 m-2.5 m/s 24 45°-4.5 m-4.5 m/s 4 25°-4.5 m-1.0 m/s 11 35°-2.0 m-2.5 m/s 18 35°-6.0 m-4.5 m/s 25 45°-6.0 m-1.0 m/s 5 25°-4.5 m-2.5 m/s 12 35°-2.0 m-4.5 m/s 19 45°-2.0 m-1.0 m/s 26 45°-6.0 m-2.5 m/s 6 25°-4.5 m-4.5 m/s 13 35°-4.5 m-1.0 m/s 20 45°-2.0 m-2.5 m/s 27 45°-6.0 m-4.5 m/s 7 25°-6.0 m-1.0 m/s 14 35°-4.5 m-2.5 m/s 21 45°-2.0 m-4.5 m/s 注：各工况以坡度(°)-路宽(m)-风速(m/s)标记。Notes: working condition is marked by slope (°)-road width (m)-wind speed (m/s). -
坡度25°时不同林道宽度对不同风速条件下林火蔓延的阻隔模拟见图4。随着风速增加，燃烧加剧，火势增强。中等强度的林火均蔓延越过宽度为2 m的林道，点燃林道另一侧可燃物的时间约70 s。林道宽度4.5、6 m均成功阻隔3种不同风速条件下中等强度的林火蔓延。
图 4 坡度25°时不同林道宽度对不同风速条件下林火蔓延的阻隔模拟
Figure 4. Blockade simulation on various widths of forest roads to fire spread under different wind speeds at slope of 25°
坡度35°时不同林道宽度对不同风速条件下林火蔓延的阻隔模拟见图5，林道宽度2 m无法满足该情况的防火隔离需求，点燃林道另一侧可燃物的时间约65 s。林道宽度4、6 m可满足风速小于4.5 m/s时中等火强度的防火阻隔需求。
图 5 坡度35°时不林道宽度对不同风速条件下林火蔓延的阻隔
Figure 5. Blockade simulation on various widths of forest roads to fire spread under different wind speeds at slope of 35°
坡度45°时不同林道宽度对不同风速条件下林火蔓延的阻隔模拟见图6。林道宽度2、4 m均无法满足此情况的防火阻隔要求，林道宽度6 m可满足风速小于4.5 m/s时该中等强度森林火灾的阻隔要求。由风速4.5 m/s的林火蔓延模拟图可看出火焰已经相当接近于上坡位的树木顶端，一旦气候环境出现极端情况，林道宽度6 m也将失去防火隔离作用。
图 6 坡度45°时林道宽度对不同风速条件下林火蔓延的阻隔
Figure 6. Blockade simulation on various widths of forest roads to fire spread under different wind speeds at slope of 45°
-
温度数据采集于上坡位树木模型顶端处探测点。由图7可看出，林道宽度2 m、坡度大于25°时，上坡位树木顶端的温度均接近或大于240 ℃，达到云南松针叶着火点，引燃上坡位树木，林区道路的防火阻隔作用失效。当林道宽度为4.5 m时，坡度小于35°情况下，温度仅达到150 ℃左右，无法达到云南松针叶的燃点，林道阻隔功能有效。当坡度45°、风速为4.5 m/s时，树端温度达到250 ℃左右，点燃上坡位云南松针叶，林道防火阻隔作用失效。林道宽度6 m时，坡度小于35°的情况即便风速达到4.5 m/s，上坡位树端测点的温度约200 ℃，无法达到云南松针叶的燃点，可实现林道的防火阻隔作用。在坡度45°、风速4.5 m/s的情况下，测点温度有一小段时间峰值达到250 ℃，若火势持续下去，极有可能出现林道的阻火功能失效。
图 7 上坡位树木模型顶端温度
Figure 7. Temperature at the top of tree model in uphill position
-
云南松林模型燃烧的HRR见图8所示。森林火灾发生到全盛时期释放的热量在3.0 × 105 ~ 9.0 × 105 kW之间，极端情况下可超过1.0 × 106 kW，HRR与坡度、风速成正比，随坡度、风速增加，蔓延速率加快，热释放速率增大。
图 8 热释放速率曲线图
Figure 8. Curves of heat release rate
-
云南松林模型可燃物载量相同，林道宽度对火场中心温度没有影响，选择坡度35°、林道宽度4.5 m时不同风速进行对比分析，如图9所示。在该中等强度的森林火灾中，火场中心温度达到700 ℃以上，极端条件时甚至高于1 200 ℃。随风速越大，火场中心温度增加明显。
图 9 火场中心温度
Figure 9. Temperature in fire center
-
云南松林模型可燃物载量相同，林道宽度对燃烧速率没有影响，选取坡度35°、林道宽度4.5 m时不同风速对比分析，如图10所示。中等强度的森林火灾1 min左右燃烧速率可达到50 kg/s，短时间内即可烧毁大量森林资源，对生态环境的破坏力非常大。从可燃物燃烧速率也证明风速可以影响森林燃烧强度，风速越大，可燃物燃烧速率增加。
图 10 燃烧速率
Figure 10. Burning rate
-
为较好地实现林道的防火阻隔功能，运用火灾动态仿真模拟软件PyroSim构建不同坡度、风速和林道宽度的云南松林微观模型，设置带状火源模拟中等强度的森林火灾进行数值仿真模拟，探究不同林道宽度的防火阻隔功能，得到以下结论。
（1）林道宽度为2 m时，上坡位树木顶端的温度均大于240 ℃，可引燃上坡位云南松针叶，林道阻火作用均失效，无法满足防火隔离需求。
（2）林道宽度为4.5 m时，坡度小于35°，风速小于4.5 m/s情况，上坡位树木顶端的温度均小于220 ℃，不能引燃上坡位云南松针叶，能满足此类情况防火阻隔需求。坡度大于45°，风速超过4.5 m/s则将出现防火隔离带失效的风险。林区道路宽度4.5 m的防火阻隔性能一般，坡度较大时容易失效。
（3）林道宽度为6 m时，坡度小于45°的上坡位树木顶端的温度大部分情况在200 ℃，林区道路拥有较强的防火隔离功能，够较好阻隔中等强度的森林火灾火势蔓延。当风速大于4.5 m/s的恶劣环境时，火势仍然可能跨越过林区道路继续蔓延。
（4）中等强度森林火灾的火场中心温度在700 ~ 1 200 ℃左右，热释放功率可达到3.0 × 105 ~ 9.0 × 105 kW，蔓延山坡10 m只需约1 min，可燃物燃烧速率60 s可达到40 kg/s以上。
-
林火蔓延受地形、气象、可燃物等因素影响，坡度、风速增加，林火蔓延速率增加。林道阻火功能能否实现，需要充分调研当地的风速，特别是防火期的最高风速和平均风速，以此设置合理的林道宽度阈值。常年处于高风速、坡度较大、树种高大、载量高的林区建议增加林道宽度至8 ~ 10 m，以便取得较好的防火效益。现阶段防火林道宽度设置多以灾后调研数据或经验确定，以PyroSim火灾模拟软件构建模型仿真模拟林道的阻火功能，可为防火林道宽度和密度设置提供一种全新思路，当需要在森林中设置林道实现阻火隔离功能时，按照地形、风速、林木可燃物燃烧性能参数建立仿真模型进行模拟燃烧，可得出满足经济效益和具备阻火隔离作用的林道宽度和密度指标。运用PyroSim构建模型模拟火灾现场，原则上模型越精准，模拟结果越可靠准确，但高精度模型需要高配置的计算机承载运行计算过程，否则不能计算出高可靠度的数据，还会因大功率火灾、复杂模型模拟计算量太大而使软件停止，根据电脑实际配置构建精确适度、数据量适中的模型进行仿真模拟所得到的结果相对真实可靠。
Study on forest road of fireproof blockade functions based on PyroSim
-
摘要:
目的 林道网络及其宽度是林业生产、交通和旅游的需求，在发挥其最大经济效益的同时，兼顾林火巡护监测、应急救援，以及发生火灾时快速运送扑救人员和装备，阻隔森林火灾蔓延和发挥防火隔离功能。 方法 运用火灾动态仿真模拟软件PyroSim分别以坡度25°、35°、45°，风速1、2.5、4.5 m/s，林道宽度2、4.5、6 m构建云南松林微观山体模型，施加长12 m、宽1 m、热释放速率4.6 × 104 kW/m2的带状火源，模拟中等强度林火研究不同林道宽度的防火阻隔功能。 结果 林道宽度2 m无法满足防火隔离需求。林道宽度4.5 m的防火阻隔功能一般，当坡度超过35°、风速大于4.5 m/s时易失效。林道宽度6 m拥有较强的防火隔离功能，能较好阻隔中等强度的林火蔓延，可满足坡度小于45°、风速小于4.5 m/s的情况，当风速大于4.5 m/s时该林道宽度的隔离功能可能失效。中等强度的森林火灾火场中心温度在700 ~ 1 200 ℃左右，热释放功率可达3.0 × 105 ~ 9.0 × 105 kW，蔓延山坡10 m只需约1 min，可燃物燃烧速率60 s可达到40 kg/s以上。 结论 云南省云南松林区道路宽度建议设置为6 m，能够取得较好的防火效益。常年处于高风速，坡度大于45°，树种高大的林区建议增加林道宽度至8 ~ 10 m。此外本研究应用PyroSim建立火场模型仿真模拟林道阻火功能是可行的，可为林区道路建设的合理设计、确定林道密度和宽度的阻火功能指标与经济效益相适应的最大限度阈值、降低森林火灾破坏力具有一定的指导意义。 Abstract:Objective Forest road network and its width is the basic demands of forestry production, transportation, and tourism. While giving maximum economic benefits, its function on forest fire patrol monitoring, emergency rescue, rapid transport of firefighters and equipment when happening fires also should be taken into account, as well as the function of block forest fire spreading and forest fire isolation. Method The micro-mountains model of Pinus yunnanensis was constructed using fire dynamic simulation software PyroSim with parameters including slopes of 25°, 35° and 45°, wind speed of 1 , 2.5 and 4.5 m/s, and forest-road widths of 2, 4.5 and 6 m. And with the application of belt fire source with 12 m in length and 6 m in width, and heat release rate of 4.6 × 104 kW/m2, which simulated moderate-intensity forest-fire, we aimed to study fire-isolating functions at different forest-road widths. Result The forest road width of 2 m failed to meet the need for fire protection and isolation. And the one with 4.5 m achieved relatively poor fire-isolating function performance. And the one with more than 35° slope and greater than 4.5 m/s wind speed was easy to lose its function. And 6 m width road got a strong fire isolation function. Comparing with the one that can better isolate medium-intensity forest-fire spreading condition can meet the one with less than 4.5 m/s wind speed and less than 45° slope. When wind speed reached more than 4.5 m/s, the fireproof blockade functions of forest-road width of 6 m may be lost. Generally, the central temperature of moderate-intensity forest fire was about 700−1 200 ℃. And corresponding heat release rate power can reach 3.0 × 105 − 9.0 × 105 kW. And it only took about 1 min to spread 10 m on the hill. And the burning rate of combustibles can reach above 40 kg/s in the 60 s. Conclusion It is suggested to build the road with 6 m width in P. yunnanensis district in Yunnan Province of southwestern China, which can achieve better fire-prevention benefits. And road with width of 8−10 m is highly recommended in the forest area, where the stands with higher tree height all year round suffering high wind speed and slope greater than 45°. It also approve that the feasibility of simulating forest-road fire-isolating function by constructing the fire-site model with PyroSim, the results in this study have certain guiding significance on the rational designs of forest-road constructions, such as the determinations of maximum threshold values of fire prevention functions of forest-road with indexes of densities and widths, so can achieve certain economic benefits, as well as reductions of destructive powers of forest fires. -
Key words:
- Pinus yunnanensis /
- forest-road width /
- heat release rate /
- blockade function /
- PyroSim
-
图 4 坡度25°时不同林道宽度对不同风速条件下林火蔓延的阻隔模拟
Figure 4. Blockade simulation on various widths of forest roads to fire spread under different wind speeds at slope of 25°
图 5 坡度35°时不林道宽度对不同风速条件下林火蔓延的阻隔
Figure 5. Blockade simulation on various widths of forest roads to fire spread under different wind speeds at slope of 35°
图 6 坡度45°时林道宽度对不同风速条件下林火蔓延的阻隔
Figure 6. Blockade simulation on various widths of forest roads to fire spread under different wind speeds at slope of 45°
表 1 4种标准t2火
Table 1. Four kinds of t2 standard fire
火灾增长类型
Fire growth type增长系数
Growth coefficient/(kW·s−2)达到1 MW需要时间
Time of reaching 1 MW/s典型可燃材料
Typical combustible material超快速 Ultra-fast speed 0.187 800 75 油池火、易燃的装饰家具、轻的窗帘
Pool fire, flammable decorative furniture, light curtains快速 Fast speed 0.046 890 150 装满东西的邮袋、塑料泡沫、叠放的木架
A mail bag full of things, plastic foam, stacked wooden shelves中速 Medium speed 0.011 270 300 棉与聚酯纤维弹簧床垫、木制办公桌
Cotton and polyester spring mattress, wooden desk慢速 Slow speed 0.002 931 600 厚重木制品
Heavy woodwork
下载: 导出CSV
表 2 云南松可燃物参数
Table 2. Forest fuel parameters of P. yunnanensis
参数
Parameter燃烧热
Combustion heat/
(kJ·g−1)参考温度
Reference
temperature/℃密实度
Density/
(kg·m−3)反应热
Reaction heat/
(kJ·g−1)热传导率
Heat conductivity/
(W·km−1)比热容
Specific heat capacity/
(J·kg−1)树干 Trunk 8 260 900 3 0.28 3.0 树叶 Leaf 13.28 260 15 1.5 0.05 3.2
下载: 导出CSV
表 3 模型试运行参数
Table 3. Model testing running parameters
工况
Working condition网格
Mesh/m3火源面积
Fire source area/m2热释放速率
Heat release rate/(kW·m−2)蔓延速度
Spreading rate/(m·min−1)1 0.5 1 4.464 × 104 0.2 2 0.5 12 2.0 × 104 1.0 3 0.5 12 4.6 × 104 8.0
下载: 导出CSV
表 4 模拟工况表
Table 4. Simulated working conditions
工况号
Working condition No.类别
Class工况号
Working condition No.类别
Class工况号
Working condition No.类别
Class工况号
Working condition No.类别
Class1 25°-2.0 m-1.0 m/s 8 25°-6.0 m-2.5 m/s 15 35°-4.5 m-4.5 m/s 22 45°-4.5 m-1.0 m/s 2 25°-2.0 m-2.5 m/s 9 25°-6.0 m-4.5 m/s 16 35°-6.0 m-1.0 m/s 23 45°-4.5 m-2.5 m/s 3 25°-2.0 m-4.5 m/s 10 35°-2.0 m-1.0 m/s 17 35°-6.0 m-2.5 m/s 24 45°-4.5 m-4.5 m/s 4 25°-4.5 m-1.0 m/s 11 35°-2.0 m-2.5 m/s 18 35°-6.0 m-4.5 m/s 25 45°-6.0 m-1.0 m/s 5 25°-4.5 m-2.5 m/s 12 35°-2.0 m-4.5 m/s 19 45°-2.0 m-1.0 m/s 26 45°-6.0 m-2.5 m/s 6 25°-4.5 m-4.5 m/s 13 35°-4.5 m-1.0 m/s 20 45°-2.0 m-2.5 m/s 27 45°-6.0 m-4.5 m/s 7 25°-6.0 m-1.0 m/s 14 35°-4.5 m-2.5 m/s 21 45°-2.0 m-4.5 m/s 注：各工况以坡度(°)-路宽(m)-风速(m/s)标记。Notes: working condition is marked by slope (°)-road width (m)-wind speed (m/s).
下载: 导出CSV
-
[1] 叶兵. 国内外森林防火技术及其发展趋势[D]. 北京: 中国林业科学研究院, 2000. Ye B. The present situation and development tendency of domestic and abroad technical level of forest fire prevention[D]. Beijing: China Academy of Forestry Research Institute of Forestry, 2000. [2] Picchio P, Pignatti G, Marchi E, et al. The application of two approaches using GIS technology implementation in forest road network planning in an Italian mountain setting[J]. Forests, 2018, 9(5): 277. doi: 10.3390/f9050277 [3] Buday E, Akay A E. Evaluation of forest road network planning in landslide sensitive areas by GIS-based multi-criteria decision making approaches in Ihsangazi Watershed, northern Turkey[J]. Šumarski List, 2019, 143(7−8): 325−336. [4] Akay O A, Akgul M, Demir M, et al. Analysis of factors associated with the amount of forest road reconstruction activity in Turkey: autoregressive distributed lag modelling approach[J]. Forest Ecology and Management, 2020, 458: 117800. [5] Keramati A, Lu P, Sobhani A, et al. Impact of forest road maintenance policies on log transportation cost, routing, and carbon-emission trade-offs: oregon case study[J]. Journal of Transportation Engineering, Part A: Systems, 2020, 146(5): 04020028. doi: 10.1061/JTEPBS.0000335 [6] Çalışkan E, Karahalil U. Evaluation of forest road network and determining wood extraction system using GIS: a case study in Anbardağ planning unit[J]. Sumarski List, 2017, 141(3−4): 163−171. doi: 10.31298/sl.141.3-4.6 [7] Petković V, Potočnik I. Planning forest road network in natural forest areas: a case study in northern bosnia and herzegovina[J]. Croatian Journal of Forest Engineering, 2018, 39(1): 45−56. [8] Tan J. Planning a forest road network by a spatial data handling-network routing system[J]. Acta Forestalia Fennical, 1992, 227: 26−41. [9] Hrůza P, Mikita T, Tyagur N, et al. Detecting forest road wearing course damage using different methods of remote sensing[J]. Remote Sensing, 2018, 10(4): 492. doi: 10.3390/rs10040492 [10] Fallahchai M M, Haghverdi K, Mojam M S, et al. Ecological effects of forest roads on plant species diversity in Caspian forests of Iran[J]. Acta Ecologica Sinica, 2018, 38(3): 255−261. doi: 10.1016/j.chnaes.2017.08.002 [11] Andrea L, Cristiano F, Fabio F, et al. Forest road planning, construction and maintenance to improve forest fire fighting: a review[J]. Croatian Journal of Forest Engineering, 2019, 40(1): 207−219. [12] Schöning D R, Bachmann D A, Allgöwer D B. GIS-based framework for wildfire risk assessment final report MINerve2[R]. Zurich:University of Zurich,1997. [13] 邱荣祖, 方金武, 詹正宜, 等. 林道网理论研究的历史回顾与展望[J]. 福建林学院学报, 2000, 20(4):370−374. doi: 10.3969/j.issn.1001-389X.2000.04.022 Qiu R Z, Fang J W, Zhan Z Y, et a1. A review and prospect for the research of forest road network[J]. Journal of Fujian College of Forestry, 2000, 20(4): 370−374. doi: 10.3969/j.issn.1001-389X.2000.04.022 [14] 罗康. 基于GIS技术的北京市森林防火资源布局的分析[D]. 株洲: 中南林学院, 2003. Luo K. GIS-based layout analysis of forest fire protection in Beijing[D]. Zhuzhou: Central South Forestry College, 2003. [15] 王国华. 森林地面防火资源调控技术研究[D]. 株洲: 中南林业科技大学, 2009. Wang G H. The technology research of ground-based resources for forest fire control and adjust[D]. Zhuzhou: Central South University of Forestry and Technology, 2009. [16] 邱荣祖, 邓志华. 林道网的森林火灾防治效益评价[J]. 森林防火, 1998(4):21−22. Qiu R Z, Deng Z H. Evaluation on forest road network of forest fire prevention benefit[J]. Forest Fire Prevention, 1998(4): 21−22. [17] 张国红. 我国森林防火阻隔系统建设现状与思考[J]. 林业资源管理, 2008(4):15−18. doi: 10.3969/j.issn.1002-6622.2008.04.005 Zhang G H. Present situation and thoughts on the development of the forest fire prevention and obstruction system in China[J]. Forest Resources Management, 2008(4): 15−18. doi: 10.3969/j.issn.1002-6622.2008.04.005 [18] Sun R Y, Jenkins M A, Krueger S K, et al. An evaluation of fire-plume properties simulated with the fire dynamics simulator (FDS) and the Clark coupled wildfire model[J]. Canadian Journal of Forest Research, 2006, 36(11): 2894−2908. doi: 10.1139/x06-138 [19] 尚超, 王克印, 黄海英, 等. 基于的火灾树木建模研究[J]. 消防科学与技术, 2013, 32(9):1030−1033. doi: 10.3969/j.issn.1009-0029.2013.09.030 Shang C, Wang K Y, Huang H Y, et al. Forest trees in fire model research based on the Pyrosim[J]. Fire Science and Tecnnology, 2013, 32(9): 1030−1033. doi: 10.3969/j.issn.1009-0029.2013.09.030 [20] 关祥毅, 刘少刚. 樟子松树冠火燃烧场数学模型的研究[EB/OL]. 北京: 中国科技论文在线(2013−01−25) [2019−09−25]. http://www.paper.edu.cn/releasepaper/content/201301-1040. Guan X Y, Liu S G. Research on the mathematical model of mongolian Scots pine crown fire burning area[EB/OL]. Beijing: Sciencepaper Online (2013−01−25) [2019−09−25]. http://www.paper.edu.cn/releasepaper/content/201301-1040. [21] 辛喆, 王顺喜, 云峰, 等. 基于火灾模拟软件的草原火灾蔓延规律数值分析[J]. 农业工程学报, 2012, 8(10):54−57. Xin Z, Wang S X, Yun F, et al. Numerical analysis on spreading laws of grassland fire based on fire dynamics simulator (FDS)[J]. Transactions of the Chinese Society of Agricultural Engineering, 2012, 8(10): 54−57. [22] 李世友. 滇中森林可燃物燃烧性及林火行为研究[D]. 北京: 中国林业科学研究院, 2014. Li S Y. Study on the combustibility of forest fuels and fire behavior in central Yunnan Province[D]. Beijing: China Academy of Forestry Research Institute of Forestry, 2014. [23] Morvan D, Hoffman C, Rego F, et a1. Numerical simulation of the interaction between two fire fronts in grassland and shrubland[J]. Fire Safety Journal, 2011, 46(8): 469−479. doi: 10.1016/j.firesaf.2011.07.008 [24] Overholt K J, Cabrera J, Kurzawski A, et a1. Characterization of fuel properties and fire spread rates for little bluestem grass[J]. Fire Technology, 2014, 50(1): 9−38. doi: 10.1007/s10694-012-0266-9 [25] Miloua H. Fire behavior characteristics in a pine needle fuel bed in northwest Africa[J]. Journal of Forestry Research, 2019, 30(3): 959−967. doi: 10.1007/s11676-018-0676-8 [26] Mell W E, Manzello S L, Maranghides A. Numerical modeling of fire spread through trees and shrubs[J]. Forest Ecology and Management, 2006, 234(1): S82. [27] Mell W E, Jenkins M A, Gould J, et a1. A physics-based approach to modeling grassland fires[J]. International Journal of Wildland Fire, 2007, 16(1): 1−22. doi: 10.1071/WF06002 [28] 程远平. 火灾过程中火源热释放速率模型及其实验测试方法[J]. 火灾科学, 2002, 11(2):70−74. doi: 10.3969/j.issn.1004-5309.2002.02.002 Cheng Y P. The model and experimental testing method of heat release rate of fuel during the development of fire[J]. Fire Safety Science, 2002, 11(2): 70−74. doi: 10.3969/j.issn.1004-5309.2002.02.002 [29] Jones W W, Forney G P, Peacock R D, et a1. A technical reference for CFAST: an engineering tool for estimating fire and smoke transport [R]. Gaithersburg: Building and Fire Research Laboratory, 2000. [30] Jones W W, Peacock R D, Forney G P, et al. Verification and validation of CFAST: a model of fire growth and smoke spread[R]. Washington: National Institute of Standards and Technology, 2004. [31] Nelson H E. An engineering analysis of the early stages of fire development-the fire at the dupont plaza hotel and casino, December 31, 1986[R/OL].(1987−05−29)[2019−01−22]. https://doi.org/10.6028/NBS.IR.87−3560. [32] 中华人民共和国林业部.林区公路路线设计规范: LYJ 113—1992 [S].北京: 中国标准出版社, 1997. National Forestry Administration. Design specification of road for forest highway:LYJ 113−1992[S]. Beijing: Standards Press of China,1993. [33] 张志斌, 杨莹, 张小平, 等. 我国西南地区风速变化及其影响因素[J]. 生态学报, 2014, 34(2):471−481. Zhang Z B, Yang Y, Zhang X P, et a1. Wind speed changes and its influencing factors in southwestern China[J]. Acta Ecologica Sinica, 2014, 34(2): 471−481. [34] 王秋华. 森林火灾燃烧过程中的火行为研究[D]. 北京: 中国林业科学研究院, 2010. Wang Q H. Study on fire behaviors in forest burning[D]. Beijing: Chinese Academy of Forestry, 2010. -
点击查看大图
图(10) / 表 (4)
计量
- 文章访问数: 194
- HTML全文浏览量: 73
- PDF下载量: 23
- 被引次数: 0
