2013年芦山Ms 7.0级地震断层参数模型反演

2013年4月20日四川芦山地区发生了Ms 7.0级地震。利用GPS三维同震形变数据获取地表形变场,基于位错模型反演芦山地震的断层几何参数及滑动分布。首先采用多峰值颗粒群算法（multiple peak particle swarm optimization,MPSO）得到断层几何参数,其中断层走向206.47°,倾角44.11°,长度21.94 km,离地表最浅处为7.66 km,最深处为17.84 km。为了反演断层面的精细滑动分布,分析地震所在的龙门山断裂西南段破裂面具有的铲状型特征,将芦山地震破裂面确立为铲状模型,即将断层的倾角预设为上陡下缓,倾角变化范围为21°~50°。结果显示,断层破裂面在不同深度区域出现了两个滑动峰值,其中最大滑动量为0.68 m,深度位于13 km。地震释放的能量为1.47×1019 N·m,对应的矩震级为Mw 6.74,与地震学的研究结果一致。     

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

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

2008年大柴旦M_w6.3级地震的InSAR同震形变观测及断层参数反演

2008年11月10日,青海省大柴旦地区发生了Mw6.3级地震。本文利用EnviSat卫星升降轨SAR数据和差分干涉测量技术提取了同震形变场,基于均匀位错模型反演确定了地震断层参数,然后利用格网迭代搜索法确定了较优断层倾角,同时基于非均匀位错模型获得了精细滑动分布。结果表明,地震使得上盘区域沿降轨、升轨视线向分别产生最大约8.5cm、10cm的抬升;较优断层倾角为47.9°;地震滑动未延伸至地表,主要发生在地下8.2~23.7km深度范围内,最大和平均滑动量为0.5m和0.19m,平均滑动角为104.9°。反演的地震矩为3.74×1018 N·m,矩震级为Mw6.35。     

基于多源SAR数据的2022年门源Ms 6.9地震同震破裂模型反演研究

2022年1月8日青海省海北州门源县发生Ms 6.9地震,震中位于青藏高原东北缘祁连-海原断裂中段,属历史地震空区,基于多源合成孔径雷达（synthetic aperture radar,SAR）遥感数据研究该地震的破裂模式对理解青藏高原东北缘构造变形机制、应变释放过程以及地震危险性评估具有重要意义。首先利用Sentinel-1数据和合成孔径雷达差分干涉测量（differential interferometry synthetic aperture radar,D-InSAR）技术获取了门源地震的同震形变场,视线（line of sight,LOS）向形变场显示此次地震造成了约20 km长的地表破裂,最大形变约0.75 m;然后基于Sentinel-2卫星数据,利用光学影像配准和相关技术获取了本次地震的东西向同震形变场,最大同震位移达2.5 m;最后基于均匀弹性半无限位错模型,以LOS向形变场为约束反演了断层的滑动分布模型。结果显示,门源地震是一次典型的左旋走滑型地震,地震破裂主要集中在0~10 km深度范围,最大滑动量3.25 m,滑动角10.44°,对应深度4.89 km;反演给出的矩震量为1.07×1019 N·m,对应矩震级Mw 6.6。结合野外考察和地质资料,初步判定发震断裂为冷龙岭断裂,并引起托莱山断裂发生同震滑动。同震库仑应力结果显示,冷龙岭断裂东段和托莱山断裂西段应力状态为加载,未来具有发生强震的风险。     

利用方差分量估计的地震同震滑动分布反演

提出了基于方差分量估计的地震同震滑动分布反演方法,该方法不仅可以确定不同数据集的相对权重,而且还能得到平滑因子的值,突破了方差分量估计方法仅用于各类数据定权的实际应用范围。模拟算例验证了此方法的有效性,从解的概率后验密度分布分析来看,方差分量估计方法得到的结果是线性和非线性混合方法反演结果的一个解。用方差分量估计方法反演了Bam地震同震滑动分布,并与线性非线性混合方法反演的结果进行了对比分析。两种方法得到的Bam地震断层滑动分布之间的空间互相关系数为0.999 9,地震同震的最大滑动量相同,都为3.04m;其深度仅有微小区别,前者为4.50km,后者为4.47km;平均滑动量前者为0.714m,后者为0.718m。但前者反演计算用时227s,后者却达到1.5×106 s,约17d,前者仅为后者的1/6 608。结果表明,本文方法相对于线性和非线性混合方法具有计算简单、计算量小和计算效率高的优点。     

