利用同震和震后位移数据联合反演2011年日本MW9.0地震同震断层滑动

基于黏弹性球体位错理论, 联合陆地和海底同震GPS数据以及日本本岛330个陆地GPS站点5~10年的震后数据, 反演了日本M_W9.0地震的断层滑动模型, 提升了断层滑动分布在细节上的合理性。 首先, 基于日本本岛330个陆地GPS站点震前2年和震后10年的连续观测数据, 获取了日本M_W9.0地震震后5~10年的年平均位移, 该时段的位移几乎完全由地幔黏弹性松弛效应引起; 接着, 利用黏弹性球体位错理论对震后5~10年的位移进行反复拟合, 确定了日本M_W9.0地震震源及周边地区的地幔黏滞性系数最优解(9.0×10¹⁸ Pa·s)。 然后, 联合同震和震后位移数据, 引入黏弹性位错格林函数, 反演了2011年日本M_W9.0地震的断层滑动分布。 结果表明, 该地震同震破裂的最大值达到了62.72 m, 同震滑动的总地震矩为4.48×10²² Nm, 相应的矩震级为M_W9.03。 由于黏弹性松弛效应引起的震后位移中包含了同震破裂的信息, 基于黏弹性球体地震位错理论, 联合同震和震后位移数据反演断层同震破裂, 有效提高了日本M_W9.0地震断层滑动分布的可靠性。 最后, 本文提出的反演方法为同震观测结果缺乏的大地震震后科考提供了理论支撑: 在大地震发生之后, 即使在同震期间没有足够的观测数据, 也可以在震后通过对震源区的加密观测积累的震后数据, 使用本文提出的反演方法优化同震断层滑动模型。     

2016年台湾美浓MW6.4地震震源参数的InSAR和GPS反演

2016年2月6日台湾西南部高雄市美浓区发生了M_W6.4地震.本文结合ALOS2卫星升降轨、Sentinel-1A升轨SAR数据,采用两轨差分干涉技术获取了该区域的同震形变场,形变结果表明震中西北部以抬升为主,最大视线向形变量约为11.2 cm.基于均匀位错模型和多峰值粒子群（MPSO）算法,利用InSAR和GPS形变数据联合反演了美浓地震的断层几何参数,结果表明震源中心位于22.920°N,120.420°E,深度约12 km,发震断层长度约15 km,走向角307°,倾角16.5°,平均滑动角为51.5°,此次地震是以逆冲倾滑兼左旋走滑的破裂模式.利用格网迭代搜索法得到最优倾角为15.7°,GPS和InSAR最优权比为18：1,最优平滑因子为0.06.基于非均匀位错模型,利用非负最小二乘方法进行线性反演,结果显示最大倾滑和走滑量分别为51.7 cm和55.3 cm,对应矩震级为M_W6.38,略小于GCMT （M_W6.4）的结果.通过与已有文献的比较和对该区域断层构造的分析,发现美浓地震的发震断层为单一断层的解释更为合理,我们推测发震断层是位于左镇、后甲里等断层之间的一条东南-西北走向往东北倾斜的盲断层,并初步推测2010年M_W6.3甲仙地震也同该断层有关.     

2015年新疆皮山MW6.4地震发震断层和滑动分布反演

本文基于Sentinel-1A卫星影像数据提取了2015年皮山M_W6.4地震的同震形变场, 震中北部以隆升为主, 最大抬升量为12.9 cm; 南部以沉降为主, 最大沉降量为5.5 cm。 采用基于单一断层滑动模型的多峰粒子群优化和蒙特卡罗算法, 以LOS向InSAR形变场为约束, 对发震断层的几何模型进行非线性反演。 在此基础上, 联合InSAR和GPS数据, 利用最速下降法反演断层滑动分布。 综合结果表明: 发震断层是顶部埋深约7.4 km的隐伏断裂, 断层面大小为48 km×35 km, 断层走向、 倾角、 断层滑动角分别为111°、 19°、 91°; 断层最大滑动量0.47 m, 位于深度为10.6 km的区域; 累计地震矩3.89×10¹⁸ N·m, 约合矩震级M_W6.33。 最后, 依据主震断层滑移量计算了主震对周围中小断裂的库仑应力扰动变化, 结果显示距离震中最近的泽普断裂受主震影响的库仑应力明显增加; 震后3年内余震集中分布在泽普断裂库仑应力增加区域, 表明皮山地震主震对余震的发生可能具有一定的应力触发作用。     

