季俊超,宮智耀,龐靚璠,張 燕,曲智林

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

?

一種非線(xiàn)性可燃物含水率模型的線(xiàn)性化方法*

季俊超,宮智耀,龐靚璠,張 燕,曲智林

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

利用泰勒展開(kāi)的的方法給出一種將非線(xiàn)性可燃物含水率模型線(xiàn)性化的方法，進(jìn)而給出了基于logistic函數(shù)可燃物含水率模型參數(shù)的估計(jì)方法，由此得到一種線(xiàn)性化的可燃物含水率預(yù)測(cè)模型.通過(guò)實(shí)驗(yàn)數(shù)據(jù)檢驗(yàn)，給出的線(xiàn)性化可燃物含水率模型的準(zhǔn)確率達(dá)到97%(相對(duì)誤差不超過(guò)5%)，說(shuō)明該非線(xiàn)性可燃物含水率模型的線(xiàn)性化方法是可行的，由此得到的線(xiàn)性化可燃物含水率預(yù)測(cè)模型可以用于預(yù)測(cè)可燃物含水率變化規(guī)律.

可燃物;含水率;預(yù)測(cè)模型

0 引言

可燃物含水率的大小可以用來(lái)判斷林火發(fā)生的難易程度，因此研究森林可燃物含水率變化規(guī)律對(duì)于控制林火蔓延速度和預(yù)報(bào)森林火險(xiǎn)都有非常重要的作用.通過(guò)建立數(shù)學(xué)模型來(lái)研究可燃物含水率變化規(guī)律是一種有效的方法，而以平衡含水率為基礎(chǔ)建立的可燃物含水率模型一般都是微分方程，模型中的平衡含水率有的是用線(xiàn)性函數(shù)描述的，如曲智林[1]等，有的是用非線(xiàn)性函數(shù)描述的，如Byram(1963)[2]，Viney(1991)[3-4]，Simard[5]，Van Wagner[5]，Nelson[6]，Anderson[7]，金森(1999)[8-9],曲智林[10]等，而用非線(xiàn)性函數(shù)描述的可燃物函數(shù)率更加真實(shí)可靠，但在實(shí)際應(yīng)用中模型的參數(shù)是很難估計(jì)的，這對(duì)于預(yù)測(cè)可燃物含水率變化帶來(lái)了困難.該文將針對(duì)基于logistic 函數(shù)的這種非線(xiàn)性可燃物含水率模型給出其線(xiàn)性化的方法，進(jìn)而給出估算模型參數(shù)的方法.這些研究對(duì)于非線(xiàn)性可燃物含水率模型的應(yīng)用是非常有意義的.

1 理論推導(dǎo)

1.1 線(xiàn)性化推導(dǎo)

根據(jù)文獻(xiàn)[5]的研究表明，森林可燃物含水率的變化模型為:

(1)

其中: Mt為t時(shí)刻的可燃物含水率(%)；E0為t時(shí)刻的可燃物平衡含水率(%)；k為模型參數(shù).

根據(jù)以往研究：可燃物的平衡含水率主要與大氣中的相對(duì)濕度、風(fēng)速和溫度有關(guān)[11],即

Et=f (Tm,Hn,Wl)

(2)

其中:T為溫度(℃)；H為大氣中的相對(duì)濕度(%)；W為風(fēng)速(級(jí)).將連續(xù)模型(1)式進(jìn)行離散化變形，得

(3)

令Δt =1，則可以得到可燃物含水率的預(yù)測(cè)模型

Mt+1=(1-k)Mt+kEt

(4)

不難發(fā)現(xiàn)，想要將非線(xiàn)性可燃物含水率模型(4)線(xiàn)性化只需將Et展為線(xiàn)性形式，于是，將(2)式在點(diǎn)(0,0,0)處展開(kāi)，有：

(5)

再代回到模型(4)中得到

(6)

至此，將非線(xiàn)性化模型化為線(xiàn)性模型.

