沈威

(遼寧省林業(yè)調(diào)查規(guī)劃監(jiān)測院，遼寧 沈陽 110122)

研究用林分平均木的材積生長率代替林分蓄積生長率可能是未來的發(fā)展趨勢。林木材積生長率加入樹高因子，通過胸徑、樹高、年齡計(jì)算林木材積生長率將成為一種可能[1]。本文旨在研究利用日本落葉松材積生長率模型計(jì)算出小班平均木的材積生長率，代替林分蓄積生長率，林分蓄積乘以林分蓄積生長率計(jì)算出林分蓄積生長量，從而可預(yù)估林分未來幾年的林分蓄積[2-5]。

1 選取樣地

收集遼寧省森林資源連續(xù)清查1990年至2015年所有樣地信息，篩選出含有日本落葉松的固定樣地30個(gè)，使用EXCEL軟件中Vlookup函數(shù)查找相同樣地號的五個(gè)連查間隔期內(nèi)同號木的檢尺類型、直徑值，組成列標(biāo)題為樣地號、樣木號，6個(gè)檢尺類型、6個(gè)直徑的EXCEL表格。

2 計(jì)算林分蓄積生長率、林木材積生長率

選取5個(gè)間隔內(nèi)，保留木齊全的固定樣地30個(gè)，查一元立木材積表計(jì)算每個(gè)固定樣地蓄積。計(jì)算樣地平均木材積。

利用普雷斯勒公式：

計(jì)算蓄積或材積的連年生長率。

式中：yt為本期蓄積或材積，yt-n為前期蓄積或材積，n為間隔期5 a。計(jì)算結(jié)果詳見表1。

表1 日本落葉松材積和蓄積生長率

由表1可知，林分蓄積生長率和平均木的材積生長率比較接近，由此可以用林分平均木材積生長率代替林分蓄積生長率，把林分和林木緊密地聯(lián)系起來，研究林木材積生長率也就是研究林分蓄積生長率[6]。

3 適用性檢驗(yàn)

為檢驗(yàn)林分平均木材積生長率代替林分蓄積生長率的可能性，使用表1數(shù)據(jù)做檢驗(yàn)樣本，檢驗(yàn)方法用總相對誤差(RS)法和F檢驗(yàn)[7]，當(dāng)RS小于3%且大于-3%，同時(shí)通過F檢驗(yàn)(F≤F0.05),則可以用林分平均木材積生長率代替林分蓄積生長率，否則不能代替。

3.1 總相對誤差

計(jì)算結(jié)果：RS=0.49%，大于-3%，小于3%。

3.2F檢驗(yàn)

根據(jù)林分蓄積生長率yi和林分平均木材積生長率xi建立回歸直線方程yi=a+bxi，若林分蓄積生長率和林分平均木材積生長率無差異時(shí)，則直線截距a=0，斜距b=1，說明林分平均木材積生長率代替林分蓄積生長率適用。反之，則需在一定可靠性條件下，將點(diǎn)列回歸直線參數(shù)a、b與理想直線參數(shù)α=0，β=1進(jìn)行差異性檢驗(yàn)[8]。將檢驗(yàn)相關(guān)數(shù)據(jù)代入方程式得回歸系數(shù)估計(jì)值。

直線參數(shù)的F值檢驗(yàn)公式：

=1.34

式中m樣本個(gè)數(shù)，Pxi林分蓄積生長率，Pci林分平均木材積生長率 ，F(xiàn)服從自由度f1=2、f2=m-2-1的F分布。

當(dāng)F>F0.05時(shí)，推翻假設(shè)，說明a不等于0，b不等于1，檢驗(yàn)樣本實(shí)際值和估計(jì)值存在顯著差異，林分平均木生長率不能代替林分蓄積生長率。

當(dāng)F≤F0.05時(shí),說明a等于0，b等于1，檢驗(yàn)樣本實(shí)際值和估計(jì)值沒有顯著差異，林分平均木生長率能代替林分蓄積生長率。

經(jīng)計(jì)算得出F值，詳見表2。由表2可知F=1.34

表2 F檢驗(yàn)值

4 結(jié)論

通過計(jì)算總相對誤差小于3%，F(xiàn)檢驗(yàn)F≤F0.05，所以用林分平均木材積生長率代替林分蓄積生長率是可行的。把對林分蓄積生長率研究，轉(zhuǎn)移到對林分平均木材積生長率研究，解析林分平均木，就能研究胸徑生長率，樹高生長率，材積生長率。