论文部分内容阅读
地表三维形变场的提取有助于地质灾害机理的反演。传统的单轨InSAR 形变测量结果只能反映真实的地表三维形变在雷达视线向(LOS)上的一维投影，难以满足地表形变监测的实际需求。理论上，在无任何先验信息情况下，重建地表真实的三维形变场(垂直向、东西向和南北向)，至少需要三个不同方向的观测量。随着新一代SAR 卫星(Sentinel 1、ALOS2 等)的发射升空，利用多颗SAR 卫星平台获取的多轨道数据，可以将InSAR 一维形变测量拓展至三维，最终实现InSAR 三维形变场估计。但是，考虑到InSAR 观测量对南北方向形变不敏感的固有缺陷，如何克服该缺陷以及更高精度的估计地表三维形变成为国内外众多学者的研究热点。本文系统研究了多源异质InSAR(包括MT-InSAR、Offset_tracking 等)观测量解算地表三维形变的函数模型，并深入分析融合多平台、多轨道和多时域InSAR 观测量求解地表三维形变场的精度和局限性。本文研究将GPS 提供的高精度的不同方向形变信息与不同LOS 向InSAR 观测量融合解算，改善InSAR 地表三维形变结果的估计精度。本文以美国南加州地区作为研究区域，由于Offset-tracking 技术的在方位向的形变监测精度较低，本文仅利用MTInSAR技术(如TCPInSAR 技术)对2015-2016 年Sentinel 1 升、降轨数据分别进行处理得到升、降轨两个LOS 向的年平均形变速率。其中，采用欧洲空间局的精确轨道数据对Sentinel 1数据进行轨道精化，地形相位的去除采用1 弧秒的SRTM DEM(空间分辨率约为3030m)，为抑制干涉图的相位噪声，距离向和方位向分别做10 视和2 视的多视处理，并采用改进的Non-Local 滤波法对干涉图进行滤波，提高干涉图的质量。由于InSAR 得到的形变观测量是形变相对值，本文依靠SCIGN 网提供的与InSAR 观测量同步的GPS 台站将InSAR LOS 向形变观测量进行绝对校正，得到真实的LOS 向形变观测量。同时，本文采用普通Kriging 法对GPS 数据进行空间插值，将插值结果与校正后InSAR 观测量进行融合解算，获得整个研究区域的高精度三维形变场。为验证分析融合不同向GPS 形变观测量对于升、降轨InSAR 观测量解算地表三维形变的精度影响，本文共分为3 种方案解算。方案1：仅将GPS 南北方向形变观测量参与三维形变融合估算；方案2：将GPS 水平方向的形变观测量参与三维形变融合估算；方案3：将GPS垂直向和水平向形变观测量参与三维形变融合估算。最终利用未参与融合解算的GPS 台站数据验证其三维形变估算精度。