基于不同水肥措施下的桉樹二元材積方程研究

黃鋒，劉文泉，孔杰君，邱權(quán)，何茜，李吉躍，蘇艷

基于不同水肥措施下的桉樹二元材積方程研究

黃鋒，劉文泉，孔杰君，邱權(quán)，何茜，李吉躍，蘇艷

(華南農(nóng)業(yè)大學(xué)林學(xué)與風(fēng)景園林學(xué)院，廣東 廣州 510520)

以不同水肥措施下的20株4年生桉樹樣木為材料，通過伐倒解析樣木收集建模資料。在干形差異分析的基礎(chǔ)上，將4個處理分為無施肥與施肥兩組處理進(jìn)行模型擬合。采用模型選優(yōu)法從6個備選模型中分別篩選出無施肥處理和施肥處理的最優(yōu)模型。分別將兩組處理最優(yōu)模型的檢驗(yàn)結(jié)果與廣東省韶關(guān)市編制的山本和藏式二元材積模型=0.000 12的模型檢驗(yàn)結(jié)果進(jìn)行比較，選出誤差最小的模型作為該試驗(yàn)地的二元材積方程。結(jié)果表明：無施肥處理的最優(yōu)模型為模型6孟憲宇式，材積方程為=?0.014 675+0.001 078?0.000 105?0.000 155+0.000 050×，施肥處理的最優(yōu)模型為模型6孟憲宇式，材積方程為=0.067 455?0.009 072+0.000 831+0.028 237?0.000 340。模型6孟憲宇式相較于韶關(guān)市編制的山本和藏式，誤差更小，精度更高，更適合于該試驗(yàn)地桉樹人工林材積和蓄積量的估算。

桉樹；單株材積；模型擬合；二元材積方程

在森林資源調(diào)查中，蓄積量調(diào)查是重要方面之一，它是反映森林原材生產(chǎn)能力的重要指標(biāo)。森林的蓄積量是由所有林分內(nèi)林木個數(shù)和林木單株材積決定的。材積是反映單株林木體積的指標(biāo)，受自身生長狀況和該地區(qū)森林生態(tài)變化以及森林經(jīng)營利用的影響而變化?？沙掷m(xù)的森林管理需要對不斷增長的蓄積量進(jìn)行估計以指導(dǎo)森林管理人員估計木材價格和確定砍伐面積，木材蓄積量可以從容易測量得到的樹木生長指標(biāo)中估計出來，最常見的方法是根據(jù)材積與胸徑和樹高等變量之間的關(guān)系使用材積方程。應(yīng)用材積方程估算單株材積降低了野外森林蓄積量調(diào)查的難度，提高了調(diào)查的效率。盡管材積方程減輕了林業(yè)工作者野外調(diào)查的負(fù)擔(dān)，但地理位置、氣候條件、環(huán)境因素、人為干擾和樹種遺傳等影響都會使樹木發(fā)生異速生長，導(dǎo)致材積方程的適用范圍變得有限，應(yīng)該針對特定條件、樹種和環(huán)境編制不同的材積方程或者把不同的變異源作為自變量包括在方程中。由于森林樹種豐富度和內(nèi)部變異性較大，通過收集涵蓋足夠數(shù)量的變異源和樹種組合的數(shù)據(jù)建立適用范圍較廣的“通用模型”使用起來不切實(shí)際，應(yīng)用性也不強(qiáng)；而建立應(yīng)用范圍較窄的相對精確的模型則可靠性更高。林場級別經(jīng)營單位的工作者們?yōu)榱擞嬎懔謭雒磕甑慕?jīng)濟(jì)效益而建立的專門適用于該單位的材積表，相較于全國級別和省市級別的材積表，這種材積表減小了誤差，提高了查表效率和精度。

