留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于度量误差方法的油松林分生长模型

王金池 冉啟香 邓华锋 黄国胜 王雪军

王金池, 冉啟香, 邓华锋, 黄国胜, 王雪军. 基于度量误差方法的油松林分生长模型[J]. 浙江农林大学学报, 2018, 35(1): 68-74. doi: 10.11833/j.issn.2095-0756.2018.01.009
引用本文: 王金池, 冉啟香, 邓华锋, 黄国胜, 王雪军. 基于度量误差方法的油松林分生长模型[J]. 浙江农林大学学报, 2018, 35(1): 68-74. doi: 10.11833/j.issn.2095-0756.2018.01.009
WANG Jinchi, RAN Qixiang, DENG Huafeng, HUANG Guosheng, WANG Xuejun. Stand growth model for Pinus tabulaeformis using measurement error method[J]. Journal of Zhejiang A&F University, 2018, 35(1): 68-74. doi: 10.11833/j.issn.2095-0756.2018.01.009
Citation: WANG Jinchi, RAN Qixiang, DENG Huafeng, HUANG Guosheng, WANG Xuejun. Stand growth model for Pinus tabulaeformis using measurement error method[J]. Journal of Zhejiang A&F University, 2018, 35(1): 68-74. doi: 10.11833/j.issn.2095-0756.2018.01.009

基于度量误差方法的油松林分生长模型

doi: 10.11833/j.issn.2095-0756.2018.01.009
基金项目: 

国家林业公益性行业科研专项 201204510

详细信息
    作者简介: 王金池, 从事森林资源监测与评价研究.E-mail:2276599635@qq.com
    通信作者: 邓华锋, 教授, 博士, 从事森林可持续经营理论与技术研究.E-mail:denghuafeng@bjfu.edu.cn
  • 中图分类号: S758.5

