廖超 陳愛斌 王丹

【摘要】我國的年均森林火災(zāi)次數(shù)和受害面積都很嚴重，森林火災(zāi)突發(fā)性強、破壞性大、處置救助較為困難，嚴重地制約著現(xiàn)代林業(yè)的可持續(xù)發(fā)展。森林火災(zāi)的預(yù)測預(yù)報是森林防火工作中最重要的一環(huán)，本文在分析林火行為的基礎(chǔ)上，闡述了國內(nèi)外幾種經(jīng)典的林火蔓延模型及其優(yōu)缺點和適用范圍，分析了國內(nèi)外林火蔓延模型的發(fā)展狀況，并著重介紹了目前最新提出的基于離散事件驅(qū)動的林火蔓延模型，最后對未來林火蔓延發(fā)展的趨勢進行了展望。

【關(guān)鍵詞】森林火災(zāi);林火蔓延模型;離散事件

1.引言

森林火災(zāi)，是指失去人為控制，在林地內(nèi)自由蔓延和擴展，對森林、森林生態(tài)系統(tǒng)和人類帶來一定危害和損失的林火行為，其突發(fā)性強、破壞性大、處置救助較為困難，森林火災(zāi)嚴重地制約著現(xiàn)代林業(yè)的可持續(xù)發(fā)展。森林火災(zāi)會產(chǎn)生大量煙霧，污染環(huán)境;大量二氧化碳與水發(fā)生化學(xué)反應(yīng)，嚴重污染水環(huán)境，嚴重破壞了全球的生態(tài)平衡，使已經(jīng)很脆弱的地球生態(tài)環(huán)境進一步惡化。據(jù)聯(lián)合國糧農(nóng)組織報告，近年來森林火災(zāi)發(fā)生的頻率和規(guī)模十分反常，已遠超過自然界的自我調(diào)節(jié)范圍，對人類、環(huán)境和經(jīng)濟造成了前所未有的危害。

森林火災(zāi)的預(yù)測預(yù)報是森林防火工作中最重要的一環(huán)，因此加強森林火災(zāi)的預(yù)測預(yù)報成為森林防火工作必須面對的一個重要問題，對各級政府做好森林消防決策，最大限度地減少森林火災(zāi)造成的損失，具有極其重要的作用。

林火蔓延是一種林火行為，與火災(zāi)預(yù)防、撲救等工作密切相關(guān)。林火行為是指森林可燃物在點燃后，所產(chǎn)生的火焰和火蔓延以及發(fā)展過程的特征，亦即是林火發(fā)生、發(fā)展，直至熄滅的全過程中著火、蔓延、能量釋放、火強度、火災(zāi)種類等特征的綜合。林火蔓延模型是林火行為研究的一個主要方面，受可燃物、火環(huán)境（氣象、地形和植被）和火源條件的制約和控制，量化的林火蔓延模型多是在各種簡化條件下，運用數(shù)學(xué)方法進行處理，導(dǎo)出蔓延速度與各種參數(shù)，諸如可燃物的理化性質(zhì)、地形、氣象因子等之間的定量關(guān)系式，供人們利用這些關(guān)系式去預(yù)測要發(fā)生的林火行為，從而指導(dǎo)撲火工作，進行日常的林火管理等。

根據(jù)模型建立的方法以及對林火蔓延本質(zhì)的認識程度，模型可分為統(tǒng)計模型、半物理模型、物理模型三類[6]。

統(tǒng)計模型不涉及任何物理機制，純粹從統(tǒng)計的角度來描述火行為，該方法把有多個變量相互關(guān)系的復(fù)雜問題，在形式上作簡單的處理，因建立在大量實際森林火災(zāi)和計劃火燒的資料基礎(chǔ)上，資料充足，有可靠的置信度，故公式計算結(jié)果與實際情況基本符合。

物理模型最先是由Fons提出來的，他的模型將燃料床理想化了，并且認為燃料達到著火的溫度即著火，使得模型與實際相距較大，且模型表達式復(fù)雜，物理參數(shù)較多又難以確定。這種模型的局限性在于它跟實際火災(zāi)情況有較大的差距，但把林火蔓延抽象成一個純物理問題來研究，為人類提供了一條認識林火蔓延一般規(guī)律的途徑。

2.幾種經(jīng)典的林火蔓延模型

目前比較經(jīng)典的林火蔓延模型有：美國的Rothermel模型[5]、澳大利亞的McArthur模型[1]、加拿大森林火險等級系統(tǒng)[4]和中國的王正非模型[2][3]。

2.1 Rothermel模型

