Excel在林地林木資源價(jià)值核算調(diào)查抽樣中的應(yīng)用

王森林，王宏年，羅軍偉，吳 可，王愛(ài)東，秦乃花，張 鵬

（山東省林業(yè)監(jiān)測(cè)規(guī)劃院，山東 濟(jì)南 250014）

2016年以來(lái)，國(guó)家林業(yè)局與國(guó)家統(tǒng)計(jì)局聯(lián)合開(kāi)展了中國(guó)森林資源核算研究，旨在通過(guò)研究推動(dòng)我國(guó)森林資源核算常態(tài)化，建立和完善符合生態(tài)文明建設(shè)需要的中國(guó)森林資源核算和林業(yè)綠色經(jīng)濟(jì)評(píng)價(jià)體系。黨的十九大強(qiáng)調(diào)堅(jiān)持人與自然和諧共生，樹(shù)立綠水青山就是金山銀山的理念，為進(jìn)一步深化森林資源核算工作提出了新任務(wù)，開(kāi)展林地林木資源核算是中國(guó)森林資源核算研究的重要內(nèi)容，林地林木資源核算包括實(shí)物量賬戶核算、價(jià)值量賬戶核算，從蓄積量核算分天然林資源、人工林資源和其它林木。林地林木資源核算，通過(guò)國(guó)民經(jīng)濟(jì)核算的方法和手段，可反映森林在一個(gè)核算期內(nèi)的存量和變動(dòng)情況，對(duì)發(fā)生變動(dòng)的原因進(jìn)行量化分析。新一輪核算對(duì)于完善生態(tài)補(bǔ)償政策和編制自然資源資產(chǎn)負(fù)債表具有重要參考價(jià)值。根據(jù)國(guó)家林業(yè)局科技司《關(guān)于開(kāi)展林地林木資源價(jià)值核算調(diào)查的函》（科綜字〔2017〕）33號(hào)要求，今年4月份林業(yè)系統(tǒng)在全國(guó)范圍內(nèi)組織開(kāi)展了林地林木價(jià)值核算調(diào)查工作，調(diào)查樣本縣、鄉(xiāng)（鎮(zhèn)）、國(guó)有林場(chǎng)確定如何，關(guān)系最終成果質(zhì)量?，F(xiàn)就調(diào)查內(nèi)容、抽樣方法做一介紹。

1 調(diào)查內(nèi)容

包括林地價(jià)值核算調(diào)查、林木價(jià)值核算調(diào)查兩個(gè)方面內(nèi)容，共涉及7個(gè)調(diào)查表格。

1.1 林地價(jià)值核算調(diào)查

林地價(jià)值核算采用年金資本化法估價(jià)，林地使用權(quán)市場(chǎng)交易價(jià)格作為參考。

表1-1：調(diào)查抽樣單元內(nèi)不同類型林地的當(dāng)年租金價(jià)格、林地面積（資源量）。

表1-2：調(diào)查交易案例的單位面積林地使用權(quán)交易價(jià)格、使用權(quán)流轉(zhuǎn)期限、林地立地質(zhì)量、流轉(zhuǎn)方式、流轉(zhuǎn)發(fā)生時(shí)間等。

1.2 林木價(jià)值核算調(diào)查

根據(jù)林木價(jià)值核算選擇的估價(jià)方法，調(diào)查內(nèi)容分為活立木交易價(jià)格調(diào)查和林木生產(chǎn)成本收益類調(diào)查。

表2-1：對(duì)調(diào)查單元內(nèi)的活立木（喬木林、經(jīng)濟(jì)林、竹林）交易案例進(jìn)行調(diào)查，調(diào)查樹(shù)種（組）、齡組、林木交易單位價(jià)格，案例單位面積蓄積。

