设计火灾中火灾热释放速率曲

1、20/20设计火灾中火灾热释放速率曲线的确定基金项目基金项目：国家“十五”科技攻关项目“建筑物性能化防火设计技术导则的研究”(2001BA803B02-04)王志刚 倪照鹏 王宗存 姜明理（公安部天津消防研究所 天津 300381）摘要：火灾的热释放速率是火灾进展过程中十分重要的参数，在运用火灾模型进行建筑物性能化消防设计与评估时，确定正确合理的火灾热释放速率曲线至关重要。本文总结了几种确定火灾热释放速率曲线和计算火灾热释放速率的方法，列举了相关计算实例。关键词：性能化消防设计；设计火灾；热释放速率1 前言在进行性能化消防设计与评估时，通常需要运用火灾模型对建筑内发生火灾时的可能场景进行模拟计

2、算，得到建筑内各区域有关火灾参数的变化进展情况，如温度、燃烧产物浓度、火焰高度、烟气层界面高度等，然后依照这些计算结果，推断不同条件下的火灾是否会对建筑内人员生命安全或结构稳定等造成危害，评价相关消防设计能否将火灾危险降低到一个合理的水平以及确定是否需要进一步完善修改或调整相关措施。火灾进展过程中的各种火灾参数直接取决于火灾的热释放速率。但当前绝大多数火灾模型都没有包含燃料的燃烧模型，其要紧缘故是燃烧是一个相当复杂的物理化学过程，对其进行数值模拟是一项极其艰巨的任务，在目前条件下还无法完全实现。鉴于这种状况，各种火灾模型都把热释放速率作为一个输入参数由用户来输入，因此输入的火灾热释放速率曲线能

3、否代表真实的火灾情况将决定火灾模拟结果的精确程度。本文依照国家“十五”科技攻关项目“建筑物性能化消防设计导则的研究”所进行的研究工作并结合一些工程实例的经验，总结了几种确定火灾热释放速率曲线的方法，并以某高架库房的设计火灾分析作了一实际算例。2 火灾的进展过程火灾的整个进展过程按时刻先后顺序分为阴燃、增长时期、充分进展时期、衰退时期直到最终熄灭，如图1所示。视燃料的燃烧特性，阴燃时期的持续时刻有所不同，关于液体和气体火灾甚至没有阴燃时期。在火灾增长时期，随着持续燃烧时刻的延长，如无外界条件干涉，会有越来越多的可燃物参与燃烧，火灾的热释放速率也相应不断增大。当火灾进展到一定程度，（如着火室空间条

4、件合适，还会出现轰燃现象）火灾进入其充分进展时期，随后的一段时刻内的火灾热释放速率差不多保持稳定。其最大热释放速率要紧取决于燃料的数量与性质、着火空间的通风条件等。火灾通过一段时刻的持续进展后，随着燃料数量的减少，火灾最终进入衰退时期，当燃料全部耗尽后火灾随之熄灭。在消防安全工程学中，通常选取火灾危害可能较小但发生火灾的可能性专门大的火灾以及发生火灾的可能性较小但危害较大的多个火灾场景组成一组，来对现有的消防设计进行评估。关于火灾的各个时期，也是如此。由于阴燃时期产生的热量较少，进展缓慢，在对常见建筑物的消防安全进行评估时通常忽略火灾的阴燃时期。火灾在衰退时期时的热释放速率呈下降趋势，而且通常

5、情况下，现在火灾差不多过了相当长的时刻。为简化处理，关于衰退时期也不予考虑。因此在实际中，人们最关怀的确实是火灾的增长时期和充分进展时期，在构建火灾的热释放速率曲线时，也要紧是考虑这两个时期。阴燃增长时期充分进展时期阴燃增长时期充分进展时期衰退时期火灾热释放速率时刻3 热释放速率曲线的确定燃料的燃烧是一个相当复杂的物理化学过程，对其进行数值模拟是一项极其艰巨的任务。当前绝大多数火灾模型都没有包含燃料的燃烧模型，在运用火灾模型进行火灾危害分析时，通常都直接采纳实验测量得到的热释放速率曲线，或采纳一些简化的特征火灾进展模型，如T平方进展火灾。以下是几种确定火源热释放速率曲线的方法：3.1 实际火灾

