基于InSAR和MAI的电离层误差校正及同震三维形变场计算与断层滑动反演.pdfVIP

和MAI的电离层误差校正及同震三维形变场计算与断

层滑动反演

地震是地球上破坏力最强的自然灾害之一,我国境内地震活动非常频繁。精

确测量地震导致的同震形变场,是了解和认识地震发生机制最为直接的手段。

结合地表形变观测和地球物理模型模拟,对揭示区域断层活动规律和评估中

远期地震风险具有重要意义。作为一种新的空间大地测量技术,合成孔径雷达差

分干涉测量(DifferentialInterferometricSyntheticApertureRadar,

DInSAR)自上世纪90年代诞生以来,已在全球地震形变测量中展现出了无可替

代的地位。

得益于其全天时和全天候成像观测能力、广覆盖范围和高分辨率等优

势,DInSAR已被应用于全球一百多个地震同震形变场的测量和模拟。然而,针对

地震导致的复杂形变场,DInSAR仍面临着局部地区相位失相干严重、大气延迟长

波信号干扰以及单轨观测难以完整揭示同震形变特征等局限性。

靠近发震断层附近的同震形变往往具有较大的形变梯度,故容易使DInSAR

测量结果由于相位失相干而产生空值。虽然利用长波段SAR数据(如L波段)恢复

大梯度形变具有独特优势,但其容易受大气电离层干扰影响,并可能导致数十厘

米的形变测量误差。

此外,DInSAR测量的地表形变其实为地表真实位移场沿卫星视线向(Lineof

Sight,LOS)投影值,因而采用单一干涉对往往无法对复杂特征的同震形变场进

行准确解译。DInSAR的这些局限性阻碍了对地震形变机制准确认识,也进一步限

制了断层参数的准确反演。

为克服DInSAR技术在同震形变测量中面临的这些局限性,本文系统研究了

方位向形变的多孔径雷达干涉技术(Multi-ApertureInSAR,MAI),并基

于MAI测量的方位向偏移研究了同震SAR干涉对电离层干扰误差校正方法,旨在

扩展基于SAR影像的多维地表形变测量能力和提高同震形变场测量精度。本文还

提出了新的地震同震三维(Three-dimensional,3D)形变场计算方法,并结合实际

震例进行了3D同震形变测量和反演,以验证本文所提方法的有效性和优势。

MAI是基于一对SAR影像进行分孔径干涉处理以获取地表方位形变的技术,

可以有效弥补常规DInSAR仅对地表视线向形变敏感的缺陷。本文系统研究了

SAR分孔径成像原理、MAI干涉模型和数据处理流程,并定量评估了MAI的测量

精度。

以2016年Mw7.1级日本熊本地震为例,MAI测量结果显示该地震导致地表产

生长约45km的地表破裂,并且由南北两盘的形变方向可清晰判断其为右旋走滑

型地震,而这些信息无法从常规DInSAR测量结果反映出来。分析表明,MAI测量

精度与相位失相干程度密切相关,并且与SAR卫星系统参数和数据处理配置参数

有关。

理论计算显示,当SAR相干系数为0.8时,基于PALSAR的MAI测量精度优于

0.08m。针对同震SAR干涉图中的电离层干扰误差,本文提出基于SAR方位偏移

量和断层模型的校正方法。

为了获取高分辨SAR方位向偏移场,提出融合MAI和像素偏移追踪(Pixel

OffsetTracking,POT)两种方法测量的方位向偏移思路,即综合利用SAR相位和

振幅信息实现两者的优势互补。为了获取电离层干扰引起的方位向偏移,提出利

用已有断层模型去除SAR方位向偏移中的地表形变贡献,从而改善电离层误差校

正精度。

轨道PALSA20080531同震干涉对为例,进行了电

离层误差校正实验。结果显示该干涉对中电离层相位最大贡献量约为17rad,对

应地表约0.32m的LOS向形变。

校正结果成功恢复了T471轨道断层南盘的三个连续干涉条纹。计算显示,

远场三个区域内干涉相位的标准差在校正后分别降低了50%、65%和84%。

针对基于DInSAR一维视线向形变难以全面揭示同震形变特征,本文提出联

合DInSAR、MAI和形变梯度张量模型(DisplacementGr