表2-2、表2-3：調(diào)查單元內(nèi)的不同樹(shù)種的營(yíng)造林成本、產(chǎn)品價(jià)格、產(chǎn)量等情況；調(diào)查喬木林不同樹(shù)種的中幼齡林營(yíng)造林成本、間伐成本收益、木材市場(chǎng)價(jià)格、林產(chǎn)品價(jià)格、采運(yùn)成本、稅費(fèi)、利潤(rùn)、單位面積產(chǎn)量、單位面積蓄積量和輪伐期等指標(biāo)。

表2-4：調(diào)查竹林的穩(wěn)產(chǎn)前期、穩(wěn)產(chǎn)期的生產(chǎn)成本和收益。

表2-5：調(diào)查經(jīng)濟(jì)林的初產(chǎn)期、盛產(chǎn)期的生產(chǎn)成本和收益。

2 *抽樣方法說(shuō)明（以江西為例）

采用二階不等概抽樣方法確定調(diào)查樣本單元。

2.1 確定樣本縣（區(qū)）級(jí)單元

此次調(diào)查全國(guó)抽取501個(gè)鄉(xiāng)（鎮(zhèn)）、國(guó)有林場(chǎng)作為調(diào)查對(duì)象。按照省級(jí)林地面積占全國(guó)林地面積比例分配省級(jí)抽取樣本數(shù)。樣本數(shù)分配結(jié)果顯示，江西省分配18個(gè)樣本單元。按照每個(gè)縣（區(qū)）抽3個(gè)的樣本點(diǎn)的要求，江西省調(diào)查需要抽選6個(gè)縣（區(qū)）（用n表示）。

M0=，按照江西省縣（市、區(qū)）正常順序?qū)魇「骺h（區(qū)）進(jìn)行編號(hào)，并根據(jù)江西省最新的森林資源二類調(diào)查數(shù)據(jù)列出所有縣（區(qū)）的林地面積和累計(jì)林地面積，詳見(jiàn)表1。

Mi代表樣本單元的林地面積，M0代表樣本單元林地累計(jì)面積，N代表總體數(shù)115。以樣本單元林地累計(jì)面積M0除以抽中樣本縣（區(qū)）數(shù)，即k=M0/n=10720217.6/6=1786703，從[1，k]中隨機(jī)抽取一個(gè)整數(shù)如 1225653，這樣 r=1225653，r+k=3012356，r+2k=4799059……r+5k=10159168，查找林地面積累計(jì)數(shù)與此6個(gè)數(shù)值接近的數(shù)，確定對(duì)應(yīng)的縣名，即修水縣、崇義縣、會(huì)昌縣、玉山縣、永豐縣、宜黃縣入樣，序號(hào)分別為 29，53，62，82，99，111。

表1 江西省縣級(jí)樣本抽樣結(jié)果 單位：hm2

2.2 確定鄉(xiāng)（鎮(zhèn)）、國(guó)有林場(chǎng)樣本單元

采用等距抽樣的方法抽取樣本鄉(xiāng)（鎮(zhèn)）和國(guó)有林場(chǎng)。即隨機(jī)起點(diǎn)等距抽樣的方法，每個(gè)鄉(xiāng)（鎮(zhèn)）、國(guó)有林場(chǎng)被抽中的概率相同。按照鄉(xiāng)（鎮(zhèn)）距離縣城的距離遠(yuǎn)近，對(duì)鄉(xiāng)（鎮(zhèn)）和國(guó)有林場(chǎng)依次進(jìn)行編號(hào)。根據(jù)總體數(shù)和樣本數(shù)，確定抽樣間隔k。

鄉(xiāng)（鎮(zhèn)）和國(guó)有林場(chǎng)個(gè)數(shù)N=20，每縣（區(qū)）抽取單位n=3，計(jì)算間隔數(shù)k=N/n=20/6=6.66，取整數(shù)k=7。在[1,k]中隨機(jī)抽取 1 個(gè)整數(shù) r，如 r=6，則 r+k、r+2k分別為13、20，對(duì)應(yīng)的代碼分別是雙明鎮(zhèn)、紫湖鎮(zhèn)大葉林場(chǎng)和國(guó)營(yíng)玉山縣懷玉山林場(chǎng)，即為抽中調(diào)查樣本。