利用考虑地形起伏的反演方法获取2017年伊拉克MW7.3地震的断层参数及滑动分布

2017年11月13日伊拉克北部地区苏莱曼尼亚省发生了M_W7.3地震, 造成了重大的人员和经济损失。 本文利用升降轨的Sentinel-1和降轨的ALOS-2卫星的SAR数据, 通过差分干涉测量技术获取了该地震的同震形变场, 联合DInSAR和MAI技术, 采用抗差最小二乘法求解该地震的同震三维形变场。 基于改进的考虑地形起伏的均匀位错模型反演确定了发震断层的断层参数, 最后基于非均匀位错模型得到了发震断层的分布式滑动分布模型。 结果显示: ALOS-2卫星降轨轨道观测到的伊拉克地震引起的LOS向地表形变最大为55.8 cm抬升和47.9 cm下沉; Sentinel-1卫星观测到的伊拉克地震引起的LOS向地表形变为: 升轨轨道最大为87.9 cm抬升和17.1 cm下沉; 降轨轨道最大为55.6 cm抬升和38 cm下沉; 相对于前人的研究成果, 本文利用改进的考虑地形起伏的反演方法得到的发震断层几何参数表明发震断层为NNW走向, 倾向角为352°, 同震破裂以逆冲为主, 同时兼有一定的左旋走滑分量。 基于均匀位错模型反演得到的断层滑动分布结果表明, 同震破裂未延伸至地表, 主要滑动量集中在12~18 km, 最大滑动量位于15 km深度, 达到4.3m, 反演得到的矩震级为M_W7.35, 与UGSG、 GCMT等机构给出的结果一致。     

基于升降轨Sentinel-1 SAR数据研究2020年新疆于田MS6.4地震震源机制

2020年6月26日新疆于田县发生M_S6.4地震。 本文先利用Sentinel-1A/B卫星升、 降轨SAR数据, 结合InSAR技术提取此次地震的同震形变场; 并以同震形变场为约束, 利用贝叶斯方法的均匀滑动模型反演发震断层的几何参数; 最后基于梯度下降法(Steepest Decent Method, SDM)来确定发震断层更为精细的滑动分布。 研究结果表明: InSAR形变场的分布主要沿北北东方向, 东西方向跨度约40 km, 南北方向跨度约30 km, 形变特征为西侧沉降与东侧隆升, 西侧最大沉降约200 mm, 东侧最大隆升约70 mm。 于田地震发震断层为一条走向187.68°、 倾角59.78°、 滑动角77.76°的隐伏断层; 于田地震的宏观震源深度为6.29 km, 距震级M_W6.19。 断层滑动分布主要集中在沿走向7~21 km与倾向4~11 km的区域内, 平均滑动量约0.2 m, 在沿倾向7~8 km处的最大滑动量约0.97 m, 同震位错主要表现为正断运动。     

2016年新疆阿克陶地震Sentinel-1 InSAR同震形变特征

2016年11月25日新疆克孜勒苏州阿克陶县发生M_W6.6地震。 本文利用合成孔径雷达差分干涉测量技术, 对Sentinel-1卫星获取的升、 降轨雷达数据进行了处理, 提取了该次地震的同震形变场, 并结合形变场特征与震源机制解, 采用梯度下降法反演发震断层的滑动分布。 结果表明, 升、 降轨LOS向同震形变场在发震断层两侧具有明显不同的形变特征, 主要形变区域分布在断层两侧, 升轨LOS向形变量可达-8.2 cm与11.2 cm, 降轨LOS向形变量可达-21.4 cm与13.1 cm; 反演的升、 降轨干涉形变场与InSAR测量值之间的残差得到有效控制, 大部分的残差介于±5 cm之间; 断层滑动分布主要集中于沿断层面深约2~18 km处, 最大滑动量位于沿断层面深约7 km处可达0.96 m; 平均滑动角约182.29°, 最大滑动处的滑动角约197.13°, 两个滑动分布中心的滑动角均接近180°, 表明阿克陶地震为一典型的右旋走滑破裂性事件; 当剪切模量取32 Gpa时, 反演的发震断层地震矩M₀可达9.75×10¹⁸, 相当于矩震级M_W6.60, 与地震波形反演结果一致。     

