研究林分位于北京市门头沟区东南部的九龙山，属北京西部山系(115°59’~116°07’E，39°54’~39°59’N)，海拔130 m，坡度20°。大陆性季风气候，年平均日照时数2 621.1h，年平均气温11.7℃，全年平均降水650.4 mm，主要集中在7、8、9月。土壤以山地棕壤为主。试验林分植于20世纪70年代。林分具体特征见表 1。

在京西山地人工侧柏试验地内设置1个大样地(图 1)，分为24个面积为15 m×12 m=180 m²的小区，依次编号(1~24号)。该林分计划施行包含4个处理的森林经营试验，分别为结构化森林经营、结构化森林经营与割灌、割灌和对照。利用完全随机优化设计方法将小区分为4组，每组6个重复。利用全站仪(Topcon_GTS600)对各小区的林木进行每木定位，测量每木的胸径、树高，并记录树种及其健康状况。各小区参数见表 2。本研究中涉及到的林分空间结构参数计算和数据处理以及模拟抽取小区重复，均使用R语言编程软件完成。

图  1  大样地各小区布置及编号示意图

Figure 1.  Arrangement and plot numbers in the sample stand

当改变小区隔离带距离或利用点抽样方法分析数据时，各小区林分状态也相应改变，因此，需要重新确定最佳隔离带距离及相应的小区重复个数、点抽样林木株数和最佳抽样点个数。因此, 本研究的数据处理方法主要分为2部分：(1)小区数据为全面调查数据，逐渐增加各小区间的隔离带距离，缩小每个小区的调查面积，找到最佳的隔离带大小；(2)只调查距离各小区中心点最近n(n=4~12)株的林木，找到最佳调查林木株数。这2种方法都涉及到小区重复，而不同小区作为重复时，得到的结果存在差异，为了避免这种不确定性，在研究小区重复时，采取了随机抽取小区的办法, 也就是在24个小区中随机抽取m个作为重复，计算各参数的均值并与样地真值比较，这个过程进行1 000次，统计这1 000次中满足条件的次数，当有90%或以上概率满足要求时，则认为结果是可靠的。

在本研究中将样地中小区设置为15 m×12 m，对调查得到的数据进行分析时忽略各小区边缘一定范围内的林木，划定隔离带，增加小区间隔距离，以减小相邻小区间距离较近的林木相互影响。为了研究小区间隔离带的最适距离，依次增大小区间的间隔距离(图 2)。分别设置为1、2、3、4 m，小区的调查面积相应减小为13 m×10 m、11 m×8 m、9 m×6 m、7 m×4 m。计算所有24个小区不同隔离距离下的基本参数(胸径D、树高H)与空间结构参数(角尺度W、大小比数U、混交度M)均值并与大样地进行比较, 其中，角尺度W的相似性以相同分布格局为准：当W>0.517时，林分为团状分布；W <0.475时为均匀分布；W处于[0.475，0.517]时为随机分布。混交度M为[0，0.25)，[0.25，0.75]，(0.75，1]时，分别为轻度混交、中度混交和强度混交标准，其它参数以相差10%以内为标准。在进行林分状态比较时，应用以上所有参数。

图  2  小区间隔离距离和调查面积示意图(图中灰色部分为小区隔离带，白色区域为进行数据与实验的核心区域)

Figure 2.  The distance of isolating neighboring plots (grey part are isolated area and the white part are core area for trial)

当固定隔离带距离时，用计算机模拟的方法，随机在24个小区中选择n(1~23)个小区作为重复，计算不同重复小区个数的林分基本参数和空间结构参数并与真值进行比较。当n取值不同时，分别进行1 000次随机模拟，以期找到在不同隔离带大小时最少的小区重复次数。

设置每个小区的中心点为抽样点，计算距离抽样点最近n(n=4，5……12)株林木的基本参数与空间结构参数的平均值并与大样地林分状态比较。图 3为调查n=4株木为例的调查方法，以调查林木不重复为原则，黑色方块为小区中心点。在该样地中, 横向为每隔12 m、纵向每隔15 m都有一个抽样点，其他株数的调查方法雷同。

图  3  调查距离中心点最近n株林木示意图(以n=4为例)

Figure 3.  Method to investigate n trees around the middle points (for example as n=4)

当随机抽取n(1~23)个小区时，为确定与样地林分状态基本一致的重复小区数，首先在24个小区中随机抽取任意n个小区(无重复)，计算这n个小区各参数平均值并与大样地林分状态比较。调查不同株林木时各重复1 000次，以期找到最佳的小区重复次数。

当隔离带由1 m增加到4 m时，小区的调查面积由130 m²逐渐减小至28 m²，得到表 3中24小区基本参数与空间结构参数均值。与大样地参数相差不大，说明当调查区域达到一定面积时，与大样地的基本参数、空间结构参数有一定相似性，尤其当隔离带为2 m时最相近。

由表 3中数据与大样地数据得到的相对误差见图 4。在图 4中，当小区间的隔离带设置为2 m或3 m时，基本生长参数和空间结构参数与大样地的相对误差均小于5%，其中，当隔离带为2 m时，也就是当小区面积为88 m²时的相对误差最小。因此，进行全面调查时，小区间的隔离带应设置为2 m。即相邻小区间的距离为4 m。

图  4  改变小区隔离带距离时基本生长参数和空间结构参数与大样地的相对误差

Figure 4.  The relative error between plots and stand of basic parameters and spatial

