• 中国中文核心期刊
  • 中国科学引文数据库(CSCD)核心库来源期刊
  • 中国科技论文统计源期刊(CJCR)
  • 第二届国家期刊奖提名奖

留言板

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

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

基于试验数据的人工用材林最优轮伐期确定的研究进展

刘林 孙洪刚 王宇华 张建国

引用本文:
Citation:

基于试验数据的人工用材林最优轮伐期确定的研究进展

    作者简介: 刘林,中国计量大学经济与管理学院副教授,博士。研究方向:林业生产效率。电话:13634119945 Email:jhliulin@cjlu.edu.cn.
    通讯作者: 孙洪刚, honggangsun@caf.ac.cn
  • 中图分类号: S752

Research Progress in the Optimal Rotation of Forest Plantations Based on Experimental data

    Corresponding author: SUN Hong-gang, honggangsun@caf.ac.cn ;
  • CLC number: S752

  • 摘要: 虽然用材林轮伐期研究长期受到国内外众多学者的关注,但有关轮伐期的内涵、轮伐期确定的标准、轮伐期模型选择以及相关变量对轮伐期的影响等方面,依然是众说纷纭,无法为构建用材林提质增效培育模式提供坚实的理论基础。本文通过梳理国内外有关轮伐期确定方面的研究成果,并以我国杉木用材林试验数据作为支撑,意在展示不同轮伐期确定准则、方法和相关影响因素对轮伐期确定的作用机理和过程,并对未来轮伐期的研究方向进行了展望。
  • 表 1  轮伐期确定的6个准则及确定条件

    Table 1.  Conditions for determining rotation under 6 criteria

    准则
    Criteria
    目标函数
    Objective function
    轮伐期确定条件
    Conditions for determining rotation
    文献
    References
    MGY max(T)
    {pQ(T)}
    pQT = 0 = QT [6]
    MSY max(T)
    {pQ(T)/T}
    pQT/pQ(T) = 1/T [7-9]
    PNW max(T)
    {pQ(T)e−rTwE}
    QT/Q(T) = r [10-11]
    LEV max(T)
    {[pQ(T)e−rTwE](1- e−rT )}
    pQT = rpQ(T) + rLEV [12-15]
    MANR max(T)
    {(pQ(T) - wE)/T}
    l/T = pQ(T)/[pQ(T)-wE] [16]
    IRR max(i)
    {[pQ(T)e−iTwE]/(1-e−iT)} = 0
    pQT = i × pQ(T) [1-3,17-20]
    注:T为林分年龄;p是活立木价格;r为折现率;i为内部收益率;wE为林地更新成本;Q(T)为林分生长量,是林分年龄的凹函数;QT为林分生长量对T的一阶导[4]
      Notes: T represents the stand age; P is the price of stumpage; r is the discount rate; i is the internal rate of return; wE is the cost of forest land renewal; Q(T) is the stand growth, which is a concave function of stand age, QT is the first derivative of Q(T) on T
    下载: 导出CSV

    表 2  6个最优化条件下求解杉木用材林最优轮伐期

    Table 2.  Optimal rotations for Chinese fir plantation under 6 criteria

    林分年龄
    Age
    Q(T)QTQ(T)/TQT/Q(T)QT/[Q(T)−wE/p]i*1/Tr/(1−e−rT)
    749.47624.9677.0680.5051.5470.0090.1430.169
    975.27824.3428.3640.3230.5800.0810.1110.138
    1199.59122.2459.0540.2230.3360.1160.0910.118
    13121.52419.7919.3480.1630.224
    0.1240.0770.105
    15140.95917.4329.3970.1240.1620.1200.0670.095
    17158.08015.3219.2990.0970.1230.1180.0590.087
    19173.15713.4899.1140.078
    0.096
    0.1120.0530.081
    21186.47011.9198.8800.0640.0790.1050.0480.076
    23198.27110.5798.6200.0530.0640.1010.0430.073
    25208.7799.4348.3510.0450.0520.0940.0400.070
    27218.1798.453
    8.0810.0390.042
    0.0910.0370.067
    29226.6287.3097.8150.0330.0350.0890.0340.065

    45326.1030.0020.0000.000
    下载: 导出CSV

    表 3  常用生长方程的数量成熟龄确定条件

    Table 3.  Condition for quantitative mature age under difference growth functions

    序号
    No.
    生长方程
    Growth function
    方程形式
    Equations
    数量成熟龄确定条件
    Condition for quantitative rotation
    1 Compertz模型 $ Q=kexp(-{e}^{a-bt}) $ $ {e}^{(a-bT)}={\left(bT\right)}^{-1} $
    2 Logistic模型 $ Q=\dfrac{k}{1 + {e}^{(a-bt)}} $ $ {e}^{(-a + bT)}=bT-1 $
    3 Mitschelich模型 $ Q=k(1-a{\mathrm{e}}^{-bt}) $ $ {e}^{bT}=abT + a $
    4 Richards模型 $ Q=A(1-b{e}^{-kt}{)}^{1/1-m} $ $ {e}^{bT}-b=\dfrac{kbT}{1-m} $
    5 Weibull模型 $ Q=A\left\{1-\mathrm{e}\mathrm{x}\mathrm{p}\left[-{\left(\dfrac{t}{b}\right)}^{c}\right]\right\} $ ${e}^{{\left(\frac{T}{b}\right)}^{c} }=1 + \dfrac{c}{b}{\left(\dfrac{T}{b}\right)}^{c-1}\cdot T$
    注:t为林龄;T 是数量成熟龄;kQ的上渐近值;a为与Q初植有关的参数;b为内禀生长率;c为形状参数;m为曲线形状及拐点位置参数
      Notes: t represents the stand age; T is the quantitative rotation; k is the upper asymptotic value; a is the parameter related with the initial Q value; b is the intrinsic growth rate; c is the shape parameter; m is the curve shape inflection point position parameter
    下载: 导出CSV

    表 4  杉木人工林不同密度下最优轮伐期

    Table 4.  Optimum rotation of Chinese fir plantation under different density

    初植密度
    Initial planting density
    (Number·ha−1)
    生长方程形式
    Functions of growth
    CompetzLogisticMitschelichRichardsWeibull
    16661617181717
    33331819191918
    50002523282425
    66663829342726
    100003829322729
    下载: 导出CSV

    表 5  对Faustmann模型进行扩展的研究

    Table 5.  Expansions for Faustmann model

    扩展方向
    Direction
    主要研究结论
    Main conclusions
    文献
    References
    考虑初植密度和间伐 初植密度对轮伐期有重要的影响,初植密度相同时,没有间伐的林分轮伐期比间伐林分的轮伐期短。 [25-33]
    考虑活立木价格变化 立木价格增加,轮伐期相应缩短。立木质量越好(如,径阶越大),价格越高,如果此时提高立木的价格可能会延长轮伐期。 [34-39]
    加入立地条件 立地条件对轮伐期的影响表现出非单调性。立地条件较差的林分,轮伐期对盈利影响不大,但是立地条件好的林分,轮伐期的确定对盈利有很大的影响。 [40-42]
    加入非木材收益 加入非木材收益后倾向于延长轮伐期(如,考虑碳存储价值)。 [43-47]
    考虑折现率的变化 当折现率较高时,林地价值对轮伐期的变化更敏感。折现率较低时,林地价值不稳定。 [8,29,33,47-50]
    加入不确定性因素 森林价值对多种风险和不确定性很敏感。木材价格的变化比其它风险(如,风灾、病虫害等)造成的经济损失更大。但价格的不确定性,并不意味着对轮伐期有显著影响。 [51-56]
    下载: 导出CSV

    表 6  不同初植密度和折现率时杉木人工林的最优经济轮伐期

    Table 6.  Optimal rotations for Chinese fir plantation using different initial planting densities and discount rates

    初植密度
    Initial planting density (Number·ha−1)
    r = 0.01r = 0.02r = 0.03r = 0.04r = 0.05r = 0.06r = 0.07r = 0.08r = 0.09
    166720181717161513118
    333319161515131211109
    5000181515141313111110
    6667191817171614131313
    10000201918181717161515
    下载: 导出CSV
  • [1] 盛炜彤, 惠刚盈, 罗云伍. 大岗山杉木人工林主伐年龄的研究[J]. 林业科学研究, 1991, 2:7-3.

    [2] 张志云, 蔡学林. 杉木人工林主伐年龄研究[J]. 江西农业大学学报, 1991, 13(4):336-343.

    [3] 周国模, 郭仁鉴, 韦新良, 等. 浙江省杉木人工林生长模型及主伐年龄的确定[J]. 浙江农林大学学报, 2001, 18(3):219-222.

    [4]

    NEWMAN D H, GILBERT C B, HYDE W F. The optimal forest rotation with evolving prices[J]. Land Economics, 1985, 6(1): 347-353.
    [5]

    OHLIN B. Concerning the question of the rotation period in forestry[J]. Journal of Forest Economics, 1995, 1: 89-114.
    [6]

    THOMSON R B. An examination of basic principles of forest valuation[M]. Bu11.6. Durham, NC: Duke University School of Forestry. 1942, 99 pp.
    [7]

    OLSCHEWSKI R, BENITEZ P C. Optimizing joint production of timber and carbon sequestration of afforestation projects[J]. Journal of Forest Economics, 2010, 16(1): 1-10. doi: 10.1016/j.jfe.2009.03.002
    [8]

    BETTINGER P, BOSTON K, SIREY J P, et al. Forest Management and Planning[M]. Amsterdam: Academic Press, 2009.
    [9]

    LUCKERT M K, WILLIAMSON T. Should sustained yield be part of sustainable forest management[J]. Canadian Journal of Forest Research, 2005, 35(2): 356-364. doi: 10.1139/x04-172
    [10]

    PETIT B, MONTAGNINI F. Growth equations and rotation ages of ten native tree species in mixed and pure plantations in the humid neotropics[J]. Forest Ecology and Management, 2004, 199: 243-257. doi: 10.1016/j.foreco.2004.05.039
    [11]

    GAFFNEY M M. Concepts of Financial Maturity of Timber and Other Assets[M]. NC. State College, Dept. Agric. Econ. , A. E. Information Series No. 62.1960.
    [12]

    ZHANG D. Faustmann in an uncertain policy environment[J]. Forest Policy and Economics, 2001, 2(2): 203-210. doi: 10.1016/S1389-9341(01)00052-1
    [13]

    CHANG S J. Forest valuation under the generalized Faustmann formula[J]. Canadian Journal of Forest Research, 2014, 44(1): 56-63.
    [14]

    TAHVONEN O, SALO S, KUULUVAINEN J. Optimal forest rotation and land values under a borrowing constraint[J]. Journal of Economic Dynamics and Control, 2001, 25: 1595-1627. doi: 10.1016/S0165-1889(99)00065-2
    [15]

    CHANG S J. A simple production function model for variable density growth and yield modelling[J]. Canadian Journal of Forest Research, 1984, 14: 783-788. doi: 10.1139/x84-139
    [16]

    BENTLEY W R, FIGHT R D. A zero interest comparison of forest rent and soil rent[J]. Forest Science, 1966, 12(4): 460.
    [17] 陈平留, 刘 健, 郑德祥. 速生丰产杉木林经济效益及伐期确定[J]. 林业科学, 2001, 37(1):35-47.

    [18] 陈少雄, 周国福, 林义辉. 尾巨桉纸浆材人工林轮伐期研究[J]. 林业科学研究, 2002, 15(4):394-398.

    [19] 黄庆丰, 孙启祥, 吴泽民, 等. 长江滩地杨树人工林主伐年龄的研究[J]. 林业科学, 2002, 38(6):154-158.

    [20] 黄和亮, 吴景贤, 许少洪, 等. 桉树工业原料林的投资经济效益与最佳经济轮伐期[J]. 林业科学, 2007, 43(6):6.

    [21] 相聪伟, 张建国, 段爱国, 等. 杉木人工林材种结构的立地及密度效应研究[J]. 林业科学研究, 2015, 28(5):54-59.

    [22] 段爱国, 付立华, 张建国. 基于长期密度试验的杉木人工林自然稀疏线研究[J]. 资源与生态学报, 2019, 10(3):315-323.

    [23] 孙洪刚, 张建国, 段爱国. 数据点选择与参数估计方法对杉木人工林自疏边界线的影响[J]. 植物生态学报, 2010, 34(4):9.

    [24] 张建国, 段爱国. 理论生长方程与直径结构模型的研究[M]. 北京: 科学出版社, 2004.

    [25]

    CHANG S J. Rotation age, management intensity, and the economic factors of timber production: do changes in stumpage price, interest rate, regeneration cost, and forest taxation matter[J]. Forest Science, 1983, 29: 267-277.
    [26]

    ZHOU W. Optimal method and optimal intensity in reforestation[D]. Doctoral thesis Swedish Univ. of Agricultural Sciences, Umea, Sweden. 1999.
    [27]

    COORDES R. Influence of planting density and rotation age on the profitability of timber production for Norway spruce in central Europe[J]. European Journal of Forest Research, 2013, 132: 297-311. doi: 10.1007/s10342-012-0675-9
    [28]

    CAO T, HYYTIAINEN K, TAHVONEN O, et al. Effects of initial stand states on optimal thinning regime and rotation of Picea abies stands[J]. Scandinavian Journal of Forest Research, 2006, 21: 388-398. doi: 10.1080/02827580600951915
    [29]

    CLARK C W, DE PREE J. A simple linear model for optimal exploitation of renewable resources[J]. Applied Mathematics and Optimization, 1979, 5(1): 181-196. doi: 10.1007/BF01442553
    [30]

    KAO C, BRODIE J D. Simultaneous optimization of thinnings and rotation with continuous stocking and entry intervals[J]. Forest Science, 1980, 26: 338-346. doi: 10.1093/forestscience/26.3.338
    [31]

    HYYTIAINEN K, TAHVONEN O. Economics of forest thinnings and rotation periods for finish conifer cultures[J]. Scandinavian Journal of Forest Research, 2002, 17: 274-288. doi: 10.1080/028275802753742945
    [32]

    HALBRITTER A, DEEGEN P. A combined economic analysis of optimal planting density, thinning and rotation for an even-aged forest stand[J]. Forest Policy and Economics, 2015, 51: 38-46. doi: 10.1016/j.forpol.2014.10.006
    [33]

    CAWRSE D C, BETTERS D R, KENT B M. A variational solution technique for determining optimal thinning and rotation schedules[J]. Forest Science, 1984, 30: 793-802.
    [34]

    TAYLOR R G, FORTSON J C. Optimum plantation planting density and rotation age based on financial risk and return[J]. Forest Science, 1991, 37: 886-902.
    [35]

    LOHMANDER P. The economically optimal number of plants, the damage probability, and the stochastic round wood market. 1994, P. 290-314 in Proc. of the Internat. Symp. on Systems analysis and management decisions in forestry, Gonzalo, L. , and V. Paredes (eds. )
    [36]

    REED W J. The decision to conserve or harvest old-growth forest[J]. Ecological Economics, 1993, 8: 45-69. doi: 10.1016/0921-8009(93)90030-A
    [37]

    WILLASSEN Y. The stochastic rotation problem: A generalization of Faustmann formula to stochastic forest growth[J]. Journal of Economic Dynamics and Control, 1998, 22: 573-596. doi: 10.1016/S0165-1889(97)00071-7
    [38]

    INSLEY M. A real options approach to the valuation of a forestry investment[J]. Journal of Environmental Economics and Management, 2002, 44(3): 471-492. doi: 10.1006/jeem.2001.1209
    [39]

    JACOBSEN B J. The regeneration decision: a sequential two-option approach[J]. Canadian Journal of Forest Research, 2007, 37: 439-448. doi: 10.1139/x06-232
    [40]

    NAKAJIMA T, SHIRAISHI N, KANOMATA H, et al. A method to maximize forest profitability through optimal rotation period selection under various economic, site and silvicultural conditions[J]. New Zealand Journal of Forestry Science, 2017, 47: 4. doi: 10.1186/s40490-016-0079-6
    [41]

    MARUTANI T. The effect of site quality on economically optimal stand management[J]. Journal of Forest Economics, 2010, 16(1): 35-46. doi: 10.1016/j.jfe.2009.05.001
    [42]

    HAIGHT R G, SMITH W D. Harvesting loblolly pine plantations with hardwood competition and stochastic prices[J]. Forest Science, 1991, 37(5): 1266-1282.
    [43]

    HARTMAN R. Harvesting decision when a standing forest has value[J]. Economic Inquiry, 1976, 14: 52-58. doi: 10.1111/j.1465-7295.1976.tb00377.x
    [44] 刘 林, 张 旭, 余素君, 等. 湿地松材脂兼用林最优轮伐期的经济分析—以江西省景德镇市枫树山林场为例[J]. 林业科学, 2022, 58(4):62-73.

    [45] 林 卓, 吴承祯, 洪 伟. 杉木人工林碳汇木材复合经济收益分析及最优轮伐期确定-基于时间序列预测模型[J]. 林业科学, 2016, 52(10):134-145.

    [46] 沈月琴, 王 枫, 张耀启, 等. 中国南方杉木森林碳汇供给的经济分析[J]. 林业科学, 2013, 49(9):140-147.

    [47]

    CHLADNÁ Z. Determination of optimal rotation period under stochastic wood and carbon prices[J]. Forest Policy and Economics, 2007, 9(8): 1031-1045. doi: 10.1016/j.forpol.2006.09.005
    [48]

    PRICE C. Optimal rotation with differently-discounted benefit streams[J]. Journal of Forest Economics, 2017, 26: 1-8. doi: 10.1016/j.jfe.2016.10.001
    [49]

    BRAZEE R J. Impacts of declining discount rates on optimal harvest age and land expectation values[J]. Journal of Forest Economics, 2018, 31: 27-38. doi: 10.1016/j.jfe.2017.06.002
    [50]

    DAVIS L S, Johanson K N, Bettinger P, et al. Forest management to sustain ecological, economic and social values[M]. New York: McGraw-Hill. 2001.
    [51]

    THORSEN B J, HELLES F. Optimal stand management with endogenous risk of sudden destruction[J]. Forest Ecology and Management, 1998, 108: 287-299. doi: 10.1016/S0378-1127(98)00233-3
    [52]

    STAUPENDAHL K, MöHRING B. Integrating natural risks into silvicultural decision models: A survival function approach[J]. Forest Policy and Economics, 2011, 13(6): 496-502. doi: 10.1016/j.forpol.2011.05.007
    [53]

    LOISEL P. Impact of storm risk on Faustmann rotation[J]. Forest Policy and Economics, 2014, 38: 191-198. doi: 10.1016/j.forpol.2013.08.002
    [54]

    KILKKI P, Vaisanen U. Determination of the optimal policy for forest stands by means of dynamic programming[J]. Acta Forestalia Fennica, 1969, 102: 100-112.
    [55]

    SOLBERG B, HAIGHT R G. Analysis of optimal economic management regimes for Picea abies stands using a stage structured optimal-control model[J]. Scandinavian Journal of Forest Research, 1991, 6: 559-572. doi: 10.1080/02827589109382692
    [56]

    COORDES R. Thinnings as unequal harvest ages in even-aged forest stands[J]. Forest Science, 2014, 60(4): 677-690. doi: 10.5849/forsci.13-004
    [57] 李子敬, 张守攻, 孙晓梅, 等. 用实物期权法确定日本落叶松纸浆林的最优轮伐期[J]. 林业科学, 2012, 48(5):61-66.

    [58]

    GARDEBROEK, HERNANDEZ. Do energy prices stimulate food price volatility? Examining volatility transmission between US oil, ethanol and corn markets[J]. Energy Economics, 2013, 40(11): 119-129.
    [59]

    NORSTROM C J. A stochastic model of the growth period decision in forestry[J]. Sweden Journal of Economics, 1975, 77: 329-337. doi: 10.2307/3438965
    [60]

    THOMSON T A. Optimal forest rotation when stumpage prices follow a diffusion process[J]. Land Economics, 1992, 68(3): 329-342.
    [61]

    HAIGHT R G, HOLMES T P. Stochastic price models and optimal tree cutting: results for loblolly pine[J]. Natural Resource Modeling, 1991, 5(4): 423-443. doi: 10.1111/j.1939-7445.1991.tb00255.x
    [62]

    CLARKE H R, REED J W. The tree-cutting problem in a stochastic environment: the case of age-dependent growth[J]. Journal of Economic Dynamics and Control, 1989, 13(4): 569-595. doi: 10.1016/0165-1889(89)90004-3
    [63]

    MANLEY B, NIQUIDET K. How does real option value compare with Faustmann value when log prices follow fractional Brownian motion[J]. Forest Policy and Economics, 2017, 85(1): 76-84.
    [64]

    INSLEY, MARGARET, ROLLINS, KIMBERLY. On solving the multirotational timber harvesting problem with stochastic prices: a linear complementarity formulation[J]. American Journal of Agricultural Economics, 2005, 87(3): 735-755. doi: 10.1111/j.1467-8276.2005.00759.x
    [65] 焦媛媛, 韩文秀, 杜 军. 确定最优轮伐期的实物期权方法[J]. 西北农林科技大学学报(自然科学版), 2003, 11:22-29.

    [66] 石海金, 宋铁英. 适应价格的用材林主伐决策模型的研究[J]. 林业科学, 1999, 35(1):17-23.

    [67]

    HYDE W F. Timber Supply, Land Allocation and Economic Efficiency[M]. Johns Hopkins Press, Baltimore. 1980.
    [68]

    HARDIE I W, Daberkow J N, McConnell K E. A timber harvesting model with variable rotation lengths[J]. Forest Science, 1984, 2: 511-523.
    [69]

    PRICE C. Optimal rotation with negative discount rates: completing the picture[J]. Journal of Forest Economics, 2017, 29: 87-93. doi: 10.1016/j.jfe.2017.08.006
    [70]

    HYYTIäINEN K, HARI P, KOKKILA T, et al. Connecting a process-based forest growth model to stand-level economic optimization[J]. Canadian Journal of Forest Research, 2004, 34(10): 2060-2073. doi: 10.1139/x04-056
    [71]

    HYYTIäINEN K, TAHVONEN O, VALSTA L. Optimum juvenile density, harvesting and stand structure in even-aged Scots pine stands[J]. Forest Science, 2005, 51(2): 120-133.
    [72]

    TEETER L D, CAULFIELD J P. Stand density management strategies under risk: effects of stochastic prices[J]. Canadian Journal of Forest Research, 1991, 21: 1373-1379. doi: 10.1139/x91-194
    [73]

    GONG P, BOMAN M, MATTSSON L. Non-timber benefits, price uncertainty and optimal harvest of an even-aged stand[J]. Forest Policy and Economics, 2005, 7(3): 283-295. doi: 10.1016/S1389-9341(03)00073-X
    [74]

    BRAZEE R J, DWIVEDI P. Optimal forest rotation with multiple product classes[J]. Forest Science, 2015, 61(3): 458-465. doi: 10.5849/forsci.13-207
    [75]

    HIRSHLEIFER J. Investment, interest, and capital[M]. Prentice-Hall, Upper Saddle River. 1970.
    [76]

    CAULFIELD J P, SOUTH D B, SOMERS G L. The influence of the price size curve on planting density decisions[C]. Sixth biennial southern silvicultural research conference, 1990:801-810, Memphis.
    [77]

    GONG P. Determining the optimal planting density and land expectation value- a numerical evaluation of decision model[J]. Forest Science, 1998, 44: 356-364.
    [78]

    LIU L, ZHANG J G, LI Y, SUN H G. Impact of initial planting density on the optimal economic rotation of Chinese Fir (Cunninghamia lanceolata (Lamb. ) Hook) in an experimental forest plantation[J]. Forests, 2019, 10(9): 713. doi: 10.3390/f10090713
    [79]

    SAMUELSON P A. Economics of forestry in an evolving society[J]. Economic Inquiry, 1976, 14(4): 466-492. doi: 10.1111/j.1465-7295.1976.tb00437.x
    [80]

    TAHVONEN O, PIHLAINEN S, NIINIMäKI S. On the economics of optimal timber production in boreal Scots pine stands[J]. Canadian Journal of Forest Research, 2013, 43: 719-730. doi: 10.1139/cjfr-2012-0494
    [81]

    HOOL J N. A dynamic programming-Markov Chain approach to forest production control[J]. Forest Science, 1966, 12(3): a0001-z0001.
    [82] 刘景芳, 童书振. 全国杉木人工林间伐表编制的研究[J]. 林业科学研究, 1996, 9(2):152-157.

    [83] 张水松, 陈长发, 吴克选, 等. 杉木林间伐强度试验20年生长效应的研究[J]. 林业科学, 2005, 41(5):56-65.

    [84]

    CAMPELL H. Land and forest economics[M]. Edward Elgar Publishing, Glensanda House. 1999.
    [85]

    PRICE C. Optimal rotation with differently-discounted benefit streams[J]. Scandinavian Journal of Economics, 2014, 43: 48-54.
    [86]

    AMACHER G S, OLLIKAINEN M, KOSKELA E. Economics of Forest Resources[M]. The MIT Press, Cambridge, MA. 2009.
  • [1] 徐宏远郑世锴 . Ⅰ-72杨人工林的投入产出及经济成熟龄. 林业科学研究, 1991, 4(5): 539-544.
    [2] 王葆芳朱灵益 . 干旱区杨树人工用材林灌溉量与土壤水分的关系*. 林业科学研究, 1995, 8(1): 94-100.
    [3] 张勇仲崇禄陈羽陈珍姜清彬韩强方发之 . 海南木麻黄无性系生长过程研究. 林业科学研究, 2017, 30(4): 588-594. doi: 10.13275/j.cnki.lykxyj.2017.04.008
    [4] 李贻铨 . 林木施肥是短轮伐期工业用材林的基础技术措施. 林业科学研究, 1992, 5(2): 214-218.
    [5] 姜景民虞沐奎童方平李江南刘光正胡世才 . 湿地松、火炬松工业用材林造林密度初步研究. 林业科学研究, 2000, 13(2): 167-176.
    [6] 马晖于卫平牛长明郭宏玲黄利江王挺张广才王培录 . 宁夏腾格里沙漠边缘杨树短轮伐期工业用材林造林技术. 林业科学研究, 2004, 17(s1): 109-113.
    [7] 梁军吕士行徐锡增 . I-69杨胶合板用材林密度及立地效应研究. 林业科学研究, 1996, 9(4): 348-353.
    [8] 胡小燕段爱国张建国杜海伦张雄清郭文福郭光智 . 南亚热带杉木人工成熟林密度对土壤养分效应研究. 林业科学研究, 2018, 31(3): 15-23. doi: 10.13275/j.cnki.lykxyj.2018.03.003
    [9] 王保平李宗然李芳东周海江周道顺 . 中龄兰考泡桐胶合板材林营养补给效应的研究. 林业科学研究, 1997, 10(2): 220-225.
    [10] 赵庆霞包维楷 . 华山松、油松不同叶龄针叶光合特性对林窗式疏伐的响应. 林业科学研究, 2018, 31(5): 104-109. doi: 10.13275/j.cnki.lykxyj.2018.05.014
    [11] 赵宝瑄 . “用材林基地立地分类、评价及适地适树研究”在北京通过成果鉴定. 林业科学研究, 1991, 4(6): 649-649.
    [12] 姚甲宝曾平生袁小平吴建国楚秀丽周志春 . 间伐强度对木荷-萌芽杉木中龄混交林生长和林分结构的影响. 林业科学研究, 2017, 30(3): 511-517. doi: 10.13275/j.cnki.lykxyj.2017.03.021
    [13] 苏上倪建伟许可耿涌杭王伟许新桥 . 香椿成熟复叶饲用价值分析评价研究. 林业科学研究, 2017, 30(5): 848-853. doi: 10.13275/j.cnki.lykxyj.2017.05.019
    [14] 汪奎宏裘福庚蔡纫秋 . 毛竹主要营林措施技术经济效果分析. 林业科学研究, 1991, 4(1): 30-37.
    [15] 徐海东苑海静熊静虞木奎成向荣 . 杉阔异龄复层林对土壤团聚体稳定性和有机碳及养分储量的影响. 林业科学研究, 2020, 33(3): 107-115. doi: 10.13275/j.cnki.lykxyj.2020.03.014
    [16] 姚小华方嘉兴杨耀仙盛能荣颜国栋 . 余甘子二龄结果母枝强修剪对一龄结果母枝的影响. 林业科学研究, 1989, 2(2): 155-162.
    [17] 王庆海武菊英滕文军刘亚平 . 种植方式和越冬保护措施对芦竹成活率的影响. 林业科学研究, 2006, 19(6): 813-815.
    [18] 梁军徐锡增吕士行 . I-69杨胶合板用材林产量及生物量结构. 林业科学研究, 2001, 14(2): 154-159.
    [19] 姜景民虞沐奎童方平 . 湿地松、火炬松工业用材林间伐制度研究. 林业科学研究, 2000, 13(4): 397-406.
    [20] 曹光球费裕翀路锦黄樱郑宏林开敏季春杉曹世江 . 林下植被不同管理措施培育杉木大径材林分土壤酶活性差异及质量评价. 林业科学研究, 2020, 33(3): 76-84. doi: 10.13275/j.cnki.lykxyj.2020.03.010
  • 加载中
