全 文 :第 33 卷第 1 期
2014 年 2 月
中南林业调查规划
CENTRAL SOUTH FOREST INVENTORY AND PLANNING
Vol. 33 No. 1
Feb. 2014
* 收稿日期: 2013 —11 —26
作者简介: 普布顿珠 ( 1973 — ) ,男,西藏昌都类乌齐人,工程师,从事森林资源保护与培育的相关研究。
通讯作者: 边巴多吉,副研究员。
西藏高原林芝云杉树高与胸径异速生长方程建模
普布顿珠1,边巴多吉1,郭其强2,3,朱雪林1,朱 莉4
( 1. 西藏自治区林业调查规划研究院,拉萨 850000; 2. 西藏高原生态研究所,西藏 林芝 860000;
3. 西藏林芝高山森林生态系统国家野外科学观测研究站,西藏 林芝 860000;
4. 西藏拉萨师范高等专科学校,拉萨 850007)
摘 要:以 150 株西藏高原林芝云杉解析木资料为样本,利用 Excel 2003 软件中的 5 种数学模型对林芝
云杉的树高与胸径异速生长的进行建模,并分析其精度;结果显示幂指数方程 H = 1. 68D0. 777拟合效果
最佳( R2 = 0. 981) ;该模型可用于西藏高原东南部快速准确地测算林芝云杉平均树高,也可为计算立木
材积、编制出材率表和评价林分的立地质量提供基础数据,对于了解林芝云杉生长状况和指导当地林业
生产具有重要作用。
关键词:林芝云杉;树高;胸径;异速生长方程
中图分类号:S 758 文献标识码:A 文章编号:1003 —6075(2014)01 —0013 —03
Height-Diameter Allometric Growth Equations for Picea likiangensis
var. linzhiensis in Tibet Plateau
Pubu Dunzhu1,Bianba Duoji1,GUO Qiqiang2,3,ZHU Xuelin1,ZHU Li4
(1. Forest Inventory and Planning Institute of Tibet Autonomous Region,Lhasa 850000 ,Tibet,China;2. Research Institute of Tibet Plat-
eau Ecology,Linzhi 860000,Tibet,China;3. National Forest Ecosystem Research Station in the Sejila Mountain of Tibet,Linzhi 860000,
Tibet,China;4. Teachers College in Lhasa,Lhasa 850007 ,Tibet,China)
Abstract:In this paper,five height-diameter allometric growth equations were evaluated by making use of Ex-
cel 2003 software,based on the data of 150 parse trees of Picea likiangensis var. linzhiensis in Tibet plateau,
and every equation’s accuracy was analyzed,the results showed that:The exponential equation (H =
1. 68D0. 777)was the best fit (R2 = 0. 981),it can be used to calculate the average tree height of Picea likian-
gensis var. linzhiensis in southeastern Tibet plateau precisely and rapidly ,and can provide the basic data for
the calculation of tree volume,compilation of merchant volume table and the evaluation of forest-site quality,
which is very important for understanding the growing status of Picea likiangensis var. linzhiensis and guiding
forestry production.
Key words:Picea likiangensis var. linzhiensis;tree height;diameter at breast height(DBH) ;allometric growth
植物的异速生长主要是研究植物个体大小和其他
属性之间的非线性关系,通常以幂指数 (Y = axb)
形式表示。在森林资源调查中树木的胸径、树高是最
重要的两项因子,二者的相互关系称为树高曲线,也
属于植物异速生长关系的一种[1]。通过应用异速生长
方程研究树高与胸径关系,一方面可以求算林分各径
·31·
DOI:10.16166/j.cnki.cn43-1095.2014.01.009
2014 年第 1 期 中南 林 业 调 查 规 划 第 33 卷
阶相对应的平均树高,尤其是林分平均胸径相对应的
林分平均高;另一方面可用于计算树干材积、编制出
材率表、评价立地质量等[2 - 3],对于指导当地林业生
产具有重要的现实意义。
林芝云杉[4] (Picea likiangensis var. linzhiensis)
为松科云杉属常绿高大乔木,为西藏亚高山暗针叶林
的主要建群种之一,主要生长在藏东南的林芝地区和
昌都地区。林芝云杉不仅是青藏高原特有的速生用材
树种,而且在提供工业用材、维护生态环境和构建高
原生态屏障等方面发挥着重要的作用。学者们针对林
芝云杉开展了群落特征与林分生产力[4 - 5]、生态[6]水
文功能和光合特性[7]等方面的研究,而对该树种在高
寒生境的生长过程及状况了解甚少。鉴于此,本研究
拟利用西藏林芝云杉解析木数据拟合胸径与树高关系
的异速生长方程,以便为当地的林业生产提供基础
资料。
1 材料与方法
1. 1 研究地概况
研究区包括西藏东南部的林芝地区和昌都地区,
位于雅鲁藏布江中游河谷地带。属高原温带半湿润季
风气候,年平均气温 6. 5 ℃,最冷月 (1 月)平均气
温 0. 3 ℃,最热月 (7—8 月份)平均气温 15. 6 ℃,
无霜期 170 d,年降水量 600 ~ 900 mm,4—10 月份
降水量占全年降水量的 85%左右;土壤类型主要为
山地棕壤和暗棕壤;亚高山针叶林区,地带性植被以
高山松林 (Pinus densata)、高山栎林 (Quercus seme-
carpifolia)、林芝云杉和急尖长苞冷杉林 (Abies geor-
gei var. smithii)为主。
1. 2 试验材料
为满足建模精度要求,林芝云杉样本按 18 个径
阶段 (分别为 2,4,6,8,10,12,14,16,18,
20,22,24,26,32,38,40,44 和 50 cm),且每
个径阶组内的样本按个体高度均匀分布,并考虑地
形、坡位、年龄和个体生长状况等因子,共选取具有
代表性的个体共 150 株。在离地面 1. 3 m处用围尺测
量胸径,再用超声波测高测距仪 (Vertex IV,瑞典
产)测定树高。
1. 3 数据处理
将林芝云杉样本按 2 cm的径阶分组并统计株数,
计算各径阶的平均胸径和树高,详见表 1。
采用 Excel 2003 软件进行数据建模和统计分析,
用 Office 2003 绘制图形。
表 1 调查各径阶林芝云杉概况
径阶 /
cm
平均胸径 /
cm
平均树高 /
m
株数 /
株
占总株数
比例 /%
2 2. 21 2. 99 5 3. 33
4 4. 21 4. 44 7 4. 67
6 6. 07 5. 34 6 4. 00
8 7. 93 6. 97 9 6. 00
10 11. 20 11. 10 3 2. 00
12 12. 50 10. 23 7 4. 67
14 13. 80 6. 77 8 5. 33
16 16. 03 12. 81 10 6. 67
18 19. 17 16. 44 12 8. 00
20 20. 14 18. 52 8 5. 33
22 21. 94 19. 32 12 8. 00
24 24. 15 21. 32 13 8. 67
26 26. 22 19. 15 12 8. 00
32 31. 85 19. 25 15 10. 00
38 38. 27 22. 43 10 6. 67
40 40. 70 24. 10 5 3. 33
44 44. 37 27. 74 5 3. 33
50 51. 53 32. 41 3 2. 00
2 结果与分析
2. 1 径阶分析
林芝云杉按胸径 < 10,10 ~ 19,20 ~ 39,≥
40 cm的株数分别统计 (结果见图 1),胸径 < 10 cm
的占 20%,10 ~ 19 cm 的占 30%,20 ~ 39 cm 的占
41. 33%,≥40 cm 的占 8. 67%。根据前人对西藏林
芝云杉林生长状况的研究[5,8 - 9]结果表明,该树种胸
径在 20 ~ 39 cm的个体处于迅速生长期。说明本研究
中的样木选择基本以中龄期和成熟期的个体为主,同
时兼顾幼龄期和老龄期个体。
图 1 林芝云杉胸径分布比例
·41·
2014 年第 1 期 普布顿珠,等: 西藏高原林芝云杉树高与胸径异速生长方程建模 第 33 卷
2. 2 建模与检验
2. 2. 1 建模
根据表 1 中各径阶实测的树高平均值,以胸径为
横坐标,树高平均值为纵坐标,绘制径阶与树高分布
的散点图 (图 2)。通过观察散点形状与分布趋势,
初步选定可能符合散点分布趋势的几种相关模型,见
表 2。根据表 1 的平均胸径与对应平均树高成对数
据,利用 Excel2003 软件拟合所选几种相关模型,得
各模型参数及复相关系数 R2。表 2 中:H 是平均树
高的回归估计值;D 是各径阶的平均胸径;a,b,c
为参数。
图 2 林芝云杉树高与胸径散点图
表 2 林芝云杉树高与胸径的模型拟合
拟合方程
参 数
a b c
R2
H = aD + b 3. 386 0. 562 — 0. 921
H = aDb 1. 468 0. 777 — 0. 981
H = a + blnD - 10. 221 9. 216 — 0. 860
H = aebD 5. 213 0. 042 — 0. 789
H = a + bD + cD2 1. 133 0. 822 - 0. 005 1 0. 936
模型 H = a + blnD 和 H = aebD拟合方程的相关系
数 R2 较小;H = aD + b,H = aDb和 H = a + bD + cD2
复相关指数 R2 均大于 0. 92,说明其拟合效果相对较
好。5 个模型中方程 H = aDb的相关系数最高,为
0. 981,说明其拟合效果最好。
2. 2. 2 显著性检验
对上述胸径与树高的回归方程 H = aDb进行方差
分析,见表 3,F = 313. 207 > F0. 01 (1,149) =
6. 83,说明方程 H = aDb能极显著地描述林芝云杉生
长过程中树高与胸径的关系。表 3 回归分析中平均树
高的差异,一方面来自于因个体对应的胸径不同引
起;另一方面是由其它随机因素引起,因此回归方差
和残差的比值较大,表明回归方程愈显著[10 - 11]。
表 3 方差分析表
项目 自由度 离差平方和
回归分析 1 7 682. 136
残差 148 3 630. 052
总计 149 11 312. 190
项目 均方 F值 Significance F
回归分析 768 2. 136 313. 207 2. 34E - 38
残差 24. 528
总计
3 结论
林芝云杉的横向和纵向生长过程与其他乔木一
样,要经历慢、快、慢的过程,幂指数的变化趋势能
够反映这一变化规律。在排除环境发生剧烈变化的前
提条件下,植物的生长量是可预见的[12],应用幂指
数生长方程完全可以对此做出预测[13]。
本研究通过对西藏高原林芝云杉不同径阶个体树
高与胸径的测定,应用数学模型拟合树高与胸径的异
速生长方程。结果表明:林芝云杉的树高与胸径关系
的最佳模型为 H = aDb,具体方程为 H = 1. 68D0. 777,
相关系数为 R2 = 0. 981。该模型可用于西藏高原东南
部快速准确地测算林芝云杉平均树高,也可为计算立
木材积[14]、编制出材率表[15]和评价林分的立地质
量[16 - 18]提供基础数据,这对指导当地林业生产具有
重要的现实意义。
参考文献:
[1]赵敏,丁慧勇 . 上海水杉树高与胸径关系模式分析 [J] . 上
海师范大学学报: 自然科学版,2009,38 ( 5) : 531 - 535.
[2]龙本楼,林顺根 . 黔东南杉木树高与胸径相关研究 [J] . 中
国新技术新产品,2011 ( 5) : 237.
[3]肖志坚 . 福建柏树高—胸径一元线性模型的改进 [J] . 宁夏
农林科技,2011,52 ( 7) : 56 - 57.
[4]罗建,方江平,王国严 . 林芝云杉群落特征的研究 [J] . 热
带亚热带植物学报,2011,19 ( 2) : 113 - 119.
[5]方江平 . 西藏南伊沟林芝云杉林生物量与生产力研究 [J] .
林业科学研究,2012,25 ( 5) : 582 - 589.
[6]方江平,项文化,刘韶辉 . 西藏南伊沟原始林芝云杉林水文学
过程的水化学特征 [J] . 林业科学,2010,46 ( 9) : 14 - 19.
[7]肖文发,韩景军,郭志华,等 . 西藏林芝云杉针叶净光合速率
对环境因子的响应 [J] . 林业科学研究,2003,16 ( 3) : 299
( 下转第 66 页)
·51·
2014 年第 1 期 中南 林 业 调 查 规 划 第 33 卷
[4]国家海洋局,中国新闻网 . 全国海洋功能区划: 2011—2020
年 [EB /OL] . ( 2012 - 04 - 25) [2014 - 01 - 16] . http: / /
www. chinanews. com /gn /2012 /04 - 25 /3846144. shtml.
[5]绿色中国 . 消逝的候鸟生命线 [EB /OL] . ( 2012 - 11 - 20 )
[2014 - 01 - 16] . http: / / greenchina. tv /news - 4190. xhtml
[6]中国经营报 . 湖泊开发亟须划定一条红线 [EB /OL] . ( 2012
- 06 - 02) [2014 - 01 - 16] . http: / / finance. sina. com. cn /
roll /20120602 /003512206065. shtml.
[7]新华社 . 全国土地利用总体规划纲要: 2006—2020 年 [EB /
OL] . 北京: ( 2008 - 10 - 24) [2014 - 01 - 16] . http: / /
www. mlr. gov. cn /xwdt / jrxw /200810 / t20081024_ 111040. htm.
[8]中国绿色时报 . 全国林地保护利用规划纲要: 2010—2020 年
[EB /OL] . ( 2010 - 08 - 26) [2014 - 01 - 16] . http: / /
www. greentimes. com /green /news / lscy /cyzc /content /2010 - 08 /
26 /content_ 102976. htm.
[9]人民网 . 致公党中央调研: 像保卫耕地红线一样保卫湿地红
线 EB /OL] . ( 2010 - 08 - 05) [2014 - 01 - 16] . http: / / cp-
pcc. people. com. cn /GB /45580 /45600 /12357540. html.
[10]人民网 . 划定湿地红线保护“地球之肾”[EB /OL] . ( 2012
- 03 - 14) [2014 - 01 - 16] . http: / /henan. people. com. cn /
news /2012 /03 /14 /602709. html.
[11]中国证卷网 . 我国湿地萎缩危及 3 亿人生存 专家建议划红线
[EB /OL] . ( 2013 - 02 - 03) [2014 - 01 - 16] . http: / /
house. focus. cn /news /2013 - 02 - 03 /2833454. html.
[12]中国网 . 执着与追求我国的湿地保护: 访国家林业局湿地办
公室主任马广仁 [EB /OL] . ( 2013 - 02 - 03) [2014 - 01 -
16] . http: / /www. china. com. cn / fangtan /zhuanti / sdxf /2012 -
05 /17 /content_ 25408971. htm.
[13]南方周未 . 候鸟迷失呼唤“湿地红线”[EB /OL] . ( 2012 -
10 - 09) [2014 - 01 - 16] . http: / /www. infzm. com /content /
81638.
[14]中国人民政府网站 . 回良玉: 推进湿地保护与恢复 推动生态
保护事业 [EB /OL] . ( 2012 -12 -14) [2014 - 01 - 16]. ht-
tp: / /money. 163. com /12 /1214 /15 /8IMP7JNN0025
櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴櫴
3B0H. html.
( 上接第 15 页)
- 305.
[8]陈克龙,刘晓宏 . 喜马拉雅冷杉和林芝云杉年轮稳定碳同位素
气候意义比较 [J] . 盐湖研究,2003,11 ( 2) : 29 - 34.
[9]韩景军,肖文发,郭泉水,等 . 西藏林芝县林芝云杉幼林更新
与物种多样性指数研究 [J] . 林业科学,2002,38 ( 5) : 166
- 168.
[10]郑传英 . 四川桤木林分各主要生长因子间相关性研究 [J] .
安徽农学通报: 上半月刊,2010,16 ( 17) : 161 - 162.
[11]明安刚,唐继新,于浩龙,等 . 桂西南米老排人工林单株生
物量回归模型 [J] . 林业资源管理,2011 ( 6 ) : 83 - 87,
93.
[12]张彩琴,杨持 . 植物生长模拟与数学模型研究 [J] . 内蒙古
大学学报: 自然科学版,2006,3 ( 4) : 435 - 440.
[13]谭美,王四清 . 观赏植物生长模拟模型研究进展 [J] . 园艺
学报,2010,37 ( 9) : 1523 - 1530.
[14]马友平,沈作奎,艾训儒,等 . 材积与胸径关系的 Richards衍
生模型研究 [J] . 生物数学学报,2010,25 ( 1) : 122 - 126.
[15]吴明山,胥辉 . 度量误差对材积模型的影响及参数估计研究
[J] . 北京林业大学学报,2008,30 ( 5) : 83 - 86.
[16]郭晋平,张浩宇,张芸香 . 森林立地质量评价的可变生长截
距模型与应用 [J] . 林业科学,2007,43 ( 10) : 8 - 13.
[17]骆期邦,吴志德,蒋菊生,等 . 用于立地质量评价的杉木标
准蓄积量收获模型 [J] . 林业科学研究,1989,2 ( 5) : 447
- 453.
[18]赵总,张兆国,谢德兵,等 . 立地条件对川滇桤木生长的影
响 [J] . 江苏农业科学,2011,39 ( 6) : 312 - 313.
·66·