嚴迎冬,談為民,蔣 彪,付世華
(湖北省當陽市林業局,湖北 當陽444100)
在處理盜伐、濫伐林木案件時,往往是盜伐、濫伐的林木已滅失,無法查證,案發現場只能調查到被伐林木的樹樁,而采伐木的胸徑和樹高難以現場求證,也就不能準確地測定被盜伐、濫伐林木的蓄積。而盜伐、濫伐蓄積的多少,是辦案的主要證據。在這種情況下,如何計算涉案林木蓄積就是個難題。以往辦理盜伐、濫伐林木案件,每一個案件都要在現場附近量測50~100株的樣木作為推算地徑和胸徑關系的依據[1],工作重復,且工作量大。如果能夠找到地徑和胸徑的關系,通過地徑計算出相對應的胸徑,那么利用一元立木材積表計算林木蓄積就給林業案件辦理提供了一個有力的科學依據,給林業部門辦理林業案件帶來極大的方便,從而降低了辦案成本。
當陽市地處大巴山脈東麓,荊山山脈以南,鄂西山地向江漢平原過渡地帶,境內地貌以丘陵、崗地為主。
當陽市樹種資源十分豐富,有木本植物76科320余種,但構成該市森林資源的主要樹種是馬尾松和櫟類、化香等闊葉樹,以馬尾松為優勢樹種的森林面積占森林總面積的33.10%,蓄積占森林總蓄積的48.25%,櫟類等闊葉樹占森林總面積的42.82%,占森林總蓄積的40.51%,這兩大優勢樹種的面積就占了森林面積的75.92%,占森林蓄積的88.76%,其它樹種所占比例很小,所以調查材料以馬尾松和闊葉樹為主。
2.2.1 外業調查
從2007年2月開始至2015年5月,對林業案件現場附近有代表性的馬尾松、闊葉樹活立木用圍尺進行了地徑和胸徑的精確測定,共調查922組馬尾松、837組闊葉樹的地徑和胸徑。調查的馬尾松最小胸徑6.2cm,最大胸徑54.4cm,調查的闊葉樹最小胸徑4.6cm,最大胸徑73.3cm,調查地點遍布全市6個鄉鎮、3個辦事處和5個國有林場,真實反映了當陽市馬尾松、闊葉樹的實際生長現狀。
2.2.2 回歸分析方法
從922組馬尾松和837組闊葉樹數據中,依檢尺順序按間隔8的周期各抽取100組數據作為馬尾松、闊葉樹回歸方程F值的檢驗數據。根據抽樣余下的822組馬尾松和737組闊葉樹地徑與胸徑的調查數據繪制散點圖,觀察散點圖分布后研判,地徑與胸徑有一元線性回歸關系,采用回歸方程式y=a+bx進行回歸分析,用Excel數據分析工具求出回歸常數項a、b的值和回歸系數r[2]。
2.2.3 回歸方程檢驗方法
對馬尾松、闊葉樹各抽取的100組實測胸徑數值與回歸分析方程式推算的對應胸徑數據,用F值對回歸方程進行檢驗。
2.2.4 地徑材積表編制
用馬尾松、闊葉樹地徑計算出胸徑,根據鄂東馬尾松一元立木材積公式v=0.000060049144×(-0.13210336+0.97887017×d)^1.8719753×(24.269237-591.97756/(24+d)^0.97180232和鄂東闊葉一元立木材積公式v=0.000050479054×(-0.21700621+0.98481055×d)^1.9085054×(17.823386-272.42014/(17+d)^0.99076507[3]計 算出當陽市馬尾松、闊葉樹地徑一元材積表。
用Excel數據分析工具求得馬尾松的回歸方程的常數項a=-1.553728048,b=0.872731373,回歸系數r=0.975173625,回歸方程式為y(胸徑)=0.872731373×x-1.553728048;闊葉樹回歸方程的常數項a=-1.027634523,b=0.884994409,回歸系數r=0.963926604,回歸方程式為y(胸徑)=0.884994409×x-1.027634523。
對馬尾松抽取的100組實測胸徑數值與回歸分析方程式推算的對應胸徑數據,用Excel數據分析進行單因素方差分析,按照可靠性0.99計算,F=19.470 73>F0.01(f1=1,f2=198)=6.76,說明馬尾松的實測胸徑與推算胸徑的關系顯著。詳見表1。

表1 馬尾松回歸方程方差分析
對闊葉樹抽取的100組實測胸徑數值與回歸分析方程式推算的對應胸徑數據,用Excel數據分析進行單因素方差分析,按照可靠性0.99計算,F=9.376 379>F0.01(f1=1,f2=198)=6.76,說明闊葉樹的地徑與胸徑的關系顯著。說明闊葉樹的實測胸徑與推算胸徑的關系顯著,詳見表2。

表2 闊葉樹回歸方程方差分析
用地徑數據計算出當陽市馬尾松、闊葉樹地徑一元材積表,見表3、表4。

表3 當陽市馬尾松地徑一元材積

表4 當陽市闊葉樹地徑一元材積
(1)馬尾松、闊葉樹的地徑與胸徑關系顯著,回歸方程式可靠性高,可用于林業生產實際工作。
(2)在處理濫伐、盜伐林木案件時,通過采伐現場測定的伐樁直徑查馬尾松、闊葉樹地徑材積表可直接計算林木蓄積,降低了辦案成本,提高了調查效率。
(3)地徑材積表編制完成之后,此表數據在用于盜伐、濫伐林木案件時,經驗證數據是較為適用的。但是為編制數表所實測的樣木是生長在自然界中的活立木,同樹種、同齡級的樹木,受立地條件的不同,其材積生長量都不相同,推算出胸徑值相對應的材積表就可能出現偏差,特殊立地或特殊樹種在引用此數據時應審慎。為解決此問題,將再選取不同立地、林分的樣木,擴大實測數據量,參與計算,使回歸方程式更接近于該市林業生產實際。
[1]何邦亮,鄧安桂,謝永勝,等 .重慶市馬尾松伐樁一元材積表編制研究[J].重慶林業科技,2012,95(5):1~7.
[2]孟憲宇 .測樹學[M].北京:中國林業出版社,2006.
[3]陳 華,龔萬鵬,王萬政,等 .宜昌市馬尾松地徑一元材積表編制[J].湖北林業科技,2007(5):25~26.