嚴迎冬,談為民,蔣 彪,付世華
(湖北省當陽市林業(yè)局,湖北 當陽444100)
在處理盜伐、濫伐林木案件時,往往是盜伐、濫伐的林木已滅失,無法查證,案發(fā)現(xiàn)場只能調(diào)查到被伐林木的樹樁,而采伐木的胸徑和樹高難以現(xiàn)場求證,也就不能準確地測定被盜伐、濫伐林木的蓄積。而盜伐、濫伐蓄積的多少,是辦案的主要證據(jù)。在這種情況下,如何計算涉案林木蓄積就是個難題。以往辦理盜伐、濫伐林木案件,每一個案件都要在現(xiàn)場附近量測50~100株的樣木作為推算地徑和胸徑關(guān)系的依據(jù)[1],工作重復(fù),且工作量大。如果能夠找到地徑和胸徑的關(guān)系,通過地徑計算出相對應(yīng)的胸徑,那么利用一元立木材積表計算林木蓄積就給林業(yè)案件辦理提供了一個有力的科學(xué)依據(jù),給林業(yè)部門辦理林業(yè)案件帶來極大的方便,從而降低了辦案成本。
當陽市地處大巴山脈東麓,荊山山脈以南,鄂西山地向江漢平原過渡地帶,境內(nèi)地貌以丘陵、崗地為主。
當陽市樹種資源十分豐富,有木本植物76科320余種,但構(gòu)成該市森林資源的主要樹種是馬尾松和櫟類、化香等闊葉樹,以馬尾松為優(yōu)勢樹種的森林面積占森林總面積的33.10%,蓄積占森林總蓄積的48.25%,櫟類等闊葉樹占森林總面積的42.82%,占森林總蓄積的40.51%,這兩大優(yōu)勢樹種的面積就占了森林面積的75.92%,占森林蓄積的88.76%,其它樹種所占比例很小,所以調(diào)查材料以馬尾松和闊葉樹為主。
2.2.1 外業(yè)調(diào)查
從2007年2月開始至2015年5月,對林業(yè)案件現(xiàn)場附近有代表性的馬尾松、闊葉樹活立木用圍尺進行了地徑和胸徑的精確測定,共調(diào)查922組馬尾松、837組闊葉樹的地徑和胸徑。調(diào)查的馬尾松最小胸徑6.2cm,最大胸徑54.4cm,調(diào)查的闊葉樹最小胸徑4.6cm,最大胸徑73.3cm,調(diào)查地點遍布全市6個鄉(xiāng)鎮(zhèn)、3個辦事處和5個國有林場,真實反映了當陽市馬尾松、闊葉樹的實際生長現(xiàn)狀。
2.2.2 回歸分析方法
從922組馬尾松和837組闊葉樹數(shù)據(jù)中,依檢尺順序按間隔8的周期各抽取100組數(shù)據(jù)作為馬尾松、闊葉樹回歸方程F值的檢驗數(shù)據(jù)。根據(jù)抽樣余下的822組馬尾松和737組闊葉樹地徑與胸徑的調(diào)查數(shù)據(jù)繪制散點圖,觀察散點圖分布后研判,地徑與胸徑有一元線性回歸關(guān)系,采用回歸方程式y(tǒng)=a+bx進行回歸分析,用Excel數(shù)據(jù)分析工具求出回歸常數(shù)項a、b的值和回歸系數(shù)r[2]。
2.2.3 回歸方程檢驗方法
對馬尾松、闊葉樹各抽取的100組實測胸徑數(shù)值與回歸分析方程式推算的對應(yīng)胸徑數(shù)據(jù),用F值對回歸方程進行檢驗。
2.2.4 地徑材積表編制
用馬尾松、闊葉樹地徑計算出胸徑,根據(jù)鄂東馬尾松一元立木材積公式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數(shù)據(jù)分析工具求得馬尾松的回歸方程的常數(shù)項a=-1.553728048,b=0.872731373,回歸系數(shù)r=0.975173625,回歸方程式為y(胸徑)=0.872731373×x-1.553728048;闊葉樹回歸方程的常數(shù)項a=-1.027634523,b=0.884994409,回歸系數(shù)r=0.963926604,回歸方程式為y(胸徑)=0.884994409×x-1.027634523。
對馬尾松抽取的100組實測胸徑數(shù)值與回歸分析方程式推算的對應(yīng)胸徑數(shù)據(jù),用Excel數(shù)據(jù)分析進行單因素方差分析,按照可靠性0.99計算,F(xiàn)=19.470 73>F0.01(f1=1,f2=198)=6.76,說明馬尾松的實測胸徑與推算胸徑的關(guān)系顯著。詳見表1。
表1 馬尾松回歸方程方差分析
對闊葉樹抽取的100組實測胸徑數(shù)值與回歸分析方程式推算的對應(yīng)胸徑數(shù)據(jù),用Excel數(shù)據(jù)分析進行單因素方差分析,按照可靠性0.99計算,F(xiàn)=9.376 379>F0.01(f1=1,f2=198)=6.76,說明闊葉樹的地徑與胸徑的關(guān)系顯著。說明闊葉樹的實測胸徑與推算胸徑的關(guān)系顯著,詳見表2。
表2 闊葉樹回歸方程方差分析
用地徑數(shù)據(jù)計算出當陽市馬尾松、闊葉樹地徑一元材積表,見表3、表4。
表3 當陽市馬尾松地徑一元材積
表4 當陽市闊葉樹地徑一元材積
(1)馬尾松、闊葉樹的地徑與胸徑關(guān)系顯著,回歸方程式可靠性高,可用于林業(yè)生產(chǎn)實際工作。
(2)在處理濫伐、盜伐林木案件時,通過采伐現(xiàn)場測定的伐樁直徑查馬尾松、闊葉樹地徑材積表可直接計算林木蓄積,降低了辦案成本,提高了調(diào)查效率。
(3)地徑材積表編制完成之后,此表數(shù)據(jù)在用于盜伐、濫伐林木案件時,經(jīng)驗證數(shù)據(jù)是較為適用的。但是為編制數(shù)表所實測的樣木是生長在自然界中的活立木,同樹種、同齡級的樹木,受立地條件的不同,其材積生長量都不相同,推算出胸徑值相對應(yīng)的材積表就可能出現(xiàn)偏差,特殊立地或特殊樹種在引用此數(shù)據(jù)時應(yīng)審慎。為解決此問題,將再選取不同立地、林分的樣木,擴大實測數(shù)據(jù)量,參與計算,使回歸方程式更接近于該市林業(yè)生產(chǎn)實際。
[1]何邦亮,鄧安桂,謝永勝,等 .重慶市馬尾松伐樁一元材積表編制研究[J].重慶林業(yè)科技,2012,95(5):1~7.
[2]孟憲宇 .測樹學(xué)[M].北京:中國林業(yè)出版社,2006.
[3]陳 華,龔萬鵬,王萬政,等 .宜昌市馬尾松地徑一元材積表編制[J].湖北林業(yè)科技,2007(5):25~26.