Rothermel模型遵循能量守恒定律，屬于半物理模型。主要研究的是火焰前鋒的蔓延，沒有考慮過火火場的持續(xù)燃燒。提出了“似穩(wěn)態(tài)”的概念，從宏觀方面研究林火行為。它假定燃料床和地形地勢在空間上是連續(xù)分布的，并且可燃物的含水量、風(fēng)速、坡度等參數(shù)是不變的。在火災(zāi)蔓延過程中，考慮到了熱傳導(dǎo)、熱對流和熱輻射的物理機制。它建立在均一的可燃物狀態(tài)下，并且不現(xiàn)實的要求90%可燃物的直徑在1/8英寸，可燃物的含水量不得超過35%，對于比較大型的可燃物的影響通常忽略。

2.2 McArthur模型

McArthur模型是Noble I.R.等人對McArthur火險尺的數(shù)學(xué)描述。它是建立在多次點燃實驗上，導(dǎo)出的林火蔓延速度與各參數(shù)之間定量的關(guān)系式，屬于統(tǒng)計模型。它的優(yōu)點是能預(yù)報火險天氣和一些重要的火行為參數(shù)。但是，它的適用范圍極其有限，只適合草地和桉樹林。我國南方森林預(yù)警系統(tǒng)可選擇性的參考應(yīng)用。

2.3 CFFCRS林火蔓延模型

加拿大森林火險等級系統(tǒng)（CFFCRS）是以加拿大林火行為為理論基礎(chǔ)進行研發(fā)的。它包含四個子系統(tǒng)：林火行為預(yù)報系統(tǒng)（FBP），林火天氣指標系統(tǒng)（FWI），林火發(fā)生預(yù)報系統(tǒng)（FOP），可燃物濕度估測系統(tǒng)（AFM）。它是通過幾百次的實驗與觀察，導(dǎo)出的林火向前蔓延速度（ROS）方程，屬于統(tǒng)計模型?？扇嘉镱愋筒煌?，其蔓延速度方程也不同，各方程都獨立地依賴于初始蔓延指標ISI。該模型不考慮任何熱傳機制。不過，它可以幫助人們比較簡單的認識林火行為的各個分過程和整個過程，以及揭示其作用規(guī)律。在各參數(shù)相似的條件下，可以比較準確地預(yù)測火行為。

2.4 王正非林火蔓延模型

王正非模型也是建立在大量的實驗基礎(chǔ)上，得出的林火蔓延速度關(guān)系式，屬于統(tǒng)計模型。此模型只適用于平地?zé)o風(fēng)、無風(fēng)上坡和順風(fēng)上坡且坡度不能超過650的情況下。對于一些粗大的可燃物，如元木和大的枝干，視為無效可燃物。模型中的各個修正系數(shù)王正非都給出了相應(yīng)的關(guān)系表，可直接由表查找，簡化了計算。

3.林火蔓延研究新進展

3.1 國內(nèi)外最新研究進展

相比與澳大利亞、加拿大和王正非模型，Rothermel模型的適用范圍最廣，近年來，有許多學(xué)者在該模型的基礎(chǔ)上，提出了許多Rothermel修正模型。

在國內(nèi)，景文峰[7]以Rothermel模型為藍本，結(jié)合交互窗口、字符型漢化技術(shù)及專家系統(tǒng)的原理，提出了用于估算特定可燃物分布狀況、變環(huán)境條件下火蔓延趨勢的專家系統(tǒng)。王海暉[8]參考Rothermel模型，將火場假定為橢圓面，將風(fēng)速和坡度進行矢量疊加并引入折合風(fēng)速的概念，通過大量室內(nèi)林火實驗，提出了一整套估算地表火行為特征量的關(guān)系式。朱啟疆[9]借助地理信息系統(tǒng)（GIS）和遙感技術(shù)，將Rothermel半經(jīng)驗?zāi)Ｐ秃徒?jīng)驗?zāi)Ｐ拖嘟Y(jié)合，分別考慮帶方向（坡度和風(fēng)）和不帶方向因素對火勢蔓延的影響，并利用迷宮算法實現(xiàn)了林火空間蔓延動態(tài)模擬。李建微[10]以Rothermel模型為基礎(chǔ)，假定點火源的火場形狀為橢圓形，利用惠更斯原理，結(jié)合粒子系統(tǒng)方法、計算機圖形和虛擬現(xiàn)實技術(shù)，在不同風(fēng)速、坡度下，對不同位置的林火火場蔓延行為進行三維模擬。蔡學(xué)理、張貴[11]以AroInfo為地理信息系統(tǒng)平臺，結(jié)合Rothermel模型，并將火場形狀假設(shè)為多邊形，提出了一種森林火災(zāi)指揮撲救方案。

