

姓名：景林木提问时间： 2007/5/14 21:48:30	邮箱：
问题：

请教吕专家,我想编制我县桉树二元材积表,但不知a ,b,c三个参数怎么计算,请吕专家给予例子说明,谢谢!

回复专家：张连翔，回复时间：2007/5/14 22:17:58

建议查阅文献并参考之：

漳州市桉树人工林(尾叶桉、巨尾桉、尾巨桉)二元材积表的编制

<<林业调查规划 >>2006年04期
黄建海

在福建省漳州市11个林场设置70块样地并实测350株桉树人工林样木的胸径、树高和材积,应用改进单纯形法拟合22个二元材积方程,并筛选最优方程编制桉树人工林二元材积表,结果表明:二元材积的通用模型V=a0(a1+a2d)a(a4+a5h)a6,拟合效果优于其他方程,其相关系数R2=0.993 8,残差平方和Q=0.003 363,利用该方程编制的桉树二元材积表能满足林业生产上的精度要求.

回复专家：吕玉奎，回复时间：2007/5/15 16:03:19

楠木二元立木材积表试编

卢昌泰

摘要　采用样木解析法收集和整理资料，探讨楠木材积与直径、树高之间的相关关系，建立楠
木最佳二元立木材积方程，并编制楠木二元立木材积表，为楠木林分蓄积量测定等提供依据。
关键词　样木解析法　楠木　二元材积方程　材积表

　　立木材积表是森林经营和森林调查中不可缺少的常用数表之一。目前我国森林调查常用数
表的建设工作正走向正规化、系统化和科学化，但对于珍稀树种和国家重点保护树种的编表工
作尚存在差距，应引起足够重视。因为这类树种需估计材积准确，从而为其保护和发展提供可
靠依据。
　　楠木(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.0004504D²_去　　　　　(1)
　　(R=0.9176004,U=0.6584922,S=0.1235740)
　　(3) 利用解析木当年带皮材积(V_带)、去皮材积(V_去)，在计算机上用回归组程序建立V_去-V_带
最优相关方程。
　　V_带=0.0005735+1.0708479V_去-0.0381490V²_去　　　　(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	4	6	8	10	12	14	16	18	20	22	24	合　计
2	34	3											37
3	40	17											57
4	41	42	6										89
5	11	32	16	1									60
6	2	35	32	6									75
7		15	32	12	2								61
8		8	35	26	9	1							79
9			15	16	14	6	2						53
10		2	16	27	26	14	1	1					87
11			6	19	17	13	7	2	2				66
12			9	28	27	26	9	3	2	1			105
13				12	14	22	13	9		2			72
14			1	15	25	26	24	14	4	1	3		113
15				6	19	20	18	21	5	4	3		96
16				2	18	12	16	25	14	6			93
17					10	9	18	20	20	5	5		87
18					6	8	9	17	21	14	3	2	80
19					3	7	1	13	16	11	9	6	66
20						1		2	4	8	6	2	23
21							1		1	4	3		9
22										1	1	2	4
23
合　计	128	154	168	170	190	165	119	127	89	57	33	12	1412

2　二元立木材积表的编制

2.1　二元材积方程的建立

　　将1412个样本资料采用含5种二元数学模型的回归程序，在计算机上进行回归参数计算与回
归显著性检验。5种数学模型如下：
　　V=a(D²H)^b　　　　　　　　(3)
　　V=aD^bH^c　　　　　　　　　(4)
　　V=a+bH+cD²　　　　　　　　(5)
　　V=a+bH+cHD²　　　　　　　(6)
　　V=a+bD²+cHD²　　　　　　　(7)
　　回归模型拟合及回归显著性检验结果见表2。

表2　回归参数及回归显著性检验

数学模型	回　归　参　数				回　硅　显　著　性　检　验
数学模型	a	b	c	r	U	Q	S	F	F_(0.05)^(m,n-m-1)
(3)	0.000073368	0.924435658		0.991565844	4449.585388	76.01734739	0.23227400	41237.07303	3.84
(4)	0.000074954	1.884652825	0.881513308	0.991581972	4449.730135	75.87260095	0.22205311	41317.08766	3.00
(5)	0.048682247	0.009481822	0.000180672	0.81387383	7.992597878	4.073693105	0.001432121	346.858799	3.00
(6)	-0.024110881	0.006030802	0.000018937	0.862293720	8.962293720	3.103997263	0.001237542	2043.131030	3.00
(7)	0.031781842	-0.000571795	0.000062852	9.852743099	9.852743099	2.213547884	0.003963615	3135.806995	3.00

由表2可以看出，楠木最佳二元材积方程为：
V=0.000074954D^1.884652825H^0.881513308　　　（8）
R=0.991581972 U=4449.730135 Q=75.87260095
S=0.22205311 F=41317.08766 F_0.05(2,1409)=3.00
　　回归关系是高度显著的。

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。
g5503.gif (2180 bytes)
　　查《F检验的临界值(Fa=0.05)表》得F_0.05(1，93)=3.94，由于F＜F_0.05(1，93)，说明实
测材积Yi与查表材积Xi间无显著差异。
　　据以上两种检验结果，本研究所编《楠木二元立木材积表》能在林场使用。

4　结语

　　(1) 研究结果表明，林场楠木二元立木材积方程式
V=0.000074954D^1.884652825H^0.881513308，回归关系显著，为最佳方程。
　　(2) 本研究编制的楠木二元立木材积表，经检验可在林场楠木林分蓄积量测定等方面使
用。根据林昌庚等关于同一树种在不同地区的干形一般差别不大的观点(4)，本表也可供其它地
区参考使用，但使用前需进行适用性检验。
　　(3) 采用样本解析法收集编表资料，具有伐木量少、工作量省且能有效保护珍稀树种的优
点，值得推广。