Stand growth model for Pinus tabulaeformis using measurement error method

  • 摘要: 为掌握北京地区油松Pinus tabulaeformis生长过程,建立相容的断面积和蓄积模型,通过引入间伐林分与未间伐林分的哑变量,分别建立林分断面积、林分蓄积生长模型,然后从相容的角度出发,建立林分断面积、蓄积量的误差变量联立方程组,并与不含哑变量的传统生长模型的误差变量联立方程组进行比较。经检验,传统误差变量联立方程组中林分断面积和林分蓄积生长模型的预测精度都在92%以上,对油松林分断面积的预测精度高达0.921 5,决定系数高达0.900 1,对油松林分蓄积量的预测精度达到了0.928 3,决定系数高达0.912 3,而引入哑变量的误差变量联立方程组中,模型的预测精度和确定系数稍高,均在93%以上,对油松林分断面积的预测精度高达0.939 8,决定系数达到了0.927 9,对油松林分蓄积量的预测精度在0.930 0以上,决定系数达0.932 8。这说明引入哑变量,一定程度上提高了模型的预测精度,而且所建模型比较合理,形式相对简单,便于应用,不仅使得林分水平上的林分断面积、蓄积量的预测结果具有相容性,同时还考虑了间伐措施对林分生长的影响,达到了林分生长与收获模型整体化研究的目的,为林分的经营管理提供了可靠依据。
  • 表  1  样地基本情况

    Table  1.   The condition of plots

    样地类型 项目 年龄/a 平均胸径/cm 平均树高/m 林木株数/(株·hm-2) 林分断面积/(m2·hm-2) 林分蓄积/(m3·hm-2)
    最大值 56 20.0 16.5 2 054 32.925 9 12.701 6
    间伐样地 最小值 19 8.4 5.4 180 1.292 7 4.287 9
    平均值 38 12.9 9.5 992 13.394 4 52.360 3
    最大值 62 22.4 15.7 2 339 46.139 5 200.689 7
    未间伐样地 最小值 15 7.0 4.2 225 0.710 2 5.082 5
    平均值 35 12.2 8.9 955 12.480 7 49.293 6
    最大值 62 22.4 16.5 2 339 46.139 5 200.689 7
    全部样地 最小值 15 7.0 4.2 180 0.710 2 4.287 9
    平均值 36 12.6 9.2 973 12.927 2 50.792 3
    下载: 导出CSV

    表  2  建模及检验数据样地、样木分布情况

    Table  2.   Distribution of plots and tree of fit data and validation data

    项目 间伐林分 未间伐林分 合计
    样地/个 样木/株 样地/个 样木/株 样地/个 样木/株
    建模数据 41 2 876 39 2 504 80 5 380
    检验数据 24 1 423 29 1 826 53 3 249
    全部样地 65 4 299 68 4 330 133 8 629
    下载: 导出CSV

    表  3  联立方程组参数统计

    Table  3.   Parameter statistics of simultaneous equations

    项目 a0/b0 a1/b1 a2/b2 a3/b3 a4/b4 a5/b5
    传统方法联立方程组 林分断面积模型 0.924 0 0.681 0 0.011 2 0.630 6 1.365 9
    林分蓄积模型 5.678 6 0.246 0 0.001 0 3.421 9 0.280 8
    哑变量方法联立方程组 林分断面积模型 0.879 8 0.867 9 0.595 7 0.015 3 0.676 0 1.456 2
    林分蓄积模型 4.433 4 4.490 5 0.314 1 0.000 2 3.497 9 0.282 9
    下载: 导出CSV

    表  4  联立方程组评价指标统计

    Table  4.   Evaluation indicators statistics of simultaneous equations

    项目 R2 DMA ERMS P
    传统模塑联立方程组 林分断面积模型 0.849 6 0.160 0 0.204 5 0.945 7
    林分蓄积模型 0.827 2 0.737 1 0.284 4 0.942 0
    哑变量模型联立方程组 林分断面积模沏 0.852 2 0.142 0 0.200 6 0.956 5
    林分蓄积模型 0.831 2 0.706 5 0.280 1 0.959 9
    下载: 导出CSV

    表  5  误差变量方程组的检验

    Table  5.   Inspection of equations with error variables

    项目 R2 DMA ERMS P
    传统方法的联立方程组 林分断面积模型 0.900 1 0.205 2 0.332 0 0.921 5
    林分蓄积模型 0.912 3 0.379 3 0.688 0 0.928 3
    含哑变量的联立方程组 林分断面积模型 0.929 7 0.201 2 0.288 1 0.939 8
    林分蓄积模型 0.932 8 0.269 3 0.589 2 0.932 6
    下载: 导出CSV
  • [1] 赵俊卉. 长白山云冷杉混交林生长模型的研究[D]. 北京: 北京林业大学, 2010.

    ZHAO Junhui. Growth Modeling for Spruce-Fir Forest in Changbai Mountains [D]. Beijing:Beijing Forestry University, 2010.
    [2] 段爱国, 张建国, 孙洪刚, 等.林分断面积生长模拟理论与技术研究[J].世界林业研究, 2013, 26(2):43-47.

    DUAN Aiguo, ZHANG Jianguo, SUN Honggang, et al. Research progress of growth simulation theories and technologies for stand basal area [J]. World For Res, 2013, 26(2): 43-47.
    [3] 张雄清, 张建国, 段爱国.基于单木水平和林分水平的杉木兼容性林分蓄积量模型[J].林业科学, 2014, 50(1):82-87.

    ZHANG Xiongqing, ZHANG Jianguo, DUAN Aiguo. Compatibility of stand volume model for Chinese fir based on tree-level and stand-level [J]. Sci Silv Sin, 2014, 50(1): 82-87.
    [4] 邓静, 陈宇拓.利用增长量分配模型的杉木林分生长预测建模[J].浙江农林大学学报, 2014, 31(6):898-904.

    DENG Jing, CHEN Yutuo. Stand growth prediction based on a growth distribution model [J]. J Zhejiang A&F Univ, 2014, 31(6):898-904.
    [5] 戎建涛, 刘殿仁, 林召忠, 等.东北过伐林区主要森林类型林分蓄积量生长模型[J].林业科技开发, 2011, 25(1):30-34.

    RONG Jiantao, LIU Dianren, LIN Zhaozhong, et al. Volume growth models of main forest types in over-logged forest region, Northeast China [J]. China For Sci Technol, 2011, 25(1): 29-33.
    [6] 冉啟香, 邓华锋, 吕常笑, 等.油松林分断面积与蓄积量生长模型研究[J].西北林学院学报, 2016, 31(5):217-223.

    RAN Qixiang, DENG Huafeng, LÜ Changxiao, et al. Stand basal area and volume growth for Pinus tabulaeformis [J]. J Northwest For Univ, 2016, 31(5): 217-223.
    [7] 高东启, 邓华锋, 程志楚, 等.蒙古栎间伐林分和未间伐林分生长模型研究[J].中南林业科技大学学报, 2014, 34(2):50-54.

    GAO Dongqi, DENG Huafeng, CHENG Zhichu, et al. Study on growth models for thinned and un-thinned stands of Quercus mongolica [J]. J Cent South Univ For Technol, 2014, 34(2): 50-54.
    [8] 蒋益, 邓华锋, 高东启, 等.用度量误差模型方法建立油松树高曲线方程组[J].东北林业大学学报, 2015, 43(5):126-129.

    JIANG Yi, DENG Huafeng, GAO Dongqi, et al. Constructing height-diameter curve equations with measurement error models for Chinese pine stands [J]. J Northeast For Univ, 2015, 43(5): 126-129.
    [9] 李际平, 郭文清, 曹小玉.基于非线性度量误差的马尾松相容性立木生物量模型[J].中南林业科技大学学报, 2013, 33(6):22-25, 32.

    LI Jiping, GUO Wenqing, CAO Xiaoyu. Compatibility single-tree biomass model for Pinus massoniana stands based on nonlinear measurement error [J]. J Cent South Univ For Technol, 2013, 33(6): 22-25.
    [10] 陈振雄, 贺东北, 丁长春.利用度量误差模型方法建立海南省桉树、木麻黄、马占相思重量与材积相容性模型[J].中南林业调查规划, 2012, 31(4):5-9.

    CHEN Zhenxiong, HE Dongbei, DING Changchun. Using error-in-variable modeling method to establish compatible weight and volume equations system for Eucalyptus, Casuarina equisetifolia and Acacia mangium in Hainan Province [J]. Cent South For Inv Plann, 2012, 31(4): 5-9.
    [11] 高东启, 邓华锋, 程志楚, 等.基于度量误差模型方法建立的林分相容性树高曲线方程组[J].西北农林科技大学学报(自然科学版), 2015, 43(5):65-70.

    GAO Dongqi, DENG Huafeng, CHENG Zhichu, et al. Measurement error model based compatible system of equations for height-DBH curve [J]. J Northwest A&F Univ Nat Sci Ed, 2015, 43(5): 65-70.
    [12] 高云昌, 张文辉, 何景峰, 等.黄龙山油松人工林间伐效果的综合评价[J].应用生态学报, 2013, 24(5):1313-1319.

    GAO Yunchang, ZHANG Wenhui, HE Jingfeng, et al. Effects of thinning intensity on Pinus tabulaeformis plantation in Huanglong Mountain Northwest China: a comprehensive evaluation [J]. Chin J Appl Ecol, 2013, 24(5): 1313-1319.
    [13] 杜纪山, 唐守正.林分断面积生长模型研究评述[J].林业科学研究, 1997, 10(6):599-606.

    DU Jishan, TANG Shouzheng. The review of studies of stand basal area growth model [J]. For Res, 1997, 10(6): 599-606.
    [14] 杜纪山, 唐守正.杉木林分断面积生长预估模型及其应用[J].北京林业大学学报, 1998, 20(4):4-8.

    DU Jishan, TANG Shouzheng. Basal area growth prediction of Chinese fir stands and its application [J]. J Beijing For Univ, 1998, 20(4): 1-5.
    [15] 卢军, 张会儒, 李凤日.大兴安岭天然林林分生长模型研究[J].林业资源管理, 2011(3):33-36.

    LU Jun, ZHANG Huiru, LI Fengri. Natural forest stand growth models for Daxing'anling Mountains [J]. For Resour Manage, 2011(3): 33-37.
    [16] 杜纪山, 唐守正, 王洪良.天然林分生长模型在小班数据更新中的应用[J].林业科学, 2000, 36(3):52-58.

    DU Jishan, TANG Shouzheng, WANG Hongliang. Application of natural stand growth models in subcomparatment data update [J]. Sci Silv Sin, 2000, 36(3): 52-58.
    [17] 贺姗姗. 北京山区油松人工林林分结构与生长模拟研究[D]. 北京: 北京林业大学, 2009.

    HE Shanshan. Research on Structure and Growth of Pinus tabulaeformis Stands Simulation in Beijing [D]. Beijing:Beijing Forestry University, 2009.
    [18] 韩文娟, 袁晓青, 张文辉.油松人工林林窗对幼苗天然更新的影响[J].应用生态学报, 2012, 23(11):2940-2948.

    HAN Wenjuan, YUAN Xiaoqing, ZHANG Wenhui. Effects of gap size on seedling natural regeneration in artificial Pinus tabulaeformis plantation [J]. Chin J Appl Ecol, 2012, 23(11): 2940-2948.
    [19] 高东启. 北京市蒙古栎、油松林分生长预估模型研究[D]. 北京: 北京林业大学, 2014.

    GAO Dongqi. The Research of Stand Growth Forecast Models for Quercus mongolica and Pinus tabulaeformis in Beijing [D]. Beijing:Beijing Forestry University, 2014.
    [20] 高东启, 邓华锋, 蒋益, 等.油松林分断面积生长预估模型研究[J].西南林业大学学报, 2015, 35(1):42-46.

    GAO Dongqi, DENG Huafeng, JIANG Yi, et al. Forecast models research of stands basal area growth for Pinus tabulaeformis [J]. J Southwest For Univ, 2015, 35(1): 42-46.
    [21] 李永慈, 唐守正.带度量误差的全林整体模型参数估计研究[J].北京林业大学学报, 2006, 28(1):23-27.

    LI Yongci, TANG Shouzheng. Parameter estimate of the whole stand model with measurement error [J]. J Beijing For Univ, 2006, 28(1): 23-27.
  • [1] 窦啸文, 吴登瑜, 张笑菁, 汤孟平.  天目山常绿阔叶林胸高断面积生长量影响因子研究 . 浙江农林大学学报, 2023, 40(5): 1063-1072. doi: 10.11833/j.issn.2095-0756.20220651
    [2] 谭灿灿, 聂稳, 刘逸夫, 王亚, 董遥, 黄睿智, 刘建锋, 史胜青, 常二梅, 赵秀莲, 江泽平, 贾子瑞.  3种栎属树种单木生长模型及树皮厚度模型研究 . 浙江农林大学学报, 2023, 40(3): 589-597. doi: 10.11833/j.issn.2095-0756.20220455
    [3] 申家朋, 陈东升, 孙晓梅, 张守攻.  基于似乎不相关回归和哑变量的日本落叶松单木生物量模型构建 . 浙江农林大学学报, 2019, 36(5): 877-885. doi: 10.11833/j.issn.2095-0756.2019.05.005
    [4] 雷泽勇, 韩艳刚, 赵国军, 周晏平, 张岩松, 于德良.  辽宁章古台樟子松生长过程分析 . 浙江农林大学学报, 2018, 35(2): 324-330. doi: 10.11833/j.issn.2095-0756.2018.02.017
    [5] 王海宾, 彭道黎, 高秀会, 李文芳.  基于GF-1 PMS影像和k-NN方法的延庆区森林蓄积量估测 . 浙江农林大学学报, 2018, 35(6): 1070-1078. doi: 10.11833/j.issn.2095-0756.2018.06.010
    [6] 刘崴, 魏天兴, 朱清科.  半干旱黄土丘陵区河北杨和油松生长季树干液流特征 . 浙江农林大学学报, 2018, 35(6): 1045-1053. doi: 10.11833/j.issn.2095-0756.2018.06.007
    [7] 罗恒春, 张超, 魏安超, 张一, 黄田, 余哲修.  云南松林分平均胸径生长模型及模型参数环境解释 . 浙江农林大学学报, 2018, 35(6): 1079-1087. doi: 10.11833/j.issn.2095-0756.2018.06.011
    [8] 刘薇祎, 邓华锋, 冉啟香, 黄国胜, 王雪军.  湖南省杉木林分相容性树高曲线方程组研究 . 浙江农林大学学报, 2017, 34(6): 1051-1058. doi: 10.11833/j.issn.2095-0756.2017.06.012
    [9] 许杰, 戚大伟.  基于特征点提取的林木生长状态无损测量方法 . 浙江农林大学学报, 2016, 33(3): 403-408. doi: 10.11833/j.issn.2095-0756.2016.03.005
    [10] 冉啟香, 邓华锋, 黄国胜, 王雪军, 陈振雄.  云南松地上生物量模型研究 . 浙江农林大学学报, 2016, 33(4): 605-611. doi: 10.11833/j.issn.2095-0756.2016.04.008
    [11] 吕常笑, 邓华锋, 王少杰, 陈振雄, 王雪军.  马尾松不同区域相容性立木材积和地上生物量模型 . 浙江农林大学学报, 2016, 33(5): 790-797. doi: 10.11833/j.issn.2095-0756.2016.05.010
    [12] 焦祥, 郑加强, 张慧春, 苏朦朦.  林木虚拟生长建模方法及建模工具研究综述和展望 . 浙江农林大学学报, 2015, 32(6): 966-975. doi: 10.11833/j.issn.2095-0756.2015.06.021
    [13] 王晓宁, 徐天蜀, 李毅.  利用ALOS PALSAR双极化数据估测山区森林蓄积量模型 . 浙江农林大学学报, 2012, 29(5): 667-670. doi: 10.11833/j.issn.2095-0756.2012.05.005
    [14] 徐文兵, 高飞, 杜华强.  几种测量方法在森林资源调查中的应用与精度分析 . 浙江农林大学学报, 2009, 26(1): 132-136.
    [15] 马丰丰, 贾黎明.  北京地区侧柏、油松带皮胸径与去皮胸径的关系 . 浙江农林大学学报, 2009, 26(1): 13-16.
    [16] 刘永红, 杨培华, 韩创举, 樊军锋, 李新会, 李安平, 杨世荣.  油松不同种源种实性状的变异分析 . 浙江农林大学学报, 2008, 25(2): 163-168.
    [17] 景丽, 朱志红, 王孝安, 郭华.  秦岭油松人工林与次生林群落特征比较 . 浙江农林大学学报, 2008, 25(6): 711-717.
    [18] 刘永红, 樊军锋, 杨培华, 韩创举.  油松单亲子代苗期生长性状遗传分析 . 浙江农林大学学报, 2005, 22(5): 513-517.
    [19] 周永学, 樊军锋, 杨培华, 高建社, 刘永红.  奥地利黑松与油松1 年生苗生长和生物量对比分析 . 浙江农林大学学报, 2003, 20(4): 438-441.
    [20] 周国模, 郭仁鉴, 韦新良, 王雪军.  浙江省杉木人工林生长模型及主伐年龄的确定 . 浙江农林大学学报, 2001, 18(3): 219-222.
  • 加载中
  • 链接本文:

    https://zlxb.zafu.edu.cn/article/doi/10.11833/j.issn.2095-0756.2018.01.009

    https://zlxb.zafu.edu.cn/article/zjnldxxb/2018/1/68

