兴安落叶松树高曲线模型研究

()文章编号:1005-5215201705-0036-04

ProtectionForestScienceandTechnology

防 护 林 科 技

,Ma2017y

(S)No.5umNo.164

兴安落叶松树高曲线模型研究

李星硕

()辽宁省本溪市明山区林业工作总站,辽宁本溪117000

摘 要 应用孟家岗林场2根据经验方程选择单木树高曲线的基011—2012年17块固定样地中85株解析木资料,本模型,应用选定的模型对不同立地条件下的落叶松人工林进行树高曲线拟合,再用参数化的方法确定参数与各个林分调查因子之间的关系,从而建立单木树高曲线模型。结果表明:Richards理论模型可作为落叶松人工林单)木树高曲线基本模型,其参数与林分的地位级指数(呈线性关系,而与年龄和密度的关系不明显;最终建立的SCI

0.05D-()树高曲线方程为H=(其优点是,当0≤H≤1.3时,D的12.38025+0.74079SCI)1-e0.66529+1.3,

取值可以为零,符合林木生长的生物学特性。关键词 兴安落叶松;单木;树高曲线模型

:/中图分类号:S791.22 文献标识码:A doi10.13601.issn.1005-5215.2017.05.012j

HeihtCurveModelofLarixgmeliniig

(,M,,,)ForestrorkStationinshanDistrictBenxiCitLiaoninrovinceBenxi117000,ChinayWgygP

LiXinshuog

17permanentsamlefrom2011to2012,Larixgmeliniplantationunderdifferentsiteconditionswereconductedp;heihtcurvefittinliedbelectedmodeltherelationshimonarameters&surveactorsofeachstandggappyspagpyf,wasdeterminedbarametricmethodandfinallthesinlecurvewasestablished.ResultshowsthattheRichardsypyg,modelcanbeusedasthebasicmodelofthesinlecurveofLarixgmeliniiplantationanditsparametersarelinearg-

Abstract Alinataobtainedfrom85basemodelofindividualtreewhichwasselectedbmiricaleuationinppygdyepq

)lelatedtothepositionalindex(SCIofstand.Therelationshietweenaeanddensitisnotobvious.Thehihyrpbgyg

,wisthatwhen0≤H≤1.3,thevalueofDcanbezerohichisconsistentwiththebioloicalcharacteristicsoftreegrowth.g

;Keords Larixgmelini;individualtreeheihtcurvemodelgyw

0.05D-:)()curvilineareuationistheadvantaeofthismodelH=(12.38025+0.74079SCI1-e0.66529+1.3;qg

其模 树高曲线在林业生产和实践中广泛应用,

型在林分材积表、径阶材种出材量预测、经济评价及林分蓄积量的计算等方面都有着重要的作用,国内外学者相继研究过树高曲线随林分和时间变化的规律,建立了以直径、树高关系为基础的林分树高曲线

1]

。建立树高-胸径曲线普遍规律的标准树高曲线[

1 材料与方法

1.1 数据来源

的方法主要有两种:一种是图解法,多采用随手绘制;另一种是数式法,选择适合的经验回归模型,配合回归得到相对最优的回归方程,以此来拟合各因子之间的相关关系。

收稿日期:2017-03-09,李星硕(男,辽宁朝阳人,大学,助理 作者简介:1985-)

工程师,现从事落叶松研究.

图1 兴安落叶松人工林单木树高与胸径的关系

Copyright©博看网 www.bookan.com.cn. All Rights Reserved. 第5期 李星硕 兴安落叶松树高曲线模型研究 固定样地内的680株样木绘制落叶松人工林单木树,高与胸径之间的关系(如图1)然后用上述经验模型对落叶松人工林单木树高生长模型进行拟合,选()看统计指标所表征的模型实际拟合效果。从生2

物学意义和统计指标两方面来选择最适合落叶松人1.2 基本树高曲线模型的选择

传统森林测树学中树高-胸径曲线称为一般树高曲线模型,常用的有以下几种:

]2

。工林单木树高生长的基本模型[

37

本研究应用孟家岗林场2011—2012年的17块

参数

1.3 基本树高曲线模型再参数化

利用上述方法选出建立落叶松人工林单木树高出发,导入适宜的林分因子建立H=fH=f(D)

3]

(。具体方法就是在每林分因子)的树高曲线[D,

块固定样地内,用落叶松人工林单木树高曲线的基础模型分别进行拟合,每块样地都得出不同的参数值,通过观察各个参数与林分各调查因子的散点图,揭示各个参数与林分各调查因子之间的关系,建立能充分表明各个参数与调查因子的模型,最终建立

]4

