基于InSAR与多源数据融合的堆石坝变形场重构方法技术

本发明专利技术公开一种基于InSAR与多源数据融合的堆石坝变形场重构方法包括：获取第一堆石坝垂直变形数据、第二堆石坝垂直向变形数据、第二堆石坝水平向变形数据；基于第一堆石坝垂直变形数据、第二堆石坝垂直变形数据得到第三堆石坝垂直变形数据；获取第一堆石坝方位向变形数据，基于第一堆石坝方位向变形数据、第二堆石坝水平向变形数据通过联合解算模型得到堆石坝的东西向变形矢量、南北向变形矢量；基于第三堆石坝垂直变形数据、东西向变形矢量、南北向变形矢量构建堆石坝外观三维变形场。本发明专利技术将InSAR观测结果与传统监测技术监测结果融合，比单一监测方式质量更高，强调了数据间的互补性，解决了对实际变形漏判或误判的问题。

【技术实现步骤摘要】
基于InSAR与多源数据融合的堆石坝变形场重构方法
本专利技术属于水利工程及岩土工程
，尤其涉及一种基于InSAR与多源数据融合的堆石坝变形场重构方法。
技术介绍
目前堆石坝安全监测技术滞后于坝工设计理论和施工技术的发展，难以满足高堆石坝变形多维度、面广、线长的特点，为了适应高堆石坝的变形监测需求，近年来通过改进现有技术、研发新技术或引进其它领域成熟的技术，在堆石坝变形监测领域出现了一些新的技术手段，比如已被广泛应用于地灾领域的InSAR技术。InSAR技术具有全天候、全天时、覆盖范围广、空间分辨率高、不受恶劣天气影响等优点，其高空间分辨率、连续实时监测的特点极大的弥补了传统监测技术(如水管式沉降仪、引张线式水平位移计、水准测量等)只能在重点部位布设分散、不连续的监测点，按照以点代面、以局部代整体的方法采集离散的变形数据的不足。然而，由于SAR传感器斜视成像的特点，InSAR技术观测的结果为地表真实三维形变在雷达视线向上的投影，难以反映地物的真实变形特征及其演化规律，而且InSAR技术易受到大气和失相关的影响，导致单点监测精度略低于传统监测技术。因此，可以看出新型和传统监测技术各具优势，且有良好的互补性，将两者进行优势互补，并利用多源数据重构堆石坝三维变形场是解决堆石坝变形监测难题的有效方案，但由于新型与传统监测技术的观测机制迥异，多源监测数据的时空分辨率、精度等差异较大，目前尚缺乏合理的三维变形场重构方法。准确评估堆石坝工作状态的需要掌握坝体整体的变形分布，目前常用的方法是进行有限元数值模拟，其主要依赖先验的知识和经验。堆石坝的变形受多因素耦合制约，特别是高堆石坝，由于其建设周期长、填筑蓄水过程复杂、高应力导致的颗粒破碎等，使得堆石坝的应力变形特性复杂且高度非线性，其本构模型参数、实际力学特性等都与先验的认识存在出入，导致有限元数值模拟结果与实际情况存在差异。
技术实现思路
本专利技术要解决的技术问题是基于堆石坝新型与传统监测技术的观测机制迥异，多源监测数据的时空分辨率、精度等差异较大的问题，如何有效地利用变形监测信息重构堆石坝变形场以进行工作性态评估，提供一种基于InSAR与多源数据融合的堆石坝变形场重构方法，利用空-天-地多层次观测体系获得多视角、多时空尺度的堆石坝变形监测数据，在顾及多源数据空间代表性差异的条件下进行最优融合，消除数据之间可能存在的冗余性和互斥性，挖掘它们之间的关联性，使其优势互补以实现堆石坝外观变形场的重构，对全面掌握堆石坝的工作性态具有重要工程应用价值。为实现上述目的，本专利技术采用如下的技术方案：一种基于InSAR与多源数据融合的堆石坝变形场重构方法，包括以下步骤：基于传统方法获取第一堆石坝垂直变形数据，基于时间序列InSAR方法获取堆石坝视线向变形数据，所述堆石坝视线向变形数据包括第二堆石坝垂直向变形数据、第二堆石坝水平向变形数据，所述第二堆石坝垂直变形数据为堆石坝视线向变形的垂直分解矢量，所述第二堆石坝水平向变形数据为堆石坝视线向变形水平分解矢量；基于所述第一堆石坝垂直变形数据、所述第二堆石坝垂直变形数据通过多源数据融合方法得到第三堆石坝垂直变形数据；基于改进的多孔径InSAR方法获取第一堆石坝方位向变形数据，基于所述第一堆石坝方位向变形数据、所述第二堆石坝水平向变形数据通过联合解算模型得到堆石坝的东西向变形矢量、南北向变形矢量；基于所述第三堆石坝垂直变形数据、所述东西向变形矢量、所述南北向变形矢量构建堆石坝外观三维变形场。优选地，基于数据同化算法将所述第二堆石坝垂直变形数据与所述第一堆石坝垂直变形数据进行融合，得到所述第三堆石坝垂直变形数据。优选地，获取第一堆石坝方位向变形数据的步骤为通过改进的多孔径InSAR方法得到时间序列的方位向变形，即得到第一堆石坝方位向变形数据。优选地，获取所述堆石坝视线向变形数据的方法为：基于SAR影像的成像几何，对于地面任意一点的视线向变形矢量分解为垂直向投影的变形矢量和水平向投影的变形矢量，所述垂直向投影的变形矢量即为第二堆石坝垂直向变形数据，所述水平向投影的变形矢量即为第二堆石坝水平向变形数据。优选地，基于所述南北向变形矢量和所述东西向变形矢量，建立有雷达视线向变形和地表三维变形的关系；基于改进的多孔径InSAR技术获取的方位向变形DAZI，位于由东西向和南北向组成的二维平面中，根据集合关系建立方位向变形和地表三维变形的关系。优选地，所述南北向变形矢量和所述东西向变形矢量为：DLOS＝DU·cosθ+DH·sinθDLOS是视线向变形矢量，DU、DH分别为垂直向变形矢量和水平向变形矢量；DE、DN为东西向和南北向的变形矢量；θ为雷达入射角，入射角是入射光线与表面法向的夹角，α为卫星飞行方位角即正北方向和卫星飞行方向的顺时针夹角，为正北方向与水平向的夹角。优选地，所述雷达视线向变形和地表三维变形的关系：优选地，所述方位向变形和地表三维变形的关系为：其中，DAZI为方位向变形矢量。优选地，所述改进的多孔径InSAR的方法具体为：结合多孔径InSAR方法和SBAS-InSAR方法，根据与SBAS方法相同的规则改进所述多孔径InSAR方法的数据处理流程，将SAR影像数据集划分为以较小的时间和空间基线为特征的干涉组合，利用多孔径InSAR方法生成MAI干涉图，并应用SBAS反演步骤得到方位向时间序列的变形信息。本专利技术的技术效果为：本专利技术将InSAR观测结果与传统监测技术监测结果实现智能融合，获取比单一监测方式更高质量的结果(例如精度、分辨率等指标)，强调了数据间的互补性，并通过联合解算模型将多维度监测数据重构堆石坝外观三维变形场以全面地掌握堆石坝的变形性态，解决了单维度监测可能造成对实际变形漏判或误判的问题。附图说明图1为本专利技术的流程示意图；图2为本专利技术改进的多孔径InSAR技术数据处理流程图；图3为本专利技术基于数据同化算法的数据融合过程线对比图；图4为本专利技术数据融合前后垂直向变形对比图，(a)数据融合前，(b)数据融合后。图5为数据融合前后精度对比图；图6为堆石坝外观三维变形场重构图。具体实施方式下面通过具体实施方式结合附图对本专利技术作进一步详细说明。参照图1-5，重构的水布垭面板堆石坝外观变形场直观的反映了坝体在2007年2月至2009年1月期间的整体变形状态，其中颜色和箭头分别表示垂直和水平年平均变化速率(图6)。此期间正值水布垭面板堆石坝初期集中蓄水阶段，坝体出现以沉降为主的明显变形，同一点处年沉降速率普遍大于年平均水平变形速率，其中沉降最大值出现在下游面板的顶部和中部，并呈现沿顺河向向下游和沿坝轴向向两岸逐渐减小的规律；水平位移沿顺河向与沉降有相似的分布规律，最大值同样出现在下游面板的顶部和中部，沿顺河向和向两岸逐渐减小，左右岸发生指向对岸的相互变形，且左岸大于右岸，这是由于右岸岸坡的较缓的原因本文档来自技高网...

