趙玉恒

(黑龍江省森林與環(huán)境科學(xué)研究院,黑龍江 齊齊哈爾 161005)

樹木的生長(zhǎng)因樹種和品種的不同，均有其自身的生長(zhǎng)規(guī)律，相同樹種也會(huì)因品種不同而差異顯著。為研究相同樹種不同品種間存在的這種差異，篩選了大樹平均木進(jìn)行樹干解析，利用SPSS曲線參數(shù)估計(jì)法模塊，分別以樹高、直徑和林齡為自變量自動(dòng)擬合11種材積曲線模型[1]，選擇最適模型建立大樹材積生長(zhǎng)曲線方程，為生產(chǎn)應(yīng)用提供科學(xué)參考。

1 樹干解析及材積生長(zhǎng)量擬合

在富錦市東風(fēng)崗林場(chǎng)30 a生試驗(yàn)林內(nèi)，我們選擇東風(fēng)崗1號(hào)楊平均木進(jìn)行樹干解析，其樹高、胸徑和材積分別為27.54 m、27.60 cm和0.651 16 m3。根據(jù)樹干解析數(shù)據(jù)，我們利用SPSS曲線估計(jì)程序中的線性、對(duì)數(shù)、反向、二次、三次、復(fù)合、冪、S、增長(zhǎng)、指數(shù)和Logistic等11個(gè)參數(shù)估計(jì)法，分別以樹高、直徑和齡階為自變量對(duì)已觀測(cè)的樹木材積進(jìn)行擬合，見表1。

表1 東風(fēng)崗1號(hào)楊解析木生長(zhǎng)進(jìn)程調(diào)查匯總表

2 大樹材積生長(zhǎng)曲線方程建立

我們分別以樹高、直徑和林齡為自變量，利用SPSS曲線估計(jì)程序?qū)Σ姆e生長(zhǎng)曲線進(jìn)行擬合，其中Cubic(立方)曲線與實(shí)測(cè)結(jié)果曲線擬合效果最好。

2.1 以樹高為自變量對(duì)立木材積進(jìn)行擬合

以樹高為自變量，對(duì)立木材積進(jìn)行擬合，擬合結(jié)果見圖1-圖2。

圖1 以樹高為自變量獲得的12種曲線擬合圖

圖2以樹高為自變量獲得的擬合度最高的材積生長(zhǎng)曲線圖

最佳生長(zhǎng)曲線的確定決定系數(shù)R2為0.993，標(biāo)準(zhǔn)誤差為0.020，F(xiàn)值為1 296.07，P值為0.000。獲得以樹高為自變量的材積生長(zhǎng)模型表達(dá)式：

Y=0.041-0.022x+0.002x2-3.133×10-5x3

式中:Y表示材積，X表示樹高。

2.2 以直徑為自變量對(duì)立木材積進(jìn)行擬合

以直徑為自變量，對(duì)立木材積進(jìn)行擬合，擬合結(jié)果見圖3-圖4。

最佳生長(zhǎng)曲線的確定決定系數(shù)R2為0.999，標(biāo)準(zhǔn)誤差為0.009，F(xiàn)值為5 810.44，P值為0.000。獲得以直徑為自變量的材積生長(zhǎng)模型表達(dá)式：

Y=-0.009+0.009x-0.001x2+6.47×10-5x3

式中，Y表示材積，X表示直徑。

圖3 以直徑為自變量獲得的12種曲線擬合圖

圖4 以直徑為自變量獲得的擬合度最高的材積生長(zhǎng)曲線圖

2.3 以林齡為自變量對(duì)立木材積進(jìn)行擬合

擬合結(jié)果見圖5-圖6。

圖5 以齡階為自變量獲得的12種曲線擬合圖

圖6 以齡階為自變量獲得的擬合度最高的材積生長(zhǎng)曲線圖

最佳生長(zhǎng)曲線的確定決定系數(shù)R2為0.997，標(biāo)準(zhǔn)誤差為0.013，F(xiàn)值為2 924.73，P值為0.000。獲得以林齡為自變量的材積生長(zhǎng)模型表達(dá)式：

Y=0.053-0.027x+0.003x2-4.493×10﹣5x3

式中，Y表示材積，X表示林齡。

根據(jù) Cubic(立方)曲線擬合結(jié)果，我們建立了東風(fēng)崗1號(hào)楊材積擬合曲線生長(zhǎng)模型(表11)，結(jié)果表明，在α=0.001時(shí)，三個(gè)方程均有統(tǒng)計(jì)學(xué)意義，其中，以直徑擬合的曲線方程決定系數(shù)和F值均最大，標(biāo)準(zhǔn)誤最小，以林齡擬合的次之，以樹高擬合的誤差相對(duì)較大。

表11 樹高、直徑和齡階與材積的擬合曲線模型

為驗(yàn)證三個(gè)方程的實(shí)用性，我們將樹干解析數(shù)據(jù)中的樹高、直徑和林齡分別代入相應(yīng)方程，計(jì)算結(jié)果表明，林齡在10 a前擬合的材積生長(zhǎng)量數(shù)值與實(shí)測(cè)數(shù)值誤差較大，10年生后的數(shù)值與實(shí)測(cè)值誤差較小，可以應(yīng)用。在生產(chǎn)中根據(jù)要求的精度選擇應(yīng)用一種方程進(jìn)行材積估算與預(yù)測(cè)，為林業(yè)生產(chǎn)提供科學(xué)參考。

3 小結(jié)

不同樹木的樹高和胸徑生長(zhǎng)都有一定的規(guī)律，因樹種和品種而不同，我們以樹高、直徑和林齡為自變量，對(duì)其材積生長(zhǎng)曲線進(jìn)行擬合，選擇最佳方案進(jìn)行材積生長(zhǎng)量預(yù)測(cè)，為立木材積計(jì)算提供了重要參考。