2021年双湖MS5.8地震InSAR同震形变场及断层滑动分布反演

基于Sentinel-1A数据,利用合成孔径雷达差分干涉(DInSAR)技术提取视线向(LOS)同震形变场,利用基于序列蒙特卡罗采样的贝叶斯方法反演发震断层几何参数与断层滑动分布。结果表明,发震断层走向约47.810°±2.218°,倾角约36.664°±2.499°;震中位置为87.715°E、34.365°N,震源深度为5.673 km。同震滑动分布主要集中在3.9~7.5 km深度,最大滑动量约0.42 m,平均滑动角约-66.598°±3.258°。本次地震以正断为主,兼具少量左旋走滑性质,反演得到的矩震级为M_W5.5,略小于USGS提供的震级。     

2023年塔吉克斯坦MW7.2地震InSAR同震形变场与滑动分布

利用Sentinel-1A卫星SAR影像数据,对2023年塔吉克斯坦M_W7.2地震开展同震形变提取,基于弹性位错模型进行断层反演,并以本文反演得到的右旋节面解为接收面,计算不同深度的静态库仑应力。同震形变结果显示,升轨LOS向最大形变量达15 cm,降轨LOS向最大形变量达16 cm。断层反演结果表明,此次地震最优发震断层走向为131.1°、倾角为85.7°,同震主滑移区分布在深度10～30 km范围内,以右旋走滑为主,最大滑移位置位于地下约20 km深度处,滑移量为3.49 m,未破裂至地表,矩震级为M_W7.16。库仑应力结果显示,该区域库仑应力符合帕米尔高原已有的应力场及地质学研究结果,随着深度增加,其影响范围以发震断层为中心向外扩张,且自5 km深度往下,应力加载区逐渐侵蚀应力卸载区,并开始以加载区为主,在约10 km深度处开始发生余震活动,与本次发震断层相邻的2条断层未来短时间内地震风险性较小。     

改则Ms6.9地震InSAR同震形变场模拟

以InSAR同震形变场为约束,采用各向同性弹性半空间位错模型,利用试错法对2008年西藏改则县扎西错\%M\%s6.9地震断层几何形状进行格网搜索寻优模拟。观测值离散化得到1 010个数据总集,并重采样80次得到80个数据子集,每个子集之间存在微小的噪声差异,采用Bootstrap方法进行80次搜索寻优。首先用单断裂矩形断层模型进行模拟,最大残差为30 cm。然后用双断裂矩形断层模型进行模拟,最大残差为15 cm。双断裂模型反映了观测形变场的基本特征,并与地震资料符合。根据80次搜索寻优得到主断裂倾滑量95%置信区间为1.6~2.4 m。模拟表明,地表变形是主、余震两次破裂的结果。实际观测和数值模拟结果表明,InSAR可以作为监测青藏高原大陆动力学的有效技术。     

2020年土耳其MW5.7地震InSAR同震形变特征与破裂滑动分布

为准确认识2020-06-14土耳其MW5.7地震的发震位置、构造特点以及地震危险性,利用D-InSAR技术对Sentinel-1A数据进行处理,基于GACOS进行大气校正,获得视线向（line of sight,LOS）同震形变场。降轨LOS向同震形变场显示,断层北侧抬升,最大抬升形变量约8.87 cm;南侧沉降,最大沉降量约-7.75 cm。以LOS向同震形变为约束,先采用贝叶斯自举优化法反演发震断层几何参数,然后使用有限断层方法反演地震破裂滑动分布。结果显示,断层走向约257.48°±0.65°、倾角约79.69°±0.98°、滑动角约154.2°±3.8°,震中位置为（40.754°E,39.389°N）,破裂区域长度约8 km、宽度约6 km,破裂的最浅埋深约0.8 km、最大埋深约8.9 km,最大滑动量约0.57 m,对应深度约4.278 km。地震释放地震矩约4.54×1017 Nm,对应矩震级MW5.7,与土耳其灾害和紧急情况管理局公布的结果一致。此次土耳其地震受近东西向的潜伏断层控制,以右旋走滑为主兼具少量的逆冲特性。地震造成部分地区的库伦破裂应力增量超过0.1 bar,这些区域未来的地震危险性值得关注。     

2022年门源MS6.9地震同震应变观测

