CN117854257B - 基于地基sar监测变形数据的次生灾害预警方法 - Google Patents

Abstract

本发明公开了一种基于地基SAR监测变形数据的次生灾害预警方法，属于地质灾害预警技术领域，包括如下步骤：通过地基SAR对次生灾害不稳定区域进行实时监测，以获取变形监测数据；构建次生灾害不稳定区域随时间产生的位移与时间曲线；设置监测周期，并基于位移与时间曲线，构建位移变化三角面积模型；基于位移变化三角面积模型构建面积与时间曲线，按监测周期分别获取面积与时间曲线的变化趋势，以及位移与时间曲线中的位移三角面积；将次生灾害发展过程划分为初始变形阶段、等速变形阶段和加速变形阶段；分别对应各阶段设置次生灾害预警等级，以分阶进行次生灾害预警；解决了无法快速全面判断识识别次生灾害发展阶段及其变形演化趋势的问题。

Description

基于地基SAR监测变形数据的次生灾害预警方法

技术领域

本发明属于地质灾害预警技术领域，尤其涉及一种基于地基SAR监测变形数据的次生灾害预警方法。

背景技术

山区地形陡峻、河谷深切、构造活跃，斜坡在长期地质营力和外动力作用下，往往容易发生突发性的、高位的且隐蔽性极强的地质灾害。在滑坡以及崩塌等高位崩滑地质灾害发生后，往往存在进一步滑塌的可能，这给首次发生滑坡以及崩塌后的救援带来了困扰，也时刻威胁着抢险现场人员和设备的安全。

当前常见的次生灾害预警判识依据主要是基于三阶段变形特征而建立的位移、速度倒数、加速度、切线角和改进切线角等预警模型，但由于滑坡的二次灾害变形具有隐蔽性和突发性，并且上述预警模型是依靠大量的滑坡实例来总结得到的，导致在使用这些预警模型进行新的次生地质灾害预警时，往往出现误报、漏报的现象，从而导致不能准确地实现次生地质灾害预警。虽然地基SAR（synthetic aperture radar，合成孔径雷达）监测数据精度较高，但在面对已发生地质灾害的危险区域时，由于监测设备无法及时架设，并且次生灾害会严重影响滑坡形变数据的获取，致使监测设备无法监测到不稳定区域发生变形破坏的完整过程，以至于获取的到的位移曲线可能没有完整的三个变形阶段，无法对滑坡的变形演化趋势进行判断，且进一步会导致预警信息不能及时发布，严重威胁到现场抢险救灾人员的生命安全。现有的地质灾害预警方法难以在应急场景下突破灾害体前兆信息，实现抢险现场滑坡以及崩塌的全局监测和实时预警。

发明内容

针对现有技术中的上述不足，本发明提供的一种基于地基SAR监测变形数据的次生灾害预警方法，以地基SAR现场监测数据为基础，通过深化处理地基SAR所获取的数据，以位移与时间的曲线面积变化趋势为基础进行分析预警，实现了更为准确和直接有效地预警次生地质灾害，解决了无法快速全面判断识识别次生灾害发展阶段及其变形演化趋势的问题。

为了达到上述发明目的，本发明采用的技术方案为：

本发明提供的一种基于地基SAR监测变形数据的次生灾害预警方法，包括如下步骤：

S1、通过地基SAR对次生灾害不稳定区域进行实时监测，以获取变形监测数据；

S2、基于变形监测数据，得到次生灾害不稳定区域随时间产生的位移曲线，并作为位移与时间曲线；

S3、设置监测周期，并基于位移与时间曲线，构建位移变化三角面积模型；

S4、基于位移变化三角面积模型构建面积与时间曲线，按监测周期分别获取面积与时间曲线的变化趋势，以及位移与时间曲线中的位移三角面积；

S5、基于面积与时间曲线的变化趋势和位移三角面积，将次生灾害发展过程划分为初始变形阶段、等速变形阶段和加速变形阶段；

S6、根据面积与时间曲线，结合划分的初始变形阶段、等速变形阶段和加速变形阶段，分别对应设置次生灾害预警等级；

S7、基于次生灾害预警等级，分阶进行次生灾害预警。

