线性回归模型在直径生长和径阶分布方面的应用.docxVIP

线性回归模型在直径生长和径阶分布方面的应用 随着林业一体化经营水平的提高和发展，我们必须及时、准确地估算森林资源的数量，客观、科学地评价森林资源的质量，然后继续开发和开发森林资源。充分发挥森林的多功能效益是森林经营领域的一个永恒课题。 林分蓄积量的增长是林分生长与枯损综合作用的结果,林分枯损率(量)的估测,将为林分蓄积量的增长和消耗提供科学方法和手段,是林业宏观决策和发展规划的基础,对于林业生产实践亦有重要的现实意义。 1 枯损模型估测技术 目前,国内外林分生长量的积累和发展规律的研究,无论是在基础理论领域,还是在实用技术方法均已经取得丰硕的成果。林分枯损方面的研究与探索,其科学理论和技术方法也具有一定的规模。60年代中期,以铃木太七为代表的日本林学家,曾经运用随机过程理论,研究和分析了林分直径分布的转移概率理论。进入80年代末期龙原哲和箕轮光博根据概率论,描述了林分各径阶生长与枯损的转移过程。美国林学家在建立国家森林资源连续清查体系的同时,利用林木直径生长率的特性建立了估测林木枯损率的直径枯损模型。 综合林分枯损率(量)的估测方法,可大致划分为两种,即固定样地法和收获表法。迄今为止,国外常用的森林枯损率(量)的估测方法不外乎采用如下技术手段:(1)线性或外线性回归分析;(2)林分直径分布模型;(3)利用理论生长曲线模拟枯损量;(4)以差分微分方程直接估测枯损量;(5)利用马尔柯夫链模型估计枯损量;(6)利用竞争指数及直径自变量估计枯损量;(7)建立存活率预估模型间接估测枯损量。我国林分枯损估测技术的研究与探索,尚在初期阶段。本研究只是在此领域所做的尝试。 2 固定样地的清查 本文研究资料均来源于辽宁省森林资源连续清查的固定样地体系。共收集了1978年至1995年期间4次清查固定样地48块,其中人工林6块,天然次生林42块,计算木5 000余株。 利用上述资料,分别样地和树种计算平均直径、直径方差及直径生长率等标志值,其具体方法和过程分述如下。 2.1 直径生长率及等级 运用普雪斯勒公式(Pressler 1857),计算各样木 5 a间直径生长率,并且将直径生长率划分成若干等级。 PD(％)=(Da?Da?n)/(Da+Da?n)×(200/5)％(1)ΡD(％)=(Da-Da-n)/(Da+Da-n)×(200/5)％(1) 2.2 直径回归模型的建立 以各年度直径数据,分别样地模拟,构成不同年度之间的一元回归模型。 y=ax+b(2) 式中:y为n年时直径值,x为n年前直径值。 2.3 平均直径、离散度和条件因子的计算 根据方差定理计算(2)式的条件方差,其计算公式表示如下: 3 枯损分布林木株数枯损模型 本研究选择直径为标志值,分别样地、树种,建立林分枯损模型。 3.1 株数枯损模型 林木株数枯损及其按径级分布的规律,必须符合以下3种假定: (1)令t时刻直径为x,t+△t时刻直径为y,则下式成立。 yi=f(Xi)+ei(5) 式中:f(Xi)为单木生长式;ei为随机误差项。 (2)随机误差项ei遵从某一概率分布。 (3)根据转移概率密度函数定义,在数值上,令(t,t+△t)期间内,当yx时林木为枯死木。其枯死概率计算如下: 式中:m(t,△t,x)为林木枯死概率;P(Y/x)为林木现实直径转移概率。 令t及t+△t时刻的直径分布概率可用Φ(x)、Φ(y)描述,根据转移概率密度函数定义: 若令t时刻直径为x的林木在t+△t时刻枯死的株数为M(t,△t,x),林分总株数为Nm(t,△t),则得: 3.1.1 生长式与枯损分布 林木不同时期的直径(t,△t)、(t+△t,y)之间,根据Prodam、大隅真一等人的研究多呈直线关系。其直径生长式表示如下:y=ax+b(10) 直径生长式的误差变动,根据实测资料,利用(10)式检验证明遵从正态分布,因此,直径转移概率函数表示如下: 故其枯损概率m(t+△t,x)则为: 3.1.2 直径分布与枯损分布 若假定林木初期直径分布Φ(x)的概率密度为正态分布,则分布函数得: 根据(11)式和(13)式,得到末期直径分布的概率密度函数: 在已知t+△t时刻直径为y的林木,作△t以前直径x的概率,根据条件概率定义可得: 将(4)式和(14)式代入(15)式可得: 若设(t,t+△t)期间林木枯死的概率为m(y),生存概率为S(y),根据假定可得: 故末期直径分布为: 式中Φ(y)为末期直径分布;N为初期立木株数。 3.1.3 株数枯损模型的建立 利用(17)、(18)式,按径阶求出相应株数枯损率,分别拟合如下模型: 式中:R为径阶株数枯损率;D1.3为径阶中值;a、b、c、d为待估参数。各样地、树种的拟合结果详见表1。 4 模型的拟合精度和适用性 本文还用直径、断面积及收获表方