本发明属于灾害区生态恢复评估,具体涉及一种历史地震单体滑坡植被恢复评估方法。
背景技术:
1、山区地震极易诱发大量滑坡,改变区内地形地貌,并对植被产生长期影响。由于历史地震诱发的古滑坡发生时间久远,无法获取到滑坡前及滑坡后一段时期区域的植被覆盖情况,所以无法采用现有滑坡区域植被恢复评估模型开展相关评估分析。目前,国内外进行滑坡区植被恢复评估的方法为通过以下评估模型计算植被恢复率vrr:
2、
3、其中,ndvi0为震前的ndvi值,ndvi1为震后的ndvi值,ndvi2为震后某一年的ndvi值,ndvi为归一化植被指数。
4、震后的植被恢复对滑坡区的修复至关重要,也能反映出滑坡的活动性,所以追踪滑坡区长时序的植被动态变化很有必要。到目前为止,世界范围内发生多起历史滑坡复活的事件,再次造成区域人员伤亡和经济损失。评估m年后滑坡区植被恢复状况,进而量化滑坡的活动性,能够及早规避历史滑坡复活造成的危害,同时也可以及时采取措施进行修复,进一步提升滑坡防灾减灾能力。
5、由于历史滑坡复活事件频发,受多种因素影响,又与植被覆盖紧密关联,因此,基于ndvi和fvc(植被覆盖度)数据评估滑坡区植被恢复率仍然是主要方法。基于此,本发明提供了一套科学合理可靠的历史地震单体滑坡植被恢复评估方法,用于计算历史地震单体滑坡的植被恢复率,量化滑坡的活动性,可预防滑坡的复活,进一步修复植被恢复率较低的滑坡区域。
技术实现思路
1、针对上述背景技术中存在的现有滑坡区植被恢复评估方法无法评估历史地震单体滑坡植被恢复的问题,本发明的目的是提供一种历史地震单体滑坡植被恢复评估方法,通过计算历史地震单体滑坡植被恢复率,定量评估滑坡的活动性,避免历史地震滑坡再次复活造成较大的人员伤亡和财产损失。
2、为了实现上述目的,本发明采用如下技术方案:
3、一种历史地震单体滑坡植被恢复评估方法,包括如下步骤:
4、步骤s1,选取历史地震单体滑坡体,设定历史地震单体滑坡植被恢复的评估时段;
5、步骤s2,在设定的评估时段m年中,计算评估时段内滑坡每个年份的植被覆盖度平均值fvci(i=1,2,……m):
6、
7、其中,ndvisoil为全裸区ndvi的值,ndviveg为全覆盖区ndvi的值,ndviveg采用滑坡区内ndvi的最大值ndvimax代替,ndvisoil采用滑坡区内ndvi的最小值ndvimin代替;
8、步骤s3,根据所述植被覆盖度平均值绘制评估时段内单体滑坡累积植被覆盖度曲线,根据曲线形态和数据拐点确定评估历史地震单体滑坡植被恢复的基准期和评估期;
9、步骤s4,分别计算基准期的植被覆盖度fvcb和评估期的植被覆盖度fvct;
10、步骤s5,根据fvcb和fvct建立历史地震单体滑坡植被恢复评估模型:
11、
12、计算得到历史地震单体滑坡评估期的植被恢复率vrr;根据vrr值划分植被恢复等级及滑坡活动性类型。
13、进一步的,步骤s3中,所述基准期和评估期的确定方法为:
14、情况1:评估时段内单体滑坡累积植被覆盖度曲线的形态起伏较大并且有明显拐点时,确定共n个拐点,选取每个拐点gj(j=1,2,……n,n≤m-2)为分割点,将整个评估时段共分为k期(k≤1,2,……n+1);第1期即为该评估的基准期,其后的任意一期为评估期;
15、情况2:评估时段内单体滑坡累积植被覆盖度曲线的形态较平滑顺直并且无明显拐点时,根据评估时段及研究特性自主划分区间,评估时段内的第1期为基准期,其后的任意一期为评估期。
16、进一步的,针对上述2种情况,所述基准期的植被覆盖度fvcb和评估期的植被覆盖度fvct计算方式如下:
17、针对情况1:为分期计算,基准期和评估期的植被覆盖度通过fvck的计算公式得到为:
18、或
19、其中,fvcl为累积植被覆盖度曲线上第k期内前一拐点对应的值,fvcl+1为累积植被覆盖度曲线上第k期内后一拐点对应的值,p为该期的年份数;fvcx为第k期内第x年份的植被覆盖度平均值;
20、针对情况2:为自主分区计算,基准期和评估期的植被覆盖度选取区间内的均值进行计算。
21、进一步的,步骤s5中,所述植被恢复等级及滑坡活动性类型划分如下:
22、
23、相比于现有技术的缺点和不足,本发明具有以下有益效果:
24、本发明能够计算历史滑坡m年后的植被恢复率,克服了由于历史地震滑坡前及滑坡后一段时期内的植被数据的缺失而无法计算vrr的问题。本发明结合基准期和评估期的数据计算评估期的植被恢复率,得到的结果与实际植被恢复情况贴近,能够精准反映历史滑坡m年后的植被恢复情况,可进一步量化滑坡的活动性。
1.一种历史地震单体滑坡植被恢复评估方法,其特征在于,包括如下步骤:
2.如权利要求1所述的历史地震单体滑坡植被恢复评估方法,其特征在于,步骤s3中,所述基准期和评估期的确定方法为:
3.如权利要求2所述的历史地震单体滑坡植被恢复评估方法,其特征在于,所述基准期的植被覆盖度fvcb和评估期的植被覆盖度fvct计算方式如下:
4.如权利要求1所述的历史地震单体滑坡植被恢复评估方法,其特征在于,步骤s5中,所述植被恢复等级及滑坡活动性类型划分如下: