伊春林区天然毛榛树高地径生长量调查研究

姜伟阳，霍冬磊，李勇，万道印

(黑龙江省带岭林业科学研究所，黑龙江 带岭 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