本发明的有益效果为：本发明提供的一种基于地基SAR监测变形数据的次生灾害预警方法，通过地基SAR实时监测到次生灾害不稳定区域的变形监测数据，并基于变形监测数据构建了次生灾害不稳定区域随时间产生的位移与时间曲线，通过对位移与时间曲线设置监测周期，构建了位移变化三角面积模型，并划分得到了各监测周期内对应的位移三角面积，本发明基于位移三角面积和面积与时间曲线的斜率将次生灾害发展过程进行了三级阶段划分，并对应各阶段进行了预警等级设置，尤其对于加速变形阶段的警报等级进一步进行了三阶划分，本发明能够全面快速地对次生灾害的发展阶段进行准确的识别，充分反映不稳定区域的变形演化趋势，同时能够为现场救灾人员提供准确的预警信号，以保障救灾人员的生命安全。

进一步地，所述S1包括如下步骤：

S11、圈定次生灾害不稳定区域；

S12、以次生灾害不稳定区域为监测目标架设地基SAR；

S13、利用地基SAR实时监测次生灾害不稳定区域，得到变形监测数据。

采用上述进一步方案的有益效果为：本发明对次生灾害不稳定区域进行圈定，并以次生灾害不稳定区域为监测目标架设地基SAR进行实时监测，能够获取到检测目标处的变形监测数据，为对变形监测数据进行深度分析以全面快速识别次生灾害发展阶段提供数据基础。

进一步地，所述S3包括如下步骤：

S31、设置地基SAR对次生灾害不稳定区域的监测周期；

S32、基于位移与时间曲线，获取次生灾害不稳定区域在各监测周期内的累计位移量；

S33、基于次生灾害不稳定区域在各监测周期内的累计位移量，构建位移增量模型；

S34、基于位移增量模型，构建位移变化三角面积模型。

采用上述进一步方案的有益效果为：本发明设置了对次生灾害不稳定区域的监测周期，以便于对位移与时间曲线进行分监测周期的分析，为从次生灾害不稳定区域的变形速率和变形加速度划分次生灾害发展阶段提供基础。

进一步地，所述位移增量模型的计算表达式如下：

，

其中，表示次生灾害不稳定区域在第n个监测周期内的位移增量，S_n表示次生灾害不稳定区域在第n个监测周期内的累计位移量，S_n-1表示次生灾害不稳定区域在第n-1个监测周期内的累计位移量。

采用上述进一步方案的有益效果为：本发明提供位移增量模型的计算方法，通过位移增量模型能够计算得到每个一个监测周期内次生灾害体的位移变形量，能够与曲线斜率一样反映次生灾害不稳定区域的变形速率，也能够为计算位移三角面积提供基础。

进一步地，所述位移变化三角面积模型的计算表达式如下：

，

其中，A_n表示第n个监测周期内的位移变化三角面积，t_n表示第n个监测周期结束的时刻，t_n-1表示第n-1个监测周期结束的时刻。

采用上述进一步方案的有益效果为：本发明提供位移变化三角面积模型的计算方法，由于位移的变化速率不同，通过位移三角面积变化能够体现次生灾害不稳定区域的变形加速度，依据位移三角面积变化规律能够为准确划分次生灾害发展阶段提供基础。

进一步地，所述S5包括如下步骤：

S51、获取各监测周期内位移与时间曲线的位移三角面积，以及位移与时间曲线的斜率，并将面积与时间曲线的斜率作为次生灾害变形速率，其中，位移与时间曲线的斜率，对应面积与时间曲线的变化趋势；

S52、若位移三角面积大于第一位移三角面积阈值，同时次生灾害变形速率变化趋势为由小变大再变小，则将位移与时间曲线以及面积与时间曲线中对应的部分作为初始变形阶段；

S53、若位移三角面积小于第一位移三角面积阈值，则次生灾害变形速率变化趋势为基本不变，并将位移与时间曲线以及面积与时间曲线中对应的部分作为等速变形阶段；

S54、若位移三角面积大于第一位移三角面积阈值，同时次生灾害变形速率的变化趋势为逐渐增大，直至渐进于竖直状态，则将位移与时间曲线以及面积与时间曲线中对应部分作为加速变形阶段。

