一種野外細(xì)小死可燃物含水率預(yù)測方法*

尹 坤，曲智林

(東北林業(yè)大學(xué))

*黑龍江省博士后科研啟動金資助項目(LBH-Q12175)

0 引言

森林可燃物含水率是影響林火發(fā)生及林火蔓延的重要因素之一，因此研究可燃物含水率的變化規(guī)律對提高林火預(yù)測預(yù)報的準(zhǔn)確性和控制林火蔓延都具有極大意義．目前研究森林可燃物含水率變化規(guī)律一般采用建立數(shù)學(xué)模型的方法，多數(shù)采用的方法是平衡含水率法和氣象要素回歸法．平衡含水率法就是利用平衡含水率建立含水率模型，如 Viney［1］模型、Simard［2］模型、Van Wagner［2］模型、Anderson［3-4］模型、Nelson［5］模型等，以及金森［6-8］、曲智林［9-10］等均利用平衡含水率法給出了可燃物含水率的預(yù)測模型，由于平衡含水率不容易估算，因此給該類模型的應(yīng)用帶來一定的難度．在實(shí)際中氣象要素回歸法顯得更加直接、有效，氣象要素回歸法就是利用回歸理論直接建立氣象要素與可燃物含水率之間的函數(shù)關(guān)系．目前已有一些研究，如杜秀文等［11］，居恩德等［12］均對幾種森林類型可燃物含水率與氣象因子關(guān)系進(jìn)行了分析．但根據(jù)野外細(xì)小死可燃物含水率的實(shí)時觀測數(shù)據(jù)直接建立氣象要素與可燃物含水率之間的關(guān)系預(yù)測誤差較大，因此，該文根據(jù)試驗(yàn)地的觀測數(shù)據(jù)，給出了一個具有延遲效應(yīng)，且分段細(xì)小死可燃物含水率預(yù)測模型，有效地提高了預(yù)測精度，從而可為準(zhǔn)確預(yù)報林火發(fā)生及林火蔓延提供理論依據(jù)．

1 數(shù)據(jù)來源與研究方法

1．1 數(shù)據(jù)來源

該文研究所用的數(shù)據(jù)為2008年4月1日到2008年4月29日國家林業(yè)局在黑龍江省韓家園子林業(yè)局(東經(jīng) 125°10'～ 127°，北緯 50°55'～52°19')設(shè)置的觀測地的觀測數(shù)據(jù)．觀測木為長20 cm、外直徑為2．6 cm、內(nèi)直徑為2．1 cm 的內(nèi)空的椴木．觀測木安裝在觀測儀器上．觀測儀器每1 h自動收錄各項數(shù)據(jù)．收錄的數(shù)據(jù)包括時間、氣溫、相對濕度、風(fēng)速、風(fēng)向、降雨量、觀測木的含水率和溫度等．

該文所涉及的可燃物的含水率為相對含水率，即

相對含水率=(濕材重量-全干材重量)/濕材重量×100%

根據(jù)試驗(yàn)地的觀測數(shù)據(jù)，統(tǒng)計出各影響因子的取值范圍見表1．

表1 各影響因子取值范圍

1．2 模型建立

由于降雨量的多少對可燃物含水率的變化有較大影響，且它們之間的關(guān)系較為復(fù)雜．因此該文研究的模型建立在多時無雨的情況下，并在固定區(qū)域和特定時段內(nèi)采集數(shù)據(jù)，因此可燃物含水率預(yù)測模型的影響因子選擇大氣的溫度、相對濕度和風(fēng)速．相關(guān)性分析表明:氣象因子對單位時間內(nèi)可燃物含水率改變量的影響有一定的滯后性，即t時刻單位時間內(nèi)可燃物含水率改變量由t時刻可燃物含水率、t-k1時刻大氣中的溫度、t-k2時刻相對濕度和t-k3時刻風(fēng)速所決定．該文選取

式中:ΔMt為t時刻單位時間內(nèi)可燃物含水率改變量(%);Mt為t時刻的可燃物含水率(%);Tt-k1為t-k1時刻溫度(℃);Ht-k2為大氣中tk2時刻的相對濕度(%);Wt-k3為t-k3時刻風(fēng)速(級);ε 為隨機(jī)誤差項;k1、k2、k3、a、b、c為待定系數(shù)．

t+1時刻的可燃物含水率為:

則得到可燃物含水率變化預(yù)測模型

其中:Mt+1為t+1時刻的可燃物含水率(%)．

利用此模型可由t-k1時刻的氣溫值、t-k2時刻的相對濕度值與t-k3時刻的風(fēng)速值以及t時刻可燃物的含水率得到t+1時刻可燃物的含水率．該文處理數(shù)據(jù)均使用STATISTICA6．0軟件完成．

2 結(jié)果與分析

2．1 相關(guān)性分析

利用試驗(yàn)地2008年4月1日至4月29日的觀測數(shù)據(jù)(剔除受降雨影響的數(shù)據(jù)，共6700000004個樣本)和偏相關(guān)系數(shù)公式:

給出了可燃物1 h內(nèi)含水率改變量與之前各時段空氣中溫度、相對濕度和風(fēng)速等各影響因子的相關(guān)性關(guān)系．圖1-3給出了不同時刻的相關(guān)關(guān)系．

從表2可以看出:該觀測物1 h含水率改變量與之前3 h的氣溫值相關(guān)性較高，并且為負(fù)相關(guān);與之前2 h的相對溫度相關(guān)性較高，并且為正相關(guān);與之前3 h的風(fēng)速值相關(guān)性最高，并且為負(fù)相關(guān);說明氣溫、相對濕度和風(fēng)速對可燃物含水率變化影響有一定的滯后性．

表2 單位時間內(nèi)含水率改變量與不同時段溫度、相對濕度和風(fēng)速的相關(guān)性

圖1 單位時間內(nèi)可燃物含水率改變量與t-k3時刻氣溫的關(guān)系

圖2 單位時間內(nèi)可燃物含水率改變量與t-k2時刻相對濕度的關(guān)系

2．2 分段可燃物含水率預(yù)測模型的建立

將選取的580個數(shù)據(jù)平均分成6段(30%～36%、37% ～43%、44% ～50%、51% ～57%、58% ～64%、65% ～71%)，根據(jù)統(tǒng)計分析理論，分別在每段上建立模型即公式(3)，并對其進(jìn)行參數(shù)估計，見表3．

圖3 單位時間內(nèi)可燃物含水率改變量與t-k3時刻風(fēng)速的關(guān)系

2．3 模型的檢驗(yàn)

利用94個檢驗(yàn)數(shù)據(jù)對分段可燃物含水率模型進(jìn)行檢驗(yàn)，分段可燃物含水率取平均值預(yù)測模型正確率到達(dá)98%(相對誤差不超過5%)．并給出了試驗(yàn)地2008年4月29日可燃物含水率的實(shí)測值與預(yù)測值見表4．

表3 分段可燃物含水率預(yù)測模中的參數(shù)估計值

表4 2008年4月29日試驗(yàn)地可燃物含水率的實(shí)測值與預(yù)測值

3 結(jié)論與討論

該文給出了在多時無雨條件下的分段可燃物含水率預(yù)測模型，利用該模型，根據(jù)某時刻的可燃物含水率和之前某時刻的氣溫、相對濕度和風(fēng)速能夠預(yù)測下一時刻的可燃物含水率，經(jīng)檢驗(yàn)?zāi)Ｐ偷臏?zhǔn)確率達(dá)到98%，說明這種分段預(yù)測方法是可行的．該文的研究是以野外細(xì)小死可燃物為研究對象，該研究方法也可以推廣到野外其它類型的可燃物．另外模型中的風(fēng)速因子使用的是級，若換成m/s預(yù)測效果將更好．該文沒有考慮降雨量的影響，以后將對這類模型繼續(xù)進(jìn)行研究．

