Excel在林地林木资源价值核算调查抽样中的应用

王森林，王宏年，罗军伟，吴 可，王爱东，秦乃花，张 鹏

（山东省林业监测规划院，山东 济南 250014）

2016年以来，国家林业局与国家统计局联合开展了中国森林资源核算研究，旨在通过研究推动我国森林资源核算常态化，建立和完善符合生态文明建设需要的中国森林资源核算和林业绿色经济评价体系。党的十九大强调坚持人与自然和谐共生，树立绿水青山就是金山银山的理念，为进一步深化森林资源核算工作提出了新任务，开展林地林木资源核算是中国森林资源核算研究的重要内容，林地林木资源核算包括实物量账户核算、价值量账户核算，从蓄积量核算分天然林资源、人工林资源和其它林木。林地林木资源核算，通过国民经济核算的方法和手段，可反映森林在一个核算期内的存量和变动情况，对发生变动的原因进行量化分析。新一轮核算对于完善生态补偿政策和编制自然资源资产负债表具有重要参考价值。根据国家林业局科技司《关于开展林地林木资源价值核算调查的函》（科综字〔2017〕）33号要求，今年4月份林业系统在全国范围内组织开展了林地林木价值核算调查工作，调查样本县、乡（镇）、国有林场确定如何，关系最终成果质量。现就调查内容、抽样方法做一介绍。

1 调查内容

包括林地价值核算调查、林木价值核算调查两个方面内容，共涉及7个调查表格。

1.1 林地价值核算调查

林地价值核算采用年金资本化法估价，林地使用权市场交易价格作为参考。

表1-1：调查抽样单元内不同类型林地的当年租金价格、林地面积（资源量）。

表1-2：调查交易案例的单位面积林地使用权交易价格、使用权流转期限、林地立地质量、流转方式、流转发生时间等。

1.2 林木价值核算调查

根据林木价值核算选择的估价方法，调查内容分为活立木交易价格调查和林木生产成本收益类调查。

表2-1：对调查单元内的活立木（乔木林、经济林、竹林）交易案例进行调查，调查树种（组）、龄组、林木交易单位价格，案例单位面积蓄积。

表2-2、表2-3：调查单元内的不同树种的营造林成本、产品价格、产量等情况；调查乔木林不同树种的中幼龄林营造林成本、间伐成本收益、木材市场价格、林产品价格、采运成本、税费、利润、单位面积产量、单位面积蓄积量和轮伐期等指标。

表2-4：调查竹林的稳产前期、稳产期的生产成本和收益。

表2-5：调查经济林的初产期、盛产期的生产成本和收益。

2 *抽样方法说明（以江西为例）

采用二阶不等概抽样方法确定调查样本单元。

2.1 确定样本县（区）级单元

此次调查全国抽取501个乡（镇）、国有林场作为调查对象。按照省级林地面积占全国林地面积比例分配省级抽取样本数。样本数分配结果显示，江西省分配18个样本单元。按照每个县（区）抽3个的样本点的要求，江西省调查需要抽选6个县（区）（用n表示）。

M0=，按照江西省县（市、区）正常顺序对江西省各县（区）进行编号，并根据江西省最新的森林资源二类调查数据列出所有县（区）的林地面积和累计林地面积，详见表1。

Mi代表样本单元的林地面积，M0代表样本单元林地累计面积，N代表总体数115。以样本单元林地累计面积M0除以抽中样本县（区）数，即k=M0/n=10720217.6/6=1786703，从[1，k]中随机抽取一个整数如 1225653，这样 r=1225653，r+k=3012356，r+2k=4799059……r+5k=10159168，查找林地面积累计数与此6个数值接近的数，确定对应的县名，即修水县、崇义县、会昌县、玉山县、永丰县、宜黄县入样，序号分别为 29，53，62，82，99，111。

表1 江西省县级样本抽样结果 单位：hm2

2.2 确定乡（镇）、国有林场样本单元

