黑龍江省樟子松人工林一元立木材積表驗證分析

李笑峰

(黑龍江省林業監測規劃院，黑龍江 哈爾濱 150080)

黑龍江省市縣林區一直使用《黑龍江省立木材積表》(黑龍江省營林局編，1981年)計量人工樟子松林的材積，該表編制至今已使用近40 年，近年來在使用過程中，發現該材積表存在偏差，甚至出現出材量大于單株木材積的現象。為了進一步分析人工樟子松材積表的估計誤差，利用近兩年收集的72 株人工樟子松樣木立木材積數據，對《黑龍江省立木材積表》中的5個人工樟子松一元立木材積表和全省樟子松一元立木材積表表進行檢驗，并對其進行差異分析，為客觀評價人工樟子松立木材積表的現狀提供科學依據，也為推動我省林木材積表的更新和修訂工作提供決策參考。

1 驗表樣木來源

驗表樣木來源：一是2016年和2017年專業調查森林資源調查標準地中的計算木和解析木；二是全國主要樹種生物量模型編制工作中所采集的樣木。樣木采集于樺南、樺川、湯原、林口、東寧、五大連池、蘿北、孫吳、嫩江、嘉蔭、海倫、賓縣等市縣區。

本次樣本采集以5.0 cm為起測徑階。樣木伐倒后，樹高小于10 m以1 m為一個區分段長，樹高大于等于10 m以2 m為一個區分段長，按中央斷面積區分求積法用圍尺量測各區分段帶皮直徑及胸高部位的帶皮直徑(精確到 0.1 cm)，用皮尺測量樹高(精確到 0.1 m)。將采集樣木的胸徑、樹高及各區分段的中央直徑，在計算機上建立數據庫，并按中央斷面積區分求積式計算各樣木的材積。將各樣木的胸徑(D)、樹高(H)及材積(V)建立數據庫，剔除異常點后，作為本次檢驗材積表的基礎數據，以檢驗材積表的實際精度和適用性。樣本基本情況見表1。

表1 樟子松樣本數據基本情況

2 現行人工樟子松材積表的檢驗

2.1 材積表的檢驗方法

設采集樣木的實測材積為yi，材積表材積為xi分析以下誤差：

(1)平均誤差(ME)：

(2)平均相對誤差(M%E)：

(3)系統誤差(SE%)：

根據《林業專業調查主要技術規定》，立木材積表的系統誤差一般不得超過 ±3% (中華人民共和國林業部，1990)。根據檢驗樣本的實際觀測值yi和材積表材積xi，計算平均相對誤差(M%E)和系統誤差(SE%)。如果M%E和SE%均不超過 ±3%，則認為原來的模型仍然是適用的；否則，就是不適用的。

2.2 材積表的預估精度及置信橢圓F檢驗

預估精度計算及置信橢圓F檢驗方法，是實測材積(y)與材積表材積(x)之間建立一元線性回歸方程：

yi=α+βxi+εi

(1)

回歸標準差(Sy.x):

(2)

(3)

誤差限(Δ):

(4)

材積表的預估精度(P%):

(5)

置信橢圓F檢驗則是在置信水平取為1-α時，對回歸模型系數α和β構造聯合置信區域。這個區域由下式給出：

(6)

從而一次完成對回歸系數α和β的顯著性檢驗。

如果材積表材積(x)與實測材積(y)完全一致，則α=0，β=1。但實際上往往不一致。因此，需要檢驗由樣本估計的a和b的值與它們的真值之間有無顯著差異。將α=0，β=1代入(6)式可得：

(7)

式中：n為檢驗樣木的株數，yi為實測材積，xi為材積表材積，a、b為回歸直線的參數。

(7)式服從自由度f1=2，f2=n-2的F分布。這個置信區域的邊界是以(a，b)為中心的橢圓。根據按(7)式計算的 F 值與α顯著水平下的理論 Fα值，來判斷材積表的適用性。當F>F0. 05時，線性回歸方程y=a+bx其參數(a，b)與 (0，1)有顯著差異，說明材積表已經存在系統誤差，已不適用；當F≤F0. 05 時，參數(a，b)與 (0，1)無顯著差異，說明不存在系統誤差，原來的立木材積表仍然是適用。

