张志才 ,陈亚文

1.石棉县自然资源和规划局，四川 雅安 625400；2.四川省林业和草原调查规划院，四川 成都 610081

地位指数作为立地质量评价与林地生产力的常用指标[1,2]，被国内外学者进行广泛研究[3-7]，北美学者早在20 世纪初期就以地位指数法来评价林地生产力[8]，同时，美国中、西部地区也将地位指数曲线用于指导实际生产[9]。随后，在对加拿大黑松的立地质量进行评价时，环境因子中的温度作为变量被引入立地质量评价模型[10-12]，以便更加准确的反映林地生产力随温度的变化情况。我国从20 世纪80 年代起，各地开始编制马尾松、落叶松、杉木、杨树、桉树[13-16]等主要树种的地位指数表，但多集中在人工林，对天然林的研究相对较少。栎类天然林在中国面积大、分布广，具有极高的经济价值、生态价值，是我国最重要的森林资源之一。已有学者对秦岭、华北、伏牛山等地的栓皮栎天然林[17-20]、湖南地区的栎类次生林[21]地位指数表进行过相关研究，横断山作为中国第一﹑第二阶梯的分界线，海拔落差大，气候垂直变化明显，孕育着丰富的植物群落，本文以横断山东部石棉县的211 块栎类天然林样地数据为基础，编制地位指数表，以期客观、科学、准确的评价横断山东部栎类天然林立地质量，为该地森林可持续经营提供依据。

1 研究区概况

研究区位于青藏高原横断山脉东部的石棉县，地处大渡河中游，东连甘洛县、汉源县，南接冕宁县、越西县，西依康定市、九龙县，北与泸定县毗邻。地理位置介于东经101°55'—102°34'，北纬28°51'—29°32'之间。境内东西跨度达60 km，南北跨度达76.5 km，全县面积2 678.2 km2。石棉县位于川西高原和四川盆地的过渡带，地貌特征为山高谷深，由极高山至中山，由高原到河谷阶地、坝子呈层状地貌发育，最低海拔790 m，最高海拔5 793 m。大渡河由泸定县流入石棉县境，由北向东，贯穿县境中部，流经10 个乡（镇）。石棉县属中纬度亚热带季风气候为基带的山地气候，孕育着丰富多样的植物群系，全县分布有高等植物208 科900 属2 468 种，其中：种子植物156 科818 属2 308 种，蕨类植物29 科51 属123 种，苔藓植物23 科31 属37 种。根据石棉县2020 年林木种质资源普查结果，县域内有壳斗科（Fagaceae）植物5 属，47 种，其中栎属（Quercus）植物20 种，主要包括麻栎（Quercus acutissima）、白栎（Quercus fabri）、高山栎（Quercus semecarpifolia）、栓皮栎（Quercus variabilis）等。

2 研究方法

2.1 数据收集与整理

本研究中用于编制地位指数表的基础数据来源于石棉县2016—2021 年设置的211 块栎类天然林样地数据，样地为25.82×25.82 m 的方形样地，记录样地各项林分指标和单木指标，主要包括海拔、坡度、坡位、坡向、土壤厚度，年龄、胸径、树高等，其中胸径采用每木检尺法，树高测量3 株优势木高取平均值。211 块栎类天然林样地以中（41 块）、幼龄林（92 块）为主，占样地总数的63%，成、过熟林分布相对较少，各样地的优势木平均高分布在5.2～19.9 m 之间（见表1）。以其中140 块样地数据作为地位指数表的建模数据，剩余71 块样地数据作为检验数据。

表1 样地年龄分布统计Tab.1 Statistics of age distribution of the sample plots

2.2 异常数据剔除

用于建模的140 块样地数据，其年龄分布在10～70 年，按照5 年一个龄阶距，共将其划分为13个龄阶，统计出分布于各龄阶样地的优势木平均树高值，并计算出标准差值，以平均树高为基准，采用3 倍标准差法（）剔除异常数据，最终得到134 组数据对（见表2），用于树高导向曲线的拟合[22]。

表2 样地林分特征统计Tab.2 Statistics of stand characteristics of the sample plots

2.3 导向曲线的选择

选取8 个常用方程对导向曲线进行拟合，其中包括对数、指数、二次函数3 个数学方程，坎派兹、考尔夫、理查德、威布尔、逻辑斯蒂5 个生长方程，各方程的表达式见表3。用SPSS 统计软件，以最小二乘法对各方程的参数进行估计，以决定系数（R2）、残差平方和（sum of squares for error，SSE）和预估精度（P）作为模型适用性检验的指标，综合考虑以上3 个指标，选取最优方程拟合导向曲线。

