金 森 褚騰飛 邸雪穎 劉礴霏 張吉利

(東北林業(yè)大學(xué)，哈爾濱，150040)

火行為是指林火發(fā)生發(fā)展過(guò)程中所表現(xiàn)出來(lái)的各種特征的總和，包括林火蔓延速度、火焰高度、駐留時(shí)間、反應(yīng)區(qū)寬度、可燃物消耗量和林火強(qiáng)度等指標(biāo)［1］。了解、掌握可燃物的林火行為，是做好森林火險(xiǎn)預(yù)報(bào)、提高森林火災(zāi)撲救效率和更好地分析、理解林火對(duì)森林生態(tài)系統(tǒng)影響的基礎(chǔ)［2］?；鹦袨槭芸扇嘉锝M成和理化性質(zhì)、可燃物床層特征及環(huán)境因素的影響，復(fù)雜性高。以火行為模型為主要任務(wù)的林火行為研究是林火科學(xué)的重要研究領(lǐng)域，國(guó)內(nèi)開(kāi)展得也相當(dāng)廣泛［3－9］?；鹦袨槟Ｐ鸵话憧煞譃槲锢?、半物理和實(shí)驗(yàn)統(tǒng)計(jì)模型3類［10］。物理模型普適性好，但要求較高的計(jì)算能力，在目前的技術(shù)條件下，還難以廣泛應(yīng)用［10］。在實(shí)際中得到廣泛應(yīng)用的都是在很大程度上或完全依賴于點(diǎn)燒試驗(yàn)的半物理模型和實(shí)驗(yàn)統(tǒng)計(jì)模型，如Rothermel模型［11］、加拿大火行為模型［12－14］、澳大利亞火行為模型［15－16］。這些半物理模型和實(shí)驗(yàn)統(tǒng)計(jì)模型，特別是后者，因使用很多來(lái)自點(diǎn)燒試驗(yàn)的方程和參數(shù)，對(duì)計(jì)算能力的要求不高，從而得到廣泛的應(yīng)用，但也恰恰是因?yàn)閷?duì)試驗(yàn)條件依賴大，影響了這些模型的推廣使用。在目前的經(jīng)濟(jì)技術(shù)條件下，對(duì)于一個(gè)地區(qū)而言，如果建立至少適用于該區(qū)域的試驗(yàn)統(tǒng)計(jì)火行為模型或半物理火行為模型，必須對(duì)該區(qū)域內(nèi)的主要可燃物類型的火行為進(jìn)行系統(tǒng)研究，了解、掌握它們?cè)诓煌扇嘉锎矊訔l件和環(huán)境條件下的火行為特征。

紅松(Pinus koraiensis)針闊混交林是溫帶森林生態(tài)系統(tǒng)的頂級(jí)群落，這些林下的可燃物由不同的針葉、闊葉組成，其火行為受這些組分的火行為影響［11，17］。了解、掌握這些組分的火行為，對(duì)于建立更準(zhǔn)確的混合可燃物模型十分必要。為此，筆者對(duì)該林分下的主要可燃物的組分進(jìn)行了實(shí)驗(yàn)室點(diǎn)燒試驗(yàn)，分析影響其火行為的主要因子。這些影響因子往往因可燃物組分的不同而有所變化［17］，需分別研究，然后才能建立相應(yīng)組分的火行為預(yù)測(cè)模型。為此，本文給出了紅松針葉床層的室內(nèi)點(diǎn)燒試驗(yàn)中有關(guān)林火蔓延速率研究的初步結(jié)果，其目的是:給出平地?zé)o風(fēng)條件下紅松針葉床層的蔓延速率的基本數(shù)據(jù);確定影響紅松針葉床層蔓延速率的主要因子，以及這些因子對(duì)林火蔓延速率的影響之間是否相關(guān);建立紅松針葉床層的林火蔓延速率模型。

1 研究方法

1.1 試驗(yàn)材料采集

