利用升降轨SAR数据提取地震三维同震形变场模型.pdf

利用升降轨SAR 数据提取地震三维同震形变场模型的完善 与实现 摘 要 合成孔径雷达干涉测量技术（Interferometry Synthetic Aperture Radar ，InSAR ） 在近二十年来得到了迅猛发展，在获取地面形变信息具有独特的优势，使得人类 从外空间对全球地表进行长时间大范围的监测成为可能。但是，InSAR 测量获取 的只是地表形变沿雷达视线向（LOS ）的分量，不能全面地反映地表真实变形， 即存在 InSAR 的视线向模糊问题。尤其是对于地震、火山喷发等大形变，需提 取地表的三维形变 （即竖直向、东西向和南北向形变）特征，以迅速获取灾害信 息、提供救援以及更加准确分析发震构造与地震机理。 目前利用InSAR 提取三维形变的方法有许多，主要方法包括以下几种： （1）方位向偏移+距离向偏移方法：该方法主要利用基于像素偏移技术获取 地表沿卫星轨道方位向上的位移和距离向的偏移，最终解算出地表位移的三维分 量。 （2 ）常规 D-InSAR+方位向偏移方法：该方法利用像素偏移技术获取方位 向的形变，再结合D-InSAR 获取的视线向的形变，根据空间几何关系，可以解 算出地表位移的三分量。 （3 ）常规D-InSAR+MAI 方法：针对使用Michel 提出的方位向像素偏移法 （AZO ）获取沿轨道方向位移分量时精度不高的问题，Bechor 提出了MAI 方法。 这种方法是基于分光束 InSAR 进行数据处理，生成前向和后向干涉图，根据两 幅干涉图上的相位差，获取沿轨道方向的位移分量。进而结合D-InSAR 获取的 视线向形变进行三维解算，获得地表位移的三个分量。 （4 ）左右视InSAR 技术：2004 年，Tim Wright 提出了使用多幅不同不同入 射角生成干涉图的方法。采用升轨左视数据、升轨右视数据和降轨左视数据、降 轨右视数据，得到四个干涉图，通过四个干涉图得到地表位移沿视线向的四个方 程，再通过最大似然法求解出地表位移的三个分量。 通过上述四种方法的介绍，可以看出，为了求得地表位移三个分量，必须有 至少三个方程组成方程组才能进行解算。而常规的 D-InSAR 技术获得的干涉图 只能提取距离向的位移，只有一个方程，即使利用升降轨数据也只能获得两个方 程；而且同一个区域有三个不同入射角、方位角且角度差较大的数据是很少见的， 所以比较常用的方法是获得方位向的偏移量，根据不同的组合可以得到至少三个 方程，进而求解地表位移三分量。 获取方位向位移的方法，有方位向像素偏移法和MAI 方法。利用MAI 方法 获得方位向的偏移的精度要明显高于方位向像素偏移法，而且与D-InSAR 处理 精度在一个数量级，进而结合 D-InSAR 技术获取的距离向位移，由几何关系， 获得四个方程，再利用最小二乘法求解出地表位移的三个分量。 在以上方法模型的基础上，本文对于右视升降轨、左视升降轨情进行了分别 讨论，得出了两点结论：①北、东方向上的位移对方位向上的形变贡献和雷达的 左视右视无关；②竖直方向上的位移对于视线向上的形变贡献相同，而南北和东 西方向上的位移对于视线向上形变的贡献则相反。针对左、右视的不同情况，提 出了不同的求解公式，进而对左、右视两种情况进行了统一，完善了利用升降轨 SAR 数据获取三维形变场的计算模型。 本文利用完善后的模型，选取意大利拉奎拉地震进行试验，选取震前升降轨、 震后升降轨ASAR 数据各两景，采用D-InSAR 的方法分别获得升轨、降轨数据 的同震形变场。进而，采用MAI 方法获得升、降轨图像对的在方位向上的干涉 条纹图，对干涉条纹图进行相位解缠，解缠方法采用区域增长法，并将解缠后的 相位图转化得到方位向上的形变图。 图1 三维同震形变场 运用完善后的计算模型，联合视线向和方位向形变获取拉奎拉地震的三维同 震形变场，结果如图1 所示。结果显示，（1）图1 中黑色实线位置为地表破裂区 即断层所在的地表位置，从三维形变图可以明显的看出断层两侧发生了位移，且 通过模型计算得到的地表破裂区和断层位置较为吻合，走向为西北—东南 （NW-SE）。（2 ）断层面西南侧在竖直方向上发生了明显的沉降运动，升降轨共 同覆盖区域内最大下降量为0.16m ，而在东北侧竖直方向上表现