以2022年门源M_S6.9地震为例,将17个台站的四分量钻孔应变仪1 sps采样率的同震应变观测值与利用QSSP软件正演计算得到的同震应变模拟值进行拟合,分析二者的相关性。结果表明,17个台站NS向和EW向的P波波形拟合较好,相关系数平均值分别为0.835和0.842,达到高度相关,表明钻孔应变观测数据可较好地反映地震破裂信息,在理论上进一步验证钻孔应变观测数据可用于地震破裂过程的研究。     

利用InSAR数据约束反演2017年Mw6.3精河地震同震破裂模型

利用Sentinel-1降轨卫星数据,研究2017-08-09新疆精河M_W6.3地震的同震形变与断层滑动破裂。考虑轨道和大气误差的影响,基于模型改正后的精细同震干涉图显示,视线向最大地表位移约为7 cm,未发现有地表破裂。此外,基于弹性半空间位错模型,采用两步法策略分别反演地震的均匀断层几何参数和同震滑动分布,结果显示,精河地震发震断层以逆冲运动为主,主破裂区位于8~17 km深度,最大滑移量为0.6 m,表明该地震是一次由逆冲型隐伏断层产生的单次破裂事件。InSAR结果给出的地震释放矩能量为2.91×10¹⁸ Nm,相当于矩震级M_W6.25,与USGS、CENC等机构给出的震级一致。     

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向断裂带上,该断裂可能属于维西-乔后-巍山断裂和红河断裂北段的次级或分支断裂。     

九寨沟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同震形变场对提高形变场精度具有一定的作用,反演断层滑动分布的结果较改正前差别不明显,发震断层的属性与虎牙断裂北段的性质基本一致。结合余震重定位结果可以推断,发震断层为虎牙断裂的北向延伸部分,此次地震事件为巴颜喀拉地块南东向扭转与华南地块碰撞的结果,中下地壳粘性流体的差异分布是导致虎牙断裂倾角变化的主要原因。     

利用球谐分析方法研究2017年九寨沟MS7.0地震前重力场变化

推导地表重力变化数据球谐展开计算公式,并根据该方法分析2017年九寨沟M_S7.0地震前震源周边区域重力场动态演化特征。结果表明：1） 2015-04～2017-04九寨沟地震震源周边区域重力场动态变化明显。2015-10～2016-04、2016-04～10和2016-10～2017-04这3个时段内,研究区内出现重力正负变化过渡带,九寨沟地震震源位于这3个时段过渡带的交会区域内部。震源附近2016-10～2017-04正负重力变化过渡带的空间走向与2015-10～2016-04和2016-04～10时段结果相比发生了约90° 旋转,九寨沟地震发生时间在过渡带空间走向旋转之后。流动重力观测数据较好地反映了九寨沟地震的前兆现象。2）采用球谐方法分析重力变化数据可以有效突出研究范围内大尺度重力场动态演化规律,抑制局部细节特征,有利于识别主要重力场变化特征。     

青海玉树地震的InSAR数据同震形变场模拟与参数反演分析

2010年4月我国青海省玉树藏族自治州的玉树县发生了Ms7.1级地震,这次地震(震中33.2°N 96.6°E,Ms=7.1)是由甘孜-玉树断层的强烈活动引起。甘孜-玉树断裂带大部分分布于青海-西藏高原,北西走向。该断裂带和鲜水河断裂带共同组成了巴颜喀拉地块的南边界,二者在甘孜附近成左阶错列分布。在区域构造背景分析和卫...     

2008年汶川MS8.0地震与2017年九寨沟MS7.0地震成因关系探讨

利用2008年汶川地震破裂模型及2017年九寨沟地震接收断层参数和断层面有效摩擦系数等不同模型参数,计算2008年汶川地震造成的库仑应力变化。结果显示,在考虑不同模型参数条件下,2017年九寨沟地震震源处的库仑应力增量（同震及震后）约为3~7 kPa,尚未达到地震触发阈值10 kPa。综合分析认为,2008年汶川地震对2017年九寨沟地震有一定的触发作用,但作用十分有限,九寨沟地震的发生主要受控于区域构造活动及地壳流变结构。分析不同参数对计算结果的影响发现,不同破裂模型的选取对2018年汶川地震计算结果的影响最大,这可能是导致前人研究结论产生差异的原因之一。     

九寨沟MS7.0地震强地震动模拟及参数敏感性分析

