融合DInSAR和MAI技术反演九寨沟地震三维形变场

宫熙雯,张正加,陈启浩,刘修国

(中国地质大学(武汉)，湖北 武汉 430070)

2017年8月8日21:19:46，四川省北部阿坝州九寨沟县发生7.0级地震(33.2°N，103.82°E)，震源深度20 km，截至2017年8月18日09:00共记录到地震总数为5 020个(余震总数5 019个)，其中最大余震为：2017年8月9日10:17四川阿坝州九寨沟县4.8级。此次地震，中国地震局徐锡伟等专家判断这次地震与巴颜喀拉块体南东向推挤有关系，发震构造推测为塔藏断裂南侧分支和虎牙断裂北段[1]。

一些学者对该地震进行了相关研究，郑绪君等采用迭代反褶积和叠加方法反演得到破裂模型，用于确定地震的破裂过程[2]，汪建军等采用库仑应力模型研究地震激发的库仑应力变化[3]，文献[4]采用多步定位法对九寨沟地震Ms7.0的主震及其余震进行重新定位，提出震源机制解决方案，并通过CAP法(Cut and Paste Method)反演了九寨沟地震的质心深度。除了对以上地震参数的反演，获取地震同震形变场对揭示地表形变、理解和分析发震构造具有重要意义。由于差分雷达干涉测量(DInSAR)可以获取全天时全天候、大范围、连续空间覆盖、高精度的地震区域同震形变图[5]，在地震监测方面应用广泛，对于本次九寨沟Ms7.0地震，赵强等利用Sentinel-1A和RadarSat-2升轨数据结合两轨法DInSAR反演了地震的同震形变场[6]；熊威等使用Sentinel-1A数据进行同震形变场反演，并确定形变区的相干系数[7]；申文豪等采用DInSAR方法利用升、降轨Sentinel-1数据反演了同震形变场[8]。Zhao基于InSAR和GPS测量，反演地震同震滑动分布和倾角，并获得地震视线向的同震形变场[9]。但由于DInSAR方法的设计本质是一种侧视雷达测量技术，得到的并不是地表真实形变，而是地表真实形变矢量通过一系列坐标系转换和几何转换后在视线向的投影[10]。……