華偉平 ，丘 甜，葉宏萌，武佳慧，池上評，張傳海，蓋新敏

（1.福建省生態(tài)產(chǎn)業(yè)綠色技術(shù)重點(diǎn)實(shí)驗(yàn)室，福建 武夷山 354300；2.武夷學(xué)院 生態(tài)與資源工程學(xué)院，福建 武夷山 354300；3.武夷學(xué)院 商學(xué)院，福建 武夷山 354300；4.福建省林業(yè)勘察設(shè)計(jì)院，福州 350001；5.寧德市林業(yè)局，福建 寧德 352100）

密控圖（SDMD），全稱林分密度控制圖，是以密度效應(yīng)模型為基礎(chǔ)，通過建立等直徑線、等樹高線、最大密度線、自然稀疏線等各種模型，確定林分密度與產(chǎn)量間的經(jīng)營圖表[1,2]。密控圖對掌握林分密度、預(yù)測林分生長、造林規(guī)劃設(shè)計(jì)、資源清查等方面提供重要的經(jīng)營管理手段[1-3]。因此，各國對重要的用材樹種均有研制相應(yīng)的密控圖，如Penner等研制了香脂冷杉（Abies balsamea）、紅云杉（Picea rubens Sarg）、黑云杉（P.mariana）和白云杉（P.glauca Voss）的密控圖[4]；Long等對新不倫瑞克省香脂冷杉的密控圖進(jìn)行了研究[5]。國內(nèi)起源于19世紀(jì)80年代末期，之后出現(xiàn)多樹種、多類型的密控圖，如杉木(Cunninghamia lanceolata)、馬尾松 (Pinus massoniana)、落葉松(Larix gmelinii)、福建柏(Fokienia hodginsii)等[1-3,6-10]。

可變密度收獲表是用來描述林分各調(diào)查（或測樹因子）的生長規(guī)律的林業(yè)數(shù)表，該表是在給定林分密度的情況下，利用各測樹因子的生長規(guī)律模擬時(shí)間上的變化情況，對指導(dǎo)營林生產(chǎn)具有重要的意義[11]。我國在80年代中期，張少昂學(xué)者利用Von Bertalanffy方程研制了興安落葉松天然林可變密度收獲表[12]，李希菲等利用Richards方程擬合地位級方程并編制了一個(gè)完整的可變密度收獲表[13]。這些研究對我國的可變密度收獲表研制打下重要的理論基礎(chǔ)。

過去的編表往往是可變密度收獲表與林分密度控制圖自成體系，從而使結(jié)果無法兼容，造成輸出的林分測樹因子不一致的問題。基于此，本次研究以黃山松(Pinus taiwanensis)為研究對象、以密度效應(yīng)模型為基礎(chǔ)模型，通過建立等樹高線、等直徑線、最大密度線、自然稀疏線等系列模型繪制黃山松密控圖，并結(jié)合地位指數(shù)模型導(dǎo)出黃山松可變密度收獲表。

黃山松為福建山區(qū)高海拔樹種，干形通直，不僅具有較高的材用價(jià)值，而且具有獨(dú)特的“黃山松文化”備受人們喜愛，目前有學(xué)者對黃山松的生態(tài)群落、生物化學(xué)、林業(yè)相關(guān)數(shù)表等方面做過研究[14-17]，但兼容性黃山松密控圖與可變密度收獲表方面的研究未見報(bào)道。通過本次研究，為制定合理的黃山松經(jīng)營密度提供依據(jù)。

1 研究材料

在福建省南平武夷山、泉州戴云山、龍巖蓮臺山、南平茫蕩山、泉州九仙山等黃山松適宜分布區(qū)，歷經(jīng)10年的黃山松樣地調(diào)查，獲取了肥沃地類、較肥沃地類、中等肥沃地類、瘠薄地類的253塊樣地?cái)?shù)據(jù)(每個(gè)面積0.067 hm2)，包括每木胸徑、樹高、空間位置、林下植被。根據(jù)樣地調(diào)查結(jié)果，伐倒1～3株樣木，利用圍徑尺測定伐倒木胸徑，皮尺測定樹干長度；按1 m為1個(gè)區(qū)分段截取原木，分別測定各區(qū)分段帶皮、去皮直徑，以及梢頭底直徑和長度，用于建立黃山松二元材積方程，估算林分蓄積量，樣地主要林分因子見表1。

