天然次生林直徑分布的研究

天然次生林直徑分布的研究

王硯峰，王鶴智

(國家林業局調查規劃設計院，北京 100714)

摘要該研究應用東北林區調查的119天然林標準地數據，通過繪制各林分直徑分布直方圖，分析林分的直徑分布規律，運用SB分布、指數分布、正態分布、Weibull分布和β分布函數對其進行擬合，建立直徑回歸模型，并采用矩解法求解參數，更精確地比較各種分布函數對實際數據的擬合效果。研究結果表明：林木株數隨著直徑的增大呈減少的趨勢；指數分布和Weibull分布的擬合效果較好，但指數分布優于Weibull分布。觀察指數分布函數的各參數的變化規律，發現參數a基本不變，而參數b在變化。繪制參數b與林分調查因子間的散點圖，發現指數分布中的參數b與林分平均直徑有顯著的相關性。通過相關性分析，選擇函數進行擬合，結果表明：冪函數的擬合效果最好，相關系數最高。通過對參數b的再參數化建立模型，最終得到最優的天然次生林林直徑分布模型。

關鍵詞天然次生林；直徑分布；林分調查因子；林分直徑

中圖分類號S757

作者簡介王硯峰(1978-)，男，山東榮口人，工程師，碩士，從事森林經理研究。

收稿日期2015-06-05

Study on DiameterDistribution of Natural Secondary Forest

WANG Yan-feng, WANG He-zhi (State Forestry Bureaue Survey Scheme Designing Institute, Beijing 100714)

AbstractUse SB distribution function, exponential distribution function, normal distribution function, Weibull distribution function and beta distribution function to fit the transformation regularities, set up parameter regression model. In order to compare various distribution function of the actual data fitting effect more precisely, the solution method using torque. Research results show: the rule of diameter distribution is that with the increasing of the diameter, the number of trees is reduced; the fitting effects of exponential distribution and Weibull distribution are good, but the exponential distribution is better.Observing each parameter variation of exponential distribution function, it can be seen that parameter a is constant while parameter b is changing. By making chat between parameter b and stand description factor, it shows that parameter b and average diameter has significant correlation. Through the correlation analysis, select the appropriate functions to fit them. The results show that the power function fits best, the correlation coefficient is the highest. Combining the relationship between parameter b and average diameter, obtain the general diameter distribution function of natural secondary distribution function.

Key wordsNatural secondary forest; Diameter distribution; Stand description factor; Stand diameter

對于林木的直徑來說，其大小與林分的蓄積量存在著非常密切的關系，總體上來說對于林分結構的合理確定是充分發揮其各種效益的基礎和條件[1-5]。因此，需要通過外業調查等方式獲得林分結構的詳細信息，充分掌握林木生長和退化的整個過程，從而對更加充分地了解林分競爭起到非常好的作用，同時也有利于充分發揮林分結構的各種特性和作用，為林分經營奠定基礎[6-8]。

在通常條件下，需要對森林的直徑結構進行合理的確定，在實踐的基礎上，需要合理地確定經營措施、制定主伐年齡[9]。林木直徑代表了森林對于經營措施的綜合反映，是一個綜合體。因此，需要不斷地通過外業調查來確定經營指標，從而為林分的正常生長提供好的條件和基礎。當然，在目前條件下，不同的國家對于林分直徑的確定方法也存在非常大的差異性，在調查方法上很難做到統一[10-13]。N·Nishimura[14]在早期的研究中，對日本溫帶地區闊葉樹種的直徑分布進行了研究，采用的方法是直方圖，結果顯示，日本溫帶地區的常綠闊葉林的直徑總體上呈現雙峰分布。但綜合國內外的研究現狀可以知道，目前對直徑分布應用最多的還是Weibull分布，這是由于Weibull分布所有的參數具有非常大的靈活性，具有明顯的生物學意義，能夠對樹木的生長等進行合理的描述和解釋，能夠進行科學的預估和評價。因此，這些年來，該分布得到了很多學者的廣泛應用，也取得了非常大的成果。該研究應用在東北林區調查的119天然林標準地數據，通過繪制各林分直徑分布直方圖，分析林分的直徑分布規律。