采用上述进一步方案的有益效果为：本发明提供基于面积与时间曲线的变化趋势和位移三角面积，将将次生灾害发展过程划分为初始变形阶段、等速变形阶段和加速变形阶段的具体方法，实现了对次生灾害发展阶段的准确三阶段划分，为进行不同紧要程度的预警提供了基础。

进一步地，所述S6包括如下步骤：

S61、根据面积与时间曲线，将初始变形阶段对应的次生灾害预警等级设置为注意级T1；

S62、根据面积与时间曲线，将等速变形阶段对应的次生灾害预警等级设置为警示级T2；

S63、根据面积与时间曲线，获取加速变形阶段中各监测周期内的位移三角面积；

S64、若加速变形阶段中位移三角面积小于第二位移三角面积阈值，则将面积与时间曲线中对应的部分作为初加速阶段，并将初加速阶段对应的次生灾害预警等级设置为第一警报级T3；

S65、若加速变形阶段中位移三角面积大于第二位移三角面积阈值且小于第三位移三角面积阈值，则将面积与时间曲线中对应的部分作为中加速阶段，并将中加速阶段对应的次生灾害预警等级设置为第二警报级T4；

S66、若加速变形阶段中位移三角面积大于第三位移三角面积阈值，则将面积与时间曲线中对应的部分作为临滑阶段，并将临滑阶段对应的次生灾害预警等级设置为第三警报级T5。

采用上述进一步方案的有益效果为：本发明提供了根据面积与时间曲线，结合划分的初始变形阶段、等速变形阶段和加速变形阶段，分别对应设置次生灾害预警等级的方法，通过分别分阶段设置对应等级的预警等级，尤其对加速变形阶段进行进一步的三阶划分，能够实现根据灾害发生严重程度，实时准确预警。

与现有技术相比，本发明具有以下优势：

本发明在监测过程中无需获取次生灾害不稳定区域发生变形过程对应的完整变形监测数据，只需根据能够获取到的或已经获取到的变形监测数据进行构建位移与时间曲线后，通过位移变化三角面积模型进行变形演化趋势判别、次生灾害发展过程划分和次生灾害预警等级设置，即可实现依据现场抢险救灾的实际具体情况，快速进入监测预警状态，实现实时的快速预警，由此，本发明能够实现阶段性快速分析识别当前次生灾害发展阶段，从而克服地基SAR在次生灾害预警监测过程中因首次灾害发生后救援时间仓促、次生灾害监测时因仪器部署时间耗时长、次生灾害突发性和隐蔽性强等特点导致难以快速获取次生灾害变形演化趋势，以供现场第一决策做出准确预判的问题。

针对于本发明还具有的其他优势将在后续的实施例中进行更细致的分析。

附图说明

为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。

图1为本发明实施例1中一种基于地基SAR监测变形数据的次生灾害预警方法的步骤流程图。

图2为本发明实施例1中位移与时间曲线的示意图。

图3为本发明实施例1中面积与时间曲线的示意图。

图4为本发明实施例2中位移与时间曲线的示意图。

图5为本发明实施例2中面积与时间曲线的示意图。

图6为本发明实施例3中位移与时间曲线的示意图。

图7为本发明实施例3中面积与时间曲线的示意图。

具体实施方式

下面将结合本发明实施例中附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。因此，以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施例。基于本发明的实施例，本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。

实施例1：

如图1所示，在本发明的一个实施例中，本发明提供一种基于地基SAR监测变形数据的次生灾害预警方法，包括如下步骤：

S1、通过地基SAR对次生灾害不稳定区域进行实时监测，以获取变形监测数据；

所述S1包括如下步骤：

S11、圈定次生灾害不稳定区域；

S12、以次生灾害不稳定区域为监测目标架设地基SAR；

S13、利用地基SAR实时监测次生灾害不稳定区域，得到变形监测数据。地基SAR在获取变形监测数据时，原理为地基SAR设备基于微波干涉技术的远程监控系统，将步进频率连续波技术、合成孔径雷达技术以及差分干涉测量技术相结合，将次生灾害不稳定区域划分为多个微小的像素单元，并通过每个微小的像素单元的连续两次电磁波的相位差获取监测目标的形变信息，以构成变形监测数据。