桉樹()因其適應(yīng)性強(qiáng)、生長迅速、產(chǎn)量高和木材工業(yè)特性好等特點(diǎn)，被廣泛用于木材加工制造業(yè)，是華南地區(qū)主要的工業(yè)用材樹種。有關(guān)桉樹材積模型的研究已見報道，如岑巨延應(yīng)用現(xiàn)代建模思想與方法，優(yōu)化模型結(jié)構(gòu)，建立了廣西桉樹二元立木材積動態(tài)模型。馮強(qiáng)等以山本材積式為基本模型，對可變參數(shù)動態(tài)模型進(jìn)行了模型研建，編制了海南桉樹立木材積表。羅明永采用多模型選優(yōu)和逐步回歸法研制了福建省尾葉桉人工林二元立木材積表。這些材積模型的研建都是基于自然環(huán)境下生長的桉樹人工林，而目前基于水肥管理下的桉樹材積模型研究仍目前處于空白。因此，本文以廣東省增城市華南農(nóng)業(yè)大學(xué)實(shí)驗(yàn)基地的桉樹人工林為研究對象，基于收集到的代表性樣木為實(shí)驗(yàn)材料，建立二元材積方程，通過與廣東省韶關(guān)市編制的材積方程進(jìn)行對比，選出材積預(yù)測誤差更小的方程作為該試驗(yàn)地的二元材積預(yù)估方程，旨在提高該試驗(yàn)地桉樹人工林林分蓄積量的估算準(zhǔn)確度。

1 試驗(yàn)地和林分概況

試驗(yàn)地位于廣東省增城市寧西鎮(zhèn)，該鎮(zhèn)位于增城市西南部(23°05′～23°37′ N，113°32′～114°00′ E)。屬丘陵地帶，地勢北高南低，南亞熱帶季風(fēng)濕潤氣候，年均氣溫21.8℃，年均降水量為2 137 mm，多集中在4—9月。試驗(yàn)地沿著山坡陽面(南面)分成較為水平的五個梯度，面積從北到南(山頂?shù)缴降?分別為599.2、1 044.23、1 231.05、1 510.4、1 282.5 m。共5 667.38 m，約0.57 hm。于2017年4月新造林，苗木均為尾巨桉(×)DH32-29無性系組培苗，株行距為2 m × 3 m，初始造林密度為1 667株·hm。

本試驗(yàn)采用隨機(jī)區(qū)組設(shè)計，共設(shè)置4個不同的水肥處理(表1)，分別為(1)無施肥條件下無旱季補(bǔ)水處理(無水無肥處理，用CK表示)；(2)無施肥條件下有旱季補(bǔ)水處理(有水無肥處理，用W表示)；(3)有施肥條件下無旱季補(bǔ)水處理(無水有肥處理，用F表示)；(4)有施肥條件下有旱季補(bǔ)水處理(有水有肥處理，用WF表示)。每個處理設(shè)置5個重復(fù)小區(qū)，4個處理共20個小區(qū)。造林后，于每年的旱季(10月—翌年3月)對試驗(yàn)地W和WF處理進(jìn)行灌溉補(bǔ)水，灌溉采用滴灌方式，灌溉頻率為每周2次，每次連續(xù)滴灌4 h。肥料采用桉樹專用基肥與追肥，施肥量均與生產(chǎn)桉樹用材林的用量一致。基肥和兩次追肥分別于2017年3月、2017年7月和2018年7月施于F、WF所有樹木中。

表1 試驗(yàn)地水肥處理設(shè)計

2 研究方法

2.1 林分調(diào)查與樣木選擇

于2020年7月12日對該試驗(yàn)地所有小區(qū)的4年生桉樹進(jìn)行每木檢尺，調(diào)查指標(biāo)包括樹高和胸徑。根據(jù)每個處理小區(qū)每木檢尺結(jié)果，計算得出每個處理的平均胸徑和平均樹高，選取各處理內(nèi)生長健康，樹形良好(無明顯受損，樹干通直不彎曲不傾斜)且胸徑和樹高最接近平均值的立木作為樣木，每個處理共選5株樣木，4個處理共20株。

2.2 樣木伐倒與材積測定

樣木選好后，下一步進(jìn)行砍伐工作。伐倒前，應(yīng)使用油漆筆準(zhǔn)確標(biāo)注樹干的根頸位置和胸徑位置，并使用噴漆沿著樹干基部至胸徑處標(biāo)記南北方向，用油鋸或斧頭在根頸處進(jìn)行伐倒后，用皮尺測出實(shí)際樹干高，用胸徑尺測量實(shí)際胸徑。單株材積采用中央斷面區(qū)分求積法求得，以2 m為區(qū)分段。具體步驟如下：以2 m為段長(第一段長度為2.6 m，其余2 m)，用噴漆在胸徑處(1.3 m)，3.6 m，5.6 m，7.6 m……，和樹梢(不足1個區(qū)分段)底部處標(biāo)記出需要截取樣木圓盤的位置和南北方向，接著對解析圓盤進(jìn)行截取，圓盤厚度一般為3～5 cm，截取時應(yīng)盡量使油鋸與干軸保持垂直，不可偏斜，同時盡量使斷面平滑。圓盤向地的一面要恰好在各分段的標(biāo)定位置上，以該面為工作面，用來量測直徑，另一面為非工作面。每個圓盤鋸下后，應(yīng)立即在當(dāng)天完成工作面直徑測量，避免標(biāo)本圓盤縮水產(chǎn)生誤差。直徑測量用游標(biāo)卡尺分別測量每個圓盤南北、東西兩個方向的去皮直徑，兩者求平均值作為該圓盤的直徑。材積計算是利用測得的圓盤直徑算出各段的體積，各段體積之和即單株材積，每段的體積為圓柱求積公式，即段長×圓盤面積，梢頂為圓錐求積公式，即1/3梢長×梢底面積，材積具體計算公式如下：

