基于d-insar的地表干涉形变场时空演化特征与震源破裂面关闭关系研究.docxVIP

基于d-insar的地表干涉形变场时空演化特征与震源破裂面关闭关系研究 0 震前地表变形场的演化过程 地震的发生和发生是由于地应力的长期积累，岩石突然失稳，能量迅速释放。能量释放后, 又会出现一定量的弹性恢复。从震前的能量积累到地震时的同震位移, 再到震后的弹性恢复, 在发震断层附近都会有不同形式和不同特征的地表形变场展布。如果我们能得到震前、同震和震后地表形变场的动态演化图像, 将对断裂带的变形演化特征、孕震过程、孕震机理、地震预测等问题有更深入的认识。 许多震例表明, 孕震区形变场在震前1～3a内有加速变形的过程, 形变量级可达每年几十毫米。例如, 位于1966年3月邢台地震极震区的CRC449水准测量点, 从1964年开始, 下降速率急剧增大, 年变化率达到数十至100mm;中强地震发生时, 其同震位错一般在几十厘米到几米之间。例如, 1998年张北地震的同震隆起变形量达25cm (Wanget al., 2001) 。玛尼地震的同震左旋水平位错达到了7.96m (Shanet al., 2004) ;地震发生后, 在一段时间内震源区又会出现地表形变微量回弹, 该量级一般在几厘米至几十厘米之间。例如,Peltzer等 (1996)收集了Landers地震后3a内的ERS-1/2 SAR资料, 进行D-InSAR处理。结果表明, 地震造成的地表破裂隆起区有大约2cm的下降回弹, 地表破裂的拉张沉降区大约有6cm的上升回弹现象。Reilinger等 (2000)利用GPS和D-InSAR测量数据研究了1999年8月土耳其7.5级地震, 结果表明同震位错量达5.7m。震后75d内, 断层恢复滑动量为0.43m。2002年Jacobs等利用D-InSAR获取的1999年美国加利福尼亚洲Hector Mine地震后4～354d内的4期干涉图像, 显示在0.5～50km的空间尺度上, 有几厘米的震后弹性恢复。这些震例表明强震前后的形变量级完全在D-InSAR的观测精度范围之内。 目前看来, 国内外学者只是研究同震和震后弹性恢复的形变场图像, 还没有用D-InSAR技术研究震前干涉形变场图像。若我们能通过1个震例, 观测到震前地表形变场的动态演化图像, 则将有助于对地震孕震过程、孕震机理、地震预测以及地震重复间隔等问题的深入认识。本文以1997年11月8日西藏玛尼MS7.9地震为例, 利用12景欧空局雷达卫星 (ERS-1/2) SAR资料, 采用D-InSAR技术获取了玛尼地震的震前干涉形变场演化特征。通过分析, 发现地震前地表形变场有明显的变化特征, 这种具有前兆意义的形变场变化, 可能为地震监测预报提供一定的科学基础。 1 地震地表变形场及断裂分布 1997年11月8日, 在西藏自治区那曲地区玛尼乡北约150km处发生了MS7.9强烈地震 (35.2°N, 87.3°E) 。玛尼地震发生在玛尔盖茶卡断层附近, 该断裂为羌塘块体北部边界, 沿该断层带屡次发生6级以上强震 (汪一鹏, 1998) 。该地震处于藏北无人区, 方圆几百公里内无形变观测台, 这些情况限制了我们对孕震区地表形变场、活动断层变形等问题的深入认识。而星载D-InSAR测量技术将不受这些因素的限制, 在了解空间形变场演化方面具有不可替代的优势。 从图1可以看出, 几个震源机制解结果非常接近 (许力生等, 1999;高原等, 2001) , 两节面分别沿NEE向和NNW向分布, 且节面较陡。美国国家地震信息中心记录的余震也具有沿近EW向分布的特征, 且均匀地散落在整个断层附近, 主震位于断层东部。据野外考察资料 (徐锡伟等, 2000) , 玛尼地震的地表同震破裂带全长约120km, 沿NEE向分布, 以左旋剪切为主, 最大左旋位移量达4.5m, 其东西两侧水平位移量约2～3m, 地表形变相当剧烈。由此看来, 玛尼地震的发震断层面为NEE—近EW方向, 断层性质主要为左旋走滑。 2 sar数据处理 研究区选在玛尼地震震中区100km×200km的范围内。通过对资料实际调研, 发现研究区恰好位于2889和2907两景区之间, 因此我们必须收集到2个景区的SAR干涉数据, 才能较为完整地解释玛尼地震的形变场。我们收集了1995年1月至2000年12月的16景欧空局ERS-1/2雷达卫星数据。通过采用“三通”、“四通”差分干涉模式, 对SAR数据进行了处理, 即由震前相差1d的两景SAR图像 (Tandem数据) 产生干涉图, 通过解缠获取高精度DEM, 来消除另一幅干涉图的地形影响, 以达到减少地形残留相位, 提高形变量观测精度的目的。由于两幅SAR图像是在相邻轨道上获取的, 入射角、星下平台方向和线采样距等卫星参数有微小的差异, 为保证干涉处理的精度, 必须统一成像