當(dāng)陽(yáng)市馬尾松、闊葉樹(shù)地徑材積表編制

嚴(yán)迎冬，談為民，蔣 彪，付世華

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

1 引言

在處理盜伐、濫伐林木案件時(shí)，往往是盜伐、濫伐的林木已滅失，無(wú)法查證，案發(fā)現(xiàn)場(chǎng)只能調(diào)查到被伐林木的樹(shù)樁，而采伐木的胸徑和樹(shù)高難以現(xiàn)場(chǎng)求證，也就不能準(zhǔn)確地測(cè)定被盜伐、濫伐林木的蓄積。而盜伐、濫伐蓄積的多少，是辦案的主要證據(jù)。在這種情況下，如何計(jì)算涉案林木蓄積就是個(gè)難題。以往辦理盜伐、濫伐林木案件，每一個(gè)案件都要在現(xiàn)場(chǎng)附近量測(cè)50～100株的樣木作為推算地徑和胸徑關(guān)系的依據(jù)［1］，工作重復(fù)，且工作量大。如果能夠找到地徑和胸徑的關(guān)系，通過(guò)地徑計(jì)算出相對(duì)應(yīng)的胸徑，那么利用一元立木材積表計(jì)算林木蓄積就給林業(yè)案件辦理提供了一個(gè)有力的科學(xué)依據(jù)，給林業(yè)部門辦理林業(yè)案件帶來(lái)極大的方便，從而降低了辦案成本。

2 材料與方法

2.1 自然概況

當(dāng)陽(yáng)市地處大巴山脈東麓，荊山山脈以南，鄂西山地向江漢平原過(guò)渡地帶，境內(nèi)地貌以丘陵、崗地為主。

當(dāng)陽(yáng)市樹(shù)種資源十分豐富，有木本植物76科320余種，但構(gòu)成該市森林資源的主要樹(shù)種是馬尾松和櫟類、化香等闊葉樹(shù)，以馬尾松為優(yōu)勢(shì)樹(shù)種的森林面積占森林總面積的33.10%，蓄積占森林總蓄積的48.25%，櫟類等闊葉樹(shù)占森林總面積的42.82%，占森林總蓄積的40.51%，這兩大優(yōu)勢(shì)樹(shù)種的面積就占了森林面積的75.92%，占森林蓄積的88.76%，其它樹(shù)種所占比例很小，所以調(diào)查材料以馬尾松和闊葉樹(shù)為主。

2.2 調(diào)查方法

2.2.1 外業(yè)調(diào)查

從2007年2月開(kāi)始至2015年5月，對(duì)林業(yè)案件現(xiàn)場(chǎng)附近有代表性的馬尾松、闊葉樹(shù)活立木用圍尺進(jìn)行了地徑和胸徑的精確測(cè)定，共調(diào)查922組馬尾松、837組闊葉樹(shù)的地徑和胸徑。調(diào)查的馬尾松最小胸徑6.2cm，最大胸徑54.4cm，調(diào)查的闊葉樹(shù)最小胸徑4.6cm，最大胸徑73.3cm，調(diào)查地點(diǎn)遍布全市6個(gè)鄉(xiāng)鎮(zhèn)、3個(gè)辦事處和5個(gè)國(guó)有林場(chǎng)，真實(shí)反映了當(dāng)陽(yáng)市馬尾松、闊葉樹(shù)的實(shí)際生長(zhǎng)現(xiàn)狀。

2.2.2 回歸分析方法

從922組馬尾松和837組闊葉樹(shù)數(shù)據(jù)中，依檢尺順序按間隔8的周期各抽取100組數(shù)據(jù)作為馬尾松、闊葉樹(shù)回歸方程F值的檢驗(yàn)數(shù)據(jù)。根據(jù)抽樣余下的822組馬尾松和737組闊葉樹(shù)地徑與胸徑的調(diào)查數(shù)據(jù)繪制散點(diǎn)圖，觀察散點(diǎn)圖分布后研判，地徑與胸徑有一元線性回歸關(guān)系，采用回歸方程式y(tǒng)=a＋bx進(jìn)行回歸分析，用Excel數(shù)據(jù)分析工具求出回歸常數(shù)項(xiàng)a、b的值和回歸系數(shù)r［2］。

2.2.3 回歸方程檢驗(yàn)方法

對(duì)馬尾松、闊葉樹(shù)各抽取的100組實(shí)測(cè)胸徑數(shù)值與回歸分析方程式推算的對(duì)應(yīng)胸徑數(shù)據(jù)，用F值對(duì)回歸方程進(jìn)行檢驗(yàn)。