计量
  • 文章访问数:  1694
  • HTML全文浏览量:  602
  • PDF下载量:  58
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-08-17
  • 录用日期:  2023-10-10
  • 网络出版日期:  2024-02-27
  • 刊出日期:  2024-04-27

基于试验数据的人工用材林最优轮伐期确定的研究进展

    通讯作者: 孙洪刚, honggangsun@caf.ac.cn
    作者简介: 刘林,中国计量大学经济与管理学院副教授,博士。研究方向:林业生产效率。电话:13634119945 Email:jhliulin@cjlu.edu.cn
  • 1. 中国计量大学 经济与管理学院,浙江 杭州 310018
  • 2. 中国林业科学研究院 亚热带林业研究所,浙江 杭州 311400
  • 3. 中国林业科学研究院 林业研究所,北京 100091

摘要: 虽然用材林轮伐期研究长期受到国内外众多学者的关注,但有关轮伐期的内涵、轮伐期确定的标准、轮伐期模型选择以及相关变量对轮伐期的影响等方面,依然是众说纷纭,无法为构建用材林提质增效培育模式提供坚实的理论基础。本文通过梳理国内外有关轮伐期确定方面的研究成果,并以我国杉木用材林试验数据作为支撑,意在展示不同轮伐期确定准则、方法和相关影响因素对轮伐期确定的作用机理和过程,并对未来轮伐期的研究方向进行了展望。