紅松松針采于東北林業(yè)大學(xué)帽兒山試驗(yàn)林場(chǎng)老山(127°30'～127°34'E，45°20'～45°25'N)紅松(Pinus koraiensis)人工林。林場(chǎng)所處地域?qū)贉貛Т箨懶詺夂?。年平均氣?.8℃;1月份最冷，氣溫日變化范圍在－18～－23℃;7月份最熱，氣溫日變化范圍在21～22℃。年平均降水量723.8 mm，雨水集中在6—8月份。林分樹(shù)齡50 a，林分平均胸徑21.1 cm，平均樹(shù)高13.4 m。林下植被以毛榛子(Corylus mandshurica)、刺五加(Acanthopanax senticosus)、東北山梅花(Philadelphus schrenkii)、苔草(Carex sp.)為主。地表可燃物主要是凋落的紅松針葉，混有少量的闊葉落葉。

1.2 室內(nèi)點(diǎn)燒試驗(yàn)

點(diǎn)燒試驗(yàn)在東北林業(yè)大學(xué)帽兒山森林防火實(shí)驗(yàn)室內(nèi)進(jìn)行。2010年春季，在研究地區(qū)的紅松林下進(jìn)行紅松針葉采集，同時(shí)測(cè)定可燃物的載量、厚度等特征，在實(shí)驗(yàn)室內(nèi)模擬這些可燃物的野外條件，鋪設(shè)不同含水率、載量和高度組合的均勻的紅松針葉床層。其中可燃物含水率設(shè)5個(gè)水平:5%、10%、15%、20%、25%。載量設(shè) 5 個(gè)水平:4 、5、6、7、8 t·hm－2?？扇嘉锎矊痈叨仍O(shè)4個(gè)水平:3、5、7、9 cm。共鋪設(shè)5×5×4=100個(gè)不同的3因素交叉水平床層，進(jìn)行了100次點(diǎn)燒試驗(yàn)?？扇嘉锖什捎萌缦路椒ㄕ{(diào)制:先將可燃物均勻浸濕，然后將其分層放置在自制的大烘箱內(nèi)(3.2 m×1.9 m×2.3 m)，控制溫度進(jìn)行烘干，使含水率逐漸下降，采用高精度快速水分測(cè)定儀AND—ML50測(cè)定含水率，直到可燃物含水率達(dá)到設(shè)定的水平。根據(jù)預(yù)定的可燃物含水率和載量，計(jì)算出每次點(diǎn)燒所需可燃物的鮮質(zhì)量，然后稱量同樣質(zhì)量的可燃物進(jìn)行鋪設(shè)。燃燒床的尺寸為2 m×2 m。燃燒床水平放置，模擬平地?zé)o風(fēng)條件下的燃燒。在燃燒床一端固定有一點(diǎn)火槽，放入乙醇后點(diǎn)燃，可迅速形成一條火線。每次點(diǎn)燒點(diǎn)燃后，火頭呈一條直線向前蔓延。當(dāng)火蔓延達(dá)到“似穩(wěn)態(tài)”(Quasisteady state)時(shí)開(kāi)始記錄燃燒時(shí)間及火焰高度，采用標(biāo)桿法測(cè)量林火蔓延速率［18－20］，用攝像機(jī)在燃燒床正面及側(cè)面同時(shí)記錄點(diǎn)燒過(guò)程。

1.3 數(shù)據(jù)分析

1.3.1 紅松針葉床層林火蔓延速率的基本特征

通過(guò)對(duì)100次點(diǎn)燒試驗(yàn)的統(tǒng)計(jì)分析，給出紅松針葉床層平地?zé)o風(fēng)條件下林火蔓延速率的基本特征，如均值、中數(shù)、最大、最小值等。

1.3.2 紅松針葉床層林火蔓延速率的影響因子分析