计量
  • 文章访问数:  3246
  • HTML全文浏览量:  491
  • PDF下载量:  367
  • 被引次数: 0
出版历程
  • 收稿日期:  2016-12-19
  • 修回日期:  2017-03-15
  • 刊出日期:  2018-02-20

基于度量误差方法的油松林分生长模型

doi: 10.11833/j.issn.2095-0756.2018.01.009
    基金项目:

    国家林业公益性行业科研专项 201204510

    作者简介:

    王金池, 从事森林资源监测与评价研究.E-mail:2276599635@qq.com

    通信作者: 邓华锋, 教授, 博士, 从事森林可持续经营理论与技术研究.E-mail:denghuafeng@bjfu.edu.cn
  • 中图分类号: S758.5

摘要: 为掌握北京地区油松Pinus tabulaeformis生长过程,建立相容的断面积和蓄积模型,通过引入间伐林分与未间伐林分的哑变量,分别建立林分断面积、林分蓄积生长模型,然后从相容的角度出发,建立林分断面积、蓄积量的误差变量联立方程组,并与不含哑变量的传统生长模型的误差变量联立方程组进行比较。经检验,传统误差变量联立方程组中林分断面积和林分蓄积生长模型的预测精度都在92%以上,对油松林分断面积的预测精度高达0.921 5,决定系数高达0.900 1,对油松林分蓄积量的预测精度达到了0.928 3,决定系数高达0.912 3,而引入哑变量的误差变量联立方程组中,模型的预测精度和确定系数稍高,均在93%以上,对油松林分断面积的预测精度高达0.939 8,决定系数达到了0.927 9,对油松林分蓄积量的预测精度在0.930 0以上,决定系数达0.932 8。这说明引入哑变量,一定程度上提高了模型的预测精度,而且所建模型比较合理,形式相对简单,便于应用,不仅使得林分水平上的林分断面积、蓄积量的预测结果具有相容性,同时还考虑了间伐措施对林分生长的影响,达到了林分生长与收获模型整体化研究的目的,为林分的经营管理提供了可靠依据。