1.2 參數(shù)估計(jì)

若己知可燃物含水率和氣象因子的實(shí)時(shí)觀測(cè)數(shù)據(jù)(Mt+1,i，Mt,i，Tm,j,Hn,j,W1,j),(i=1,2,…,n),m時(shí)刻氣溫，n時(shí)刻的相對(duì)濕度和l時(shí)刻的風(fēng)速來(lái)預(yù)測(cè)t+1時(shí)刻的可燃物含水率.令

(7)

根據(jù)線(xiàn)性最小二乘原理，由正規(guī)方程組的解就可作為模型(6)中參數(shù)的估計(jì)值.

2 模型的驗(yàn)證

2.1 數(shù)據(jù)的選取

該文進(jìn)行研究所利用的觀察數(shù)據(jù)的來(lái)源為2009年3月20日至2009年4月20日黑龍江海林市石河子林場(chǎng)試驗(yàn)地(北緯43°28＇～45°69＇,東經(jīng)128°3＇～130°47＇ ).采集的數(shù)據(jù)包括時(shí)間、溫度、風(fēng)速、相對(duì)濕度、可燃物含水率、風(fēng)向以及降水量等數(shù)據(jù).

通過(guò)統(tǒng)計(jì)林場(chǎng)試驗(yàn)地的觀測(cè)數(shù)據(jù)，得到了各影響因子的取值范圍，見(jiàn)表 1.

處理數(shù)據(jù)均使用STATISTICA6．0軟件完成.

表1各影響因子取值范圍

2.2 基于logistic曲線(xiàn)的可燃物含水率線(xiàn)性化模型

根據(jù)文獻(xiàn)[12]的研究表明：平衡含水率的變化規(guī)律遵循logistic曲線(xiàn)，即溫度、相對(duì)濕度等氣象條件決定可燃物平衡含水率最高接近極限值Emax,最小接近或者可以達(dá)到0.因此選擇的平衡含水率為

(8)

其中：Mt為t時(shí)刻的可燃物含水率(%)；Et為t時(shí)刻的可燃物平衡含水率(%)；k為模型參數(shù)；Emax為平衡含水率的最大值，該文Emax=100,A,a,b,c為待定參數(shù).將上式代入到公式(4)中，即可以得到可燃物含水率預(yù)測(cè)模型

(9)

由1.1中線(xiàn)性化推導(dǎo)過(guò)程可知，為了將上述方程(2)展開(kāi)成線(xiàn)性關(guān)系,只需對(duì)Et進(jìn)行泰勒展開(kāi)為線(xiàn)性形式，

≈(A-2A2)(aTt-k1+bHt-k2+cWt-k3)+A2-A+1.

(10)

將該結(jié)果代入(9)，則得到非線(xiàn)性可燃物含水率模型線(xiàn)性化后的含水率模型為：

Mt+1≈(1-k)Mt+kEmax[(A-2A2)(aTt-k1+bHt-k2+cWt-k3)+A2-A+1]

(11)

利用此模型，t+1時(shí)刻的可燃物含水率可由t時(shí)刻的可燃物含水率、t-k1時(shí)刻的氣溫、t-k2時(shí)刻的相對(duì)濕度和t-k3時(shí)刻的風(fēng)速算得.

2.3 模型中參數(shù)的估算方法

若己知可燃物含水率和氣象因子的實(shí)時(shí)觀測(cè)數(shù)據(jù)(Mt+1,i，Mt,i，Tt-k1,i,Ht-k2,i,Wt-k3,i),(i=1,2,…,n),t-k1時(shí)刻氣溫，t-k2時(shí)刻的相對(duì)濕度和t-k3時(shí)刻的風(fēng)速來(lái)預(yù)測(cè)t+1時(shí)刻的可燃物含水率.根據(jù)非線(xiàn)性最小二乘原理，令

得到方程組