研究［21］表明，影響林火蔓延的因子間存在一定的交互作用。一些因子對(duì)林火蔓延的影響往往被其他因子的作用所掩蓋。只有在保持其他因子不變而只有目的因子變化時(shí)，才能準(zhǔn)確揭示目的因子對(duì)林火的作用。因此，文中的分析包括兩個(gè)方面:一是全部點(diǎn)燒數(shù)據(jù)放在一起，即各因子都在變化時(shí)，分析某一因子對(duì)可燃物蔓延速率的影響，以判斷該因子對(duì)蔓延速率影響的強(qiáng)度;另一個(gè)是其他因子不變時(shí)，單獨(dú)剝離出該因子的作用。對(duì)兩者進(jìn)行比較，以分析該因子的影響與其他因子間的交互關(guān)系。

1.3.2.1 含水率對(duì)林火蔓延速率的影響

用全部點(diǎn)燒數(shù)據(jù)繪制紅松針葉含水率和林火蔓延速率的散點(diǎn)圖，并計(jì)算其相關(guān)系數(shù)。然后以載量(5個(gè)水平)和高度(4個(gè)水平)為分類變量，將100次點(diǎn)燒試驗(yàn)數(shù)據(jù)分為5×4=20組，每組點(diǎn)燒試驗(yàn)的可燃物床層的載量和高度相同而含水率不同，以揭示含水率對(duì)蔓延速率的獨(dú)立影響。分別繪制這20組數(shù)據(jù)的紅松針葉含水率和林火蔓延速率的散點(diǎn)圖，并計(jì)算其相關(guān)系數(shù)。比較分類數(shù)據(jù)的散點(diǎn)圖和全部數(shù)據(jù)的散點(diǎn)圖所表現(xiàn)出來(lái)的含水率和蔓延速率之間的關(guān)系，分析紅松針葉含水率對(duì)林火蔓延速率的影響及這種影響與針葉床層載量和高度的關(guān)系。

1.3.2.2 床層載量對(duì)林火蔓延速率的影響

與1.3.2.1同法，只是采用紅松針葉含水率(5個(gè)水平)和床層高度(4個(gè)水平)作為分類變量進(jìn)行數(shù)據(jù)分組，共20組分類數(shù)據(jù)。

1.3.2.3 床層高度對(duì)林火蔓延速率的影響

與1.3.2.1同法，只是采用紅松針葉含水率(5個(gè)水平)和床層載量(5個(gè)水平)作為分類變量進(jìn)行分組，共25組分類數(shù)據(jù)。

1.3.2.4 床層壓縮比對(duì)林火蔓延速率的影響

可燃物床層壓縮比對(duì)林火蔓延速率有影響，但壓縮比與可燃物床層載量和高度是聯(lián)系在一起的，在試驗(yàn)中無(wú)法單獨(dú)剝離處理。在1.3.2.2和1.3.2.3中，紅松針葉床層載量和高度對(duì)蔓延速率的影響分析中就隱含著壓縮比對(duì)林火蔓延速率的影響，且兩者的趨勢(shì)是一致的。為此，只用全部點(diǎn)燒試驗(yàn)數(shù)據(jù)繪制可燃物床層壓縮比和林火蔓延速率的散點(diǎn)圖，并計(jì)算相關(guān)系數(shù)，以綜合分析壓縮比對(duì)林火蔓延速率的影響。壓縮比采用計(jì)算公式為:β=L/(Dρp)。式中:β為壓縮比;ρp為可燃物顆粒密度(g·cm－3)，試驗(yàn)測(cè)定;L為床層載量(g·cm－2);D為床層高度(cm)。

1.3.3 建立紅松針葉床層林火蔓延速率的預(yù)測(cè)模型

式中:xi為第i個(gè)建模因子;fi為因子xi對(duì)林火蔓延速率的影響函數(shù);n為所采用的建模因子數(shù)量。

本文的加式模型采用線性回歸形式，要求各因子之間不能存在強(qiáng)烈的相關(guān)，因此，載量和壓縮比、高度和壓縮比不能同時(shí)出現(xiàn)在加式模型中。同時(shí)要考慮各因子之間的交互作用，若交互作用顯著，則將其加入模型中。對(duì)于所有的加式和乘式模型，都進(jìn)行100次的交叉驗(yàn)證，即每次選擇1次點(diǎn)燒數(shù)據(jù)作為驗(yàn)證數(shù)據(jù)，用其余99次點(diǎn)燒數(shù)據(jù)建模。計(jì)算這100次驗(yàn)證的平均絕對(duì)誤差(EM，A)和平均相對(duì)誤差(EM，R)。誤差的計(jì)算公式為:

