免费文献传递   相关文献

湖北桂花林场檫木次生林单木生长模型的研究



全 文 :Vol. 33 No. 4
Apr. 2013
第 33 卷 第 4期
2013年 4月
中 南 林 业 科 技 大 学 学 报
Journal of Central South University of Forestry & Technology
收稿日期:2012-12-03
基金项目:国家林业局行业公益性项目“南方集体林区次生林抚育间伐与高效利用技术研究”(201004032);湖南省“十一五”重
点学科建设计划资助项目“森林经理学科”(2006-028-0015)
作者简介:胡焕香(1985-),男,河南信阳人,硕士研究生,主要从事林业信息工程和森林可持续经营方 面的研究
通讯作者:佘济云(1966-),男,湖南邵东人,博士,教授,博导,主要从事森林经理和林业资源管理方面的教学和科研工作;
E-mail:shejiyun@126.com
檫木 Sassafras tzumu,别名檫树,为樟科檫木
属落叶乔木,主要分布在长江以南地区,广泛作为
用材林和风景林,其木材浅黄色,材质优良,细
致,耐久,多用于水车、造船及上等家具;檫木春
开黄花,且先花后叶,叶形奇特,秋季变红,花、
叶均具有较高的观赏价值,多用于庭园、公园栽植
或用作行道树、山区造林绿化,是行道绿化或城郊
风景林的理想树种 [1-2]。此外,檫木的根和树皮入
药,具有活血散瘀、祛风去湿等功能。因此,系统、
全面地研究檫木对于更加科学合理地开发、利用和
保护这一树种具有重要意义。目前关于檫木的研究
已有很多,但大多是关于生物量 [3-5]、林分结构 [6-7]、
育苗与栽培 [8-9]等方面的研究,对其单木生长模型 [10]
的研究还比较罕见。直径 [11-12]是林分内部最重要、
最基本的测树因子,其测量简单、方便而且精确,
因此,本研究拟对湖北国营桂花林场的檫木次生林
直径生长过程进行分析与探讨,以期为研究区檫木
的经营管理提供一定的参考依据。
1 材料与方法
1.1 试验地概况
实验地点位于湖北省国营桂花林场,林场建
湖北桂花林场檫木次生林单木生长模型的研究
胡焕香 1,佘济云 1,李 俊 2,张 敏 1,孟 伟 1,程玉娜 1
(1中南林业科技大学,湖南 长沙 410004;2广东省林业调查规划院,广东 广州 510520)
摘 要:为了更好地开发和利用檫木资源并为其经营管理提供参考依据,利用大小比数并结合多模型选优法对
湖北省国营桂花林场檫木次生林的单木直径生长过程进行了分析与探讨 . 结果表明:不论是竞争优势木还是竞争
劣势木,用 Richards函数来拟合檫木的直径生长过程效果最好,其次是 Schumacher函数,最差的是 Gauss函数。
但是 Richards函数的参数偏多,并且通过 Richards函数预估的直径生长量与檫木的实际生长状况不符合,因此,
Schumacher函数是描述湖北桂花林场檫木次生林直径生长过程的最佳模型。
关键词:檫木次生林;生长模型;大小比数;直径
中图分类号:S758.1 文献标志码:A 文章编号:1673-923X(2013)04-0061-05
Study on individual tree growth model of Sassafras tsumu secondary
forest on Osmanthus Forest Farm in Hubei Province
HU Huan-xiang1, SHE Ji-yun1, LI Jun2, ZHANG Min1, MENG Wei1, CHENG Yu-na1
(1.Central South University of Forestry and Technology, Changsha 410004, Hunan, China; 2.Forestry Surveying and Designing Institute
of Guangdong Province, Guangzhou 510520, Guangdong, China)
Abstract: In order to develop and use Sassafras albidum better and provide reference basis for its operation and management, the neighborhood
comparison combined with multi-model selection methods were used to analyze and discuss individual tree diameter growth progress of S.
tsumu secondary forest in Osmanthus forest farm in Hubei province. The results show that whether it is a dominant tree or oppressed tree,
Richards function was the best to fi t the growth process of S. albidum diameter, better than Schumacher function, while Gauss function’s fi tting
effect was the worst. However, the parameters of Richard s function were too much, and the amount of diameter growth of S. albidum estimated
by Richards function was not incompatible with the actual growth conditions. Therefore, Schumacher function was the best model to describe
the diameter growth process of S. tsumu secondary forest in Osmanthus forest farm in Hubei province.
Key words: Sassafras tsumu secondary forest; growth model; neighborhood comparison; diameter; Hubei Province
胡焕香,等:湖北桂花林场檫木次生林单木生长模型的研究62 第 4期
于 1964年,位于湘鄂赣边陲的湖北省崇阳县桂
花泉镇,距离崇阳县城大约 16 km。林场地处幕
阜山北坡低山丘陵与江汉平原的过渡地带,属低
山丘陵地貌 [13]。整个林场森林覆盖面积约有 1万
hm2,活立木蓄积42万m3,是全国“十佳国有林场”,
是鄂东南的森林氧吧,是崇阳县的一道绿色屏障。
桂花林场属亚热带季风气候,气候温暖,四季分明,
雨量充沛、湿热同步,年平均降雨量 1 636.2 mm,
年平均气温 15.5℃。林场内土层深厚,土壤肥沃,
十分适宜众多植物的生长与繁衍,生物多样性丰
富。山上生长着大量的人工次生林,较常见的乔
木树种有檫木、南酸枣 Choerospondias axillaria、
苦槠 Castanopsis sclerophylla、杉木 Cunninghamia
lanceolata、 马 尾 松 Pinus massoniana、 枫 香
Liquidambar formosana 等。
1.2 试验材料
2011年 5月与 7月,在湖北省国营桂花林场
中选择具有代表性的檫木次生林林分类型,用罗
盘仪分别设置 6块面积大小为 400 m2,规格为 20
m×20 m的正方形标准地,对其中胸径≥ 5 cm的
檫木进行每木检尺,用围尺测量胸径,用测高仪
测量树高,并做好记录。依据标准地每木检尺数据,
寻找 2~ 3株与林分平均直径和平均高相接近(一
般要求相差在±5%以下)且干形中等的林木作为
平均标准木,每块标准地选取 2棵标准木做树干
解析。样地基本信息见表 1。
表 1 样地基本信息
Table 1 Basic information of the sample plots
样地
编号
地形地貌 主要测树因子
海拔
/m 坡向 坡位
坡度
/(°)
平均胸
径 /cm
平均树
高 /m
平均冠
幅 /m
平均年
龄 /a
1 203.8 东南 中 24 9.2 7.8 2.2 15
2 390.3 南 下 16 11.8 10.9 2.6 21
3 387.8 南 下 16 10.2 7.0 2.4 18
4 292.4 北 中 21 15.0 9.9 2.2 24
5 291.6 北 中 19 15.8 10.1 2.1 25
6 292.9 南 中 9 16.4 11.8 3.1 25
1.3 研究方法
1.3.1 林木竞争态势的研究方法
本研究利用大小比数 [6,13]对檫木的竞争态势
进行分类,从而科学、合理地表达现实林木的生
长状态。大小比数是用来描述相邻木与参照木之
间在胸径、树高、冠幅等某一指标上的个体优势
程度,本研究定义为胸径大于参照木的相邻木占
4株最近相邻木的株数比例,其计算公式如下:



