嚴迎冬，談為民，蔣 彪，付世華

（湖北省當陽市林業(yè)局，湖北 當陽444100）

1 引言

在處理盜伐、濫伐林木案件時，往往是盜伐、濫伐的林木已滅失，無法查證，案發(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è)案件帶來極大的方便，從而降低了辦案成本。

2 材料與方法

2.1 自然概況

當陽市地處大巴山脈東麓，荊山山脈以南，鄂西山地向江漢平原過渡地帶，境內(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 調(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］計 算出當陽市馬尾松、闊葉樹地徑一元材積表。

3 結(jié)果與分析

3.1 地徑與胸徑的關(guān)系

用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。

3.2 回歸方程F值檢驗

對馬尾松抽取的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 闊葉樹回歸方程方差分析

3.3 地徑材積表編制

用地徑數(shù)據(jù)計算出當陽市馬尾松、闊葉樹地徑一元材積表，見表3、表4。

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

表4 當陽市闊葉樹地徑一元材積

4 結(jié)語

（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.