2022年青海门源地震震源机制与同震滑动分布研究

采用Sentinel-1A卫星提供的升降轨雷达影像数据研究门源地震同震形变及震源机制。首先利用合成孔径雷达差分干涉(D-InSAR)技术获得门源地震的同震形变场;然后以升降轨同震形变场为源数据,利用弹性半空间位错模型进行反演,确定地震的断层几何参数和滑动分布;最后基于同震滑动模型对升降轨同震形变场进行正演。结果表明,沿雷达视线方向的升轨和降轨同震形变场最大抬升形变量分别为39 cm和58 cm,最大沉降形变量均为56 cm。此次门源地震为左旋走滑型地震事件,发震断层方向为NWW-SEE、走向为109°、倾角为86°,主要集中在地下2～6 km处,最大滑动量为4.2 m,释放的地震距为8.22×10¹⁸Nm(M_W6.6)。正演结果表明,本文滑动分布模型准确可靠。     

2021年云南漾濞地震InSAR监测及断层滑动分布反演

采用双轨D-InSAR技术获取2021-05-21云南漾濞6.1级地震的升降轨同震形变场,并以此为约束,基于Okada均匀弹性半空间位错模型反演发震断层的几何参数与同震滑动分布,最后通过正演模型预测地震的同震三维形变场,分析地震的形变特征及发震断层的几何形态和构造运动。结果显示,在InSAR形变场中,升轨LOS向最大形变量约10 cm,降轨LOS向最大形变量约13 cm;利用Okada模型反演得出,发震断层长约12.53 km,宽约6.44 km,位于地下3～9 km处,断层走向约136.34°,倾角约83.79°,滑动角约-175.46°;断层的同震滑动主要分布在沿走向9～21 km与沿倾向向下3～9 km范围内,最大滑动位移出现在地下4 km深度处,滑移量达到0.64 m;通过正演模型的地表三维形变预测发现,垂直方向最大抬升约3.22 cm,最大沉降约4.04 cm。总体而言,此次地震的发震断层类型为右旋走滑,且带有正断分量,初步判断地震发生在一条未知的NW-SE向断裂带上,该断裂可能属于维西-乔后-巍山断裂和红河断裂北段的次级或分支断裂。     

Sentinel-1 SAR数据约束的2018年谢通门MW5.8地震同震破裂模型

通过Sentinel-1卫星升降轨数据获取谢通门地震的同震形变场,并基于均匀弹性半无限位错模型反演地震的同震滑动分布模型。InSAR同震形变场表明,升降轨视线向最大形变量分别为0.049 m和0.051 m,形变场长轴大致呈南北方向,位于甲岗-定结断裂西侧。通过对倾角和倾向进行格网搜索发现,西倾节面更可能为该地震的发震节面。反演结果表明,滑动分布主要位于2~10 km深度范围内,平均滑动量为0.02 m,最大滑动量为0.10 m,发震断层倾角为47°,平均滑动角为-81.60°,显示该地震以正倾滑动为主。大地测量数据约束的该地震震中为30.27°N、87.75°E,震源深度为6.58 km,释放地震矩为5.056×10¹⁷ Nm,对应矩震级为M_W5.7,与GCMT、USGS公布的震级基本一致。综合分析震中位置和滑动机制认为,甲岗-定结断裂的分支断层为本次谢通门地震的发震断层。     

基于升降轨Sentinel-1 SAR影像研究精河MS6.6地震震源机制

采用Sentinel-1A/B卫星升降轨SAR数据获取了2017-08-09新疆精河MS6.6地震的同震形变场，并以InSAR形变场为约束，利用均匀滑动模型和分布式滑动模型反演了发震断层的滑动分布。结果表明，升降轨InSAR形变场均显示为隆升，形变影响范围约30 km × 40 km；发震断层走向约为83.8°，倾角约为40.6°，滑动角约为89.1°，震源深度约为20.7 km，矩震级为MW6.35，同震位错以逆冲运动为主兼有极少量左旋走滑分量。精河地震的发震断层为库松木楔克山前断裂，该次地震是新疆北天山地区逆冲断裂带深部发生错动的结果。     

汶川地震断层滑动模型反演及地震区域构造的大地测量约束

以汶川地震前后GPS观测得到的同震形变场为约束,利用带约束的最小二乘反演方法,分别反演考虑与不考虑滑脱层的最优同震断层滑动模型,并在此基础上正演了同震重力变化。结果表明,含滑脱层的断层滑动模型模拟的GPS台站的同震水平形变和重力台站的同震重力变化都更接近实际的观测结果,证实龙门山的中央断裂带南段存在大规模的近水平的滑脱层参与了汶川地震的同震滑动。     

九寨沟Ms7.0地震的InSAR观测及断层破裂研究