式中:Ri和Ri'分別為實(shí)測(cè)和預(yù)測(cè)的紅松針葉床層的林火蔓延速率(m·min－1)。

2 結(jié)果與分析

2.1 紅松針葉床層林火蔓延速率的特征

表1給出了100次點(diǎn)燒試驗(yàn)的基本情況。其中紅松針葉床層的林火蔓延速率為0～0.399 m·min－1，平均值為0.2 m·min－1。從 75%的區(qū)間值看，多數(shù)紅松針葉點(diǎn)燒試驗(yàn)的林火蔓延速率小于0.3m·min－1。紅松針葉床層的壓縮比從 0.014 到 0.084。在100次點(diǎn)燒試驗(yàn)中，有9次沒(méi)有持續(xù)燃燒，其中6次試驗(yàn)的紅松針葉含水率為25%，其余3次的含水率為20%。

2.2 紅松針葉床層林火蔓延速率的影響因子

2.2.1 可燃物含水率的影響

圖1給出了用100次點(diǎn)燒試驗(yàn)數(shù)據(jù)繪制的紅松針葉床層林火蔓延速率和含水率之間的關(guān)系。從圖1中可見(jiàn)，含水率對(duì)紅松針葉床層林火蔓延速率具有顯著的阻滯作用。圖2給出了相同載量和高度時(shí)紅松針葉床層含水率與林火蔓延速率之間的關(guān)系，其中各小圖的擬合直線的斜率比較相似，表明紅松針葉床層林火蔓延速率和含水率之間的關(guān)系在不同的載量和高度的紅松針葉床層間是一致的，說(shuō)明紅松針葉含水率對(duì)林火蔓延速率的影響與紅松針葉床層的載量、高度無(wú)明顯關(guān)系。圖2中各小圖的壓縮比不同，這也說(shuō)明紅松針葉含水率對(duì)林火蔓延速率的影響?yīng)毩⒂诖矊訅嚎s比，這在2.3的建模過(guò)程中得到了驗(yàn)證。圖1和圖2的比較研究表明，紅松針葉含水率對(duì)紅松針葉床層林火蔓延速率的影響在其他床層特征一致和混合床層特征兩個(gè)條件下是一致的，即都表現(xiàn)為較強(qiáng)的線性關(guān)系。

表1 點(diǎn)燒試驗(yàn)和蔓延速率的統(tǒng)計(jì)數(shù)據(jù)(樣本數(shù)為100)

圖1 林火蔓延速率與可燃物含水率之間的關(guān)系(全部100次點(diǎn)燒試驗(yàn)數(shù)據(jù))

在本文的研究中，可燃物含水率對(duì)林火蔓延速率的影響是線性的，這與 Burrows［22］和 Baeza［23］的點(diǎn)燒試驗(yàn)結(jié)果相似。在Rothermel模型［11］中，可燃物含水率對(duì)林火蔓延速率的影響是用含水率的三次多項(xiàng)式和一次多項(xiàng)式的商來(lái)表示的，其他研究則用冪函數(shù)［12，16］或指數(shù)函數(shù)［24，15］來(lái)描述，與本文所得結(jié)果不同。

2.2.2 紅松針葉床層載量的影響