在國外，以Rothermel模型為基礎(chǔ)，所提出的應(yīng)用比較廣泛的林火蔓延模型有：FARSITE[12][13]、BehavePlus[14][15]、和HFire[16]。FARSITE和BehavePlus集成了現(xiàn)有的地表火、樹冠火、飛火和火加速等子模型，它們都采用Rothermel模型來計算地表火蔓延速度。不過，BehavePlus是基于橢圓形來模擬二維林火蔓延趨勢，而FARSITE是采用基于惠更斯波動理論來模擬二維林火蔓延趨勢。HFire是一種基于柵格的空間直觀模型，它也采用Rothermel模型來計算林火蔓延速率，運用惠更斯理論來模擬二維林火蔓延趨勢。

3.2 基于離散事件驅(qū)動的林火蔓延研究

最近，出現(xiàn)了采用離散事件來模擬林火蔓延。其中，比較成熟和完善的林火蔓延模型是DEVS-FIRE模型[17]。它將森林看成是由一個個大小相同的矩形單元組成的空間，空間的尺寸由GIS來確定。每個單元代表一個DEVS原子模型。一旦著火，每個單元都利用Rothermel模型來計算一維火蔓延速度和方向，然后根據(jù)橢圓形將其分解為八個方向的二維火蔓延。在林火蔓延的過程中，采用元胞自動機模型，分解的八個方向剛好與八個相鄰的元胞方向一致。元胞之間信息的交流采用預(yù)時間表（pre_schedule）策略[18]，也就是說，一旦點燃，元胞便立即同時發(fā)送八條信息給相鄰的元胞。在模擬過程中，有可能存在大量的細胞單元，然而只有一小部分細胞單元處于活動狀態(tài)。

基于內(nèi)存空間問題的考慮，DEVS-FIRE采用了動態(tài)結(jié)構(gòu)細胞空間模型，這使得在仿真過程中可以根據(jù)實際情況動態(tài)地創(chuàng)建和移除細胞單元，從而節(jié)省了大量的內(nèi)存空間。由于林火蔓延過程是動態(tài)的、復(fù)雜的、非線性和非高斯的，對于數(shù)據(jù)的獲取不可能是準確無誤的，林火蔓延計算的抽象性也會帶來誤差，故DEVS-FIRE引入了動態(tài)數(shù)據(jù)驅(qū)動應(yīng)用系統(tǒng)（DDDAS）[20]，對仿真結(jié)果進行動態(tài)地調(diào)整，并結(jié)合粒子濾波算法（PF）[19]對從傳感器獲取的實時數(shù)據(jù)進行數(shù)據(jù)同化來降低這些誤差。從傳感器獲取的測量數(shù)據(jù)具有高度的空間相關(guān)性，如果將其孤立的看待，將導(dǎo)致原始粒子權(quán)重偏差較大，對系統(tǒng)狀態(tài)分布的估計正確率降低。為此，DEVS-FIRE采用一種傳感器信息相關(guān)的方法。

4.林火蔓延研究展望

目前林火蔓延模型都是針對林火的始發(fā)階段，今后對森林大火和特大火的模型研究將成為重點[1]。在DEVS-FIRE模型的基礎(chǔ)上，我們可以嘗試用其他方法對數(shù)據(jù)進行同化來降低估計誤差。虛擬現(xiàn)實地理信息系統(tǒng)技術(shù)（VR-GIS）的發(fā)展，為林火蔓延模型的研究向三維、立體可視化的方向發(fā)展提供了有力的技術(shù)保障。物聯(lián)網(wǎng)技術(shù)的興起，為林火行為的預(yù)測、監(jiān)控等，向智能化方向發(fā)展提供了保障。WebGIS和OpenGIS的發(fā)展，使得林火蔓延仿真向網(wǎng)格化方向發(fā)展，成為一種多用戶端的系統(tǒng)[21]。地理信息系統(tǒng)（GIS）、衛(wèi)星遙感（RS）和全球定位系統(tǒng)（GPS）的集成為獲取更全面準確的數(shù)據(jù)信息提供了保障。為此，我們在進行林火行為模擬數(shù)據(jù)采集過程中，引入其他類型的實時數(shù)據(jù)（如航空影像）。

5.結(jié)論