2008年于田MW7.2地震滑动分布反演及三维同震形变场解算

基于ENVISAT ASAR升降轨数据, 利用InSAR获取2008年于田M_W7.2地震同震形变场; 采用SDM反演本次地震断层滑动分布; 使用PSCMP正演获取于田M_W7.2地震南北向模拟形变量, 并结合升降轨同震形变场, 解算三维同震形变场。 同震形变场分析表明, 2008年于田M_W7.2地震以正断为主, 且带有走滑运动特征, 破裂带走向为NNE向。 同震滑动分布反演结果显示, 断层沿走向被分为4段F1、 F2、 F3、 F4, 其滑动分布集中在0~14 km区间, 以F2、 F3段为主, 最大滑动量约5.31 m, 位于F2段深部2.76 km处; 沿破裂带走向, 左旋走滑位移与垂直位移比值有增大的趋势; 反演获得的地震矩M₀=5.58×10¹⁹N·m, 相当于矩震级M_W7.1。 三维同震形变场解算结果显示, 断层上盘整体表现为沉降, 断层下盘整体表现为隆升, 且沉降量明显大于隆升量, 表明地震以正断破裂为主; 除靠近断裂带中上部表现为向东南运动外, 上盘整体上表现为向西南运动; 断层下盘则整体表现为向东北运动, 证明破裂兼有左旋走滑运动。 滑动分布反演、 正演与三维同震形变场解算结果皆表明, 于田M_W7.2地震破裂以正断为主, 且带有一定的左旋走滑。     

基于InSAR和远场地震波联合反演2008年MW6.3大柴旦地震震源破裂过程

利用InSAR同震形变升、降轨数据和远场地震波数据,基于均方根最小与标量地震矩最小双重约束下的模拟退火方法,联合反演2008年11月10日大柴旦M_W6.3地震震源破裂过程.结果表明,2008年大柴旦地震震源破裂过程整体表现为沿倾向方向从深部破裂起始点处向上往地表传播,且破裂未到达地表;在前7 s内,滑动沿西北和东南两个方向传播,7 s后主要沿东南方向传播,破裂过程时间持续约为11 s,同震滑动分布主要集中在地下10~20 km范围内,最大滑动量达-0.71 m;反演结果揭示本次地震为西南倾高角度兼具微量走滑分量的逆冲破裂事件,反演矩张量为3.96×10¹⁸N·m,矩震级约M_W6.37.通过大柴旦地震发震断层和破裂机制综合分析,初步判断发震断层为西南倾向的大柴旦—宗务隆山断裂.     

汶川地震同震形变场的GPS和InSAR邻轨平滑校正与断层滑移精化反演

为克服InSAR观测汶川地震同震形变场的邻轨不连续问题,提出联合GPS观测值与邻轨平滑约束的同震位移校正方法,采用GPS观测形变去除PALSAR轨道误差引入的残留平地相位,基于形变平滑条件校正邻轨干涉相位的不连续性.ALOS/PALSAR干涉处理结果表明,校正后同震形变场的准确度与平滑性得以显著提高,InSAR高相干点残差达3.6 cm,校正后精度提高约60%,低相干点精度提高约40%,校正后形变场的邻接平滑因子标准差减小约33%,验证了轨道误差校正与邻轨平滑约束方法的准确性与可靠性.进一步基于弹性半空间位错模型的断层滑移反演结果表明,断层滑移主要分布于映秀、北川和青川地区,集中于地壳深度0~16 km范围,最大滑动量（位于北川县城）约为9.0 m,GPS反演模型残差为5.5 cm,InSAR反演模型残差达9.2 cm,InSAR反演精度约有30%的显著提高,由模型反演计算得到的地震矩为8.0469×10²⁰ N·m.     

