吳勝義，王 飛，閔志強(qiáng)

（國家林業(yè)和草原局西北調(diào)查規(guī)劃設(shè)計院 旱區(qū)生態(tài)水文與災(zāi)害防治國家林業(yè)和草原局重點實驗室，陜西西安 710048）

森林資源規(guī)劃設(shè)計調(diào)查（簡稱二類調(diào)查）是以國有林業(yè)局（場）、自然保護(hù)區(qū)和森林公園等森林經(jīng)營單位或縣級行政區(qū)域為調(diào)查單位，為滿足森林經(jīng)營方案、總體設(shè)計、林業(yè)區(qū)劃與規(guī)劃設(shè)計需要而進(jìn)行的森林資源調(diào)查[1]。二類調(diào)查是全面掌握縣域森林資源狀況的重要手段，林分蓄積量是森林資源監(jiān)測最基礎(chǔ)的信息之一[2]，其精度直接影響二類調(diào)查的成效。

森林調(diào)查是獲取森林資源信息并進(jìn)行分析與決策的一項基礎(chǔ)性工作。為了獲得森林資源數(shù)據(jù)，通常使用的方法為全林調(diào)查法，該方法費時費力，且成本高[3]。20世紀(jì)70年代以來，林業(yè)調(diào)查的主管部門根據(jù)我國森林資源調(diào)查的經(jīng)驗與實踐，提出了抽樣控制總體法，將小班調(diào)查法與抽樣調(diào)查法相結(jié)合，小班區(qū)劃調(diào)查可將資源數(shù)據(jù)落實到小班（山頭地塊），系統(tǒng)抽樣控制總體蓄積精度和小班調(diào)查蓄積，既有抽樣精度作保證，又發(fā)揮了各自的長處。該法在二類調(diào)查中得到了廣泛的應(yīng)用[4-6]。林分蓄積量是評價森林?jǐn)?shù)量、質(zhì)量及經(jīng)營效果的重要指標(biāo)，它的準(zhǔn)確估測有利于及時把握森林資源動態(tài)，是制定森林經(jīng)營計劃與采伐的重要依據(jù)[7-8]。本研究在西昌市森林資源現(xiàn)狀的基礎(chǔ)上，對西昌市總體蓄積抽樣設(shè)計方式進(jìn)行分析與比較，為總體蓄積控制和抽樣設(shè)計提供技術(shù)指導(dǎo)。

1 研究區(qū)概況

1.1 研究范圍

西昌市實測抽樣總體抽樣范圍為西昌市行政范圍內(nèi)的森林、林木和林地，包括西昌市管轄的37個鄉(xiāng)（鎮(zhèn)）和7 個國有林場（包括飛播站和森林經(jīng)營所）。

1.2 森林資源狀況

西昌市二類調(diào)查經(jīng)營范圍內(nèi)，土地總面積為265 570.73 hm2，有、疏林地面積為131 633.08 hm2。有、疏林地中，可抽樣對象面積為128 690.69 hm2，占全市有、疏林地面積的97.8%，占全市森林經(jīng)營總面積的48.5%[9]。

1.3 森林資源特征

西昌市森林植被類型復(fù)雜多樣，地帶性植被以殼斗科（Fagaceae）常綠闊葉林和云南松（Pinus yunnanensis）林為主。由滇青岡（Cyclobalanopsis glaucoides）、黃毛青岡（C.delavayi）和高山錐（Castanopsis delavayi）等組成常綠闊葉林，伴隨有少量落葉和硬葉的櫟屬（Quercus）或冬青屬（Ilex）植物，常與其他樹種形成塊狀分布的針闊混交林、常綠闊葉混交林和常綠落葉闊葉混交林。由于長期受人類活動的影響，常綠闊葉林呈不斷減少趨勢，逐漸被耐干旱和貧瘠的云南松林和華山松（P.armandii）林代替。主要闊葉樹種（組）有滇青岡、高山錐、桉樹、榿木（Alnus cremastogyne）、樺木屬（Betula）和楊屬（Populus）等；針葉樹種（組）有云南松、冷杉（Abies fabri）、華山松、柏木（Cupressus funebris）、杉木（Cunninghamia lanceolata）、柳杉（Cryptomeria japonicavar.sinensis）、云杉（Picea asperata）和鐵杉（Tsuga chinensis）等。

2 研究方法

2.1 抽樣總體設(shè)計

