楊錦昌，尹光天，馮昌林，李榮生，鄒文濤

（1. 中國林科院熱帶林業(yè)研究所，廣東 廣州 510520；2. 中國林科院熱帶林業(yè)實驗中心，廣西 憑祥 532600）

林分生長和收獲模型是預(yù)測林分的生長、產(chǎn)量及生長潛力的科學依據(jù)，也是有效制訂管理措施、正確評價經(jīng)營成效和合理利用森林資源的重要工具[1]。近年來，隨著科技的發(fā)展和研究的不斷深入，林分生長和收獲模型在林業(yè)生產(chǎn)上的應(yīng)用空前廣泛，為實現(xiàn)林業(yè)管理的科學性、定量性、準確性和持續(xù)性提供了重要的參考[2～9]。相比之下，因棕櫚藤具刺、攀緣等生物學特性的限制，國內(nèi)外關(guān)于棕櫚藤生長收獲模型的研究偏少。曾炳山等選擇冪函數(shù)描述了黃藤生長規(guī)律并模擬了莖長分配結(jié)構(gòu)[10]，Lee利用Richards和Logistic等幾種方程研究了西加省藤和疏刺省藤的生長和收獲量[11～12]；近年來，楊錦昌等通過比較 5種生長方程建立了單葉省藤（Calamus simplicifolius）和黃藤（Daemonorops margaritae ）人工林生長模型[13]?？傮w上看，棕櫚藤生長收獲模型的研究基礎(chǔ)比較薄弱，尤其莖長分布模型的研究亟待加強。本文以我國華南地區(qū)重要商品藤種單葉省藤為對象，通過理論生長方程的篩選和采用兩次曲線擬合法來研制單葉省藤莖長分布模型，為準確預(yù)測藤林生長收獲和科學經(jīng)營藤林提供參考。

1 材料與方法

1.1 試驗地概況

試驗地位于廣西憑祥市熱帶林業(yè)實驗中心英陽試驗站，22° 07′ N，106° 44′ E，屬南亞熱帶季風氣候，干濕季節(jié)交替明顯。年平均氣溫21.6℃，最冷月1月平均氣溫13.3 ℃，極端最低氣溫0.1℃，最熱月7月平均氣溫27.7℃，極端最高氣溫37.9℃，≥10℃年積溫7 596℃；年平均降水量1 388 mm，蒸發(fā)量為1 275 mm，相對濕度81.75%；年日照1 512 h，風速0.84 m/s。土壤為花崗巖發(fā)育而成的磚紅壤性紅壤，土層厚度一般在2 m以上，土壤呈強酸性，pH值4.4 ～ 4.9，表土層腐殖質(zhì)含量26.1 ～ 40.2 g/kg，全N、全P和全K含量分別為1.024 ～1.374、1.032 ～ 1.30和 1.68 ～ 2.23 g/kg。試驗地海拔310 ～ 370 m，上層樹種為石梓（Gmelina arborea），1982年造林，密度為1 650株/hm2。

1.2 試驗設(shè)計與觀測

試驗設(shè)計分為2個部分，即生長監(jiān)測和采收設(shè)計。通過典型選取，在石梓林分下設(shè)置2塊15 m×12 m的單葉省藤固定標準地以監(jiān)測藤林生長動態(tài)，株行距為2 m×3 m，藤叢數(shù)為30；采收設(shè)計根據(jù)隨機區(qū)組設(shè)計安排了4個區(qū)組，每區(qū)組5個小區(qū)，小區(qū)水平排列，大小為20 m×3 m，其形狀為長帶形，小區(qū)內(nèi)沿坡度按2 m株距單行種植10株，并加以標記和編號。

固定標準地的調(diào)查時間為第3、第4、第5、第6、10第和第13年，調(diào)查時均實測標準地內(nèi)所有藤叢的植株數(shù)和莖長。采收試驗設(shè)計中，在第5至第9年每年對單葉省藤試驗地各區(qū)組的藤叢進行全面調(diào)查，調(diào)查因子包括母莖長、萌莖長及萌蘗株數(shù)；同時分別于第5、第6、第7、第8、第9年對各區(qū)組中的第1、第2、第3、第4和第5小區(qū)依次采收，采收時選擇株高5 m以上的植株，然后除去葉片和葉鞘，最后再實測長度。植株的長度在未采收時用測高桿測量，而采收后的植株則用皮尺測量；根據(jù)植株長度采前測量值與采后實測值之間的相關(guān)關(guān)系，對測高桿測量6 m以上植株長度的測量值進行校正。所有數(shù)據(jù)在計算機上利用SPSS12.0、Excel2003軟件進行相關(guān)的處理和分析。