基于InSAR形变观测反演2015年皮山MS6.5地震滑动分布

2015年7月3日在新疆皮山县发生了M_S6.5地震, 该地震使当地遭受了巨大的经济损失。 本文利用欧空局提供的Sentinel-1A卫星差分干涉数据对该地震的震源机制情况进行了反演研究, 首先运用两轨法对卫星雷达影像进行差分干涉处理, 获取了覆盖皮山地震震区的同震形变场, 然后利用弹性半空间的均匀滑动模型反演获取了发震断层的几何参数, 并对原始观测数据进行降采样处理, 在此基础上运用分布式滑动模型反演获取了更为精细的断层滑动分布, 结果显示分布式滑动模型与观测结果有很高的拟合度。 反演结果表明发震断层是一个以逆冲为主兼有极少量左旋走滑的盲断层, 此次地震断层面的同震活动分布主要集中在7~15 km深度范围内, 同震的地震矩为6.28×10¹⁸N·m, 矩震级为M_W6.46, 与前人的研究结果非常一致。     

基于InSAR和GPS观测数据的尼泊尔地震发震断层特征参数联合反演研究

利用日本ALOS-2和欧空局Sentinel-1A卫星获得的尼泊尔地震同震形变场,结合GPS同震位移数据,联合反演了断层滑动分布特征和空间展布.结果表明:尼泊尔地震的同震形变场主要集中在150km×100km的范围内,且分为南北两个相邻的形变中心,南形变中心的视线向抬升量约为1.2m,北形变中心的视线向沉降量约为0.8m,均位于发震断层上盘.位于形变抬升区的KKN4和NAST两个GPS站,抬升量和南向运动量均达到了m级,而远离震区的其他GPS台水平和垂直观测量均在1cm以内.联合反演得到的断层位错分布主要集中在沿走向150km,沿倾向70km的范围内,最大滑动量为5.59m,平均滑动量为0.94m.断层面倾角在浅部约为7°,随着深度增加,倾角逐渐变大,到垂直深度20km时倾角接近12°;5月12日MW7.2级余震位于主震破裂区的"凹"型滑动缺损区域;主震破裂区的上边界与MBT空间位置十分吻合,主震破裂区主要集中的MBT以北50~60km处,垂直深度为8~9km,倾角为9°,继续向北时主震破裂面以10°~12°的倾角向深延伸,在18~20km可能与MHT交汇.因此,初步判定MBT为此次地震的发震断层.     

2017年8月8日四川九寨沟7.0级地震InSAR同震形变场及断层滑动分布反演

基于InSAR技术,利用欧空局升降轨Sentinel-1A/IW宽幅数据,获取了2017年8月8日四川九寨沟7.0级地震InSAR同震形变场,并以升降轨InSAR观测结果为约束,反演了断层滑动分布,基于三种不同接收断层计算了同震库仑应力变化.结果表明,同震形变场发生在塔藏断裂、岷江断裂和虎牙断裂交汇的三角地带,升降轨干涉位移均显示本次地震的形变场影响范围约为50 km×50 km,形变场长轴方向为NW向,升降轨观测的形变量相反,反映断层运动性质以走滑运动为主,升降轨数据观测得到的最大LOS （Line of Sight,视线向）形变量分别为~22 cm和~14 cm.非对称形变场反映出断层两侧的运动差异.反演结果显示,最大滑动量约为1 m,平均滑动角为-9°,矩震级为M_W6.5,地震破裂主要集中在地下1~15 km深度范围内,但整体而言本次地震破裂较为充分,基本将该区域1973年及1976年4次 > M_W6.0地震的破裂空区完全破裂.考虑到塔藏断裂和虎牙断裂的运动性质,可初步判定发震断层为虎牙断裂北侧延伸分支.基于三种不同接收断层模型的同震库仑应力变化计算结果反映出该区域以应力释放为主,进一步触发较大走滑型余震的可能性不大.     