表1 樣地主要林分因子Tab.1 Main stand factors of sample plots

2 研制思路

2.1 林分密度效應(yīng)模型

林分密度效應(yīng)模型是用于制作林分密控圖的林分模型。理論基礎(chǔ)由吉良龍夫提出，用于描述個(gè)體與群體密度的關(guān)系式，見式(1)。

式中：y與植物群落有關(guān)的林木個(gè)體產(chǎn)量；A、B為模型參數(shù)；N為林木株數(shù)；公式含義下同。

之后，許多學(xué)者將式(1)用于描述多年生的森林植被。但安滕·貴等提出樹干材積不僅受年齡影響，也受立地條件影響。因此，對吉良龍夫的密度效應(yīng)模型進(jìn)行了改進(jìn)，改進(jìn)后的林分密度效應(yīng)模型為：

式中：Hu為林分優(yōu)勢高；a1、a2、b1、b2為模型參數(shù)。

江希鈿等[18]研究發(fā)現(xiàn)，吉良龍夫提出的林分密度效應(yīng)模型存在林木生長最大速度(vmax)在L/2之處到達(dá)，其中L為林木個(gè)體因子產(chǎn)量y的上限值。但林木個(gè)體或者林分的vmax可能在L/2之前或之后，導(dǎo)致密控圖出現(xiàn)精度較低的現(xiàn)象。因此，本次采用下式作為黃山松密度效應(yīng)機(jī)制模型。

式中：c1、c2為模型參數(shù)。當(dāng)式(3)中的時(shí)，即為式(2)的密度效應(yīng)模型；當(dāng)不等于1時(shí)，克服了式(2)存在的問題。因此，本次將式(3)作為黃山松密度效應(yīng)模型一般式。當(dāng)y表示林分平均胸徑(D)時(shí)，則林分平均胸徑模型式為：

式中：M為林分蓄積量，m3。式(4)、式(5)不僅是研制黃山松林分密度控制圖的等樹高線和等直徑線，也是編制可變密度收獲表的重要模型。本次采用改進(jìn)單純形法[19,20]求解林分密度效應(yīng)模型參數(shù)。

2.2 最大密度法則

當(dāng)現(xiàn)實(shí)林分密度很大時(shí)，林木間由于競爭強(qiáng)烈而出現(xiàn)自然稀疏現(xiàn)象，密度效應(yīng)模型將不適用。這種情況可用最大密度法則（又稱3/2法則）進(jìn)行描述，方程如下：

式中，α、β為待求參數(shù)。

前人研究表明β約為3/2，所以最大密度法則式(6)又稱3/2法則(或二分之三法則)。根據(jù)林分密度效應(yīng)模型和最大密度法則式(6)即可研制密控圖。

2.3 林分密度指數(shù)

為了確定林分不同胸徑的株數(shù)，需要建立林分密度指數(shù)模型。林分密度指數(shù)（SDI）是將任意林分單位面積上的株數(shù)通過林分平均胸徑換算成標(biāo)準(zhǔn)平均胸徑時(shí)所具有的株數(shù)，是賴內(nèi)克在1933年提出，認(rèn)為單位面積株數(shù)與林分平均胸徑之間存在冪函數(shù)關(guān)系，用于評定林分密度的一種指標(biāo)；丹尼爾（1979）年在賴內(nèi)克的基礎(chǔ)上進(jìn)一步認(rèn)為，該指標(biāo)受林分年齡和立地質(zhì)量的影響很小。該指標(biāo)在評定林分競爭情況方面簡單直觀，并易于測算。本研究選用該指標(biāo)作為編制黃山松可變密度收獲表的林分密度指標(biāo)，其林分最大密度方程見下式。

式中：N是每公頃最大株數(shù)；d1、d2為待求參數(shù)，其中d2大于為林分最大株數(shù)對應(yīng)的平均胸徑。通過林分密度指數(shù)的定義，可得到林分密度指數(shù)公式為：