式中：為單株材積(m)，為各區(qū)分段斷面積(m)，為各區(qū)分段長度(m)，g為梢頭斷面積(m)，l為梢頭長度(m)。

2.3 胸高形數(shù)計算

胸高形數(shù)是林木單株材積與胸高斷面積乘以樹高的比值，它是反映林木干形的常用指標(biāo)。其計算公式如下：

式中：為立木的單株材積，為胸高形數(shù)，為胸高直徑，為樹高。

2.4 數(shù)據(jù)處理與分析

使用WPS excel軟件對實(shí)測的樹高、胸徑進(jìn)行分類整理，用韶關(guān)市編制的材積方程得出單株材積理論值，用中央斷面區(qū)分求積法公式計算出各解析木單株材積實(shí)測值，并用公式計算出其胸高形數(shù)。使用SPSS 22.0軟件對各解析木胸高形數(shù)進(jìn)行單因素方差分析和鄧肯式多重比較，并對解析木的樹高、胸徑和材積進(jìn)行非線性回歸分析。

3 結(jié)果與分析

3.1 干形差異分析

若建模樣本(表2)之間干形差異較大，將樣本總體擬合會降低最終方程的預(yù)估精度。因此，在對4個處理的桉樹胸高形數(shù)進(jìn)行顯著性差異分析的基礎(chǔ)上(表3)，再決定是否應(yīng)該分處理進(jìn)行建模。CK處理與W處理之間胸高形數(shù)無顯著差異，F(xiàn)處理和WF處理之間胸高形數(shù)無顯著差異，CK與F和WF處理之間胸高形數(shù)無顯著差異，W與F和WF處理之間胸高形數(shù)無顯著差異。施肥處理與不施肥處理之間干形具有顯著差異，說明施肥影響了桉樹的干形生長。因此，為了提高材積方程的預(yù)估精度，需要將有干形差異的桉樹分組建模，即分為施肥與無施肥兩組處理進(jìn)行模型擬合。

3.2 模型回歸分析

為了保證研究的準(zhǔn)確性以及適應(yīng)性，從我國林業(yè)行業(yè)標(biāo)準(zhǔn)《二元立木材積表編制技術(shù)規(guī)程》(LY/T 2012—2013)推薦的12種二元立木材積模型中選擇擬合桉樹二元材積式常用的6個回歸模型(表4)作為備選模型。以胸徑和樹高的實(shí)測值作為模型的自變量，以通過中央斷面分區(qū)求積法得到的材積作為模型的因變量，利用SPSS 22.0軟件的“非線性回歸”功能進(jìn)行模型的回歸分析，將參數(shù)a、b、c、d、e初始化為0，選擇Levenberg-Marquardt估計方法對各樣本數(shù)據(jù)進(jìn)行回歸分析，最終求得a、b、c、d、e的值。

表2 建模樣本資料

決定系數(shù)越大越接近1，說明胸徑和樹高與材積之間相關(guān)性更高，擬合效果更好。無施肥處理中6個模型的決定系數(shù)均在97%以上(表5)，施肥處理的6個模型的決定系數(shù)均在90%以上(表6)，說明這6個模型均能較好地擬合無施肥與施肥處理的材積與樹高和胸徑的關(guān)系。在無施肥處理中，模型6的為0.978，稍大于其他5個模型，施肥處理中，模型6的為0.962，明顯大于其它5個模型。因此，根據(jù)擬合結(jié)果得出，無施肥處理和施肥處理的最優(yōu)模型都是模型6(孟憲宇式)。

表3 處理間胸高形數(shù)差異分析

