跨时间尺度联合的古滑坡复活风险分级评估方法和装置

1.本发明一般涉及地理空间数据处理领域，并且更具体地，涉及一种跨时间尺度联合的古滑坡复活风险分级评估方法和装置。


背景技术：

2.古滑坡是斜坡长期复杂演化过程的产物，物质组成复杂、岩土力学性质特殊，具有极强的隐蔽性和扰动敏感性。古滑坡泛指曾历经一次或多次滑动，在受到外界因素(强降雨、库水位变化)扰动时会再次滑动或具有滑动趋势，不稳定的古滑坡是人类工程活动的潜灾体。我国西南地区因其特殊的地理位置、地形地貌特征和极端气候条件等因素，致使区域内古滑坡失稳及复活所造成的地质灾害日益频发。西南地区地理地质环境(即长时间尺度的因素)与库水、降雨等内外因素耦合诱发的古滑坡失稳复活，其发育、变形影响因素众多、形变破坏机理复杂，随着大量基础设施的建设和外界极端因素(如强降雨、库水位变化等)的影响，滑坡地质灾害呈现日趋严重的趋势。
3.由于古滑坡复活灾害带来的广泛性、威胁性、损失性等恶劣影响，其失稳复活风险精准定量评估如今已成为灾前科学预警和主动防灾减灾的必要手段和迫切需求。野外勘测已无法满足对古滑坡复活灾害隐患风险定量评估的需求，利用数据驱动模型进行滑坡易发性分析成为主流。其中，基于机器学习方法的滑坡易发性分析非线性关系数据处理能力强，结合sbas-insar技术得到的地表形变速率数据，能够跟踪古滑坡失稳复活全生命周期的形变破坏细部信息，寻找二次灾害与灾害链的临界触发点，再进行短时间尺度的动力学模型局部微观分析，提高古滑坡失稳复活隐患风险评估的精准性与可靠性及科学的监测预警能力。
4.使用传统的机器学习方法依赖于大量有代表性的数据样本，处理小样本数据易导致风险评估分析结果的欠拟合和过度拟合的问题，同时容易忽略古滑坡复活受外界极端因素(如强降雨、地震等)的影响，评估未与短期动力学触发过程建立动态关联，无法应对诱灾动力学触发因素突发性强、短期触发过程动态随机等情况，且古滑坡复活灾害时空演变信息量较为片面滞后、时效性不足，导致评估结果的随机性、偏移性和盲目性难以避免。
5.在长时间尺度古滑坡发育演化过程中对潜在失稳复活特征机理信息进行深入挖掘的同时，古滑坡复活灾害还面临着诱灾动力学触发因素突发性强、短期触发过程动态随机等问题。目前，对于古滑坡滑坡的复活触发过程，尤其是动力学过程研究关注较少。风险评估难以与短期动力学触发过程建立动态关联，以至于机器学习无法准确评估复活风险等级，面向海量待预测潜灾体无法事先精准评估滑坡灾害隐患并实施科学的监测预警措施，最终造成防灾减灾预警和应急时效性不足，损失惨重等痛点问题。


技术实现要素：