1 研究地概況及方法

1.1 研究地區概況及數據收集帽兒山實驗林場位于我國黑龍江省哈爾濱市，距離哈爾濱市區僅有100 km左右。該地區的氣候屬于溫帶季風氣候，非常干燥，冬季寒冷、時間長，夏季炎熱、時間短，季節之間差別非常明顯。帽兒山林場年均氣溫為2.8 ℃，年降水量為723.8 mm，主要集中在6～8月份。

穆棱林業局位于黑龍江省東南部地區，屬于長白山地區的一部分，與老爺嶺相鄰，位于我國黑龍江省的西南部地區。地理坐標為129°46′55″～130°32′27″ E，43°49′57″～44°35′27″ N。穆棱林業局有8個林場。在整個黑龍江省來說，穆棱地區的氣溫稍微偏高，氣候類型屬于中溫帶大陸性季風氣候，夏季干燥炎熱，總是會出現很高的溫度，但炎熱持續的時間比較短暫。在該地區，降水時間主要集中在6～8月份。

為了使數據比較有代表性，在研究區域內不同立地條件下進行標準地的設置，對于標準地內的所有林木都要進行每木檢尺，測定每株林木的胸徑、樹高、枝下高等林木因子。依據所有林木的測量情況進行標準木的選取。該研究所用數據來自于帽兒山林場和穆棱林場，其中包括帽兒山2007年調查的45塊標準地和2008年調查的6塊標準地以及穆棱2008年調查的68塊標準地。

1.2研究方法

1.2.1直徑分布模型。根據以往的研究成果，正態分布、指數分布、Weibull分布、β分布以及SB分布都在一定范圍內得到了廣泛的應用，對所有數據進行擬合和檢驗，對直徑分布模型采用矩解法進行參數估計，個別異常點對參數擬合會造成很大的影響，會使模型產生離散，因此，參數估計之前把一些異常點都剔除。通過研究發現，Weibull分布和指數分布的通過率是非常大的，在該研究中選擇4塊標準地進行分析。

采用矩法估計對其進行研究：

(1)正態分布

正態分布的概率密度函數為：

(Ⅰ)

(2)Weibull分布

其概率密度函數為：

(Ⅱ)

式中：a為位置參數，一般指林分的最小直徑Dmin；b為比度參數；c為形狀參數；D為林木直徑。

(3)β分布

其概率密度函數為：

(Ⅲ)

式中：a,b分別為徑級的上限和下限；α,β為形狀參數。

(4)SB分布

其概率密度函數為：

(Ⅳ)

而：

式中：ε為徑級分布下限，ε=Dmin；λ為范圍參數，λ=Dmax-Dmin；σ為峰度參數；γ為偏度參數。

(5)指數分布

二參數負指數分布概率密度函數為：

(Ⅴ)

式中：a為位置參數(直徑分布最小徑階下限值)；b為尺度參數。

(Ⅵ)

1.2.2參數的再參數化。分別研究每一模型參數與數據的關系，選出最優的模型，分別探究各個參數的變化規律，通過描述其散點圖進而探究其變化規律。因此，模型的再參數化是構建預估模型的一個過程。

表1　各種分布函數擬合檢驗結果

2結果與分析

2.1直徑分布各模型的擬合檢驗結果將5種函數對各樣地直徑分布擬合檢驗結果列于表1。通過對5種函數進行χ2檢驗，可以得出：指數分布和Weibull分布接收率較大，擬合較好。指數分布優于Weibull分布。部分樣地的指數分布和Weibull分布擬合效果見圖1～4。