采用等距抽样的方法抽取样本乡（镇）和国有林场。即随机起点等距抽样的方法，每个乡（镇）、国有林场被抽中的概率相同。按照乡（镇）距离县城的距离远近，对乡（镇）和国有林场依次进行编号。根据总体数和样本数，确定抽样间隔k。

乡（镇）和国有林场个数N=20，每县（区）抽取单位n=3，计算间隔数k=N/n=20/6=6.66，取整数k=7。在[1,k]中随机抽取 1 个整数 r，如 r=6，则 r+k、r+2k分别为13、20，对应的代码分别是双明镇、紫湖镇大叶林场和国营玉山县怀玉山林场，即为抽中调查样本。

表2 抽样的乡（镇）国有林场结果 单位：km

3 excel抽样（以山东为例）

3.1 确定样本县（区）级单元

按照2.1确定样本县（区）级单元抽样步骤，编制县（区）级单元抽样表，见表3。现对表格的栏目加以说明。

1）序号栏：B7：B143，1—137，指山东省 137 个县（区、市）顺序号，有省略，下同；

2）统计单位栏：C7：C143，历下区—东明县，指山东省济南市、青岛市、淄博市、枣庄市、东营市、烟台市、潍坊市、济宁市、泰安市、威海市、日照市、莱芜市、临沂市、德州市、聊城市、滨州市、菏泽市17个市137个县（区、市）具体名称；

3）林地面积栏：D7：D143，2644.49—9521.88，指山东省137个县（区、市）2016年度森林资源动态监测林地面积，单位为hm2；

4）累计栏：E7：E143，2644.49—3510693.7，指山东省137个县（区、市）2016年度森林资源动态监测林地面积依次累计，单位为hm2；

5）E3栏：3510693.7，指山东省137个县 （区、市）2016年度森林资源动态监测林地面积依次累计，单位为公顷，以E2栏字母M0表示，也是E7:E143 栏最大值，即“=MAX（E7：E143）”；

6）F4 栏：6，指某省抽中县（区）级单元数，按照国家林业局抽样要求，最多抽13个单位，最少抽1个单位；

7）E6 栏：585116，指抽样数据间隔数，是 E3 栏数据除以F4栏数据取整，即 “=ROUND （$E$3/$F$4，0）”，以 E4 栏字母 k 代表，$E$3、$F$4 是绝对引用E3栏、F4栏数据；

8）H2 栏：指在（1,k）之间随机选取的整数，以公式 “=RANDBETWEEN（1，$E$5)）”表示，把鼠标放在公式窗口，回车，即出现随机数值，以I2栏字母r代表；

9）K2栏：423166，可以设置为与H2栏数据相等的随机数r，也可以是手工填写的任意数，后面一系列计算及抽中县（区）单元都与此数有直接关系，如此数是423166时抽中的县（区）单元即是平度市、长岛市、青州市、环翠区、兰陵县、沾化县；

10）F6 栏：累计-r，F7：F143，指林地面积累计数与随机数之差绝对值，如F7栏数据420521.51=2644.49-423166，以公式“=ABS（E7-$K$2）”表示；

11）F5栏：4403.14，指林地面积累计数与随机数之差绝对值最小的数，也就是说林地面积累计数与随机数最接近的数，以公式“=MIN（F7:F143）”；

12）B3 栏：平度市，指通过 LOOKUP 在 F7：F143栏找到与F5栏一致的数据，然后再到C7：C143中找到对应的县（区），以公式“=LOOKUP（1，0/（F7：F143=F5），C7：C143）”表示；

13）C3、D3、B5、C5、D5 栏：长岛市、青州市、环翠区、兰陵县、沾化县，指分别通过LOOKUP在G7：G143、H7：H143、I7：I143、J7：J143、K7：K143 栏找到分别与 G5、H5、I5、J5、K5 栏一致的数据，然后再到 C7：C143中找到对应的县 （区），分别以公式“=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）序号栏：1-44，指抽中县级单元组数，多次抽取便于比选，有省略，下同，见表4；