English Abstract

  • 我国政府从上世纪80年代初“六.五”林业科技支撑计划项目开始,围绕杉木(Cunninghamia lanceolata (Lamb.) Hook.)、马尾松(Pinus massoniana Lamb.)和杨树(Populus)等主要用材树种,逐步构建了以“遗传控制、立地控制、密度控制、苗木质量控制和轮伐期控制”为核心的用材林培育技术体系,林分最优轮伐期的确定是这一培育体系中的重要组成部分。在确定用材林轮伐期的采伐年龄时,主要依据林分数量成熟龄、工艺成熟龄和经济成熟龄。上世纪90年代,为满足国民经济建设对木材的需求,国内众多学者提出,要在考虑立地潜能的基础上,依据“以工艺成熟为基限,重点考虑经济成熟,适当兼顾数量成熟”的原则来确定用材林的主伐年龄[1-3]。进入21世纪,为了响应全球加快造林增汇以应对气候变化的要求,我国林业发展目前致力于提高森林质量,发挥森林多功能效益,这就需要建立以森林经营方案为基础的森林多功能培育技术体系。如何科学确定用材林的轮伐期又成为培育技术体系中的关键问题。

    虽然国内外有关轮伐期的研究数量众多,但有关轮伐期的内涵、轮伐期确定的标准、不同类型轮伐期模型的选择以及相关变量对轮伐期的影响等方面,依然是众说纷纭,无法为构建用材林提质增效培育模式提供坚实的理论基础。本文通过梳理国内外有关轮伐期确定方面的研究成果,并以我国杉木用材林试验数据作为支撑,意在展示轮伐期确定的不同准则、方法和相关影响因素对轮伐期确定时的影响机理和过程,并对未来轮伐期的研究方向进行展望。

    • 根据不同的林分产出目标,目前主要有6 个用材林最优轮伐期的确定准则[4]:(1)最大化总产出(Maximum Gross Yield,MGY)准则;(2)最大化可持续产出(Maximum Sustainable Yield,MSY)准则;(3)单一轮伐期收入净现值最大(Present Net Worth ,PNW)准则;(4)林地期望价值(Land Expectation Value,LEV)最大准则,也被称为无限轮伐期收入净现值最大准则或无限折现收益准则[5];(5)年净收入最大(Maximum Annual Net Revenue,MANR)准则,以及(6)资本收益率最大(Internal Rate of Return,IRR)准则。前两个准则以林分生长量为要素,不考虑林分培育过程中产生的成本和林地价值。而后4个准则兼顾林分生长和经济因素。本文基于上述6个准则建立最优化目标函数,通过数理解析得到在6个准则下轮伐期确定的最优条件(见表1),并以40年生杉木密度试验林 1为例,比较不同准则下轮伐期的异同 2(见表2)。

      表 1  轮伐期确定的6个准则及确定条件

      Table 1.  Conditions for determining rotation under 6 criteria

      准则
      Criteria
      目标函数
      Objective function
      轮伐期确定条件
      Conditions for determining rotation
      文献
      References
      MGY max(T)
      {pQ(T)}
      pQT = 0 = QT [6]
      MSY max(T)
      {pQ(T)/T}
      pQT/pQ(T) = 1/T [7-9]
      PNW max(T)
      {pQ(T)e−rTwE}
      QT/Q(T) = r [10-11]
      LEV max(T)
      {[pQ(T)e−rTwE](1- e−rT )}
      pQT = rpQ(T) + rLEV [12-15]
      MANR max(T)
      {(pQ(T) - wE)/T}
      l/T = pQ(T)/[pQ(T)-wE] [16]
      IRR max(i)
      {[pQ(T)e−iTwE]/(1-e−iT)} = 0
      pQT = i × pQ(T) [1-3,17-20]
      注:T为林分年龄;p是活立木价格;r为折现率;i为内部收益率;wE为林地更新成本;Q(T)为林分生长量,是林分年龄的凹函数;QT为林分生长量对T的一阶导[4]
        Notes: T represents the stand age; P is the price of stumpage; r is the discount rate; i is the internal rate of return; wE is the cost of forest land renewal; Q(T) is the stand growth, which is a concave function of stand age, QT is the first derivative of Q(T) on T

      表 2  6个最优化条件下求解杉木用材林最优轮伐期

      Table 2.  Optimal rotations for Chinese fir plantation under 6 criteria

      林分年龄
      Age
      Q(T)QTQ(T)/TQT/Q(T)QT/[Q(T)−wE/p]i*1/Tr/(1−e−rT)
      749.47624.9677.0680.5051.5470.0090.1430.169
      975.27824.3428.3640.3230.5800.0810.1110.138
      1199.59122.2459.0540.2230.3360.1160.0910.118
      13121.52419.7919.3480.1630.224
      0.1240.0770.105
      15140.95917.4329.3970.1240.1620.1200.0670.095
      17158.08015.3219.2990.0970.1230.1180.0590.087
      19173.15713.4899.1140.078
      0.096
      0.1120.0530.081
      21186.47011.9198.8800.0640.0790.1050.0480.076
      23198.27110.5798.6200.0530.0640.1010.0430.073
      25208.7799.4348.3510.0450.0520.0940.0400.070
      27218.1798.453
      8.0810.0390.042
      0.0910.0370.067
      29226.6287.3097.8150.0330.0350.0890.0340.065

      45326.1030.0020.0000.000

      表2可见,在不同准则下计算出的最优轮伐期明显不同。MGY准则下的轮伐期最长,由此确定的轮伐期其实就是用材林林分生长曲线的最高点,此后林分产出随林分年龄增加而逐步下降。在生产实践中很少依据该准则确定轮伐期,原因在于轮伐周期太过漫长,无法满足木材即时需求,尤其在生长缓慢的欧洲北部森林,据此确定的轮伐期会超过100年。IRR准则下的轮伐期最短,因为依据这一准则计算得到的资本收益率要与市场收益率相比较,市场收益率越高,依据该准则所确定的轮伐期越短[3]。MSY准则和MANR准则在该例中结果接近,在其它研究中,依据这两个准则计算得到的轮伐期也较为接近[4]。基于PNW准则的计算结果要与折现率相比较,折现率越高,依此准则确定的轮伐期越短,本例中将折现率设置为0.05时,轮伐期为23年。

      在生产实践中具体采用哪个准则,要根据用材林的培育目标来确定。就目前用材林轮伐期研究来看,森林培育方向的研究者大多采用MSY准则计算用材林的数量成熟龄,林业经济的研究者主要采用LEV准则计算用材林的经济成熟龄。因此,下面主要讨论这两个准则下的轮伐期确定方法及其影响因素。

    • 数量成熟龄主要根据MSY准则来确定主伐年龄。基于林分最大可持续产量理论,数量成熟龄确定的传统方法是,将林分材积的年平均生长量(Mean Annual Increment,MAI)等于连年生长量(Cumulative Annual Increment,CAI)时的林分年龄作为主伐年龄[10]。CAI函数和MAI函数均是从林分生长方程中推导出来的,CAI函数为生长方程的一阶导数,MAI为生长方程Q除以林分年龄t。虽然基于不同林分生长方程推导出的CAI和MAI函数形式不同,但都是依据CAI = MAI来确定数量成熟龄。从表达式可以看出,林分生长方程的选择对数量成熟龄的确定有决定作用。表3给出常用生长方程的形式以及利用该生长方程的数量成熟龄确定条件。

      表 3  常用生长方程的数量成熟龄确定条件

      Table 3.  Condition for quantitative mature age under difference growth functions

      序号
      No.
      生长方程
      Growth function
      方程形式
      Equations
      数量成熟龄确定条件
      Condition for quantitative rotation
      1 Compertz模型 $ Q=kexp(-{e}^{a-bt}) $ $ {e}^{(a-bT)}={\left(bT\right)}^{-1} $
      2 Logistic模型 $ Q=\dfrac{k}{1 + {e}^{(a-bt)}} $ $ {e}^{(-a + bT)}=bT-1 $
      3 Mitschelich模型 $ Q=k(1-a{\mathrm{e}}^{-bt}) $ $ {e}^{bT}=abT + a $
      4 Richards模型 $ Q=A(1-b{e}^{-kt}{)}^{1/1-m} $ $ {e}^{bT}-b=\dfrac{kbT}{1-m} $
      5 Weibull模型 $ Q=A\left\{1-\mathrm{e}\mathrm{x}\mathrm{p}\left[-{\left(\dfrac{t}{b}\right)}^{c}\right]\right\} $ ${e}^{{\left(\frac{T}{b}\right)}^{c} }=1 + \dfrac{c}{b}{\left(\dfrac{T}{b}\right)}^{c-1}\cdot T$
      注:t为林龄;T 是数量成熟龄;kQ的上渐近值;a为与Q初植有关的参数;b为内禀生长率;c为形状参数;m为曲线形状及拐点位置参数
        Notes: t represents the stand age; T is the quantitative rotation; k is the upper asymptotic value; a is the parameter related with the initial Q value; b is the intrinsic growth rate; c is the shape parameter; m is the curve shape inflection point position parameter

      利用5个生长方程,拟合不同初植密度下40年生杉木用材林生长进程并求解轮伐期(表4)。5个方程下的轮伐期均随着林分初植密度的增大而增加,不同方程下轮伐期之间的差异随着初植密度的增大而增加。当初植密度较小时,生长方程的形式对确定轮伐期影响不大。当初植密度达到6 666 株·ha−1时,用Weibull方程求算的轮伐期要比用Mitschelich方程计算的轮伐期短8年,比用Compertz方程计算的轮伐期短12年。因此,在利用数量轮伐期来确定用材林主伐年龄时,造林密度因素是必须要考虑的因素。相聪伟等[21]和段爱国等[22]对杉木密度效应进行研究的结果也证明了这一观点。

      表 4  杉木人工林不同密度下最优轮伐期

      Table 4.  Optimum rotation of Chinese fir plantation under different density

      初植密度
      Initial planting density
      (Number·ha−1)
      生长方程形式
      Functions of growth
      CompetzLogisticMitschelichRichardsWeibull
      16661617181717
      33331819191918
      50002523282425
      66663829342726
      100003829322729

      利用不同生长方程求算的轮伐期结果之间差异很大。这说明数量成熟龄的确定与生长方程形式有关。在生产实践中确定最优轮伐期时,对生长方程的选择要慎重。另外,生长方程参数估计值与采用的估计方法也有关系,如何选择参数估计方法参见孙洪刚等[23]的相关研究,这里不再赘述。理论上,生长方程的选择应遵循如下原则:要充分反映林分真实生长进程;对不同培育措施(如,林分密度、立地质量和抚育间伐等)具有相容性;模型结构要尽可能简单,且模型参数具有明显的生物学意义[24]

    • 经济成熟龄主要依据LEV准则来确定轮伐期。1849年Faustmann在假定林地无限期连续轮伐作业的基础上,将木材收益增加值与木材延迟采伐所增加成本相等时的林分年龄,作为林分最优轮伐期,此时生产的边际收入等于投入的边际成本[25]。Faustmann模型的基本形式是:

      $ \mathit{L}\mathit{E}\mathit{V}=\underset{\mathit{t}}{\mathit{Max}}\frac{{\mathit{e}}^{-\mathit{r}\mathit{t}}\mathit{p}\mathit{V}\left(\mathit{t}\right)-\mathit{c}}{1-{\mathit{e}}^{-\mathit{r}\mathit{t}}} $

      上式中, p 是活立木价格;V为林分材积,它是林分年龄t的函数;c 是更新成本;r是折现率。在Faustmann模型框架下,最优轮伐期就是林地期望价值最大时所对应的林分年龄。由于该模型包含林分机会成本和土地机会成本,被认为是“计算最优轮伐期的最准确的模型”[8]。因此,林地期望价值模型成为众多研究最优轮伐期的理论基础。Faustmann假设活立木价格、林分材积、更新成本和折现率在各轮伐期之间都是相等的,由此导致了在不同轮伐期主伐年龄相同的问题。更为重要的是,这些假设也隐含了林地所有者期望林地价值永远不变的思想[25]。事实上,由于上述假定不变的因素在同一轮伐期和不同轮伐期都会发生变化,对土地期望价值和轮伐期都会产生重要的影响。为此,研究人员对Faustmann模型进行了扩展,主要通过加入影响经济成熟龄的各类影响因素,如林分生长因子、培育措施因子、经济因素以及各种风险或者不确定性因素。表5总结了近些年国内外学者对Faustmann模型的扩展方向以及相关主要研究结论。

      表 5  对Faustmann模型进行扩展的研究

      Table 5.  Expansions for Faustmann model

      扩展方向
      Direction
      主要研究结论
      Main conclusions
      文献
      References
      考虑初植密度和间伐 初植密度对轮伐期有重要的影响,初植密度相同时,没有间伐的林分轮伐期比间伐林分的轮伐期短。 [25-33]
      考虑活立木价格变化 立木价格增加,轮伐期相应缩短。立木质量越好(如,径阶越大),价格越高,如果此时提高立木的价格可能会延长轮伐期。 [34-39]
      加入立地条件 立地条件对轮伐期的影响表现出非单调性。立地条件较差的林分,轮伐期对盈利影响不大,但是立地条件好的林分,轮伐期的确定对盈利有很大的影响。 [40-42]
      加入非木材收益 加入非木材收益后倾向于延长轮伐期(如,考虑碳存储价值)。 [43-47]
      考虑折现率的变化 当折现率较高时,林地价值对轮伐期的变化更敏感。折现率较低时,林地价值不稳定。 [8,29,33,47-50]
      加入不确定性因素 森林价值对多种风险和不确定性很敏感。木材价格的变化比其它风险(如,风灾、病虫害等)造成的经济损失更大。但价格的不确定性,并不意味着对轮伐期有显著影响。 [51-56]

      为了展示利用Faustmann模型确定的经济成熟龄,本文使用40年杉木密度试验林定期观测数据,在计算时考虑不同径阶杉木活立木价格差异,求解5个初植密度下不同折现率时的最优经济轮伐期(见表6)。结果表明,随着初植密度增加,最优经济轮伐期的变化呈非单调性。初植密度较小时,随着初植密度增加,最优经济轮伐期缩短;初植密度较大时,随着初植密度增加,最优经济轮伐期反而会延长。当折现率较低时,不同初植密度林分的最优经济轮伐期差异较小。

      表 6  不同初植密度和折现率时杉木人工林的最优经济轮伐期

      Table 6.  Optimal rotations for Chinese fir plantation using different initial planting densities and discount rates

      初植密度
      Initial planting density (Number·ha−1)
      r = 0.01r = 0.02r = 0.03r = 0.04r = 0.05r = 0.06r = 0.07r = 0.08r = 0.09
      166720181717161513118
      333319161515131211109
      5000181515141313111110
      6667191817171614131313
      10000201918181717161515
    • 表4表6可以看出影响用材林最优轮伐期的因素很多,本文对学者们讨论的较多的影响因素进行梳理。

    • Faustmann模型要求所有变量必须是时间的确定性函数,从而得出活立木价格增加将会缩短轮伐期的结论。而现实的经济社会中,活立木的价格具有不确定性,是随机波动的。因此,有学者将研究活立木未来价格的不确定性以及与价格波动相关的实物期权理论应用到林业生产研究中[57-58]。利用实物期权理论分析林地价值的研究可以分为两类,一类视活立木价格是平稳的[59-60];另一类将活立木价格视为非平稳的[36,61-62],并用分数布朗运动[63]、自回归条件异方差效应[64]、几何布朗运动[65]、自相关随机过程[66]来描述价格的变化。但Jacobsen[39]指出,利用实物期权理论计算得到的最优轮伐期不能同Faustmann模型计算出的轮伐期进行比较。因为如果不考虑活立木价格服从何种数据生成过程,单独考虑价格对轮伐期的影响,会得出当活立木价格增加,轮伐期将缩短的结论。但是如果考虑与价格相关的其他因素时,价格与轮伐期之间的关系就不能确定了。比如,Hyde[67]就发现,由于价格水平和价格增长率之间的相互影响,价格提高所带来的效应是不确定的。Hardie,Daberkow和Mc Connel[68]在考虑活立木价格和生产成本都随时间变化时,假设允许林地在未来可以另作它途,发现最优轮伐期与价格和折现成本之间的关系也是不确定的。

      在确定最优轮伐期时,价格和木材质量之间的关系也被认为是制定主伐和培育模式的重要依据[69]。在多数早期的研究中,活立木不同径阶(木材质量)间的价格差异被忽视了[34],直到Zhou[26]分析欧洲赤松(Pinus sylvestris L.)木材质量对木材价格的影响时,木材质量才在最优轮伐期的研究中考虑进来。Hyytiäinen等人 [70-71] 在单木生长模型中加入了木材质量对木材价格的影响,证明了初植密度、木材质量、木材价格之间是密切相关的。如果考虑了木材质量,价格与最优轮伐期的关系与标准的Faustmann模型的结论是相反的[72-73]。比如,Brazee和Dwivedi[74]发现不同规格木材的价格会直接影响最优轮伐期,尤其是木材径阶越大,价格越高,而木材价格提高会延长最优经济轮伐期。由此可见,有关价格和轮伐期之间的结论都是在一定假设条件下的,条件不一样,结论也不一样。

    • 尽管初植密度对林分生长动态有非常重要的影响,但是有关初植密度对林地价值的影响却鲜见报道[68]。初植密度不仅仅决定木材生产利润,还可能决定其它林分培育措施应用的规模和范围。如何反映初植密度对轮伐期的影响也就成为不可回避的重要问题之一。Chang[25] 在Faustmann和Hyde[67]以及 Hirshleifer[75]这些早期研究的基础上,为后来的学者建立了数量分析框架,并得到了广泛的认可和应用[52,76-77]。比如,Liu 等人 [78]在Chang的框架下提出确定初植密度对最优经济轮伐期的作用时要把握两个关键:一是林分材积增长率与折现率的大小关系,如果林分材积增长率高于折现率时,初植密度对最优经济轮伐期具有负向影响,反之为正向;另一个关键是初植密度对木材质量的影响,低密度林分会随林分年龄增加而立木价格差异大,对最优轮伐期的影响也大,并认为在人工成本日益增加的情况下,采用皆伐生产木材的林分应采用低密度造林,但最低密度如何确定还有待于进一步的研究。

      间伐是林分密度调控的重要手段。虽然Faustmann在最初的LEV模型中考虑了间伐因素,但Samuelson[79]为了方便对其模型进行经济分析,并没有考虑间伐因素。此后,多数经济学家在研究轮伐期时也将间伐因素摒弃。Chang[25]分析初植密度和轮伐期之间的关系时也没有涉及间伐。因为他怀疑任何对初植密度与轮伐期模型的扩展,比如加入间伐因子,将会导致模型更高的复杂性,而这种复杂性将可能不会得到有意义的结果。但他同时指出,间伐对最优轮伐期的影响是非常重要的,仅对林地价值的贡献就可达到40%[80]。为此,关于初植密度、间伐和轮伐期三者关系的研究再次受到重视[28]。目前,此类研究大都应用动态规划[32]或者庞特里亚金(Pontryagin)的最大化原则[29,31,33,56],比如,Hool[81]应用动态规划马尔柯夫链研究木材间伐与主伐问题时,在假设木材生长具有随机性情况下,提出了最优间伐和皆伐策略。Clark 和De Pree[29]以及Cawrse 等人[33]在利用Faustmann-Pressler-Ohlin理论比较有无间伐对轮伐期的影响时,假设两者的初植密度是相同的,并推断出没有间伐的轮伐期比有间伐林分的轮伐期短。而从树种类型来说,有研究认为,针叶树应该进行高密度造林,然后通过多次间伐,获得较高的间伐收益[77,52-53,82]。而有些研究结论恰好相反,认为针叶树应该直接用最终皆伐密度进行低密度造林[26,78]。事实上,随着劳动力价格逐年增加,间伐的人工成本越来越高,甚至超过间伐收入,因此,在中国的很多杉木林区,出现了即使林分已经郁闭甚至发生自疏也不进行抚育间伐的现象。张水松等通过杉木抚育间伐强度试验,认为间伐可以促进杉木直径生长和单株材积生长,但不能提高林分总出材量[83]。因此,在林分总出材量稳定的情况下,提高林业生产效率的唯一途径就是通过减少劳动力投入以节约生产成本。

    • 森林从苗木繁育、整地造林、除草施肥、抚育间伐、日常管护到主伐整个营造林生产过程所支出的直接费、管理费和固定资产折旧费等,都应计入森林培育成本。因此,不同树种的用材林,从育苗造林到木材采伐,几乎年年都产生生产费用,构成生产成本。关于生产成本,Campell H.[84]指出:直接成本不会影响最优轮伐期,因为它是固定投入。但是,它会影响裸地造林的经济可能性。所有未来轮伐收入的现值必须高于直接成本,造林才会有收益,生产单位才会从事营造林生产。关于用材林培育的直接成本变化也是很复杂的。有些成本与初植苗木数量呈正相关,比如初植密度越大,则造林成本越高,但高密度造林会降低抚育成本;间伐强度越大,人工成本越高,但木材质量越好,相应木材收益也越大。因此,培育措施对成本的影响在确定最优经济轮伐期时应该予以考虑。

      另外,折现率或者利率对最优轮伐期和土地期望值的影响也是不容忽视的[31,48-49,54,85]。在典型的Faustmann模型下,采伐年龄在各轮伐期是不变的。Amacher 等人[86]在此条件下利用比较静态分析证明了折现率提高会缩短轮伐期;更新成本增加会延长轮伐期。但是这些结论都是在假设成本与折现率是确定性的条件下得出的,事实上,林业生产周期长,经济因素随着时间发生变化,必然会影响森林培育决策,经济因素的不确定性将使得对最优轮伐期的研究更加复杂。目前,多数学者先使用确定的成本和折现率计算轮伐期,然后在敏感性分析中将成本、折现率的变化进行主观设置,进而分析其变化对轮伐期的影响。如果把用材林培育视为一种长期投资,那么成本、折现率这些经济因素的不确定性如何体现在最优轮伐期的理论模型中,将值得关注。

    • 通过上述理论分析和实证结果可以看出,最优轮伐期是特定立地条件、经济环境以及定向培育措施共同作用下的结果。传统上,最优轮伐期的应用仅限于在外部给定的条件下来确定主伐年龄。然而,当进行初始造林时,森林所有者通常面临着初植密度大小、成本投入高低的不同选择,并且随着林龄的变化,活立木价格、利率等市场情况经常变化,导致最优轮伐期的变化规律极为复杂,未来仍然需要在以下几方面进行探索。

      (1)最优初植密度的确定。用材林最优轮伐期的研究在向林分密度调控方向进行扩展时,一个关键的因素是初植密度。初植密度的重要性表现在:它决定了林分首次间伐时的林分年龄,密度越大,首次间伐时的林分年龄越小;它通过对木材质量的影响使得木材价格在不同径阶上的分布不同;初植密度越大,初植成本越高。因此,确定最优初植密度要在这三个方面实现平衡。简单地说,初植密度就是造林的初始投入,轮伐期就是何时主伐获得产出收益,这种投入产出关系决定了林业生产效率的高低。因此,对由初植密度和轮伐期的变化所导致的森林管理生产应用和理论分析都非常重要。

      (2)不确定性因素对轮伐期的影响。用材林分从造林到主伐,一般要经过10~50年的培育周期。在漫长的生长过程中,林分要经受病虫害和风、雪、火灾等自然灾害,以及由于管理不当所引起的毁林和无法成林等问题,林业投资不能回收的可能性很大,这就要求按复利计算立木价格,以保证林业生产的经济可持续性。另外,在林木不同生长阶段,培育林木的投资也年年不断发生,所形成的利息也年年相应增加。在传统的用材林培育模型中,立木价格、林分培育成本和折现率都假定为时间的确定性函数,但用于实际培育决策则常常是不适宜的,因为市场的变动有很多的不确定性。如何分析这种不确定性对用材林培育决策的影响,制定适应培育策略,对森林培育者是一个很重要的问题。

参考文献 (86)

目录

    /

    返回文章
    返回