【技术保护点】
1.一种基于InSAR与多源数据融合的堆石坝变形场重构方法，其特征在于，包括以下步骤：/n基于传统监测方法获取第一堆石坝垂直变形数据，基于时间序列InSAR方法获取堆石坝视线向变形数据，所述堆石坝视线向变形数据包括第二堆石坝垂直向变形数据、第二堆石坝水平向变形数据，/n所述第二堆石坝垂直变形数据为堆石坝视线向变形的垂直分解矢量，所述第二堆石坝水平向变形数据为堆石坝视线向变形的水平分解矢量；/n基于所述第一堆石坝垂直变形数据、所述第二堆石坝垂直变形数据通过多源数据融合方法得到第三堆石坝垂直变形数据；/n基于改进的多孔径InSAR方法获取第一堆石坝方位向变形数据，/n基于所述第一堆石坝方位向变形数据、所述第二堆石坝水平向变形数据通过联合解算模型得到堆石坝的东西向变形矢量、南北向变形矢量；/n基于所述第三堆石坝垂直变形数据、所述东西向变形矢量、所述南北向变形矢量构建堆石坝外观三维变形场。/n

【技术特征摘要】
1.一种基于InSAR与多源数据融合的堆石坝变形场重构方法，其特征在于，包括以下步骤：
基于传统监测方法获取第一堆石坝垂直变形数据，基于时间序列InSAR方法获取堆石坝视线向变形数据，所述堆石坝视线向变形数据包括第二堆石坝垂直向变形数据、第二堆石坝水平向变形数据，
所述第二堆石坝垂直变形数据为堆石坝视线向变形的垂直分解矢量，所述第二堆石坝水平向变形数据为堆石坝视线向变形的水平分解矢量；
基于所述第一堆石坝垂直变形数据、所述第二堆石坝垂直变形数据通过多源数据融合方法得到第三堆石坝垂直变形数据；
基于改进的多孔径InSAR方法获取第一堆石坝方位向变形数据，
基于所述第一堆石坝方位向变形数据、所述第二堆石坝水平向变形数据通过联合解算模型得到堆石坝的东西向变形矢量、南北向变形矢量；
基于所述第三堆石坝垂直变形数据、所述东西向变形矢量、所述南北向变形矢量构建堆石坝外观三维变形场。


2.如权利要求1所述的基于InSAR与多源数据融合的堆石坝变形场重构方法，其特征在于，
基于数据同化算法将所述第二堆石坝垂直变形数据与所述第一堆石坝垂直变形数据进行融合，得到所述第三堆石坝垂直变形数据。


3.如权利要求1所述的基于InSAR与多源数据融合的堆石坝变形场重构方法，其特征在于，
获取第一堆石坝方位向变形数据的步骤为通过改进的多孔径InSAR方法得到时间序列的方位向变形，即得到第一堆石坝方位向变形数据。


4.如权利要求3所述的基于InSAR与多源数据融合的堆石坝变形场重构方法，其特征在于，获取所述堆石坝视线向变形数据的方法为：基于SAR影像的成像几何，对于地面任意一点的视线向变形矢量分解为垂直向投影的变形矢量和水平向投影的变形矢量，所述垂直向投影的变形矢量即为第二堆石坝垂直向变形数据，所述水平向投影的变形矢量即为第二堆石坝水平向变形数据。...

【专利技术属性】
技术研发人员：马刚，郭承乾，周伟，程翔，常晓林，
申请(专利权)人：武汉大学，
类型：发明
国别省市：湖北;42