[发明专利]一种耦合InSAR和数值模拟分析城区滑坡形变的方法

权利要求书

1.一种耦合InSAR和数值模拟分析城区滑坡形变的方法，其特征在于，包括：

获取滑坡体的二维横截面视图，并基于所述二维横截面视图构建一定大小的有限元网格，且初始化所述滑坡体的材料属性及定义用于表征所述滑坡体的真实载荷的边界条件；

基于所述有限元网格、所述材料属性及所述边界条件，建立用于模拟所述滑坡体形变的数值模型，并利用所述数值模型模拟所述滑坡体的形变，得到所述滑坡体的形变模拟值；

获取所述滑坡体在InSAR下的形变测量值，并在所述滑坡体的形变模拟值等于对应的形变测量值的目标条件下调整所述滑坡体的材料属性值，最终得到所述滑坡体的最佳材料属性值；

基于所述最佳材料属性值对应的数值模型，模拟所述滑坡体的真实物理运动过程，得到所述滑坡体的全方位运动。

2.如权利要求1所述的耦合InSAR和数值模拟分析城区滑坡形变的方法，其特征在于，在所述滑坡体的形变模拟值等于对应的形变测量值的目标条件下调整所述滑坡体的材料属性值，最终得到所述滑坡体的最佳材料属性值的过程，包括：

计算目标函数并判断所述目标函数是否达到最小化要求；其中，(p₁,p₂,…,p_m)为所述滑坡体的材料属性集，为所述滑坡体上第i个InSAR点的形变模拟值，U_i为所述滑坡体上第i个InSAR点的形变测量值，n为所述滑坡体上所有InSAR点的数量；

若是，则将当前调整的所述滑坡体的材料属性值作为最佳材料属性值；

若否，则通过遗传算法调整所述滑坡体的材料属性值，并返回执行计算目标函数的操作。

3.如权利要求2所述的耦合InSAR和数值模拟分析城区滑坡形变的方法，其特征在于，所述滑坡体的材料属性集包括有效弹性模量、水力属性、激活孔隙水压力及泊松比。

4.如权利要求1所述的耦合InSAR和数值模拟分析城区滑坡形变的方法，其特征在于，所述边界条件包括降水引起的水力作用、体载荷及建筑施工引起的力的作用。

5.如权利要求1所述的耦合InSAR和数值模拟分析城区滑坡形变的方法，其特征在于，获取所述滑坡体在InSAR下的形变测量值的过程，包括：

基于InSAR，分别获取所述滑坡体上各InSAR点在升轨轨道和降轨轨道上的形变观测值，并将各InSAR点在升轨轨道和降轨轨道上的形变观测值进行二维形变速度分解，得到各InSAR点的东西形变分量和垂直形变分量；其中，所述InSAR点包括所述滑坡体上的PS点和DS点；

分别将各InSAR点的东西形变分量和垂直形变分量进行合并，得到各InSAR点的形变合并量，并各InSAR点的形变合并量对应作为所述滑坡体上各InSAR点的形变测量值。

6.如权利要求5所述的耦合InSAR和数值模拟分析城区滑坡形变的方法，其特征在于，分别获取所述滑坡体上各InSAR点在升轨轨道和降轨轨道上的形变观测值的过程，包括：

获取所述滑坡体的SAR图像，对所述SAR图像进行数据配准，并根据数据配准后的SAR图像生成差分干涉图；其中，所述SAR图像包括升轨图像和降轨图像；

基于所述差分干涉图构建用于监测所述滑坡体的部分PS点形变的第一层网络，并基于所述第一层网络监测所述部分PS点的形变值，得到所述部分PS点的形变观测值；

对所述差分干涉图进行相位优化处理，并基于相位优化处理的差分干涉图构建用于监测所述滑坡体的剩余PS点和所有DS点形变的第二层网络，且基于所述第二层网络监测所述剩余PS点和所有DS点的形变值，得到所述剩余PS点和所有DS点的形变观测值。