2 結(jié)果與分析

藤叢由年齡各異、莖長參差不齊的植株組成，這種特殊性使得藤林既區(qū)別于同齡純林，又與異齡林不盡相同。藤林由藤叢組成，而藤叢又由母莖和若干萌莖構(gòu)成；若將母莖從各藤叢中提取出來則可看成同齡純林，而萌莖則構(gòu)成異齡林[14]；因此，單葉省藤人工林可看作由同齡純林和異齡林組成的復(fù)合體。基于上述分析，將母莖和萌莖區(qū)分開來，選擇在林分直徑結(jié)構(gòu)研究中應(yīng)用最多的3種生長方程[8,14]，采用兩次曲線擬合法建立單葉省藤人工林莖長結(jié)構(gòu)模型。

表1 單葉省藤莖長劃分標準Table 1 Division of stem length of C. simplicifolius

2.1 長度級的劃分及生長方程的表達

為研究莖長分布規(guī)律，根據(jù)單葉省藤生物學特性以及藤條采收（5 m以上）和收購要求，主要按2 m為莖長間距劃分長度級；同時，考慮到1.5 m以下尤其是不足0.5 m的莖長特別多，特將1.5 m以下的莖長按0.5 m和1.0 m的莖長間距進行劃分，最后形成表1的莖長劃分標準。

鑒于Logistic、Weibull和Richards方程已廣泛應(yīng)用于描述人工林直徑分布規(guī)律，且擬合精度高[14～15]；為此，選擇這三個方程對單葉省藤人工林的莖長分布進行擬合，其數(shù)學表達式及解析性見文獻[8]。利用累加生成、標準化數(shù)據(jù)處理方法將林分莖長分布初始數(shù)據(jù)轉(zhuǎn)化成（0, 1] 區(qū)間數(shù)列，這樣Logistic、Weibull和Richards分布的累積分布函數(shù)則分別簡化為：

式中，F(xiàn)(x)為相對累計頻率，x為莖長，b和c為待求解參數(shù)。

2.2 生長方程的擬合效果

基于觀測數(shù)據(jù)，利用SPSS軟件之非線性回歸法求解方程中的參數(shù)，結(jié)果如表2和表3所示?？傮w上，3種方程對單葉省藤母莖長和萌莖長分布的擬合效果比較理想，相關(guān)指數(shù)在0.9以上，但3種方程的擬合精度因莖長類型而有所差異。擬合母莖長分布時，Logistic和Weibull方程均優(yōu)于Richards方程；前二者擬合相關(guān)指數(shù)均在0.99以上，高者達0.999 9，而后者擬合母莖長分布的相關(guān)指數(shù)變化范圍為0.975 0 ～ 0.999 8；擬合萌莖長分布時，Logistic和Richards方程的相關(guān)指數(shù)均在0.98以上，其擬合效果優(yōu)于Weibull方程，特別是隨著莖長和株數(shù)的增加，Weibull方程模擬萌莖長分布的相關(guān)指數(shù)有下降趨勢，這種優(yōu)越性更加突出。從擬合精度來看，Logistic方程擬合莖長分布規(guī)律的效果最好，而Weibull和Richards方程則分別適合于描述母莖長和萌莖長的分布特征。

表2 理論生長方程對母莖長分布的擬合效果Table 2 Fitting effect of growth equations with length distribution of mother stem

2.3 參數(shù)預(yù)測模型