基于InSAR分析2015年尼泊尔MW7.8地震形变场特征及断层滑动分布反演

2015年4月25日尼泊尔发生了M_W7.8地震, 本文基于震前、 震后两景Sentinel-1A雷达影像, 采用D-InSAR两轨差分干涉法提取了此次地震的同震形变场。 结果显示, 同震形变场位于喜马拉雅造山带—主边界逆冲断裂(MBT)和主前锋逆冲断裂(MFT)附近, 形变场整体表现为自西北向往东南方向延伸近150 km的纺锤形包络状, 以大面积隆起抬升形变为主, 视线向最大隆升形变达1.18 m, 抬升区北侧存在一小沉陷区, 以InSAR观测值定位同震最大形变中心。 基于均匀介质弹性半空间模型(Okada模型)与InSAR观测数据反演了断层滑动分布。 反演结果表明该地震属于典型逆冲型地震, 发震断层为主喜马拉雅逆冲断裂(MHT), 同震破裂从主喜马拉雅逆冲断裂(MHT)向上沿着主前锋逆冲断裂(MFT)传递。 基于InSAR同震形变场局部形变细节, 结合震区地质背景、 断裂分布及断层运动特征, 获得了同震破裂拟出露地表迹线。     

联合强震记录和InSAR/GPS结果的四川九寨沟7.0级地震震源滑动分布反演及其地震学应用

2017年8月8日我国四川九寨沟发生里氏7.0级地震.本研究利用基线校正方法获得距震中100km范围内9个强震台站同震位移,基于Sentinel-1卫星干涉SAR影像对获取了InSAR同震形变场.结合GPS形变数据,本研究进行了震源滑动模型联合反演,结果显示此次地震整体以走滑运动为主,释放地震矩约为7.60×1018 N·m(～MW6.52).通过对比模拟形变场和观测值显示,联合反演结果优于单独基于InSAR形变场的反演结果.静态应力变化计算结果显示断层平均静态应力降为1.07MPa.反演滑动模型沿走向和倾角方向拐角波数值分别为0.99×10-4和1.10×10-4.同震静态库仑应力变化计算结果显示共有83.6%的余震位于库仑应力增加的区域,被主震所触发的余震占总数的77.9%,主震对后续余震具有显著触发作用.强地面运动模拟结果显示模拟结果在烈度分布范围和等级方面与调查烈度符合度很高,模拟结果能够很好地反映断层破裂的方向性效应等特征.本研究计算结果显示九寨沟地震无论是平均静态应力降还是拐角波数均低于同类型地震的平均水平,这可能是造成本次地震强地震动水平相对不高的原因.     

2014年新疆于田MS7.3地震同震位移及位错反演研究

利用于田震中300 km范围内的1个GPS连续站和12个GPS流动站数据,解算得到了2014年新疆于田M_S7.3地震地表同震位移,并反演了发震断层滑动分布,探讨此次地震对周边断裂的影响.地表同震位移结果显示,GPS观测到的同震位移范围在平行发震断裂带的北东-南西向约210 km,垂直发震断裂带的北西-南东方向约为120 km,同震位移量大于10 mm的测站位于震中距约120 km以内;同震位移特征整体表现为北东-南西方向的左旋走滑和北西-南东方向的拉张特征,其中在北东-南西方向,I069测站位移最大,约为32.1 mm,在北西-南东方向,XJYT测站位移最大,约为28.1 mm;位错反演结果表明,最大滑动位于北纬36.05°,东经82.60°,位于深部约16.6 km,最大错动量为2.75 m,反演震级为M_W7.0,同震错动呈椭圆形分布,以左旋走滑为主并具有正倾滑分量,两者最大比值约为2.5：1,同震错动延伸至地表,并向北东方向延伸,总破裂长度约50 km,地表最大错动约1.0 m;同震水平位移场模拟结果显示贡嘎错断裂、康西瓦断裂和普鲁断裂等不同位置主应变特征具有差异性,这种差异特征是否影响断裂带以及周围区域的应力构造特征,值得关注.