從圖3中可見(jiàn)，紅松針葉層載量對(duì)林火蔓延速率具有非常弱的影響，隨載量的增加，紅松針葉床層的林火蔓延速率略有增加，但未達(dá)到0.05的統(tǒng)計(jì)顯著水平。圖4給出了相同床層高度和含水率時(shí)載量與紅松針葉床層林火蔓延速率之間的關(guān)系，其中各小圖所表現(xiàn)出的床層載量與林火蔓延速率之間的關(guān)系比較復(fù)雜，在20個(gè)小圖中，只有4個(gè)達(dá)到了0.05以上的顯著水平，4個(gè)達(dá)到了0.10的顯著水平，而其他小圖的紅松針葉載量與林火蔓延速率沒(méi)有表現(xiàn)出線性關(guān)系。圖3和圖4的比較研究表明，對(duì)于紅松針葉床層而言，從統(tǒng)計(jì)上講，載量對(duì)林火蔓延速率沒(méi)有顯著的影響。但在一定的可燃物床層條件下(如圖4(d)、(l)、(t)蔓延速率與載量線性關(guān)系顯著)，林火蔓延速率隨載量增大而增加。這表明可燃物載量對(duì)林火蔓延速率的影響比較復(fù)雜:一種可能是兩者之間的關(guān)系依賴于其他條件，如含水率或床層高度或壓縮比等;另一種可能是兩者之間沒(méi)有關(guān)系，其中一些顯著的相關(guān)是由于試驗(yàn)誤差導(dǎo)致。

圖2 可燃物床層載量和高度相同條件下的林火蔓延速率與可燃物含水率之間的關(guān)系

圖3 林火蔓延速率與可燃物床層載量之間的關(guān)系(全部100次點(diǎn)燒試驗(yàn)數(shù)據(jù))

2.2.3 紅松針葉床層高度的影響

從圖5中可見(jiàn)，紅松針葉床層高度對(duì)林火蔓延速率具有較強(qiáng)的影響，隨可燃物床層高度的增加，紅松針葉床層的林火蔓延速率增加，顯著水平超過(guò)0.0001。圖6給出了相同可燃物載量和含水率時(shí)床層高度與紅松針葉床層林火蔓延速率之間的關(guān)系。其中，各小圖所表現(xiàn)出的可燃物床層高度與林火蔓延速率之間的關(guān)系也比較復(fù)雜:在25個(gè)小圖中，6個(gè)達(dá)到了0.05以上的顯著水平;5個(gè)達(dá)到了0.10的顯著水平;其他小圖中的可燃物高度雖與林火蔓延速率沒(méi)有表現(xiàn)出簡(jiǎn)單的相關(guān)關(guān)系，但從相關(guān)系數(shù)上看，其關(guān)系要遠(yuǎn)比載量與林火蔓延速率之間的關(guān)系(圖4)密切。因此，盡管從一定的可燃物床層條件下紅松針葉床層林火蔓延速率與高度關(guān)系不緊密，但從統(tǒng)計(jì)或大量試驗(yàn)數(shù)據(jù)上看，高度對(duì)林火蔓延速率具有顯著的影響。從圖6中還可以看出，當(dāng)可燃物含水率為5%和15%時(shí)，床層高度對(duì)林火蔓延速率的影響顯著，含水率為10%時(shí)則不顯著，而含水率為20%～25%則基本沒(méi)有影響。這表明，床層高度對(duì)林火蔓延速率的作用可能與含水率有關(guān)，含水率越高，床層高度的影響越不明顯，或者影響模式從線性向非線性過(guò)渡。

圖4 可燃物含水率和床層高度相同條件下的林火蔓延速率與可燃物載量之間的關(guān)系(圖4(e)為不持續(xù)燃燒數(shù)據(jù))

圖5 林火蔓延速率與可燃物床層高度之間的關(guān)系(全部100次點(diǎn)燒試驗(yàn)數(shù)據(jù))

2.2.4 紅松針葉床層壓縮比的影響

圖7給出了紅松針葉床層壓縮比和林火蔓延速率的關(guān)系。從圖7中可見(jiàn)，壓縮比對(duì)林火蔓延速率具有顯著影響。在可燃物含水率不變的條件下，當(dāng)載量和高度兩個(gè)變量中的一個(gè)保持不變而另一個(gè)變化時(shí)(即2.2.2 和 2.2.3 中的試驗(yàn))，可燃物床層壓縮比也在變化。因此，圖4和圖6也可以用來(lái)研究剝離紅松針葉床層壓縮比對(duì)林火蔓延速率的影響，只是其變化與可燃物床層載量或高度同步。紅松針葉床層壓縮比與林火蔓延速率之間的相關(guān)系數(shù)與高度和林火蔓延速率之間的相關(guān)系數(shù)在絕對(duì)值上接近，因此，對(duì)于紅松針葉床層而言，壓縮比和高度對(duì)林火蔓延速率的影響是相似的，高度對(duì)林火蔓延的影響應(yīng)通過(guò)改變可燃物床層的壓縮比而發(fā)生作用。