在森林資源規(guī)劃設(shè)計調(diào)查中，要以經(jīng)營單位或縣級行政單位為總體進(jìn)行總體活立木蓄積量抽樣控制[10]。抽樣調(diào)查是一門應(yīng)用性很強(qiáng)的學(xué)科,實施抽樣調(diào)查工作的要點有獲得正確的信息、省時省力省錢和充分利用已知信息提高抽樣精度[11]。根據(jù)西昌市國有林和集體林森林資源分布特點，結(jié)合上期二類調(diào)查總體劃分情況，將西昌市作為1 個調(diào)查總體，設(shè)置4 個副總體，分別為巴汝副總體（巴汝森林經(jīng)營所）、西寧-鹽中副總體（西寧飛播站和鹽中森林經(jīng)營所）、其他國有副總體（四合林場、大箐林場、瀘山森林經(jīng)營所和石嘉森林經(jīng)營所）和集體副總體（全市集體林區(qū)）。

2.2 抽樣方法選擇

系統(tǒng)抽樣是應(yīng)用較廣泛的抽樣控制方法。該方法操作簡便，樣本單位的獲取十分方便，樣本定位容易且分布均勻，對于總體有更好的代表性，能保證較好的抽樣效果，抽樣估計效率較高[12]。該方法從含有N個單元的總體中,隨機(jī)確定起點,按照預(yù)先設(shè)計的固定間隔模式抽取n 個樣本單元，組成樣本，估計總體[13]。綜合考慮森林資源特點和抽樣精度，西昌市采用隨機(jī)布點和機(jī)械抽取的系統(tǒng)抽樣方法。

2.3 精度與可靠性指標(biāo)

以公益林為主的經(jīng)營單位或縣級行政單位抽樣控制精度為75% ～85%，以幼中林和疏林地為主的精度為75%，以近、成和過熟林為主的精度為85%[14]。為保持全市調(diào)查精度的一致性，4個副總體抽樣精度均取85%，可靠性指標(biāo)取95%。

將總體抽樣調(diào)查方法計算的總體蓄積量與各小班目測蓄積量統(tǒng)計匯總的總體蓄積量進(jìn)行比較。當(dāng)兩者相差（累偏）小于±1 倍標(biāo)準(zhǔn)誤差時，認(rèn)為符合精度要求，蓄積統(tǒng)計分析直接使用小班調(diào)查數(shù)據(jù)；累偏在±1倍標(biāo)準(zhǔn)～±3倍標(biāo)準(zhǔn)誤差之間時，應(yīng)進(jìn)行小班調(diào)查檢查，找出誤差較大因素，對需要進(jìn)行現(xiàn)地核實的小班進(jìn)行野外調(diào)查，綜合分析后，對小班蓄積進(jìn)行修正；累偏大于±3 倍標(biāo)準(zhǔn)誤差時，應(yīng)對小班蓄積量進(jìn)行重新調(diào)查。

2.4 樣本單元面積

各總體抽樣調(diào)查選擇20 m × 20 m 的正方形地塊作為樣地。

2.5 樣本單元數(shù)估算

標(biāo)準(zhǔn)差有兩種計算方法。一是在調(diào)查區(qū)域內(nèi)借助現(xiàn)存樣本資料直接進(jìn)行標(biāo)準(zhǔn)差計算，此方法一般由于資料的缺乏較難完成；二是使用全距法計算標(biāo)準(zhǔn)差，此方法可通過預(yù)估樣本的大值和小值完成計算，在實際工作中應(yīng)用廣泛[13]。本研究采用全距法進(jìn)行計算。

2.5.1 蓄積平均值預(yù)估

根據(jù)西昌市上期二類調(diào)查資料結(jié)合本期小班調(diào)查結(jié)果，采用平均蓄積抽樣理論全距法進(jìn)行蓄積預(yù)估。各副總體預(yù)估平均蓄積見表1，各單位面積蓄積分布見圖1。

表1 各副總體預(yù)估平均蓄積Tab.1 Estimated average volume of vice overall areas (m3/hm2)

圖1 各副總體單位面積蓄積分布Fig.1 Distribution of volume unit area of vice overall areas

2.5.2 預(yù)估參數(shù)計算

各抽樣實測總體蓄積預(yù)估樣地數(shù)量均按100個進(jìn)行R值（全距比值）計算。其近似計算公式為：

式中，e為自然常數(shù)，N為樣本單元數(shù)，計算出R值為5.13。

各抽樣實測總體標(biāo)準(zhǔn)差（S）和變動系數(shù)（C）計算公式為：

式中，Y為單位面積蓄積。計算結(jié)果見表2。

根據(jù)可靠性指標(biāo)95%，查得閾值T0.05=2。抽樣精度（P）為85%，即E值（1-P）為15%。

表2 各副總體標(biāo)準(zhǔn)差與變動系數(shù)Tab.2 Standard deviation and coefficient of variation of vice overall areas

2.5.3 樣本單元數(shù)估算

樣本單元數(shù)（N）估算公式為：

結(jié)果見表3。

2.6 抽樣風(fēng)險分析

