14实验十四、全林分生长模型的构造-课件

实验十四东北林业大学测树学实验创新性实验

主要内容一、实验目的二、实验工具与资料三、模型的相关信息四、建模的方法步骤五、思考题六、实验报告下一页上一页实验目的

1．掌握用多元回归技术建立全林分生长模型的方法2．掌握用微机编制林分可变密度收获表的方法3．理解全林分模型的特点及林分生长规律

1．掌握用多元回归技术建立全林分生长模型的方法2．掌握用微机编制林分可变密度收获表的方法下一页上一页实验工具与资料计算机及其统计软件（Statistics6.0）。大兴安岭东部地区1990年～2000年复测两次的150块兴安落叶松固定标准地数据（文件名LYSDT.xls）。下一页上一页可变密度生长模型特征因子分类及特点概念全林分生长模型三、模型的相关信息下一页上一页方法步骤

1.建模资料的收集2．资料的整理3.可变密度收获模型4.实验数据下一页上一页全林分生长模型：用以描述全林分总量(如断面积、蓄积量)及平均单株木的生长过程(如平均直径的生长过程)的生长模型称为全林分生长模型（WholeStandModel），亦称第一类模型或全林分模型。下一页上一页三、模型的相关信息特征因子：此类模型是应用最广泛的模型，其特点是以林分总体特征指标为基础，即将林分的生长量或收获量作为林分特征因子如：年龄(A)，立地(SI)，密度(D)及经营措施等的函数来预估整个林分的生长和收获量。这类模型从其形式上并未体现经营措施这一变量，但经营措施是通过对模型中的其它可控变量(如密度和立地条件)的调整而间接体现。这一过程主要通过增加一些附加的输入变量(如间伐方案及施肥等)来调整模型的信息。下一页上一页全林分模型又可分为可变密度的生长模型及正常或平均密度林分的生长模型。林分生长与收获预估模型就是基于这四个因子（年龄(A)，立地(SI)，密度(D)及经营措施）采用生物统计学方法所构造的数学模型。所以，林分生长量或收获量预估模型一般表达式为

式中：Y——林分每公顷的生长量或收获量，A——林分年龄；SI——地位指数或其它立地质量指标；SD——林分密度指标。下一页上一页以林分密度为主要自变量反映平均单株木或林分总体的生长量和收获量动态的模型，称为可变密度的全林分模型。林分密度常用林分断面积（G）、每公顷株数（N）、林分密度指数（SDI）、树冠竞争因子（CCF）等来表示。可变密度的林分生长和收获模型可以预估各种密度林分的生长过程，所以它是合理经营林分的有效工具。由于林分密度随林分年龄而变化，并且林分密度对林分生长的影响又比较复杂。对于林分生长过程，很难找出一个形式简单的模型进行准确地描述。通常采用先拟合含林分密度自变量的林分收获量方程，再依此导出相应的林分生长量方程。但是，随着全林分模型研究的不断深入，在模型系统中同时包含林分生长模型和收获模型，并保证了模型所预估的林分生长量和收获量的一致性。下一页上一页构建林分生长和收获模型所需资料的收集方法，根据收集资料时标准地性质的不同可有一下三种方法：（1）固定标准地长期观测法

对某一树种（组）、某一地域分别不同年龄（t）、不同密度(SD)、不同立地条件（SI）设置符合要求的规定标准地，按一定间隔期（一般未5或10年）进行重复测定，一直到主伐为止，从而获得单木和林分准确的生长过程数据。这是建模收集数据的最佳方法，但这种方法所需时间太长、花费高，所以在实际工作中基本上无法采用这种方法。下一页上一页1.建模资料的收集方法步骤（2）临时标准地短期观测法（一次测定法）

在规定的建模地域范围内，分别树种设置大量临时标准地。临时标准地分布于不同年龄、不同密度和不同立地条件的林分中，实测林分的各调查因子。该法提供资料迅速、花费较少，但是不能合理地反映林木或林分的生长规律和动态信息。下一页上一页

建模时将取自不同林分相同立地条件的标准地予以归类，修匀后作为该立地条件下的林分发育过程。这样做的结果是人为地将不相关的林分进行了组合来反映林分生长过程，只能说明实际林分发育过程的表面现象（平均结果），而很难从本质上揭示林分生长的内在规律，更甚者会得到错误的结论。我国固定标准地少，以往多采用此方法。下一页上一页（3）固定样地和临时样地相结合的综合法

