胡小月 蔣海平 濮厚平

摘要：為了解決林業(yè)生產(chǎn)經(jīng)營活動和行政執(zhí)法工作中只知樹木重量不知樹木材積的難題，筆者研究了馬鞍山市江東片區(qū)國外松樹干材積與重量的相關(guān)關(guān)系，獲得國外松材積（Y）與重量（X）的關(guān)系模型，以0.95的可靠性檢驗了其顯著性，并對該關(guān)系模型進(jìn)行了驗證，驗證精度達(dá)95%。從求算的結(jié)果可以看出，馬鞍山市江東片區(qū)國外松樹干材積與濕重總體線性回歸的關(guān)系是顯著性的。

關(guān)鍵詞：國外松；樹干重量與材積；關(guān)系模型

文章編號：1674-7437（2022）12-0166-03? ? ? 中國圖書分類號：S791.246? ? ?文章標(biāo)識碼：B

如何從樹木重量得出樹木材積是在林業(yè)工作尤其是行政執(zhí)法的過程中會時常遇到的問題。例如，根據(jù)相關(guān)司法解釋，盜伐濫伐林木的立案、查處、量刑定罪，是以所盜伐濫伐林木的材積為依據(jù)的[1]；有時，滿足某個木材加工企業(yè)對林木的需求依據(jù)是樹木材積而非重量。從樹木重量得出樹木材積，國家和安徽省均沒有相關(guān)的數(shù)表可供查找，也沒有可以通用的經(jīng)驗公式可供求算，這就給林業(yè)生產(chǎn)經(jīng)營活動和行政執(zhí)法工作帶來一定的困難。因此，采用某種科學(xué)方法，在調(diào)查的區(qū)域隨機抽取一定樣地，并在樣地上隨機抽取若干樣木組成樣本，獲得高精度的樹干材積與重量的關(guān)系模型，并在此模型公式基礎(chǔ)上編制出樹干濕重與材積對應(yīng)值數(shù)表，進(jìn)而得出的研究結(jié)果既能為研究區(qū)域林政執(zhí)法提供一定的科學(xué)依據(jù)，也會為該區(qū)域林業(yè)生產(chǎn)經(jīng)營管理提供一定的基礎(chǔ)數(shù)據(jù)?；谝陨峡紤]，結(jié)合工作實際，筆者在馬鞍山市江東的一縣三區(qū)，開展了國外松樹干重量與材積的關(guān)系研究。

1? ?研究區(qū)域概況

馬鞍山市地處長江下游、安徽東部，下轄三縣三區(qū)，國土總面積4 049km2，屬北亞熱帶季風(fēng)氣候，四季分明，雨量充沛，植被豐富。據(jù)馬鞍山市第三次全國國土調(diào)查主要數(shù)據(jù)公報顯示，全市林地73 653.19hm2；其中，喬木林地49 775.04hm2，占比67.58%；竹林地

1 204.42hm2，占比1.64%；灌木林地5 983.97hm2，占比8.12%；其它林地16 689.76hm2，占22.66%[2]。馬鞍山市長江以東（即江東片區(qū)）下轄當(dāng)涂縣、花山區(qū)、雨山區(qū)及博望區(qū)，主要喬木樹種有杉、松、硬闊、軟闊、國外松、楊樹等。其中，國外松是江東片區(qū)重要造林樹種，也是構(gòu)成江東片區(qū)森林資源的主要樹種。

2? ?研究樣本的抽取與調(diào)查

2.1? ?樣地、樣木的隨機抽取

在馬鞍山市江東片區(qū)的花山區(qū)、雨山區(qū)、博望區(qū)及當(dāng)涂縣國外松林分中，隨機抽取了15塊樣地，并在這15塊樣地上隨機抽取了40株國外松（火炬松拉丁學(xué)名：Pinus taeda L，濕地松拉丁學(xué)名：Pinus elliottii）樣木組成樣本，用于本次國外松樹干重量與材積的關(guān)系研究。其中，幼齡林樣木8株、中齡林樣木1株、近熟林樣木21株、成熟林樣木10株。按地域分，花山區(qū)樣地4塊、樣木11株，雨山區(qū)樣地3塊、樣木15株，博望區(qū)樣地2塊、樣木3株，當(dāng)涂縣樣地6塊、樣木11株。

2.2? ?樣木樹干的材積測算和濕重稱量

由3名林業(yè)技術(shù)人員組成的調(diào)查組，在當(dāng)?shù)貒辛謭龊袜l(xiāng)鎮(zhèn)、村的協(xié)助下，對隨機抽取的每一塊樣地上的每一株樣木，在測量其胸徑之后將其伐倒，伐樁高度均控制在3cm以下。將伐倒木打去枝椏后，測量其長度（即樹高）、底端直徑。按照1m一個區(qū)分段，從底端起將伐倒木區(qū)分為若干段，分別測量每一段的中央直徑，并測量不夠分段的稍頭長度及稍底直徑[3]。最后用磅秤秤量伐倒木的重量[4]（為濕重，這里用“Xi”表示）。

在室內(nèi)，根據(jù)現(xiàn)場記錄的伐倒木測量結(jié)果，按照中央斷面積區(qū)分求積法，對所有伐倒木的材積進(jìn)行求算，得出每株樣木的帶皮材積（這里用“Yi”表示）。

對在馬鞍山市江東片區(qū)國外松林分中隨機抽取的樣木樹干濕重秤量及材積[5]測算結(jié)果進(jìn)行了記錄（見表1）。

3? ?總體均方直線回歸模型研究與精度檢驗