S2、基于变形监测数据，得到次生灾害不稳定区域随时间产生的位移曲线，并作为位移与时间曲线；本实施例中，基于变形监测数据，将地基SAR获取的监测目标在不同时刻的累计位移S与时刻T的数据，绘制呈关系曲线，即位移与时间曲线；位移与时间曲线的构建方法包括如下步骤：

A1、基于变形监测数据，将监测目标的累计位移S与对应的时间T导入表格，通过表格将时间T与累计位移S相关联；

A2、以时间T为横轴，累计位移S为纵轴，绘制具有平滑曲线和数据标记的散点图，并作为初始位移与时间曲线图；

A3、对初始位移与时间曲线图进行去离散处理，优化得到位移与时间曲线；

位移与时间曲线体现次生灾害不稳定区域在一个时间段内的位移变化，由于架设设备时间以及次生灾害发生的突发性，实际获取到的位移与时间曲线可能无法完整体现出次生灾害不稳定区域完整的变形阶段。

S3、设置监测周期，并基于位移与时间曲线，构建位移变化三角面积模型；

所述S3包括如下步骤：

S31、设置地基SAR对次生灾害不稳定区域的监测周期；本实施例中，次生灾害不稳定区域的监测周期设置为30min，即1800s；

S32、基于位移与时间曲线，获取次生灾害不稳定区域在各监测周期内的累计位移量；

S33、基于次生灾害不稳定区域在各监测周期内的累计位移量，构建位移增量模型；

所述位移增量模型的计算表达式如下：

，

其中，表示次生灾害不稳定区域在第n个监测周期内的位移增量，S_n表示次生灾害不稳定区域在第n个监测周期内的累计位移量，S_n-1表示次生灾害不稳定区域在第n-1个监测周期内的累计位移量。

S34、基于位移增量模型，构建位移变化三角面积模型。

所述位移变化三角面积模型的计算表达式如下：

，

其中，A_n表示第n个监测周期内的位移变化三角面积，t_n表示第n个监测周期结束的时刻，t_n-1表示第n-1个监测周期结束的时刻。

S4、基于位移变化三角面积模型构建面积与时间曲线，按监测周期分别获取面积与时间曲线的变化趋势，以及位移与时间曲线中的位移三角面积；

如图2所示，在每个监测周期内，以每个监测周期的位移变化量与时间所形成的三角形面积作为该监测周期对应的位移三角面积，在连续相邻的监测周期内，由于次生灾害不稳定区域位移的变化速率不同，其所形成的位移三角面积存在明显区别，利用位移三角面积随时间变化的规律能够划分且准确识别出次生灾害不稳定区域的变形速率以及变形演化趋势，图中纵坐标为位移，单位为mm，横坐标为监测时间，单位为s。

S5、基于面积与时间曲线的变化趋势和位移三角面积，将次生灾害发展过程划分为初始变形阶段、等速变形阶段和加速变形阶段；

所述S5包括如下步骤：

S51、获取各监测周期内位移与时间曲线的位移三角面积，以及位移与时间曲线的斜率，并将面积与时间曲线的斜率作为次生灾害变形速率，其中，位移与时间曲线的斜率，对应面积与时间曲线的变化趋势；

S52、若位移三角面积大于第一位移三角面积阈值，同时次生灾害变形速率变化趋势为由小变大再变小，则将位移与时间曲线以及面积与时间曲线中对应的部分作为初始变形阶段；在本实施例中，第一位移三角面积阈值为0.1mm·s，即表示该监测周期内次生灾害不稳定区域的时空间隔达到了0.1mm·s；

S53、若位移三角面积小于第一位移三角面积阈值，则次生灾害变形速率变化趋势为基本不变，并将位移与时间曲线以及面积与时间曲线中对应的部分作为等速变形阶段；

S54、若位移三角面积大于第一位移三角面积阈值，同时次生灾害变形速率的变化趋势为逐渐增大，直至渐进于竖直状态，则将位移与时间曲线以及面积与时间曲线中对应部分作为加速变形阶段。