注：表中不同數(shù)字后小寫字母表示<0.05。

表4 備選模型

表5 無施肥處理擬合結(jié)果

表6 施肥處理擬合結(jié)果

3.3 模型精度檢驗(yàn)

利用模型6(孟憲宇式)和韶關(guān)式(山本和藏式)得出的材積理論值與實(shí)際值計算出單株材積的誤差絕對值，以平均誤差絕對值作為模型精度檢驗(yàn)的指標(biāo)(表7～8)。無施肥處理下，模型6(孟憲宇式)的平均誤差絕對值為1.87%，韶關(guān)市編制的二元材積式(山本和藏式)的平均誤差絕對值為3.67%；施肥處理下，模型6(孟憲宇式)的平均誤差絕對值為1.67%；韶關(guān)市編制的二元材積式(山本和藏式)的平均誤差絕對值為12.03%。無論是無施肥處理還是施肥處理，模型6(孟憲宇式)的平均誤差絕對值都小于韶關(guān)市編制的材積模型。

表7 無施肥處理兩個模型精度檢驗(yàn)結(jié)果

表8 施肥處理兩個模型精度檢驗(yàn)結(jié)果

4 結(jié)論與討論

本研究在干形差異分析的基礎(chǔ)上對4個水肥措施下的桉樹林分進(jìn)行分組建模，最終確定模型6(孟憲宇式)為該試驗(yàn)地施肥處理和不施肥處理的二元材積模型。模型6擬合回歸得到桉樹無施肥處理和施肥處理的兩個二元立木材積方程的決定系數(shù)均在96%以上，擬合效果良好，方程能很好的擬合樹高和胸徑的關(guān)系。材積方程分別為無施肥處理：=?0.014 675+0.001 078?0.000 105?0.000 155+0.000 050×，施肥處理：=0.067 455?0.009 072+0.000 831+0.028 237?0.000 340。

從兩個模型的精度檢驗(yàn)結(jié)果可知，相較于廣東省韶關(guān)市編制的山本和藏式材積方程，本研究篩選出的最優(yōu)模型孟憲宇式提高了對該試驗(yàn)地的桉樹材積的預(yù)估精度。韶關(guān)式編制的材積方程是基于自然條件生長的桉樹而編制，未考慮施肥對桉樹干形生長的影響，因此在施肥處理中，用韶關(guān)市編制的材積式的預(yù)測值普遍偏小，誤差絕對值均超過了10%。這也應(yīng)證了土壤養(yǎng)分特征變化會對樹木干形產(chǎn)生顯著影響，應(yīng)用相同的材積式估算單株材積可能會導(dǎo)致誤差較大。而通過模型6孟憲宇式計算出的兩組處理的平均誤差絕對值均不超過3%，這表明用孟憲宇式對單株材積的預(yù)估準(zhǔn)確度可達(dá)97%以上，這說明分組建?？梢詼p少組內(nèi)異速生長給模型帶來的影響，是建立更準(zhǔn)確模型的可行方法。

模型6(孟憲宇式)具有良好的全面切合性能，可以應(yīng)用于該試驗(yàn)地桉樹人工林單株材積和林分蓄積的估算，但由于本研究樣本量較少，未對最優(yōu)模型(孟憲宇式)進(jìn)行適用性檢驗(yàn)，因此，本研究結(jié)果僅對該試驗(yàn)地的單株材積和林分蓄積的估計有效，能否適用于廣東省其他地區(qū)的桉樹林，有待進(jìn)一步考證。

[1] MAUYA E W, MUGASHA W A, ZAHABU E, et al. Models for estimation of tree volume in the miombo woodlands of Tanzania[J]. Southern forests, 2014, 76(4): 209-219.

[2] 李奇,朱建華,馮源,等.中國森林喬木林碳儲量及其固碳潛力預(yù)測[J].氣候變化研究進(jìn)展,2018,14(3):287-294.

[3] AKINDELE S O, LEMAY V M. Development of tree volume equations for common timber species in the tropical rain forest area of Nigeria[J]. Forest Ecology and Management. 2006, 226(1/3): 41-48.

[4] 何誠,馮仲科,劉柯珍,等.云杉二元材積與蓄積量精度分析模型研究[J].林業(yè)調(diào)查規(guī)劃,2010,35(4):1-6.

[5] 張偉,岑巨延,秦旭東,等.廣西杉木胸徑和樹高二元立木材積模型研究[J].廣西林業(yè)科學(xué),2020,49(3):397-402.