为提高九寨沟2017-08-08 MS7.0地震的InSAR同震形变场精度，厘清发震断层的构造形态及形成机制，对形变区相干系数进行统计分析并确定相干性阈值，利用GACOS对形变场进行大气校正，再根据余震分布和地质背景确定发震断层的基本形状，最后基于Okada弹性半空间位错模型反演发震断层的滑动分布。改正后的InSAR同震形变场显示，视线向最大下沉为25 cm，最大抬升为10 cm，分别位于震中西北和东南，形变长轴为北西向,主要形变区位于发震断层西部。改正后的InSAR形变场残差均方根较改正前小，最大滑动量为0.9 m，平均滑动角为-0.5°，破裂主要集中在地下1~20 km范围内，矩震级为6.5，与USGS和GCMT等机构的结果一致。研究表明，利用GACOS改正九寨沟InSAR同震形变场对提高形变场精度具有一定的作用，反演断层滑动分布的结果较改正前差别不明显，发震断层的属性与虎牙断裂北段的性质基本一致。结合余震重定位结果可以推断，发震断层为虎牙断裂的北向延伸部分，此次地震事件为巴颜喀拉地块南东向扭转与华南地块碰撞的结果，中下地壳粘性流体的差异分布是导致虎牙断裂倾角变化的主要原因。     

同震变形分析中不同地球结构模型差异性研究

顾及地球结构的非均匀性、分层和球体特性，定量研究了球体分层模型、球体均匀模型、半空间均匀模型分别作用于地震断层（走滑断层，逆冲断层）产生的地表同震变形的差异。计算结果显示，相对实际地球模型，利用半空间均匀位错模型计算的同震变形差异很大，大大超出现有高精度测量技术的精度。如果利用实测资料（例如GPS，VIBL，InSAR等高精度测量数据）进行反演研究，必须根据实际地壳构造情况考虑地壳分层结构的影响。     

2021年漾濞MW6.1地震同震形变提取与断层源模型反演研究

基于D-InSAR技术,收集Sentinel-1卫星升、降轨及精密轨道数据,结合外部DEM,成功提取2021-05-21云南省大理州漾濞县M_W6.1地震的同震形变场。结果显示,升轨卫星视线向最大形变量级约为6.0 cm,降轨卫星视线向最大形变量级约为7.9 cm。以升、降轨InSAR观测结果为约束,对漾濞地震的断层几何参数及滑动分布进行联合反演发现,此次地震的最优发震断层走向为136.6°、倾角为83.1°,断层破裂主要集中在地下2~12 km深度处,最大滑动量约为0.45 m,位于地下7 km深度处,断层在近地表未出现大面积显著滑移现象,表明此次地震未破裂至地表。     

基于升降轨Sentinel-1数据分析2017-11-12伊拉克MW7.3地震震源参数

利用Sentinel-1卫星提供的多个轨道SAR数据对2017-11-12伊拉克M_W 7.3地震进行研究，通过处理不同轨道数据并比较，获取震区可靠的同震形变场。利用降轨同震形变场，采用数据分辨率约束的四叉树方法降采样后作为约束条件，基于弹性半空间模型（Okada）反演断层几何参数及滑动分布。结果表明，本次地震位于扎格罗斯褶皱冲断带，产生的破裂带长度为44.5 km、宽度为16.1 km；发震机制为走滑逆冲断层，断层滑动角为140.68°，倾角为8°，最大滑动量为3.83 m，地震造成视线向最大抬升和最大下沉为50 cm和40 cm；此次地震震源深度为15.5 km，累计释放地震矩0.97×1020 Nm。     

2015年西藏定日MW5.7地震震源参数估计和静态应力触发研究

采用InSAR形变监测资料计算2015年西藏定日M_W5.7地震同震形变场,反演发震断层几何参数和滑动分布。在此基础上,研究尼泊尔M_W7.8主震对定日地震的静态库仑应力触发影响。综合分析地震滑动机制和构造特征,认为定日断层为西倾隐伏断层。反演结果表明,地震破裂相对集中,主要深度在6~9 km,破裂以正断滑动为主。发震断层走向约178°,倾角约 48°,破裂区长约5 km,宽约5 km,最大滑动量约0.2 m,释放的地震矩约3.7×10¹⁷ N·m,对应矩震级M_W5.6。尼泊尔主震同震库仑应力在定日地震震源处约为0.2 bar,造成藏南申扎-定结拉张地堑应变加载。     

荣县地震InSAR同震与时间序列形变初步分析