式中：ym为 树高实际值；为树高理论值；为预测值的标准差；t0.05为置信水平α=0.05 时的t分布值。

3 结果与分析

3.1 导向曲线拟合结果

拟合结果显示（见表4），8 个方程的确定系数都在0.960 以上，其中逻辑斯蒂（Logistic）方程的确定系数最大（R2=0.992），并且，其精度（P=98.23%）最高、残差平方和（SSE=0.848）最小，将逻辑斯蒂方程作为导向曲线的拟合公式，即：

表4 曲线方程计算结果统计Tab.4 Statistical analysis of curve calculation results

3.2 基 准年龄（A0）与地位指数级 距（C）的确定

林分的树高生长在基准年龄时应该处于稳定状态，同时不同立地条件的树高生长值应该存在明显的差异，以134 块标准地的树高-年龄数据为基础，通过公式可以计算出不同龄阶时树高标准差产生变化的幅度（ΔS）以及变动系H数值的变化幅度（ΔCH），以年龄为横坐标，分别以ΔSH和ΔCH值绘制变化幅度图。从图1 可以看出，在40 年以前ΔSH的变化范围为0.76～0.96，ΔCH的变化范围为1.31～1.09，二者的变化幅度相对较大，在40 年以后，ΔSH的值稳定在0.97～0.98，ΔCH的值稳定在1.06～1.08 左右，二者的变化幅度逐渐趋于稳定，且其变化值无限接近于1，表明在40 年以后树高的生长趋势比较稳定，这时的年龄值（40 年）可以确定为横断山栎类天然林地位指数表编制的基准年龄[18-19,22]。

图1 树高标准差、变异系数变化幅度Fig.1 Variation range of standard deviation and coefficient of variation of tree height

在基准年龄（A0=40a）时栎类天然林林分树高的变化范围为7.2～17.4m，即ΔH=10.3m，本研究以2m 为一个指数级距，通过公式n=ΔH/C确定n为6，即6 个地位指数级，最终确定为8、10、12、14、16、18 m。

3.3 地位指数表的编制

编制地位指数表常用的方法有3 种：变动系数调整法、标准差调整法以及相对优势高法。针叶树种地位指数表的编制多用变动系数法，阔叶树种的地位指数表编制多用标准差调整法，且标准差调整法相对其余2 种方法来说更为简便好用[17]，因此本研究以标准差调整法来展开横断山栎类天然林地位指数表。

以各龄阶的平均年龄值为自变量，各龄阶树的高标准差值为因变量，利用对数函数SH=a+bln（A）式进行拟合，即可得到各龄阶树高标准差理论方程为：

以导向曲线为基础，根据基准年龄（A0=40a）时的树高值以及地位指数级距（C=2 m），采用阔叶树种常用的地位指数表编制方法标准差调整法展开横断山东部栎类天然林地位指数曲线簇。将各龄阶年龄（A）代入最终选定的导向曲线方程，则可以得到各龄阶导向曲线的树高理论值（Hmk）。

因在基准年龄A0时，通过导向曲线计算出的理论树高值通常与地位指数级的数值不相等，所以应根据H0和S0的大小，采用下式进行相应调整：

式中：Hmn为第m 龄阶第n 指数级调整后的树高；Kn为调整系数；Hmk为第m 龄阶导向曲线的树高；H0k为基准年龄时导向曲线的树高值；H0n为基准年龄时第n 指数级的树高；SAm为第m 龄阶树高标准差理论值；SA0为基准年龄所在龄阶树高标准差理论值。

将各龄阶年龄值（A）代入（5）式，就可以得到各龄阶的树高标准差理论值（SAm）。本研究中，在基准年龄40 年时，导向曲线上树高的理论值为11.40 m，与划定的地位指数有一定差距，而与其最为接近的地位指数级数为H=12 m，则应按照式（6）进行相应的调整。首先计算出基准年龄40 年时导向曲线上的树高值（11.40 m）和指数级12 m 的差值d=H0n-H0k=0.60 m，然后计算调整系数值Kn=（H0n-H0k）/SA0=0.227，最后将计算出的调整系数（Kn）与各龄阶树高标准差理论值（SAm）相乘，就可以得到各龄阶的树高调整值（Kn·SAm）。将各龄阶导向曲线上的树高值（Hmk）与各龄阶调整值相加（减）就可以得到以12 m 指数级为基准的各龄阶调整后的树高值（H0）[22,23]，结果见表5。