当隔离带大小不同时，需要不同个数的重复。当隔离带越小、核心调查面积越大时，需要的重复越少。当隔离带为1 m时，任意抽取5个小区作为重复即可达90%以上的概率与样地林分状态基本一致；当隔离带为2 m时也只需要6个重复；隔离带达到最大即4 m时，最多只需要11个重复。当重复个数达到19个或更多时，不论小区间的隔离带多大，均有100%的概率与林分状态相似。表 4给出了不同隔离带时，1 000次模拟随机抽取n个小区作为重复达到与样地林分状态相似的个数。当重复大于19个时，1 000次随机抽取的结果都与大样地林分状态相似的。

当改变每个小区点抽样附近株数时，得到的基本参数与空间结构参数也随之改变。表 5为调查抽样点最近4株至12株林木24个小区均值的变化。将不同抽样结果的参数与大样地参数比较得到其相对误差。当选取的抽样点附近的树木过少时，误差相对较大，可达到40%。当抽样点达到4个时，除了混交度(M)误差为-12.24%，其余基本参数与空间结构参数其误差基本控制在15%以内。当抽样点达到5个时，误差可控制在10%以内，其中，以6株木为佳，其误差均在5%内。当改变小区中心点最近调查林木株数时，得到的基本参数与空间结构参数也随之改变。

随后研究了当调查株数不同时需要的重复。随机抽取n(1~23)个小区作为重复，统计模拟1 000次时其小区均值与大样地相似的概率(图 5)。随着调查株数的增多，需要的重复逐渐减少。当调查株数达到6株时，只需要9个小区即可与大样地林分状态相近；当达到11株或12株时，只需要6个小区即可。当重复的小区个数达到22个或更多时，不论每个小区调查的株数多少，都有100%的概率可以满足与大样地参数基本相似的要求。

图  5  距中心点最近林木株数和小区重复次数改变时林分状态与大样地相似的次数

Figure 5.  The number of times of the investigating results similar to the stand when changing the number of neighboring tree around the middle point or repeat times

调查不同株数时，林木代表的面积不同。一般来说，调查株数越多代表的面积越大，调查株数距离中心点最远的1株林木距离边缘也越近。研究统计了调查不同株数时距离各小区边界的距离，从表 6可见：当调查株数逐渐增加时，距小区边界逐渐减小。当调查4株时，第4株林木距小区边缘的平均距离约为4 m，当调查林木增加至10株以上时，最后1株林木到小区边缘的平均距离约为3 m。这从另一个方面反应出小区间间隔的程度，说明2个小区间的隔离带只要小于3即可。在上文中将小区之间的隔离带设置为2 m。

将小区的调查面积设置为88 m²时，其基本参数与空间结构参数与大样地林分状态最为接近，其与大样地的误差均小于2%，其中, 角尺度(W)的相对误差只有-0.14%。将各小区调查面积设置为88 m², 以保证其与大样地林分状态基本相似，小区重复后的各参数与大样地相似度可达到98%以上，且只需要6个重复小区就可达到90%概率与样地林分状态相似。因此，在该林分经营试验中将小区间的隔离带设置为2 m，也就是相邻的小区间距为4 m。

当只对各小区中心点附近林木进行调查比较时发现，选取林木株数过少则误差相应增大，当选取树木达到4株时，小区各参数与大样地林分相似度达85%；当抽样点附近树木达到5个时，小区各参数与大样地林分相似度达90%, 其中，以6株木为佳，小区各参数与大样地林分相似度达95%。当选取6株木为调查对象时，逐渐减少重复小区，开始相对误差变化较小，表明与林分状态仍很接近，直到重复小区减少到11块时，林分基本参数与空间结构参数与大样地误差均保持在5%以内；当重复小区继续减少时，误差开始明显增大，一些参数的相对误差可超过10%。当减少为7块小区时，误差可达到20%以上。若调查林木选取为11株以上时，只需要6个重复即可保证与大样地林分特征一致。

研究首先尝试了当设置小区之间的隔离带分别为1、2、3、4 m时，其各项参数与不设置隔离带的各小区的参数比较。显然，当隔离带逐渐增大、小区调查面积逐渐减小时，与其本身参数相近的小区个数越来越少。当隔离带为1 m时，只有1/3的小区偏离了原来的取值，而当隔离带增大到4 m时，只有1/3的小区没有偏离原来的取值。这一结果说明设置隔离带后并不能准确的表达小区原本的林分状态，而应设置重复。研究这个问题涉及到2个变量：隔离带的最佳距离和最少重复个数。本研究的结果也证实当小区隔离带过小为1 m时误差相对更大，2 m误差最小，重复6次最为适宜；而对距小区中心点最近林木的研究发现，当选取6个小区作为重复时，需要11株或更多林木，而距小区中心点最远的林木到小区边缘的平均距离约为3 m，这也从另一个方面说明，小区间的隔离带不能大于3 m。因此, 对于该样地来说，应将小区间的隔离带设置为2 m，相邻小区间的间距为4 m。

在抽样调查中，分别研究了4株至12株的调查方法。从理论上分析，株数越多，调查覆盖面积越大，精度应越高，但从计算结果来看，7株至12株的效果并不比6株好，这可能是由于样圆过大时林木分布不均，林间空地、林道等多方面对其产生了影响。在调查较少株数时这种情况较少，而增大样圆面积无疑会增加这种情况发生的概率。这个结果与刘凤阁等^[20]的6株调查优于20株调查效果的结论相似。

(1) 将小区的调查面积设置为88 m²时，其基本参数与空间结构参数与大样地林分状态最为接近，其与大样地的误差均小于2%，且只需要6个重复小区就可达到90%概率与样地林分状态相似。因此，在该林分经营试验中可将小区间的隔离带设置为2 m，也就是相邻的小区间间距为4 m。

(2) 在对该样地进行抽样调查时，调查12个小区内距中心点最近的4株木即可满足调查的要求，或6个小区内各调查距中心点最近的11株林木。