2.3 紅松針葉床層林火蔓延速率的預(yù)測(cè)模型

根據(jù)前面的分析，在加式模型(線性模型)中選擇兩個(gè)模型。第一個(gè)模型為林火蔓延速率R=b0+b1M+b2D(m·min－1)，其自變量為紅松針葉含水率M、紅松針葉床層高度D。由于2.2.1中的研究已表明紅松針葉含水率對(duì)林火蔓延速率的影響與床層的其他特征無(wú)關(guān)，所以，含水率和床層高度的交互作用沒(méi)有包括在模型中。另一個(gè)線性模型的自變量為紅松針葉含水率M、可燃物床層壓縮比B。與上同理，含水率與壓縮比的交互作用也沒(méi)有使用，故模型形式也是R=b0+b1M+b2D。實(shí)際上，在包括交互作用的上述模型中，交互作用的回歸系數(shù)不顯著，也需要從模型中剔除，這也進(jìn)一步證實(shí)了2.2.1中關(guān)于紅松針葉含水率對(duì)林火蔓延速率的影響?yīng)毩⒂诖矊犹卣鞯姆治?。乘式模型也選擇兩個(gè):第一個(gè)模型的自變量為可燃物含水率和床層高度，形為R=(b0+b1M)(b2+b3D):另一個(gè)模型的自變量為可燃物含水率和床層壓縮比，形為R=(b0+b1M)(b2+b3B)。

圖6 可燃物含水率和床層載量相同條件下的林火蔓延速率與可燃物床層高度之間的關(guān)系

圖7 林火蔓延速率與可燃物床層壓縮比之間的關(guān)系(全部100次點(diǎn)燒試驗(yàn)數(shù)據(jù))

表2給出了4個(gè)模型的估計(jì)參數(shù)、誤差和統(tǒng)計(jì)檢驗(yàn)。其中，參數(shù)和統(tǒng)計(jì)檢驗(yàn)為用100次點(diǎn)燒試驗(yàn)數(shù)據(jù)估計(jì)的結(jié)果，而誤差為100次交叉驗(yàn)證的結(jié)果。從表2中的決定系數(shù)和平均絕對(duì)誤差上看，同樣形式的模型，以可燃物含水率和床層高度為預(yù)測(cè)因子的模型要比以含水率和床層壓縮比為自變量的模型要好。采用同樣預(yù)測(cè)因子時(shí)，加式模型比乘式模型好。其中，效果最好的模型是以可燃物含水率和床層高度為預(yù)測(cè)因子的加式模型，盡管其平均相對(duì)誤差比相應(yīng)的乘式模型略大一些。圖8給出了該模型的預(yù)測(cè)值和實(shí)測(cè)值的對(duì)比，從中可見(jiàn)，該模型對(duì)紅松針葉床層林火蔓延速率的預(yù)測(cè)在高估和低估方向上的比例基本是一致的，因此，模型是無(wú)偏的。但預(yù)測(cè)偏低的幅度要比預(yù)測(cè)偏高的幅度大。

圖8 以含水率和床層高度為自變量的林火蔓延速率線性預(yù)測(cè)模型的預(yù)測(cè)值和實(shí)測(cè)值對(duì)比

表2 4個(gè)林火蔓延速率模型的參數(shù)和誤差

3 結(jié)論與討論