表2 抽樣的鄉(xiāng)（鎮(zhèn)）國(guó)有林場(chǎng)結(jié)果 單位：km

3 excel抽樣（以山東為例）

3.1 確定樣本縣（區(qū)）級(jí)單元

按照2.1確定樣本縣（區(qū)）級(jí)單元抽樣步驟，編制縣（區(qū)）級(jí)單元抽樣表，見(jiàn)表3?，F(xiàn)對(duì)表格的欄目加以說(shuō)明。

1）序號(hào)欄：B7：B143，1—137，指山東省 137 個(gè)縣（區(qū)、市）順序號(hào)，有省略，下同；

2）統(tǒng)計(jì)單位欄：C7：C143，歷下區(qū)—東明縣，指山東省濟(jì)南市、青島市、淄博市、棗莊市、東營(yíng)市、煙臺(tái)市、濰坊市、濟(jì)寧市、泰安市、威海市、日照市、萊蕪市、臨沂市、德州市、聊城市、濱州市、菏澤市17個(gè)市137個(gè)縣（區(qū)、市）具體名稱；

3）林地面積欄：D7：D143，2644.49—9521.88，指山東省137個(gè)縣（區(qū)、市）2016年度森林資源動(dòng)態(tài)監(jiān)測(cè)林地面積，單位為hm2；

4）累計(jì)欄：E7：E143，2644.49—3510693.7，指山東省137個(gè)縣（區(qū)、市）2016年度森林資源動(dòng)態(tài)監(jiān)測(cè)林地面積依次累計(jì)，單位為hm2；

5）E3欄：3510693.7，指山東省137個(gè)縣 （區(qū)、市）2016年度森林資源動(dòng)態(tài)監(jiān)測(cè)林地面積依次累計(jì)，單位為公頃，以E2欄字母M0表示，也是E7:E143 欄最大值，即“=MAX（E7：E143）”；

6）F4 欄：6，指某省抽中縣（區(qū)）級(jí)單元數(shù)，按照國(guó)家林業(yè)局抽樣要求，最多抽13個(gè)單位，最少抽1個(gè)單位；

7）E6 欄：585116，指抽樣數(shù)據(jù)間隔數(shù)，是 E3 欄數(shù)據(jù)除以F4欄數(shù)據(jù)取整，即 “=ROUND （$E$3/$F$4，0）”，以 E4 欄字母 k 代表，$E$3、$F$4 是絕對(duì)引用E3欄、F4欄數(shù)據(jù)；

8）H2 欄：指在（1,k）之間隨機(jī)選取的整數(shù)，以公式 “=RANDBETWEEN（1，$E$5)）”表示，把鼠標(biāo)放在公式窗口，回車，即出現(xiàn)隨機(jī)數(shù)值，以I2欄字母r代表；

9）K2欄：423166，可以設(shè)置為與H2欄數(shù)據(jù)相等的隨機(jī)數(shù)r，也可以是手工填寫的任意數(shù)，后面一系列計(jì)算及抽中縣（區(qū)）單元都與此數(shù)有直接關(guān)系，如此數(shù)是423166時(shí)抽中的縣（區(qū)）單元即是平度市、長(zhǎng)島市、青州市、環(huán)翠區(qū)、蘭陵縣、沾化縣；

10）F6 欄：累計(jì)-r，F(xiàn)7：F143，指林地面積累計(jì)數(shù)與隨機(jī)數(shù)之差絕對(duì)值，如F7欄數(shù)據(jù)420521.51=2644.49-423166，以公式“=ABS（E7-$K$2）”表示；

11）F5欄：4403.14，指林地面積累計(jì)數(shù)與隨機(jī)數(shù)之差絕對(duì)值最小的數(shù)，也就是說(shuō)林地面積累計(jì)數(shù)與隨機(jī)數(shù)最接近的數(shù)，以公式“=MIN（F7:F143）”；