4
14
1
i
iji kU (1)
式中:Ui表示大小比数; kij表示一个离散型
的变量,当相邻木的胸径比参照木大时 kij=1,反
之 kij=0。
1.3.2 单木生长模型的选择
以林分中各单株林木为基本单位,以与其相
邻木之间的竞争关系为基础,模拟林分中单株树
木生长过程的模型,称为单木生长模型 [14]。用来
描述单木生长的模型有很多,本研究根据实际情
况来选取以下 5个逻辑性强,适应性广的数学模
型对湖北省国营桂花林场的檫木次生林的直径生
长过程进行研究与分析。并利用 DPS软件来处理
与分析数据。
(1)Logistic函数:

)(exp1 bta
ky  。 (2)
(2)Gauss函数:
y=a[1-bexp(-kt2)]。 (3)
(3)Richards函数:

cktb
ay 1
])(-exp1[ 
 。 (4)
(4)Gompertz函数:
y=aexp[-bexp(-ct)]。 (5)
(5)Schumacher函数:
)(-exp
tc
bay  。 (6)
式中:y表示直径生长量;k、a、b、c为随机
参数;t表示树木的年龄。
1.3.3 模型的求解与检验
本研究利用麦夸特法来求取各模型的参数,
并通过方差分析来对模型的拟合效果进行显著性
检验。在拟合过程中,可根据数据拟合图中实际
观测值和拟合曲线的拟合情况、残差平方和、决
定系数及显著水平的高低来选择最优的模型。残
差平方和越小,决定系数越大,显著水平越低,
参数越少,则拟合的效果越理想。
2 结果与分析
2.1 标准木竞争态势的划分
由于不同林分的生长条件不同,其竞争状态
也存在差异,因此对同树种不同竞争态势的林木
63第 33卷 中 南 林 业 科 技 大 学 学 报
建立相应的生长模型是更加科学、合理的。依据
大小比数可以将竞争态势分为两类,当大小比数
≥ 0.5时,为竞争优势木,反之为竞争劣势木。经
统计分析,湖北省国营桂花林场檫木次生林 12株
标准木的大小比数分布情况如图 1所示。由图 1
可知,竞争优势木(大小比数≥ 0.5)与竞争劣势
木(大小比数< 0.5)的标准木株树均为 6株,说
明檫木次生林的竞争态势比较均衡。
2.2 单木生长模型的拟合与检验
由单木生长模型的拟合与检验结果见表 2。可
知,不论是竞争优势木还是竞争劣势木,5种数学
模型对檫木次生林生长过程的拟合效果都比较接
近,而且效果都非常理想,直径生长量与树木年
龄回归关系均表现为极其显著。就单木(无论是
竞争优势木还是竞争劣势木)而言,残差平方和
最小,决定系数最大,对现有标准木数据拟合效
果最好的都是 Richards函数,其次是 Schumacher
函数,Gauss函数的拟合效果最差。残差平方和的
大小顺序为:Richards函数< Schumacher函数<
Gompertz函数< Logistic函数< Gauss函数;决定
系数大小排序均为:Richards函数> Schumacher
函数> Gompertz函数> Logistic函数> Gauss函
数。说明 Richards函数用来拟合湖北省国营桂花
林场檫木次生林的单木生长过程效果最佳,其次
是 Schumacher函数,Gauss函数的拟合效果最差。
表 2 单木生长模型的拟合与检验结果
Table 2 Fitting and testing results of individual tree growth model
竞争态势 模型名称
方差分析
决定系数
参数
F值 显著水平 残差平方和 k a b c
优势木
Logistic函数 153.074 3 0.000 0 380.476 7 0.843 0 25.155 3 2.175 4 -0.115 4 ——
Gauss函数 151.843 3 0.000 0 383.073 6 0.842 0 0.001 7 24.279 8 0.871 3 ——
Richards函数 105.577 7 0.000 0 364.189 5 0.849 8 0.010 9 74.983 9 -0.995 6 -0.939 9
Gompertz函数 157.928 0 0.000 0 370.580 8 0.847 1 —— 30.053 2 2.802 3 0.062 7
Schumacher函数 160.152 1 0.000 0 366.201 9 0.848 9 —— 65.754 4 52.382 7 13.125 1
劣势木
Logistic函数 285.962 0 0.000 0 70.885 9 0.899 4 14.598 6 2.120 6 -0.146 8 ——
Gauss函数 279.491 6 0.000 0 72.358 4 0.897 3 0.002 8 14.093 0 0.862 5 ——
Richards函数 213.393 1 0.000 0 63.104 4 0.910 4 0.005 9 82.550 1 -0.997 2 -1.021 5
Gompertz函数 308.879 6 0.000 0 66.120 2 0.906 1 —— 17.687 0 2.710 1 0.078 0
Schumacher函数 320.402 2 0.000 0 63.958 3 0.909 2 —— 39.030 9 42.671 4 11.152 1
2.3 直径生长量预估分析
由直径生长量预估分析结果(见图 2和图 3)
可知,从整体而言,无论是竞争优势木还是竞争
劣势木,不论是直径总生长量还是直径连年生长
量,5种数学模型对现有标准木数据的拟合效果都
比较接近,并且重合度很高。在总生长量方面,
竞争优势木基本上都是从 40年左右开始明显分
化,而竞争劣势木则是在 30年左右开始明显分化,
Richards函数预估的直径总生长量最大,Gauss
函数预估的直径总生长量最小,其大小排序为:
Richards函数> Schumacher函数> Gompertz函数
> Logistic函数> Gauss函数。
在连年生长量方面,竞争优势木在 5年左右
出现第一个明显高峰,随后出现一个 5年左右时
间的缓慢生长期,之后生长速率一路上升至 20年
左右到达连年生长第二高峰,随后又开始进入缓
慢生长期直至衰老;竞争劣势木在 5年左右出现
一个第 1个明显的高峰,随后出现一个 5年左右
时间的短时期缓慢生长期,之后生长速率一路上
升至 15年左右到达连年生长第 2个高峰,随后又
开始进入缓慢生长期直至衰老。Gauss函数的收
敛速度最快,Richards函数的收敛速度最慢,5中
数学模型收敛速度的快慢排序为:Gauss函数>
Logistic函数> Gompertz函数> Schumacher函数
> Richards函数。
就成熟期和生长停止期而言,在竞争优势木
方面:Logistic 函数、Gauss 函数和 Gompertz 函
数 3种模型的成熟期和生长停止期都出现的较
图 1 檫木次生林 12株标准木的大小比数分布
Fig.1 Neighborhood comparison of 12 standard wood in
S. tsumu secondary forest
胡焕香,等:湖北桂花林场檫木次生林单木生长模型的研究64 第 4期
早,在 20年左右出现成熟期,在 55年左右出现
生长停止期,与檫木次生林生长的现实状况不符
合。Richards函数在 5年左右到达高峰后一路缓
慢走低,因此无法判断成熟期,5~ 100年一直为
缓慢生长期,100年左右也未出现生长停止期;
Schumacher函数的成熟期出现的比较适中,20年
左右出现成熟期,而在 20~ 100年出现缓慢生
长期,与檫木次生林的现实生长状况比较符合。
在竞争劣势木方面:Logistic函数、Gauss函数和
Gompertz函数 3种模型的成熟期和生长停止期都
出现的较早,在 15年左右出现成熟期,在 50年
左右出现生长停止期,不符合檫木次生林生长的
实际状况。Richards函数在 5年左右到达高峰后一
路缓慢走低,因此无法判断成熟期,5~ 100年一
直为缓慢生长期,100年左右也未出现生长停止期;
Schumacher函数的成熟期出现的适中,15年左右
出现成熟期,而在 15~ 100年出现缓慢生长期,
较符合檫木次生林的实际生长状况。
3 结论与讨论
本研究利用 5个生长模型对湖北省国营桂
花林场的檫木次生林的单木直径生长过程进行了
拟合,得出对现有标准木数据拟合效果最好的
是 Richards函数,其次是 Schumacher函数,而
Gauss函数的拟合效果最差。但是考虑到 Richards
函数的参数较多(有 4个参数),而拟合效果与
其相近的 Schumacher函数的参数适当偏少(有
3个参数),从减少模型复杂程度方面考虑,选择
Schumache函数作为檫木次生林的单木直径生长模型
最为合适。此外,根据直径生长量预估分析的结果
可知,无论是竞争优势木还是竞争劣势木,Richards
函数都无法判断檫木直径生长过程的成熟期,而通
过 Schumache函数来判断成熟期得出的结果最为理
想,与檫木次生林的实际生长情况十分吻合。因
此,不论是竞争优势木还是竞争劣势木,Schumache
函数均是描述檫木次生林直径生长过程的最理想模
型,竞争优势木与竞争劣势木的直径生长模型表达
式分别为:y=65.754 4exp[-52.382 7/ (13.125 1+t)];
y=39.030 9exp[-42.671 4/(11.152 1+t)]。
檫木是湖北省国营桂花林场的主要树种之一,
研究结果可以用来直接判定檫木各单株木的直径
生长状况和生长潜力,这对于研究区檫木林分的
集约经营、生长与收获预估等方面具有一定的参
考价值。但是由于条件限制,本研究所选取的样
地和标准木数量偏少,这对研究结果会造成一定
的影响。如果今后的研究能对这些方面加以补充
和完善,并增加对树高、冠幅、断面积等因子进
行单木生长模型探讨,研究成果将会更加具有应
图 2 竞争优势木直径生长量预估
Fig. 2 Forecast amount of diameter growth of dominant tree
图 3 竞争劣势木直径生长量预估
Fig.3 Forecast amount of diameter growth of oppressed tree
65第 33卷 中 南 林 业 科 技 大 学 学 报
用与参考价值。
参考文献:
[1] 田大伦 ,张昌剑 ,罗中甫 ,等 .天然擦木混交林的生物量及营
养元素分布Ⅰ .生物生产量及生产力 [J].中南林学院学报 ,
1990, 10(2):121-127.
[2] 范海兰 ,洪 滔 ,吴承祯 .檫木种群基面积增长规律研究 [J].
江西农业大学学报 ,2005,27(5):744-747.
[3] 沈 燕,田大伦,闫文德,等 .湖南沅陵天然次生白栎+檫
木+杉木群落生物量及其分布规律 [J].中南林业科技大学学
报,2011,31 (5):44-50.
[4] 沈 燕,田大伦,项文化,等 .天然次生檫木枫香混交林生
物量及生产力研究 [J].中南林业科技大学学报,2011,31(5):
26-30.
[5] 张 琼,洪 伟,吴承祯,等 .天然更新檫木林的能量分析
[J].植物资源与环境学报,2004,13(2):36-39.
[6] 曾思齐,李东丽,宋武刚,等 .檫木次生林空间结构的研究
[J].中南林业科技大学学报,2012,32 (3):1-6.
[7] 李 俊 .南方集体林区典型林分类型结构特征及生长模型研
究 [D].长沙:中南林业科技大学,2012.
[8] 黄宝祥,朱培林,符树根 .檫木的组织培养 [J].江西林业科
技 , 2010,(4):11-12.
[9] 丁小飞,陈红林,曹 健 .檫木三个群体的遗传结构初探 [J].
湖北林业科技 ,2006,(5):1-2.
[10] 刘 平 ,王玉涛 ,杨 帆 .基于单木生长模型的森林动态模拟
系统研究进展 [J].世界林业研究 ,2011,24(5):25-29.
[11] 李 俊,佘济云,胡焕香,等 .昌化江流域天然林直径结构
研究 [J].中南林业科技大学学报,2012,32 (3):37-43.
[12] 胡焕香,佘济云,刘立武,等 .海南阿陀岭 4种加勒比松林
分直径分布规律研究 [J].中国农学通报 ,2012,28(22):16-21.
[13] 曾思齐,李 俊,李东丽,等 .南方集体林区南酸枣次生林
林分结构研究 [J].中南林业科技大学学报,2012,32 (4):1-6.
[14] 孟宪宇 .测树学 [M].北京:中国林业出版社,2008:254-293.
[本文编校:吴 彬 ]