Visual Simulation of Thinning of Chinese Fir Based on Stand Spatial Vertical Structure Characteristics

研究区位于湖南省株洲市攸县境内的黄丰桥国有林场，113°04′～113°43′ E，27°06′～27°04′ N，属于亚热带季风湿润气候，年平均气温17.8℃，无霜期达229 d，年均降水量为1 410.8 mm，森林覆盖率为86.24%，树种主要以杉木（Cunninghamia lanceolata（Lamb.）Hook.）、马尾松（Pinus massoniana Lamb.）为主，其中杉木面积3 197.6 hm²，全部为人工造林。

在研究区选取杉木人工纯林典型样地4块（分别命名为1、2、3、4。地位指数：18，坡向：南坡，坡度：5～10°，海拔：270～320 m）。利用全站仪测量得到杉木单木根部相对坐标（X, Y, Z），以常规外业调查方式调查胸径（DBH）、树高（H）、冠幅（CW）、活枝下高（UBH）、生长状态等，按1年时间间隔采集2012—2017年调查数据。其核心区（距边界5 m）属性分布如表1所示。

为了表达对象木与其竞争林木之间的大小差异，引用参考文献[13]李思佳研究提出的垂直空间结构参数$A\_Pv$。

本研究中将文献[13]式中的空间结构单元数n取4，当$A\_P{v_i}$越接近−1时，整个样地的平均竞争状况越好，可以适当减小间伐强度；当$A\_P{v_i}$越接近1时，说明样地的平均竞争状况处于劣势地位，应加大间伐强度来改善样地内的竞争状况。

结合垂直空间结构参数$A\_P{v_i}$和森林经营方案原定间伐强度，经过多次实验认为平均$A\_P{v_i}$对间伐强度的影响系数为0.1较为合理，提出考虑林分垂直结构的间伐强度指标。公式如下：

式中：S为经过改进后的抚育间伐强度；P为森林经营方案中原定间伐强度，本研究中为27%；$A\_P{v_i}$为林分平均垂直结构参数。当样地平均竞争状况良好时，$A\_P{v_i}$为负值，则经过调整后的间伐强度较原间伐强度更小；当样地平均竞争状况处于劣势时，$A\_P{v_i}$为正值，则经过调整后的间伐强度较原间伐强度更大，以期对样地的竞争状况获得更好的调整效果。

以黄丰桥国有林场的森林经营方案为依据，确定4块样地所在小班的间伐强度均为27%，采用机械抽样的方式选取间伐木。

通过上文提到的考虑林分垂直结构的间伐强度指标来对原有的间伐强度进行改善，对样地进行模拟间伐，并与采用27%强度间伐的样地进行对比。

为了检验考虑林分垂直结构的间伐强度与原定间伐强度相比对林分空间结构的改善情况，本研究从林分空间分布格局和林木个体大小差异两个方面来描述林分空间结构，选取了大小比数和角尺度作为表达林分空间结构的参数。

基于Unity3D引擎实现间伐影响下的林分生长动态可视化模拟。主要采用了高亮显示技术、公告板技术、对空间结构单元中树木的信息、距离等进行可视化模拟，能更直观的表达林木的生存空间。采用Unity中最新的UI系统UGUI对林分中大小比数、角尺度、垂直结构参数$A\_P{v_i}$进行可视化显示，采用读取间伐表的方式对样地进行模拟间伐。

使用系统的林分加载功能将样地调查表导入到可视化系统中，通过样地调查表中的树高、胸径、冠幅、枝下高数据，模拟样地真实三维场景。通过读取间伐表的方式模拟采伐，再运行林分结构分析功能，可以直接显示样地的大小比数和角尺度的分布频率，平均大小比数和平均角尺度，以及平均垂直空间结构参数$A\_Pv$，如图1所示（以3号样地为例，下同）。

Figure 1.  3rd sample plot before the thining 3D virtual scene simulation