借助Sentinel-1A卫星升降轨数据，利用合成孔径雷达差分干涉测量(DInSAR)和Okada均匀滑动模型反演获得四川荣县地区三维同震形变场及荣县地震震源参数。结果表明，此次地震造成的断裂现象大体沿NS向，两侧沉降区域中间伴随抬升现象，最大垂直形变量约18 mm；荣县地震发震断层符合逆断层特性，断层长约9.5 km，宽约2.1 km。采用时间序列分析技术获取该地区2018-11～2019-03各时间段的累计形变，结果发现，该区域在震前形变波动较小，地震是造成累计形变较大的主要原因，震后地表变化趋于稳定。     

菲律宾棉兰老岛2019年4次MW>6.0地震断层运动模型及库仑应力传输研究

基于雷达干涉测量技术,利用ALOS-2、Sentinel-1卫星升降轨雷达影像,获得2019-10~12发生在菲律宾棉兰老岛的4次M_W>6.0地震的同震形变场,并以此形变结果为约束,反演得到4次地震的断层运动模型。综合分析发现,此次地震序列由3条断裂的破裂引起,其中2019-10-16和2019-10-31的2次地震为同一发震断裂,2019-10-31地震断层破裂区域位于2019-10-16地震断层破裂的东北延伸段,最大滑动量约为1.1 m,约为2019-10-16地震最大滑动量的2倍。2019-10-29地震由一条独立断层破裂引起,断层最大滑动量约为2.0 m。2019-12-15地震由一条东北向倾斜断层破裂引起,断层最大滑动量约为3.0 m。此外,2019-10-16地震引起2019-10-29地震显著滑动区明显的正向库仑应力传输;而2019-10-29地震显著增加了2019-10-31地震震源区域的库仑应力;前3次地震对2019-12-15地震孕震断层的库仑应力传输总和为负值,说明静态库仑应力传输可能不是此次地震触发的主要诱因。     

2015年塔吉克斯坦MS7.4地震的InSAR观测及其构造意义

利用欧空局提供的Sentinel-1A卫星差分干涉数据为约束，对2015-12-07塔吉克斯坦MS7.4地震的震源机制进行反演。首先运用两轨法对卫星雷达影像进行差分干涉处理，获取覆盖塔吉克斯坦地震震区的同震形变场，然后运用分布式滑动模型反演获取较精细的断层滑动分布。结果显示，分布式滑动模型与观测结果有很高的拟合度。反演结果表明，发震断层以左旋走滑运动为主，此次地震断层的最大滑移量为4.42 m，同震的地震矩为73.71×1018 Nm，矩震级为MW7.19，与震后USGS、GCMT、IPGP等机构给出的震源机制解结果一致，是在具有左旋走滑运动特征的显著构造背景下的一次正常、必然的破裂事件。     

新疆于田MS7.3地震同震与震后形变机制研究

利用InSAR技术获取2008-03-21新疆于田MS7.3地震的同震和震后形变场。同震分布式滑动反演结果表明,同震断层最大滑动量达5.4 m,主要分布在南部断层的0~5 km深度附近,地震以正断错动为主,兼有左旋走滑分量。震后形变结果表明,发震断层北段两侧存在差异性运动,最大累积差异形变在震后782 d达15 cm。进一步分析表明,震后断层余滑可能是震后形变的主要机制。余滑反演结果表明,震后2 a断层余滑量相对较小,滑移区范围明显减小且均位于浅部区域,北部断层能量释放较彻底,南部仍存有少量能量,整体能量基本释放完。     

基于Sentinel-1A数据反演2015年新疆皮山MW6.5地震断层滑动分布

使用Sentinel-1A卫星观测的2015年新疆皮山地震震前震后一对SAR影像,提取该地震的同震形变场。结合Okada位错模型,在考虑渐变入射角和解缠基准偏差的情况下,反演本次地震的断层滑动分布。结果表明,皮山地震是小倾角逆冲型地震,断层破裂开始于地表以下约5 km处,最大滑动量为0.76 m,高滑区集中于深度9～14 km,累计释放地震矩达5.14×10¹⁸,相当于矩阵级M_W6.47。     

基于Sentinel-1A的2021年双湖MW5.7地震同震破裂模型

利用Sentinel-1A SAR数据提取2021年西藏双湖县M_W5.7地震同震形变场及2.5D形变场,反演断层滑动分布模型。计算不同节面解为接收断层产生的库仑应力变化差异确定发震构造,并结合余震分布信息评估未来地震风险性。结果表明,地震震中为34.37°N、87.71°E,震源深度6.51 km,发震断层倾向东、走向33°、倾角50°、平均滑动角-74°,以倾滑为主兼有少量左旋走滑分量,最大滑动量0.26 m。短时间内,震区南部地震风险较小,北部则需要结合更多资料进一步分析。本次地震是在羌塘块体持续向东扩张的背景下,受EW向拉伸作用使得黄水湖正断层发生的一次弥散型变形活动,SN向地堑得到进一步扩张。