2.2.4 地徑材積表編制

用馬尾松、闊葉樹(shù)地徑計(jì)算出胸徑，根據(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］計(jì) 算出當(dāng)陽(yáng)市馬尾松、闊葉樹(shù)地徑一元材積表。

3 結(jié)果與分析

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

用Excel數(shù)據(jù)分析工具求得馬尾松的回歸方程的常數(shù)項(xiàng)a=-1.553728048，b=0.872731373，回歸系數(shù)r=0.975173625，回歸方程式為y（胸徑）=0.872731373×x-1.553728048；闊葉樹(shù)回歸方程的常數(shù)項(xiàng)a=-1.027634523，b=0.884994409，回歸系數(shù)r=0.963926604，回歸方程式為y（胸徑）=0.884994409×x-1.027634523。

3.2 回歸方程F值檢驗(yàn)

對(duì)馬尾松抽取的100組實(shí)測(cè)胸徑數(shù)值與回歸分析方程式推算的對(duì)應(yīng)胸徑數(shù)據(jù)，用Excel數(shù)據(jù)分析進(jìn)行單因素方差分析，按照可靠性0.99計(jì)算，F(xiàn)=19.470 73＞F0.01（f1=1，f2=198）=6.76，說(shuō)明馬尾松的實(shí)測(cè)胸徑與推算胸徑的關(guān)系顯著。詳見(jiàn)表1。

表1 馬尾松回歸方程方差分析

對(duì)闊葉樹(shù)抽取的100組實(shí)測(cè)胸徑數(shù)值與回歸分析方程式推算的對(duì)應(yīng)胸徑數(shù)據(jù)，用Excel數(shù)據(jù)分析進(jìn)行單因素方差分析，按照可靠性0.99計(jì)算，F(xiàn)=9.376 379＞F0.01（f1=1，f2=198）=6.76，說(shuō)明闊葉樹(shù)的地徑與胸徑的關(guān)系顯著。說(shuō)明闊葉樹(shù)的實(shí)測(cè)胸徑與推算胸徑的關(guān)系顯著，詳見(jiàn)表2。

表2 闊葉樹(shù)回歸方程方差分析

3.3 地徑材積表編制

用地徑數(shù)據(jù)計(jì)算出當(dāng)陽(yáng)市馬尾松、闊葉樹(shù)地徑一元材積表，見(jiàn)表3、表4。

表3 當(dāng)陽(yáng)市馬尾松地徑一元材積

表4 當(dāng)陽(yáng)市闊葉樹(shù)地徑一元材積

4 結(jié)語(yǔ)

（1）馬尾松、闊葉樹(shù)的地徑與胸徑關(guān)系顯著，回歸方程式可靠性高，可用于林業(yè)生產(chǎn)實(shí)際工作。

（2）在處理濫伐、盜伐林木案件時(shí)，通過(guò)采伐現(xiàn)場(chǎng)測(cè)定的伐樁直徑查馬尾松、闊葉樹(shù)地徑材積表可直接計(jì)算林木蓄積，降低了辦案成本，提高了調(diào)查效率。

（3）地徑材積表編制完成之后，此表數(shù)據(jù)在用于盜伐、濫伐林木案件時(shí)，經(jīng)驗(yàn)證數(shù)據(jù)是較為適用的。但是為編制數(shù)表所實(shí)測(cè)的樣木是生長(zhǎng)在自然界中的活立木，同樹(shù)種、同齡級(jí)的樹(shù)木，受立地條件的不同，其材積生長(zhǎng)量都不相同，推算出胸徑值相對(duì)應(yīng)的材積表就可能出現(xiàn)偏差，特殊立地或特殊樹(shù)種在引用此數(shù)據(jù)時(shí)應(yīng)審慎。為解決此問(wèn)題，將再選取不同立地、林分的樣木，擴(kuò)大實(shí)測(cè)數(shù)據(jù)量，參與計(jì)算，使回歸方程式更接近于該市林業(yè)生產(chǎn)實(shí)際。

［1］何邦亮，鄧安桂，謝永勝，等 .重慶市馬尾松伐樁一元材積表編制研究［J］.重慶林業(yè)科技，2012，95（5）：1～7.

［2］孟憲宇 .測(cè)樹(shù)學(xué)［M］.北京：中國(guó)林業(yè)出版社，2006.

［3］陳 華，龔萬(wàn)鵬，王萬(wàn)政，等 .宜昌市馬尾松地徑一元材積表編制［J］.湖北林業(yè)科技，2007（5）：25～26.