根据黄丰桥国有林场森林经营方案，对样地进行间伐强度为27%的定量间伐，采用机械抽样的方法选取间伐木，采用系统读取间伐表的方式进行模拟间伐。系统运行后的样地情况和各空间结构参数如图2所示。具体间伐后结果见表2：

Figure 2.  3rd sample plot after 27% intensity thinning 3D virtual scene simulation

可以看到，经过27%间伐强度间伐后的样地，与间伐前样地相比，4块样地的平均$A\_Pv$均出现了较大幅度的增长，样地内单株木与其4棵相邻木的竞争依然处于不平衡的状态；平均角尺度变化不大，依旧维持在0.5左右，处于随机分布状态；在平均大小比数方面，4块样地的平均大小比数与间伐前相比，都出现了较大幅度的上升，说明单株木在与其最近4株相邻木的竞争中普遍处于优势地位。

根据4块样地在间伐前的垂直结构参数A_Pv对原有的间伐强度进行改进，根据(4)式可以计算出改进后的间伐强度及间伐株数，具体结果如表3：

根据表3的间伐强度及间伐株数计算结果，使用林分三维场景可视化模拟系统对4块样地进行模拟间伐，可以得到经过考虑林分垂直结构的间伐强度间伐后的林分三维场景，运行林分结构分析功能，可以得到经过考虑林分垂直结构的间伐强度间伐后的样地内的大小比数、角尺度以及垂直空间结构参数$A\_Pv$（图3），具体间伐后结果如表4：

Figure 3.  3rd sample plot after improved intensity thinng 3D virtual scene simulation

可以看出，经过考虑林分垂直结构的间伐强度间伐后4块样地与采用27%间伐强度间伐后的样地相比，平均$A\_Pv$更加趋近于0，其中2、3和4号样地更加明显，说明这3块样地的单株木与其最近4棵相邻木的竞争中整体上处于均势；4块样地的平均角尺度与间伐前的样地和经过27%间伐强度间伐后的样地相比，均无明显变化，整体均维持在0.5左右，属于随机分布；4块样地中，除了3号样地外，平均大小比数均出现了较大幅度的减少。

本研究使用考虑林分垂直结构的间伐指标与森林经营方案原定间伐强度相比，对$A\_Pv$改善效果明显，使林分中的竞争更加平均。但是对角尺度的改善无明显效果，可能是样地均为人工纯林的原因；而对大小比数的影响无规律可循，可以在今后的研究中加大模拟样地间伐的数量，探寻其规律。本研究使用的林分三维可视化模拟系统对样地真实场景的模拟效果较好，开发的森林结构功能模块对间伐活动及间伐后样地的各空间结构参数分析准确，为森林经营者提供了新的决策手段。

在今后研究中，一方面需要考虑大场景森林可视化模拟优化技术，实现大规模、高细节的场景模拟，以及在林分空间结构分析功能模块中加入聚集指数、开敞度等其他空间结构参数，以直接反映林分的空间结构特征；另一方面应继续改进考虑林分垂直结构的间伐强度，合理考虑更多的变量，如聚集指数、林分拥挤度等空间结构参数，以期相较原定间伐强度而言，对林分空间结构的改善效果更加显著。

本研究以杉木人工林为研究对象，利用林分三维场景可视化模拟系统中的森林结构功能模块，通过读取杉木样地调查表中的每木坐标、胸径、树高、冠幅和枝下高数据，可以在虚拟场景中模拟样地的真实情况，分析样地的大小比数、角尺度及$A\_Pv$的平均值，为研究间伐影响下的林分空间结构变化提供了可视化手段。经过考虑林分垂直结构的间伐强度间伐后的4块样地与原定间伐强度相比，平均$A\_Pv$更加趋近于0，说明经过考虑林分垂直结构的间伐强度间伐后样地内林木竞争状况更加均匀；平均角尺度无大幅度变化，4块样地都维持在0.5左右，呈随机分布状态；4块样地中，除了3号样地外，平均大小比数均出现了较大幅度的减少。