抽樣調(diào)查是在假設(shè)的基礎(chǔ)上，根據(jù)樣本的數(shù)量特征推斷總體的數(shù)量特征，通過樣本測定估算總體平均值，存在抽樣誤差和一定的風(fēng)險，可通過技術(shù)手段和措施把誤差和風(fēng)險降至最小，確保抽樣調(diào)查的精度。本次抽樣設(shè)計采用系統(tǒng)抽樣方法，減小了樣本單元周期性分布風(fēng)險。西昌市總體抽樣對象呈不規(guī)則分布，也大大降低了樣本單元周期性分布的可能性。同時，對各抽樣實測副總體采用不同方案進(jìn)行風(fēng)險控制，采取2 種方案對樣地數(shù)量進(jìn)行分析與確定。

方案1：樣地數(shù)量變動，正常計算樣地數(shù)量增加，即+10%、+15%保險系數(shù)。

方案2：平均蓄積的變動為主要因素，采用平均蓄積增減，即±5%、±10%、±15%蓄積進(jìn)行分析。

為切實保證抽樣方案安全可行，降低總體蓄積控制風(fēng)險，提高調(diào)查效率，確保調(diào)查質(zhì)量，分別對各副總體按方案1和方案2進(jìn)行分析，最終確定各抽樣實測副總體理論布設(shè)樣地采用方案及樣地數(shù)量（表4）。

表3 各副總體樣本單元數(shù)估算Tab.3 Estimation on number of sample plots in vice overall areas

表4 各抽樣實測副總體理論布設(shè)樣地方案及數(shù)量Tab.4 Scheme and quantity of theoretical sample plots of vice overall areas

2.7 樣本布設(shè)方案

2.7.1 抽樣范圍選取

在西昌市二類調(diào)查小班數(shù)據(jù)庫中，以坡度≤45°的中齡林以上（含中齡林）和平均直徑≥5 cm 的幼齡林的有林地和疏林地為條件選擇抽樣小班（表5）。

2.7.2 抽樣間隔與樣本單元確定

各總體應(yīng)用統(tǒng)一的抽樣間隔，以全市西北角公里坐標(biāo)（477 000，3 119 000）為起始值，按照1 km ×1 km、2 km ×2 km 和4 km ×4 km 間距進(jìn)行排序，然后導(dǎo)入到ARCGIS 中，利用添加XY 數(shù)據(jù)功能，生成抽樣布點圖，導(dǎo)出圖層數(shù)據(jù)。

將抽樣布點圖與抽樣范圍進(jìn)行疊加分析，其中2 km×2 km 布點方案與預(yù)估數(shù)量最接近，再將坐標(biāo)位置與抽樣范圍進(jìn)行比對，將樣地范圍內(nèi)蓄積面積較小的排除，剩余點為抽樣單元（表6）。

表5 各抽樣實測副總體布設(shè)樣地方案及數(shù)量Tab.5 Layout scheme and quantity of sample plots of vice overall area

表6 各抽樣實測副總體實際布設(shè)樣地方案及數(shù)量Tab.6 Layout scheme and quantity of actual sample plots of vice overall areas

2.7.3 樣本單元（樣地）布設(shè)

以小班地類為背景，添加行政區(qū)劃界線，標(biāo)注樣地號和總體號，經(jīng)整飾形成樣地布置圖。

2.8 補(bǔ)點方案

如樣地調(diào)查后達(dá)不到抽樣精度要求，則需補(bǔ)充樣地數(shù)量，并添加到原抽樣系統(tǒng)中。先計算間隔樣地數(shù)（已實測樣地數(shù)/需要補(bǔ)充樣地數(shù)），對間隔樣地數(shù)四舍五入取整形成間隔樣地組，第一個間隔樣地組和第二個間隔樣地組中間位置添加1 個樣地，以此類推。

2.9 樣地定位與周界測量

2.9.1 樣地定位

采用GPS 導(dǎo)航進(jìn)行樣地定位，定位坐標(biāo)為樣地西南角點。到達(dá)樣地位附近（30 m 左右），尋找一個明顯地物（大樹或獨立石頭等）為引點，連測引點至樣地。記入“羅盤儀測量記錄”。

2.9.2 樣地周界測量

用羅盤儀定方向，用測繩量距。樣地按坐標(biāo)方位角0°、90°、180°和270°進(jìn)行樣地周界測量?，F(xiàn)地周界閉合誤差≥1/200。方位角允許誤差不得大于30'。在角點埋設(shè)標(biāo)志，并書寫編號；標(biāo)志直徑大于2 cm，高度不低于10 cm，埋入地下部分不少于10 cm。

3 研究結(jié)果

3.1 實測樣地數(shù)

根據(jù)西昌市蓄積實測抽樣結(jié)果統(tǒng)計,全市共布設(shè)樣地242 個，實際收回樣地234 個，放棄樣地8 個（表7）。

3.2 實測抽樣精度