要建立良好的參數(shù)預(yù)測模型，從而客觀預(yù)測莖長分布動態(tài)，除考慮生長方程擬合莖長分布的精度之外，也要特別重視方程中參數(shù)的變化規(guī)律。為此，分別選擇 Logistic、Weibull方程中的母莖長分布參數(shù)及 Logistic、Richards方程中的萌莖長分布參數(shù)與藤林因子建立相關(guān)關(guān)系，以相關(guān)指數(shù)為入選標準，最后分別建立 Weibull方程中的b、c參數(shù)與母莖長相關(guān)關(guān)系的參數(shù)預(yù)測模型以及Richards方程中的b、c參數(shù)與萌莖長和萌蘗數(shù)相關(guān)關(guān)系的參數(shù)預(yù)測模型。

Weibull方程中參數(shù)預(yù)測模型：

式（4）至式（7）中：Wb、Wc、Rb、Rc分別代表Weibull和Richards方程中的參數(shù)b和c；Nss為藤叢平均萌蘗數(shù)；Lms、Lss分別為母莖均長和萌莖均長。

為確定林分莖長的上下限，建立了下列莖長極值模型。

式（8）至式（10）中：Lmin、Lmax、maxL′ 分別表示母莖長下限、母莖長上限和萌莖長上限。

表3 理論生長方程對萌莖長分布的擬合效果Table 3 Fitting effect of growth equations with length distribution of sucker stem

2.4 林分生長模型

建立參數(shù)預(yù)測模型后，為預(yù)測不同年齡莖長分布動態(tài)，引用了已建的林分生長模型[14]。

式（11）至式（13）中：Lms為母莖均長，TLs為叢萌莖總長，Nss為叢萌蘗數(shù)，t為年齡，Lss= TLs/Nss。

建立了莖長分布模型后，只要母莖均長、萌莖均長和萌蘗數(shù)已知，相應(yīng)方程中的參數(shù)b和c便可求出；然后將參數(shù)b和c代入式對應(yīng)的累積分布函數(shù)，再通過確定莖長的上下限，求出母莖長和萌莖長中某一長度級的相對累計頻率和該長度級的株數(shù)，最后將母莖長中各長度級的株數(shù)和萌莖長中各長度級的株數(shù)匯總，就可實現(xiàn)藤林莖長分布的模擬和預(yù)測。各長度級的株數(shù)可通過下式計算：

2.5 莖長分布模型的檢驗

表4 單葉省藤人工林莖長分布實際值與預(yù)測值的比較及檢驗Table 4 Comparison between practical and predicted values of stem length distribution and its test

為了檢驗參數(shù)預(yù)測模型的適用性，選擇了未參與建模的5組單葉省藤標準地數(shù)據(jù)進行驗證。同時，為了說明采用兩次曲線擬合法建立莖長分布模型（方法 2）是否提高不同長度級株數(shù)的預(yù)測精度，同未區(qū)分母莖和萌莖而通過經(jīng)三個理論方程篩選后直接擬合而建立的莖長分布模型（方法1）進行比較，其結(jié)果如表4所示。由表4可知，不論是從預(yù)測株數(shù)與實際株數(shù)的接近程度來考慮，還是選擇株數(shù)剩余標準差或卡方統(tǒng)計量來衡量，方法2的預(yù)測精度明顯高于方法1。另一方面，通過參數(shù)預(yù)測模型所預(yù)測的各長度級株數(shù)與實際值比較吻合，每組數(shù)據(jù)均通過卡方檢驗，說明莖長分布模型預(yù)測精度高，可在實際中應(yīng)用。

3 結(jié)論與討論

三種生長方程應(yīng)用于描述單葉省藤人工林莖長分布規(guī)律時總體上表現(xiàn)出了良好的模擬性能，但擬合精度因不同生長方程而異；Logistic方程適于描述母莖長和萌莖長結(jié)構(gòu)規(guī)律，而Weibull和Richards方程則分別適合于擬合母莖長和萌莖長的分布特征。多個生長方程的比較和分析克服了以往選擇單一方程建立林分模型的缺陷，為準確描述藤林生長規(guī)律提供了理論依據(jù)。

Logistic方程在描述莖長分布結(jié)構(gòu)時優(yōu)于Weibull和Richards方程，但在建立參數(shù)預(yù)測模型時又遜于后兩個方程，這可能與方程的解性特性有關(guān)[8]；建立莖長分布模型，不僅要考慮生長方程擬合莖長分布的精度，也要特別重視方程中參數(shù)的變化規(guī)律。因而，基于多個方程的篩選和比較，確定最適宜建立林分莖長分布模型且性能較好的生長方程是可能的，也是非常必要的。