12）B3 欄：平度市，指通過(guò) LOOKUP 在 F7：F143欄找到與F5欄一致的數(shù)據(jù)，然后再到C7：C143中找到對(duì)應(yīng)的縣（區(qū)），以公式“=LOOKUP（1，0/（F7：F143=F5），C7：C143）”表示；

13）C3、D3、B5、C5、D5 欄：長(zhǎng)島市、青州市、環(huán)翠區(qū)、蘭陵縣、沾化縣，指分別通過(guò)LOOKUP在G7：G143、H7：H143、I7：I143、J7：J143、K7：K143 欄找到分別與 G5、H5、I5、J5、K5 欄一致的數(shù)據(jù)，然后再到 C7：C143中找到對(duì)應(yīng)的縣 （區(qū)），分別以公式“=LOOKUP（1，0/（G7：G143=G5），C7：C143）”、“=LOOKUP （1，0/（H7：H143=F5），C7：C143）”、 “=LOOKUP （1，0/（I7：I143=I5），C7：C143）”、“=LOOKUP （1，0/（J7：J143=J5），C7：C143）”、 “=LOOKUP （1，0/（K7：K143=K5），C7：C143）”表示；

14）序號(hào)欄：1-44，指抽中縣級(jí)單元組數(shù)，多次抽取便于比選，有省略，下同，見(jiàn)表4；

15）r欄：隨機(jī)數(shù)欄，直接從與表3中引用，以公式 “=$K$2”表示，r數(shù)值是隨機(jī)變化的，如在r=423166時(shí)抽中的縣（區(qū)）單元即是平度市、長(zhǎng)島市、青州市、環(huán)翠區(qū)、蘭陵縣、沾化縣，如果把抽中這一組單元復(fù)制粘貼為數(shù)值，下一行就出現(xiàn)r=495562，對(duì)應(yīng)的縣（區(qū)）即是張店區(qū)、萊陽(yáng)市、高密市、文登市、費(fèi)縣、曹縣；

16）k欄：間隔數(shù)，直接從表3中引用，以公式“=$E$5”表示；

17）抽中單位1-6欄：直接從表3中引用，分別以公式“=$B$3”、“=$C$3”、“=$D$3”、“=$B$5”、“=$C$5”、“=$D$5”表示；

18）r+k、r+2k、r+3k、r+4k、r+5k 欄：根據(jù)隨機(jī)數(shù) r與間隔數(shù)k推算數(shù)據(jù)，根據(jù)此數(shù)據(jù)查找林地面積接近的數(shù)，進(jìn)一步確定抽取縣（區(qū)）級(jí)調(diào)查單元。

3.2 確定鄉(xiāng)（鎮(zhèn)）、國(guó)有林場(chǎng)樣本單元

按照2.2確定鄉(xiāng)（鎮(zhèn)）、國(guó)有林場(chǎng)樣本單元抽樣步驟，以臨沂市沂水縣為例編制鄉(xiāng)（鎮(zhèn)）、國(guó)有林場(chǎng)樣本單元抽樣表，見(jiàn)表4?，F(xiàn)對(duì)表格的欄目由左至右順序加以說(shuō)明。

1）序號(hào)欄：B7：B28，1-22，指依據(jù)沂水縣 22 個(gè)鄉(xiāng)（鎮(zhèn)）、國(guó)有林場(chǎng)距離縣政府的近遠(yuǎn)順序編號(hào)；

2）統(tǒng)計(jì)單位欄：C7：C28，指依據(jù)沂水縣 22 個(gè)鄉(xiāng)（鎮(zhèn)）、國(guó)有林場(chǎng)具體名稱；

3）離縣城距離：D7：D28，指依據(jù)沂水縣 22 個(gè)鄉(xiāng)（鎮(zhèn)）、國(guó)有林場(chǎng)距離縣城公里數(shù)值；

4）E6 欄：7，指抽樣數(shù)據(jù)間隔數(shù)，全縣鄉(xiāng)（鎮(zhèn)）、國(guó)有林場(chǎng)總數(shù)除以全縣抽取數(shù)，以公式“=ROUND(COUNT(D7:D28)/3,0)”表示，以F6欄字母k代表；