6、实验通过实际的火灾实验，能够得到火灾的热释放速率曲线。依照实验的规模，火灾实验可分为实验室规模实验、中型实验以及实尺寸火灾实验。其中实验室规模的实验要紧通过锥形量热器测量单一可燃物的热释放速率，试样的大小约为100mm100mm。中型实验除了能够测量单一可燃物的热释放速率外，还能够测量几种可燃物组成的可燃组件的热释放速率，试样的大小最大可达1m1m。大型实尺寸火灾实验模拟建筑的实际尺寸，依照实际的可燃物种类、火灾荷载及摆放方式进行实验，事实上验结果与真实火灾较为接近。但由于这类实验的花费较大，目前此类实验的相关数据较少。相对较多的是单一可燃组件的火灾实验数据。三种实验方法中，最有有用价值的是实

7、尺寸火灾实验，其次是中型火灾实验，而实验室规模实验结果则常用于分析火灾引燃源。3.2 类似实验由于已做的火灾实验有限，多数情况下找不到有待分析的可燃组件的实验数据。在这种情况下，能够找其他具有类似的燃料类型、燃料布置及引燃场景的火灾实验数据。因此，所采数据的实验条件与实际要考虑的情况越接近越好。例如：在考虑会展中心中的一个展位发生火灾时，因缺少展位起火的实验数据，能够采纳一个办公家具组合单元的火灾试验数据1，实验的照片见图2。实验中的办公家具组合单元包括两面办公单元的分隔板、组合书架、软垫塑料椅、高密度层压板办公桌以及一台电脑，还有98kg纸张和记事本等纸制品。该办公家具组合单元中包含了展览中

8、较为常见的可燃物，物品的摆放形式也差不多与展位的布置相同，且其尺寸与一个展位相当。因此，在缺少展位火灾实验数据的情况下，能够用如此一个办公家具组合单元的火灾试验数据来替代，相关实验获得的火灾热释放速率曲线见图3。图2 一个办公家具组合单元火灾实验的燃烧情形图3 办公家具组合单元火灾实验热释放速率实测曲线3.3 依照火灾蔓延进行计算火灾在初期仅有少量的可燃物参与燃烧，之后通过热辐射将相邻区域内的可燃物引燃，再通过一段时刻后，被引燃的可燃物产生的热辐射又将邻近的可燃物引燃。假如火源周围存在足够的可燃物，随着时刻的推移，不难设想卷入火灾的可燃物将会成倍增加，火势也随之不断增长。如图4所示，在距火源中

9、心半径为R的空间范围内，火源对该区域内可燃物的热辐射为2：(1)其中：q对目标可燃物的单位辐射热流，kW/m2；P火焰的总辐射热流，kW；R与目标可燃物之间的距离，m；Xr热辐射效率，依照不同的燃料类型在0.20.6范围内取值；Q火源的总热释放速率，kW。关于一般的可燃物，公式(1)中的热辐射效率Xr取1/3，即火源1/3的能量以热辐射的方式传播出去。因此，公式(1)可转化为公式(2)：.(2)火源半球目标可燃物火源半球目标可燃物假如明白了可燃物被引燃的热辐射q以及其距火源的距离R，就能够通过公式(2)计算出可燃物被引燃时火源的热释放速率。若火源的热释放速率曲线已知，就能够确定火源能够达到引燃

10、可燃物程度的时刻，也就相应能确定可燃物被引燃的时刻。可燃物被引燃后，作为次生火源与起始火源的功率叠加在一起，又会向周围的可燃物发射更强的辐射热流。运用同样的方法计算下一个可燃物被引燃的时刻，如此进行下去，就能够得到在一定时刻内的火灾热释放速率的变化情况。关于特定的可燃物，它被引燃的热辐射通量能够通过实验手段测定，也能够查相关的实验数据获得。但在工程计算中通常依照被引燃的难易程度将可燃物分为三类3，参见表1：表1 可燃物被引燃难易程度的粗略分类可燃物类不单位面积可燃物表面引燃所需要的辐射热流(kW/m2)易引燃10一般可引燃20难引燃40上述计算方法适用于物品摆放形式规则一致，可燃物种类单一的情