從已有的林火模型研究來看，目前要建立一個理想的林火蔓延模型是比較困難的，不管是物理模型還是統(tǒng)計模型，都要對可燃物、氣象條件和地形等因子進行統(tǒng)計意義上的估計[1]。而這些因子隨地域的不同而存在多樣性和不確定性。由于林火行為的復(fù)雜性，目前還沒有一個真正具有普遍適用性的林火蔓延模型誕生。我國對林火蔓延的研究多是建立在美國的Rothermel模型的基礎(chǔ)上，并對其進行不斷的改進，對其存在的局限性無實質(zhì)性的突破。結(jié)合我國森林生態(tài)系統(tǒng)的實際情況和自然特點，通過大量的實地觀察與建模，獲取更全面真實的參數(shù)信息，構(gòu)建適合我國的林火蔓延模型，將是今后研究的重要方向。

參考文獻

[1]唐曉燕.林火蔓延模型及蔓延模擬的研究進展[J].北京林業(yè)大學(xué)學(xué)報，2002，24（1）：87-91.

[2]王正非.山火初始蔓延速度測算法[J].山地研究，1983，

1（2）：42-51.

[3]王正非.通用森林火險等級系統(tǒng)[J].自然災(zāi)害學(xué)報，1992，1（3）：

39-44.

[4]羅光裕.加拿大林火行為預(yù)報系統(tǒng)進展[J].森林防火，1994，1（40）：44-46.

[5]Richard C.Rothermel a mathematical model for predicting fire spread in wildland fuels[J].USDA：Forest Service Research Paper INT-115，1972.

[6]袁春明.林火行為研究概況[J].世界林業(yè)研究，2000，13（6）：27-30.

[7]景文峰.林火蔓延預(yù)測計算機圖形顯示系統(tǒng)[J].火災(zāi)科學(xué)，1993，2（2）：44-51.

[8]王海暉.森林地表火行為估算的數(shù)學(xué)模型[J].火災(zāi)科學(xué)，1994，3（1）：33-41.

[9]朱啟疆.GIS支持下森林火場蔓延的空間模擬[J].環(huán)境遙感，1995，10（2）：81-87.

[10]李建微.虛擬森林景觀中林火蔓延模型及三維可視化表達[J].應(yīng)用生態(tài)學(xué)報，2005，16（5）：838-842.

[11]蔡學(xué)理，張貴.基于ArcInfo技術(shù)的森林火災(zāi)指揮撲救決策方案[J].湖南林業(yè)科技，2005，32（2）：17-19.

[12]吳志偉.FARSITE火行為模型的原理、結(jié)構(gòu)及其應(yīng)用[J].生態(tài)學(xué)雜志，2012，31（2）：494-500.

[13]Mark A.Finney FARSITE：Fire Area Simulator-Model Development and Evaluation[R].Research Paper RMRS-RP-4，US Department of Agriculture，F(xiàn)orest Service，Rocky Mountain Research Station，Ogden，UT，1998.

[14]Patricia L.Andrews BehavePlus FIRE MODELING SYSTEM：PAST，PRESENT，AND FUTURE US Forest Service[R].Rocky Mountain Research Station，Missoula，Montana.

[15]Andrews PL.BEHAVE：Fire Behavior Prediction and Fuel Modeling System-BURN Subsystem[R].Part 1.General Technical Report INT-194，US Department of Agriculture，F(xiàn)orest Service，Intermountain Research Station，Ogden，UT，1986.

[16]Morais M.Comparing Spatially Explicit Models of Fire Spread Through Chaparral Fuels：A New Model Based Upon the Rothermel Fire Spread Equation[D].MA Thesis，University of California，Santa Barbara，CA，2001.

[17]Ntaimo.DEVS-FIRE：Towards an Integrated Simulation Environment for Surface Wildfire Spread and Containment[J].SIMULATION，2008，84（4）：137-155.

[18]Xiaolin Hu.DEVS-FIRE：design and application of formal discrete event wildfire spread and suppression models[J].SIMULATION，2012.

[19]H.Xue，F(xiàn).Gu，X.Hu，Data Assimilation Using Sequential Monte Carlo Methods in Wildfire Spread Simulation[J].The ACM Transactions on Modeling and Computer Simulation（TOMACS），Vol.22，No.4，Article No.23，2012.

[20]H.Xue，X.Hu，Exploiting Sensor Spatial Correlation for Dynamic Data Driven Simulation of Wildfire[C].Proc.ACM/IEEE/SCS Workshop on Principles of Advanced and Distributed Simulation（PADS），2012.

[21]張運生.林火蔓延模擬研究進展[J].安徽農(nóng)業(yè)科學(xué)，2010，

38（32）：18208-18209.

項目資助：國家948項目（2014-09）;湖南省科技計劃項目（2013WK3011）。

作者簡介：

廖超（1988—），女，湖南望城人，碩士研究生，研究方向：林火蔓延仿真。

陳愛斌（1971—），男，湖南攸縣人，博士，教授，研究方向：圖像處理，林業(yè)信息化，森林防火。