5）E4欄：指在（1,k）之間隨機(jī)選取的整數(shù)，以公式 “=RANDBETWEEN（1，$E$6)）”表示，把鼠標(biāo)放在公式窗口，回車，即出現(xiàn)隨機(jī)數(shù)值，以F4欄字母r代表；

表3 縣級(jí)單元抽樣表 單位：hm2

表4 縣級(jí)單元抽樣表

6）E5欄：6，可以設(shè)置為與E4欄數(shù)據(jù)相等的隨機(jī)數(shù)r，也可以是手工填寫的任意數(shù)，后面一系列計(jì)算及抽中鄉(xiāng)（鎮(zhèn)）、國(guó)有林場(chǎng)樣本單元都與此數(shù)有直接關(guān)系；

7）B3 欄：道托鎮(zhèn)，指通過(guò) LOOKUP 在 B7：B28欄找到與E5欄一致的數(shù)據(jù)，然后再到C7：C28中找到對(duì)應(yīng)的鄉(xiāng) （鎮(zhèn)）、國(guó)有林場(chǎng)，以公式“=LOOKUP（1，0/（B7：B28=E5），C7：C28）”表示；

表5 鄉(xiāng)（鎮(zhèn)）、國(guó)有林場(chǎng)樣本單元抽樣表

8）C3、D3欄：泉莊鎮(zhèn)、沂山林場(chǎng)，指通過(guò)LOOKUP 在 B7：B28 欄找到分別與 （E5+E6）、（E5+2*E6）一致的數(shù)據(jù)，然后再到C7：C28中找到對(duì)應(yīng)的鄉(xiāng)（鎮(zhèn)）、國(guó)有林場(chǎng)，分別以公式“=LOOKUP（1，0/（B7：B28=（E5+E6 欄 ）），C7：C28）”、 “=LOOKUP（1，0/（B7：B28=（E5+2*E6）），C7：C28）”表示；

9）序號(hào)欄：1-22，指抽中鄉(xiāng)（鎮(zhèn)）、國(guó)有林場(chǎng)單元組數(shù)，多次抽取便于比選；

10）r欄：隨機(jī)數(shù)欄，r數(shù)值是隨機(jī)變化的，如在r=6時(shí)抽中的鄉(xiāng)（鎮(zhèn)）、國(guó)有林場(chǎng)樣本單元即是道托鎮(zhèn)、泉莊鎮(zhèn)、沂山林場(chǎng)，如果把抽中這一組單元復(fù)制粘貼為數(shù)值，下一行就出現(xiàn)r=3，對(duì)應(yīng)的鄉(xiāng)（鎮(zhèn)）、國(guó)有林場(chǎng)樣本單元即是龍家圈鎮(zhèn)、高橋鎮(zhèn)、馬站鎮(zhèn)；

11）抽中單位1-3欄：直接從表左面引用，分別以公式“=$B$3”、“=$C$3”、“=$D$3”表示；

12）r+k、r+2k欄：根據(jù)隨機(jī)數(shù)r與間隔數(shù)k推算數(shù)據(jù)，根據(jù)此數(shù)據(jù)查找與此數(shù)對(duì)應(yīng)的鄉(xiāng)（鎮(zhèn)）、國(guó)有林場(chǎng)序號(hào)，確定抽取鄉(xiāng)（鎮(zhèn)）、國(guó)有林場(chǎng)調(diào)查單元。

4 討論與建議

4.1 本文僅編寫了6個(gè)縣（區(qū)）級(jí)調(diào)查單元抽取方法，望讀者能夠舉一反三編制更多的縣（區(qū)）級(jí)調(diào)查單元抽取方法。

4.2 林地林木價(jià)值核算很有必要，建議山東林業(yè)主管部門領(lǐng)導(dǎo)，近期按照抽樣有關(guān)要求，在全省舉行林地林木價(jià)值核算全面調(diào)查。