利用STATISTIC6.0軟件解出此正規(guī)方程組的解就可作為模型(11)中的參數(shù) 的估計(jì)值.

2.4 模型的建立

2.4.1 可燃物單位時(shí)間內(nèi)含水率改變量與各影響因子的相關(guān)性分析

為了確定模型中的反應(yīng)時(shí)間k1、k2、k3，可以先對(duì)試驗(yàn)地各個(gè)觀測(cè)變量進(jìn)行相關(guān)性分析，該文利用試驗(yàn)地3月20日至3月28日的觀測(cè)數(shù)據(jù)(共200個(gè)樣本)并結(jié)合偏相關(guān)系數(shù)公式：

(13)

確定了單位時(shí)間內(nèi)可燃物含水率變化量與之前各個(gè)時(shí)段內(nèi)空氣中溫度、相對(duì)濕度和風(fēng)速等影響因子的相關(guān)性關(guān)系.

由表2中可得：該可燃物1h內(nèi)的含水率改變量與前2 h的氣溫值相關(guān)性較高，并呈負(fù)相關(guān)；與前1 h的相對(duì)濕度相關(guān)性較高,并呈正相關(guān)；與前3 h的風(fēng)速值相關(guān)性最高，并呈負(fù)相關(guān).由此可以說(shuō)明，氣象要素中的氣溫、相對(duì)濕度和風(fēng)速對(duì)可燃物含水率變化的影響有一定的滯后性，因此，針對(duì)該文的研究所要選擇的氣象因子為前3 h的氣溫、前1 h的相對(duì)溫度和前3 h的風(fēng)速，即公式(11)中k1=2,k2=1,k3=3.

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

2.4.2 模型參數(shù)估算與模型驗(yàn)證

該文依據(jù)林場(chǎng)試驗(yàn)地的200個(gè)觀測(cè)數(shù)據(jù)建立了模型，并根據(jù)線(xiàn)性最小二乘原理，按照該文1.2中方法步驟，對(duì)模型(11)的參數(shù)進(jìn)行估算(見(jiàn)表3).以上模型均通過(guò)t檢驗(yàn)(α=0.05).

表3 模型中參數(shù)的估計(jì)值

即得到可燃物含水率預(yù)測(cè)模型：

Mt+1=0.926Mt+0.74×100×[0.044×

(-0.043Tt-2+0.023Ht-1-0.653Wt-3)+0.752]

(14)

模型(14)表示可以用t時(shí)刻的可燃物含水率、t-2時(shí)刻氣溫、t-1時(shí)刻的相對(duì)濕度和t-3時(shí)刻風(fēng)速來(lái)預(yù)測(cè)t+1時(shí)刻的可燃物含水率.

為了說(shuō)明模型的精度，對(duì)模型(14)進(jìn)行了檢驗(yàn).檢驗(yàn)發(fā)現(xiàn)，選取200個(gè)試驗(yàn)地的觀測(cè)數(shù)據(jù)，其正確率達(dá)到97%(相對(duì)誤差不超過(guò)5%).并給出了林場(chǎng)試驗(yàn)地2009年3月28日可燃物含水率的實(shí)測(cè)值、預(yù)測(cè)值和誤差值(見(jiàn)表4)

由表4中的數(shù)據(jù)可見(jiàn)，該文建立的模型能夠較準(zhǔn)確地對(duì)試驗(yàn)地的可燃物含水率進(jìn)行預(yù)測(cè)，說(shuō)明該模型是可行的，能夠在實(shí)際中進(jìn)行應(yīng)用.

3 結(jié)束語(yǔ)