11、况。较为典型的例子有图书馆内的书库、高架仓库内的货架以及展览中心内的展位等。算例4：分析一个物流仓库内的火灾蔓延情况及确定火灾的热释放速率曲线。图5为某物流高架库房内货架商品的摆放情况示意图。图中两个货架背靠背布置形成一个货架列，物资储存在标准规格的自动化存取包装箱内，放置在货架内的货架托盘上。本文将通过分析货架列内的一个货架单位着火时，火灾沿货架列的传播蔓延情况，以讲明如何应用上述计算方法确定火灾的热释放速率曲线。货架物资货架托盘货架物资货架托盘图5 货架商品摆放示意图（1）分析火灾蔓延的差不多情形为便于分析，将货架的正面视图用图6抽象表示。由图6可见，货架共有6层，其中每个格子对应上图中的

12、一个货架格（以下简称为货架单位）。由于货架是背靠背布置，因此每个货架单位中放置两个货架托盘和两个装满物资的箱子。假定起火点位于货架底层位置为A处的货架单位，起火之后，火灾会沿垂直和水平方向分不向上部及其两侧蔓延。按货架单位被引燃的先后顺序将货架单位分组，火焰蔓延的先后顺序则可如图6所示表现为A-B-C-D-E-F-G-H。由此能够看出，火灾向上蔓延首先呈金字塔状向上蔓延，同时火灾还向货架两侧蔓延，使得火灾蔓延过程将逐步呈梯形状进展。HGFGHHGFEFGHHGFEDEFGHHGFEDCDEFGHHGFEDCBCDEFGHHGFEDCBABCDEFGH图6 货架单位A起火后，火灾沿货架的传播蔓延

13、情形简化示意（2）确定首先起火的一个货架单位A的热释放速率进展情况由于货箱内的物资不确定，本文采纳木制货架托盘的燃烧实验结果作为分析的基础。该实验中参与燃烧的共有六个托盘（实验照片见图7），将其火灾载荷等价于一个货架单位内托盘及其上面物资的火灾载荷。图8为实验得到的热释放速率曲线，并在计算中将其视为一个货架单位着火时的热释放速率曲线。货架托盘的尺寸为1m1m，两个托盘之间的间隙约为30mm，因此一个货架单位起火后，其火源中心距相邻货架单位的距离为0.5+0.03=0.53（m）。计算时，货架托盘及商品按一般可引燃物品计算(见表1，其引燃所需要的辐射热流为20kW/m2)。依照公式(2)，得到引

14、燃相邻货架单位所需的热释放速率为：（kW）图7 木制货架托盘燃烧实验情形图8 一个货架单位着火时的热释放速率曲线从一个货架单位的热释放速率曲线（图8）能够查出，火源热释放速率约在180s（3min.）时达到211kW，即表示3min后相邻的货架单位将被引燃。相邻的货架单位被引燃后，也按相同的热释放速率曲线进展。依照相同的方法再计算其他相邻的货架单位被引燃的时刻。图9为按被引燃的先后顺序将一个货架列中不同货架单位组的火灾热释放率曲线绘制而成。一个货架列的火灾热释放速率一般以时刻为步长，综合各货架单位顺序起火后某一时刻的热相互作用来获得，但目前这方面的研究尚不完善。因此，本文假设一个货架列的火灾热

15、释放速率差不多等同于将各组热释放速率曲线的合成结果。如此，就能够得到一个货架列中货架单位A-H组起火时的火灾热释放速率随时刻的变化曲线（如图10所示）。图9 货架底层起火时，同列其他货架单位起火后的热释放速率曲线图10 一个货架列中货架单位A-H组起火时的合成热释放速率曲线3.4 特征火灾曲线当确实没有其他相关信息来得到火灾的热释放速率曲线时，还能够进行适当的假设，将火灾的进展特征化。一种假设是认为火灾将按稳态火灾进展，即火源的热释放速率始终保持一个恒定的值。另一种假设是认为火灾按非稳态火灾进展，但火灾的进展具有一定的规律性，例如T平方特征火灾。本文要紧介绍后一种类型的火灾非稳态火灾。3.4.

16、1 稳态火灾关于稳态火灾，在其整个进展过程中，火源的热释放速率始终保持一个恒定的值。火灾进展过程中的充分进展时期能够近似看成是稳态火灾。通常在进行排烟系统的计算时一般保守地设定火灾为稳态火灾，从而计算排烟系统的排烟速率和排烟量。稳态火灾的火灾热释放速率可采纳公式(3)计算：.(3)其中：Q稳态火灾的热释放速率，kW；燃料的质量燃烧速率，kg/s；hc燃料的燃烧值，kJ/kg。3.4.2 T平方特征火灾T平方特征火灾是非稳态火灾之一，消防安全工程学中一般采纳该模型来描述非稳态火灾热释放速率随时刻的变化。该模型假定火灾热释放速率与时刻的平方成正比，见公式(4)：(4)其中：火灾增长系数，kW/s2