15）r栏：随机数栏，直接从与表3中引用，以公式 “=$K$2”表示，r数值是随机变化的，如在r=423166时抽中的县（区）单元即是平度市、长岛市、青州市、环翠区、兰陵县、沾化县，如果把抽中这一组单元复制粘贴为数值，下一行就出现r=495562，对应的县（区）即是张店区、莱阳市、高密市、文登市、费县、曹县；

16）k栏：间隔数，直接从表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 栏：根据随机数 r与间隔数k推算数据，根据此数据查找林地面积接近的数，进一步确定抽取县（区）级调查单元。

3.2 确定乡（镇）、国有林场样本单元

按照2.2确定乡（镇）、国有林场样本单元抽样步骤，以临沂市沂水县为例编制乡（镇）、国有林场样本单元抽样表，见表4。现对表格的栏目由左至右顺序加以说明。

1）序号栏：B7：B28，1-22，指依据沂水县 22 个乡（镇）、国有林场距离县政府的近远顺序编号；

2）统计单位栏：C7：C28，指依据沂水县 22 个乡（镇）、国有林场具体名称；

3）离县城距离：D7：D28，指依据沂水县 22 个乡（镇）、国有林场距离县城公里数值；

4）E6 栏：7，指抽样数据间隔数，全县乡（镇）、国有林场总数除以全县抽取数，以公式“=ROUND(COUNT(D7:D28)/3,0)”表示，以F6栏字母k代表；

5）E4栏：指在（1,k）之间随机选取的整数，以公式 “=RANDBETWEEN（1，$E$6)）”表示，把鼠标放在公式窗口，回车，即出现随机数值，以F4栏字母r代表；

表3 县级单元抽样表 单位：hm2

表4 县级单元抽样表

6）E5栏：6，可以设置为与E4栏数据相等的随机数r，也可以是手工填写的任意数，后面一系列计算及抽中乡（镇）、国有林场样本单元都与此数有直接关系；

7）B3 栏：道托镇，指通过 LOOKUP 在 B7：B28栏找到与E5栏一致的数据，然后再到C7：C28中找到对应的乡 （镇）、国有林场，以公式“=LOOKUP（1，0/（B7：B28=E5），C7：C28）”表示；

表5 乡（镇）、国有林场样本单元抽样表

8）C3、D3栏：泉庄镇、沂山林场，指通过LOOKUP 在 B7：B28 栏找到分别与 （E5+E6）、（E5+2*E6）一致的数据，然后再到C7：C28中找到对应的乡（镇）、国有林场，分别以公式“=LOOKUP（1，0/（B7：B28=（E5+E6 栏 ）），C7：C28）”、 “=LOOKUP（1，0/（B7：B28=（E5+2*E6）），C7：C28）”表示；

9）序号栏：1-22，指抽中乡（镇）、国有林场单元组数，多次抽取便于比选；

10）r栏：随机数栏，r数值是随机变化的，如在r=6时抽中的乡（镇）、国有林场样本单元即是道托镇、泉庄镇、沂山林场，如果把抽中这一组单元复制粘贴为数值，下一行就出现r=3，对应的乡（镇）、国有林场样本单元即是龙家圈镇、高桥镇、马站镇；

11）抽中单位1-3栏：直接从表左面引用，分别以公式“=$B$3”、“=$C$3”、“=$D$3”表示；

12）r+k、r+2k栏：根据随机数r与间隔数k推算数据，根据此数据查找与此数对应的乡（镇）、国有林场序号，确定抽取乡（镇）、国有林场调查单元。

4 讨论与建议

4.1 本文仅编写了6个县（区）级调查单元抽取方法，望读者能够举一反三编制更多的县（区）级调查单元抽取方法。

4.2 林地林木价值核算很有必要，建议山东林业主管部门领导，近期按照抽样有关要求，在全省举行林地林木价值核算全面调查。