如图2所示，在位移与时间曲线图中，初始变形阶段对应AB段部分。初始变形阶段在坡体初期变形时逐渐产生，坡体中出现明显的裂缝，在受到外力和自重作用下，裂缝逐渐压密，位移与时间曲线最初表现出较大的斜率，随时间持续，位移与时间曲线的AB段部分中存在一个拐点，即在该拐点处的斜率最大，这对应于次生灾害不稳定区域中的次生灾害体由于受到外力及自重作用，具有相对大的变形速率，此后裂缝压密完成，变形速率又逐渐减小，且位移与时间曲线的斜率逐渐减小，表现出减速变形的特征，即在初始变形阶段次生灾害变形速率的变化趋势呈现为由小变大再变小，位移与时间曲线呈现为上凸形态，直至到达B点时，达到初始变形阶段最小斜率。

在位移与时间曲线图中，等速变形阶段对应BC段部分。在等速变形阶段中，次生灾害不稳定区域的坡体，在重力作用下以等速发展的趋势继续变形，但由于外力及自身重力作用产生下滑力，同时次生灾害不稳定区域底部又受到滑床部位的摩擦阻力，且外力及自身重力所产生的下滑力与摩擦阻力大小相当，导致次生灾害体的变形呈现为等速变形，在位移与时间曲线中呈现为一条倾斜的直线且直线斜率总体不变，次生灾害变形速率也基本保持不变，变形的加速度几乎为0，直至到达C点。

在位移与曲线图中，加速变形阶段对应CF段部分：随次生灾害不稳定区域变形的持续，当次生灾害体底部受到滑床部位的摩擦力小于外力及自身重力导致的下滑力时，次生灾害体则逐渐加速变形，此过程中，次生灾害变形速率逐渐增大，呈现出不断增长的趋势，位移与时间曲线呈现为下凹的形态，尤其在临近次生灾害不稳定区域破坏的时候，曲线渐进于直立状态，次生灾害变形速率极大，变形的加速度也极大。本实施例中，在加速变形阶段根据加速度量级（即位移三角面积）的不同，可将该阶段进行进一步的划分，划分结果分别为初加速阶段、中加速阶段和临滑阶段，其中，初加速阶段对应CD段部分，中加速阶段对应DE段部分，临滑阶段对应EF段部分。

S6、根据面积与时间曲线，结合划分的初始变形阶段、等速变形阶段和加速变形阶段，分别对应设置次生灾害预警等级；

如图3所示，图中纵坐标为位移三角面积，单位为mm·s，横坐标为时间，单位为s，0~T₁对应初始变形阶段，次生灾害变形速率的变化趋势为由小变大再变小时，在面积与时间曲线中体现为曲线呈下降趋势，上凹型；T₁~T₂对应等速变形阶段，次生灾害变形速率的变化趋势为基本不变时，在面积与时间曲线中体现为曲线平稳，基本保持不变的直线型；T₂~T₅对应加速变形阶段，次生灾害变形速率的变化趋势为逐渐增大，直至渐进于竖直状态时，在面积与时间曲线中体现为曲线呈上升趋势，下凹型；在临滑阶段中，根据位移三角面积的不同，将加速变形阶段进一步划分为了，T₂~T₃对应的初加速阶段、T₃~T₄对应的中加速阶段和T₄~T₅对应的临滑加速阶段。

所述S6包括如下步骤：

S61、根据面积与时间曲线，将初始变形阶段对应的次生灾害预警等级设置为注意级T1；

S62、根据面积与时间曲线，将等速变形阶段对应的次生灾害预警等级设置为警示级T2；

S63、根据面积与时间曲线，获取加速变形阶段中各监测周期内的位移三角面积；

S64、若加速变形阶段中位移三角面积小于第二位移三角面积阈值，则将面积与时间曲线中对应的部分作为初加速阶段，并将初加速阶段对应的次生灾害预警等级设置为第一警报级T3；本实施例中，第二位移三角面积阈值为0.5mm·s，即表示该监测周期内次生灾害不稳定区域的时空间隔达到了0.5mm·s；