在不同年龄、立地的同类林分中设置一定数量的固定标准地，每块标准地进行短期（3～5次）重复测定，并结合临时标准地一次测定结果来建立生长和收获模型。近30年来，我国的森林资源连续清查体系在全国设置了40多万固定样地，并且各地结合科学研究设置了一些固定标准地。此方法比较适合于我国实际，建议采用该方法来建立林分生长合收获模型。下一页上一页收集资料前，应拟定计划其内容有：1）确定地域和树种；2）确定标准地的条件、数量；3）确定标准地调查内容与方法。（1）标准地设置标准地应分布于不同年龄、不同立地和不同密度，其数量应在200块以上。（2）测定项目：标准地所需测定项目包括：1）林分各调查因子：林分年龄（t）、每木检尺，树高、枝下高、冠幅。建立单木生长模型还需要通过定株观测，测定每株树木的直径、树高、冠幅、冠长和树木的相对位置；2）记载标准地的地形、地势、海拔、植被；3）做土壤剖面进行土壤调查；4）详细记载林分经营历史，尤其是间伐次数、间伐时间及间伐强度等。下一页上一页(1)资料的整理首先应进行检查，检查标准地的设置、测定因子、记录是否符合规定的要求，是否有漏测或漏记因子、数据是否有误。将各种调查数据建立计算机数据库，包括各标准地因子库和每木检尺库，并利用微机分别活立木、枯立木计算标准地各测树因子：林分年龄（t）、林分平均直径（）、平均高（）优势木平均高（HT）、林分断面积（G/hm2）、株数（N/hm2）、蓄积（M/hm2）、枯损量等。将所收集的全部标准地数据,大致按4：1（75%和25%）的比例分成两组独立样本：建模样本和检验样本，分别用于构建和检验林分生长和收获模型。下一页上一页2．资料的整理（2）资料的分析—取舍建模的数据是总体中的一组样本，如有个别过大或过小的异常数据混杂进去，缩减模型的精度会受到影响。为此，必须剔除异常点（测错、计算误差、异常数据）以提高模型的预测精度。下一页上一页异常数据的剔除过程分两步进行：首先，用计算机绘制林分（或林木）各调查因子之间（如年龄(t)和优势木平均高(HT)、林分断面积和蓄积等）的散点图，通过肉眼观察确定出明显远离样点群的数据并删除，这类数据是属于因登记、计算等错误而引起的异常值；其次是根据候选的模型初步拟合建模数据，并绘制模型预估值()与标准化残差()之间的残差图。在残差图中，超出±2倍标准差以外的数据作为极端观测值予以剔除。下一页上一页含有林分密度的收获预估模型主要用于现实收获量的直接预测。根据建模方法的不同可划分为以下三种。下一页上一页3.可变密度收获模型（1）基于多元回归技术的经验方程

典型的可变密度收获预估模型Schumacher(1939)收获模型：

(14－1)

式中M—单位面积林分蓄积量；t—林分年龄；SI—地位指数；G—林分断面积；b0～b4为方程待定参数。下一页上一页这一研究开创了定量分析林分生长和收获量的先河，类似的研究方法沿用至今。之后，许多研究者采用多元回归技术来预测林分生长或收获量。这类可变密度收获模型的基础模型为Schumacher(1939)蓄积收获曲线：

(14－2)下一页上一页

基于（14－2）式构造的可变密度收获模型的一般形式为：

(14－3)

式中f(SI)——地位指数（SI）的函数;g(SD)——林分密度(SD)的函数;a0～a1以及β0～β4为方程参数.

(14－3)式称作Schumacher收获模型。下一页上一页迄今为止许多学者均采用(14－3)式这一模型形式，构建了不同树种的全林分可变密度收获模型，如：（1）美国火炬松天然林（Clutter和Sullivan，1972)（英制单位）

(14－4)（2）台湾二叶松人工林（冯丰隆和罗绍麟，1986)

(14－5)（3）大兴安岭兴安落叶松天然林（蒋伊尹和李凤日，1989)

(14－6)

下一页上一页(2)基于理论生长方程的林分收获模型由于理论生长方程具有良好的解析性和适用性，近30年来，各国倾向于将稳定性较强的林分密度指标引入适用性广的理论生长方程，来建立林分生长和收获预估模型。常用的理论生长方程有：逻辑斯谛(Logistic)方程、单分子(Mitcherlich)式、坎派兹(Gompertz)方程、考尔夫(Korf)方程和理查德(Richards)方程等。下一页上一页许多研究者使用这些理论方程拟合林分生长量和收获量，都取得较好的结果，这也说明这些方程具有较强的通用性和稳定性。从70年代开始，许多研究者开始研究这些方程中的参数与林分密度或单木竞争之间的关系，并将林分密度指标引入这些方程之中，预估各种不同密度林分的生长过程，这样建立的收获模型具有较好的预估效果，使模型也具有更强的通用性。下一页上一页现以理查德方程为例说明利用这种方法建模的基本思路，理查德生长方程的基本形式为：

(14－7)下一页上一页首先分析方程（14－7）中各参数A、k和c与地位指数（SI）和林分密度(SD)之间的关系并建立函数关系，比如将最大值参数A作为立地的函数：；而生长速率参数k主要受林分密度的影响，与SI相关不紧密，故；关于形状参数c与立地条件和林分密度的关系尚无定论。最后，根据所建立的函数关系，采用再次参数化的方法将地位指数（SI）和林分密度(SD)变量引入（14－7）中来构造林分生长和收获预估模型。下一页上一页以这种方法成功地建立可变密度收获模型的实例如下：1)兴安落叶松天然林断面积预估模型（张少昂，1986)假设林分平均单株直径生长潜力的发挥程度随其所占有的林地空间的减少而下降，在此基础上，将林分密度指数(SDl)引入了Richards生长方程，利用临时标准地资料分别地位指数（SI）建立了兴安落叶松天然林林分断面积预估模型：

(14－8)式中G—林分每公顷断面积（SI＝22m时，基准年龄为80年）；S——表示每株林木平均占有的林地面积积，SDI—林分密度指数，。

下一页上一页2）马尾松人工林断面积预估模型（唐守正，1991)在（14－8）基础上，通过分析理查德(Richards)方程中的渐近值(A)与地位指数（SI）的关系，将立地因子引入方程，建立了马尾松人工林断面积预估模型：

(14－9)式中：SI—地位指数（基准年龄20年）；SDI—林分密度指数，。

下一页上一页3）大兴安岭北部地区兴安落叶松天然林断面积生长模型（李凤日，2003）利用SDI作为密度指标，在几个假设的条件下将SDI引入Schumacher方程中，利用688块固定标准地数据建立了落叶松天然林断面积生长模型（14－10）。图14－1为SDI=600时，不同立地（SI=12m～24m）条件下落叶松天然林林分断面积生长曲线。