3 材積表的檢驗結果與分析

《黑龍江省立木材積表》(黑龍江省營林局編，1981年)中樟子松材積表共有5個，它們分別是：

(1)《牡丹江地區樟子松人工林圍徑一元立木材積表》

(2)《松花江地區樟子松人工林圍徑一元立木材積表》

(3)《合江地區樟子松人工林圍徑一元立木材積表》

(4)《黑河地區樟子松人工林圍徑一元立木材積表》

(5)《嫩江平原地區樟子松人工林圍徑一元立木材積表》

據調查發現，當前我省現存的樟子松人工林，特別是林分平均胸徑超過20 cm的林分，其林相與樟子松天然比較接近，立木干型相似，因此，利用本次采集的樣本數據，對樟子松天然林的一元材積表進行誤差分析和適用性檢驗。

利用采集的檢驗樣本數據(見表1)，對現行的5個樟子松一元立木材積表進行了材積誤差分析、預估精度及置信橢圓F檢驗，檢驗結果見表2、表3。

表2 樟子松材積表誤差分析表

表3 樟子松材積表檢驗結果 =0.05

表3 樟子松材積表檢驗結果 =0.05

材積表樣本數量 N回歸標準差Sy.x回歸標準誤 S`y誤差限D預估精度p%FFa牡丹江930.065 50.006 80.013 394.531 037.483.98松花江930.063 40.006 60.012 994.711 114.693.98合江860.007 90.000 90.001 799.161 046.863.99黑河930.062 20.006 50.012 794.811 160.293.98松嫩平原630.034 90.004 40.008 692.27362.464.06

根據分析和檢驗的結果，黑龍江省樟子松人工林目前所使用的5個一元立木材積表，其相對誤差(M%E)和系統誤差(SE%)分別大于±36.6%和±43.0%，遠遠超過了±3%的要求，材積表估計的材積偏低很多。

所有材積表的置信橢圓 F 值均大于臨界值F0.05，推翻原假設，即推翻實測材積(y)與材積表材積(x)完全一致的假設，二者之間差異顯著。

通過檢驗得出《黑龍江省立木材積表》對人工樟子松材積的估計不具有適用性。因此，應盡快編制符合我省樟子松人工林現實情況的一元或二元材積表。

4 誤差產生的原因

通過精度計算和實用性檢驗，我省現行的5個樟子松人工林一元立木材積表，精度均超過允許范圍，不具有實用性。

4.1 林分結構以及單株木材積發生變化

《黑龍江省立木材積表》于1981年編制完成并頒發使用，其編表樣本多采集于1978年，當時我省樟子松人工林林分年齡一般不超過30年，且多為幼中齡林。根據歷史資料，1986年黑龍江省市縣林區人工林幼(a≤20)、中(a=20～30)、近成過(a>30)以上面積分別占總面積的45.8%、10.6%和0.3%。當時樟子松人工林的林分特點是初植密度大多為1.5 m×1.5 m或1 m×1.5 m，成林后幼齡林保留株數多3 500株以上，中齡林多在3 000株左右，林分密度大，林分年齡小，其林分平均胸徑多在12 cm以下，林分內以小徑木為主，中徑級以上林木不多。而現今樟子松人工林林分與80年代林分相比發生了較大變化，首先是林分齡組占比從以幼齡林為占比最大轉變為以中齡林為主，近熟林以上林分大幅增加，其二是林分密度降低、林分質量提高， 2000年以后營造的人工樟子松林，初植密度多為2 m×1.5 m，成林時林分年齡達到20年之間，保存株數多3 000株左右，中齡林經撫育后公頃株數多在1 500株左右，近熟林以上林分經過多次的生產經營和撫育以上株數多為800株左右。根據過往研究成果，林分密度變化可促使林木的削度和胸徑與樹高的線性關系發生變化，特別是經過撫育的林分與未經過撫育的林分比較產生的變化更大。而在胸徑相同的條件下，立木材積取決于立木的樹干削度和樹高，因此現今林分與20世紀80年代前林分比較，其林木材積發生變化是必然的。

4.2 歷史條件所限

現行樟子松人工林一元立木材積表，除了嫩江平原外都是以按當時的行政范圍為單位編制的，這種按人為區域劃分編制材積表的方法，本身就存在弊端，因為自然因素對林木的生長是有影響的，而行政區劃一般不能影響林木的生長。

由于編表時樟子松人工林幼齡林占絕對優勢，林木徑級偏小，材積表徑級范圍小，5個材積表的徑級范圍分別為3～38、3～34、3～27、3～32和3～20。現今林分超出材積表徑級的林木不在少數，甚至有超出徑級范圍20 cm以上的林木，超出范圍的林木雖然可采用插入法查其材積，但是所查材積一般都沒有經過檢驗，精度如何更是無從談起。

4.3 材積測量工具不同

由《黑龍江省立木材積表》使用說明：“三、圍尺立木材積表的編制程序(1)用圍尺徑與輪尺徑原始材料求出經驗公式并列出各區域各樹種的圍尺、輪尺換算表，查出所對應的輪尺直徑。(2)把輪尺直徑代入《黑龍江省立木材積表》中相應的一元或二元材積式即得出圍尺材積。”可知， 黑龍江省立木材積表編表樣本原始數據，為輪尺數據，即其胸徑和材積均采用輪尺測量得出，其數學模型是以輪尺胸徑與輪尺材積的相關性為基礎建立的。而《黑龍江省立木材積表》在編表時僅考慮了圍尺直徑與輪尺胸徑之間的測量誤差，如利用經驗公式D輪= -0.256 357 203 5+1.002 957 973×D圍，將樟子松人工林圍尺直徑轉為輪尺直徑，轉換后圍尺直徑比輪尺胸徑降低0.2 cm，即15.2 cm 的圍尺胸徑等于15.0 cm的輪尺胸徑，然后代入輪尺測量數據構建的數學模型，得出相對應的材積。這種只考慮圍尺直徑與輪尺直徑誤差，而不考慮圍尺材積與輪尺材積差異，直接套用模型的編表方法值得探討(研究點為圍尺胸徑d對應的數學模型的值與其圍尺材積值的差異)。此外，分別按胸徑為15.0 cm和15.2 cm，查樟子松人工林立木材積表二者材積相差3.3%，若徑級范圍由6～32 cm，直徑相差0.2 cm，查材積表其材積差由胸徑從小到大遞減，材積差幅度為1%～8%。

4.4 材積表誤差產生的影響

森林蓄積的多少是評價森林經濟和生態指標的重要因素，也是判定森林質量的重要因子。林木材積表是為計算森林蓄積量服務的，它是林業重要的基礎計量數表。如果材積表誤差過大，就會導致森林蓄積量估計產生偏差，無法獲得真實的森林資源現狀及森林資源的消長動態情況，導致森林資源調查、森林資源監測、森林資源評估、森林撫育、采伐設計、征占用林地補償等成果的精度降低和出現系統性偏差，也會對森林經營方案編制、采伐限額制定產生重大影響。甚至能夠影響到森林資源管理部門和森林資源經營部門的決策，導致林業宏觀政策如林業發展規劃、產業政策、經營方針制定出現偏差。

5 結論和建議

在森林資源管理和森林資源經營中，材積表具有不可替代的作用。目前我省市縣林區除柞樹、山楊、白樺、人工落葉松以及人工楊類以外，一直使用1981年省營林局頒布的《黑龍江省立木材積表》，使用時間已近40年。這期間我省市縣林區森林的林分質量發生了變化，特別是人工林的林相變化更大。現行材積對林分發生了變化的林分蓄積量估計的精度和適用性如何應引起高度重視。因此，應盡快對《黑龍江省立木材積表》進行檢驗或者修訂。此外，應盡快重新編制用材林中主要樹種特別是人工林樹種的一元立木材積表和二元立木材積表。

由于材積表從編制至頒布周期較長，此段時間作為過渡期，要采用科學的方法，為誤差較大無法使用的材積表確定一個替代材積表。