系統(tǒng)設(shè)計(jì)的室內(nèi)點(diǎn)燒試驗(yàn)結(jié)果表明，平地?zé)o風(fēng)條件下紅松針葉床層的林火蔓延速率不超過(guò)0.4 m·min－1。紅松針葉含水率、床層高度和壓縮比對(duì)林火蔓延速率具有顯著的影響。其中，含水率對(duì)林火蔓延速率的影響與可燃物床層的高度、載量等無(wú)顯著關(guān)系，而可燃物床層高度和壓縮比對(duì)林火蔓延的影響可能受可燃物床層特征的影響，需進(jìn)一步研究確定。床層載量對(duì)紅松針葉床層的林火蔓延速率影響不大。以可燃物含水率和床層高度為預(yù)測(cè)因子的線性林火蔓延速率預(yù)測(cè)模型能夠解釋85%的林火蔓延速率變差，模型的平均絕對(duì)誤差為0.03 m·min－1，平均相對(duì)誤差不超過(guò)13%。

研究中，有些點(diǎn)燒試驗(yàn)在20%的含水率時(shí)不燃，有些在25%的含水率時(shí)不燃。這與 Rothermel［11］研究中采用的固定熄滅含水率不同。滅絕含水率本身就是一個(gè)比較模糊的概念，并且許多環(huán)境因素和可燃物本身的屬性都可能影響滅絕含水率的取值［25］(本文中可能與可燃物床層的結(jié)構(gòu)有關(guān))。Wilson［26］的研究表明，可燃物含水率對(duì)林火蔓延的持續(xù)性影響與可燃物床層結(jié)構(gòu)有關(guān)。本文的研究與其相似，但由于不燃點(diǎn)燒試驗(yàn)的次數(shù)只有9次，還無(wú)法具體分析其中的聯(lián)系，需進(jìn)一步研究，以便更好地理解可燃物含水率對(duì)火行為的影響。

可燃物床層的載量、高度和壓縮比對(duì)林火蔓延速率的影響比較復(fù)雜，Anderson［27］的研究表明，林火蔓延速率與可燃物載量正相關(guān);Catchpole［28－29］認(rèn)為，林火蔓延速率與可燃物床層高度關(guān)系不大;而Wilson［26］則認(rèn)為，無(wú)風(fēng)時(shí)林火蔓延速度正比于床層高度的平方根。本文的研究與Wilson［26］的結(jié)論相似，與其他研究相反。一些野外點(diǎn)燒試驗(yàn)，如Anderson［27］，受點(diǎn)燒時(shí)條件限制，一些可燃物床層特征變量的變化幅度可能較小，所得結(jié)果無(wú)法體現(xiàn)該變量對(duì)林火蔓延速率的影響。Catchpole［28－29］的研究與本研究試驗(yàn)條件相似，但所得結(jié)論與本研究相反，需進(jìn)一步研究。

［1］鄭煥能.森林防火［M］.哈爾濱:東北林業(yè)大學(xué)出版社，1992:24.

［2］舒立福，田曉瑞，馬林濤.林業(yè)科學(xué)研究［J］.林火生態(tài)的研究與應(yīng)用，1999，12(4):422－427.

［3］柴瑞海，趙雨森，杜秀文.大興安嶺林區(qū)草甸火順風(fēng)蔓延模型的研究［J］.東北林業(yè)大學(xué)學(xué)報(bào)，1998，16(4):90－93.

［4］王海暉，朱霽平.森林地表火行為估算的數(shù)學(xué)模型［J］.火災(zāi)科學(xué)，1994，3(1):33－41.

［5］袁春明，文定元.林火行為研究概況［J］.世界林業(yè)研究，2000，12(6):27－31.

［6］鐘占榮，周建軍，鄒樣輝，等.山地林火蔓延模型的研究［J］.火災(zāi)科學(xué)，2001，10(2):83－87.

［7］唐曉燕，孟憲宇，易浩若.林火蔓延模型及蔓延模擬的研究進(jìn)展［J］.北京林業(yè)大學(xué)學(xué)報(bào)，2002，24(1):87－91.

［8］田曉瑞，王明玉，殷麗，等.大興安嶺南部春季火行為特征及可燃物消耗［J］.林業(yè)科學(xué)，2009，45(3):90－95.