2.2指數分布與林分調查因子間的關系由指數分布的數據處理結果(表2)可以看出：指數分布的2個參數a，b中，a的值(除YH-7外)均為5，只有b的值在變化，通過對b值與林分調查因子之間的散點圖可知b值與林分平均直徑(DG)相關性顯著，而與其他因子均無明顯的相關性。選擇對數函數、指數函數、冪函數擬合b值和林分平均直徑(DG)關系，擬合效果圖(圖5)和擬合結果如下：

對數函數：b=-0.151ln(DG)+0.536

(Ⅶ)

其中：R2=0.921 4。

指數函數：b=0.410 3e-0.080 5DG

(Ⅷ)

其中：R2=0.935

冪函數：b=5.671 2DG-1.141 5

(Ⅸ)

其中：R2=0.957 3

通過對圖5分析可知：冪函數分布的參數b與林分平均直徑的相關系數最大，擬合效果最優。又因為直徑分布的指數函數為：f(D)=be-b(D-a)D>a，所以可得直徑分布函數為：

f(D)=5.671 2DG-1.415exp(-5.671 2DG-1.415(D-5))，D>5

(Ⅹ)

式中：D為林木的胸徑；DG為林分的平均直徑。

3 結論

該研究利用帽兒山2007年調查的45塊標準地和2008年調查的6塊標準地以及穆棱2008年調查的68塊標準地，選擇5種直徑分布函數進行擬合。

(1)該研究為了探究直徑分布的規律，對數據進行非常詳細的分析，對一些異常的數據進行了剔除。總體上，該文采用5種分布函數進行研究，這5個分布函數分別是指數分布、Weibull分布、SB分布、β分布和正態分布。結果表明，指數分布和Weibull分布都比較合適，但總體上指數分布的擬合效果比較好。

表2　指數分布擬合的參數 a、 b值

(2)在指數分布中，a值總體上不變，為5，只有參數b的數值會不斷變化。

(3)對參數b與林分平均直徑的散點圖進行分析，分別選擇對數函數、指數函數、冪函數進行分析，結果顯示對數函數比較好。

參考文獻

[1] 胡文力.長白山過伐林區云冷杉針闊混交林分結構的研究[D].北京：北京林業大學，2003.

[2] 李毅,孫雪新,康向陽.甘肅胡楊林分結構的研究[J].干旱區資源與環境,1994,8(3):88-95.

[3] 陳東來,秦淑英.山楊天然林林分結構的研究[J].河北農業大學學報,1994,17(1):36-43.

[4] 孟憲宇.測樹學[M].北京:中國林業出版社,1996.

[5] SANNIKOVA N S.Microecosystem analysis of the structure and functions of forest biogeocenoses[J].Russian Journal of Ecology, 2003,34(2):80-85.

[6] 邵青還.第二次林業革命——“接近自然的林業”在中歐興起[J].世界林業研究,1991,4(4):8-15.

[7] 沈國舫.現代高效持續林業——中國林業發展道路的抉擇[J].林業經濟,1998(4):1-8.

[8] 邵青還.德國的林業保護政策及其評價[C]//中國林業如何走向21世紀.北京:中國林業出版社,1995：53-60.

[9] 李芳東,李宗然,周道順,等.蘭考泡桐林分結構規律研究[J].林業科學研究，1996(2):114-120.

[10] 黃清麟.福建青岡萌芽林分結構及生產力的研究[J].福建林學院學報,1995,15(2):107-111.

[11] 陳昌雄,陳平留.閩北天然異齡林林分結構規律的研究[J].福建林業科技,1997,24(4):1-4.

[12] WAYAN C.Zipperer Species composition and structure of regenerated and remnant forest patches within an urban landscape[J].Urban Ecosystems,2002,6:271-290.

[13] CHOKKALINGAM U，WHITE A.Structure and spatial patterns of trees in old-growth northern hardwood and mixed forests of northern Maine[J].Plant Ecology,2001，156(2)：139-160.

[14] NISHIMURA N，HARA T，MIURA M.Tree competition and species coexistence in a warm-temperate old-growth evergreen broad-leaved forest in Japan[J].Plant Ecology,2002,164:235-248.