留言板

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

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

带约束条件的森林防火最优路径算法

姜广宇 徐爱俊 黄小银 邵建龙

王剑武, 季碧勇, 王铮屹, 等. 浙江省丽水市亚热带森林景观格局对森林碳密度的影响[J]. 浙江农林大学学报, 2024, 41(1): 30-40. DOI: 10.11833/j.issn.2095-0756.20230205
引用本文: 姜广宇, 徐爱俊, 黄小银, 等. 带约束条件的森林防火最优路径算法[J]. 浙江农林大学学报, 2013, 30(1): 76-82. DOI: 10.11833/j.issn.2095-0756.2013.01.011
WANG Jianwu, JI Biyong, WANG Zhengyi, et al. Impact of subtropical forest landscape pattern on forest carbon density in Lishui City of Zhejiang Province[J]. Journal of Zhejiang A&F University, 2024, 41(1): 30-40. DOI: 10.11833/j.issn.2095-0756.20230205
Citation: JIANG Guangyu, XU Aijun, HUANG Xiaoyin, et al. Optimal path algorithm for forest fire prevention with constraints[J]. Journal of Zhejiang A&F University, 2013, 30(1): 76-82. DOI: 10.11833/j.issn.2095-0756.2013.01.011

带约束条件的森林防火最优路径算法

DOI: 10.11833/j.issn.2095-0756.2013.01.011
详细信息
    通信作者: 徐爱俊

Optimal path algorithm for forest fire prevention with constraints

