伊春林區天然毛榛樹高地徑生長量調查研究

姜偉陽，霍冬磊，李勇，萬道印

(黑龍江省帶嶺林業科學研究所，黑龍江 帶嶺 153106)

伊春林區天然毛榛樹高地徑生長量調查研究

姜偉陽，霍冬磊，李勇，萬道印*

(黑龍江省帶嶺林業科學研究所，黑龍江 帶嶺 153106)

通過對伊春林區天然毛榛生長情況進行調查，擬合了毛榛生物量預測方程，探討天然毛榛林合理經營周期。研究結果表明：天然毛榛榛樹樹齡(A)與地徑(D)、樹高(H)相關程度較高，冪函數模型是地徑(D)、樹高(H)生長最佳模型，天然毛榛林合理經營周期為10年。

毛榛；生長量；經營周期

毛榛(CorylusmandshuricaMaxim.)為樺木科榛屬落葉灌木或小喬木，榛樹木材堅硬、耐腐。近幾年，隨著封山育林的進行，袋裝食用菌培養料來源受到限制，伊春林區經過試驗，“架條木耳” 培育取得成功，榛樹椴木栽培黑木耳技術得到廣泛應用，并獲得較好的經濟效益，拓寬了食用菌培養料來源。但椴木栽培黑木耳要求榛樹樹干胸徑達到2cm以上，小于2cm的樹干和樹枝不能被利用。為充分利用毛榛木質資源，本研究通過對伊春林區天然毛榛生長情況進行調查，擬合了毛榛生物量預測方程[1-2]，探討天然毛榛林合理經營周期，為伊春林區毛榛林開發利用提供參考依據。

1 研究方法

1.1 樣木的采集及生長量測定

表1 天然榛樹生長情況調查統計表

為便于測定天然毛榛樹高和地徑生長量，在2013年秋季落葉后對榛樹生長情況及生物量進行了調查。在帶嶺林業局和翠巒林業局選擇4塊不同立地條件集中連片天然榛林，其中帶嶺林業局2塊榛林是經過多次反復割采形成的萌生林，翠巒林業局2塊榛林未見割采痕跡，無法確定林分起源。在所選林分中隨機選取一定數量的榛樹，從貼近地面部位伐倒后，查數年輪確定樹齡，測量地徑、樹高、冠幅。本次共采集榛樹樣木121株，調查數據統計結果見表1。

1.2 建模方法

榛樹地徑、樹高生長方程以樹齡(A)為自變量，以地徑D、樹高H、生物量W為因變量，選擇冪函數、對數、多項式、指數函數、Bertalanffy、Compertz、Schumacher、logistic經驗回歸模型，根據所測樣木的各調查因子數據，采用SAS9.1軟件包分別擬合出樹木地徑、樹高生長方程。

1.3 模型的評價指標

采用殘差平方和(SSE)、均方根誤差 (RMSE) 和確定系數 (R2)等擬合統計量作為比較和評價備選模型的標準。根據所計算的方程擬合統計量，選擇其中殘差平方和(SSE)、均方根誤差(RMSE)相對較小而確定系數(R2)相對較大的方程作為最優的回歸方程。

2 結果與分析

利用8種理論模型，擬合的榛樹地徑、樹高生長模型及統計量見表2～表5。

表2 地徑生長模型

表3 地徑生長模型擬合結果

表4 樹高生長模型

表5 樹高生長模型擬合結果

從表2～表5可知，用 8種理論生長模型分別擬合榛樹地徑(D)、樹高(H)生長過程效果很好。從統計指標R2看，樹齡(A)與地徑(D)、樹高(H)相關程度較高。根據殘差平方和 ( SSE) 、均方根誤差 ( RMSE) 較小比較好而確定系數 (R2) 相對較大比較好，冪函數模型是地徑(D)、樹高(H)生長最佳模型。

3 天然榛林合理采伐利用時間的確定

采用前面擬合的榛樹地徑、樹高生長模型，編制了天然榛樹生長過程表，繪制了榛樹生長曲線，見表6。

表6 天然毛榛榛樹生長過程表

從表6可以看出，榛樹幼齡時期地徑、樹高生長較快，這是與多數喬木樹種生長過程明顯的區別，可能與榛樹起源多是天然萌生有關。從第10年起地徑、樹高連年生長量開始放緩，平均地徑達到26.47mm，已經能夠滿足食用菌生產的需要，工藝成熟。因此，建議天然毛榛林合理經營周期為10年。

4 小結

4.1 天然毛榛榛樹樹齡(A)與地徑(D)、樹高(H)相關程度較高，冪函數模型是地徑(D)、樹高(H)生長最佳模型。

4.2 由于榛樹生命周期短，地徑每年生長量小，樹干解析困難，本次使用同一樹齡多株樣木的平均數來擬合榛樹一般生長過程，沒區分林分起源、經營密度和立地條件，模型的適用性還存在一定的局限性，有待于今后在實踐中不斷完善。

4.3 天然毛榛林合理經營周期為10年。

[1]王曉林，郭斌.柞樹樹高與胸徑相關關系的研究[J].森林工程，2012,28(6)：18-21.

[2]馬武，雷相東.蒙古櫟天然林單木生長模型的研究[J].西北農林科技大學學報，2015,43(3)：83-90.

2017-09-30

黑龍江省財政資助項目(2011-01)

姜偉陽(1988-)，女，助理工程師，主要從事林業生態方向研究；*通訊作者：萬道印(1986-)，男，高級經濟師，主要從事林業資源開發與利用方向科研工作，E-mail：mrs_jiangjiang@163.com。

DOI.:10.13268/j.cnki.fbsic.2017.06.015

S664.4

A