利用樣地實(shí)測數(shù)據(jù)估算，得到d2為1.248 8，所以黃山松林分密度方程為：

3 研制結(jié)果

3.1 林分密度控制圖的研制

3.1.1 等直徑線與等樹高線的建立

在林分密度控制圖研究中，采用本次密度效應(yīng)倒數(shù)模型一般式，并應(yīng)用改進(jìn)單純形法對191塊樣地求解模型參數(shù)，得到等直徑線模型：

由林分密度方程可得林分單位面積株數(shù)與林分密度指數(shù)（SDI）的關(guān)系式為：

將上式代入到等直徑線模型中，得到林分平均胸徑與林分密度指數(shù)SDI關(guān)系式為：

由式(12)可知，在給定林分密度的情況下，并利用地位指數(shù)確定林分不同年齡的優(yōu)勢高，通過迭代即可得到林分平均胸徑D。

選用191塊樣地實(shí)測數(shù)據(jù)作為建模樣本，經(jīng)改進(jìn)單純形法估算，建立的林分蓄積量模型為：

由式(13)可知，確定了林分不同年齡的優(yōu)勢高和株數(shù)即可得到林分蓄積量。分別給定優(yōu)勢高，將不同的株數(shù)密度代入式(13)，求出相應(yīng)的蓄積量，并在雙對數(shù)坐標(biāo)上點(diǎn)出，即可連成等樹高線。

3.1.2 最大密度線的建立

最大密度線是林木生長過程中，M達(dá)到最大、對應(yīng)的N最多時(shí)的曲線。根據(jù)最大密度法則，采用式(6)作為林分最大密度線模型。而確定最大密度線模型參數(shù)α、β值的常用方法為提升法，但該方法理論依據(jù)不夠充分，受樣地?cái)?shù)據(jù)和主觀因素影響。因此，本次張惠光利用具有理論依據(jù)技術(shù)確定福建柏最大密度線的方法[9]確定黃山松的最大密度線模型參數(shù)。根據(jù)樣地的最大蓄積量和最大株數(shù)，采用回歸分析技術(shù)建立最大密度線模型為：

3.1.3 自然稀疏線

林分生長過程中，由于空間資源限制，林木之間出現(xiàn)競爭而導(dǎo)致立木株數(shù)減少的現(xiàn)象，即自然稀疏過程，本次用下式描述黃山松的自然稀疏規(guī)律。

式中：N0為單位面積林分初植密度；N為任意階段林分密度。

將最大密度線模型式(14)中已知的參數(shù)值α、β代入式(15)，給定不同單位面積林分初植密度N0和任意階段林分密度N，得到對應(yīng)的林分蓄積量M，即單位面積林分任意初植密度的自然稀疏線。

3.1.4 等疏密度線

編制黃山松密控圖過程中，還需要繪制等疏密度線(又稱等比產(chǎn)量線)。繪制該線的方法為：根據(jù)等樹高線，按1.0、0.9、0.8、…、0.1的最大蓄積量計(jì)算下降量，并在密控圖中，將相同比值點(diǎn)連接起來，得到與最大密度線基本平行的等疏密度線。

3.1.5 精度檢驗(yàn)

密控圖包括等樹高線、等直徑線、最大密度線、自然稀疏線、等疏密度線，其中等樹高線、等直徑線決定密控圖的精度。因此，將未參加建模的62塊數(shù)據(jù)對等樹高線和等直徑線進(jìn)行相對誤差、精度、F檢驗(yàn)。結(jié)果顯示，等樹高線相對誤差5.23%，精度為95.37%，F(xiàn)=1.32＜F0.05(2,60)=3.07。等直徑線相對誤差5.92%，精度為94.24%，F(xiàn)=1.87＜F0.05(2,60)=3.07。說明黃山松林分密度控制圖的理論值和實(shí)際值差異不顯著，實(shí)際使用誤差小、精度高，可用于編制可變密度收獲表。

3.2 地位指數(shù)模型的建立

通過林分密度控制圖導(dǎo)出可變密度收獲表或者通過建立方程編制可變密度收獲表均需要建立地位指數(shù)。本次采用地位指數(shù)模型來描述優(yōu)勢高生長情況，結(jié)合不同立地質(zhì)量等級的黃山松樣地?cái)?shù)據(jù)，本次構(gòu)建如下公式作為黃山松地位指數(shù)模型。