[6] FELDPAUSCH T R, BANIN L, PHILLIPS O L, et al. Height-diameter allometry of tropical forest trees[J]. Biogeosciences Discussions, 2010,7(5):7727-7793.

[7] 曹忠.立木材積無損精測與建模方法研究[D].北京:北京林業(yè)大學(xué),2015.

[8] 國家林業(yè)和草原局.中國森林資源報告(2014—2018)[M].北京:中國林業(yè)出版社,2019.

[9] 岑巨延.廣西桉樹人工林二元立木材積動態(tài)模型研究[J]. 華南農(nóng)業(yè)大學(xué)學(xué)報,2007,28(1):91-95.

[10] 馮強(qiáng),方良.海南桉樹立木材積模型研建[J].中南林業(yè)調(diào)查規(guī)劃,2014,33(3):47-51.

[11] 羅明永.尾葉桉人工林二元材積表的研制[J].福建林業(yè)科技,2009,36(4):170-175.

[12] 王海華.增城區(qū)鄉(xiāng)土樹種人工林生長表現(xiàn)及生態(tài)效益評價[D].廣州:華南農(nóng)業(yè)大學(xué),2016.

[13] 唐偉東.紅松材積與生長性狀相關(guān)分析及二元回歸方程的擬合[J].林業(yè)勘查設(shè)計,2017(1):42-43.

[14] 吳志華,尚秀華,王睿,等.基于干曲線方程的桉樹立木單株材積測算研究[J].桉樹科技,2018,35(3):1-8.

[15] 全國森林資源標(biāo)準(zhǔn)化技術(shù)委員會.二元立木材積表編制技術(shù)規(guī)程:LY/T 2102-2013[S].北京:中國標(biāo)準(zhǔn)出版社,2013.

[16] 劉治平,曹宗英,張旭華,等.甘肅省人工云杉林二元立木材積方程的編制研究[J].林業(yè)資源管理,2020(1):172-176.

[17] PHIRI D, PHIRI E, KASUBIKA R, et al. The implication of using a fixed form factor in areas under different rainfall and soil conditions forin Zambia[J]. Southern forests,2016,78(1):35-39.

[18] MUGASHA W A, BOLLANDSAS O M, EID T, et al. Relationships between diameter and height of trees in natural tropical forest in Tanzania[J]. Southern Forests: a Journal of Forest Science,2013,75(4):221-237.

Study on Binary Volume Equation ofBased on Different Water and Fertilizer Measures

HUANG Feng, LIU Wenquan, KONG Jiejun, QIU Quan, HE Qian, LI Jiyue, SU Yan

()

Data were obtained from 20 felled 4-year-oldsample trees, that had been grown under different water and fertilizer regimes. Based on stem shape analyses, the four treatments were divided into two groups: no fertilization and fertilization for model fitting. The selection methods used to select the optimal model for no fertilization treatment and the fertilization treatment involved examination of six alternative models. The results of the two groups of optimal models were compared with the model test results of Yamamoto and Tibetan binary volume model=0.000 12compiled by Shaoguan City, Guangdong Province, respectively, and the model with the least error was selected as the binary volume equation of the test site. The results showed that the optimal model of no fertilization treatment was model 6 Meng Xianyu comprising the following:=?0.014 675+0.001 078?0.000 105?0.000 155+0.000 050×. For the fertilization treatment, the optimal model was model 6 Meng Xianyu comprising the following:=0.067 455?0.009 072+0.000 831+0.028 237?0.000 340. Model 6 Meng Xianyu volume model has smaller error and higher precision than Yamamoto and Tibetan volume model compiled by Shaoguan City, whilst also being more suitable for the estimation of individual tree volume and also stand volumes ofplantation in this experimental site.

; individual tree volume; model fitting; binary volume equation

S758.2

A

10.13987/j.cnki.askj.2021.03.003

國家重點(diǎn)研發(fā)計劃“桉樹、馬尾松人工林生產(chǎn)力形成的生理”（2016YFD0600201;2016YFD060020102）；國家自然科學(xué)基金“桉樹人工林生產(chǎn)力對旱季水分輸入的變化”(31800527)

黃鋒(1997— )，男，在讀碩士，主要從事林木生長模型研究，E-mail:249211084@qq.com

李吉躍(1959— )，男，博士，教授，主要從事森林培育研究，E-mail:564857527@qq.com