楠木二元立木材积表试编
卢 昌 泰
摘要 采用样木解析法收集和整理资料,探讨楠木材积与直径、树高之间的相关关系,建立楠木最佳二元立木材积方程,并编制楠木二元立木材积表,为楠木林分蓄积量测定等提供依据。关键词 样木解析法 楠木 二元材积方程 材积表
立木材积表是森林经营和森林调查中不可缺少的常用数表之一。目前我国森林调查常用数表的建设工作正走向正规化、系统化和科学化,但对于珍稀树种和国家重点保护树种的编表工作尚存在差距,应引起足够重视。因为这类树种需估计材积准确,从而为其保护和发展提供可靠依据。 楠木(Phoebe Zhennan S.Lee et F.N.Wei),系樟科楠木属树种,材质优良,为我国珍贵用材树种,也是国家二级重点保护树种。我省楠木林分蓄积量测定在目前仍沿用《四川和滇西北阔叶树二元立木材积表》,测定误差很大。为解决楠木林分调查用表问题,笔者采用样本解析法(1)收集编表资料,编制了林场《楠木一元材积表》(2)(1994),本文在此基础上,建立了林场楠木二元立木材积方程,编制了《楠木二元立木材积表》。
1 资料收集与数据预处理
1.1 资料收集
详见文献(2)。
1.2 数据预处理
(1) 计算解析木生长过程,取得各龄阶树高(H),去皮直径(D去)、去皮材积(V去)及当年带皮材积(V带)。 (2) 利用解析木各圆盘当年带皮直径(D带)与去皮直径(D去)求算其树皮系数(KB),KB=(D带)/D去,在计算机上用回归组程序建立D去,-KB最优相关方程。 KB=1.1459542-0.0150898D去+0.0004504D2去 (1) (R=0.9176004,U=0.6584922,S=0.1235740) (3) 利用解析木当年带皮材积(V带)、去皮材积(V去),在计算机上用回归组程序建立V去-V带最优相关方程。 V带=0.0005735+1.0708479V去-0.0381490V2去 (2) (R=0.9989950,U=1.0856784,S=0.0021854) (4) 将各解析木各龄阶D去代入(1)式求相应KB,再根据D带=KBD去求相应的D带;将V去代入(2)式求相应的V带。 这样就从对一株样木进行树干解析,经数据预处理得到了龄阶数个编表样本,即D带-H-V带对应值。 利用95株解析木共整理得编表样本1412个,样本量按径阶、树高分配情况见表1。
表1 样本按径阶、树高阶分配表
2.1 二元材积方程的建立
将1412个样本资料采用含5种二元数学模型的回归程序,在计算机上进行回归参数计算与回归显著性检验。5种数学模型如下: V=a(D2H)b (3) V=aDbHc (4) V=a+bH+cD2 (5) V=a+bH+cHD2 (6) V=a+bD2+cHD2 (7) 回归模型拟合及回归显著性检验结果见表2。
表2 回归参数及回归显著性检验
2.2 二元立木材积表的编制
将选定的最佳二元材积方程(即8式)编制程序输入计算器,代入相应的径阶(D)、树高(H)即得相应材积(V)。按径阶、树高阶整理编制成楠木二元立木材积表(表略)。
3 二元立木材积表的检验
为检验所编楠木二元立木材积表是否适用,利用林场楠木林分中收集的95株伐倒木材料,计算实测材积与回归估计值之间的系统误差,并进行F检验。
3.1 适用精度
根据收集的检验样本,分别径阶进行检验。经计算样本实测材积与回归估计值之间的系统误差为-2.57%,在《林业专业调查主要技术规定》要求范围内(±3%)(3)。
3.2 F检验
配制实测材积(Yi)与查表材积(Xi)的回归直线方程Yi=a+bXi。通过计算统计,结果为: Yi=0.0067687+1.0182804Xi (9) (R=0.9861005) 根据计算统计得到的采用下式进行F值的计算,得F=3.55。 查《F检验的临界值(Fa=0.05)表》得F0.05(1,93)=3.94,由于F<F0.05(1,93),说明实测材积Yi与查表材积Xi间无显著差异。 据以上两种检验结果,本研究所编《楠木二元立木材积表》能在林场使用。
4 结语
(1) 研究结果表明,林场楠木二元立木材积方程式V=0.000074954D1.884652825H0.881513308,回归关系显著,为最佳方程。 (2) 本研究编制的楠木二元立木材积表,经检验可在林场楠木林分蓄积量测定等方面使用。根据林昌庚等关于同一树种在不同地区的干形一般差别不大的观点(4),本表也可供其它地区参考使用,但使用前需进行适用性检验。 (3) 采用样本解析法收集编表资料,具有伐木量少、工作量省且能有效保护珍稀树种的优点,值得推广。