应用线性联立方程组方法建立相容性生物量模型研究.docxVIP

应用线性联立方程组方法建立相容性生物量模型研究* 张会儒 赵冇贤 王学力 王柱明 摘要：森林牛物量是森林生态系统的最基木数量特征，生物量数据是 研究许多林业问题和生态问题的基础，因此，准确测定生物量十分重耍。 建立生物量模型是生物量估测的主要手段。以往所建模型，存在一个 严重的缺陷，即总量与分量模型不相容。如何解决相容性问题，一直 是生物量估计领域所面临的一个难题。本文引入线性联立方程组方法, 较好地解决了相容性问题，并使模型的预佔精度和适用性较以往生物量 模型有了显箸提高。 关键词: 关键词: 生物量模型；相容性；线性联立方程组 中图分类号：S711 文献标识码：A 众所周知，森林生物量测定非常困难，而且耗时费力，因此确定一 种行之有效而又准确的调查方法是I ?分重要的。牛物量模型估计法是 目前比较流行的方法，它是利用林木易测因子來 推算难于测定的林木 生物量，从而减少测定生物量的外业工作。但是以往所建立的生物量模 型，都存在一个严重不足，即模型间不相容。换句话说，就是木材、树 皮、树枝和树叶的模型估计值之和不等于总量模型估计值。例如，在 《东北主要林木生物量手册》⑴屮，胸径10cm紫檄，地上部分总生 物量为16. 392kg,而木材、树皮、枝、叶四部分干重之和为20. 154kg, 二者误差达23%。由于生物量模型的不相容性，从而导致在一次调查 结果屮，出现两个不完全相等的数字，这将给森林资源调查及统计工 作造成数据上的不一致。如何解决相容性问题，一直是生物量估计领 域所面临的一个难题。为此，木研究首次引入线性 联立方程组，探讨 解决这一难题的可行方法。 1数据来源及整理 本次研究的范围只限于立木地上部分。研究的数据来源于吉林省红 石林业局的51株冷杉(Abies fabri (Mast. )Graib.)样木。每株样木 实测的因子有：年龄，直径、树 高、平均冠幅、冠长、材积等。样木 概况见表1。 表1 样木概况表 样木数(株) 胸径(cm) 年龄 1树高1 合计 全称重法 密度法 范围 平均 51 4 47 7. 6?69. 7 37.2 25 ?186 8. 8-25. 3 树干鲜重采用两种方法获得，对于幼树和小树，采用全称重法；对 于大树，采用材积密度法。枝、叶的鲜重测定采用全枝称重法。同时在 野外对干材、树皮、枝和叶抽取样品，在室内通过测定样品的含水率或 密度而推算岀各部分的干重。 2建模方法 在生物量测定吋，总生物量是由各分量(木材、树皮、树枝和树叶) 的合计得到的，因此他们之间是百分之百吻合的。而模型估计结果出现 不相容现象完全是由于各维量独立建模所致。在林业数表领域，一般通 用性回归模型的建立都是针对单一方程而言的，唯生物量模型例夕卜，各 维量模型的建立必须受到“各分量之和等于总量”这一约束条件的制 约。因此，各维量必须联合建模，联立求解，这正好符合线性联立方程 组的思想。线性联立方程组是曲两个或两个以上方程组成的模型系统， 系统中任一个方程的参数估计必须考虑这个系统中其他方 程提供的信 息才能作出⑵。因此，本研究引入线性联立方程组方法建立相容性生物 量模型。 2. 1四个基本概念 为了便于理解线性联立方程组方法，冇必要先清楚以下四个基本概 外生变量：即通常情况下的自变量，本研究屮是指直径D、树 高H、冠幅3、冠KC1、材积V及其组合形式。外主变量在模型拟合 吋，用实际观测值。 内生变量：即通常情况下的因变量，本研究中是指总量及木材、 树皮、枝、叶的干重W 1, W2, W3, W4, W5o在模型拟合时，内生变量 作为因变量时用观测值，作为自变量时用估计值。 简化模型：就是通常所见到的模型形式，即将每个内生变量仅 仅表示为外生变量的关系式。 结构模型：是一组描述系统内部结构关系的方程组，此时，一 个内生变量用其他内生变量和外生变量來表示。 2.2简化模型形式设计 多元线型模型的一般结构形式为： y 二 5+c ixz+c2x2+??? +CpXp (1) 对于生物量模型，因变量y为某维量的干重W,自变量X,的基本形 式有D, II, Cw, Cl, V五个，但其组合形式却很多。本研究经过多次筛 选，确定的自变量系统为:X!=D2, x2=DH, x3=Cw2, xt=Vo由于线性模型 的形式多种多样，各维量模型选型工作量很大，因此，各维量模型均采 用了统一形式。具体形式如下： Wi=Cio+CiiD2+Ci2DH+Ci3Cw2+Ci4V (2) 式中i=l, 2, 3, 4, 5,分别代表总量、木材、树皮、枝、叶各维 量，胜表示各维量的生物量。 2.3结构模型设计 在深入分析各内生变量之间的关系的基础上，并要保证所设计的结 构模型是口J识别的，经过多次筛选，确定的结构模型见下式，此模型为 恰当方程，即正