17、；t火灾进展时刻，s。超快火快速火中速火慢速火热释放速率(KW)超快火快速火中速火慢速火热释放速率(KW)图11 各类T方特征火灾的热释放速率曲线表2 对应于T方火灾的各类可燃物火灾增长分级典型的可燃材料火灾增长系数(kW/s2)热释放速率达到1MW的时刻3,10（s）慢速火0.002931（0.00277）3600中速火棉质/聚酯垫子0.01127（0.0111）300快速火装满的邮件袋、木制货架托盘、泡沫塑料0.04689（0.0444）150超快火池火、快速燃烧的装饰家具、轻质窗帘0.1878（0.177）75从图11能够看出，T平方特征火灾实际描述的是火灾的增长时期，那么如何确定火灾何

18、时达到它的充分进展时期呢？或者讲火灾在充分进展时期的热释放速率能达到多大规模？本文按以下两种情况分不进行分析：（1）火灾为无限制地自由进展状态在这种状态下，火灾无限制地自由进展到其最大热释放速率(见图12中曲线)，这种状态也能够认为是自动灭火系统失效的情况。火灾的最大热释放速率可由公式(5)计算6：.(5)其中：Qmax最大热释放速率，kW；q单位面积的热释放速率，又称热释放速率密度，kW/m2；Ac可燃物所占的面积，m2。（2）火灾受灭火系统的操纵灭火系统动作并有效作用的情况下，火灾将受到抑制，热释放速率不再接着增长，继而会呈下降趋势。但在消防安全工程学中，为安全起见，通常保守地假定灭火系统

19、动作后热释放速率不再增长，同时接着维持在这一水平(见图12中曲线)。火灾的最大热释放速率能够按自动喷水灭火系统动作时的热释放速率考虑。自动喷水灭火系统的动作时刻可运用美国国家标准及技术研究院（NIST）开发的软件工具包FPETool中的DETECT-QS工具计算7。热释放速率热释放速率时刻QcQ时刻QcQmax图12 火灾分不在自由进展和受到灭火系统操纵两种状态下的热释放速率曲线关于火灾的热释放速率计算方法，许多学者尝试从各个角度进行了多方面研究，得出了一些经验公式8,9，限于篇幅不再赘述。4 结束语在运用火灾模拟模型进行性能化消防设计与评估时，需要输入火灾的热释放速率曲线。该曲线能否代表火灾

20、的真实情况直接阻碍性能化消防设计与评估的可靠性及其应用。热释放速率曲线可直接通过火灾实验获得，但由于实尺寸火灾实验的费用较大，此类可用的实验数据较少，而较多的是中型火灾实验与实验室规模的火灾实验数据。当无法找到有待考虑的可燃组件的实验数据时，能够采纳类似的火灾实验数据替代，因此实验采纳的燃料类型、布置及引燃场景与实际情况越接近越好。本文还介绍了另外一种确定热释放速率曲线的方法，即在确定了火灾引燃源之后，通过计算火灾的蔓延情况来得到火灾的热释放速率曲线，并举了一个计算实例。在没有相关数据确定火灾的热释放速率曲线时，可近似采纳特征火灾曲线。火灾的热释放速率是进展火灾过程中的最为重要的参数，在运用火

21、灾模型进行计算时，确定正确合理的火灾热释放速率曲线至关重要。本文仅对一些常用的有用计算方法进行了总结。在今后的研究工作中，一方面需要广泛进行各类场所内各种典型火灾场景的火灾实验，丰富补充目前的火灾实验数据库，从而为推广与完善性能化消防设计提供完备的数据基础。另一方面也要进行燃料燃烧过程数值模拟的研究，并尽早将其纳入火灾模型中，尽管迄今为止这仍是火灾科学界的一大难题。参考文献1 NIST网站: .2 Edward K.Budnick, Harold E.Nelson, David D.Evans., “Simplified Fire Growth Calculations”, Fire Protection Handbook, 18th Edition, National Fire Protection Association, Batterymarch Park, Quincy, MA, 1997.3 “Fire Engineering Guidelines” First Edition-March 1996, Fire Code Reform Center