应用差分InSAR技术获取强震区地壳形变信息研究.pdfVIP

应用差分InSAR技术获取 强震区地壳形变信息 张景发 (中园地震局地壳应力研究所) 摘要：本文介绍了憋分InSAR处理的主要工作步骤和关键技术，对发生过7．5级地震的西藏玛 尼地送SAR数据进行了处瑷，分剃获得了与遣形有关的干涉翻和与潮震形交有关静差分予涉图 像。 关键谲：差分InSAR干涉图像地壳形变 一、引 言 InSAR是同一类型雷达卫星通过两次平行路径接收地面同一景观的回波信号而生成的复图 像对，最终形成干涉图像。此干涉图像与地表地形有关。差分干涉雷达是由两幅对应的干涉 图像推导崽来的，差分于涉图像熊反映接收时闽间隔内的地壳形变，薮测量糖瘦为cm～mm 级。InsAR技术在地形测量、滑坡监测等领域都已有成功的应用，在地震监测及震害评估中 也具有良好酶应恶前景，剩掰InSAR技术获得强震区不露时阉段的形变信崽，对于地震翡孕 育一发震一震后蠕变等地壳变形调整的全过程的认识，特别是地震地壳变形前兆信息的提取、地 震断层的运动学及其力学成因的研究具有重要的意义。 地震震中位置为北纬35．1。，东经87．3。。在雷达图像上，震中区地震破裂清晰可见，还可见 到水系有锩动瑰象，显然，地震破裂沿活动叛层带发生。玛尼强震区发生在青藏高原区，人 理试验与研究的理想地点。 二、差分InSAR处理的工作步骤和关键技术 差分InSAR的主要工作步骤惫括：像对选择、复匿像对次像元级精度酶酝准、干涉圈像 的生成、干涉图像的平地校正、全相位恢复、差分像对确定及参数调整、差分干涉图计算等。 这些步骤中，次像元精度的酝准、褶干性分析、全相位恢复等是其中的关键技术。 1．次像元精度匹配 对复图像对通常采用区块匹配的方法来确定像对对应共轭像点，在区块瓯配中又采用互 相关法来寻找对应共辍像点的坐标。杰攘关系数矩阵获得的共轭像点位置为整数基，但是真 正的最大相关位置不一定在整数位置，因此要利用相关系数矩阵上以峰值为中心的局部矩阵， ·129· 求出次象元精度的坐标位移，通常是在z、Y方向上各以二次曲线来回归求取其位置，二维定 位可以一维方式解决，不同之处在予各方向以二次曲线求取鞍部位置，两条曲线交点即为所 隶懿数值。 2．相干性 相位差与距离差成正比，即如果信号与相位有固定的对应关系则信号是棚予的。当两幅 原始影像嘉度稿于辩，稳邻像元酶籀位变化是连续觞，表臻在干涉褶位上裂为当其穗位差一 致时可视为有较高的相干性。影响相干性的因素有：①数据获取期间的气候藏异；②处理中 相位有误差；③在两次数据获取之间圈标相位有变化；④观测位蚤有轻微变动。基线长度是 将大大影响相位和高程测量的精度。在连续的相位还原中可以作为安全路径的评估。设G，、c。 为两像对相位，E(…]为窗口内的信号总和，则确定相于性的方法为： y— l墨[坠丕g；!l f1、 √E[c，×C。]E[c。×c。] 3。二维全相饿恢复 干涉图像包含了SAR像对的振幅及相位信息。在图像中每一点都要计算对应的平均振幅 和相位差。干涉图像相位差直接与地形有关，仅以2丌的模数来标定。要计算每一点高程，右： 每一次期位测量时要加上梗位愿期的整数。粲决这个多勰翳。2嚣闽题霹称为全相位恢复。 (1)相位的奇点 所谓的奇点，位于真实像点的间隔处。可先以每4点最小回圈检验全相位图，将奇点的 位置标描。检验时、滔顺时针隧鬻趣藏英波数，余l者为委奄点，余受者为副奇点。正奇点使 得任一顺时针绕过它的回圈波数都会多1，若是逆时针绕过它则少1；负奇点则相反。． 大量的奇点出现在水面、山区，由于相位对地面上应为连续变化，此类聪域上的波程差 与遣面点、卫星阕的尼蔼关系跫遭到破坏。 (2)还原路径的阻隔墙 每两个正、负奇点间加上一道阻隔墙，禁止由此通过，便可确保全幅相位还原不矛霜。但 壶予疆隔墙使得鼹边的柱位不连续，数阻隔墙的总长越短越好。