对2017年九寨沟M_S7.0地震进行强地面运动模拟,应用复合震源模型建立特定圆形子震分布的断层破裂运动学模型,对地震矩、应力降、破裂速度、子断层尺寸和震源深度等参数进行敏感性分析,得到一组合理的震源参数,然后基于复合震源模型计算得到该地区地震动,最后与观测的峰值加速度(PGA)、反应谱和PGA的空间分布特征进行对比。结果表明,该模型整体上可以再现各台站地震动强度特征,验证了复合震源模型计算强地面运动的有效性,同时敏感性规律也为复合震源模型计算缺乏强震记录区域的地震动参数的变化区间提供了依据。     

2017年四川九寨沟MS7.0地震震前地壳密度变化特征

基于紧凑约束最小二乘优化反演算法,利用2014~2017年地表重力变化观测成果,反演得到2017年九寨沟MS7.0地震前震中及周边地区的地壳三维密度变化。结果表明：1）在震中及周边地区出现大范围的重力变化,并沿张掖-西宁-九寨沟一带形成最大值达9 μGal的重力变化梯度带;2） 地壳密度变化量随深度增加而增大,其中在45 km附近达到最大,为0.2 g/cm3;3） 震中及邻区范围内形成密度变化的正负四象限分布（2014-04~2017-04）,其表现的物质运移方向呈左旋走滑式分布,与震源机制解基本一致。据此推断,在震前产生了与地震孕育相关的深部构造运动,而强震易发生在物质运移汇聚地方。     

2017-08-08四川九寨沟MS7.0地震对周边区域应力影响

基于2017年九寨沟M_S7.0地震破裂模型及区域均匀分层粘弹性地壳模型,计算此次地震引起的周边地区同震及震后库仑应力变化,讨论地震对周边若干活动断层的影响。结果表明,此次地震增加了虎牙断裂北段发震破裂区两端的库仑应力,最大同震库仑应力增量为148 kPa,在很大程度上超过地震应力触发阈值10 kPa,震后50 a内破裂区两端的库仑应力有所增强,但增量较小,仅在3~5 kPa左右;虎牙断裂南段同震库仑应力增量较小,在3 kPa左右,震后50 a内应力基本无变化;塔藏断裂带西北段的同震应力增量在5～33 kPa之间,在大部分区域超出应力触发阈值,震后应力有所加强,但是增量较小(约为3 kPa)。库仑应力加载有可能增加以上活动断层的地震危险性,需加强监测。地震对雪山断裂带西段(以虎牙断裂与雪山断裂交叉点为界)的影响较小,其同震库仑应力最大增量仅在6.15 kPa左右,未超过地震应力触发阈值。同时,地震还减小了雪山断裂带东段、白龙江断裂舟曲段、岷江断裂北部以及龙日坝断裂北部地区的库仑应力,可能短时间内减小了这些地区的地震危险性;而对距离震中较远(150 km以外)的地区(如青川断裂)影响十分微弱,其同震与震后应力变化可忽略不计。     

2016年杂多MS6.2地震的InSAR形变及断层滑动模型

利用欧空局Sentinel-1A卫星的升轨观测,基于InSAR技术获取2016-10-17青海杂多地区MS6.2地震的同震形变场,并分别采用最速降法和三角元法反演地震的位错分布模型。结果表明,本次地震发生在杂多上拉秀断裂带,走向北东东,地震破裂以正断倾滑为主,兼具少许左旋走滑运动分量,最大滑移量0.17 m,震源深度为9 km,平均滑动角约为-76°,最大形变中心位于32.85°N、94.86°E附近,地震位错主要分布在深度8~12 km,与余震的深度分布较为一致,反演给出的矩震级为MW5.9,与USGS和CENC测定的矩震级和面波震级基本一致。     

九寨沟MS7.0地震前跨断层水准前兆异常与震后变化

利用2017-08-08九寨沟MS7.0地震前震区北侧青藏块体东北缘跨断层短水准观测资料,以及震后8月、9月两期加密复测结果,综合观测曲线、异常时空演化特征及特征强度指标,分析九寨沟地震前断层形变异常及震后变化。结果显示,九寨沟地震前1 a多西秦岭至六盘山构造区异常比较集中;震前1~3个月西秦岭构造区及其附近(震区周边300 km范围内)短期异常增多,异常特性以逆断性突跳为主,对该次地震有中短期-短期前兆反映;震后1个月西秦岭构造区多数异常调整恢复。