3.1? ?樹干材積與濕重相關(guān)關(guān)系分析

將馬鞍山市江東片區(qū)的所有國外松看作一個總體。由數(shù)理統(tǒng)計[6]知曉，對于這個總體中的每株國外松樹干濕重，都有樹干材積的一個概率分布與其相對應(yīng)，即樹干材積與樹干濕重存在著回歸關(guān)系，這種回歸關(guān)系可以確定為一個樹干材積與樹干濕重的均方直線回歸，用直線回歸方程Yi=A+BXi來表示。其中，Y代表因變量-樹干材積，X代表自變量-樹干濕重。

在馬鞍山市江東片區(qū)國外松林分中，隨機抽取40株國外松樣木組成研究樣本，對總體均方直線回歸方程進(jìn)行估計，求算出A、B的估計值a、b，從而導(dǎo)算出所估計的江東片區(qū)國外松樹干材積與濕量的均方直線回歸模型。

3.2? ?國外松材積與濕重均方直線回歸模型導(dǎo)算

依據(jù)表1的相關(guān)數(shù)據(jù)，按照均方直線回歸方程參數(shù)求算公式，可以得出估計值a、b，見公式（1）（2）：

因此，依據(jù)上述樣本所估計的馬鞍山市江東片區(qū)國外松樹干濕重與材積之間均方直線回歸方程，見公式（3）：

y=0.0244+0.0011x

3.3? ?國外松材積與濕重直線回歸模型顯著性F檢驗

用國外松樹干材積與濕重均方直線回歸模型導(dǎo)算樣木所屬的馬鞍山市江東片區(qū)國外松總體，滿足線性、正態(tài)、等方差條件，抽樣方式滿足獨立的條件。因此，可以采用線性回歸顯著性的F檢驗方法，對國外松樹干材積與濕重總體線性回歸的顯著性進(jìn)行F檢驗（即方差分析）。根據(jù)上述導(dǎo)算國外松樹干材積與濕重均方線性回歸模型的樣本資料，以0.95的可靠性，求算出線性回歸模型顯著性F檢驗所需數(shù)值，并列出方差分析表（見表2）。

從求算的結(jié)果可以看出，F(xiàn)>F0.05 ，所以，從0.95的可靠性來檢驗，上述馬鞍山市江東片區(qū)國外松樹干材積與濕重總體線性回歸的關(guān)系是顯著性的。

3.4? ?模型驗證

隨機抽取未參加建模的10株國外松單株濕重，代入獲得的國外松材積與濕重關(guān)系模型，比較分析回歸材積與單株實際材積，驗證樹干重量與材積關(guān)系模型的適用性，以下為驗證結(jié)果（見表3）。

表3 模型驗證表

方差分析的相對誤差RS值在林業(yè)上的精度要求范圍之內(nèi)。

4? ?結(jié)論與討論

4.1? ?結(jié)論

綜上所述，隨機抽取樣地以及在樣地上隨機抽取若干樣木，組成馬鞍山市江東片區(qū)國外松樹干材積與濕重均方直線回歸模型研究樣本，其每株樣木的樹干濕重、材積均是獨立、隨機、正態(tài)和均方差的。因此，因變量-樹干材積（y）對于自變量-樹干濕重（x），存在著均方直線回歸關(guān)系。

根據(jù)樣本數(shù)據(jù)求算的樹干材積（y）與樹干濕重（x）的直線回歸方程，見公式（4），編制出了《馬鞍山市江東片區(qū)國外松濕重與材積對應(yīng)值數(shù)表》，供馬鞍山市江東片區(qū)的林業(yè)生產(chǎn)經(jīng)營活動和行政執(zhí)法工作參考（見表4）。

y=0.024 4+0.001 1x

4.2? ? 討論

本次調(diào)查研究了馬鞍山市江東片區(qū)國外松伐倒木的樹干濕重與材積之間的關(guān)系模型，范圍不夠全面，對當(dāng)?shù)亓謽I(yè)生產(chǎn)經(jīng)營活動和行政執(zhí)法工作的作用還十分有限。下一步需要開展如下研究：一是增加樣本數(shù)量，使導(dǎo)算的直線回歸方程精度更高；二是擴大研究范圍，在包括和縣、含山縣在內(nèi)的馬鞍山市全境內(nèi)進(jìn)行，使研究成果可在全市范圍內(nèi)發(fā)揮作用；三是擴展研究樹種，從國外松逐步擴展至馬鞍山市全境內(nèi)的硬闊、軟闊、楊樹、杉類、其它松類、柏類等；四是分林木不同伐倒時段、不同立地條件等進(jìn)行調(diào)查研究。

參考文獻(xiàn)：

[1]石兆明. 宜賓主要用材樹種（組）樹干重量與木材材積的關(guān)系模型研究[D].成都：四川農(nóng)業(yè)大學(xué)，2012.

[2]馬鞍山市人民政府.馬鞍山市第三次全國國土調(diào)查主要數(shù)據(jù)公報[EB/OL].（2022-03-10）[2022-03-10].https：//www.mas.gov.cn/xxgk/openness/detail/content/622ab263886688bb2d8b4570.html.

[3]李鳳日.測樹學(xué)[M].北京：中國林業(yè)出版社，2019.

[4]曾偉生.全國立木生物量方程建模方法研究[D].北京：中國林業(yè)科學(xué)研究院，2011.

[5]謝華.國外松一元立木材積表的編制[J].安徽農(nóng)業(yè)科學(xué)，2015，43（20）：206-207+235.

[6]賈乃光.數(shù)理統(tǒng)計[M].北京：中國林業(yè)出版社，1979.