多视线向InSAR形变场约束的改则地震滑动分布反演

利用3种不同视线向LOS(Line Of Sight)的ENVISAT ASAR数据进行干涉处理,提取多视线向(Multi-LOS)的同震形变场;结合同震形变场特征与震源机制解,构建了改则地震双断层破裂模型;利用四叉树采样后的多视线向同震形变场进行约束,通过梯度下降法(Steepest Descent Method,SDM)与Crust2.0地壳分层模型反演了改则地震的同震滑动分布特征。结果表明:反演的形变残差得到有效控制,基本介于0±10 cm之间;主震断层的滑动量主要位于断层面2—16 km深部,最大滑动量可达1.34 m,位于断层面6.4 km深处;余震断层滑动量主要位于断层面2—6 km深部,最大滑动量可达0.90 m,位于断层面3.52 km深处;主震断层与余震断层均以正断为主,但主震断层还具有一定的左旋走滑分量,而余震断层的左旋走滑不明显;当剪切模量μ取3.2×1010Pa时,反演获得的主震与余震地震矩M0分别为6.34×1018N·M与1.20×1018N·M,分别相当于矩震级MW6.47与MW5.98。     

利用断层自动剖分技术的2008年青海大柴旦Mw6．3级地震InSAR反演研究

2008-11-10青海大柴旦地区发生了Mw6．3级地震,其发震断层位于青藏高原东北缘的大柴旦一宗务隆山断裂带。利用欧空局Envisat／ASAR卫星雷达影像数据,采用二通差分干涉技术获得了地震的同震地表形变场,基于1D协方差函数估计InSAR同震形变场的中误差为0．52cm,方差一协方差衰减距离为5．9km。在此基础上,采用弹性半空间矩形位错模型进行断层几何参数反演,并利用断层自动剖分技术确定了地震的最佳同震滑动分布模型。结果表明,该地震的震源机制解为走向107．19°,倾角56．57°,以逆冲为主兼具少量右旋走滑分量;滑动分布主要发生在10-20km深度范围内,最大滑动量为0．51m,释放的能量为4．3×10^18Nm。     

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数据对震源机制进行反演、模拟和分析,获取了高精度的震源参数,可以为分析地震危险性和断层破裂参数等提供数据支持。     

日本2011 Tohoku-Oki Mw 9.0级地震的同震形变及其滑动分布反演:GPS和InSAR约束

利用覆盖整个日本东北部的3个条带的Envisat/ASAR降轨数据和6个条带的ALOS/PALSAR升轨数据,通过二通差分干涉处理,首先得到了2011年日本Tohoku-Oki Mw 9.0级地震的初步同震形变干涉结果,然后,利用GPS同震观测值对InSAR结果进行校正和基准统一。在此基础上,采用弹性半空间矩形位错模型,联合GPS和InSAR观测结果对同震滑动分布进行反演,获取了发震断层的断层滑动分布。结果显示,此次地震的滑动分布主要发生在40～50km深度范围内,最大滑移量为50.3m,释放的能量为3.20×1022 N·m(相当于Mw 8.94级)。     

门源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;区域断裂带的深部几何形态以花状构造为特征,整体倾向南西,门源地震发震断裂为花状构造中未出露地表的盲断层。相关成果能够为研究区域地壳运动与变形、活动断裂与地震孕育发生等提供参考。     

2015年尼泊尔地震同震滑动及震后余滑的三角位错模型反演

2015年尼泊尔Mw 7.9级地震发生在印度板块向欧亚板块低角度俯冲的喜马拉雅断裂带上。对该地震的滑动模型进行精化,对于理解喜马拉雅断裂带的孕震模式具有重要意义。采用三角形位错元构建主喜马拉雅断裂“双断坡”几何模型,联合全球定位系统（Global Positioning System,GPS）和合成孔径雷达干涉测量（interferometric synthetic aperture radar,InSAR）资料反演2015年尼泊尔地震同震滑移及震后余滑。结果表明,尼泊尔地震最大同震滑移达到7.8 m,深度为15 km,位于中地壳断坡和浅层断坪的接触部位。不考虑中地壳断坡结构会使反演的最大滑移量偏低。震后余滑主要分布在同震破裂区北侧,释放的地震矩为1.02×1020 N·m,相当于一次Mw 7.3级地震,约占主震释放地震矩的12%。同震库伦应力变化和震间断层闭锁分布均表明,尼泊尔地震破裂区南部宽约60 km的区域仍具有较高的地震危险性。     

