森林抚育间伐后林分的总生产量最大的间伐优化设计方法
【专利摘要】本发明公开了一种森林抚育间伐后林分的总生产量最大的间伐优化设计方法，属于森林抚育间伐技术领域，根据本发明中提出的一系列模型，以及所建立线性规划目标函数和约束方程，利用OFFICE的EXCEL中的“规划求解”选项，求解拟开展抚育间伐林分的各树种各径级林木的间伐株数和保留株数，并结合设计经验和生产要求以及间伐效果，可通过所提供的部分约束值的人?机交互，修正方案，使得间伐后林分结构最好、间伐间隔期末的林分总生产量最大，获取满意的间伐设计方案。据此发明得到的优化设计方案，用于开展森林抚育间伐生产，能使森林的结构最好、产量最大，从而最大限度地发挥森林的生态功能和经济效益。
【专利说明】
森林抚育间伐后林分的总生产量最大的间伐优化设计方法
技术领域
[0001 ]本发明涉及森林抚育间伐技术领域，特别涉及一种森林抚育间伐后林分的总生产 量最大的间伐优化设计方法。
【背景技术】
[0002] 森林结构和生长、功能和质量都影响人类生存环境(气候、水源、土壤、空气等）。在 林业的森林经营与管理中，森林抚育间伐是改善森林结构、促进森林生长、强化森林多种功 能、提高森林质量的最重要的手段。
[0003] 然而，我国的森林抚育间伐技术还很落后，没有得到应有的重视。特别是抚育间伐 设计技术，简单、粗放，多凭生产经验设计，没有效果预测。具体体现在：
[0004] (1)某森林是否有必要开展抚育间伐?间伐效果如何?靠简单判断。
[0005] (2)间伐木是哪些？间伐强度多大?凭经验粗放设计。
[0006] (3)间隔期是多少?间隔期末产量有多大?不能定量预测。
[0007] (4)间伐后郁闭度怎样控制?如何保证间伐后总生产量最大?不能优化计算。
[0008] (5)如何将优化法和经验法相结合，进行人机交互快速优化设计，有效解决上述问 题?亟待发明。
[0009] 与本发明相关的现有技术一是：温带寒温带森林抚育间伐方法，
【发明人】：何玉琦 (1996年）。
[0010] 该发明属于森林抚育间伐技术类，是通过森林分别类型的特点和各种树种的生物 特性，本着培育优质高价林的原则，对由多种类树种组成的整体森林，按幼、中龄林各种不 同的林相、林层结构依据所划分的条件，将森林分成同龄单层林组、同龄近似单层林组、同 龄复层林相、异龄复层林组四个林相组和林相类型、林分起源、林层郁闭度等因子，进行抚 育间伐措施，最终达到森林最终目标林层结构。
[0011] 现有技术一的缺点
[0012] 该发明，只注重分层间伐，人为判断确定间伐对象，调整林分树种组成，没有考虑 如何使间伐后林分总生产量保持最大的问题，这也是本发明要解决的主要问题。另外，该技 术属"人为主观判断"的"定性"抚育间伐，而本发明属"人机结合优化"的"定量"抚育间伐。
[0013] 与本发明相关的现有技术二是:一种天然次生林抚育间伐对象木选择方法(专利 编号：201310559470,发明人：吕勇；刘发林；曾思齐；肖化顺；熊露桥；吕飞舟;刘玉;罗鹏飞， 中南林业科技大学，专利类型:发明专利，申请日：2013.11.12,【公开日】：2014.02.19)。
[0014] 该发明公开了 "一种天然次生林抚育间伐对象木选择方法"，它结合次生林特点及 经营目标，选择影响次生林中林木被确定为抚育间伐对象概率大小的定量指标，如自由度、 混交度、健康指数、目的树种特性指数、空间密度指数和开阔比数。应用层次分析法求算6个 定量指标的权重。根据构建林木的间伐指数ICIi的指标的生物学意义，分析各指标(作为自 变量)与间伐指数(作为因变量)之间的相关关系，以林木的间伐指数ICIi由大至小排序来 确定优先抚育间伐对象木。本发明方法使用方便、操作简单、效率高、精度满足要求。
[0015] 现有技术二的缺点
[0016] 该发明只提供一种如何确定抚育间伐对象的方法(即ICIi指数法），用于确定间伐 哪些林木(或保留哪些林木），同样没有解决"如何抚育间伐才能保证间伐后林分的总生产 量最大"的问题。
[0017] 本发明的理论基础、依据及其缩略语和关键术语定义
[0018] (1)林分
[0019] 在森林中，某一区域地段的林木特征或环境特征明显不同于其相邻地区域段，则 该区域地段就称为是一个林分，其相邻的其它不同的区域地段就是其它不同的林分。一片 森林是由若干林分组成。不同的林分，其特征不同，应该采取不同的经营措施，包括抚育间 伐方案也不同。
[0020] (2)森林抚育间伐
[0021] 指自森林郁闭至成熟之前的幼龄林、中龄林、近熟林等生长阶段，进行定期、定量 而重复伐除部分林木的森林采伐。目的是促进保留林木生长、改善森林结构、强化森林多种 功能。也可称抚育间伐、抚育采伐、抚育伐、间伐等。
[0022] (3)抚育间伐方法
[0023]①下层间伐法
[0024]是指在主林层下层选取树高低、胸径细、生长落后的部分林木做为砍伐对象的一 种抚育间伐方法。目的是消除生长不良林木，改善林内卫生、通风、光照条件，提高保留木生 长空间、促进保留木生长。
[0025]②上层间伐法
[0026]是指在主林层上层选取高大、粗壮、生长优势、抑制其它林木生长的部分林木做为 砍伐对象的一种抚育间伐方法。目的是消除优势霸王木，解放中、下层大部分林木，促进其 生长。
[0027]③综合间伐法
[0028] 是指既在主林层下层选取生长劣势林木、又在上层选取生长优势林木做为砍伐对 象的一种抚育间伐方法。目的是促进中庸木的生长。
[0029] (4)人机交互
[0030]本发明的"人机交互"是指抚育间伐优化方案设计过程中的人机互动。即：抚育间 伐设计者在给定一组"设计方案"约束项的数值后，EXCEL进行优化的"规划求解"，根据 EXCEL "求解"结果，设计者根据专业知识，结合设计经验和生产要求，再改变某一项或某几 项约束值，EXCEL再优化"求解"，设计者通过若干次改变约束项的数值，根据EXCEL若干次 "求解"所得的"优化"方案，从中选择"满意"的抚育间伐设计方案。从而发挥计算机EXCEL的 规划优化功能和设计者的专业实践知识，此过程可称为"人机交互"式抚育间伐设计。
[0031] (5)规划求解
[0032] "规划求解"是系统工程学的量化、优化技术，是OFFICE的EXCELL中的加载项。用于 对用户在开发有关应用模型的基础上，进行优化求解。本发明的方案设计指标的"优化"是 利用EXCEL的"规划求解"加载项功能得以实现的。
【发明内容】
[0033] 本发明针对现有技术的不足，提供一种森林抚育间伐后林分的总生产量最大的间 伐优化设计方法，旨在解决"如何进行优化设计，使间伐后的林分保证结构良好、总生产量 最大"的问题。
[0034] 为解决以上问题，本发明采用的技术方案如下：一种森林抚育间伐后林分的总生 产量最大的间伐优化设计方法，包括以下步骤：
[0035] S1.通过野外样地调查获取林分的径级株数分布特征数据；
[0036] S2.根据S1获取的特征数据，自动生成抚育间伐优化设计所需的样地统计表、株数 分布图和间伐设计表;其中，间伐设计表包括下层间伐法优化方案设计表、综合间伐法优化 方案设计表和上层间伐法优化方案设计表；
[0037] S3.结合S2的各种表项选择至少一种间伐方法，间伐方法包括下层间伐法、综合间 伐法和上层间伐法，利用OFFICE的EXCEL中的规划求解选项，根据模型建立线性规划目标函 数和约束方程，求解拟开展抚育间伐林分的各径级林木的间伐株数；
[0038] S4.如果对EXCEL所得的最优解不满意，再进行人机交互，修改相应的约束项，利用 EXCEL的规划求解加载项重新优化求解，获得满意的优化方案为止。
[0039] 作为优选，S1所述的野外样地调查，采用无边界样圆法，调查工具为林业常用的测 高器、轮尺、围尺和测绳。
[0040] 作为优选，S3所述的模型包括林木冠幅胸径模型：
[0041 ] Cw= 1.1047-0.03113*Ht+0.16662*D
[0042] Cw-林木冠幅直径，Ht-林分优势高，D-林木胸径。
[0043] 作为优选，S3所述的模型还包括林冠覆盖度模型：
[0044] P = N*Cw~2*3.14/4/10000
[0045] P-林分林冠覆盖度，Cw-冠幅直径，10000-公顷平米。
[0046] 作为优选，S3所述的模型还包括间伐蓄积、间伐期末蓄积和总生产量模型：
[0047] (1)间伐蓄积
[0048]为林分间伐木材积之和：
[0049] 帥司=N间* V间=Σ Ni间V i间
[0050] N间-间伐木株数，V间-间伐木平均单株材积，Ni间-第i径级间伐木株数，Vi间-第i径级 间伐木单株材积；
[0051] (2)间伐间隔期末林分蓄积
[0052]为间伐林分间伐间隔期末保留木蓄积之和：
[0053] Μ留=N留 * V留=Σ Ni留V i留
[0054] Nir间伐间隔期末保留木株数，Vir间伐间隔期末保留木平均单株材积，Ni末-伐间 隔期末第i径级间保留木株数，Viir间伐间隔期末第i径级保留木单株材积；
[0055] (3)总生产量
[0056] 为林分间伐蓄积与间伐林分间伐间隔期末林分蓄积之和：
[0057] 啦=司+M留。
[0058] 作为优选，S3所述的线性规划目标函数如下：
[0059] 决策变量为间伐林分各径级的保留株数附留，目标函数为间伐林分间隔期末总生 产量最大：
[0060] max IVfe=Mfsj+M留=N间V间+N留V留=Σ Ni间Vi间+ Σ Ni留Vi留 [0061 ] = ? ((NitrNi留)Vi间)+ Σ (Ni留Vi留）
[0062] Nitr间伐前林分第i径级株数。
[0063] 作为优选，S3所述的线性规划约束条件方程包括：
[0064] (1)各径级保留株数约束之一
[0065] 各径级保留株数要小于或等于相应径级最大保留株数：
[0066] Ni 留彡 Ni 前*Pii〇/100
[0067] Pi大-第i径级最大保留株数强度，由人机交互根据需要确定，以实现更好地控制径 级结构；
[0068] (2)各径级保留株数约束之二
[0069] Ni 留彡 Ni 前*Pi 小/100
[0070] Pi小-第i径级最小保留株数强度，由人机交互根据需要确定，以实现更好地控制径 级结构；
[0071 ] (3)间伐后林分保留总株数约束
[0072] 间伐林分总保留株数要小于或等于间伐间隔期末林分郁闭度1.0的林分株数：
[0073] ΣΜ 留彡
[0074] N?i = 10000/Cw'2/0.5
[0075] Cw = 1.1047-0.03113*(Ht+(DtD后））+0.16662*D末
[0076] Cw-间隔期末冠幅直径，Ht-林分优势平均高，D末-间隔期末林分平均胸径，D后-间伐 后林分平均胸径；
[0077] (4)间伐前后林分平均断面积约束
[0078] 间伐后林分保留木平均断面积要大于或等于间伐前林分林木平均断面积，即间伐 后林分平均胸径要大于等于间伐前林分平均胸径：
[0079] Σ (Ni留*DΓ 2) /Ν?ι 彡 Σ (Ni前*DΓ 2) /%?
[0080 ] N；wi-间伐间隔期末林分郁闭度1.0的林木株数，ΝιΜ_间伐前林分的林木株数；
[0081] (5)间伐株数强度约束
[0082] 间伐株数强度要小于或等于最大间伐株数强度：
[0083] (%总-ΣΝ? 留)/%总*100 彡 P株
[0084] Ρ株-最大间伐株数强度，可根据生产需要和设计技术通过人机交互确定输入值； [0085] (6)间伐材积强度约束
[0086] 间伐材积强度要小于或等于最大间伐材积强度：
[0087] Σ ((NitrNi留)*Vi)/2 (Ni前*Vi)$P材
[0088] Nitr间伐前第i径级株数，Ni留-第i径级保留株数，Vi-第i径级单株材积;P扩最大 间伐材积强度，可根据生产需要和设计技术通过人机交互确定输入值；
[0089] (7)间伐间隔期末林分蓄积量约束
[0090] 间伐间隔期末林分蓄积量要大于等于间伐前蓄积，即保证蓄积量不能越伐越少： [0091 ] V末* ΣΜ留彡 Σ (Ni前*Vi)
[0092] 间伐间隔期末林分平均单株材积，Niir第i径级保留株数，Nifr间伐前第i径级 株数，Vi-第i径级单株材积；
[0093] (8)间伐后郁闭度约束
[0094] 间伐后保留木树冠覆盖度要小于等于规定值P:
[0095] Σ(Ni留*Cwi~2*3.14159/4)/10000彡P
[0096] Cwi-第i径级冠幅直径，P-间伐后保留木树冠覆盖度限制值。
[0097] 作为优选，S4所述的人机交互修改约束各项数值时，要根据如下约束原则进行修 改：
[0098] 株数强度<约束株数强度保留株数 < 最大保留株数
[0099] 材积强度 < 约束材积强度保留株数多最小保留株数 [0100]伐后冠覆<约束伐后冠覆期末蓄积彡伐前蓄积
[0101] 保留株数 < 期末疏1.0株数伐后胸径多= < 伐前胸径
[0102] 在以上原则约束下，适当增加或减小约束项的数值，修改后再进行规划求解，即得 到修改方案参数后的优化方案。
[0103] 本发明的有益效果：本发明根据林业经营单位的森林现状调查数据（即小班调 查），可进行林业生产中的抚育间伐方案优化设计，利用OFFICE的EXCEL中的"规划求解"选 项，根据本发明提出的一系列模型，建立线性规划目标函数和约束方程，求解拟开展抚育间 伐林分的各径级林木的间伐株数，并结合生产经验，通过人机交互，修正方案，使得间伐后 林分结构最好、间伐期末的林分总生产量最大。据此发明得到的优化设计方案，用于开展森 林抚育间伐生产，能使森林的结构最好、产量最大，从而最大限度地发挥森林的生态功能和 经济效益。经多方案模拟研究结果表明：与现行的抚育间伐设计技术相比，本发明的森林抚 育间伐方案优化设计的设计效率可提高50 %以上，使森林总生产量提高10 %以上。
【附图说明】
[0104] 图1为本发明的流程方法示意图；
[0105] 图2为间伐前后各树种径级株树分布图；
[0106] 图3为无边界样圆法示意图。
【具体实施方式】
[0107] 为使本发明的目的、技术方案及优点更加清楚明白，以下参照附图并举实施例，对 本发明做进一步详细说明。
[0108] S1.抚育间伐野外样地调查，若要进行抚育间伐方案设计，必须先开展抚育间伐林 分的野外样地调查，获取设计林分的胸径(径级)株数分布等特征数据。
[0109] 具体操作如下：
[0110] a)打印野外调查表(表1、表2)，用于抚育间伐方案设计的野外样地调查。
[0111] b)野外样地调查，采用无边界样圆法（图3)。调查工具为林业常用的测高器、轮尺、 围尺、测绳等。
[0112] 图3说明：
[0113] 1)林分边界(虚线)与样圆分布(十字交叉）
[0114] 2)样圆半径3~5米，密度大，半径小。
[0115] 3)样圆间距5~10米，林分大，间距大。
[0116] 4)在样圆内每木检尺胸径和5株优势树高，填入样地调查表。
[0117] c)将调查结果分别填入表1和表2。
[0118] d)室内设计前再野外调查表内的数据填入EXCEL电子表格中的表1和表2。
[0119] 注:如使用iPAD，可直接在野外将调查数据录入EXCEL表格中的表1和表2。
[0120] S2.获取"样地统计表、株数分布图、间伐设计表"
[0121 ]根据填写的EXCEL表1和表2，自动生成抚育间伐优化设计所需的"样地统计表、株 数分布图、间伐设计表"。其中，"间伐设计表"包括"下层间伐法优化方案设计表、综合间伐 法优化方案设计表、上层间伐法优化方案设计表"（EXCEL表3~表6、图2)。
[0122] S3.选择间伐方法，输入表中所需的各项约束值
[0123] ①选择间伐方法
[0124] 选择间伐方法是通过"人机交互"，选择"满意"间伐设计方案的基础。
[0125] 设计者结合林分特征（表1~表3)、生产实际和设计经验选择间伐方法（表4~表 6)，可以选择其中之一进行间伐设计，以寻求"满意"间伐设计方案;也可选择其中之二或三 种方法进行间伐设计，通过"人机交互"比较，选择其中的一个"满意"间伐设计方案提交上 级部门批准、交付抚育间伐生产使用。
[0126] ②输入表中所需的各项约束值
[0127] 在表4或表5或表6的"约束Γ和"约束2"各项中输入相应的约束要求数值，就是开 始"人机交互"间伐设计。具体方法和要求是：
[0128] a)如"约束Γ中，在"落松最大保留％"的"6"径级单元格内输入"50"（即落松这个 树种6cm径级的间伐后最大保留株数强度约束限制数值为50%，亦即最多可以保留50%的 6cm径级的落叶松株数），"落松最小保留％"的"6"径级单元格内输入"0"（即落松这个树种 6cm径级的间伐后最小保留株数强度约束限制数值为0%，亦即可以全部间伐砍掉）；
[0129] b)再如"约束2"中，在"期末径级"中输入"16"（即抚育间伐间隔期末预期想要达到 的径级厘米约束数值为16cm) 〇
[0130] c)其中，表4~表6的约束2中"伐后冠覆"的约束输入值为0.5~1.5(伐后郁闭度越 小，"伐后冠覆"的约束输入值越小）；表4 "下层间伐法优化方案设计表"中的"约束2"中的 "调整系数"约束值输入范围为1.1~1.5(密度越大、被压木越多的林分，其"调整系数"约束 值越大）。其它如举例。
[0131] d)实际设计中，如果生产方面没有特殊要求，各约束值可不需人为输入，使用系统 提供的默认值即可开展优化求解，如果对系统所得的最优解不满意，再进行"人机交互"，根 据经验和需求修改相相应的约束项，然后重新优化求解。
[0132] S4.间伐设计方案优化"规划求解"
[0133] 输入并确认上述各约束项的具体数值后，便可实施"规划求解"，即获取间伐设计 优化方案。方法如下：
[0134] a) "规划求解"有解情况下的"人机交互"获得满意设计方案
[0135] 选择EXCEL菜单栏"数据"中的"规划求解"选项卡（在屏幕的右上角），再点击出现 选项卡窗口中的"求解"选项。当EXCEL求解后弹出的窗口中有"规划求解找到一解，可满足 所有的约束和最优状况"提示时，说明已求出相应约束条件下的间伐设计最优指标值，点击 "确定"返回到设计窗口，可查看优化设计结果，如满意(Y)就确定该设计方案为优化的满意 方案。否则(N)，通过"人机交互"修改各项"约束"数值，再"规划求解"，寻求其它的满意设计 方案。
[0136] b) "规划求解"无解情况下的"人机交互"获得满意设计方案
[0137] 当EXCEL "规划求解"后弹出的窗口中有"规划求解找不到有用的解，规划求解找不 到满足所有约束的点"提示时，说明不能求出相应给定约束条件下的间伐设计最优指标值， 点击"确定"返回到设计窗口，可查看求解结果，通过"人机交互"修改各项"约束"数值，再 "规划求解"，寻求其它的"可求解"的满意设计方案(下一步）。
[0138] c)如何"人机交互"修改"约束"各项数值？
[0139] "人机交互"修改"约束"各项数值时，要根据如下约束原则进行修改：
[0141] 在以上原则约束下，适当"增加"或"减小"约束项的"数值"，修改后再进行"规划求 解"，可得到修改方案参数后的"优化方案"，根据专业知识和生产要求，比较各"优化方案"， 从中选择"满意方案"为抚育间伐生产设计方案。
[0142] 其中，要特别注意"伐后冠覆"和"期末径级"两个约束项的"增加"或"减小"。这两 个约束项的数值对方案优化的"规划求解"的"有解"或"无解"影响较大，二者的约束值如果 "太小"，可能造成"无解"，这就需要设计者根据生产要求和专业知识进行"人机交互"具体 确定，以体现科学性、灵活性、适用性。
[0143] 此外，注意"间伐效果"中的"增产％ "值，如果出现负值(说明所产生的方案不增 产，反而减产），则需调整有关"约束"值，如:可增加"期末胸径"等约束项的值，再进行"规划 求解"，求得其它的优化方案。
[0144] 本发明只要确定各项约束值后，"规划求解"很迅速。
[0145] 兴安落叶松抚育间伐方案优化设计技术指标主要模型说明：
[0146] 1、林木分级标准
[0147] 采用林分竞争动态林木分级法对全林分林木进行分级，以确定自然稀疏对象和间 伐采伐对象。
[0148] 定义:林木直径或径级为d，林分平均胸径为D，R = d/D。则
[0149] I级（优势木)R3 1.336
[0150] Π 级（亚优势木）1.336>R会 1.026
[0151] ΙΠ 级（中等木）1.026>R 会 0.712
[0152] IV 级（被压木)0.712>R3 0.383
[0153] V 级(濒死木)R〈0.383
[0154] (据王立明，1989)
[0155] 2、单木材积与林分蓄积模型
[0156] v = 〇 · 0008701*D~2-0 · 0060536*D+0 · 0106236
[0157] v-单株材积，D-林木胸径
[0158] (据内蒙古大兴安岭落叶松一元材积表）
[0159] M = N*v
[0160] Μ-林分蓄积，N-单位林木株数，v-平均单株材积
[0161] 3、林木冠幅胸径模型
[0162] Cw= 1.1047-0.03113*Ht+0.16662*D
[0163] Cw-林木冠幅直径，Ht-林分优势高，D-林木胸径
[0164] 4、林冠覆盖度模型
[0165] P = N*Cw~2*3.14/4/10000
[0166] P-林分林冠覆盖度(简称"冠覆"），Cw_冠幅直径，10000-公顷平米
[0167] 5、伐前和间伐间隔期末的疏密度1.0株数密度模型
[0168] 疏密度1.0株数密度N:
[0169] N=184007*D'-1.7028
[0170] R = 0.9998
[0171] D-林分平均胸径(为伐前林分平均胸径和人机交互选择约束输入的"期末径级"期 望目标值），D={8,22}。
[0172] (据林业部调查设计局，森林调查内业资料汇编，中国林业出版社，1956)
[0173] 6、间伐间隔期末郁闭度1.0株数密度模型 [0174]间伐间隔期末郁闭度1.0株数：
[0175] N獅= 10000/Cw~2/0.5
[0176] Cw-间隔期末冠幅直径
[0177] Cw= 1 · 1047-0 · 03113*(Ht+(DtD后））+0 · 16662*D末
[0178] Ht-林分优势平均高，D末-间隔期末林分平均胸径(为人机交互选择"约束输入"的 "期末径级"目标值），D后-间伐后林分平均胸径（由间伐对象径级和株数计算得出）
[0179](据内蒙古林学院等，1990;王立明，1996)
[0180] 7、间伐间隔期末林冠覆盖度模型
[0181] 间伐间隔期末林冠覆盖度：
[0182] P末=N*Cw~2*3.14159/4/10000
[0183] Cw-间伐间隔期末平均冠幅直径
[0184] Cw= 1 · 1047-0 · 03113*(Ht+(DtD后））+0 · 16662*D末
[0185] Ht-林分优势平均高，D末-间隔期末林分平均胸径(为人机交互选择"约束输入"的 "期末径级"目标值），D后-间伐后林分平均胸径（由间伐对象径级和株数计算得出）
[0186] (据内蒙古林学院等，1990;王立明，1996)
[0187] 8、间伐蓄积、间伐期末蓄积、总生产量模型
[0188] (1)间伐蓄积
[0189] 为林分间伐木材积之和：
[0190] M|司=N间 *V间=ΣΝ?间 Vi间
[0191] N间-间伐木株数，V间-间伐木平均单株材积。-第i径级间伐木株数，-第i径级 间伐木单株材积。
[0192] (2)间伐间隔期末林分蓄积
[0193] 为间伐林分间伐间隔期末保留木蓄积之和：
[0194] Μ留=N留*V留=ΣΝ溜V溜
[0195] Nir间伐间隔期末保留木株数，Vir间伐间隔期末保留木平均单株材积。化末-伐间 隔期末第i径级间保留木株数，V溜-间伐间隔期末第i径级保留木单株材积。
[0196] (3)总生产量
[0197] 为林分间伐蓄积与间伐林分间伐间隔期末林分蓄积之和：
[0198] Μ,^=Μ|'β]+Μ^
[0199] 9、线性规划目标函数
[0200] 决策变量为间伐林分各径级的保留株数Nlgi，目标函数为间伐林分间隔期末总生 产量最大：
[0201 ] max M^=Mfg]+M留=N间V间+N留V留=Σ Ni间Vi间+ Σ N溜V溜
[0202] = Σ ((NitrN溜）Vi间）+ Σ (N溜V溜）
[0203] Nltr间伐前林分第i径级株数。
[0204] 10、线性规划约束条件方程
[0205] (1)各径级保留株数约束之一
[0206]各径级保留株数要小于或等于相应径级最大保留株数：
[0207] N 溜彡 Ni 前*Ρ?〇/100
[0208] Pi大-第i径级最大保留株数强度（％ )，由人-机交互根据需要确定，以实现更好地 控制径级结构。
[0209] (2)各径级保留株数约束之二
[0210] N 溜彡 Ni 前*Pi 小/100
[0211] P,小-第i径级最小保留株数强度（％ )，由人-机交互根据需要确定，以实现更好地 控制径级结构。
[0212] (3)间伐后林分保留总株数约束
[0213] 间伐林分总保留株数要小于或等于间伐间隔期末林分郁闭度1.0的林分株数：
[0214] ΣΝ溜
[0215] N獅= 10000/Cw~2/0.5
[0216] Cw= 1 · 1047-0 · 03113*(Ht+(DtD后））+0 · 16662*D末
[0217] Cw-间隔期末冠幅直径，Ht-林分优势平均高，D末-间隔期末林分平均胸径(为人机 交互选择"约束输入"的"期末径级"目标值），D后-间伐后林分平均胸径（由间伐对象径级和 株数计算得出）
[0218] (4)间伐前后林分平均断面积约束
[0219] 间伐后林分保留木平均断面积要大于或等于间伐前林分林木平均断面积，即间伐 后林分平均胸径要大于等于间伐前林分平均胸径：
[0220] Σ (N溜*Di~2)/Ν?!彡 Σ (Ni前*Di~2)/%?
[0221 ] N；wi-间伐间隔期末林分郁闭度1.0的林木株数，ΝιΜ_间伐前林分的林木株数。
[0222] (5)间伐株数强度约束
[0223] 间伐株数强度（％ )要小于或等于最大间伐株数强度：
[0224] (%总-ΣΝ溜)/_奸100 彡 P株
[0225] P株-最大间伐株数强度（％ )。可根据生产需要和设计技术通过人-机交互确定输入 值(默认为1〇〇%)。
[0226] (6)间伐材积强度约束
[0227] 间伐材积强度（％ )要小于或等于最大间伐材积强度：
[0228] Σ ((NitrN溜）*Vi)/X (Ni前*ViKP材
[0229] Nltr间伐前第i径级株数，^留-第i径级保留株数，Vi-第i径级单株材积。P?-最大间 伐材积强度（％ )，可根据生产需要和设计技术通过人-机交互确定输入值(默认为100% )。 [0230] (7)间伐间隔期末林分蓄积量约束
[0231] 间伐间隔期末林分蓄积量要大于等于间伐前蓄积，即保证蓄积量不能越伐越少：
[0232] V末* Σ N溜彡 Σ (Ni前*Vi)
[0233] V末-间伐间隔期末林分平均单株材积(根据人-机交互给定的"期末径级"值，由单 木材积模型计算），N溜-第i径级保留株数，Nltr间伐前第i径级株数，^-第i径级单株材积。 [0234] (8)间伐后郁闭度约束
[0235] 间伐后保留木树冠覆盖度要小于等于规定值P:
[0236] Σ (N溜*Cwi~2*3.14159/4)/10000彡P
[0237] CWl-第i径级冠幅直径，P-间伐后保留木树冠覆盖度限制值（由人-机交互给定）。 [0238](以上目标函数和约束条件方程为本设计方法发明而设计）
[0239]以兴安落叶松抚育间伐优化设计为实施例：
[0240] S1.抚育间伐野外样地调查，根据调查，填写表1和表2。
[0241 ]表1为样地概况调查表，需选设样地后首先调查，获得数据，用于"抚育间伐方案" 设计的其它数据(表3~表7)统计与计算，以及写入"设计书文本"之需要。
[0242] 表2为开展抚育间伐林分的样地每木调查表，需采用无边界样圆法(附件2)调查后 填写。此表用于表3的统计，以便开展间伐设计。
[0243] S2.获取"样地统计表、株数分布图、间伐设计表"
[0244]完成第一步（即表1和表2填写完成)后，系统会自动生成表3 (样地统计表）、图2 (径 级株数分布图）、表4~表6(下层间伐设计表、综合间伐设计表、上层间伐设计表）。其中，表3 用于了表4~表6的伐前密度统计，图2可用于确定"约束Γ的各项约束数值（8卩"人机交互" 时参考）。
[0245] S3.选择间伐方法，输入表中所需的各项约束值
[0246] 完成第二步后，便可根据林分特征和生产要求选择间伐方法。供选择的间伐方法 有：
[0247] a)下层间伐法
[0248] 在"表4下层间伐优化方案设计"中输入表中所需的各项约束值。
[0249] b)综合间伐法
[0250]在"表5下层间伐优化方案设计"中输入表中所需的各项约束值。
[0251] c)上层间伐法
[0252] 在"表6下层间伐优化方案设计"中输入表中所需的各项约束值。
[0253] 以下各间伐方法举例中：
[0254] "约束1"中的各约束项的输入值："最大保留％"的值都是"1〇〇"（即保留株数最大 可100%保留），"最小保留％"的值都是"Ο"（即保留株数最小可为0%，即全部伐除）。
[0255] "约束2"中的各约束项的输入值："材积强度"约束和"株数强度"约束均为"100" (BP100%可以伐除），"伐后冠覆"约束为"1"，"期末径级"约束为"16"，即间伐间隔期末（BP 下次再开始间伐时）的林分平均胸径要达到16cm(间伐前为11 .lcm，间伐后为9.9cm)〇
[0256] 以上各约束项的约束值，设计者都可根据林分特征和生产要求及技术规程，通过 "人机交互"进行改变，以便寻求"满意"的森林抚育间伐优化方案。
[0257] ①下层间伐法优化方案设计
[0258] 表4为"下层间伐优化方案设计"表。设计者可根据自动统计得出的伐前林分各树 种、各径级的株数密度及其最大密度和林木分级（用于表明林分竞争状况），结合专业知识 和生产要求(保留树种、保留径级等要求），输入"约束Γ和"约束2"的各项指标约束值。之 后，选择EXCEL的"数据"菜单的"规划求解"选项，得到"优化结果"。
[0259] ②综合间伐法优化方案设计
[0260]表5为"综合间伐优化方案设计"表。设计者可根据自动统计得出的伐前林分各树 种、各径级的株数密度及其最大密度和林木分级（用于表明林分竞争状况），结合专业知识 和生产要求(保留树种、保留径级等要求），输入"约束Γ和"约束2"的各项指标约束值。之 后，选择EXCEL的"数据"菜单的"规划求解"选项，得到"优化结果"。
[0261 ]③上层间伐法优化方案设计
[0262] 表6为"上层间伐优化方案设计"表。设计者可根据自动统计得出的伐前林分各树 种、各径级的株数密度及其最大密度和林木分级（用于表明林分竞争状况），结合专业知识 和生产要求(保留树种、保留径级等要求），输入"约束Γ和"约束2"的各项指标约束值。之 后，选择EXCEL的"数据"菜单的"规划求解"选项，得到"优化结果"。
[0263] S4.间伐效果对比及其优化方案决策
[0264] 表7是根据表4~表6自动产生的"间伐效果比较"表。用于选择、决策"间伐优化方 案"。分析表7,根据专业知识、间伐技术规程和生产实际要求，可选择间伐方法和间伐方案。
[0265] 本例中，具有专业知识和设计经验的设计者明显要选择"上层间伐法"，其方案的 具体设计指标优化结果为：间伐前林分胸径11. lcm，密度2046株/公顷，蓄积111.位方米/ 公顷，树种组成7落3杨桦；间伐后林分胸径9.9cm，密度1639株/公顷，蓄积62.6立方米/公 顷，树种组成8落2杨桦；间伐间隔期末林分胸径16cm(即当间伐林分生长到16cm时，可开始 下一次抚育间伐），间伐间隔期为6.6cm(以胸径生长代表年龄增长），间伐间隔期末林分密 度1639株/公顷，蓄积223.73立方米/公顷，间伐蓄积量48.7立方米/公顷，总生产量272.位 方米/公顷，树种组成8落2杨桦;该方案的总生产量比未开展抚育间伐的林分总生产量增长 23.5%，其总生产量的增长率是其它两种方法的1.5~2.5倍，间伐间隔期末蓄积量是间伐 后蓄积量的3.57倍。间伐前冠覆率为1.01，郁闭度为0.9,间伐后冠覆率为0.68,郁闭度为 0.68。株数间伐强度(株强)为20%，蓄积间伐强度(蓄强)为44%。
[0266] 实际中，设计者也可在本举例选择的"上层间伐法"中，输入"约束Γ中的"6径级" 的"最大保留％"和"最小保留％"的值均为"〇"，即将各树种最小的6径级的被压木也全部伐 除，这样就更加科学。即间伐方案可以确定为:在"上层间伐法"为主的条件控制下，也适当 伐除下层将要被自然稀疏掉的被压木(濒死木）。
[0267] 通过此抚育间伐优化方案，可以看出：树种组成得到调整（由原来的7落3白调整为 8落2杨桦），伐除了上层霸王木和下层的濒死木，使林内的光照和通风环境改善，必将提高 林下植物多样性，并扩大了保留木的地上和地下生长的营养空间，促进了全林分生长，增强 了森林的健康性。实现了生态效益和经济效益双赢。
[0268] 表1到表7如下：
【主权项】
1. 一种森林抚育间伐后林分的总生产量最大的间伐优化设计方法，其特征在于，包括 以下步骤：51. 通过野外样地调查获取林分的径级株数分布特征数据；52. 根据S1获取的特征数据，自动生成抚育间伐优化设计所需的样地统计表、株数分布 图和间伐设计表;其中，间伐设计表包括下层间伐法优化方案设计表、综合间伐法优化方案 设计表和上层间伐法优化方案设计表；53. 结合S2的各种表项选择至少一种间伐方法，间伐方法包括下层间伐法、综合间伐法 和上层间伐法，利用OFFICE的EXCEL中的规划求解选项，根据模型建立线性规划目标函数和 约束方程，求解拟开展抚育间伐林分的各径级林木的间伐株数；54. 如果对EXCEL所得的最优解不满意，再进行人机交互，修改相应的约束项，利用 EXCEL的规划求解加载项重新优化求解，获得满意的优化方案为止。2. 根据权利要求1所述的森林抚育间伐后林分的总生产量最大的间伐优化设计方法， 其特征在于，S1所述的野外样地调查，采用无边界样圆法，调查工具为林业常用的测高器、 轮尺、围尺和测绳。3. 根据权利要求1或2所述的森林抚育间伐后林分的总生产量最大的间伐优化设计方 法，其特征在于，S3所述的模型包括林木冠幅胸径模型： Cw= 1.1047-0.03113*Ht+0.16662*D Cw-林木冠幅直径，Ht-林分优势高，D-林木胸径。4. 根据权利要求3所述的森林抚育间伐后林分的总生产量最大的间伐优化设计方法， 其特征在于，S3所述的模型还包括林冠覆盖度模型： P = N*Cw~2*3.14/4/10000 P-林分林冠覆盖度，Cw-冠幅直径，10000-公顷平米。5. 根据权利要求4所述的森林抚育间伐后林分的总生产量最大的间伐优化设计方法， 其特征在于，S3所述的模型还包括间伐蓄积、间伐期末蓄积和总生产量模型： (1) 间伐蓄积 为林分间伐木材积之和： Μι司=N间* V间=Σ N i间V i间 N间-间伐木株数，V间-间伐木平均单株材积，Ni间-第i径级间伐木株数，Vi间-第i径级间伐 木单株材积； (2) 间伐间隔期末林分蓄积 为间伐林分间伐间隔期末保留木蓄积之和： Μ留=N留*V留=ΣΝ?留Vi留 N留-间伐间隔期末保留木株数，Vir间伐间隔期末保留木平均单株材积，Ni末-伐间隔期 末第i径级间保留木株数，Viir间伐间隔期末第i径级保留木单株材积； (3) 总生产量 为林分间伐蓄积与间伐林分间伐间隔期末林分蓄积之和： fe=M|司+M留。6. 根据权利要求5所述的森林抚育间伐后林分的总生产量最大的间伐优化设计方法， 其特征在于，S3所述的线性规划目标函数如下： 决策变量为间伐林分各径级的保留株数Ni留，目标函数为间伐林分间隔期末总生产量 最大： max =Mfg]+M留=Nfg^、g]+N留V留=Σ N i 间V i 间+ Σ N i 留V i 留 =Σ ((NitrNi留)Vi间)+Σ (Ni留Vi留） Nitr间伐前林分第i径级株数。7.根据权利要求6所述的森林抚育间伐后林分的总生产量最大的间伐优化设计方法， 其特征在于，S3所述的线性规划约束条件方程包括： (1) 各径级保留株数约束之一 各径级保留株数要小于或等于相应径级最大保留株数： Ni 留彡 Ni 前*PiiyiOO Pi大-第i径级最大保留株数强度，由人机交互根据需要确定，以实现更好地控制径级结 构； (2) 各径级保留株数约束之二 Ni留彡Ni前*Pi小/100 Pi小-第i径级最小保留株数强度，由人机交互根据需要确定，以实现更好地控制径级结 构； (3) 间伐后林分保留总株数约束 间伐林分总保留株数要小于或等于间伐间隔期末林分郁闭度1.0的林分株数： XNigSN細 N?i = 10000/Cw'2/0.5 Cw= 1.1047-0.03113*(Ht+(DtD后））+0.16662*D末 Cw-间隔期末冠幅直径，Ht-林分优势平均高，D末-间隔期末林分平均胸径，D后-间伐后林 分平均胸径； (4) 间伐前后林分平均断面积约束 间伐后林分保留木平均断面积要大于或等于间伐前林分林木平均断面积，即间伐后林 分平均胸径要大于等于间伐前林分平均胸径： Σ (Ni留*Di~2)/N末郁1 彡 Σ (Ni前*Di~2)/N·、 间伐间隔期末林分郁闭度1.0的林木株数，ΝιΜ_间伐前林分的林木株数； (5) 间伐株数强度约束 间伐株数强度要小于或等于最大间伐株数强度： (%总-Σ Ni留)/NiM* 100彡卩株 P株-最大间伐株数强度，可根据生产需要和设计技术通过人机交互确定输入值； (6) 间伐材积强度约束 间伐材积强度要小于或等于最大间伐材积强度： Σ ((NitrNi留）*Vi)/Σ (Ni前*Vi KP材 Ni前-间伐前第i径级株数，Ni留-第i径级保留株数，Vi-第i径级单株材积;P?-最大间伐 材积强度，可根据生产需要和设计技术通过人机交互确定输入值； (7) 间伐间隔期末林分蓄积量约束 间伐间隔期末林分蓄积量要大于等于间伐前蓄积，即保证蓄积量不能越伐越少： V末*ΣΝ?留彡 Σ (Ni前*Vi) V末-间伐间隔期末林分平均单株材积，Niir第i径级保留株数，Ni前-间伐前第i径级株 数，Vi-第i径级单株材积； (8)间伐后郁闭度约束 间伐后保留木树冠覆盖度要小于等于规定值P: Σ(Ni留*Cwi~2*3·14159/4)/10000彡P Cwi-第i径级冠幅直径，P-间伐后保留木树冠覆盖度限制值。8.根据权利要求1所述的森林抚育间伐后林分的总生产量最大的间伐优化设计方法， 其特征在于，S4所述的人机交互修改约束各项数值时，要根据如下约束原则进行修改： 株数强度<约束株数强度保留株数 < 最大保留株数 材积强度 < 约束材积强度保留株数多最小保留株数 伐后冠覆<约束伐后冠覆期末蓄积多伐前蓄积 保留株数 < 期末疏1.0株数伐后胸径多= < 伐前胸径 在以上原则约束下，适当增加或减小约束项的数值，修改后再进行规划求解，即得到修 改方案参数后的优化方案。
【文档编号】G06F19/00GK106096305SQ201610465204
【公开日】2016年11月9日
【申请日】2016年6月22日
【发明人】张秋良, 王立明, 王飞, 张秀丽, 王冰, 韩胜利, 铁牛
【申请人】张秋良, 王立明