S65、若加速变形阶段中位移三角面积大于第二位移三角面积阈值且小于第三位移三角面积阈值，则将面积与时间曲线中对应的部分作为中加速阶段，并将中加速阶段对应的次生灾害预警等级设置为第二警报级T4；本实施例中，第三位移三角面积阈值为1.3mm·s，即表示该监测周期内次生灾害不稳定区域的时空间隔达到了1.3mm·s；

S66、若加速变形阶段中位移三角面积大于第三位移三角面积阈值，则将面积与时间曲线中对应的部分作为临滑阶段，并将临滑阶段对应的次生灾害预警等级设置为第三警报级T5。

S7、基于次生灾害预警等级，分阶进行次生灾害预警。

本实施例中，对于预警等级设置为注意级T1时，应在处于次生灾害威胁范围内的村庄和公路边树立“滑坡危险区域，观察通过”的标牌，以进行告警标识；对于预警等级设置为警示级T2时，应加强群防群监测，时刻监测次生灾害不稳定区域的降雨量、墙体裂缝、地裂缝的情况；对于预警等级设置为第一警报级T3时，应通过救灾抢险中心发出蓝色警报，警戒受次生灾害威胁对象，做好撤离逃生准备；对于预警等级设置为第二警报级T4时，应通过救灾抢险中心发出橙色警报，警戒受次生灾害威胁对象，时刻准备撤离；对于预警等级设置为第三警报级T5时，应通过救灾抢险中心发出红色警报，警戒受次生灾害威胁对象，立即撤离。

在本实施例中，若获取的变形监测数据不足以演化完整的次生灾害不稳定区域变形情况时，采用实际获取到的变形监测数据依次构建位移与时间曲线、三角面积模型，以及面积与时间曲线，仍能够通过S5和S6对实际情况下的次生灾害不稳定区域根据时间阶段进行准确的分级预警，从而保障救灾人员的生命安全。

实施例2：

如图4所示，在本发明的一个实用实例中，在滑坡后，以滑坡灾害抢险救援现场中的地基SAR实时监测的变形监测数据为数据支撑，构建了该滑坡情况下的位移与时间曲线图，图中纵坐标为位移，单位为mm，横坐标为监测时间，单位为s。

如图5所示，以30分钟为一个监测周期，基于位移增量模型获取滑坡不稳定区域的累计位移量；并通过位移变化三角模型，得到每个监测周期内位移量与时间构成的三级面积，以得到面积与时间曲线，图中纵坐标为位移三角面积，单位为mm·s，横坐标为时间，单位为s。由图5能够发现，从滑坡不稳定区域的变形监测数据开始获取时，滑坡已进入到等速变形阶段，在监测周期达到4.8×10⁵s（8000分钟）时，滑坡进入初级加速阶段，此时预警级别到达警报级，具体为第一警报级T3，在第一警报级T3下，救援人员应该迅速撤离灾害发生区域，防止次生灾害的发生影响救援人员的生命安全。

实施例3：

如图6所示，在某高位牵引式滑坡灾后，以地基SAR实时监测的变形监测数据为支撑，构建了该滑坡情况下的位移与时间曲线图，图中纵坐标为位移，单位为mm，横坐标为监测时间，单位为s。

如图7所示，与实施例1和实施例2中对变形监测数据和位移与时间曲线的处理方式相同，设置30分钟为一个监测周期，构建位移与面积曲线，图中纵坐标为位移三角面积，单位为mm·s，横坐标为时间，单位为s。由图7能够发现，滑坡自地基SAR开始监测后已进入加速变形阶段，且在监测周期达到1.18×10⁷s（19710分钟）时，滑坡进入临滑阶段，此时次生灾害极易发生，应当立即执行第三警报级T5的警报。

以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。

Claims (6)

1.一种基于地基SAR监测变形数据的次生灾害预警方法，其特征在于，包括如下步骤：

S1、通过地基SAR对次生灾害不稳定区域进行实时监测，以获取变形监测数据；

S2、基于变形监测数据，得到次生灾害不稳定区域随时间产生的位移曲线，并作为位移与时间曲线；

S3、设置监测周期，并基于位移与时间曲线，构建位移变化三角面积模型；

S4、基于位移变化三角面积模型构建面积与时间曲线，按监测周期分别获取面积与时间曲线的变化趋势，以及位移与时间曲线中的位移三角面积；