采用兩次曲線擬合法建立的莖長分布模型的預(yù)測精度明顯優(yōu)于直接擬合法，從而克服了直接擬合法在描述莖長分布結(jié)構(gòu)時按預(yù)測精度偏低的問題；將母莖長和萌莖長相區(qū)分，采用參數(shù)預(yù)測法建立莖長分布模型為動態(tài)模擬單葉省藤林分莖長分布規(guī)律和預(yù)測藤林產(chǎn)量提供了良好的途徑。

方程擬合精度的高低取決于方程表達形式和林分數(shù)據(jù)的差異；對于相同的方程而言，不同的林分數(shù)據(jù)擬合所得到的精度不盡相同；反之亦然[8]。本文僅以特定條件下單葉省藤人工林為研究對象建立了莖長分布模型，由于藤林莖長分布結(jié)構(gòu)與立地條件、種植密度和上層樹種有關(guān)，因而模型的應(yīng)用具有一定的局限性，適用于不同條件下的莖長分布模型有待于深入研究。另一方面，本文所建立的莖長分布模型在統(tǒng)計學上獲得滿意的預(yù)測結(jié)果，但忽略了一些大的長度級株數(shù)的表達，這成為今后研制莖長分布模型時必須充分考慮的重要因素。

[1] 馬豐豐，賈黎明. 林分生長和收獲模型研究進展[J]. 世界林業(yè)研究，2008，21（3）：21－27.

[2] 江傳陽. 天然闊葉林生長收獲動態(tài)預(yù)估模型的研究[J]. 林業(yè)勘查設(shè)計，2007（2）：9－12.

[3] 江希鈿，莊晨輝，陳信旺，等. 免疫進化算法在建立地位指數(shù)曲線模型中的應(yīng)用[J]. 生物數(shù)學學報，2007，22（3）：515－519.

[4] Kokkila T，M?kel? A，F(xiàn)ranc A，et al. Comparison of distance-dependent and distance-independent stand growth models-Is perfect aggregation possible [J]. For Sci，2006，52（6）：623－635.

[5] Bermejo I，Ca?ellas I，Miguel S A. Growth and yield models for teak plantations in Costa Rica[J]. For Ecol Manag，2004（189）：97－110.

[6] Valkonena S，Valsta L. Productivity and economics of mixed two－storied spruce and birch stands in Southern Finland simulated with empirical models[J]. For Ecol Manag，2001（140）：133－149.

[7] 胡曉龍. 長白落葉松林分斷面積生長模型的研究[J]. 林業(yè)科學研究，2003，16（4）：449－452.

[8] 段愛國，張建國，童書振. 6種生長方程在杉木人工林林分直徑結(jié)構(gòu)上的應(yīng)用[J]. 林業(yè)科學研究，2003，16（4）：423－429.

[9] 胡喜生，周新年，蘭樟仁，等. 人工林桉樹胸徑分布模型的研究[J]. 福建林學院學報，2008，28（4）：314－318.

[10] 曾炳山，許煌燦，尹光天，等. 黃藤藤叢結(jié)構(gòu)和生長的研究[J]. 林業(yè)科學研究，1993，6（4）：414－422.

[11] Lee Y F. Some models for estimating rattan growth and yield[J]. J Trop For Sci，1994（6）：346－355.

[12] Lee Y F. Genetic and ecological studies relevant to the conservation and management of some Bornean Calamus species[D]. Aberdeen：University of Aberdeen，1995.

[13] 楊錦昌，尹光天，李榮生，等. 5種生長方程在建立2種藤林生長模型中的應(yīng)用[J]. 福建林學院學報，2007，27（3）：217－221.

[14] 張建國，段愛國. 理論生長方程對杉木人工林林分直徑結(jié)構(gòu)的模擬研究[J]. 林業(yè)科學，2003，39（6）：55 61.

[15] 陳永富. 短周期桉樹人工林直徑分布模型研究[J]. 林業(yè)科學研究，林業(yè)科學，2008，21（增刊）：50－54.