基于多源InSAR数据的三维地表形变解算方法研究

王志伟

1. 山东建筑大学测绘地理信息学院，山东 济南 250101; 2. 山东科技大学测绘科学与工程学院，山东 青岛 266590

随着不同成像参数SAR卫星不断发射升空以及新的InSAR技术不断出现，同一地区可以获取多尺度、多源、异质的InSAR地表形变监测数据。将这些多源InSAR地表形变监测数据作为观测量，实际三维地表形变作为未知参数，建立模型并解算，实现东西、南北和垂直向地表形变监测，为地质灾害预报及监测提供更全面的基础资料。本文以InSAR三维地表形变解算为核心研究内容，针对大梯度、轨道和地形、粗差等因素对解算精度的影响，利用现代测量数据处理理论和方法，消除或抑制上述因素导致的各类噪声，实现多源InSAR数据的三维地表形变高精度解算。主要包括以下几个方面的研究内容：

(1) 系统研究了SAR成像及InSAR干涉相位技术原理，特别对视线向及方位向地表形变监测技术进行了深入细致的分析与研究，结合技术原理给出主要的数据处理流程。在此基础上，获取了巴姆地震方位向地表形变信息。相对于常规InSAR技术获取的一维视线向地表形变而言，方位向地表形变可以提供更丰富的地表形变信息。

(2) 针对大梯度、失相干等对干涉相位技术的影响，提出一种多种InSAR技术联合的大梯度三维地表形变解算方法，该方法基于升、降轨SAR影像，通过联合D-InSAR和Offset-tracking获取视线向地表形变，联合MAI和Offset-tracking获取方位向地表形变，在此基础上，针对两个方向地表形变监测精度不同，采用Helmert方差分量估计的方式确定权重，建立三维地表形变模型并解算。本文通过分析可监测形变梯度、MAI干涉图标准差与相干性系数的关系，描述了以相干阈值为基准的大梯度形变区域识别流程。以山东济宁某煤矿为研究区，利用上述方法对该区域进行三维地表形变监测，获取包含大梯度、失相干区域在内的三维地表形变场。研究结果表明，该方法可以突破InSAR技术受失相干、可监测形变梯度小等因素限制的问题，实现大梯度、失相干三维地表形变监测。

(3) 针对轨道、地形和噪声等误差对多视线向D-InSAR三维地表形变解算的影响，提出了一种基于GPS约束定权的多视线向D-InSAR三维地表形变解算方法。该方法以同期高精度GPS形变观测量对多视线向D-InSAR观测量进行精度评价，顾及轨道、地形和噪声等误差对多视线向D-InSAR观测量影响的特点，建立与像元位置和高程相关的联合拟合模型，完成逐像元精度评价，实现对多视线向D-InSAR观测量定权。本文以分别模拟和真实数据试验对该定权方法的可行性和精度进行验证。结果表明，所提出的解算方法提高了解算精度，解决了多源InSAR数据联合解算权矩阵确定中存在的定权难、适用性差、精度低等问题。

(4) 针对粗差对多源InSAR数据三维地表形变解算的影响，基于Helmert方差分量估计，发展了多源InSAR数据三维地表形变抗差解算方法。该方法利用最小二乘方法进行三维地表形变解算，获取预平差观测量的残差，建立最小二乘残差与观测量单位权方差的函数关系，并通过计算出的单位权方差对InSAR地表形变观测量进行定权。根据选权迭代法等价权原理对最小二乘残差的判断，进而对观测量进行保权、降权和淘汰处理，实现三维地表形变的抗差解算。论文采用模拟试验数据和巴姆地震实测数据对所提方法的可行性和精度进行了验证。结果表明，上述抗差解算方法可以有效消除粗差的影响。