式中：to為林分基準(zhǔn)年齡；t為林分年齡。

其中：

式中：a1、a2、a3、a4、b1、b2、b3、b4、c1、c2、c3、c4為待求參數(shù)；k1、k2、k3、k4為4種立地質(zhì)量等級，取值為0或1，即k1取1時(shí)，k2、k3、k4均為0，以此類推。

為提高擬合精度，采用混合蛙跳算法[21]來估計(jì)模型參數(shù)，結(jié)合樣地實(shí)測數(shù)據(jù)，求得各方程的參數(shù)見表2，相關(guān)系數(shù)為0.988。

表2 地位指數(shù)導(dǎo)向曲線擬合結(jié)果Tab.2 Fitting results of site index guidance curve

3.3 可變密度收獲表

在給定林分密度指數(shù)的情況下，根據(jù)上述建立的地位指數(shù)模型、林分密度指數(shù)、林分平均胸徑模型、林分蓄積量模型，即可編制可變密度收獲表。副林木資料記錄較少，對副林木將不做進(jìn)一步探討。具體步驟如下：

（1）給定立地質(zhì)量等級和地位指數(shù)，通過地位指數(shù)模型得到不同年齡的優(yōu)勢高；

（2）給定林分密度指數(shù)，利用林分平均胸徑模型通過迭代法得到不同林分平均高的平均胸徑；

（3）通過步驟（2）得到不同年齡下的平均胸徑后，利用林分密度指數(shù)得到不同年齡的林分單位面積株數(shù)；

（5）最后林分蓄積量模型確定不同年齡的林分蓄積量。

現(xiàn)給定立地質(zhì)量等級為肥沃，地位指數(shù)為16 m，密度指數(shù)為1 400株/hm2，根據(jù)上述步驟即可編制黃山松可變密度收獲表，具體編制結(jié)果見下表。

表3 黃山松可變密度收獲表Tab.3 Variable density harvest of Pinus taiwanensis

4 結(jié)論

為了合理利用林地生產(chǎn)力，經(jīng)營好黃山松，在前人研究的基礎(chǔ)上，根據(jù)植物種群產(chǎn)量增長模式和最終產(chǎn)量恒定理論，本次采用的密度效應(yīng)機(jī)制模型克服了安滕·貴等提出改進(jìn)密度效應(yīng)法模型存在的問題。利用樣地?cái)?shù)據(jù)，建立了等樹高線、等直徑線、最大密度線模型，編制了精度較高的黃山松密控圖，從而為掌握黃山松經(jīng)營密度、定量間伐、林分生長預(yù)測、資源清查、規(guī)劃造林設(shè)計(jì)等森林經(jīng)營管理提供科學(xué)依據(jù)。

以本研究建立的黃山松林分效應(yīng)模型為基礎(chǔ)，通過建立以立地質(zhì)量等級為啞變量的地位指數(shù)模型得到不同年齡的優(yōu)勢高，并在給定立地質(zhì)量等級為肥沃、地位指數(shù)為16 m、密度指數(shù)為1400株/hm2，根據(jù)編表步驟，得到黃山松可變密度收獲表，保證了林分密控圖與可變密度收獲表一致性，解決了不兼容的問題。

收獲表屬于林業(yè)數(shù)表中的經(jīng)營類數(shù)表，是全林分收獲模型的重要內(nèi)容之一，包括按“標(biāo)準(zhǔn)林分”編制標(biāo)準(zhǔn)收獲表、“平均密度”編制的現(xiàn)實(shí)收獲表及可變密度收獲表等3類，而可變密度收獲表與其他2類林業(yè)數(shù)表相比，既能反應(yīng)不同林分密度下的蓄積量，又不受標(biāo)準(zhǔn)收獲表的“標(biāo)準(zhǔn)林分”限制。因此，通過研制黃山松林分密控圖導(dǎo)出可變密度收獲表，不僅能將黃山松林分控制在合理的經(jīng)營密度范圍內(nèi)，同時(shí)，對制定合理的黃山松經(jīng)營措施具有重要意義。