基于ALOS PALSAR数据进行雷达差分干涉提取玉树地震同震地表形变场

本文针对2010年4月14日玉树地震引起的地表形变,使用日本ALOS卫星PALSAR L波段雷达影像数据,应用两轨雷达差分干涉(DInSAR)处理得到了以玉树为中心11 000km2范围内的同震形变场,空间分辨率为8m,并在此基础上对玉树地震的震源机制和发震机理进行了分析。研究结果表明L波段雷达数据适合在地形起伏较大的地区进行DInSAR形变探测。该同震形变场信息可为玉树地震的同震形变反演提供参考数据。该研究进一步证实DInSAR技术在大规模地表形变探测和地学研究领域具有广阔的应用前景。     

2016年意大利阿马特里切Mw 6.2地震震源机制InSAR反演

2016年8月24日，意大利中部阿马特里切（Amatrice）地区发生Mw 6.2地震。采用ALOS-2条带模式和SENTINEL-1A宽幅模式的合成孔径雷达（synthetic aperture radar，SAR）数据分别进行SAR差分干涉测量处理，获取了该地震的同震形变场。结果显示，本次地震造成意大利中部地区发生明显的地壳形变，在雷达视线向最大沉降量达19.6 cm。基于合成孔径雷达干涉测量（interferometry synthetic aperture radar，InSAR）和GPS同震形变场数据对此次地震的发震断层进行联合反演，通过改进倾角和平滑系数获取方法，得到了最优滑动分布模型。通过使用单断层模型和双断层模型进行反演可知，双断层模型反演结果优于单断层反演结果，两种模型下反演模型相关系数分别为0.85和0.89，发震断层走向分别为160°和158°，倾角分别为44°和46°，倾滑分布主要位于地下5~7 km，平均倾滑角为-80°，最大倾滑量0.9 m位于地壳深度5 km处，该发震断层是亚平宁冲断带的一部分，为NW-SE向延伸的正断层，断层长约20 km。综合使用地震同震形变场和GPS数据对震源机制进行反演、模拟和分析，获取了高精度的震源参数，可以为分析地震危险性和断层破裂参数等提供数据支持。     

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

差分雷达干涉测量(DInSAR)在地震形变监测中具有重要作用,但DInSAR仅能获取雷达视线向的形变量,无法准确反映地表三维形变.为了获取真实的地表形变,文中采用一种联合DInSAR和多孔径干涉技术(MAI)的地震三维形变场反演方法,以2017年九寨沟Ms7.0级地震为例,融合升降轨Sentinel-1A数据,联合解算获取研究区地表在竖直、南北、东西向上的形变,构建地震三维同震形变场.实验结果表明,本次地震的断裂带为西北—东南(NW-SE)走向,地震断裂面西北侧(下降盘)错动位移较大,向东南偏东方向滑动;东南侧(上升盘)错动位移较小且3个方向位移的贡献较为平均,大致向正西北方向滑动.地表运动以水平形变为主,纵向起伏较小,东西方向对本次地震的同震形变贡献最大,存在向东最大0.25 m、向西最大0.12 m的形变,断层性质为左旋走滑断层.     

青海玉树地震差分干涉雷达同震形变测量

2010-04-14青海玉树发生7.1级地震后,作者利用震前和震后获取的日本ALOS卫星PALSAR遥感数据,开展了差分干涉雷达(D-InSAR)地震同震形变测量与分析。结果表明:玉树地震引起较大范围地表变形,地震变形沿玉树—甘孜断裂带向南东东方向扩展,在N33.7°,E96.81°附近达到最大形变量,D-InSAR监测到雷达视向上的最大形变量为35cm。地表形变特征对于评价玉树地震破坏程度、推断断层性质、研究地震形变和地震孕育特征具有重要的参考价值。     

基于GPS观测数据的汶川地震断层形变反演分析

通过高精度数据处理方法获取了2008年5月12日汶川地震(Mw=8.O)可靠的同震形变位移,较大形变主要在汶川、青川和北川县,其中最大同震形变不在震源处,而在北川县,水平方向达到2.3 m,垂向为0.7 m.并进一步采用分布滑动模型,利用弹性半空间均匀位错理论反演得到同震断层滑动分布.研究结果表明,不约束断层滑动特性的反演结果显示汶川地震断层不仅仅是右旋、逆断层,还伴随着少许左旋走滑及正断层变化.其地震矩张量为2.38X1021N·M,基于反演的断层滑动分布正演的区域同震形变场主要变形区域分别在都江堰-汶川县,北川县,青川县等一带,地表最大变形达到东方向4.5m,北方向0.8m,具有极强的右旋走势,均与野外地质调查结果相符合.     