English Abstract

王金池, 冉啟香, 邓华锋, 黄国胜, 王雪军. 基于度量误差方法的油松林分生长模型[J]. 浙江农林大学学报, 2018, 35(1): 68-74. doi: 10.11833/j.issn.2095-0756.2018.01.009
引用本文: 王金池, 冉啟香, 邓华锋, 黄国胜, 王雪军. 基于度量误差方法的油松林分生长模型[J]. 浙江农林大学学报, 2018, 35(1): 68-74. doi: 10.11833/j.issn.2095-0756.2018.01.009
WANG Jinchi, RAN Qixiang, DENG Huafeng, HUANG Guosheng, WANG Xuejun. Stand growth model for Pinus tabulaeformis using measurement error method[J]. Journal of Zhejiang A&F University, 2018, 35(1): 68-74. doi: 10.11833/j.issn.2095-0756.2018.01.009
Citation: WANG Jinchi, RAN Qixiang, DENG Huafeng, HUANG Guosheng, WANG Xuejun. Stand growth model for Pinus tabulaeformis using measurement error method[J]. Journal of Zhejiang A&F University, 2018, 35(1): 68-74. doi: 10.11833/j.issn.2095-0756.2018.01.009
  • 研究林分生长与收获预估模型可以了解森林的生长变化规律和预估林分生长量、收获量[1]。在林分生长与收获预估模型的研究中,林分断面积是林分生长收获模型系统的核心部分[2],而林分蓄积量的大小又标志着林地生产力的高低及经营措施的效果[3]。在目前的研究中,许多学者都是分别建立林分断面积与蓄积生长模型的[2-6],但理论上,林分断面积与林分蓄积量的大小密切相关,是一个常用的林分密度指标,有学者认为选用林分断面积作为密度指标来建立林分蓄积模型时,模型的预估效果会更好[7]。在林分生长与收获模拟预测体系中林分断面积既可以作为重要的自变量来预估蓄积,又可以作为因变量被模拟预测。通常我们认为只有因变量的观测值含有误差,但实际上,无论是自变量还是因变量都会存在来自抽样或测量等方面的误差。我们把这些误差都称为度量误差[8]。当自变量和因变量都有误差时,传统模型方法就不再适用于模型的拟合,因此需要采用度量误差模型来建立误差联立方程组,使模型之间具有相容性。最近几年,度量误差模型在林业上的应用研究越来越深入。李际平等[9]利用非线性度量误差模型建立了南方马尾松Pinus massoniana地上生物量与树干、树冠的相容性联立方程组模型;陈振雄等[10]利用度量误差模型方法建立了海南省桉树Eucalyptus spp.,木麻黄Casuarina equisetifolia,马占相思Acacia mangium的树干去皮材积与质量相容性联立方程组;高东启等[11]以Richards方程为基础模型,利采用度量误差模型对蒙古栎Quercus mongolica建立了林分相容性树高曲线方程组;蒋益等[8]也用该方法建立了油松Pinus tabulaeformis林分相容性树高曲线方程组。但目前,利用度量误差模型方法建立林分蓄积、断面积的误差联立方程组的研究尚鲜有报道。油松根系发达,有较强的适应性和抗逆性,是中国北方温带针叶林中分布最广的群落[12],因此,本研究在前人研究的基础上,利用北京市油松一类清查数据,在保证各因子之间既相容又能预测精度的前提下,用传统方法和引入以间伐措施为特征的哑变量方法为基础模型来建立相容性的林分断面积、蓄积联立方程组,使断面积、蓄积2个因子预测时保持相容性和一致性,达到在优化模型结构的同时减少外业调查误差对预测结果影响的目的,使模型更具生物学意义和解释性。

    • 本研究所用数据为国家林业局调查规划设计院1996,2001,2006年的北京市油松一类清查数据,各个样地面积为0.066 7 hm2。记录林木胸径、林分平均高、林分年龄、林分蓄积、采伐蓄积、郁闭度、海拔高度等林分调查因子。调查期内对部分林分进行了抚育间伐。剔除数据缺失、记录不详及明显有误的样地及疏林地样地,最终选出133块油松人工林样地用以建模和检验,包括65块间伐样地(共4 299株样木)和68块未间伐样地(共4 330株样木)。对133块样地进行随机抽样,其中80块样地用于建模,剩余53块样地用于检验。统计样地基本情况及随机抽样的样地、样木分布情况分别如表 1表 2所示。

      表 1  样地基本情况

      Table 1.  The condition of plots

      样地类型 项目 年龄/a 平均胸径/cm 平均树高/m 林木株数/(株·hm-2) 林分断面积/(m2·hm-2) 林分蓄积/(m3·hm-2)
      最大值 56 20.0 16.5 2 054 32.925 9 12.701 6
      间伐样地 最小值 19 8.4 5.4 180 1.292 7 4.287 9
      平均值 38 12.9 9.5 992 13.394 4 52.360 3
      最大值 62 22.4 15.7 2 339 46.139 5 200.689 7
      未间伐样地 最小值 15 7.0 4.2 225 0.710 2 5.082 5
      平均值 35 12.2 8.9 955 12.480 7 49.293 6
      最大值 62 22.4 16.5 2 339 46.139 5 200.689 7
      全部样地 最小值 15 7.0 4.2 180 0.710 2 4.287 9
      平均值 36 12.6 9.2 973 12.927 2 50.792 3

      表 2  建模及检验数据样地、样木分布情况

      Table 2.  Distribution of plots and tree of fit data and validation data

      项目 间伐林分 未间伐林分 合计
      样地/个 样木/株 样地/个 样木/株 样地/个 样木/株
      建模数据 41 2 876 39 2 504 80 5 380
      检验数据 24 1 423 29 1 826 53 3 249
      全部样地 65 4 299 68 4 330 133 8 629
    • 近年来,对林分断面积、蓄积生长模型的研究越来越多,但其中以Richards和Schumacher模型2种应用较为广泛[13-14]。本研究在初步的试验以及前人研究的基础上发现Richard方程能够较好地拟合油松的蓄积、断面积生长模型。利用地位指数来反映林分的立地质量[15-16],选取林分密度指数和林分断面积作为密度指标来拟合林分断面积和蓄积生长模型,采用Richard方程对油松的蓄积和断面积生长模型进行拟合。公式如下:

      $$ G = {a_1}I_{\rm{S}}^{{a_2}}{\{ 1- \exp [-{a_3}{({I_{{\rm{SD}}}}/1\;000)^{{a_4}}}A]\} ^{{a_5}}}; $$ (1)
      $$ M = {b_1}I_{\rm{S}}^{{b_2}}{[1-\exp (-{b_3}{G^{{b_4}}}A)]^{{b_5}}}。 $$ (2)

      式(1)和式(2)中:GM分别为林分的断面积和蓄积,IS为地位指数[17]ISD为林分密度指数[17](根据公式ISD=N(Dg/D0)b获得,其中N为单位面积株数,Dg为林分平均胸径,D0为20 cm,b为自稀疏率1.384),A为林分年龄,a1~a5b1~b5为模型的预估参数。

    • 间伐是最重要的营林技术措施之一[18]。现实中的林分由于抚育间伐的影响,间伐与未间伐林分在生态环境、林分结构和生长规律上可能存在差异。一般来说,需要分开建模,但间伐和未间伐林分之间又存在一些联系[19],因此如何把不同的林分结合起来建立统一的模型并在保证模型预测值可靠性下使模型简化,是值得研究的问题。哑变量模型方法为解决这一问题提供了可能途径[20]。本研究在前人研究的基础上,在模型中引入哑变量,用定性代码来表示间伐和未间伐林分,以整合构建具有相容性的统一模型。具体过程是将第i个类型的林分编号为Si,将定性数据Si转化为0和1,即:

      $$ {S_i} = \left\{ \begin{array}{l} 1, \;当{S_i}为第i个类型的林分时\\ 0, \;否则 \end{array} \right.。 $$

      其中:i= 1,2,S1S2分别为间伐林分和未间伐林分的定性代码。其形式如下:

      $$ G = ({a_0}{S_1} + {a_1}{S_2})I_{\rm{S}}^{{a_2}}\{ 1- \exp [-{a_3}{({I_{{\rm{SD}}}}{\rm{/}}1\;000)^{{a_4}}}A)]{\} ^{{a_5}}}; $$ (3)
      $$ M = ({b_0}{S_1} + {b_1}{S_2})I_{\rm{S}}^{{b_2}}{[1-\exp (-{b_3}{G^{{b_4}}}A)]^{{b_5}}}。 $$ (4)

      式(3)和式(4)中:GM分别为林分的断面积和蓄积,IS为地位指数,ISD为林分密度指数,A为林分年龄,S1S2为间伐林分和未间伐林分的定性代码,a0~a5b0~b5为模型的预估参数。

    • 野外调查中,无论是自变量还是因变量都存在一定误差,且误差来源于许多方面,如抽样误差、测量误差等。这些随机误差都被称为度量误差[8]。采用度量误差模型来建立误差联立方程组,使模型之间更具有相容性。公式如下[21]

      $$ \left\{ \begin{array}{l} f({y_i}, {x_i}, c) = 0\\ {Y_i} = {y_i} + {e_i}, i = 1, 2, \cdots, n\\ E({e_i}) = 0, {\rm{cov}}({e_i}) = {\sigma ^2}\psi \end{array} \right.。 $$ (5)

      式(5)中:xiq维无误差变量的观测数据,yip维误差变量的观测数据,fm维向量函数,Yiyi的未知真值,误差的协方差矩阵记为Φ=σ2ψψei的误差结构矩阵,σ2为估计误差。

      以传统模型和哑变量模型为基础,利用非线性度量误差方法分别联立式(1),式(2)和式(3),式(4)构成联立方程组,以保证MG方程之间的相容性和一致性。基于传统方法的联立方程组形式如下:

      $$ \left\{ \begin{array}{l} G = {a_1}I_{\rm{S}}^{{a_2}}\{ 1- \exp [-{a_3}{({I_{{\rm{SD}}}}{\rm{/}}1\;000)^{{a_4}}}A)]{\} ^{{a_5}}}\\ M = {b_1}I_{\rm{S}}^{{b_2}}{[1-\exp (-{b_3}{G^{{b_4}}}A)]^{{b_5}}} \end{array} \right.。 $$ (6)

      基于哑变量方法的联立方程组如下:

      $$ \left\{ \begin{array}{l} G = ({a_0}{S_1} + {a_1}{S_2})I_{\rm{S}}^{{a_2}}\{ 1- \exp [-{a_3}{({I_{{\rm{SD}}}}{\rm{/}}1\;000)^{{a_4}}}A)]{\} ^{{a_5}}}\\ M = ({b_0}{S_1} + {b_1}{S_2})I_{\rm{S}}^{{b_2}}{[1-\exp (-{b_3}{G^{{b_4}}}A)]^{{b_5}}} \end{array} \right.。 $$ (7)
    • 利用Forstat 2.0软件求解参数,采用Excel,SPSS进行数据计算和模型检验。对所建传统模型和哑变量模型进行t检验,并通过平均绝对偏差(DMA),均方根误差(ERMS),决定系数(R2)和预估精度(P)等几个指标对模型进行检验和比较。计算表达式如下:

      $$ {D_{{\rm{MA}}}} = \sum\limits_{i = 1}^n {\left| {\frac{{{y_i}-{{\hat y}_i}}}{n}} \right|} ; $$ (8)
      $$ {E_{{\rm{RMS}}}} = \sqrt {\frac{{\sum\limits_{i = 1}^n {{{({y_i}-{{\hat y}_i})}^2}} }}{n}} ; $$ (9)
      $$ {R^2} = 1 - \frac{{\sum\limits_{i = 1}^n {{{({y_i} - {{\hat y}_i})}^2}} }}{{\sum\limits_{i = 1}^n {{{({y_i} - {{{\hat{\bar{y}}}}_i})}^2}} }}; $$ (10)
      $$ P = \left( {1 - \frac{{{t_{0.05}} \cdot \sqrt {\frac{{\sum\limits_{i = 1}^n {{{({y_i} - {{\hat y}_i})}^2}} }}{{n(n - p)}}} }}{{{\hat{\bar{y}}}}}} \right) \times 100\% 。 $$ (11)

      式(8)~式(11)中:yi为断面积、蓄积的实测值,${{\hat y}_i}$为模型预估值,${{\hat{\bar{y}}}}$为模型预估值的平均值,n为样本数,p为模型参数个数,t0.05为置信水平为0.05时的t分布值。

    • 通过ForStat软件中的非线性度量误差联立方程组求解林分蓄积、断面积生长模型的参数如表 3所示,统计各评价指标结果如表 4所示。

      表 3  联立方程组参数统计

      Table 3.  Parameter statistics of simultaneous equations

      项目 a0/b0 a1/b1 a2/b2 a3/b3 a4/b4 a5/b5
      传统方法联立方程组 林分断面积模型 0.924 0 0.681 0 0.011 2 0.630 6 1.365 9
      林分蓄积模型 5.678 6 0.246 0 0.001 0 3.421 9 0.280 8
      哑变量方法联立方程组 林分断面积模型 0.879 8 0.867 9 0.595 7 0.015 3 0.676 0 1.456 2
      林分蓄积模型 4.433 4 4.490 5 0.314 1 0.000 2 3.497 9 0.282 9

      表 4  联立方程组评价指标统计

      Table 4.  Evaluation indicators statistics of simultaneous equations

      项目 R2 DMA ERMS P
      传统模塑联立方程组 林分断面积模型 0.849 6 0.160 0 0.204 5 0.945 7
      林分蓄积模型 0.827 2 0.737 1 0.284 4 0.942 0
      哑变量模型联立方程组 林分断面积模沏 0.852 2 0.142 0 0.200 6 0.956 5
      林分蓄积模型 0.831 2 0.706 5 0.280 1 0.959 9

      表 4可以看出:在传统方法的联立方程组中,林分断面积模型的决定系数为0.849 6,林分蓄积模型的决定系数为0.827 2,两者的预估精度都在0.940 0以上,说明2个模型的拟合效果均较好;而在含哑变量模型的联立方程组中,林分断面积生长模型的决定系数为0.852 2,预估精度为0.956 5,林分蓄积生长模型的决定系数为0.831 2,预估精度为0.959 9,相比较于传统方法而言,含哑变量方法的林分断面积、蓄积生长模型的确定系数和预估精度都有所提高,而且各项误差也相应减少了。这在一定程度上说明含哑变量方法的模型拟合效果比传统模型好。

      利用检验数据分别对2种方法的断面积、蓄积联立方程组进行检验,检验结果如表 5所示。从表 5可以看出:传统方法与哑变量方法的预测精度都较高,超过了0.900 0,就林分断面积而言,传统方法的林分断面积模型的平均绝对偏差为0.205 2,均方根误差为0.332 0,决定系数为0.900 1,预测精度为0.921 5;含哑变量方法的林分断面积模型的平均绝对偏差为0.201 2,均方根误差为0.288 1,决定系数为0.929 7,预测精度为0.939 8;对于林分蓄积,传统方法的林分蓄积模型的平均绝对偏差为0.379 3,均方根误差为0.688 0,决定系数为0.912 3,预测精度为0.928 3;含哑变量方法的林分蓄积模型的平均绝对偏差为0.269 3,均方根误差为0.589 2,决定系数为0.932 8,预测精度为0.932 6。以上数据表明各项误差较小,说明所建的林分断面积、蓄积量联立方程组是比较合理的。

      表 5  误差变量方程组的检验

      Table 5.  Inspection of equations with error variables

      项目 R2 DMA ERMS P
      传统方法的联立方程组 林分断面积模型 0.900 1 0.205 2 0.332 0 0.921 5
      林分蓄积模型 0.912 3 0.379 3 0.688 0 0.928 3
      含哑变量的联立方程组 林分断面积模型 0.929 7 0.201 2 0.288 1 0.939 8
      林分蓄积模型 0.932 8 0.269 3 0.589 2 0.932 6

      对模型进行成对t检验,在0.05的显著水平上回归效果显著,实测值与预测值之间无显著差异,说明林分断面积、蓄积量的联立方程组预测效果较好,采用度量误差方法建立的林分断面积、蓄积量误差联立方程组是比较合理的。进一步分析传统模型及哑变量模型联立方程组的断面积、蓄积残差可知:绝大多数残差点都落在各自“残差平均值± 2倍标准差”之间,残差分布无明显异质性,说明模型的拟合效果较好,而且在一定程度上也说明引入哑变量,能适当提高模型的精度并能解决模型单独建模模型不相容的问题。

    • 为使林分断面积、蓄积量的预测结果保持一致,预估精度更高,本研究从模型相容的角度出发,以北京地区油松一类清查数据为例,采用度量误差方法分别建立了基于传统方法的林分断面积与林分蓄积联立方程组和基于哑变量方法的林分断面积、蓄积联立方程组,并将两者进行对比。经检验,2种模型方法的林分断面积、蓄积模型的预测精度都较高,超过了90%。在林分断面积联立方程组中,含哑变量方法的林分断面积模型的决定系数为0.929 7,预测精度为0.939 8;传统方法的林分断面积模型的决定系数为0.900 1,预测精度为0.921 5。在林分蓄积量联立方程组中,传统方法的林分蓄积模型的决定系数为0.912 3,预测精度为0.928 3;含哑变量方法的林分蓄积模型的决定系数为0.932 8,预测精度为0.932 6。

      从林分水平生长模型相容的角度出发,引入间伐林分与未间伐林分的哑变量建立了林分断面积与蓄积的误差变量联立方程组。结果表明:含哑变量方法的联立方程组的预测精度和确定系数较传统方法都稍高,各项预测误差均稍小,说明模型较合理,使得林分断面积、蓄积量模型具有相容性和一致性。所建模型不仅考虑了抚育间伐措施对林分断面积、蓄积量生长的影响,同时还考虑了林分断面积和林分蓄积之间的度量误差,优化了模型结构,减小了外业调查中的误差对预测结果的影响,为从林分水平预测油松林分断面积和蓄积生长模型提供了参考和依据。

      所建立的林分断面积与蓄积的误差变量联立方程组形式相对简单,便于在林业实践中具体应用,可以用来预测北京市油松人工林断面积、蓄积量的生长规律。但是,本研究所使用的数据来源于3期不同时期的数据,在建模过程中并未考虑时间效应的影响,而且间伐与未间伐林分的实验数据有限。在后续研究中,可收集更多实验数据,同时考虑区域、时间和其他经营管理措施等因素,采用混合模型的方法建模,使得建模结构更加合理,从而为油松林分的经营管理提供参考和依据。

参考文献 (21)

目录

    /

    返回文章
    返回