根據(jù)計算,抽樣精度滿足調(diào)查要求（表8）。

表7 各抽樣實測副總體樣地數(shù)Tab.7 Number of actual sample plots of vice overall areas

表8 各抽樣實測副總體蓄積特征值Tab.8 Volume characteristic values of vice overall areas

續(xù)表8 Continued

3.3 小班調(diào)查蓄積與總體蓄積對比

各小班目測蓄積量統(tǒng)計匯總的總體蓄積量為各小班目測蓄積量的累加值；總體抽樣調(diào)查方法計算的總體蓄積量為實測樣地平均值乘以可抽樣對象面積計算出的總體蓄積量。通過比較與修正，各抽樣實測副總體累偏均在±1倍標(biāo)準(zhǔn)誤差之內(nèi)。

4 結(jié)論與討論

4.1 總體劃分應(yīng)綜合考慮資源分布特征與樣地數(shù)量

森林資源分布特征是總體劃分的重要條件。西昌市集體林森林資源分布廣泛且零散，經(jīng)營水平基本相同，可將集體林劃為一個副總體。西昌市的國有林7個經(jīng)營單位也可按照自然地理特征和森林資源特點進(jìn)一步劃分副總體。從地形上看，巴汝森林經(jīng)營所、西寧飛播站和鹽中森林經(jīng)營所、四合林場、大箐林場、瀘山森林經(jīng)營所和石嘉森林經(jīng)營所縱向分列于西昌市西、中和東部。西部的牦牛山，是安寧河與雅礱江的分水嶺，其東邊是西寧飛播站和鹽中森林經(jīng)營所，西邊則是巴汝森林經(jīng)營所；四合林場、大箐林場、瀘山森林經(jīng)營所和石嘉森林經(jīng)營所均位于安寧河?xùn)|側(cè)的螺髻山脈。由于所處自然地理條件的不同，國有森林資源地帶性分異明顯，可將國有林劃分為3個副總體類型?？傮w劃分也應(yīng)考慮樣地數(shù)量，總體數(shù)量劃分越多，理論上抽樣精度越容易控制，但同時將極大增加外業(yè)調(diào)查工作，降低調(diào)查工作效率，因此樣地數(shù)量也是總體劃分需要考慮的重要因素。

4.2 合理設(shè)置蓄積預(yù)估值最大和最小值

總體蓄積預(yù)估一般參考的是上期二類調(diào)查成果或者本次二類調(diào)查小班目測調(diào)查結(jié)果，因此預(yù)估的平均蓄積均存在一定誤差。為了提高精度，降低誤差，應(yīng)對蓄積預(yù)估最大和最小值進(jìn)行分析，不能簡單地將小班中公頃蓄積最大和最小的設(shè)為預(yù)估最大和最小值。巴汝副總體蓄積中有290 ～305 m3/hm2區(qū)間的小班，此區(qū)間面積占比僅0.08%，因此預(yù)估最大蓄積取285 m3/hm2較合理；集體副總體蓄積中185 ～210 m3/hm2區(qū)間的小班，面積占比也較小，因此預(yù)估最大蓄積取185 m3/hm2。

4.3 標(biāo)準(zhǔn)地目測調(diào)查訓(xùn)練十分重要

通過對比小班目測蓄積與總體樣地實測蓄積，少量林班存在目測調(diào)查蓄積整體偏小的情況。由此對這部分林班調(diào)查情況展開了進(jìn)一步分析，發(fā)現(xiàn)這部分林班由1 個調(diào)查工組完成。經(jīng)仔細(xì)分析，造成調(diào)查結(jié)果存在系統(tǒng)性偏差的主要原因為該工組標(biāo)準(zhǔn)地目測調(diào)查訓(xùn)練時未全面掌握技術(shù)標(biāo)準(zhǔn)，導(dǎo)致目測蓄積整體偏小。由此可見，標(biāo)準(zhǔn)地目測調(diào)查訓(xùn)練在二類調(diào)查工作中十分重要，只有熟練且準(zhǔn)確地掌握小班目測調(diào)查方法和技術(shù)標(biāo)準(zhǔn)，才能保證調(diào)查精度。

4.4 采用隨機(jī)布點和機(jī)械抽取減少人為干擾因素

采用隨機(jī)布點和機(jī)械抽取的抽樣布點方案，排除了人為干擾因素的影響，樣本單元均由軟件隨機(jī)布設(shè)，符合樣地機(jī)械設(shè)置原則。通過分析小班目測蓄積與樣地實測蓄積，實測樣地基本覆蓋了各個總體低、中、高蓄積區(qū)間。隨機(jī)布點和機(jī)械抽取的系統(tǒng)抽樣布點方法符合總體蓄積抽樣設(shè)計原則，能達(dá)到總體蓄積控制的目的和要求。