在基于 logistic曲線(xiàn)建立可燃物含水率模型的基礎(chǔ)上，利用泰勒展開(kāi)實(shí)現(xiàn)非線(xiàn)性可燃物含水率模型的線(xiàn)性化，并利用最小二乘法對(duì)模型參數(shù)進(jìn)行估計(jì)，經(jīng)檢驗(yàn)發(fā)現(xiàn)，模型的準(zhǔn)確精度達(dá)97% (相對(duì)誤差不超過(guò)5%)，說(shuō)明在多時(shí)無(wú)雨條件下，利用該模型可以預(yù)測(cè)可燃物含水率的變化,并且一定程度上解決了平衡含水率的參數(shù)估算的困難，證明了所采用的方法是可行的.但對(duì)模型參數(shù)進(jìn)行估計(jì)時(shí)，在求解正規(guī)方程組時(shí)可能會(huì)出現(xiàn)多個(gè)解，因此在利用STATISTICA6.0軟件求解方程組時(shí)應(yīng)給定好初始值.另外，應(yīng)指出模型建立在多時(shí)無(wú)雨的情況下，具有一定的局限性；選用的影響因素只有溫度、相對(duì)濕度和風(fēng)速，對(duì)其他因素沒(méi)有考慮；同時(shí)，該文是針對(duì)特定可燃物進(jìn)行研究，對(duì)于不同類(lèi)型可燃物的含水率模型是不同的，但方法是可以借鑒的.

表4 2009年3月28日試驗(yàn)地可燃物含水率的實(shí)測(cè)值、預(yù)測(cè)值和誤差值

[1] 曲智林，吳娟，閔盈盈.具有時(shí)滯的可燃物含水率預(yù)測(cè)模型[J]．東北林業(yè)大學(xué)學(xué)報(bào)，2012，40(3):120-122．

[2] Byram G B.An analysis of the drying process in forest fuel material[J].Unpublished report,USDA Forest Service,Southern Forest Fire Laboratory,Macon,GA,USA,1963.

[3] Viney N R,Catchpole E A.Estimating fuel moisture response time from field observations [J].International of Wildland Fire,1991,1(4):211-214.

[4] Viney N R.A Review of Fine Fuel Moisture Modeling[J].International Journal of Wildland Fire,1991,1:215-234.

[5] 劉曦，金森.平衡含水率法預(yù)測(cè)可燃物含水率的研究進(jìn)展[J]．林業(yè)科學(xué),2007,43(12):126-133．

[6] 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．

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

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

[9] 金森，李亮.時(shí)滯和平衡含水率直接估計(jì)法的有效性分析[J]．林業(yè)科學(xué)，2010,46(2):96 -102．

[10] 曲智林，李昱燁，閔盈盈.可燃物含水率實(shí)時(shí)變化的預(yù)測(cè)模型[J]．東北林業(yè)大學(xué)學(xué)報(bào),2010,38(6):66-67.

[11] Nelson R M. Prediction of diurnal change in 102hour fuel moisture content[J]. Canadian Journal of Forest Research，2000，30:1071-1087.

[12] 尹坤.基于logistic曲線(xiàn)的森林可燃物含水率模型[D].東北林業(yè)大學(xué)，2015.

(責(zé)任編輯：季春陽(yáng))

A Linearization Method for Nonlinear Fuel Moisture Content Model

Ji Junchao,Gong Zhiyao,Pang Liangfan,Zhang Yan,Qu Zhilin

(Northeast Forestry University)

In this paper, a linearization method for nonlinear fuel moisture content model is presented by Taylor’s theory，and then the method of estimating parameters from the fuel moisture content model based on the logistic function is given, thus a linear model for predicting the fuel moisture content is obtained.Through experimental data testing,the established linear model had an average precision of 97%(the relative error is no more than 5%),which means that the linearization method of the nonlinear fuel moisture content model is feasible,and the resulting linear fuel moisture content prediction model can be used to predict the changing patterns of fuel moisture content.

Forest fuels; Moisture content; Prediction model

2016-01-21

*東北林業(yè)大學(xué)大學(xué)生創(chuàng)新創(chuàng)業(yè)訓(xùn)練計(jì)劃項(xiàng)目(201510225159)

S762.1

A

1000-5617(2016)02-0065-04