基于曲面断层结构的2007年所罗门群岛地震同震滑动分布反演

2007年4月1日,南太平洋岛国所罗门群岛的俯冲带地区发生了Mw 8.1级地震.利用日本宇宙航空研究开发机构(Japan Aerospace Exploration Agency,JAXA)的ALOS/PALSAR卫星雷达影像数据,采用二通差分干涉技术获得了覆盖所罗门群岛地震震区的同震地表形变场,根据四叉树方法对原始观测数据进行降采样,利用弹性半空间三角形位错模型进行断层几何参数反演;并结合断层深度剖分图的特征,构建更为合理的上下两个不同倾角连接的断层面模型,通过提出的断层剖分技术确定了最佳同震滑动分布模型.为了测试断层倾角的最佳估计,对不同断层倾角进行测试,对比合成孔径雷达干涉测量(interferometic synthetic aperture radar,InSAR)技术获取的同震观测数据和不同模型预测的InSAR干涉图,发现上下断层面倾角分别为10°和23°的模型更为合理.据此模型得到该断层同震滑动分布主要发生在距地表40 km深度范围内;最大同震滑移量为5.227 m,深度为2 km.     

门源Mw5.9级地震形变InSAR观测及区域断裂带深部几何形态

利用Sentinel-1A卫星升轨、降轨合成孔径雷达影像数据，提取了2016年门源Mw5.9级地震的高精度合成孔径雷达干涉同震形变场，利用单纯形法和非负最小二乘法反演确定了地震断层几何和滑动分布，并构建了区域断裂带的深部几何形态模型。结果表明，门源Mw5.9级地震同震形变以地表抬升为主，沿升轨、降轨视线向的最大值分别为5.3 cm、7.1 cm；地震断层走向、倾角分别为133°、43°；地震滑动以逆冲为主，主要发生在地下6.14~12.28 km处，最大滑动量约0.5 m，平均滑动角为66.85°，地震矩为1.0×1018 N·m（Mw5.94）；形变观测拟合残差均方根为0.36 cm；区域断裂带的深部几何形态以花状构造为特征，整体倾向南西，门源地震发震断裂为花状构造中未出露地表的盲断层。相关成果能够为研究区域地壳运动与变形、活动断裂与地震孕育发生等提供参考。     

基于D-InSAR的长宁地震形变场提取与模拟

2019年6月19日四川长宁县发生Ms6.0地震，提取地震同震形变场，反演震源机制，对于地震破裂分析、指导救灾具有重要意义。本文利用两轨差分干涉处理覆盖长宁地震影响区域的两景Sentinel-1A影像，在对D-InSAR关键技术和影像特征研究的基础上，配置处理方法和参数，提取了此次地震的同震形变场，采用Okada模型和正向建模对发震断层的几何参数及形变场进行了反演和模拟。结果表明，长宁地震形变场呈西北—东南走向的次级断层控制的不规则椭圆形，断层的两侧区域特征差异明显，断层左下侧为沉降区，右上侧为隆升区，两者的最大视线向形变分别为8、6 cm。断层的运动主要以左旋走滑为主，平均滑动距离约0.38 m，平均滑动角约55°，正向模拟的形变场与观测结果相符，这表明观测结果较可靠，同时也提高了低相干形变区的观测精度。     

利用PALSAR数据获取拉奎拉地震同震形变场

2009年4月6日发生的意大利拉奎拉Mw6.3地震造成了地表的严重破坏,为了解地震引起的地表变化情况,本文利用ALOS/PALSAR影像数据,采用二轨差分干涉测量技术提取了该次地震的地表同震形变场。通过分析可知:地表在西北-东南方向上发生错动,形变主要发生在40×30 km的范围内,西南部出现了沿视线方向的下降,东北部出现了沿视线方向的抬升,最大值分别为0.337 m和0.122 m。结果表明,获得的同震形变场与地震地质调查的结果一致。     

升降轨InSAR数据约束下的2007年阿里地震反演分析

2007年5月5日,西藏阿里地区发生了Mw6.1级地震。本文采用Envisat卫星的升、降轨SAR数据获取地震的精确同震地表位移,然后基于弹性半空间位错模型,分别采用MPSO非线性和最小二乘线性反演算法确定了断层几何参数和滑动分布模型。结果表明,分布式滑动分布模型能较好地解释观测到的InSAR地表形变场。地震发震断层为走向158°、倾角43°的西南倾断层,主要的滑动量集中在7~12km的深度,最大滑动量约0.3m,位于9km的深度。反演给出的地震矩为1.24×1018 Nm,与地震学结果相一致。