联合GNSS和InSAR观测位移反演2008年汶川大地震断层位错模型参数

利用现代空间大地测量技术,尤其是卫星合成孔径雷达干涉测量,能够获取高精度、高空间分辨率的同震和孕震形变,为地震断层形变和破裂机制研究提供了前所未有的机遇。本文介绍了利用大地测量观测数据反演地震断层位错模型参数的贝叶斯反演方法。联合运用2008汶川大地震前后GNSS和InSAR技术观测获得的同震位移,反演了地震断层的几何参数和滑动位错分布。研究结果表明,汶川地震的断层滑动主要集中在倾角较陡的浅部,同时包含逆冲和右旋走滑,其中最大逆冲6.1m,最大右旋6.5m。根据断层滑动分布正演计算得到的上盘同震位移明显小于下盘,预示该断层两侧孕震形变可能存在较大的不对称性。     

基于敏感性迭代拟合算法反演汶川Mw 8.0地震InSAR同震形变场

基于单一断层模型,运用敏感性迭代拟合算法反演了汶川地震InSAR同震形变场,获得了断层滑动分布及部分震源参数.结果表明,倾角线性变换的单断层模型模拟的同震形变场与InSAR形变场吻合较好,且残差较小,平均残差为0.11 m;反演的滑动分布主要集中于地下0～20km,其中,汶川地区在震源深度附近有较大的滑动量,10～15 km深度最高可达7 m,地表滑动量却相对较小,平均值仅为2.5 m,平均滑动角均值约为121°;北川地区最大滑动量可达到10 m,平均滑动角均值约为109°;而青川地区10～15 km深度滑动量最高也达8 m,平均滑动角均值约为135°;滑动矢量图显示,沿SW-NE走向断层面的滑动方向以强烈的逆冲为主,兼有一定右旋走滑分量.反演矩张量为1.0×1021 N·m,矩震级达Mw8.0.     

基于InSAR同震形变场反演汶川M_w7.9地震断层滑动分布

通过综合分析2008年5月12日汶川地震野外地震地质考察的地表破裂带空间分布及分段资料,结合InSAR干涉形变场资料,构建了五段断层儿何结构模型,该模型与野外地震地质考察结果在多数分段上基本一致;基于此五段断层模型,运用敏感性迭代拟合算法反演了汶川地震InSAR同震形变场,获得了断层滑动分布及部分震源参数.结果表明,基于余震精定他获得的地震断层倾角模型模拟的同震形变场与InSAR形变场吻合较好,且残差较小;反演的滑动分布主要集中于地下0～20 km,最大滑动量分别位于北川及青川等地区,最大可达到10 m;沿SW-NE走向,断层面的滑动方向主要以右旋兼逆冲形式为主,在汶川及都江堰地区以强烈的逆冲为主兼有一定右旋走滑分量,在北川及映秀地区以逆冲兼右旋运动为主,在平武及青川等地区则逐渐过渡为以右旋运动为主兼有一定的逆冲分量,其中汶川地区的平均滑动角为97°,北川地区的平均滑动角为119°,青川地区平均滑动角为138°.反演矩张量为7.7×10~(20)N·m,矩震级达M_w7.9.     

Source parameters of the 2009 L'Aquila earthquake,Italy from Envisat and ALOS satellite SAR images

2009年4月6日意大利L’Aquila地区发生了Mw6.3级地震,该地震造成了300余人的人员死亡.本文联合不同波长、不同入射倾角的升降轨Envisat和ALOS卫星的差分干涉数据对该地震进行震源机制解的反演研究.研究首先对卫星雷达影像进行二通差分干涉处理,获取了覆盖L' Aquila地震震区的完整InSAR同震形变场,然后结合四叉树和均匀采样方法对原始观测数据进行降采样.在此基础上,联合GPS形变观测数据,利用弹性半空间矩形和三角位错模型,以及断层自动剖分技术对断层面进行最优离散剖分,反演获取了发震断层的精确几何参数和最优断层滑动分布,结果显示分布式三角位错滑动模型能够很好地解释观测到的地表形变场.反演结果表明发震断层是一个以正倾滑为主兼有少量右旋走滑的盲断层;基于观测数据最优确定的断层单元的最短边长为0.4 km,最长边长为6.3 km;此次地震的滑动分布主要发生在5～14 km深度的范围内,最大滑移量为1.07 m,释放的能量为3.43×1018N·m(Mw6.32),与地震学的研究结果非常一致.     

基于InSAR同震形变场反演汶川Mw7.9地震断层滑动分布

通过综合分析2008年5月12日汶川地震野外地震地质考察的地表破裂带空间分布及分段资料,结合InSAR干涉形变场资料,构建了五段断层几何结构模型,该模型与野外地震地质考察结果在多数分段上基本一致;基于此五段断层模型,运用敏感性迭代拟合算法反演了汶川地震InSAR同震形变场,获得了断层滑动分布及部分震源参数.结果表明,基于余震精定位获得的地震断层倾角模型模拟的同震形变场与InSAR形变场吻合较好,且残差较小;反演的滑动分布主要集中于地下0～20 km,最大滑动量分别位于北川及青川等地区,最大可达到10 m;沿SW-NE走向,断层面的滑动方向主要以右旋兼逆冲形式为主,在汶川及都江堰地区以强烈的逆冲为主兼有一定右旋走滑分量,在北川及映秀地区以逆冲兼右旋运动为主,在平武及青川等地区则逐渐过渡为以右旋运动为主兼有一定的逆冲分量,其中汶川地区的平均滑动角为97°,北川地区的平均滑动角为119°,青川地区平均滑动角为138°.反演矩张量为7.7×10²⁰ N·m,矩震级达M_w7.9.     