S5、基于面积与时间曲线的变化趋势和位移三角面积，将次生灾害发展过程划分为初始变形阶段、等速变形阶段和加速变形阶段；

所述S5包括如下步骤：

S51、获取各监测周期内位移与时间曲线的位移三角面积，以及位移与时间曲线的斜率，并将面积与时间曲线的斜率作为次生灾害变形速率，其中，位移与时间曲线的斜率，对应面积与时间曲线的变化趋势；

S52、若位移三角面积大于第一位移三角面积阈值，同时次生灾害变形速率变化趋势为由小变大再变小，则将位移与时间曲线以及面积与时间曲线中对应的部分作为初始变形阶段；

S53、若位移三角面积小于第一位移三角面积阈值，则次生灾害变形速率变化趋势为基本不变，并将位移与时间曲线以及面积与时间曲线中对应的部分作为等速变形阶段；

S54、若位移三角面积大于第一位移三角面积阈值，同时次生灾害变形速率的变化趋势为逐渐增大，直至渐进于竖直状态，则将位移与时间曲线以及面积与时间曲线中对应部分作为加速变形阶段；

S6、根据面积与时间曲线，结合划分的初始变形阶段、等速变形阶段和加速变形阶段，分别对应设置次生灾害预警等级；

S7、基于次生灾害预警等级，分阶进行次生灾害预警。

2.根据权利要求1所述的基于地基SAR监测变形数据的次生灾害预警方法，其特征在于，所述S1包括如下步骤：

S11、圈定次生灾害不稳定区域；

S12、以次生灾害不稳定区域为监测目标架设地基SAR；

S13、利用地基SAR实时监测次生灾害不稳定区域，得到变形监测数据。

3.根据权利要求1所述的基于地基SAR监测变形数据的次生灾害预警方法，其特征在于，所述S3包括如下步骤：

S31、设置地基SAR对次生灾害不稳定区域的监测周期；

S32、基于位移与时间曲线，获取次生灾害不稳定区域在各监测周期内的累计位移量；

S33、基于次生灾害不稳定区域在各监测周期内的累计位移量，构建位移增量模型；

S34、基于位移增量模型，构建位移变化三角面积模型。

4.根据权利要求3所述的基于地基SAR监测变形数据的次生灾害预警方法，其特征在于，所述位移增量模型的计算表达式如下：

，

其中，表示次生灾害不稳定区域在第n个监测周期内的位移增量，S _n 表示次生灾害不稳定区域在第n个监测周期内的累计位移量，S _n-1表示次生灾害不稳定区域在第n-1个监测周期内的累计位移量。

5.根据权利要求4所述的基于地基SAR监测变形数据的次生灾害预警方法，其特征在于，所述位移变化三角面积模型的计算表达式如下：

，

其中，A _n 表示第n个监测周期内的位移变化三角面积，t _n 表示第n个监测周期结束的时刻，t _n-1表示第n-1个监测周期结束的时刻。

6.根据权利要求1所述的基于地基SAR监测变形数据的次生灾害预警方法，其特征在于，所述S6包括如下步骤：

S61、根据面积与时间曲线，将初始变形阶段对应的次生灾害预警等级设置为注意级T1；

S62、根据面积与时间曲线，将等速变形阶段对应的次生灾害预警等级设置为警示级T2；

S63、根据面积与时间曲线，获取加速变形阶段中各监测周期内的位移三角面积；

S64、若加速变形阶段中位移三角面积小于第二位移三角面积阈值，则将面积与时间曲线中对应的部分作为初加速阶段，并将初加速阶段对应的次生灾害预警等级设置为第一警报级T3；

S65、若加速变形阶段中位移三角面积大于第二位移三角面积阈值且小于第三位移三角面积阈值，则将面积与时间曲线中对应的部分作为中加速阶段，并将中加速阶段对应的次生灾害预警等级设置为第二警报级T4；

S66、若加速变形阶段中位移三角面积大于第三位移三角面积阈值，则将面积与时间曲线中对应的部分作为临滑阶段，并将临滑阶段对应的次生灾害预警等级设置为第三警报级T5。