［1］ Viney N R ，Catchpole E A．Estimating fuel moisture response time from field observations［J］．International of Wildland Fire，1991，1(4):211-214．

［2］ 劉曦，金森．平衡含水率法預(yù)測死可燃物含水率的研究進(jìn)展．林業(yè)科學(xué)，2007，43(12):126-133．

［3］ Anderson H E，Schuetle R D，Mutch R W．Timelag and equilibrium moisture content of Ponderosa pine needles［M］．Minneapolis:Dept．of Agriculture ，F(xiàn)orest Service ，Intermountain Forest and Range Experiment Station ，1978．8-14．

［4］ Anderson H E．Moisture diffusivity and response time in fine forest fuels［J］．Canadian Journal of Forest Research，1990，20(3):315-325．

［5］ Nelson R M ，Jr．Prediction of diurnal change in 10-h fuel stick moisture content［J］．Canadian Journal of Forest Research，2000，30(7):1071-1087．

［6］ 金森，姜文娟，孫玉英．用時滯和平衡含水率準(zhǔn)確預(yù)測可燃物含水率的理論算法［J］．森林防火，1999(4):12-14．

［7］ 金森，李銘堯，李有祥．幾種細(xì)小可燃物失水過程中含水率的變化規(guī)律［J］．東北林業(yè)大學(xué)學(xué)報，2000，28(1):35-38．

［8］ 金森，李亮．時滯和平衡含水率直接估計法的有效性分析［J］．林業(yè)科學(xué)，2010，46(2):96-102．

［9］ 曲智林，李昱燁，閔盈盈．可燃物含水率實(shí)時變化的預(yù)測模型［J］．東北林業(yè)大學(xué)學(xué)報，2010，38(6):66-71．

［10］曲智林，吳娟，閔盈盈．具有時滯的可燃物含水率預(yù)測模型［J］．東北林業(yè)大學(xué)學(xué)報，2012，40(3):120-122．

［11］杜秀文，李茹秀，王英杰．幾種森林類型可燃物含水率與氣象因子關(guān)系的分析［J］．東北林業(yè)大學(xué)學(xué)報，1988，16(5):88-90．

［12］居恩德，陳貴榮，王瑞君．可燃物含水率與氣象要素相關(guān)性的研究［J］．森林防火，1993，1:17-21．