2016年4月16日熊本地震滑动分布反演与分析

综合利用强震数据、GPS数据和InSAR数据基于双断层模型反演熊本地震滑动分布,通过选择合理的介质模型和平滑因子,分别对数据进行单独反演和联合反演。从结果分析可以看出:三种数据联合反演的结果最优,最终滑动模型为:断层1走向为236°,倾角65°,滑动角-150.6°,最大滑动量为6m;断层2走向为206°,倾角72°,滑动角-155°,最大滑动量为4m。基于K-net和Kik-net获取永久位移快速反演得到的滑动分布结果与基于GPS数据,Sentinel-1A InSAR数据反演甚至联合反演得到滑动分布结果比较一致,表明大震后利用高密度强震动台网后快速获取滑动分布用于震后应急响应和灾害评估是切实可行的,同时认为此次地震发震断层为右旋走滑的断层系统。     

2016年新疆阿克陶MS6.7地震同震形变特征和断层滑动分布反演

本文针对2016年11月25日在新疆阿克陶发生的地震, 用差分干涉测量技术(D-InSAR)对3种不同观测模式的升、 降轨数据进行处理, 提取了多视线向的同震形变场; 根据不同模式的LOS向形变量, 构建形变分解模型, 将其分解为垂直向形变量和沿断层走向形变量; 结合同震形变场特征与震源机制解, 采用单断层模型, 利用梯度下降法(SDM), 以Multi-LOS向形变进行约束, 反演了阿克陶地震的同震滑动分布特征。 研究结果表明, 升、 降轨LOS向同震形变场在发震断层两侧具有不同的形变特征, 发震断层走向近EW向; LOS向形变量分解表明, 此次地震破裂以右旋走滑为主; 滑动分布反演的形变残差介于0~5 cm之间, 发震断层的滑动量主要位于2~16 km深部, 最大滑动量可达1.02 m, 位于断层面深部5.83 km处, 最大滑动量处的滑动角为185.24°; 平均滑动角为181.32°, 平均滑动量为0.12 m; 滑动分布反演也证明该地震为右旋走滑破裂事件, 与LOS向形变分解结果一致; 当剪切模量μ=3.2 GPa时, 反演得到的地震矩震级约M_W6.6。     

利用InSAR获取2020年新疆于田Mw6.3地震同震形变场与断层滑动分布反演

本文基于InSAR技术,利用欧空局Sentinel-1A/B升降轨SAR数据,提取了2020年6月26日新疆于田Mw6.3地震的同震形变场.利用升、降轨同震形变场约束,分别采用MPSO算法和Bayesian方法反演此次地震发震断层均匀滑动的几何参数,并进行对比.然后采用SDM方法获得发震断层非均匀滑动分布,并分析了同震...     

利用 Yabuki＆Matsu'ura 反演方法计算2011年日本东北地区太平洋海域Mw9.0级地震同震滑动分布

Yabuki & Matsu'ura反演方法是利用ABIC最佳模型参数选取方法和平滑的滑动分布作为约束条件,由形变观测数据计算发震断层滑动分布.本文基于日本列岛同震GPS观测数据和发震断层曲面构造模型,利用Yabuki&Matsu'ura反演方法计算2011年日本东北地区太平洋海域Mw9.0级地震的发震断层同震滑动分布.反演结果表明,断层面上的最大滑动量为35 m,较大滑动分布在浅于30 km的震源中心上部,最大破裂集中在20 km深度的地方,其地震矩约为3.63×1022N·m,对应的矩震级为Mw9.0.模拟结果显示Yabuki&Matsu'ura反演方法更适用于倾角低于40°的断层模型反演.最后,本文基于上述方法获得的发震断层滑动模型,利用地球体位错理论正演计算该地震在中国及其邻区产生的远场形变,正演计算结果基本可以解释由中国GPS陆态网络观测到的同震形变.     

Sentinel-1A及ALOS2数据约束下的阿克陶地震震源参数反演及发震构造研究

2016年11月25日,在新疆阿克陶县发生了M_W6.6地震.本文利用Sentinel-1A宽幅数据和ALOS2精细数据获得了同震形变场,干涉形变场沿木吉断裂展布,显示本次地震破裂长度可达70 km,在地表形成两个形变中心,且震中东部形变场条纹密集而西部稀疏、影响范围广,跨断层剖面显示视线向最大形变量可达12 cm.利用一种结合先验知识的多视角最小二乘迭代分解法求解出地表的三维形变场,结果显示震中东部形变中心垂向向下运动最大可达20 cm,木吉断裂南侧西向运动可达10 cm,断层表现为右旋走滑兼具正断作用.采用前向模拟的方法确定二段式分段断层模型能够较好地恢复观测形变场,进而以InSAR观测数据为约束,基于弹性半空间形变模型采用两段非均匀断层滑动模型来反演断层面上的精细滑动值.反演结果显示本次地震可能为两次地震事件,发生在断层西段的第一次事件以右旋走滑为主,走向103°,倾角76°,滑动角-167°,震源深度10.1 km,累计地震矩为7.2×10¹⁸N·M;东段的第二次事件为右旋走滑兼具正断作用,走向109°,倾角略缓约55°,滑动角-160°,震源深度5.3 km,累计地震矩7.76×10¹⁸N·M.本次地震是一次发生在公格尔拉张系的拉张环境下的构造地震.