樹(shù)木生長(zhǎng)與降水量擬合方程平衡精度的研究

李斐，張?chǎng)危瑥埦福瑥埖氯诹栾w，張鵬

(1.山東省招遠(yuǎn)市自然資源和規(guī)劃局，山東 煙臺(tái) 265701；2.龍口市林業(yè)技術(shù)推廣站，山東 煙臺(tái) 265701；3.山東省林業(yè)監(jiān)測(cè)規(guī)劃院，山東 濟(jì)南 250014)

2010年以來(lái)，各種經(jīng)驗(yàn)方程在林業(yè)科學(xué)實(shí)踐中得到了廣泛應(yīng)用，鞏延蘋(píng)，張德全，矯興杰等人先后利用經(jīng)驗(yàn)方程對(duì)樹(shù)木與降水量的關(guān)系、碳儲(chǔ)量動(dòng)態(tài)、樹(shù)木成熟齡、光照對(duì)樹(shù)木生長(zhǎng)影響及生長(zhǎng)因子剔除等進(jìn)行了大量研究，取得了較好的效果[1-7]，但是至于試驗(yàn)精度的問(wèn)題，一直是大家十分困惑的問(wèn)題，筆者已對(duì)時(shí)間因子的試驗(yàn)精度進(jìn)行了研究，下面再對(duì)降水量試驗(yàn)精度進(jìn)行研究。

1 數(shù)據(jù)來(lái)源

數(shù)據(jù)資料來(lái)源于在 2014年做了牟平區(qū)6棵赤松以1 a為齡階的解析木材料，關(guān)于解析木的做法，筆者在以前的文獻(xiàn)中已作過(guò)介紹，在這里不再累述。降水量數(shù)據(jù)由牟平區(qū)氣象局提供。

2 研究方法

張德全等人曾經(jīng)用經(jīng)驗(yàn)方程對(duì)樹(shù)木生長(zhǎng)過(guò)程做過(guò)多次研究，用經(jīng)驗(yàn)方程y=exp(a+bx)(其中y為樹(shù)木的各種生長(zhǎng)量指標(biāo)，如樹(shù)徑、樹(shù)高、材積等，為因變量，x為影響因變量變化的各項(xiàng)因子值，即時(shí)間、光照、氣溫、降水等，a，b為待求系數(shù)，由于本函數(shù)為成長(zhǎng)函數(shù)，樹(shù)木一般不會(huì)出現(xiàn)負(fù)增長(zhǎng)，故b一般為正值)研究樹(shù)木生長(zhǎng)快漫，及各項(xiàng)因子及其交互作用的影響程度，用經(jīng)驗(yàn)方程y=exp(a-b/x)(方程中的符號(hào)含義同上一經(jīng)驗(yàn)方程)來(lái)研究因變量的平均變化速度和即時(shí)變化速度，即x=b時(shí)，因變量平均變化速度最快，當(dāng)x=b/2時(shí)即時(shí)變化速度最快(通過(guò)對(duì)方程的一階求導(dǎo)和二階求導(dǎo)可得)，用以追求效益最大化。先用經(jīng)驗(yàn)方程y1=exp(a1+b1p)和y2=exp(a2-b2/p)(為了加以區(qū)別，在函數(shù)和待求系數(shù)加了下標(biāo))建立樹(shù)木生長(zhǎng)量與年降水量的回歸關(guān)系，求得待求系數(shù)a1、b1、a2、b2，在這里y1、y2為以降水量p為自變量、樹(shù)木實(shí)際生長(zhǎng)量(實(shí)測(cè)值)為因變量而得到的經(jīng)驗(yàn)方程擬合值。

如果令y3=exp(a3+b3p)和y4=exp(a4-b4/p)，y3、y4為對(duì)y1、y2而言， 將上述的樹(shù)木生長(zhǎng)量(實(shí)測(cè)值)替換為y1^m×y2^n，其中m、n為小于1的純小數(shù)，滿足m+n=1的條件，并且使得方程y3、y4的方程擬合精度完全相等，m、n用電子表格計(jì)算采用逐步接近的辦法解決求算問(wèn)題，筆者將求算精度定為10億分之一，即1E-9。

3 研究結(jié)果與分析

表1的徑階起始年齡為35 a,對(duì)于6齡階來(lái)說(shuō)，如果起始年齡增大，其平衡精度反而降低(時(shí)間方程增大)，其余均表現(xiàn)為起始年齡越大。在齡階數(shù)相同的前提下，其試驗(yàn)精度則越高。對(duì)于42齡階來(lái)說(shuō)，則表現(xiàn)出了在部分起始年齡不同，但是平衡精度卻相同，只有極個(gè)別不一致現(xiàn)象,表現(xiàn)異常的均為65年生赤松(前面得到的6齡階異常也是如此)，由于樹(shù)木高大，采取2 m區(qū)段(其他為1 m區(qū)段)，但是對(duì)于這棵65年生解析木來(lái)說(shuō)，其胸徑、樹(shù)高和材積項(xiàng)目則表現(xiàn)正常。