6.根据本发明的实施例，提供了一种跨时间尺度联合的古滑坡复活风险分级评估方案。本方案能够加强古滑坡发育演化机理与动力学触发因素在时域上的动态关联，解决单
一尺度评估精准性与可靠性不足的问题，实现孕灾场景复杂以及跨时间尺度内外因素耦合关联约束的古滑坡复活风险分级评估。
7.在本发明的第一方面，提供了一种跨时间尺度联合的古滑坡复活风险分级评估方法。该方法包括：
8.获取古滑坡区域历史观测数据，构建古滑坡复活灾害特征评价因子体系，计算古滑坡长期发育演进过程复活灾害发生的时空概率，根据所述时空概率划分长时间尺度古滑坡复活灾害的易发性等级；
9.获取sentinel-1影像数据集，通过sbas-insar方法计算出沿雷达视线方向的地表形变速率；对所述沿雷达视线方向的地表形变速率进行二维形变转换，得到沿坡度方向的地表形变速率，并划分地表形变速率等级；以所述地表形变速率等级对应更新所述长时间尺度古滑坡复活灾害的易发性等级，得到长时间尺度古滑坡复活风险性等级；所述长时间尺度古滑坡复活风险性等级包括高风险等级；
10.将所述高风险等级对应的区域作为古滑坡发育演化的关键节点，在所述关键节点上分析诱发古滑坡复活的动力学触发因素，构建短时间尺度古滑坡复活的动力学方程，并根据短时间尺度古滑坡复活的动力学方程计算出二次灾害下的滑带土残余抗剪强度，更新所述长时间尺度古滑坡复活风险性等级，得到古滑坡复活风险分级评估模型；
11.将待分析的古滑坡区域历史观测数据和sentinel-1影像数据集输入所述古滑坡复活风险分级评估模型，输出古滑坡复活风险分级评估结果。
12.进一步地，所述构建古滑坡复活灾害特征评价因子体系，包括：
13.从所述古滑坡区域历史观测数据中提取深部地质条件特征因子与地形地貌条件特征因子；
14.对所述深部地质条件特征因子与地形地貌条件特征因子进行归一化；
15.计算归一化后的特征因子的滑坡频率比和信息熵权值。
16.进一步地，所述通过sbas-insar方法计算出沿雷达视线方向的地表形变速率，包括：
17.对所述sentinel-1影像数据集进行干涉像对的配对，得到若干对配对后的干涉像对；
18.利用所述若干对配对后的干涉像对生成干涉图，对所述干涉图进行去平处理，对去平后的干涉图进行自适应滤波，并生成相干系数，得到优化后的干涉图；
19.将所述优化后的干涉图进行sbas第一次反演，得到第一次估算的形变速率和残余地形；再对第一次估算的形变速率和残余地形进行sbas第二次反演，得到时间序列上的位移；再对所述时间序列上的位移进行地理编码，得到沿雷达视线方向的地表形变速率。
20.进一步地，所述对所述雷达视线方向的地表形变速率进行二维形变转换，得到沿坡度方向的地表形变速率，包括：
[0021][0022][0022]
其中，v
slope
为沿坡度方向的地表形变速率；v
los
为沿雷达视线方向的地表形变速
率；αs为方位向和正北方向的夹角；α为斜坡坡向；β为视线与斜坡的夹角；θ为入射角；为斜坡坡度。
[0023]
其中，v
slope
为沿坡度方向的地表形变速率；v
los
为沿雷达视线方向的地表形变速率；αs为方位向和正北方向的夹角；α为斜坡坡向；β为视线与斜坡的夹角；θ为入射角；为斜坡坡度。
[0024]
进一步地，所述以所述地表形变速率等级对应更新所述长时间尺度古滑坡复活灾害的易发性等级，包括：
[0025]
若目标古滑坡区域对应地表形变速率等级的级别高于所述目标古滑坡区域对应古滑坡复活灾害的易发性等级的级别，则将所述目标古滑坡区域对应古滑坡复活灾害的易发性等级的级别更新为所述目标古滑坡区域对应地表形变速率等级的级别；将更新后的古滑坡复活灾害的易发性等级作为长时间尺度古滑坡复活风险性等级；
[0026]
所述地表形变速率等级与古滑坡复活灾害的易发性等级的级别数相同。
[0027]
进一步地，所述根据短时间尺度古滑坡复活的动力学方程计算出二次灾害下的滑带土残余抗剪强度，包括：
[0028]
在古滑坡失稳过程中，当滑体处于间歇期和滑动期的临界触发状态时，计算滑带土的抗剪强度：
[0029]
τf＝c
′
re c
+{σ-[u0+p(zm，t
diff
)]}tanφ
′
re c
[0030]
其中，τf为滑带土的抗剪强度；c
′
re c
为经强度恢复后滑带土的有效黏聚力；σ为总应力；u0为滑带土的初始孔隙水压力；p(zm,t
diff
)为降雨或库水位变动产生的水头压力扩散到滑动面处引发的滑带土孔隙的水压力增量,zm为滑坡体厚度；t
diff
为降雨或库水位变动产生的水头压力扩散到滑动面的时间；
[0031]
选择所述滑体中单位宽度的滑坡块体作为研究对象，在所述在古滑坡失稳过程中，所述滑坡块体不发生相对位移且运动速率大小和方向均相同，当所述滑坡块体处于间歇期和滑动期的临界触发状态时，计算下滑力τ：
[0032][0033]
当所述滑坡块体进入滑动期，开始沿滑动面滑动时，获得初始加速度a0，计算滑动期滑带土的残余抗剪强度：
[0034][0035]
其中，τ为下滑力；为所述滑动块体的平均重度；g为重力加速度；a0为突破临界状态后滑坡块体的初始加速度；a为滑动期滑坡块体的加速度；τv为滑动期滑带土的残余抗剪强度。
[0036]
进一步地，所述更新所述长时间尺度古滑坡复活风险性等级，包括：
[0037]
若所述滑带土残余抗剪强度的剪切速率效应为正速率效应，则所述关键节点对应的长时间尺度古滑坡复活风险性等级的级别与所述滑带土残余抗剪强度的剪切速率效应呈正相关，并根据正相关的相关程度分级更新所述关键节点对应的长时间尺度古滑坡复活
风险等级的级别；
[0038]
其中，所述滑带土残余抗剪强度的剪切速率效应为正速率效应，则：
[0039]
h＝τ
v-τr＞0
[0040]
h为剪切速率效应系数，表示正相关的相关程度，若h越大，则正相关的相关程度越高；τv为滑动期滑带土的残余抗剪强度；τr为间歇期滑带土的残余抗剪强度。
[0041]
在本发明的第二方面，提供了一种跨时间尺度联合的古滑坡复活风险分级评估装置。该装置包括：
[0042]
第一分级模块，用于获取古滑坡区域历史观测数据，构建古滑坡复活灾害特征评价因子体系，计算古滑坡长期发育演进过程复活灾害发生的时空概率，根据所述时空概率划分长时间尺度古滑坡复活灾害的易发性等级；
[0043]
第二分级模块，用于获取sentinel-1影像数据集，通过sbas-insar方法计算出沿雷达视线方向的地表形变速率；对所述沿雷达视线方向的地表形变速率进行二维形变转换，得到沿坡度方向的地表形变速率，并划分地表形变速率等级；以所述地表形变速率等级对应更新所述长时间尺度古滑坡复活灾害的易发性等级，得到长时间尺度古滑坡复活风险性等级；所述长时间尺度古滑坡复活风险性等级包括高风险等级；
[0044]
等级更新模块，用于将所述高风险等级对应的区域作为古滑坡发育演化的关键节点，在所述关键节点上分析诱发古滑坡复活的动力学触发因素，构建短时间尺度古滑坡复活的动力学方程，并根据短时间尺度古滑坡复活的动力学方程计算出二次灾害下的滑带土残余抗剪强度，更新所述长时间尺度古滑坡复活风险性等级，得到古滑坡复活风险分级评估模型；
[0045]
评估模块，用于将待分析的古滑坡区域历史观测数据和sentinel-1影像数据集输入所述古滑坡复活风险分级评估模型，输出古滑坡复活风险分级评估结果。
[0046]
在本发明的第三方面，提供了一种电子设备。该电子设备至少一个处理器；以及与所述至少一个处理器通信连接的存储器；所述存储器存储有可被所述至少一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器能够执行本发明第一方面的方法。
[0047]
在本发明的第四方面，提供了一种存储有计算机指令的非瞬时计算机可读存储介质，所述计算机指令用于使所述计算机执行本发明第一方面的方法。
[0048]
应当理解，发明内容部分中所描述的内容并非旨在限定本发明的实施例的关键或重要特征，亦非用于限制本发明的范围。本发明的其它特征将通过以下的描述变得容易理解。
附图说明
[0049]
结合附图并参考以下详细说明，本发明各实施例的上述和其他特征、优点及方面将变得更加明显。在附图中，相同或相似的附图标记表示相同或相似的元素，其中：
[0050]
图1示出了根据本发明的实施例的跨时间尺度联合的古滑坡复活风险分级评估方法的流程图；
[0051]
图2示出了根据本发明的实施例的跨时间尺度联合的古滑坡复活风险分级评估装置的方框图；
[0052]
图3示出了能够实施本发明的实施例的示例性电子设备的方框图；
[0053]
其中，300为电子设备、301为cpu、302为rom、303为ram、304为总线、305为i/o接口、306为输入单元、307为输出单元、308为存储单元、309为通信单元。
具体实施方式
[0054]
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的全部其他实施例，都属于本发明保护的范围。
[0055]
另外，本文中术语“和/或”，仅仅是一种描述关联对象的关联关系，表示可以存在三种关系，例如，a和/或b，可以表示：单独存在a，同时存在a和b，单独存在b这三种情况。另外，本文中字符“/”，一般表示前后关联对象是一种“或”的关系。
[0056]
本发明中，能够加强古滑坡发育演化机理与动力学触发因素在时域上的动态关联，解决单一尺度评估精准性与可靠性不足的问题，实现孕灾场景复杂以及跨时间尺度内外因素耦合关联约束的古滑坡复活风险分级评估。
[0057]
图1示出了本发明实施例的跨时间尺度联合的古滑坡复活风险分级评估方法的流程图。
[0058]
该方法包括：
[0059]
s101、获取古滑坡区域历史观测数据，构建古滑坡复活灾害特征评价因子体系，计算古滑坡长期发育演进过程复活灾害发生的时空概率，根据所述时空概率划分长时间尺度古滑坡复活灾害的易发性等级。
[0060]
作为本发明的一种实施例，古滑坡是曾历经一次或多次滑动，在受到外界因素，例如强降雨、库水位变化等扰动时会再次滑动或具有滑动趋势的地质环境。我国西南地区因其特殊的地理位置、地形地貌特征和极端气候条件等因素，致使区域内古滑坡失稳及复活所造成的地质灾害日益频发。西南地区地理地质环境(即长时间尺度的因素)与库水、降雨等内外因素耦合诱发的古滑坡失稳复活，其发育、变形影响因素众多、形变破坏机理复杂，随着大量基础设施的建设和外界极端因素(如强降雨、库水位变化等)的影响，滑坡地质灾害呈现日趋严重的趋势，作为本实施例的研究对象。
[0061]
作为本发明的一种实施例，所述构建古滑坡复活灾害特征评价因子体系，包括：
[0062]
首先，从所述古滑坡区域历史观测数据中提取深部地质条件特征因子与地形地貌条件特征因子。所述深部地质条件特征因子，例如断层、岩性、土壤侵蚀、裂缝发育情况等。所述地形地貌条件特征因子，例如坡度坡向、剖面曲率、植被覆盖等。
[0063]
其次，对所述深部地质条件特征因子与地形地貌条件特征因子进行归一化。
[0064]
古滑坡复活灾害特征评价因子呈现多源异构的特点，其中存在连续型和离散型两种类型，需要进行空间化预处理，对离散型评价因子进行分类量化，对连续型评价因子采取分级离散化处理，统一模型输入。
[0065]
在本实施例中，归一化处理如下：
[0066]
[0067]
其中，y
*
为归一化处理前的样本数据；y为归一化处理前的样本数据；max为样本数据的最大值，min为样本数据的最小值。
[0068]
综合提取形成各评价因子分级量化，并生成滑坡各评价因子专题图。
[0069]
最后，计算归一化后的特征因子的滑坡频率比和信息熵权值。
[0070][0071]
其中，ni、si分别表示分布在评价因子ui中某一分级类别中的古滑坡单元数和分级区域面积；n、s为研究区域古滑坡单元总数和研究区总面积；fr
ij
、p
ij
分别为第i个评价因子指标下第j个分级类别的滑坡频率比和分布概率，fr
ij
＞1表明因子该类别的状态条件利于古滑坡复活发生；fr
ij
＜1表明该状态条件不利于古滑坡复活发生；k为常数项，由每个古滑坡灾害环境评价因子的分级类别数ki取值决定；wei为熵权值，ei为信息熵，熵权值与信息熵呈反比，熵权值越大，特征评价因子数据中的信息量不确定性越小，该评价因子对古滑坡复活灾害发育贡献度越大。
[0072]
通过计算滑坡频率比及信息熵权值，量化古滑坡复活灾害特征评价因子对古滑坡灾害环境特征规律及贡献度的影响，结合古滑坡知识机理规律与研究区基础地理概况综合归纳区域古滑坡发育演变的作用机制和空间统计规律。
[0073]
作为本发明的一种实施例，所述计算古滑坡长期发育演进过程复活灾害发生的时空概率，根据所述时空概率划分长时间尺度古滑坡复活灾害的易发性等级，包括：
[0074]
以所述归一化且分级量化的特征评价因子研究区历史专题数据为输入，结合bagging和随机选择特征建立并融合多个决策树构建随机森林模型，分级刻画古滑坡长期发育演进过程复活灾害发生的时空概率。所述随机森林模型(random forest)指的是利用多棵树对样本进行训练并预测的一种分类器；在机器学习中，随机森林是一个包含多个决策树的分类器，并且其输出的类别是由个别树输出的类别的众数而定。
[0075]
在本实施例中，根据所述时空概率划分长时间尺度古滑坡复活灾害的易发性等级可以划分为多个等级，例如划分为五个等级，分别为极低易发性、低易发性、中易发性、高易发性、极高易发性。
[0076]
作为本发明的一种实施例，在划分易发性等级后，利用自然断点法将连续性古滑坡易发性图层对应重分类为5个离散的古滑坡易发性等级分区图。
[0077]
s102、获取sentinel-1影像数据集，通过sbas-insar方法计算出沿雷达视线方向的地表形变速率；对所述沿雷达视线方向的地表形变速率进行二维形变转换，得到沿坡度
方向的地表形变速率，并划分地表形变速率等级；以所述地表形变速率等级对应更新所述长时间尺度古滑坡复活灾害的易发性等级，得到长时间尺度古滑坡复活风险性等级；所述长时间尺度古滑坡复活风险性等级包括高风险等级。
[0078]
sentinel-1影像数据集为通过sentinel-1卫星采集到的影像数据。sentinel-1卫星是欧洲航天局哥白尼计划(gmes)中的地球观测卫星，由两颗卫星组成，载有c波段合成孔径雷达，可提供连续图像(白天、夜晚和各种天气)。
[0079]
作为本发明的一种实施例，所述通过sbas-insar方法计算出沿雷达视线方向的地表形变速率，包括：
[0080]
首先，对所述sentinel-1影像数据集进行干涉像对的配对，得到若干对配对后的干涉像对。
[0081]
在本实施例中，对sentinel-1影像数据集进行干涉像对的配对，像对结果图表方式输出，输入n景数据，能得到的最大配对数是(n*(n-1))/2，生成连接图工具会选择最优的组合方式进行配对。程序自动选择超级主影像，在整个处理中，超级主影像作为参考影像，所有的像对都会配准到超级主影像上。
[0082]
其次，利用所述若干对配对后的干涉像对生成干涉图，对所述干涉图进行去平处理，对去平后的干涉图进行自适应滤波，得到优化后的干涉图。
[0083]
具体地，对所有的配对的干涉像对进行干涉处理，生成干涉图，需要设置多视的视数，从传感器的原始采样间隔和入射角计算得出，能让最后的结果尽可能的得到方位向和距离向一致的分辨率。注：视数大于1：1。
[0084]
具体地，去平处理：去除间隔相同，由平地效应造成的干涉条纹；保留因地形起伏、高程起伏所导致的干涉条纹。
[0085]
具体地，自适应滤波：对上一步去平后的干涉图进行滤波，去掉由平地干涉引起的位相噪声，让干涉条纹变得更加平滑。同时，生成干涉的相干图(描述位相质量)和滤波后的主影像强度图,增加干涉图的信噪比，提供了更可靠的相干性，为sbas反演做好数据准备。
[0086]
最后，将所述优化后的干涉图进行sbas第一次反演，得到第一次估算的形变速率和残余地形。所述sbas第反演，使用线性模型，在第一次估算速率图中，可以看到明显的沉降和抬升，sbas反演的下一步就会得到优化的速率估算值，在这一步中确定了形变的位置和空间范围。再对第一次估算的形变速率和残余地形进行sbas第二次反演，得到时间序列上的位移。sbas第二次反演是在第一次估算形变速率之后，在这一步用滤波的方法去除大气相位，得到最终的时序上的形变量。大气滤波起到了对时间序列上的形变平滑的作用，大气的去除是通过低通滤波和高通滤波实现的。这一步的最后，会根据第一次反演选择的形变模型计算出最终的形变速率，是根据多项式计算出来的纯净的时序上的形变量。这里还会生成模型计算出的平方值。再对所述时间序列上的位移进行地理编码，得到沿雷达视线方向的地表形变速率。其中，形变结果和形变速率会重投影到用户定义的方向，此处已知直接投影到沿雷达视线方向。
[0087]
作为本发明的一种实施例，由于滑坡多沿斜坡面进行滑动，雷达视线方向的形变信息无法准确反映斜坡面的真实形变情况，考虑雷达视线方向和坡度方向的几何关系，假设运动沿单位矢量指定的方向发生，采用下式将视线方向的形变速率转化为坡度方向的形变速率。
[0088][0089]vslope
＝v
los
/cosβ
[0090]
上式中，cosβ表示为：
[0091][0092]
其中，v
slope
为沿坡度方向的地表形变速率；v
los
为沿雷达视线方向的地表形变速率；αs为方位向和正北方向的夹角；α为斜坡坡向；β为视线与斜坡的夹角；θ为入射角；为斜坡坡度。
[0093]
通过sbas-insar方法计算出沿雷达视线方向的地表形变速率；对所述沿雷达视线方向的地表形变速率进行二维形变转换，得到沿坡度方向的地表形变速率，旨在减少单一使用机器学习模型可能出现的假阴性数，精准更新易发性等级，进一步建立深度精准复活风险分级标准，从而宏观确定古滑坡的危险性等级及危险区域范围。
[0094]
在本实施例中，经sbas-insar方法处理得到实际地表形变速率数据，根据地表形变速率阈值范围，与所述长时间尺度古滑坡复活灾害的易发性等级划分相同的等级，即同样划分五个速率等级。
[0095]
作为本发明的一种实施例，以所述地表形变速率等级对应更新所述长时间尺度古滑坡复活灾害的易发性等级，得到长时间尺度古滑坡复活风险性等级，具体包括：
[0096]
若目标古滑坡区域对应地表形变速率等级的级别高于所述目标古滑坡区域对应古滑坡复活灾害的易发性等级的级别，则将所述目标古滑坡区域对应古滑坡复活灾害的易发性等级的级别更新为所述目标古滑坡区域对应地表形变速率等级的级别；将更新后的古滑坡复活灾害的易发性等级作为长时间尺度古滑坡复活风险性等级。
[0097]
在本实施例中，如表1表示将古滑坡复活灾害的易发性等级更新为长时间尺度古滑坡复活风险性等级：
[0098][0099]
表1
[0100]
表1中，列表示将地表形变速率等级划分为5个级别，行表示将易发性等级划分为5个级别，表中元素表示更新后的长时间尺度古滑坡复活风险性等级。例如，第一行第二列表
示易发性等级为1级，地表形变速率等级为2级时，由于地表形变速率等级的级别高于易发性等级的级别，则需要将易发性等级的级别提升1级，作为对应风险性等级的级别。
[0101]
所述长时间尺度古滑坡复活风险性等级与所述地表形变速率等级以及古滑坡复活灾害的易发性等级的级别数均相同。例如，在本实施例中，所述长时间尺度古滑坡复活风险性等级划分为5级，分别为极低风险级别、低风险级别、中风险级别、高风险级别和极高风险级别。将高风险级别和极高风险级别作为高风险等级。即高风险等级为古滑坡复活风险较高的级别。
[0102]
s103、将所述高风险等级对应的区域作为古滑坡发育演化的关键节点，在所述关键节点上分析诱发古滑坡复活的动力学触发因素，构建短时间尺度古滑坡复活的动力学方程，并根据短时间尺度古滑坡复活的动力学方程计算出二次灾害下的滑带土残余抗剪强度，更新所述长时间尺度古滑坡复活风险性等级，得到古滑坡复活风险分级评估模型。
[0103]
在一些实施例中，诱发古滑坡复活的动力学触发因素，例如季节性强降雨导致的降雨入渗引起的周期性的库水位变化。
[0104]
在一些实施例中，构建短时间尺度古滑坡复活的动力学方程，将降雨和库水位变动在滑坡边界产生的作用可简化为具有某一水头的压力，该压力以水力扩散的形式影响滑坡各部位的孔隙水压力。假设滑坡体物质是均质的土体，则根据均匀多孔介质中的一维竖向瞬时达西流的规律特征，降雨或库水位变动产生的水头压力在滑坡体中的扩散过程可用下式表示：
[0105][0106]
其中，p(z,t)为滑动块体深度z位置在t时刻的孔隙水压力增量；d为滑坡体的水力扩散系数。
[0107][0108]
其中，τr为间歇期滑带土的残余抗剪强度；c
′r为滑带土的有效残余黏聚力；σ为总应力；u0为滑带土的初始孔隙水压力；为有效残余内摩擦角。
[0109]
当滑体处于间歇期和滑动期的临界状态时，滑带土的抗剪强度τf可表示为：
[0110]
τf＝c
′
re c
+{σ-[u0+p(zm，t
diff
)]}tanφ
′
re c
[0111]
其中，τf为滑带土的抗剪强度；c
′
re c
为经强度恢复后滑带土的有效黏聚力；σ为总应力；u0为滑带土的初始孔隙水压力；p(zm,t
diff
)为降雨或库水位变动产生的水头压力扩散到滑动面处引发的滑带土孔隙的水压力增量,zm为滑坡体厚度；t
diff
为降雨或库水位变动产生的水头压力扩散到滑动面的时间。
[0112]
取单位宽度滑坡块体为研究对象，在间歇性复活的缓慢滑动过程中，块体本身不发生相对位移，那么在滑坡块体上任一点(包括滑动面处)的运动速率大小和方向均相同。假定滑坡块体为刚体，其受到上部块体的推力和下部块体的阻力大小相等且方向相反，因此，当滑坡块体处于间歇期和滑动期的临界状态时，有：
[0113]
[0114]
其中，τ为下滑力；为所述滑动块体的平均重度；g为重力加速度；a0为突破临界状态后滑坡块体的初始加速度；a为滑动期滑坡块体的加速度；τv为滑动期滑带土的残余抗剪强度；m为单位宽度滑坡块体的质量；s为单位宽度滑坡块体的底面积。
[0115]
在获得初始加速度a0后，进入滑动期的滑坡体沿滑动面开始滑动，此时有：
[0116][0117]
其中，f(v)为滑坡运动速率v的函数。
[0118]
在一些实施例中，所述更新所述长时间尺度古滑坡复活风险性等级，包括：
[0119]
若所述滑带土残余抗剪强度的剪切速率效应为正速率效应，则所述关键节点对应的长时间尺度古滑坡复活风险性等级的级别与所述滑带土残余抗剪强度的剪切速率效应呈正相关，并根据正相关的相关程度更新所述关键节点对应的长时间尺度古滑坡复活风险等级的级别。
[0120]
其中，所述滑带土残余抗剪强度的剪切速率效应为正速率效应，则：
[0121]
h＝τ
v-τr＞0
[0122]
h为剪切速率效应系数，表示正相关的相关程度，若h越大，则正相关的相关程度越高；τv为滑动期滑带土的残余抗剪强度；τr为间歇期滑带土的残余抗剪强度。
[0123]
具体地，当τv＞τr时，表明滑带土的残余强度具有正速率效应，即随剪切速率的增大而增大；反之，当τv＜τr时，滑带土的残余强度具有负速率效应；而当τv＝τr时，则无效应。
[0124]
所述根据正相关的相关程度更新所述关键节点对应的长时间尺度古滑坡复活风险等级的级别，具体包括：
[0125]
若剪切速率效应系数高于高风险系数阈值范围，则对应提升所述关键节点对应的长时间尺度古滑坡复活风险等级的级别至相应最高风险等级。
[0126]
综上，即为描述库岸古滑坡短时间尺度触发复活的动力学方程。
[0127]
s104、将待分析的古滑坡区域历史观测数据和sentinel-1影像数据集输入所述古滑坡复活风险分级评估模型，输出古滑坡复活风险分级评估结果。
[0128]
根据本发明的实施例，以短时间尺度动力学触发因素与长时间尺度古滑坡发育进化机理定量关联的联合机制为核心，克服单一尺度对古滑坡复活风险进行分析研究，在系统分析古滑坡隐式成因机理和显式表征的基础上，深度分析短期动力学触发因素和长期滑坡形变位移演化过程存在的典型时空关联性特征，寻找二次灾害与灾害链的临界触发点，建立短时间尺度动力学触发因素及长时间尺度古滑坡发育演化机理定量关联的联合响应机制。在机器学习框架下，建立古滑坡演化机理模型和动力学触发模型之间的高层映射，多尺度多阶段实时动态风险跟踪，在易发性分析基础之上实现跨尺度联合的古滑坡失稳复活风险精准定量评估。最终，够实现孕灾场景复杂，跨尺度内外因素耦合关联约束的古滑坡复活风险定量测算，能够加强古滑坡发育演化机理与动力学触发因素在时域上的动态定量关联，解决单一尺度评估精准性与可靠性不足的问题，能够支持多尺度多阶段内外因素耦合关联约束的风险定量测算，做到事先精准评估古滑坡复活隐患并实施科学有效的监测预警。
[0129]
需要说明的是，对于前述的各方法实施例，为了简单描述，故将其都表述为一系列
的动作组合，但是本领域技术人员应该知悉，本发明并不受所描述的动作顺序的限制，因为依据本发明，某些步骤可以采用其他顺序或者同时进行。其次，本领域技术人员也应该知悉，说明书中所描述的实施例均属于可选实施例，所涉及的动作和模块并不一定是本发明所必须的。
[0130]
以上是关于方法实施例的介绍，以下通过装置实施例，对本发明所述方案进行进一步说明。
[0131]
如图2所示，装置200包括：
[0132]
第一分级模块210，用于获取古滑坡区域历史观测数据，构建古滑坡复活灾害特征评价因子体系，计算古滑坡长期发育演进过程复活灾害发生的时空概率，根据所述时空概率划分长时间尺度古滑坡复活灾害的易发性等级。
[0133]
第二分级模块220，用于获取sentinel-1影像数据集，通过sbas-insar方法计算出沿雷达视线方向的地表形变速率；对所述沿雷达视线方向的地表形变速率进行二维形变转换，得到沿坡度方向的地表形变速率，并划分地表形变速率等级；以所述地表形变速率等级对应更新所述长时间尺度古滑坡复活灾害的易发性等级，得到长时间尺度古滑坡复活风险性等级；所述长时间尺度古滑坡复活风险性等级包括高风险等级。
[0134]
等级更新模块230，用于将所述高风险等级对应的区域作为古滑坡发育演化的关键节点，在所述关键节点上分析诱发古滑坡复活的动力学触发因素，构建短时间尺度古滑坡复活的动力学方程，并根据短时间尺度古滑坡复活的动力学方程计算出二次灾害下的滑带土残余抗剪强度，更新所述长时间尺度古滑坡复活风险性等级，得到古滑坡复活风险分级评估模型。
[0135]
评估模块240，用于将待分析的古滑坡区域历史观测数据和sentinel-1影像数据集输入所述古滑坡复活风险分级评估模型，输出古滑坡复活风险分级评估结果。
[0136]
所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，所述描述的模块的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
[0137]
本发明的技术方案中，所涉及的用户个人信息的获取，存储和应用等，均符合相关法律法规的规定，且不违背公序良俗。
[0138]
根据本发明的实施例，本发明还提供了一种电子设备和一种可读存储介质。
[0139]
图3示出了可以用来实施本发明的实施例的电子设备300的示意性框图。电子设备旨在表示各种形式的数字计算机，诸如，膝上型计算机、台式计算机、工作台、个人数字助理、服务器、刀片式服务器、大型计算机、和其它适合的计算机。电子设备还可以表示各种形式的移动装置，诸如，个人数字处理、蜂窝电话、智能电话、可穿戴设备和其它类似的计算装置。本文所示的部件、它们的连接和关系、以及它们的功能仅仅作为示例，并且不意在限制本文中描述的和/或者要求的本发明的实现。
[0140]
设备300包括计算单元301，其可以根据存储在只读存储器(rom)302中的计算机程序或者从存储单元308加载到随机访问存储器(ram)303中的计算机程序，来执行各种适当的动作和处理。在ram 303中，还可存储设备300操作所需的各种程序和数据。计算单元301、rom 302以及ram 303通过总线304彼此相连。输入/输出(i/o)接口305也连接至总线304。
[0141]
设备300中的多个部件连接至i/o接口305，包括：输入单元306，例如键盘、鼠标等；输出单元307，例如各种类型的显示器、扬声器等；存储单元308，例如磁盘、光盘等；以及通
信单元309，例如网卡、调制解调器、无线通信收发机等。通信单元309允许设备300通过诸如因特网的计算机网络和/或各种电信网络与其他设备交换信息/数据。
[0142]
计算单元301可以是各种具有处理和计算能力的通用和/或专用处理组件。计算单元301的一些示例包括但不限于中央处理单元(cpu)、图形处理单元(gpu)、各种专用的人工智能(ai)计算芯片、各种运行机器学习模型算法的计算单元、数字信号处理器(dsp)、以及任何适当的处理器、控制器、微控制器等。计算单元301执行上文所描述的各个方法和处理，例如方法s101～s104。例如，在一些实施例中，方法s101～s104可被实现为计算机软件程序，其被有形地包含于机器可读介质，例如存储单元308。在一些实施例中，计算机程序的部分或者全部可以经由rom 302和/或通信单元309而被载入和/或安装到设备300上。当计算机程序加载到ram 303并由计算单元301执行时，可以执行上文描述的方法s101～s104的一个或多个步骤。备选地，在其他实施例中，计算单元301可以通过其他任何适当的方式(例如，借助于固件)而被配置为执行方法s101～s104。
[0143]
本文中以上描述的系统和技术的各种实施方式可以在数字电子电路系统、集成电路系统、场可编程门阵列(fpga)、专用集成电路(asic)、专用标准产品(assp)、芯片上系统的系统(soc)、负载可编程逻辑设备(cpld)、计算机硬件、固件、软件、和/或它们的组合中实现。这些各种实施方式可以包括：实施在一个或者多个计算机程序中，该一个或者多个计算机程序可在包括至少一个可编程处理器的可编程系统上执行和/或解释，该可编程处理器可以是专用或者通用可编程处理器，可以从存储系统、至少一个输入装置、和至少一个输出装置接收数据和指令，并且将数据和指令传输至该存储系统、该至少一个输入装置、和该至少一个输出装置。
[0144]
用于实施本发明的方法的程序代码可以采用一个或多个编程语言的任何组合来编写。这些程序代码可以提供给通用计算机、专用计算机或其他可编程数据处理装置的处理器或控制器，使得程序代码当由处理器或控制器执行时使流程图和/或框图中所规定的功能/操作被实施。程序代码可以完全在机器上执行、部分地在机器上执行，作为独立软件包部分地在机器上执行且部分地在远程机器上执行或完全在远程机器或服务器上执行。
[0145]
在本发明的上下文中，机器可读介质可以是有形的介质，其可以包含或存储以供指令执行系统、装置或设备使用或与指令执行系统、装置或设备结合地使用的程序。机器可读介质可以是机器可读信号介质或机器可读储存介质。机器可读介质可以包括但不限于电子的、磁性的、光学的、电磁的、红外的、或半导体系统、装置或设备，或者上述内容的任何合适组合。机器可读存储介质的更具体示例会包括基于一个或多个线的电气连接、便携式计算机盘、硬盘、随机存取存储器(ram)、只读存储器(rom)、可擦除可编程只读存储器(eprom或快闪存储器)、光纤、便捷式紧凑盘只读存储器(cd-rom)、光学储存设备、磁储存设备、或上述内容的任何合适组合。
[0146]
为了提供与用户的交互，可以在计算机上实施此处描述的系统和技术，该计算机具有：用于向用户显示信息的显示装置(例如，crt(阴极射线管)或者lcd(液晶显示器)监视器)；以及键盘和指向装置(例如，鼠标或者轨迹球)，用户可以通过该键盘和该指向装置来将输入提供给计算机。其它种类的装置还可以用于提供与用户的交互；例如，提供给用户的反馈可以是任何形式的传感反馈(例如，视觉反馈、听觉反馈、或者触觉反馈)；并且可以用任何形式(包括声输入、语音输入或者、触觉输入)来接收来自用户的输入。
[0147]
可以将此处描述的系统和技术实施在包括后台部件的计算系统(例如，作为数据服务器)、或者包括中间件部件的计算系统(例如，应用服务器)、或者包括前端部件的计算系统(例如，具有图形用户界面或者网络浏览器的用户计算机，用户可以通过该图形用户界面或者该网络浏览器来与此处描述的系统和技术的实施方式交互)、或者包括这种后台部件、中间件部件、或者前端部件的任何组合的计算系统中。可以通过任何形式或者介质的数字数据通信(例如，通信网络)来将系统的部件相互连接。通信网络的示例包括：局域网(lan)、广域网(wan)和互联网。
[0148]
计算机系统可以包括客户端和服务器。客户端和服务器一般远离彼此并且通常通过通信网络进行交互。通过在相应的计算机上运行并且彼此具有客户端-服务器关系的计算机程序来产生客户端和服务器的关系。服务器可以是云服务器，也可以为分布式系统的服务器，或者是结合了区块链的服务器。
[0149]
应该理解，可以使用上面所示的各种形式的流程，重新排序、增加或删除步骤。例如，本发明中记载的各步骤可以并行地执行也可以顺序地执行也可以不同的次序执行，只要能够实现本发明的技术方案所期望的结果，本文在此不进行限制。
[0150]
上述具体实施方式，并不构成对本发明保护范围的限制。本领域技术人员应该明白的是，根据设计要求和其他因素，可以进行各种修改、组合、子组合和替代。任何在本发明的精神和原则之内所作的修改、等同替换和改进等，均应包含在本发明保护范围之内。