2022年泸定Mw 6.6地震InSAR同震形变与滑动分布

2022-09-05,青藏高原东缘的鲜水河断裂上发生了泸定Mw 6.6地震,该地震是鲜水河断裂上40年来发生的最大地震,研究该地震的运动学和同震破裂模式对理解青藏高原东缘构造形变机制和评估鲜水河断裂以及安宁河断裂的地震危险性具有重要意义。利用Sentinel-1和ALOS-2卫星雷达影像,采用合成孔径雷达干涉技术获取了泸定地震的同震形变场,进而基于弹性半空间的位错模型,确定了本次地震发震断层的几何参数和滑动分布。结果表明,泸定地震是一次典型的左旋走滑事件,发震断层西倾,倾角约为72°,走向沿NNW-SSE方向,约为167°;断层破裂主要集中在0~10 km深度,最大滑动发生在约5.8 km深度,约为2.23 m;同震释放的地震矩约为8.74×1018 N·m,相当于矩震级Mw 6.59。通过对震后光学影像解译,发现此次地震诱发的滑坡多集中分布在发震断层西侧,该现象与余震主要集中在断层西侧的结果相一致,可认为是地震上盘效应的体现。     

InSAR数据约束下2016年和2022年青海门源地震震源参数及其滑动分布

中国青海省门源县于2016年和2022年分别发生了Mw 5.9和Mw 6.7地震,相距不足40 km。利用欧洲空间局Sentinel-1A升降轨雷达影像,采用合成孔径雷达干涉测量(interferometric synthetic aperture radar, InSAR)技术分别获取两次地震的同震地表形变场,进而利用弹性半空间的位错模型确定上述事件的震源参数,基于分布式滑动模型反演确定两次地震断层面上的滑动分布,并探讨2016年门源地震对2022年门源地震的发震影响及触发机制。结果表明,2016年门源地震为逆冲型地震,并未破裂到地表,升、降轨同震形变场沿视线向的最大形变量分别为6.7 cm和7.0 cm,断层的最大滑动量为0.53 m,主要集中在地下4~12 km区域滑动。2022年门源地震同震形变场沿NWW-SEE向破裂,降轨影像最大视线向地表形变量为78 cm,断层的最大滑动值达到3.5 m,处于地下4 km左右,断层滑动分布模型揭示此次地震为左旋走滑型地震;结合冷龙岭断裂的运动性质和几何特征,可初步判定发震断层主要为冷龙岭断裂的西段、且极有可能破裂到了其西北端西侧的托莱山断裂。静态库仑应力触发关系显示,2016年门源地震对2022年门源地震的发生有一定的促进作用。     

唐山7.8级地震三维地表形变和断层运动特征

通过系统收集整理1976年唐山7.8级地震前后区域水准和边角网观测资料,在研究发展动态平差程序的基础上,解算得到了唐山地震震前和同震三维形变定量结果与误差分布,并反演得到了断层面同震滑动分布。结果显示,震前阶段丰润-唐山-唐海-昌黎-迁安一带整体隆升,隆升高值区位于唐山附近,震中区的水平形变量不大,唐山菱形地块相对于周围地块逆时针扭转;同震阶段唐山断裂以右旋错动为主,并呈现西升东降的垂向变形特征;反演结果表明多段断层参与了唐山7.8级地震的应变释放,主震发震断层的南段破裂主要发生在深度约6~18 km、长轴约50 km范围,北段破裂主要发生在深度约7~17 km、长轴约30 km范围,滦县地震断层也发生了少量滑动;主震断层最大水平位错为7.82 m、最大垂直位错为2.04 m,等效震级为Mw 7.58。挖掘得到的唐山震前和同震三维形变场定量结果为深入研究该地震的孕育发生机制及地震影响提供了珍贵资料。