表1 各齡階平衡精度表 %

從表1中可以看出，平衡精度從3齡階開(kāi)始，隨著齡階的上升，其平衡精度表現(xiàn)出逐步增加的態(tài)勢(shì)，但19齡階達(dá)到峰值99.8670%，20齡階又略有下降，至42齡階降到最低，至于42齡階以上的表現(xiàn)，需要另外收集更多年份的氣象資料進(jìn)行研究了。從表中可以看出，對(duì)于同一齡階來(lái)說(shuō)，降水量的平衡精度則比時(shí)間平衡精度要高，因此可以說(shuō)降水量的表現(xiàn)更為敏感，因此降水量的年際波動(dòng)對(duì)樹(shù)木生長(zhǎng)則有較為深刻的影響，由于時(shí)間的運(yùn)動(dòng)是均一的，無(wú)差別的，它表現(xiàn)只能是累加效應(yīng)。

如果將樹(shù)木生長(zhǎng)量替換為y1^m/y2^n，建立y4=exp(a4+b4p)和y5=exp(a5-b5/p)的擬合經(jīng)驗(yàn)方程，則y4和y5的方程擬合精度非常接近，理論上是一致的，由于計(jì)算是用試算的辦法產(chǎn)生的，涉及小數(shù)進(jìn)位的問(wèn)題，故出現(xiàn)了細(xì)微偏差，但試驗(yàn)精度大大降低，相關(guān)系數(shù)降低40%以上，這就從另一方面表明y1^m與y2^n具有同等重要的作用，即達(dá)到平衡狀態(tài)，那么將系數(shù)指數(shù)m、n進(jìn)行互換處理，將樹(shù)木生長(zhǎng)量替換為y1^n×y2^m則為樹(shù)木生長(zhǎng)的理想實(shí)際狀態(tài)，再次建立樹(shù)木生長(zhǎng)擬合經(jīng)驗(yàn)方程y6=exp(a6+b6p)和y7=exp(a7-b7/p)，即為樹(shù)木生長(zhǎng)的理想實(shí)際降水生長(zhǎng)方程和降水生長(zhǎng)阻力方程。

4 結(jié)論與討論

對(duì)于經(jīng)驗(yàn)方程y=exp(a+bp)與y=exp(a-b/p)來(lái)說(shuō)，對(duì)于完整齡階，19個(gè)齡階(即樣本數(shù))即達(dá)到最高精度要求，并非樣本數(shù)(齡階數(shù)量)越多，精度越高，這樣利用數(shù)據(jù)的累加效應(yīng)，現(xiàn)代計(jì)算機(jī)技術(shù)的應(yīng)用，使得任意年齡作為齡階年齡成為可能，可對(duì)解析木的總年齡先行進(jìn)行預(yù)判，通過(guò)走訪調(diào)查樹(shù)木年齡，可以大大減少樹(shù)木年輪判讀工作量，起到事半功倍的作用。

經(jīng)過(guò)研究發(fā)現(xiàn)，對(duì)于經(jīng)驗(yàn)方程y=exp(a+bx)的演算，當(dāng)x的所有數(shù)據(jù)同時(shí)加上或減去同一數(shù)值，試驗(yàn)結(jié)果只是a值變化，b值和用來(lái)檢驗(yàn)的t值、F值、R值及試驗(yàn)精度保持不變，因此通過(guò)這一方程，來(lái)實(shí)現(xiàn)以往沒(méi)有實(shí)際記載數(shù)據(jù)的重建成為可能，而保持現(xiàn)有試驗(yàn)的精度不變。試驗(yàn)表明，只要是樣本數(shù)目相同，解析木區(qū)段長(zhǎng)度相同，齡階間隔時(shí)間長(zhǎng)度一致，其平衡精度是完全一致的。但是對(duì)于2 m區(qū)段來(lái)說(shuō)，只有胸徑和樹(shù)高、材積表現(xiàn)和1 m區(qū)段一樣，其余均表現(xiàn)異常。研究表明，起始年齡越大，其平衡精度越高(2 m區(qū)段同1 m區(qū)段相比也表現(xiàn)異常)。至于不同區(qū)段表現(xiàn)不一致，胸徑和樹(shù)高、材積三個(gè)項(xiàng)目表現(xiàn)又一致，究其是什么原因，需待以后的研究中加以解決。研究表明相同條件下，降水量的平衡精度要高于時(shí)間要素的平衡精度，說(shuō)明降水量的年際變化，對(duì)樹(shù)木的影響更為敏感。