More Information
    Corresponding author: XU Aijun
  • 摘要: 以森林防火中的路径分析算法为研究对象,针对森林防火中对路径分析算法的特殊要求,以Dijkstra最短路径算法为基础,对它们进行改进与优化。结合林区交通网络自身的特点,引入道路宽度、最大限速、必经地和忽略地等森林防火路径分析约束因子,并以此作为道路权值确定的依据,提出了带约束条件的森林防火最优路径算法,并给出详细的算法设计与实现。该算法符合森林防火的实际情况,能够在森林防火应急指挥与决策中发挥作用。图2表2参17
  • 气候变化是人类面临的全球性问题。到2050年,为了将全球升温控制在1.5 ℃以内,必须增加约1亿 hm2森林[1]。森林生态系统是全球陆地生态系统中最大的碳库[2]。森林碳密度是衡量森林生态系统固碳能力的重要指标,其影响机制越来越受到重视,已围绕树种组成[3]、气候驱动变化[4]、森林管护水平[5]、道路网络[6]和城市扩张[7]等自然与人为因素开展了富有成效的研究。

    实际上,景观格局反映了自然与人为因素在不同时空尺度上作用的最终结果,并能够体现各生态过程在不同时空尺度上的相互作用关系[8]。森林景观格局反映了森林景观的组成、空间分布、数量、结构和功能等信息[9],森林景观规模异质性的增加,改变了森林的结构、功能和过程[10],关注森林景观异质性能更好地衡量碳循环等森林生态功能[11]。森林景观如何影响森林储量受到越来越多关注[12],关于森林景观格局对森林碳储量的影响机制也进行了研究。森林景观格局优化可以促进森林碳储量提高[13],扩大斑块面积从而提高森林景观连通性,促进次生林碳储量提升[14],而森林破碎化则会减少森林边缘的碳储量[15]。张丹[16]研究了不同尺度城市森林景观格局与碳储量的相关性,并指出增加森林斑块面积和连结度,减小斑块的平均邻近距离与周长面积比均可增加城市森林碳储量。宋洁[17]分析了祁连山森林景观格局对森林碳储量的相关性,认为增加森林景观面积、景观形状复杂度、景观聚集度和连通性能够提升森林碳储量。吕海亮[18]研究了城市森林碳密度与景观格局关系,认为平均斑块面积与碳密度呈显著正相关,而景观形状指数则与其呈显著负相关。也有学者聚焦森林景观格局对森林碳密度的影响研究,通过耦合协调度模型探究碳密度与森林景观不同组成、分布等特征的关系[19]。然而,量化森林景观格局对森林碳密度影响的研究还不多,且以上研究虽然解析了森林景观格局对森林碳密度的相关性,但很少从空间上解释森林景观对森林碳密度的影响规律。实际上森林碳密度具有空间自相关性[20-21],忽略地理因素的空间非平稳性,易导致结果有偏差或估计效率低下。

    目前,地理加权回归模型(GWR)、空间计量模型、地理探测器等方法已被用于分析研究对象与各类驱动因素之间的空间关联性。其中,GWR解释模型方差比其他模型更有效,且相比于普通最小二乘模型(OLS),能较好地揭示解释变量的空间异质性,但GWR模型是基于恒定带宽,忽略了解释变量尺度作用差异。为此,FOTHERINGHAM等[22]提出了多尺度地理加权(MGWR)模型,解决了不同解释变量在不同尺度同一带宽的问题,即允许在多个空间尺度上建立自变量和因变量的关系模型。MGWR模型已成功应用于生境质量[23]、生态效率[24]、森林火灾[25]等的驱动因素分析,但在森林景观研究中的应用仍处于探索阶段。

    本研究选取了浙江省丽水市为研究区域,采用MGWR模型探究亚热带森林景观格局对森林碳密度的影响机制和空间非平稳性,为研究区域森林碳密度时空特征提供框架,以期为制定国土绿化、森林质量提升政策,因地制宜开展森林保护与管理提供依据。

    丽水市位于浙江省西南部,27°25′~28°57′N和118°41′~120°26′E,地貌以丘陵、中山为主,市域面积为1.73万km2。丽水市属中亚热带季风气候,温暖湿润,雨量充沛,具有典型的山地气候。作为南方亚热带重要的集体林区,及浙江省森林资源最为丰富的设区市,截至2019年,丽水市森林面积为142.14万hm2,森林覆盖率为82.27%,森林蓄积量为0.96亿m3,居浙江省前列。近年来,随着快速的城市化和社会发展,该地区土地利用发生了相当大的变化,这对森林提供生态系统服务能力产生了重要影响[26]。同时该地区国家公园和自然保护地建设,以及森林管理活动也对森林碳循环产生了重大影响[27],增强了森林景观格局和森林碳密度的空间异质性。

    ①土地利用矢量数据来源于浙江省森林资源年度变更成果。将研究区的土地利用分为乔木林地、竹林地、灌木林地、未成林造林地、苗圃地、迹地、宜林地和非林地等8种类型。为更细致刻画森林景观特征,进一步把乔木林地细分为针叶林、阔叶林、针阔混交林等3种类型。地类和树种分类按照DB33/T 640—2017《森林资源规划设计调查规程》[28]执行。②样地数据来源于丽水市716个森林资源连续清查样地(以下简称连清样地)的现场调查,全市以公里格网4 km×6 km间距布设样地,单个样地面积为0.08 hm2,形状为正方形。③数字高程模型来源于地理空间数据云(http://www.gscloud.cn/),分辨率为30 m。④年降水量和年平均气温数据根据资源环境科学与数据中心下载的中国756个气象站观测数据,通过普通Kriging法进行空间插值形成栅格数据,分辨率为30 m。以上数据均取2012和2019年的数据。

    ①利用2012和2019年连清样地实测的样木数据,将每株样木树种归类至表1所列的树种类型。②使用树高-胸径曲线模型计算单株树(竹)高;③使用冠长模型计算单株树(竹)冠长;④使用单株立木(竹)生物量模型计算各样地内的单株活立木(竹)的地上和地下部分生物量,再通过与含碳率的乘积测算单株活立木(竹)的碳储量;⑤将样地内各株活立木(竹)的碳储量累加形成样地林分水平的碳储量,然后除以样地面积得到森林碳密度。⑥利用ArcMap 10.8的球面半变异函数模型的普通Kriging法,形成丽水市森林碳密度空间分布图。树高-胸径曲线模型、冠长模型、生物量模型见表1,含碳率见表2

    表 1  树高-胸径曲线模型、冠长模型和生物量模型
    Table 1  Height-diameter curves model, crown length model and biomass model
    树种类型树高-胸径曲线模型冠长模型生物量模型
    松类 $H = 78.711\,0 + \dfrac{ { - 10\,051.620\,0} }{ {0.984\,8 \; D_{\rm{BH}} + 129.719\,5} }$ $L = 0.908\,2 \; {H^{0.741\,3} }$ ${B_{松类 } } = {B_1} + {B_2} + {B_3}$
    ${B_1} = 0.060\,0 \; {H^{0.793\,4} } \; {D_{\rm{BH} }^{1.800\,5} }$
    ${B_2} = 0.137\,7 \; {L^{0.405\,2} } \; {D_{\rm{BH} }^{1.487\,3} } $
    ${B_3} = 0.041\,7 \; {H^{ - 0.078\,0} } \; {D_{\rm{BH} }^{2.261\,8} }$
    杉类 $H = 119.583\,9 + \dfrac{ { - 24\,448.214\,0} }{ {0.988\,4 \; D_{\rm{BH}} + 205.692\,4} }$ $L = 0.487\,0 \; {D_{\rm{BH}}^{0.170\,7} } \; {H^{0.897\,1} }$ ${B_{ 杉类 } } = {B_1} + {B_2} + {B_3}$
    ${B_1} = 0.064\,7 \; {H^{0.895\,9} } \; {D_{\rm{BH} }^{1.488\,0} }$
    ${B_2} = 0.097\,1 \; {L^{0.034\,6} } \; {D_{\rm{BH} }^{1.781\,4} }$
    ${B_3} = 0.061\,7 \; {H^{ - 0.103\,7} } \; {D_{\rm{BH} }^{2.115\,3} }$
    硬阔类Ⅰ $H = 58.208\,2 + \dfrac{ { - 6\,994.739\,0} }{ {0.984\,8 \; D_{\rm{BH}} + 127.719\,5} }$ $L = 0.631\,6 \; {H^{1.180\,1} } \; {{\rm{e}}^{ - 0.051\,1 \; H} }$ ${B_{ 硬阔类{\text{Ⅰ}} } } = {B_1} + {B_2} + {B_3}$
    ${B_1} = 0.056\,0 \; {H^{0.809\,9} } \; {D_{\rm{BH} }^{1.814\,0} }$
    ${B_2} = 0.098\,0 \; {L^{0.461\,0} } \; {D_{\rm{BH} }^{1.648\,1} }$
    ${B_3} = 0.054\,9 \; {H^{0.106\,8} } \; {D_{\rm{BH} }^{2.095\,3} }$
    硬阔类Ⅱ ${B_{ 硬阔类{\text{Ⅱ}}} } = {B_1} + {B_2} + {B_3}$
    ${B_1} = 0.080\,3 \; {H^{0.781\,5} } \; {D_{\rm{BH} }^{1.805\,6} }$
    ${B_2} = 0.286\,0 \; {L^{0.945\,0} } \; {D_{\rm{BH} }^{1.096\,8} }$
    ${B_3} = 0.247\,0 \; {H^{0.174\,5} } \; {D_{\rm{BH} }^{1.795\,4} }$
    软阔类 $L = 0.441\,3 \; {H^{1.377\,0} } \; {{\rm{e}}^{ - 0.060\,3 \; H} }$ ${B_{ 软阔类 } } = {B_1} + {B_2} + {B_3}$
    ${B_1} = 0.044\,4 \; {H^{0.719\,7} } \; {D_{\rm{BH} }^{1.709\;5} }$
    ${B_2} = 0.085\,6 \; {L^{0.397\,0} } \; {D_{\rm{BH} }^{1.226\;6} }$
    ${B_3} = 0.045\,9 \; {H^{0.106\,7} } \; {D_{\rm{BH} }^{2.024\,7} }$
    乔木经济
     树种类
    $L = 0.618\,9 \; {H^{1.204\,8} } \; {{\rm{e}}^{ - 0.038\,2 \; H} }$ 使用硬阔类Ⅰ的公式
    竹类 $H = 24.557\,0 + \dfrac{ { - 233.809\,9} }{ {D_{\rm{BH}} + 8.434\,3} }$ $L = 0.705\,6 \; {H^{1.174\,8} } \; {{\rm{e}}^{ - 0.060\,1 \; H} }$ ${B_{ 竹类 } } = {B_1} + {B_2} + {B_3}$
    ${B_1} = 0.039\,8 \; {H^{0.577\,8} } \; {D_{\rm{BH} }^{1.854\,0} }$
    ${B_2} = 0.280\,0 \; {L^{0.274\,0} } \; {D_{\rm{BH} }^{0.835\,7} }$
    ${B_3} = 0.371\,0 \; {H^{0.135\,7} } \; {D_{\rm{BH} }^{0.981\,7} }$
      说明:H是树高(m),DBH是胸径(cm),L是冠长(m),B1B2B3分别是单株立木(竹)的树干生物量(kg)、树冠生物量(kg)和树根生物量(kg)。松类包括马尾松Pinus massoniana、湿地松P. elliottii、黄山松P. taiwanensis等树种;杉类包括杉木Cunninghamia lanceolata、水杉Metasequoia glyptostroboides、池杉Taxodiun distichum var. imbricatum等树种;硬阔类Ⅰ包括木荷Schima superba、红楠Machilus thunbergii、樟Camphora officinarum等树干木材密度小于0.7 g·cm−3的树种;硬阔类Ⅱ包括青冈Quercus glauca、苦槠Castanopsis sclerophylla、栎类Quercus等树干木材密度大于0.7 g·cm−3的树种;软阔类包括桤木Alnus cremastogyne、槭树Acer miyabei、檫木Sassafras tzumu等树种;乔木经济树种类包括栗Castanea mollissima、枇杷Eriobotrya japonica、香榧Torreya grandis ‘Merrillii’等树种;竹类包括刚竹Phyllostachys sulphurea、毛竹Ph. edulis、早竹Ph. violascens等。
    下载: 导出CSV 
    | 显示表格
    表 2  树种含碳率
    Table 2  Carbon contents of species
    树种含碳率树种含碳率树种含碳率
    冷杉Abies fabri 0.499 9 柳杉Cryptomeria japonica var. sinensis 0.523 5 杨树Populus 0.495 6
    云杉Picea asperata 0.520 8 水杉Metasequoia glyptostroboides 0.501 3 硬阔类 0.483 4
    铁杉Tsuga chinensis 0.502 2 樟树Cinnamomum camphora 0.491 6 软阔类 0.495 6
    柏木Cupressus funebris 0.503 4 楠木Phoebe zhennan 0.503 0 针叶混交林 0.510 1
    黑松Pinus thunbergii 0.514 6 栎类Quercus 0.500 4 阔叶混交林 0.490 0
    华山松Pinus armandii 0.522 5 桦木类Betula 0.491 4 针阔混交林 0.497 8
    油杉Keteleeria fortunei 0.499 7 椴树Tilia tuan 0.439 2 竹类 0.504 2
    马尾松Pinus massoniana 0.459 6 檫木Sassafras tzumu 0.484 8 乔木经济树种 0.483 4
    高山松Pinus densata 0.500 9 桉树Eucalyptus 0.525 3 其他树种 0.500 0
    杉木Cunninghamia lanceolata 0.520 1 木麻黄Casuarina equisetifolia 0.498 0
    下载: 导出CSV 
    | 显示表格

    为验证碳密度空间分布结果精度可靠性,基于已有数据,先将空间分布结果由栅格转为矢量面要素的格网,并裁切提取出丽水市行政范围内的格网;再将每个格网的面积与其碳密度相乘并汇总加和得到丽水市森林碳储量;然后按照文献[29]所述的系统抽样统计方法,基于连清样地计算丽水市森林碳储量的估计中值、估计精度和估计区间;最后根据文献[30],基于空间分布结果计算的森林碳储量如落入估计区间,则说明该结果是有精度保证的,森林碳储量是碳密度乘以丽水市面积得到,也可说明森林碳密度空间分布结果是有精度保证的。

    选择5个常用且相关性较低的景观指标:香农多样性指数(SHDI)、最大斑块指数(LPI)、斑块密度(PD)、蔓延度指数(CONTAG)、景观形状指数(LSI)。各指标的含义和计算公式见文献[31]。另外,本研究选择了年平均降水量(PREP)、年平均气温(TEMP)、海拔(ELEV)作为控制变量。气温和降水量是影响森林固碳能力最常见的气候要素,其影响是双向的。气温会影响植被光合作用和净初级生产力,从而影响森林生长[32],但高温引起的呼吸增加可能使碳储量随着温度升高而减少[33]。一定范围内的水资源会促进森林的生长,水分过多则会抑制土壤呼吸,影响植物生长[34]。使用Fragstats 4.2软件,参照已有对市级尺度的景观研究[35-36],选取5 km×5 km为窗口大小,采用移动窗口法生成各景观指标栅格图,基于ArcMap 10.8生成4 km×4 km的格网(研究区共计1 204个网格),对每个格网的景观格局指标数值和控制变量数值分区统计。

    地理学第一定律表明,在涉及具有地理空间的数据时,事物之间距离越近关联性越大[37]。采用Global Moran’ s I和Local Moran’ s I来描述森林碳密度的全域和局部聚类特征。计算公式可见文献[38]。Moran’ s I介于−1到1之间,>0表示存在正相关性,<0表示存在负相关性, 0表明不存在空间自相关。使用GeoDa 1.12绘制LISA (local indicators of spatial association)图,反映森林碳密度的空间格局,LISA图包括4个类别,“高—高”和“低—低”表示空间正相关,“高—低”和“低—高”表示空间负相关。Global Moran’ s I也将用于量化MGWR模型的残差空间自相关,如果残差存在显著的空间自相关,则表明模型中缺失了关键解释变量。

    采用MGWR模型探究森林景观格局对森林碳密度的多尺度影响。MGWR模型通过不断寻找各解释变量的最优带宽并不断更新参数估计值,直至系数迭代收敛,从而形成最终回归结果[39]。该模型如下:

    $$ {y_i} = \sum\limits_{j = 1}^m {{\beta _{{bw}j}}\left( {{u_i},{v_i}} \right)} {x_{ij}} + {\varepsilon _i} 。 $$

    其中,yi是格网i的被解释变量;xij是格网i的第j个解释变量;$ {\beta _{{bw}j}}\left( {{u_i},{v_i}} \right) $为格网i的第j个解释变量的局部回归系数,(ui, vi)为格网i的空间位置,bwj为第j个解释变量回归系数所使用的带宽bwm为解释变量的个数;εi为误差项。本研究采用MGWR 2.2进行模型构建,MGWR形式选择Gaussian,空间Kernel类型选择Adaptive Bisquare函数,带宽搜寻方法采用Golden Section,模型最优化准则采用AICc准则。

    2012—2019年,研究区森林碳密度平均值从23.19 t·hm−2提高到31.96 t·hm−2,年均增量为1.25 t·hm−2。根据研究区森林资源连续清查结果,森林面积和森林蓄积量分别从2012年的132.81万hm2和0.65亿m3增加到142.14万hm2和0.96亿m3,年均分别增加1.33万hm2和0.04亿m3,森林面积和森林蓄积量逐年上升促进了森林碳密度增加。森林碳密度空间分布格局存在显著差异,低森林碳密度主要在研究区的北部、西部和西南部的森林碳密度较高(图1)。2012和2019年研究区的Global Moran’ I分别为0.889和0.891 (P<0.01),表明森林碳密度较高或较低区域趋于聚集。由图2发现:①“高—高”区域主要分布在西南部和中南部地区。西南部是瓯江源头,拥有钱江源-百山祖国家公园百山祖园区等多处自然保护地,得益于较少的人为干扰、公益林保护和森林生态保护修复,森林质量高。②“低—低”区域主要分布在北部,北部地区包含松谷平原和碧湖平原,这些地区耕地、建筑用地较多,受人类活动影响频繁。③“低—高”区域分布极少,说明出现周围森林碳密度高而本地森林碳密度低的空间结构概率较低。

    图 1  丽水市森林碳密度空间分布图
    Figure 1  Spatial distribution map of forest carbon density in Lishui City
    图 2  丽水市森林碳密度LISA分布特征图
    Figure 2  Spatiotemporal characteristics of LISA map of forest carbon density in Lishui City

    基于丽水市连清样地,系统抽样统计结果(表3)表明:2012年全市森林碳储量估计中值为4 023.48万t,估计区间为(3 729.05~4 317.90)万t,估计精度为92.68% (P<0.05)。2019年全市森林碳储量估计中值为5 541.19万t,估计区间为(5 192.46~5 889.91)万t,估计精度93.71%为(P<0.05)。2012和2019年基于空间分布结果计算的森林碳储量分别为4 023.08万t和5 543.42万t,均位于相应年份的估计区间内,说明森林碳密度空间分布结果是有精度保证的。

    表 3  分别基于连清样地、森林碳密度空间分布结果计算的丽水市森林碳储量
    Table 3  Forest carbon storage in Lishui City based on the continuous forest inventory and spatial distribution of forest carbon density, respectively
    年份基于空间分布结果的
    计算值/万t
    样地森林碳储量
    均值/t
    基于连清样地的
    计算值/万t
    基于连清样地计算值的
    估计区间/万t
    估计精度/%
    20124 023.081.864 023.48(3 729.05, 4 317.90)92.68
    20195 543.422.575 541.19(5 192.46, 5 889.91)93.71
    下载: 导出CSV 
    | 显示表格
    2.2.1   变量尺度效应

    表4可见:2012和2019年的MGWR模型的Adj-R2分别达0.898和0.892。模型的带宽包含44、79、992和1 204,由于研究区被划分为1 204个网格,因此从行政边界角度看,带宽992和1204的空间尺度可定义为市级,44和79可定义为县级。LSI带宽为44,作用尺度较小,表明景观形状复杂度对森林碳密度的影响在空间上存在较大差异。CONTAG、LPI和SHDI的带宽都较大,说明森林景观连通性、面积优势度和景观多样性对森林碳密度影响的空间差异相对较小,也表明上述3个影响因素回归系数空间异质性较不明显。PD带宽从79升至1 204,作用尺度在扩大,意味着景观破碎化程度对森林碳密度影响的空间差异在变小。

    表 4  MGWR模型性能和回归参数描述性统计
    Table 4  MGWR model performance and descriptive statistics for regression parameters
    年份变量带宽均值标准差最小值中位数最大值正值/%负值/%
    2012 CONTAG 992 −0.107 0.001 −0.108 −0.107 −0.106 0 100
    PD 79 −0.045 0.129 −0.322 −0.077 0.493 23.84 76.16
    LPI 1 204 −0.016 0.067 −0.212 −0.021 0.206 38.04 61.96
    LSI 44 −0.019 0.054 −0.198 −0.017 0.139 27.74 72.26
    SHDI 1 204 −0.069 0.001 −0.071 −0.069 −0.068 0 100
    PREP 44 −0.013 0.573 −1.661 −0.010 1.793 48.59 51.41
    TEMP 44 −0.619 1.335 −7.906 −0.704 1.557 37.54 62.46
    ELEV 44 0.150 0.133 −0.126 0.132 0.554 90.28 9.72
    常数项 44 0.683 0.393 −0.063 0.629 1.990 99.58 0.42
    局部R2 0.916 0.774 0.116 0.178 0.792 0.959
    Adj-R2 0.898
    AICc 968.668
    残差平方和 101.026
    有效参数数量 218.323
    残差Global Moran’ I 0.019 3 (P = 0.092)
    2019 CONTAG 1 204 −0.047 0.001 −0.049 −0.047 −0.045 0 100
    PD 1 204 −0.092 0.099 −0.465 −0.095 0.135 16.69 83.31
    LPI 1 204 −0.004 0.082 −0.190 −0.013 0.309 43.02 56.98
    LSI 44 0.041 0.023 0.007 0.041 0.079 100 0
    SHDI 1 204 −0.055 0.006 −0.066 −0.055 −0.045 0 100
    PREP 44 0.113 0.618 −1.400 0.011 1.376 51.41 48.59
    TEMP 44 0.196 0.550 −1.636 0.210 2.421 65.95 34.05
    ELEV 44 0.152 0.129 −0.162 0.154 0.517 86.88 13.12
    常数项 44 −0.256 0.444 −1.122 −0.196 0.692 34.88 65.12
    局部R2 0.911 0.783 0.116 0.221 0.810 0.966
    Adj-R2 0.892
    AICc 1 018.145
    残差平方和 107.759
    有效参数数量 208.759
    残差Global Moran’ I 0.048 8 (P= 0.002)
      说明:SHDI为香农多样性指数;LPI为最大斑块指数;PD为斑块密度;CONTAG为蔓延度指数;LSI为景观形状指数;PREP为年平均降水量;TEMP为年平均气温;ELEV为海拔。
    下载: 导出CSV 
    | 显示表格
    2.2.2   回归结果

    表4可见:从景观格局指标的回归参数2期均值的绝对值看, CONTAG对森林碳密度影响最大,其次是SHDI和PD,LPI和LSI相对影响较小。从图2的局部R2看,研究区大部分区域拟合优度较高,模型可以揭示研究区2012年森林碳密度总变化的17.8%~95.9%和2019年的22.1%~96.6%。从各影响因素beta系数的正负值占比看,除CONTAG和SHDI外,其余景观格局的参数估计值正负情况都存在,即对森林碳密度表现不同效应影响,说明MGWR模型能将系数的局部特征呈现出来。

    2.2.3   回归结果的空间非平稳性

    图3可见:2012和2019年,CONTAG对森林碳密度的负向影响均从研究区的西北向东南下降,但对森林碳密度的负向驱动力在研究期间逐渐减弱。这可能是因为西北部包含了松古平原,同时遂昌县境内分布有若干小盆地,斜坡平缓,茶叶等经济林种植面积广,应当通过加强农田林网建设和平原绿化增加片状森林面积,降低经济树种的景观连通度,提高区域森林碳密度。PD系数在研究区范围内以负值为主,并且负值区域在逐渐增多,说明降低森林破碎化程度一定程度上可减少对森林碳密度的负向影响。LPI系数以负值为主,意味着大部分区域需要减少人为干扰,维持植被自然生长状态。从2012年到 2019年,LPI系数正值区域有所增加,在西北部和中西部等区域与森林碳密度呈正相关关系,说明这些地方通过适当增加森林抚育、补植造林等措施,能够对提高森林碳密度产生积极影响。从2012年到2019年,LSI对森林碳密度的驱动力逐渐增强。说明在一定程度上,随着斑块形状的复杂度增加,森林碳密度也会随之增加。2019年,LSI对森林碳密度的正向影响由西向东下降,这表明相对于研究区东部,西部森林斑块复杂程度的增加对森林碳密度的正向影响更大。这可能是因为森林景观斑块形状复杂,能够促进斑块边缘的植被形成高效率的能量共生网络[40],增强森林的边缘效应,而斑块边缘通常具有较高的初级生产力[41]。从2012年至2019年,SHDI对森林碳密度的负向驱动力逐渐减弱,形成负向影响由西向东增加的格局,说明研究区的林分种类分布比较均匀。年平均降水量系数总体上呈现由东部大于西部转变为西部大于东部的趋势,但2期正向与负向影响区域的数量差异较小。年平均气温对森林碳密度的正向影响区域逐渐增多,而海拔系数空间分布差异不大。2019年,年平均降水量主要在西部和西南部与森林碳密度呈正相关关系,其余区域主要呈负相关关系。年平均气温系数空间分布与年平均降水量有一定相似,主要在西南部和西北部与森林碳密度呈正相关关系。研究区的绝大部分区域的海拔与森林碳密度呈正相关关系,并且中部区域的正相关关系更强。年平均降水量、年平均气温系数正值分布区域包含了钱江源-百山祖国家公园百山祖园区、九龙山国家级自然保护区、括苍山省级森林公园,充分利用丰富的水热资源,加强自然保护地内森林生态系统保护修复,有助于提高森林碳密度。

    图 3  MGWR模型系数和局部R2空间格局
    Figure 3  Spatial pattern of MGWR model coefficient and local R2

    本研究发现:森林碳密度具有显著的空间自相关性,考虑影响因素的空间尺度对模型的结果会产生重大影响。MGWR模型拟合结果表明,斑块密度、最大斑块指数、年平均降水量、年平均气温和海拔等变量在不同区域的回归系数既有正值也有负值,因MGWR模型使用自适应带宽,可以更精确地探索驱动因素对森林碳密度的作用机制和空间异质性影响,并明确各驱动因素的具体影响尺度。

    森林景观格局对森林碳密度影响表现出不同大小的驱动力和多尺度空间特征,不同尺度需要关注的影响因素是不同的,因此林业管理者可以根据这些特征制定宏观与微观的政策。在市级层面,林业管理者需要加强农田林网和生态廊道建设,降低纯林的景观连通度,同时对生态保护修复进行全区域尺度的管理,将森林恢复和减少森林破碎化战略纳入空间规划中,减少山区森林生态系统中人为干扰和森林破碎化,以提高森林碳密度的整体水平。在县级尺度上,通过合理配置森林景观组成可以有效增强森林固碳量[42],并发挥森林组成成分间相互作用效能,发挥森林边缘效应优势。

    区域森林碳密度计算是一项复杂的研究任务且有多种计算方法,如采用CASA模型计算植被净初级生产力(NPP),再将其转换为固碳量,或采用InVEST模型计算固碳量。森林碳密度受各种因素的影响,本研究是在系统布设的固定样地计算结果基础上,采用空间插值形成区域森林碳密度,计算方法在某些方面可能仍需改进。由于缺少灌木林实测数据,本研究未测算样地内的灌木林碳储量。随着尺度增加,森林碳密度的影响因素也愈加复杂,由于景观格局具有尺度依赖性,MGWR模型作为线性回归模型,对愈加复杂的影响因素解释能力需进一步探讨。此外,本研究使用边长4 km正方形的网格尺度来创建景观格局指标,并探究其与森林碳密度的关系,这个尺度可能更适用于建模,对实际森林管理会有不便。在未来研究中可以以乡镇、村级行政单位为建模单元,进一步增强模型结果的实际应用性。鉴于森林碳密度的空间自相关性,可以考虑使用空间计量模型进一步探究森林碳密度对森林经营管理和经济社会发展的响应机制。

    本研究以浙江省丽水市为研究区,使用MGWR模型探讨了森林景观格局对森林碳密度的影响。 2012和2019年,森林碳密度的空间自相关度分别为0.889和0.891,呈显著的空间自相关性。森林景观格局对森林碳密度的影响在空间上存在非平稳性,CONTAG和SHDI存在负向影响,PD、LPI、LSI在空间上正向和负向影响并存。森林景观格局也表现出不同尺度效应,CONTAG、LPI、SHDI的作用尺度较大,LSI作用尺度较小。基于研究结果,建议进一步重视森林经营管理的空间决策支持,因地制宜实施自然演替和人工干预相结合的森林生态系统保护修复措施,促进森林生态系统碳汇功能提升。

  • [1] 沈鹏, 李功权.  基于生态位因子模型的湖北省松材线虫病风险评估 . 浙江农林大学学报, 2021, 38(3): 560-566. doi: 10.11833/j.issn.2095-0756.20200365
    [2] 叶碧欢, 李海波, 叶卫邦, 陈友吾.  锈色粒肩天牛幼虫的COⅠ、COⅡ、Cytb和28S基因片段序列分析 . 浙江农林大学学报, 2020, 37(2): 303-310. doi: 10.11833/j.issn.2095-0756.2020.02.015
    [3] 褚燕琴, 牛树奎, 陈锋, 芮淑君, 王凯.  火干扰及环境因子对油松林林下植被的影响 . 浙江农林大学学报, 2017, 34(1): 96-103. doi: 10.11833/j.issn.2095-0756.2017.01.014
    [4] 王凯, 牛树奎.  基于Rothermel模型的北京鹫峰国家森林公园潜在火行为 . 浙江农林大学学报, 2016, 33(1): 42-50. doi: 10.11833/j.issn.2095-0756.2016.01.006
    [5] 石晶晶, 张超, 余树全, 茅史亮, 李修鹏.  浙江省防火期内森林火险天气指标动态 . 浙江农林大学学报, 2014, 31(2): 231-237. doi: 10.11833/j.issn.2095-0756.2014.02.011
    [6] 冯雪, 吴达胜, 骆文建.  森林防火多出救点应急调度算法 . 浙江农林大学学报, 2013, 30(2): 257-262. doi: 10.11833/j.issn.2095-0756.2013.02.016
    [7] 秦飞飞, 唐丽华.  基于数字高程模型的森林火灾远程视频同步跟踪算法 . 浙江农林大学学报, 2012, 29(6): 917-922. doi: 10.11833/j.issn.2095-0756.2012.06.017
    [8] 徐爱俊, 黄小银, 姜广宇, 邵建龙.  南方林区森林火险精准区划方法 . 浙江农林大学学报, 2012, 29(1): 65-71. doi: 10.11833/j.issn.2095-0756.2012.01.012
    [9] 黄小银, 徐爱俊, 张晨辉.  森林防火监测型资源配置评价及优化方法 . 浙江农林大学学报, 2011, 28(4): 583-589. doi: 10.11833/j.issn.2095-0756.2011.04.010
    [10] 信晓颖, 江洪, 周国模, 余树全, 王永和.  加拿大森林火险气候指数系统(FWI)的原理及应用 . 浙江农林大学学报, 2011, 28(2): 314-318. doi: 10.11833/j.issn.2095-0756.2011.02.023
    [11] 徐华潮, 骆有庆.  松材线虫入侵对森林生态系统的影响 . 浙江农林大学学报, 2010, 27(3): 445-450. doi: 10.11833/j.issn.2095-0756.2010.03.020
    [12] 仇金宏, 沈明霞, 丛静华, 李龙国.  基于单目视觉的森林火点实时定位方法 . 浙江农林大学学报, 2010, 27(5): 651-658. doi: 10.11833/j.issn.2095-0756.2010.05.003
    [13] 李世友, 罗文彪, 舒清态, 马长乐, 马爱丽, 张桥蓉.  昆明地区25种木本植物的燃烧性及防火树种筛选 . 浙江农林大学学报, 2009, 26(3): 351-357.
    [14] 夏其表, 李光辉, 尹建新.  渗透理论在森林火灾模拟系统中的仿真研究 . 浙江农林大学学报, 2009, 26(2): 233-238.
    [15] 陈培金, 徐爱俊, 邵香君, 刘爱君.  基于GIS的森林火灾灾后评估算法的设计与实现 . 浙江农林大学学报, 2008, 25(1): 72-77.
    [16] 楼雄伟, 方陆明, 徐爱俊, 叶祥根.  多源数据融合下的森林防火本底数据库建立 . 浙江农林大学学报, 2008, 25(4): 405-410.
    [17] 王义平, 于振东, 吴鸿.  林木昆虫演变为重大害虫的主要环境因子 . 浙江农林大学学报, 2007, 24(6): 752-757.
    [18] 陈清林.  毛竹林林分因子与节肢动物群落的关系 . 浙江农林大学学报, 2006, 23(5): 565-569.
    [19] 李兰英, 高岚, 温亚利, 黄水灵, .  森林病害的环境影响分析 . 浙江农林大学学报, 2006, 23(5): 491-496.
    [20] 李土生, 杨幼平, 贾伟江, 翁卫松, 王同新, 华均健.  浙江省森林防火地理信息指挥系统的设计与开发 . 浙江农林大学学报, 2002, 19(3): 273-276.
  • 期刊类型引用(1)

    1. 童红卫,吴文娟,吕胤凯,王伟,吴初平,刘子琪,焦洁洁. 亚热带3种阔叶林造林模式植被层碳储量及其分配特征. 浙江林业科技. 2025(01): 8-15 . 百度学术

    其他类型引用(1)

  • 加载中
  • 链接本文:

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

    https://zlxb.zafu.edu.cn/article/zjnldxxb/2013/1/76

计量
  • 文章访问数:  3636
  • HTML全文浏览量:  437
  • PDF下载量:  1464
  • 被引次数: 2
出版历程
  • 收稿日期:  2012-02-16
  • 修回日期:  2012-06-03
  • 刊出日期:  2013-02-20

带约束条件的森林防火最优路径算法

doi: 10.11833/j.issn.2095-0756.2013.01.011
    通信作者: 徐爱俊

摘要: 以森林防火中的路径分析算法为研究对象,针对森林防火中对路径分析算法的特殊要求,以Dijkstra最短路径算法为基础,对它们进行改进与优化。结合林区交通网络自身的特点,引入道路宽度、最大限速、必经地和忽略地等森林防火路径分析约束因子,并以此作为道路权值确定的依据,提出了带约束条件的森林防火最优路径算法,并给出详细的算法设计与实现。该算法符合森林防火的实际情况,能够在森林防火应急指挥与决策中发挥作用。图2表2参17

English Abstract

王剑武, 季碧勇, 王铮屹, 等. 浙江省丽水市亚热带森林景观格局对森林碳密度的影响[J]. 浙江农林大学学报, 2024, 41(1): 30-40. DOI: 10.11833/j.issn.2095-0756.20230205
引用本文: 姜广宇, 徐爱俊, 黄小银, 等. 带约束条件的森林防火最优路径算法[J]. 浙江农林大学学报, 2013, 30(1): 76-82. DOI: 10.11833/j.issn.2095-0756.2013.01.011
WANG Jianwu, JI Biyong, WANG Zhengyi, et al. Impact of subtropical forest landscape pattern on forest carbon density in Lishui City of Zhejiang Province[J]. Journal of Zhejiang A&F University, 2024, 41(1): 30-40. DOI: 10.11833/j.issn.2095-0756.20230205
Citation: JIANG Guangyu, XU Aijun, HUANG Xiaoyin, et al. Optimal path algorithm for forest fire prevention with constraints[J]. Journal of Zhejiang A&F University, 2013, 30(1): 76-82. DOI: 10.11833/j.issn.2095-0756.2013.01.011

目录

/

返回文章
返回