Weibull函數在杉木人工林直徑分布中的研究

陳晨+劉光武

摘要：用Weibull 分布模型擬合林分直徑分布，用最大似然估計法求其模型參數，再求出各徑階的理論株數值，并進行了x2檢驗，擬合精度達到84.91。研究結果表明：Weibull分布能較好地模擬林分直徑分布。

關鍵詞：杉木；人工林；Weibull函數；直徑分布

中圖分類號：S791.28文獻標識碼：A文章編號：16749944（2017）03005902

1引言

林分直徑分布是指在林分內各種大小直徑林木按徑階分配的狀態，它是基本的林分結構[1]。林分直徑分布影響樹高、干形、材積、材種及樹冠等因子的變化，同時，研究林分直徑分布規律對森林經營技術、編制經營數表、預估林分生長與收獲模型及林分調查都有十分重要的意義。目前，關于林分直徑分布規律的研究有很多，但普遍都是針對油松[2，3]、落葉松[4，5]、馬尾松[6，7]等樹種的，而對杉木的研究還很少。

在擬合直徑分布時常根據直徑分布的具體形態特征和變化規律，選用不同的概率分布函數，常用的分布函數有正態分布、對數正太分布、Sb分布等，本研究用Weibull函數來擬合杉木人工林的直徑分布。

2材料與方法

2.1材料來源

供研究的材料來源于河南省信陽市雞公山自然保護區某固定標準地，該保護區位于豫鄂兩省交界處，面積287 km2，同時也是我國淮河與長江兩大水系的分水嶺，中國南北氣候的分界線，北亞熱帶和暖溫帶的過渡地帶。特殊的地理位置與氣候特征，造就了此地的生物多樣性。

2.2數據篩選

研究結果表明，不論樹種、年齡、密度和立地條件如何，其林分平均直徑（Dg）在株數累計分布曲線上的位置大致在55%～64%，一般在60%處[1]。從表1所示資料可以得出，相對直徑等于1.0（即林分平均直徑）所對應的株數累計百分數在57.0%處，說明所用固定標準地林分未遭受重大自然災害、病蟲害或人工撫育破壞，樣本數據符合人工林直徑分布規律，可以用來建立模型。

2.3研究方法

筆者用Weibull分布模型擬合林分直徑分布。Weibull分布的密度函為：

式中：a為位置參數；b為尺度參數；c為形狀參數，x為徑階中值。參數a一般定為林分直徑最小徑階的下限值，b和c的求解以最大似然估計法最精確。這樣，根據最大似然原理得到：

[∑ki=1fi·xci·In（xi）/（∑ki=1fi·xci）]-1c=1n[∑ki=1fi·In（xi）]（2）

b=[1n（∑1nfi·xci）1/c]（3）

式中：n=∑kifi；k為徑階個數；xi為第i徑階中值；fi為第i徑階內林木株數；n為總株數。

利用最大似然估計法求出三參數a、b、c的值之后，再利用下式求出標準地林木直徑的Weibull分布各徑階的理論株數值，并進行X2檢驗。

ni=N·W·（ab）（xi-ab）c-1exp[-（xi-ab）c]（4）

式中：ni為第i徑階內理論株數值；N為林木總株數；W為徑階距；xi為第i徑階中值，其它符號同前。

3擬合結果

由表2可知， Weibull分布模型的擬合精度為84.91%，這說明Weibull函數具有較好的擬合效果，能4結論

用雞公山林場固定標準地直徑分布數據為研究材料，利用Weibull分布函數建立模型，再求出標準地林木直徑的Weibull分布各徑階的理論株數值，并進行了x2檢驗。擬合精度達到84.91，研究證明Weibull分布模型在擬合林分直徑分布時，靈活性高、適應性強，能較好地擬合實直徑分布。

參考文獻：

[1]孟憲宇.測樹學（第二版）[M].北京：中國林業出版社，1996：66～98.

[2]刁軍，雷相東，國紅，等.北京西山油松人工林結構特征研究[J].河北林果研究，2009，24（3）：225～230.

[3]羅瑞平，亢新剛.黃龍山林區天然油松林直徑結構規律[J].林業調查規劃，2006，31（3）：86～88.

[4]王方，鄭小賢，蔣桂娟.不同密度落葉松人工林直徑結構研究[J].內蒙古林業調查設計，2011，34（5）：9～11.

[5]王香春，張秋良，春蘭，等.大青山落葉松人工林直徑分布規律的研究[J].山東農業大學學報（自然科學版），2011，42（3）：349～355.

[6]黃家榮，孟憲宇，關毓秀.馬尾松人工林直徑分布神經網絡模型研究[J].北京林業大學學報，2006，28（1）：28～31.

[7]黃慶豐.馬尾松林分直徑結構規律的研究[J].安徽農業大學學報，1998，25（1）：23～25.

收稿日期：20170111

作者簡介：賴寶乾（1964—），男，工程師，主要從事森林資源調查規劃設計技術研究。2017年2月Journal of Green Science and Technology第3期