［9］王明玉，李濤，任云卯，等.森林火行為與特殊火行為研究進(jìn)展［J］.世界林業(yè)研究，2009，22(2):45－49.

［10］Sullivan A L.Wildland surface fire spread modeling，1990－2007.2:Empirical and quasi-empirical models［J］.International Journal of Wildland Fire，2009，18(4):369－386.

［11］Rothermel R C.A mathematical model for predicting fire spread in wildland fuels［R］.USDA:Forest Service Research Paper，1972:115.

［12］Stocks B J，Lynham T J，Lawson B D，et al.The Canadian system of forest fire danger rating system:An overview［J］.The Forestry Chronicle，1989，65(6):450－457.

［13］Forestry Canada Fire Danger Group.Development and structure of the Canadian forest fire behavior prediction system［R］.Ottawa:Forestry Canada，Science and Sustainable Development Directorate Information Report，1992.

［14］Taylor S W，Alexander M E.Science，technology，and human factors in fire danger rating:the Canadian experience［J］.International Journal of Wildland Fire，2006，15(1):121－135.

［15］Cheney N P，Gould J S，Catchpole W R.Prediction of fire spread in grasslands［J］.International Journal of Wildland Fire，1998，8(1):1－13.

［16］Gould J S，McCaw W L，Cheney N P，et al.Project vesta:Fire in dry eucalypt forest:Fuel structure，fuel dynamics and fire behaviour［M］.［s.l.］:Csiro Publishing，2007.

［17］Viegas D X，Almeida M，Miranda A I，et al.Linear model for spread rate and mass loss rate for mixed-size fuel beds［J］.International Journal of Wildland Fire，2010，19(5):531－540.

［18］崔文彬，喬啟宇，牛樹(shù)奎，等.地表火燃燒參數(shù)及其火頭前影響區(qū)的實(shí)驗(yàn)研究［J］.北京林業(yè)大學(xué)學(xué)報(bào)，1998，20(5):22－27.

［19］Wotton B M，McAlpine R S，Hobbs M W.The effects of fire front width on surface fire behaviour［J］.International Journal of Wildland Fire，1999，9(4):247－253.

［20］Balbi J H，Santoni P A，Dupuy J L.Dynamic modeling of fire spread across a fuel bed［J］.International Journal of Wildland Fire，1999，9(4):275－284.

［21］Viegas D X，Matos P L P，Palheiro P.Slope and wind effects on fire spread［M］.Rotterdam:Forest Fire Research＆ Wildland Fire Safety，2002.

［22］Burrows N，Ward B，Robinson A.Fire behaviour in spinifex fuels on the Gibson Desert Nature Reserve，Western Australia［J］.Journal of Arid Environments，1991，20(2):189－204.

［23］Baeza M，De Luís M，Raventós J，et al.Factors influencing fire behaviour in shrublands of different stand ages and the implications for using prescribed burning to reduce wildfire risk［J］.Journal of Environmental Management，2002，65(2):199－208.

［24］McCaw L.Predicting fire spread in Western Australian mallee-heath shrubland［D］.Canberra:University of New SouthWales，1997.

［25］Jolly W M.Sensitivity of a surface fire spread model and associated fire behaviour fuel models to changes in live fuel moisture［J］.International Journal of Wildland Fire，2007，16(4):503－509.

［26］Wilson R A，Jr.Reexamination of Rothermel’s fire spread equations in no-wind and no-slope conditions［R］.USDA:Forest Service Research Paper，1990:434.

［27］Anderson H E，Brackebusch R P，Mutch RW，et al.Mechanisms of fire spread research progress report no.2［R］.USDA:Forest Service Research Paper，1966:28.

［28］Catchpole W R，Catchpole E A，Butler B W，et al.Rate of spread of free-burning fires in woody fuels in a wind tunnel［J］.Combustion Science and Technology，1998，131(1/6):1－37.

［29］Catchpole E A，Catchpole W R，Rothermel R C.Fire behavior experiments in mixed fuel complexes［J］.International Journal of Wildland Fire，1993，3(1):45－57.