森林资源系统的计算机仿真与现代化管理
森林资源系统的计算机仿真与现代化管理 李 昀** （北京林业大学信息学院，100083） 崔赛华 （北京林业大学信息学院，100083） 郑治刚 （国家林业局规划院，100714） 摘要 本文论述了森林资源系统的计算机仿真理论和方法，它是森林资源现代化管理的理 论基础和最主要的系统方法之一。 关键词 森林资源 模型 计算机仿真 现代化管理 1 引言 森林资源系统的计算机仿真，就是借助描述森林资源系统的数学模型，在计算机上对 森林资源系统的动态变化规律和生灭过程进行各种模拟、试验和计算分析。由于计算机 仿真具有模拟、试验、协调、控制、优化、运筹、决策等功能，使得其成为现代化管理 理论的依据和最主要的系统方法之一。因此，森林资源系统的计算机仿真是林业科学工 作和系统理论工作者研究的热点问题。 目前，无论是国外还是国内所进行的森林资源的计算机仿真，所使用的数学模型不是 经验模型就是统计模型，并多为常微分方程。它们都不是根据森林资源的动态变化规律 和生灭过程而推出的数理方程，这些模型都不能从空间和时间状态同时对森林资源的“生 ”、“长”、“死”及动态变化规律作出深刻而准确的描述和表达。因此，用这样的数学模型 进行的计算机仿真不可能对森林资源的动态变化规律和生灭过程作深刻、精确的描述和 表达，不能对森林资源进行准确的模拟、试验和分析，用这样的模型对森林资源进行计 算机仿真很难建立起森林资源现代化管理的计算机仿真理论和技术方法。 本文论述的森林资源系统计算机仿真所使用的数学模型是我们在已完成的国家自然科 学基金项目中所推导的描述森林资源动态变化和生灭过程的数理方程，这些模型是维数 不等的偏微分方程，描述了森林资源的状态空间和时间空间内的动态变化规律、生灭过 程和运动轨迹。因此，森林资源系统的计算机仿真是属于分布参数系统的计算机仿真。 2 森林资源系统仿真的数学模型 森林资源是森林的数量和质量的统一标度和表达，是森林具有财产和价值的属性部分 。目前，森林资源一般是指森林的面积资源和立木蓄积资源。森林资源是由森林各种各 样维数不等的数量分布结构组成，森林的这种数量分布结构称为森林资源分布结构。森 林资源是森林数量分布结构的总称，是各种特定（具体）的森林资源分布结构的集合体 ，每一种特定的森林资源分布结构表示了一种特定（具体）的森林资源。 根据森林资源是由许许多多一维、二维……n维等数量分布结构组成，并依据森林资源 分布结构的动态变化规律和生灭过程推导出了描述森林资源分布结构动态变化规律和生 灭过程的数理方程。 (1)同龄纯林的林龄面积分布结构变化方程： [pic] (1 ) (2)林木直径分布结构结构变化方程： [pic] (2 ) (3)林木株数的直径·树高分布结构变化方程： [pic](3 ) (4)混交林的林木直径分布结构变化方程： [pic] (4 ) (5)混交异龄林的林木株数的直径·林龄分布结构变化方程： [pic] (5 ) (6)混交异龄林的面积按混交系数·平均年龄分布的分布结构变化方程： [pic] (6 ) (7)森林资源分布结构变化方程的一般形式（通式方程） [pic] (7 ) 以上均为方程的连续型模型。森林资源的计算机仿真就是用上述方程的离散型模型。 离散型模型一般可用向量的形式或矩阵的形式表示，由于篇幅关系上述方程的离散型模 型从略。 在上述方程中，（1）至（6）式是描述森林资源某些具体分布结构的动态变化规律和 生灭过程的数理方程；（7）式是描述森林资源各种连续型、阶梯型分布结构的动态变化 规律和生灭过程的数理方程，是通式方程，（1）至（6）式是（7）式的具体化。由于森 林资源是由各种分布结构组成，因此森林资源就得到了描述。 上述所有模型都把自己所描述的系统（分布结构）在t时的状态（P(t)），起始状态 （p(o)）、t时刻生的状态（φ(t)）、死的状态（μ(t)）、长的状态（ψ(t)）等有机的连 系在一起，充分地表示了它们之间的数量关系，描述了系统与个体、个体与个体之间的 相互数量关系及生灭过程的数量关系，深刻地揭示了森林资源系统的动态变化规律和生 灭过程，描述了森林资源在时间空间上运动的轨迹。同时这些模型都有解，且解唯一， 并具有连续依赖变量而变化。这些都为森林资源的计算机仿真提供了可靠的数学模型， 坚实的理论依据和数理方法。 3 森林资源计算机仿真计算的主要内容 森林资源系统计算机仿真的关键是对偏微分方程进行求解，即如何将数学模型变为仿 真模型。由仿真模型可以研制适合于森林资源特点的分布参数系统仿真软件，从而在计 算机上实现仿真计算。本文对如何具体编制仿真软件不作讨论，只对如何通过仿真实现 各种仿真目标进行讨论。 3.1森林资源仿真计算的主要目标 (1)模拟森林资源的动态变化和生灭过程的仿真计算 描述森林资源分布结构动态变化和生灭过程的连续型模型离散化变为离散型模型。将 其任意时间t的起始（初值）状态p(0)、生的状态φ(t)、长的状态Ψ(t)、死的状态μ(t)等 值代入方程，可得到系统（森林资源分布结构）在各个不同时间，不同的生、长、死状 态下的状态值，从而在计算机上实现森林资源分布结构动态度化和生灭过程的模拟仿真 计算。 ①森林资源的监测与预测计算。监测就是通过各种调查、统计、遥感等方法和手段， 即时获得各个时间段上系统的造林更新率β(t)、保存率γ(t)、生长率ψ(t)、死亡率μ(t) 等参数代入方程，计算出各时间段上森林资源（系统）的状态值，进行动态监测。利用 模型进行预测计算，其实质就是对模型的基本函数γ(t)、β(t)、ψ(t)、φ(t)、μ(t)等进 行外推，外推值代入方程计算，实现对某系统未来的预测。 ②森林资源现状分析计算。 ③森林资源再现计算。有了森林资源现在的状态值S（k）和k—1时到K时的ψ（k- 1）、φ(k-1)、μ(k-1)的值就可以计算过去（k- 1）的森林资源的状态值，实现森林资源的再现。 ④计算指标参数体系，对森林资源进行评价。由于描述森林资源分布结构的数学模型 是数理方程即理论模型，因而所有描述森林资源质量和数量的指标参数都可以通过它的 状态函数而计算出来。通过模拟仿真计算，只要计算出它任意时间后状态值，就可以计 算出相应时间的各种指标参数，形成指标参数体系。可以实现随时随地对森林资源进行 科学的评价。 (2)对经营方案，措施等进行优化的仿真计算 模型中的γ(t)、β(t)、φ(t)、Ψ(t)、μ(t)等基本函数值的大小，在实际生产中和森林 资源经营和管理活动中的经营措施和经营方案等有密切的关系。而这些基本函数的值在 已知t时的状态值P(t)，起始状态（初值）的值P(0)及其它基本函数值的情况下，无论是 通过方程的连续型模型还是离散型模型都可以求出。在实际生产中就是在已知经营森林 资源者所需要的状态或要达到的状态P(t)，以及目前的状态即起始状态P(0)，如何从目 前的状态P(0)达到所需要或要达到的状态P(t)，应采取什么经营措施，即通过采（μ）多 少，造（φ）多少，长（φ）多少才能达到所要求的状态，即通过什么样的经营措施和方 案才能达到。这些都可以通过方程分别求得每一个基本函数的值，该值就是应该采取的 措施值，这一完整的经营方案，经营措施的优化过程，通过系统的仿真模型编制成仿真 软件，在计算机上通过仿真计算完全可以实现。并且可以对多个方案，各种可能的经营 措施进行模拟、试验、对比等仿真计算，寻求最优或最佳方案和措施。这是计算机仿真 的一大优点，人工计算是不可想象的。 (3)控制的仿真计算 在控制论中，许多复杂系统的控制问题都是利用系统的仿真计算而实现可控计算和讨 论，把仿真计算作为实现可控计算的工具和方法之一。同样森林资源的结构调整以及森 林资源的最优控制等都可以通过描述森林资源的数理方程在计算机上进行仿真计算来实 现。 森林资源的结构调整其实质是森林资源的分布结构调整。就是由现在的一种不合理（ 不理想）的分布结构状态出发，在一定时间内或在最短的时间内通过不断的调整控制变 量，平稳过渡到理想（合理）的结构状态，是一个控制问题。在林业实际生产中，森林 资源分布结构的调整主要是通过采伐和造林更新达到调控的目的，在目前的科技水平下 林木的生长速度还无法做到控制。 森林资源的控制就是通过可控制变量的取值范围和大小进行调控，在使得路径最短（ 最优）或时间最短（最优）或结构最佳等的情况下由目前的状态过渡到理论或要求的状 态。这一讨论和分析过程同样也可以通过仿真计算而实现。例如森林资源的可持续发展 的定量描述只有通过描述森林资源的数理方程而建立的控制理论和方法、仿真理论和方 法才能建立和完成。同样，森林资源与环境和人类社会关系的协调、平衡和发展及其实 现也只有通过森林资源的控制理论和方法与森林资源的仿真理论和方法才能建立和实现 。 (4)运用森林资源的计算机仿真替代传统的经营方案编制 森林经理学（森林的经营和管理的科学）的主要研究内容之一是编制施业案或称经营 方案等。过去编制经营方案是根据森林的现状（静态）和林分的状况及地形、土壤、气 候等因素，并且根据编制人员理论和经验而确定森林资源今后5年或10年或更长时间对经 营的目的（要达到目标），应采取的措施如抚育、间伐等编制一个实施方案。这样编制 出的方案是一个静态、一成不变的东西，是凭人的某种主观意识和经验而编制的。而森 林资源及其环境时刻都在变化，其中很多因素是人无法控制和预知的。因此，用一个既 定的方案去要求一个分复杂而不断变化着的系统以实现预先的目标，这是很难想象的。 事实上生产中几乎所有编制的经营方案都无法实践，都行不通。借助森林资源的数理模 型，在计算机实现森林资源系统动态变化和模拟仿真，从而达到制定各种最佳方案并加 以实施，做到时时掌握森林资源系统的变化情况以及发展趋势。因此，用建立森林资源 计算机仿真计算和决策支持系统代替经营方案及其编制完全可能和可行。 3.2 森林资源仿真的辅助计算 为了实现各种仿真计算的目标，需要进行仿真的辅助计算，以满足仿真计算所需的参 数，指标、数据及处理数据的技术方法，使得仿真计算有较好的效果。其内容主要有基 本参数的提取，参数的延伸、指标参数体系的计算，将n阶转移矩阵变为一阶转移矩阵的 计算等。 (1)基本参数的提取的计算 在进行森林资源的计算机仿真计算时，需要基本参数即郁闭成林率φ(t)，生长率ψ(t )，死亡率μ(t)才能进行各种模拟计算。这些参数的获取主要是通过森林资源的各种调查 资料如连续定位调查的材料；各种图面材料；遥感的卫片或航片；森林资源的各种档案 等，从它们当中进行提取、统计而获得。因此形成了较为完整的参数提取方法。 (2)基本函数值外延的计算 利用数学模型对森林资源进行预测，其实就是对模型的基本函数μ(死)、φ（生）、ψ （长）等进行外推。根据这些基本函数的物理意义、特点和规律可以从许多外推的数学 方法中选取适合于基本函数的特点和需要的外推方法进行外推计算。其方法一般有平均 法、曲线拟合法、参数辩识法、时间序列法等等。 (3)指标参数体系的计算 在森林资源仿真计算时，不仅需要计算出森林资源的状态值而且需要计算出许多评价 森林资源质量和数量的各种指标参数，而每一项指标参数是从某一方面和某一角度定量 的刻画系统的动态变化规律和特点。而所有的各种指标参数构成了指标参数体系，指标 参数体系是从各个方面较完整地描述了系统数量和质量的特点。对各种指标参数的计算 给出公式或计算方法并编成软件，在森林资源仿真计算时可随时调用。 (4)将n阶转移矩阵变为一阶转移矩阵的计算 在进行森林资源系统分析中，特别是在林业资源预测或监测中，怎样将步长为n年的 转移频率矩阵变为一年的平均转移频率矩阵，以便达到逐年对森林资源进行系统分析以 及预测、监测、结构调整和其它用途的目的，是一项经常遇到的重要实际问题和技术方 法。目前我国的森林资源调查都有一定的间隔期，一般是每五年进行一次或更长时间， 由此得到的各种转移矩阵其步长为5年或更长。而生产中需要对森林资源的变化进行每年 的分析，这就需要如何将步长为n年的转移矩阵变为步长为1年的转移矩阵。有了这样的 计算公式及其软件就能顺利进行森林资源计算机仿真中的各种数值的年度计算。 4 森林资源的计算机仿真在森林资源中的作用 由于森林资源仿真具有模拟、试验、控制、协调、优化、决策等功能和作用，因此， 在森林资源的现代管理中具有重要的理论意义和十分广泛的用途。其作用如下： (1)森林资源的现代化管理提供了理论依据 通过森林资源的计算机...
森林资源系统的计算机仿真与现代化管理