。落叶松人工林单木树高曲线模型[

(择模型的原则:看模型及其参数的生物学意义;1)

生长模型的基本模型以后,试图用再参数的方法从

D2

;()()抛物线H=H=3a+bD+cD2;42

(a+bD)

bDc0.05Dc--);)Richards式H=a(1-eH=a(1-ebD-);()柱体屈曲式H=aa(1-e7D3;;式中:H-树高(m)D-胸径(cm)a、b-待定

2D;()()直线H=双曲线H=a1a+bD;2(D+b)

2 结果与分析

;()5Loistic式H=g

a;()单分子式H=6cD()1+be-

2.1 基本树高曲线模型的选择

本研究对经验树高曲线,即直线方程、双曲线方程、抛物线方程、单分子式Richards式、Loistic式、g以及柱体屈曲式分别进行了拟合,结果见表1。

模型

a表1 兴安落叶松人工林单木树高曲线的拟合

b0.32014713.404990.1752640.08948

———

—————————0.00527590.39094

———

cH=a+bDaDH=(

D+b)H=

14.3066534.454980.89646851.706122.661978

剩余残差平方和(RSS)

112.07989112.37139113.54495105.34179109.05248205.20407

复相关系数R20.57141460.56252230.55795350.57544330.2011117

D2

(a+bD)2

H=a+bD+cD2

0.05D)c-H=a(1-eH2D3=a16.7200780.5898896

从相关系数和残差综合分析看,抛 由表1可知,

物线理论模型的拟合效果最好,考虑到树高的生长不可能随着直径的生长而不断地增大,树高的生长应该随胸径的生长而生长,而生长的速率应该是逐渐减小的,所以从生物学规律方面考虑,本文选择了Richards理论模型作为落叶松人工林树高单木生长

0.05Dc-),的基本模型,即H=而这个方程有a(1-e这样,即可以保持曲线成“形,又能满足当01.3,S”

[]

≤H≤1.3时,D的取值可以为零的基本要求5。

选定了落叶松人工林单木树高生长曲线的基本模型以后,本文分别对不同年龄、不同立地条件下的落叶松人工林单木树高进行拟合,并分析参数与年2。

龄及各地各龄阶地位级指数的关系。其结果见表

0.05Dc-)+所以,本文将基本方程转换为H=a(1-e一个缺点,即当0≤H≤1.3时,D的取值不能为零,

样地一二三四五六七八九十

7156499647861086081

N104

50505351434541373747

T表2 落叶松人工林不同立地条件下树高曲线的拟合结果

SCI24.9026427.0072326.8046525.0189223.6860424.8813226.2248526.8818528.7245427.83869

18.36618.45121.27918.36117.10219.92820.65519.22219.03720.124

a0.500180.587370.389000.516780.540890.344150.528880.745330.822620.65199

c剩余残差平方和(RSS)

108.59993144.9948661.281685122.8834466.7430192.61858129.9160958.7323568.6510894.97381

0.49357910.65440560.40376140.51214820.64794360.39189130.55591550.74833430.74566710.6228527

相关系数RCopyright©博看网 www.bookan.com.cn. All Rights Reserved.38

十一十二十三十四十五十六十七

102856720917852

45484139384039

18.62318.55918.50920.43618.22719.19117.620

防 护 林 科 技 2017年

26.2597625.7006125.3684725.5222722.5884025.5761722.99497

0.648910.654750.543250.714630.480290.535810.57764

119.72279104.7110387.1118213.7915375.97018127.8737944.51092

0.67049710.68245440.58160920.65676360.66641050.56514660.6231239

从模型拟合的相关系数来看,除了 由表2可知,

标准地六的效果稍差外,其余林分拟合的效果都比较好,模型的相关系数都在0.从剩余残差403以上,平方和的情况看也比较良好。

2.2 单木树高曲线模型的建立

利用表中的a与地位级指数(的数据做散SCI)点图,如图2,发现参数a与地位级指数(呈正SCI)长越快。

同时从图2中可知参数与地位级指数(可SCI)线性:a=bbCI0+1×Sb2指数型:a=bCI1+S分别对这种关系进行拟合,其拟合结果如表3。

),相关(说明立地质量越好,树高曲线增R=0.5531

能的关系有:

2

()抛物线性:a=bbCI+bCI31+2×S3×S()1()2

图2 参数a与地位级指数的关系

模型

)表3 参数a与地位级指数关系模型的拟合结果(n=17

b1

b2

10.033790.947865

0.8152315.152579.22809

a=bbCI0+1×Sa=bbCI+bCI1+2×S3×Sb2a=bCI1+S2

b3

———

———

剩余残差平方和RSS32.0650732.0401127.83251

相关系数R0.5530770.5535650.630488

发现只有在拟合线性方程( 在拟合的过程中,2

)时,所有的参数通过T检验和P检验,而指数-1

型模型和抛物线型模型虽然在相关系数和残差平方检验和P检验。所以,只有线性方程适合于两者的关系。与此同时,在研究的过程中也发现两个参数与年龄和密度指数的关系都不是非常紧密,说明在和上略优于线性模型,但其参数没有完全的通过T

-127.49677-0.37245

而年龄和竞争因子对于树高生长的影响并不是非常突出。

最终本研究所得到的落叶松人工林单木树高曲线模型为:

0.05Dc-)+1.()H=(bbSCI)(1-e340+1

本文用2011-2012两年的孟家岗林场17块固

落叶松人工林中,立地质量对树高生长的影响最大,

参数b0b1c

估计值12.380250.740790.66529

准标差

定样地内的680株样木用Statistic6.0统计软件包对模型进行拟合,其结果见表4。

剩余残差平方和(RSS)

1384.396973

剩余标准差(Sx)y.

1.310943

复相关系数(R2)0.60350415

)表4 落叶松人工林单木树高曲线模型的统计量(n=680

12.2346315.060628.8252

1.0119030.0491870.02308

t-值

p-值

0.000

0.0000.000

落叶松人工林单木树高的生长方程 由此得出,

为:

得到每株样木的树高估计值,从而可由二元表来推算材积。其应用的步骤为:以胸径和树高为目标变量的估计得到各样木或径阶树高估计值后,即可根据胸径和树高来估计以此为解释变量的各类目标变量,如借助二元立木材积表就可以估计材积,利用二元生物量表可以估计生物量,利用二元出材量表就可以估计出材量等。为孟家岗森林资源变档和规划设计提供了依据,也为提高森林资源经营决策提供了技术保障。

H=(12.38025+0.74079SCI)

0.05D-()1-e0.66529+1.3

()5

应用这个树高曲线模型,利用森林资源二类调

查数据,建立落叶松人工林树高生长过程表,应用于孟家岗林场森林经营方案编制和林业调查规划设计等林业生产及科研的实践中,它达到直接的目的就是可以根据样地每木检尺记录及平均树高的测定,

Copyright©博看网 www.bookan.com.cn. All Rights Reserved. 第5期 李星硕 兴安落叶松树高曲线模型研究

在研究落叶松人工林单木树高生长曲线的时候,由于数据的局限性,并没有考虑地域的差异,所以对于模型的适用性并没有很好的验证,而且只有所以对于树高曲线中年龄的考17块固定样地数据,虑并不是全面的,对于模型的适用范围仍需要进一步的研究,在以后的研究中应加大样本容量,这样可

变量H

平均误差(ME)0.05947529

平均绝对误差(MAE)

1.050137

平均相对误差(M%E)

-0.0856015

39

以更好的考虑模型中各个因子相互之间的关联。

3 模型的检验

在α结果见表5。=0.05的水平下进行独立性检验,

平均相对误差绝对值(MA%E)

0.04945907

预估精度(P)0.9942004

本研究用其余300个固定样地样木数据对模型

)表5 落叶松人工林单木树高曲线模型检验统计结果(n=300

4.945907%,模型的预估精度(P)达到了

模型的实际值与预估值之间的差异不99.42004%,显著,说明模型的拟合效果良好。

(为1.平均相对误差(为MAE)050137,M%E)平均相对误差绝对值(为-8.56015%,MA%E)

树高(的模型拟合检验结果为: 由表5可知,H)

平均误差(为0.平均绝对误差ME)05947529,

的散点呈均匀分布,没有发散的现象,且残差都控制在3倍之内,没有出现异常值,说明预估的效果良好。由图4可知,残差的分布近似正态,说明拟合的优度良好,模型可用。

由对预测值和残差所做的分布图3可知,残差

图4 单木树高生长模型残差直方图

0.40136268

F值

显著性不显著

3 结论

3.1 本研究用参数法建立了落叶松人工林单木树

高曲线模型,选择Richards式作为树高曲线的基本模型,根据参数与立地级指数呈线性关系,建立了单

0.05D-()1-e0.66529+1.3。

落叶3.2 从落叶松人工林树高曲线模型建立可知,

木树高曲线模型:H=(12.38025+0.74079SCI)

松人工林树高的生长与胸径呈正相关,胸径越大,树高越高,同时树高与立地质量亦相关,在其他条件相同的条件下立地质量越好,单木的树高越高。而树高的生长与年龄和密度的关系不明显。参考文献:

图3 单木树高生长模型残差的分布

[]惠刚盈,盛炜彤,等.杉木人工林收获模型系统的研1GadowKV,[]刘平,马履一,贾黎明,等.油松人工林单木树高生长模型研究2

[]孙圆,]王万江.江苏省杨树树高曲线模型的研制[应用研究,3J.[]吕勇.]林木树高曲线模型研究[中南林学院学报,4J.1997,17[]王明亮,]李希菲.非线性树高曲线模型的研究[林业科学研5J.

():究,2000,13175-79():486-89():2005,19531-34

[]():林业资源管理,J.2008,10550-55]():究[林业科学研究,J.1994,74353-358

Copyright©博看网 www.bookan.com.cn. All Rights Reserved.