表5 不同年龄的树高调整值统计Tab.5 Statistics of tree height adjustment values at different ages

以调整后的导向曲线（12 指数级）为准，按指数级距C=2 m，逐龄阶倒算出各地位指数级曲线上的树高值，其余指数级的调整系数Kn为：

本研究中Kn=2/2.64=0.76，各龄阶内相邻指数级间的调整值为Kn·SAm=0.76SAm，然后按式（4）计算出各龄阶内各地位指数级的树高值，最后展开可得地位指数表和地位指数曲线簇（见表6、图2）。

图2 横断山栎类天然林地位指数曲线簇Fig.2 Site index curve cluster of Quercus natural forest in Hengduan Mountains

表6 横断山栎类天然林地位指数表Tab.6 Site index of Quercus natural forest in Hengduan Mountains

4 地位指数表检验

对于所编制的地位指数表，需进行检验，确保其有较高的精度和适用性，方能投入实际应用。本文以拟合显著性检验和落点检验等常用的检验方法对所编制的横断山栎类天然林地位指数表进行检验。

4.1 树高值 χ2检验

根据71 个检验样地的年龄与树高值，查所编制的横断山栎类天然林地位指数表，可得出每块样地对应的地位指数值，求出各样地的树高理论值，并与实际树高进行 χ2检验，在SPSS 软件中计算得出χ2=0.405 6，远小于表明用地位指数表查的树高理论值与样地实测的树高值之间差异不显著，即研究所得的地位指数曲线能够较好地反映横断山栎类天然林树高的生长趋势。

4.2 落点检验

统计71 块样地的年龄—树高值，将每块样地的树高值绘制在地位指数曲线簇上（见图3），从其落点分布可以看出，有69 块样地的树高值落在地位指数曲线簇范围内，落点精度高达97.18%。表明本研究所编制的地位指数表能够准确反映出横断山栎类天然林的立地质量。

图3 地位指数曲线簇与落点检验Fig.3 Site index cluster and placement test

5 地位指数表的应用

已知样地林分年龄和优势木平均高，查所编制的横断山东部栎类天然林地位指数表，可得到样地所对应的地位指数值。将该方法应用到未参与建模的71 块样地中，可得到每块样地的地位指数值（见表7）。

表7 各检验样地的地位指数值Tab.7 Site index value of each sample plot

（续表 7）

71 块样地的地位指数值在8～18 m 范围内均有分布，其中8～14 指数级分布的样地数较多，分别为13、19、13、16 块，占样地总数的85.92%，18 指数级的样地仅有1 块，占样地总数的1.41%。所有样地的地位指数平均值为11.77 m，与基准年龄时导向曲线上的树高值接近，再次表明本研究所编横断山东部栎类天然林地位指数表在研究区内适用性强。

6 结论与讨论

本研究选用横断山东部石棉县的211 块栎类天然林样地调查数据，以其中140 块样地数据为建模数据，剩余71 块样地数据为检验样本。选择8 个常用的树高生长方程作为拟合年龄与树高的数学模型，得出树高生长曲线，最终展开得到横断山栎类天然林地位指数表。拟合结果显示，逻辑斯蒂方程的确定系数（R2=0.992）、拟合精度（P=98.23%）最高，残差平方和最小（SSE=0.848），最终选定其作为导向曲线：H=13.894 17/（1+4.065 53e-0.07257A），以基准年龄40 年，指数级距2 m，运用标准差调整法展开形成了横断山栎类天然林地位指数表。χ2检验和落点检验结果显示，所编的地位指数表精度高、适用性强，且能准确反映研究区栎类天然林的地位质量。

表8 各检验样地地位指数分布情况Tab.8 Distribution of site index of each sample plot

本次参与编表的样地在5 个龄组中均有分布，其中以幼、中龄林为主，说明编制的地位指数表能够客观反映研究地不同生长阶段栎类天然林的地位质量，并可广泛运用于该地区栎类天然林的地位质量评价。

研究区栎类天然林的平均地位指数在12 m 左右，树高生长相对缓慢，这是由于该地区的森林遭受到较为严重的人为干扰，许多林分为多代萌生的次生林，其生长过程与原始实生林存在较大的差异。但人为干扰为可控因素，且天保工程实施以来，对栎类等天然林的保护起到了重要作用，该地区的成、过熟林数量将逐渐增加，后期将增加成、过熟林样地数据，进一步提高地位指数表的适用范围，更好的运用于生产实践。