利用PS-InSAR技术获取长时间地表形变历史.pdfVIP

其他空间技术地震应用专题 利用PS．InSAR技术获取长时间地表形变历史 宋小刚，单新建，屈春燕 (中国地震局地质研究所。“地震动力学”国家重点实验室北京100029) 传统D-InSAR技术进行长时间地表微小形变监测中，往往受到时问和几何去相干以及大 气干扰等因素的影响。研究人员对大量SAR图像的研究发现，在长时间范围内图像上仍然存 在相位和幅度变化稳定的点，PS～InSAR技术能够利用这些稳定点上的相位特征，实现长时 间尺度上的地表形变分析。本文利用PS—InSAR技术对某个研究区域内横跨1992到2003年 的56幅影像进行了PS-InSAR处理，获取了该地区l万多个“稳定点”的形变历史。通过对 某些点的形变过程分析，我们发现有一些偏离线性形变速率曲线太大的值可能是相位解缠错 误的结果，也可能是有大的非线性成分。而对于有的值，如果加上或减去2．8cm(半个波长)， 那么它就会很好的符合在线性曲线或非线性曲线上，这些点就有可能是解缠误差；我们还可 以观察到一些特定的点只有在一段时间内具有高质量，特别是在时间序列的开始观察到了一 些点的质量不够好，它对应了一些点在某个特定的时间点上才开始变为一个永久性点散射 体，而另外一些点在时间序列的后期出现了质量降低，它说明了该点只有在某个特定的时间 之前是一个永久散射体，基础设施建设和修改可以解释这种情形。 另外，对整个PS-InSAR流程中的几个关键步骤我们采用了特殊处理，以保证结果的可 取的形变才是可靠的。本文把PS点的筛选贯穿于整个PS—InSAR处理流程中，确保选取PS 点的准确性。首先基于单个SLC数据的低光谱相位变化和低强度变化来识别PS候选点列， 然后在后续一次次的回归分析，一次次的模型精化中，利用回归分析的相位标准差作为质量 测度，来探测和拒绝那些不适合做PS分析的点。而后在分离掉大气和非线性形变后，用噪 声相位计算了时间维上的相干性，然后利用时间相干性作为标准，再一次去除低质量的点。 最后把精确求解的模型扩展到更多的点上，乃至整个影像区域上，来识别遗漏的PS点。如 果其他的点可以为它们找到解，那么在现有解的基础上去检查其他另外的这些点是有可能 的，对于另外的点，一个已经接受的点可以在积分步骤中被用作局部参考点，而使用一个局 部参考点有低大气和非线性形变相位的优势。2)采用分块回归分析的方法。时间上的差分 干涉相位分析是PS分析中的重要一环。更明确的是，我们分析的是点目标对之间差分干涉 相位的空间差。先建立观测相位与垂直基线和时间的线性模型，然后进行回归分析。对于原 始的差分干涉图堆栈，回归分析只有在那些距离比较近的点对之间进行结果才可靠。因为随 着点对距离的增加，大气相位成分、与残余轨道误差趋势相关的基线误差、相对比较高的形 变速率会导致单点偏离同归面。在一定的噪声水平下，对于缠绕相位，回归分析不再能够可 靠地求解。所以，对初始的差分干涉图，一定要分块同归。每一块选一参考点，对该块内的 其他点做回归分析。最后在各块参考点之间做回归分析。3)时间解缠与空问解缠相结合。 对回归分析中解缠不正确的干涉图(即时间维解缠)，采用空间解缠的方法对残余相位进行 解缠，来改正相位模糊度误差。4)基线精化模型是基于干涉图数据本身。是通过最小二乘 意义上模型相位和计算相位之差的最小值来精化基线。这个结果不是从物理意义上去改正基 线，因为其他一些大尺度的相位趋势(像大气路径延迟)也通过基线精化被部分地补偿了。 5)模型的多次精化。在基线精化后，大气相位去除后以及对差分干涉图进行时间上的滤波 后都要重新进行回归分析，然后更新模型，使模犁在一次次精化中得到完善。 在分析中我们发现，在一些植被覆盖密集的地方，没有PS点存在，需要布设人工角反 射体来弥补这样的不足。所以PS—InSAR和CR—InSAR的结合是我们下一步研究的内容。 43