全 文 ：书第 39 卷 第 9 期 东 北 林 业 大 学 学 报 Vol． 39 No． 9
2011 年 9 月 JOURNAL OF NORTHEAST FORESTRY UNIVERSITY Sep． 2011
1)国家“十一五”科技支撑课题(2006BAD03A0802)和林业公益
性行业科研专项(201004002)。
第一作者简介:马武，男，1986 年 12 月生，中国林业科学研究院
资源信息研究所，硕士研究生。
通信作者:雷相东，中国林业科学研究院资源信息研究所，研究
员。Email:xdlei@ caf． ac． cn。
收稿日期:2011 年 2 月 18 日。
责任编辑:张建华。
基于 CAPSIS平台开发落叶松云冷杉林矩阵生长模型1)
马 武 沈琛琛 雷相东 Samuel Dufour Franois de Coligny
(中国林业科学研究院资源信息研究所，北京，100091) (法国农科院 AMAP实验室)
摘 要 为了提高生长模型的开发效率，扩大生长模型的影响和应用，根据落叶松云冷杉林矩阵生长模型的
参数特点及各种算法，利用计算机 java语言将这些参数和算法转移到 CAPSIS共享森林建模平台上，结合 CAPSIS
本身具有的各种适用于单木或林分的模块属性，基于 CAPSIS 平台开发了落叶松云冷杉林矩阵生长模型，该模型
可用来预测和模拟自然生长及采伐对经济和生态目标的影响。建议将来的生长模型选择类似 CAPSIS 通用的模
拟器来实现，以避免重复开发和扩大模型的应用。
关键词 CAPSIS;共享森林建模平台;径阶生长模型;JAVA语言
分类号 S757． 1
Development of Matrix Growth Model for Larch-Spruce-Fir Forest Based on CAPSIS Platform/Ma Wu，Shen
Chenchen，Lei Xiangdong(Institute of Forest Resource Information Techniques，Chinese Academy of Forestry，Beijing
100091，P． R． China) ;Samuel Dufour，Franois de Coligny(INRA，UMR AMAP，Montpellier，F-34000 France)/ /
Journal of Northeast Forestry University． － 2011，39(9)． － 112 ～ 115
Computer algorithms and parameters of the matrix growth model for larch-spruce-fir forest were transferred to CAPSIS
shared forest-modeling platform by Java language，in order to improve the efficiency of the growth model and expand the
impact and application of growth models． A matrix growth model for larch-spruce-fir forest was developed based on CAPSIS
platform according to a variety of module properties of CAPSIS which are adapted to single trees and stands． The model can
predict and simulate the impact of natural growth and cutting on economic and ecological goals． It suggests that the growth
models should be developed by CAPSIS general purpose simulator in order to avoid the repeated development and expand
the application of the models．
Keywords CAPSIS;Shared forest-modeling platform;Diameter growth models;JAVA language
由于林业的长周期性，为了有效地经营森林，林业工作者
需要用模型来预测林分生长并制定经营决策。许多生长模型
已经用来模拟不同尺度下的林分生长动态［1］。然而，这些单
独开发的具体模型由于它们的系统结构和技术解决方案有很
大的差异，有时候很难被其他人使用，也很难相互交流。为了
在各种情况下都可以重复地使用一些工具，一些团队选择建
立一个更通用的模拟器来实现模型。这些模拟器使用共同的
数据结构和一些功能(如生长、竞争、环境经营等) ，可以用于
多个树种或森林类型。目前只有少数软件如 LMS［2］、
MUSE［3］、Symfor、CAPSIS等提出了一个更加通用的方法，允
许使用、比较和建立不同的生长模型。其中，法国农科院研究
人员开发的 CAPSIS最具有代表性［4 － 8］。本文基于 CAPSIS平
台开发了落叶松云冷杉林径阶生长模型，以预测和模拟自然
生长及采伐对经济和生态目标的影响，为利用通用模拟器实
现具体的生长模型提供了方法和参考。
1 CAPSIS共享森林建模平台简介
CAPSIS是一个共享的森林模拟器，开发于 1994 年［10］。
该项目的目标是建立一个长期、开放和共享的建模平台来实
现:①促进模型的发展及评价;②共享工具和方法;③比较不
同模型得出的结果;④将模型传授给森林经营管理人员并作
为教材进行培训。目前最新版本为 CAPSIS4，已经包括了 50
个不同的模型，如单木、径阶和林分生长模型，木材质量模型，
种子扩散模型等，它提供的森林经营管理工具可以设计和比
较不同的森林经营方案。
CAPSIS围绕内核建立并提供公共服务和通用数据结构
(图 1)。它有一个稳定的内核和一些应用技术库，可以接受
具有不同数据结构、模拟步长和模拟算法的生长模型。通常，
生长模型作为一个模块集成到 CAPSIS 中，这些模块包含的
算法可以应用到随时间变化的模型预测的结构中。CAPSIS
也提供了数据结构库和每个建模者都能使用的处理工具。这
种结构最容易受面向对象编程概念的影响，尤其是类继承。
因此，通用模式可以被模块中描述子类的内核超类来描述，而
且后者(子类)可以继承前者(父类)的属性(变量和函数)并
进行修改和添加。CAPSIS中包含了许多这样的父类，建立模
型时可以大量继承父类中的已有属性(变量和函数)。
图 1 CAPSIS原型结构
为了最大限度地节省人力和提高效率，CAPSIS采用了共
同开发(Co － development)的策略。由于采用 JAVA 语言，很
容易实现调整和改进。通过约定加入 CAPSIS 组织后，首先
对建模人员进行培训，然后由开发人员和建模人员共同进行
设计。建模人员通过远程共享所有代码，并和开发人员交流
得到软件开发需要的帮助。
2 利用 CAPSIS平台开发落叶松云冷杉林矩阵生长模型
2． 1 落叶松云冷杉林矩阵生长模型
矩阵生长模型结构简单、数据容易获得、易校检和检验、易与
优化模型相结合，已得到广泛应用，并用以模拟不同经营方案的
生态和经济效果。矩阵生长模型的关键是生长、进界和枯损 3个
子模型及保留转移概率的确定。落叶松云冷杉林矩阵生长模型
(LSFMGM)以我国东北过伐林区吉林省汪清林业局近天然落叶
松云冷杉为对象，利用 20块样地 4次的调查数据(每次调查间隔
5 a)，建立的落叶松、红松云冷杉(红松、云杉、冷杉)、慢阔(色木、
水曲柳、椴树和枫桦)和中阔(白桦、榆树和杂木)的非线性矩阵
生长模型，可以模拟自然生长和采伐对木材产量、树种和大小多
样性、树木地上碳贮量等多个目标的影响。本文仅侧重于用
CAPSIS平台实现模型的运行及结果输出。
2． 2 模块设计与开发
CAPSIS模型结构的设计可以确保不同种类模型的整合。建
模者可以自由地选择数据结构、经验公式或处理流程、模拟步长、
详细程度和尺度。建模者还可以整合模型，这种模型既可为林分
计算全局变量，也可以描述单木，并可以考虑局部竞争环境。
2． 2． 1 模型实现的基本原理
CAPSIS模型是对一系列数据结构的描述，它的处理过程
相互关联，主要依赖于内核描述。模型中两个主要元素为林
分类和模块类。后者涵盖功能模拟和预测，可以模拟林分的
连续状态，并且可以在新一轮的拟合过程将其添加到项目中。
生长模型作为独立的模块集成，这些模块包含它们自身的数
据结构，是从内核中选择得出的一些模块的延伸。它们还包
含了林分预测模型中的算法，在模块中可以自行设计经营管
理方案或在外部平台上进行相关的扩展，如:间伐、修枝等。
每个模拟过程拥有初始步骤，该步骤与初始林分相对应。该
模拟主要用于处理连续生长的林分，通过模块来预测各个生
长阶段及其扩展中的各项处理措施。
将落叶松云冷杉林矩阵生长模型结合到 CAPSIS 平台
上，主要是通过计算机语言 java 来实现。根据落叶松云冷杉
林矩阵生长模型的参数特点及各种算法，利用 java 语言将这
些参数和算法均转移到 CAPSIS 平台上，结合 CAPSIS 本身具
有的各种适用于单木或林分的模块属性，最后可以精确地将落
叶松云冷杉林矩阵生长模型建立在 CAPSIS平台上，并实现落
叶松云冷杉林矩阵生长模型的各项原有功能和新增的功能。
2． 2． 2 模型实现的基本流程
在对模型分析的基础上进行模块设计，在每个模块中充
分利用 CAPSIS提供的类进行开发实现。模型实现的基本流
程如图 2 所示。
图 2 利用 CAPSIS平台实现落叶松云冷杉林径阶生长模型的流程
2． 2． 3 模块设计与实现
利用计算机语言 java 将落叶松云冷杉林矩阵生长模型
结合到 CAPSIS平台上，设计了初始化、预测和输出 3 个模块。
①初始化模块。运行 CAPSIS 并初始化 LSFMGM 工程，
首先得到的是初始化模块，初始化模块主要是提供需要输入
的初始变量，包括:起测径阶(Minimum diameter)、径阶宽度
(Diameter class width)、采伐周期(Calibration time step)、坡度
(Slope angle)、海拔(Elevation)、坡向(Aspect) ，以及初始直径
分布数据(Diameter distribution file)和选择是否分树种(By
species)进行模拟(图 3)。另外，也可以增加和减少树种名，
修改对应的枯损率模型(Mortality model)、进阶率模型(Up-
growth model)、进界率模型(Ingrowth model) ，以及生物量(Bi-
omass model)和材积方程(Volume model) (图 4)。Capsis mod-
ule classes的 stand class为结合的模型提供了初始化模块，初
始化模块可以根据模型的需要添加变量。完成以上操作后，
点击确定进入预测模块。
②预测模块:Capsis module classes 中的 model class 在预
测模块中为与 capsis结合后的模型提供了预测功能，模块类
包含了功能模拟(function simulation)和预测(evolution with a
loop) ，可以计算林分的连续状态，也可以计算林分模型中的
转移概率。因此，在预测模块中可以添加特定的采伐方案和
进行转移概率计算。根据落叶松云冷杉林矩阵生长模型中的
采伐方案，在预测模块中需要添加的变量包括:模拟期(Evo-
lution duration)、采伐周期(Cutting frequency)、采伐强度(Total
percentage)、最小采伐径阶(Minimun cut diameter class) ，以及
可选参数是否分树种(selectSpeciesMethod)、采伐类型(Cut-
ting type)、采伐方法(Cut Tree selection method) (图 5 ～图 7)。
按要求完成以上操作后，点击确定进入输出模块，在该过程中
根据初始化模块中输入的枯损率模型(Mortality model)、进阶
率模型(Upgrowth model)、进界率模型(Ingrowth model)、生物
311第 9 期 马 武等:基于 CAPSIS平台开发落叶松云冷杉林矩阵生长模型
量(Biomass model)和材积方程(Volume model) ，同时进行转
移概率计算。
图 3 初始化模块窗口
图 4 各树种的转移概率模型及生物量和材积式窗口
图 5 预测模块窗口 1
③输出模块:CAPSIS 提供的林分浏览器和数据提取器
(Stand viewers and data extractors)可以以曲线图、柱状图和表
格等形式输出模拟结果，并允许浏览模拟结果。根据落叶松
云冷杉林矩阵生长模型的输出要求，结合输出模块可以输出
落叶松云冷杉林矩阵生长模型中的各项指标随时间、树种、径
阶变化的曲线或柱状图(图 8)。
图 6 预测模块窗口 2
图 7 预测模块窗口 3
图 8 模型输出模块
3 个模块设计好之后，在每个模块中充分利用 CAPSIS 提
供的类进行开发实现。首先在计算机 D 盘上安装 capsis，然
后建立一个工程 lsfmgm，在 D:\capsis4 \bin \ lsfmgm新建 model
包，在 model包中建立 8 个类，分别为:LSFMInitialParameters、
LSFMEvolutionParameters、LSFMFormula、LSFMMethodProvider、
LSFMModel、LSFMSpecies、LSFMStand、LSFMStandLSFMSpec-
iesTransitionProba，这些类的功能属性如表 1 所示。
通过以上 3 个模块的开发和实现，可以方便地模拟不同
经营方案的经营效果，如模拟 40 a内 4 种不同采伐方案下的
411 东 北 林 业 大 学 学 报 第 39 卷
生物量的变化(图 9)。方案 a:不采伐;方案 b:不分树种，按
株数从大到小采伐，采伐周期 5 a，采伐强度 30%，最小采伐径
阶 20 cm;方案 c:不分树种，按蓄积成比例采伐，采伐周期 10
a，采伐强度 20%，最小采伐径阶 15 cm;方案 d:分树种，按断
面积随机采伐，采伐周期 10 a，采伐强度 20%，最小采伐径阶
15 cm。
表 1 利用 CAPSIS平台实现落叶松云冷杉林径阶生长模型建立的类
结合模型新建的类 新建的类继承或实现 Capsis的类或接口 新 建 类 的 功 能 属 性
LSFMInitialParameters AbstractSetting类 InitialParameters接口、Automatable接口 定义初始化模块中各种初始变量以及不同树种转移概率计算的方程
LSFMEvolutionParameters EvolutionParameters接口、Automatable接口 定义预测模块中各种变量、可选参数(是否分树种、采伐类型、采伐方法)
LSFMStand SimpleScene类 定义预测模块中的各种变量以及参数的运算法则
LSFMMethodProvider MethodProvider接口 调用 CAPSIS中提供的方法进行相关运算
LSFMModel GModel类 对选择不同参数的各种情况进行相关运算
LSFMSpecies 无继承或实现 确定转移概率计算法则(枯损率模型、进阶率模型、进界率模型等)
LSFMSpeciesTransitionProba Cloneable接口、Serializable接口 对转移概率计算公式的 8 个变量赋值
LSFMFormula 无继承或实现 调用转移概率计算公式并进行计算
图 9 4 种不同方案 40 a林分生物量的变化
2． 3 与可视化的连接
CAPSIS可以通过与其他模拟器(如 SVS － Stand Visualization
System)的连接，输出林分的三维可视化较低，可使研究人员能以
更直观的方式发现隐藏在数据中的科学规律。图10为利用Cap-
sis模拟林分生长的可视化结果。通过对林分生长的模拟结果可
视化，可很大程度地方便研究人员检验、观察模拟结果。
图 10 基于 CAPSIS的林分可视化
3 结论与讨论
基于 CAPSIS共享森林建模平台开发了一个落叶松云冷
杉林矩阵生长模型，为利用通用模拟器实现具体的生长模型
提供了方法和实例。一方面提高了开发效率，另外可以扩大
模型的影响和应用，为模型的比较检验提供了好的媒介。建
议将来的生长模型选择类似 CAPSIS 通用的模拟器来实现，
以避免重复开发和扩大模型的应用。
经过多年的实践，CAPSIS项目的成功证实了这个共享森
林建模平台的优势。首先，该共享森林建模平台为形成一个
林业建模团体提供了很好的机会。CAPSIS 项目汇集了法国
与森林建模有关的主要研究机构，并建立了 70 多个合作伙伴
关系，可以促进建模者和林业工作人员的规范并进行调整。
目前该项目也对其他国家的同行进行开放。其次，这样一个
共享森林建模平台能够共享工具和模型。在 CAPSIS 项目
内，定期组织培训班用来促进科学家之间相互交流，优化
CAPSIS的内核和解决 CAPSIS 中常见的技术问题，其中建模
者主要负责模型的结合。由于 CAPSIS 是一个免费的软件
(通用公共许可证) ，它很容易再次使用 CAPSIS 单元内的模
块或内核中的子模型。例如，不少模型通过简单修改相同类
型的已存在模型的演化规律，可以很容易地在新模块中集成。
库也为再次使用建模工具提供了很好的机会。例如，本文建
立的 LSFMGM模型可为将来类似的径阶生长模型的建立提
供了工具。第三，共同森林建模平台是一个非常实用的工具，
可以用来比较不同的模型和造林方案，不同的模型可以整合
到一个单一平台上，并且数据(如初始状态)是兼容的。因
此，在相同的初始状态下很容易比较两个模型的模拟结果，或
在相同的模型中比较多种不同的方案。这些对评价模型有很
大的作用。最后，这些优势确实改变了建模者的设计习惯，并
建议扩大这种共享森林建模平台的作用。
参 考 文 献
［1］ Pretzsch H，Biber P，Dˇursk J，et al． Recommendations for standard-
ized documentation and further development of forest growth simulators
［J］． Forstwissenschaftliches Centralblatt，2002，121(3):138 －151．
［2］ McCarter J B，Wilson J S，Baker P J，et al． Landscape manage-
ment through integration of existing tools and emerging technologies
［J］． Journal of Forestry，1998，96(6) :17 － 23．
［3］ Gignoux J，Menaut J C，Noble I R，et al． A spatial model of savanna
function and dynamics:model description and preliminary results［C/
OL］． Newbery D M，Prins H H T，Brown N D． Dynamics of tropical
communities:the 37th symposium of the British Ecological Society，
Cambridge University，1996，1998． 361 －383．［2010 －12 －10］． http:/ /
www． cabdirect． org /abstracts /19980612809． html?freeview = true#)．
［4］ de Coligny F，Ancelin P，Cornu G，et al． CAPSIS:Computer-ai-
ded projection for strategies in silviculture:advantages of a shared
forest-modelling platform:Reality，models and parameter estima-
tion-the forestry scenario 2 － 5 June 2002［C］． Sesimbra，Portu-
gal:［s． n．］，2002:1 － 4．
［5］ de Coligny F，Ancelin P，Cornu G，et al． CAPSIS:Computer-ai-
ded projection for strategies in silviculture:Open architecture for a
shared forest-modelling platform:Fourth Workshop IUFRO S5． 01．
04，British Columbia，Canada-Septeber 8 － 15，2002［C］． Colum-
bia，Canada:Harrison Hot Springs，2002:371 － 380．
［6］ de Coligny F． Efficient building of forestry modelling software with the
capsis methodology［C］/ /Fourcaud T，Zhang X P． Proceedings of the
second international symposium on plant growth modelling，simulation，
visualization and applications，PMA06，13 －17 /11 /2006，Beijing． Los
Alamitos:IEEE Computer Society，2007:216 －222．
［7］ Courbaud B，Goreaud F，Dreyfus Ph，et al． Evaluating thinning
strategies using a tree distance dependent growth model:some ex-
amples based on the CAPSIS software“uneven-aged spruce forests”
module［J］． Forest Ecology and Management，2001，145(1 /2):15 －28．
［8］ Cucchi V，Meredieu C，Stokes A，et al． Modelling the windthrow risk
for simulated forest stands of Maritime pine(Pinus pinaster Ait．)［J］．
Forest Ecology and Management，2005，213(1 /3):184 －196．
511第 9 期 马 武等:基于 CAPSIS平台开发落叶松云冷杉林矩阵生长模型