第45卷第4期燕山大学学报Vol.45No.42021年7月Journal of Yanshan UniversityJuly 2021㊀㊀文章编号:1007-791X (2021)04-0352-05树冠火生长蔓延模型的改进与实时仿真吕梦雅1,2,闫㊀瑾1,2,任喜亮1,2,唐㊀勇1,2,∗,周莉莎1,2(1.燕山大学信息科学与工程学院,河北秦皇岛066004;2.河北省计算机虚拟技术与系统集成重点实验室,河北秦皇岛066004)㊀㊀收稿日期:2019-12-19㊀㊀责任编辑:孙峰基金项目:河北省自然科学基金资助项目(F2018203060)㊀㊀作者简介:吕梦雅(1965-),女,辽宁黑山人,博士,教授,主要研究方向为计算机图形学㊁多媒体技术;∗通信作者:唐勇(1964-),男,四川遂宁人,博士,教授,博士生导师,主要研究方向为计算机动画㊁计算机图形学㊁虚拟现实技术及其应用,Email:tangyong@㊂摘㊀要:现有林火研究多偏重于火焰蔓延的数值模拟以及简单可视化,存在火焰细节丢失和蔓延过程不精确的问题,本文通过深入分析地表火向树冠火的转换过程,完成树冠火生长蔓延模型的改进,并进行实时仿真㊂首先,为解决火焰蔓延过程中火焰运动缺乏真实感及细节丢失问题,引入脉动风场重新构建树冠火生长蔓延模型,丰富火焰细节,提高火焰运动真实性;其次,应用能量守恒定律,建立树冠火随温度场动态变化模型,提高地表火向树冠火生长转换过程精确性;最后,通过纹理映射技术模拟烧焦现象,采用Huygens 原理结合树冠火生长蔓延模型对火焰蔓延进行优化仿真,保证实时性㊂实验对比和分析表明,能够实时仿真出树冠火生长蔓延过程㊂关键词:树冠火生长蔓延;火焰模拟;脉动风场;温度场;Huygens 原理中图分类号:TP391.9㊀㊀文献标识码:A㊀㊀DOI :10.3969/j.issn.1007-791X.2021.04.0080 引言近年来,各地林火火灾高发,山东威海㊁四川凉山㊁山西沁源等陆续发生的森林火灾都带来重大损失,变化多端的林火蔓延尤其给火灾预防和扑救带来巨大的困难㊂其中树冠火蔓延对森林的摧毁性最大,不仅能烧毁针叶㊁树枝和地被物等,而且燃烧产热量大㊁蔓延速度飞快,难以扑灭㊂因此研究树冠火的发生条件和蔓延规律,并对树冠火的生长转换过程进行实时仿真,有着重要的学术价值和现实意义㊂1946年,W.R.Fons 首次提出林火蔓延数学模型,此后,许多学者也都开始利用数学工具对林火蔓延进行建模并以此来分析预测林火的蔓延过程㊂目前,用于刻画林火蔓延的数学模型主要有McArthur 模型㊁加拿大林火蔓延模型㊁Rothermel 模型和王正非林火蔓延模型,此外还有很多基于上述模型修正演变而来的各种蔓延模型㊂其中树冠火蔓延模型是广受关注的重点之一㊂此类模型多以地表火蔓延模型为研究基础,重点开展地表火向树冠火的生长转换和树冠火自身蔓延两部分的研究㊂1972年,美国农业部Rothermel [1]采用了加权平均法获得可燃物的参量,从宏观尺度来描述林火蔓延,该模型是一个半经验模型,没有考虑温度对火焰蔓延模型精确性的影响;1977年,加拿大森林研究所Wagner [2]提出了树冠火蔓延模型,但是他偏向于数值模拟和简单可视化,模拟过程复杂,同时没有考虑交互的影响,效果不够真实;1996年,中国科学技术大学吴清松等人[3]对燃烧过程中水分蒸发等物理条件对树冠火蔓延速率的影响进行了研究,并建立了物理模型从而控制树冠火蔓延,该模型是一维模型,适用范围小,同时没有对转换过程进行描述;2012年,天津大学刘世光等人[4]提出一种基于Rothermel 模型的地表火随时间变化的建模方法,并对火焰蔓延进行真实感绘制,但实现的火焰细节不够丰富;2013年,燕山大学唐勇等人[5]充分考虑了火焰燃烧过程中的物理特性,并利用基于物理的N-S 方程对火焰运. All Rights Reserved.第4期吕梦雅等㊀树冠火生长蔓延模型的改进与实时仿真353㊀动进行建模,并对火焰进行光线投影渲染后移植到GPU 进行计算,极大提高了仿真速度,但是加入烟雾和阴影之后仿真真实性有待提高;2015年,丁伟龙等人[7]为了在视觉上模拟单一植物的燃烧现象,提出了一种基于物理的多参数植物燃烧模拟方法,同时引入了分形和粒子系统来增加植物模型燃烧的真实性,但该方法没有考虑自然环境因素对火焰运动的影响,实现的火焰运动不够真实,仿真效果不好;2017年,斯坦福大学Sören Pirk等人[8]通过改变植物的物理化学属性,实现与火焰的交互,生成真实的植物燃烧现象,不过他们的研究局限于单株树木,未考虑大场景多株树木下树冠火生长蔓延运动的情况;2018年,Miguel GCruz 等人[11]分析了5种不同燃料类型的火势蔓延模型的准确性,并与其他模型进行对比,但是模拟精确性较低,误差较大;2019年,奥胡斯大学Michael B.Nielsen 等人[12]通过物理方程建模的方式模拟物体燃烧,视觉效果上很逼真,但缺乏实时性㊂因此,本文提出一种改进的树冠火生长蔓延模型,通过引入脉动风场丰富火焰细节,提高火焰蔓延真实性㊂再次,应用能量守恒定律,建立树冠火随温度场动态变化模型,增加树冠火转换过程的精确度㊂最后,运用Huygens 原理和纹理映射技术对火焰蔓延进一步优化仿真,保证实时性㊂1㊀树冠火生长蔓延模型1.1㊀传统的树冠火生长蔓延模型本文采用Wagner 提出的树冠火模型进行计算,此模型主要考虑了树木自身因素对火焰蔓延过程及火焰运动的影响,包括树叶的含水率M 和林木枝下高H ,计算公式如下:I O =[0.01H ㊃(460+25.9M )]32,(1)I b =I R60㊃12.6Rσ,(2)其中,I O 是树冠火发生临界状态时地表火蔓延的火焰强度阈值,I b 是地表火蔓延时火焰强度,I R 是地表火蔓延的反应强度,R 是地表火火焰蔓延速率值,σ是蔓延时可燃物的表面积与体积比㊂1.2㊀改进的树冠火生长蔓延模型传统树冠火蔓延模型没有考虑外力项因素对蔓延过程及火焰运动的影响,难以用于复杂条件下火焰蔓延的模拟㊂真实环境中,火焰蔓延会受到风力㊁地形等外部因素的影响㊂其中,风速作为主要影响因子,其对火焰蔓延的影响远超其他因素㊂因此,本文引入脉动风场重新计算火焰强度,对传统树冠火生长蔓延模型进行改进,以增加火焰的细节,提高火焰蔓延模拟的真实性㊂改进后计算公式为I O =[0.01H ㊃(460+25.9M )]32㊃14.2e0.1547v (t ),(3)式中,v (t )为脉动风速㊂火焰的蔓延会受到风场的影响,将脉动风场看作是一个典型的非完全均匀时空随机场㊂设m 个点空间相关脉动风速时程列向量的(Autoregressive model)模型可表示为v (t )=-ðpK =1ψK v (t -K Δt )+N (t ),(4)其中,p 是AR 模型的阶数值;Δt 是模拟风速时间步长;ψK 是AR 模型的自回归系数矩阵;v (t -K Δt )是t 时刻之前K 个时刻的脉动风速;N (t )是独立随机过程向量,N (t )=L ㊃n (t ),n (t )是正态随机过程;L 是m 阶下三角矩阵㊂将速度的表达式代入式(3)得到火焰强度I O值,通过引入脉动风场改进火焰强度的计算,增加火焰的细节,提高火焰运动真实性㊂2㊀树冠火随温度场动态变化过程在火焰蔓延过程中,当地表火蔓延火焰强度I b 大于阈值I O 时,就会引发地表火向树冠火的转换;反之,当I b 小于等于I O 时,表示地表火不会向树冠火发生转换㊂为了提高模拟的准确性,本文以任一时刻温度变化状态出发,应用能量守恒定律,建立树冠火随温度场动态变化生长过程㊂结合Wagner 提出的树冠火生长蔓延模型,得到火的蔓延速度计算式:R =R O +αT ,(5)其中,R O 是火的初始蔓延速度,T 是温度,α是比例系数,值为0.05㊂由能量守恒定律推导出温度场随时间变化的微分方程表达式:. All Rights Reserved.354㊀燕山大学学报2021ρC P ∂T ∂t =λ∂2T ∂x 2+∂2T ∂y 2+∂2T ∂z 2()+∂L∂t ,(6)其中,T 为火场温度;t 为火焰蔓延时间;ρ为火源密度;C P 为比热;λ是导热系数;L 是潜热㊂为简化计算,忽略潜热的计算项,即温度场的表达式为ρC P∂T∂t =λ∂2T ∂x 2+∂2T ∂y 2+∂2T ∂z2()㊂(7)㊀㊀使用Jacobi 迭代求解出每一时刻T 值,将其代入式(5)得到火焰蔓延的速度R ㊂通过引入温度场动态变化函数改进火焰蔓延速度的计算,使树冠火生长过程更加准确㊂3㊀树冠火生长蔓延的实时仿真3.1㊀基于Huygens 原理的实时仿真为了使树冠火生长蔓延更具有真实感和实时性,本文采用基于Huygens 的波动传播模型对火焰蔓延区域进行实时仿真,原理是通过计算火场边界上每一个点的过火区形状来描述火场边界,过火区由随时间动态变化的连续扩展多边形表示,为保证实时性和精度要求,将选择多边形的一点设置为一个独立的初始着火点,本文将多边形顶点的数量控制在1~624,在计算中,一定时间间隔内参数值设置是近似不变的,火场边界是通过这些着火点依次引燃邻近未燃区来完成蔓延的㊂根据Huygens 波动传播原理进行蔓延区域的实时仿真可分三步进行:首先,选取一个着火点,蔓延形成一系列着火点㊂然后,应用连续扩展的多边形顶点判别的算法,计算逐个着火点的位置关系,定义3点P 1(x 1,y 1),P 2(x 2,y 2),P (x 3,y 3)的坐标,其行列式形式如下:det(P 1,P 2,P 3)=x 1y 11x 2y 21x 3y 31㊂(8)㊀㊀最后,设图1中各个连续的顶点为P 1,P 2, ,P i ,利用行列式结果的正负值判断着火点向量的位置关系㊂重复上述过程可以建立多个着火区域,通过实验可证明基于Huygens 原理的仿真可以优化火焰蔓延过程,保证实时性和计算精确性,使蔓延过程更加真实㊂图1㊀多边形顶点示意图Fig.1㊀Polygon vertices3.2㊀ 烧焦 纹理仿真基于物理的烧焦效果真实但实时性差,本文采用设定燃料类型的方法,通过模拟过火前后地表形态的变化,获得较为真实且实时的烧焦效果㊂被燃烧纹理是主要燃料,对于每个被燃烧纹理,在操作面板中设定对应的燃料值,如果纹理不易燃,则唯一需要设定的变量是 纹理ID ,它是纹理顺序中纹理的数组索引:0~n ,本文设定不同的烧焦标记纹理,将表示随时间变化不同的烧焦程度㊂由实验结果可知,使用纹理映射技术可以模拟不同的烧焦程度,并且在树冠火生长蔓延实时仿真中依然可以保持很高的帧率㊂4㊀实验结果与分析本文实验基于Windows 系统,处理器为IntelCore i594002.90GHz,内存为16G,显卡为NVDIAGeForce GTX 1650㊂软件环境采用Unity3D 虚拟引擎进行系统开发,并结合C #以及CG /Unity Surface Shader 编程语言㊂图2为引入动态温度场和脉动风场的火焰模拟,a)为文献[6]通过传统的粒子系统模拟的火焰效果,可以看出火焰的边缘较为模糊,缺乏真实感,b)为文献[11]通过物理方法模拟的火焰效果,火焰的细节比较真实,但是做不到实时,c)为本文引入动态温度场和脉动风场进行改进后的火焰模拟效果,火焰细节更为明显,更具真实感㊂图3为本文改进的树冠火蔓延模型模拟效果,a)为文献[7]基于生物学-数学公式和支持建模的方式联合模拟植物树模型燃烧效果,b)为文献[6]基于物理定律的方法模拟植物燃烧过程的. All Rights Reserved.第4期吕梦雅等㊀树冠火生长蔓延模型的改进与实时仿真355㊀效果,可以看出二者火焰运动不够真实,细节不够丰富,c)为本文引入脉动风场后树木燃烧的蔓延模拟,火焰运动更具有随机性,效果更真实㊂图2㊀火焰模拟细节对比Fig.2㊀Comparison of flame simulation details图3㊀树冠火燃烧模拟对比Fig.3㊀Comparison of crown fire combustion simulation㊀㊀图4为本文通过调整风系数,不同风向不同风速对火焰的蔓延方向以及蔓延大小产生不同的影响,根据大气运动原理可知,风的方向和火焰蔓延方向一致㊂a 为无风情况下火焰的蔓延效果,b 引入风场后的火焰蔓延效果且风为东风㊁风速为2,c 为风速为5下火焰的蔓延情况㊂图4㊀引入风场对比Fig.4㊀Comparision of introducing wind field㊀㊀图5为引入动态温度场树冠火的生长转换过程㊂a 为I b ɤI O 时效果图,此时地表火不会向树冠火发生转换,只是地表火蔓延,b 为I b >I O 时效果图,此时地表火强度大于阈值,发生地表火向树冠火的转换㊂通过应用能量守恒定律得到随时间动态变化的火焰蔓延速度,该方法使转换过程更加精确,同时满足火焰模拟的实时性要求,大大加强了真实感㊂图6为采用Huygens 原理和纹理映射技术对烧焦效果进行实时仿真㊂a 为文献[4]的林火蔓延烧焦效果图,火焰的细节和效果都不是很好,同时做不到实时㊂b,c,d 为本文烧焦效果图,通过不同的纹理ID,展示不同时间的烧焦程度,可以看到烧焦效果更好,更具有真实感,同时保证实时性㊂图5㊀地表火向树冠火转换Fig.5㊀Conversion of surface fire to crown fire图6㊀不同烧焦程度对比Fig.6㊀Comparison of different scorch levels㊀㊀为检测改进树冠火生长蔓延模型的仿真效率,表1列出了本文实验以及部分文献数据的对比状况㊂其中,可以明显看出本文模型增加了火焰模拟细节的同时,确保了树冠火生长转换的精确度,保证了实时性㊂表1㊀不同方法的实验性能对比Tab.1㊀Comparison of experimentalperformance of different methods实验场景方法粒子数帧率/fps 火焰细节文献[11]3000023.7模拟实验本文方法3000065.2树冠火燃烧文献[6]3000014.6模拟实验本文方法3000072㊀不同烧焦程文献[4]3800023.6度模拟实验本文方法3800069.25 结论本文完成了对树冠火生长蔓延模型的改进㊂首先,引入脉动风场改进树冠火生长蔓延模型,增强了火焰细节,使火焰模拟更加真实;然后,应用. All Rights Reserved.356㊀燕山大学学报2021能量守恒定律建立动态温度场模型,使树冠火转换过程更加精确;最后,将Huygens原理与树冠火生长蔓延模型相结合,并采用纹理映射技术的方法改善了实时仿真的问题㊂但是该模型的适用范围还存在一定问题,没有考虑不同可燃物以及更加多变的环境因子下火焰的蔓延情况,场景过大时系统的运算速率偏低,这些将是我们未来工作的重点㊂参考文献1ROTHERMEL R C.A mathematical model for predicting fire spread in wildland fuels INT-115 R .Ogden UT U.S. Department of Agriculture Intermountain Forest and Range Experiment Station 1972.2WAGNER C E.Conditions for the start and spread of crown fire J .Canadian Journal of Forest Research 1977 7 1 23-34. 3吴清松童盛周建军.树冠火蔓延模型和数值分析 J .火灾科学1996 5 2 29-34.WU Q S TONG S ZHOU J J.Crown fire propagation model and numerical analysis J .Fire Safety Science 19965 2 29-34. 4刘世光傅应坻蒋浴芹.地表火蔓延的真实感仿真 J .自然灾害学报2012 21 2 180-186.LIU S G FU Y D JIANG Y Q Realistic simulation of surface fire spread J .Journal of Natural Disasters 2012 21 2 180-186. 5唐勇刘一君吕梦雅.NS方程和体渲染结合的火焰实时模拟 J .小型微型计算机系统2013 34 12 2845-2848. TANG Y LIU Y J LÜM Y.Real-time fire simulation combined with NS equation and volume rendering J .Journal of ChineseComputer Systems 2013 34 12 2845-2848.6唐勇张利卫吕梦雅等.自由场景下烟雾动态追踪的实时仿真 J .燕山大学学报2015 39 6 541-545. TANG Y ZHANG L W LÜM Y.Real-time simulation of smoke dynamic tracking in free scene J .Journal of Yanshan University 2015 39 6 541-545.7DING W L ZHENG L XU L F.Visual simulation of plant interacting with burning flame based on the law of physics C// IEEE International Conference on Multimedia Big Data Beijing 2015 362-365.8PINK S MICHAL J TORSTEN H et al.Interactive wood combustion for botanical tree models J .ACM Transactions on Graphics 2017 36 6 1-12.9张菲菲解新路.一种改进的林火蔓延模型及其实现 J .测绘与空间地理信息2012 35 2 50-53. ZHANG F F XIE X L.An improved forest fire spread model and its realization J .Geomatics&Spatial Information Technology 2012 35 2 50-53.10唐玉丽毛行辉陈崇成.基于FARSITE的林火蔓延三维可视化模拟 J .自然灾害学报2015 24 2 221-227. TANG Y L MAO X H CHEN C C.Three-dimensional visual simulation of forest fire spread based on FARSITE J .Journal of Natural Disasters 2015 24 2 221-22711CRUZ M G ALEXANDER M E SULLIVAN A L et al. Assessing improvements in models used to operationally predict wildland fire rate of spread J .Environmental Modelling& Software 2018 105 54-63.12NIELSEN M B STAMATELOS K BOJSEN H M et al. Physics-based combustion simulation in Bifrost C//ACM SIGGRAPH2019Talks.LosAngeles 2019 1-2.Improvement and real-time simulation of crown fire growth spread modelLÜMengya1 2YAN Jin1 2REN Xiliang1 2TANG Yong1 2ZHOU Lisha1 21.School of Information Science and Engineering Yanshan University Qinhuangdao Hebei066004 China2.Hebei Key Laboratory of Computer Virtual Technology and System Integration Qinhuangdao Hebei066004 ChinaAbstract Existing forest fire research mostly focuses on the numerical simulation and simple visualization of flame spreading andthere are problems with the loss of flame details and the inaccuracy of the spreading process.In this paper by analysis of thetransition from surface fire to crown fire in-depth the growth and spread model of crown fire was improved and real-timesimulation is performed.Firstly in order to solve the problem of lack of realism and loss of details during the flame spread processa pulsating wind field was introduced to reconstruct the growth and spread model of the crown fire to enrich the flame details andimprove the authenticity of the flame movement.Secondly the law of conservation of energy was used to establish a model of thedynamic change of the crown fire with temperature which improved the accuracy of the process of conversion from surface fire tocrown fire growth.Finally scorch phenomenon was simulated by texture mapping technology and Huygens principle was combinedwith crown fire growth and propagation model to optimize the simulation of flame spread to ensure real-time performance.Theexperimental results show that the crown fire growth and spread process can be simulated in real time.Keywords crown fire spread flame simulation pulsating wind field temperature field Huygens principle. All Rights Reserved.。
