2019年6月17日四川长宁MS6.0地震序列震源机制解与发震构造分析

2019年6月17日,在青藏高原东缘四川盆地南缘宜宾市长宁县发生M_S6.0地震,其后5天内相继发生了珙县M_S5.1、长宁M_S5.3和珙县M_S5.4强余震;7月4日,在珙县珙泉镇再次发生M_S5.6地震.因灾害叠加,本次地震序列导致13人死亡,200多人受伤,大量房屋受损,造成了重大的人员伤亡和财产损失.本文基于四川区域地震台网提供的地震资料,采用多阶段定位方法,对长宁M_S6.0地震序列早期（2019年6月17日至22日）余震进行了重新定位,同时,利用CAP波形反演方法,获得了序列中截止至7月4日的16次M_S≥3.6地震的震源机制解与震源矩心深度,对该序列的发震构造进行了初步分析.长宁M_S6.0地震序列重新定位后的610次M_L≥1.5地震分布显示余震区呈NW-SE向展布,长约25 km,宽5 km;序列震源深度在0~10 km区间,深度均值约3.2 km,但空间上呈西深东浅的分布特征.长宁M_S6.0地震位于余震区的东南端,具单侧破裂特征.CAP波形反演结果显示长宁M_S6.0地震序列以逆冲和逆冲兼走滑型地震为主;16次M_S≥3.6地震的震源矩心深度在1~7 km范围,平均深度3.5 km,与定位结果一致,揭示本次长宁地震序列发生在上地壳浅部.根据序列空间分布、震源机制解及震区构造特征,推测本次长宁M_S6.0地震序列的发生可能与长宁—双河复式大背斜中白象岩—狮子滩背斜和双河场褶皱及其伴生断层活动有关,位于余震区西北段的6月17日珙县M_S5.1、22日珙县M_S5.4及7月4日珙县M_S5.6地震应为6月17日长宁M_S6.0地震触发白象岩—狮子滩背斜伴生断层活动所致.序列发震构造整体呈NE-SW向挤压为主、兼具一定NW-SE向拉张分量的构造变形特征,与南侧2018年12月16日兴文M_S5.7和2019年1月3日珙县M_S5.3地震所呈现的NW-SE向挤压、NE-SW向拉张构造变形特征具有显著差异,揭示四川盆地南缘地带处于构造变形模式的转换区域,所处构造环境的变化导致本次长宁地震序列震源区及附近区域发震构造变形特征具有复杂性.     

基于升降轨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, 同震位错主要表现为正断运动。     

2022年芦山MS 6.1、马尔康MS 6.0地震井水位、水温同震响应特征

基于全国地震地下流体台网数据,分析了芦山M_S 6.1、马尔康M_S 6.0地震引起的地下流体井水位、水温同震响应特征。结果表明,对于芦山M_S 6.1地震水位同震响应观测井数量较多,以上升变化为主,同震变化幅度较大;而对于马尔康M_S 6.0地震水位同震响应观测井数量较少,以振荡为主,变化幅度较小;2次地震引起的水位同震响应能力均强于水温测项,水温记录同震响应的前提是同井能记录到水位同震变化;2组地震引起的同震响应特征差异主要与震源机制解、台站分布密度、同震响应机理不同等有关。     

基于Sentinel-1A数据的四川长宁MS6.0地震同震形变场分析及断层滑动分布反演

利用大地测量数据研究2019年6月17日四川长宁M_S6.0地震同震形变场特征和发震断层参数, 基于DInSAR技术处理升降轨Sentinel-1A数据获取干涉相位图, 并考虑大气折射效应和余震形变误差实现同震形变场改正。四叉树采样后的形变数据作为反演数据源, 采用弹性半空间分层模型反演发震断层几何面滑动分布。结果表明本次地震发震机制为兼具逆冲和左旋走滑, 矩震级为M_W5.9, 断层破裂尺度达28 km×20 km, 震源深度约9.4 km。升降轨视线向同震形变场在断层两侧呈现形变特征差异, 最大沉降量分别是8.34 cm(升轨)和4.23 cm(降轨), 最大抬升量分别是5.5 cm(升轨)和7.5 cm(降轨); 发震断层走向为302°, 倾角为43°, 平均滑动角为50°, 断层面最大滑动量达到0.28 m。     

2020年2月18日长清MS4.1地震发震构造研究

本文利用基于波形互相关的双差定位方法对2020年2月18日长清M_S4.1地震序列进行了精定位计算, 共得到33个地震事件的精定位结果。 结果显示, 地震序列主要沿NW向分布, 在水平方向上具有自NW向SE迁移, 在深度上具有由浅向深迁移的特征; 序列震源深度主要集中在2~7 km, 其中, 主震的震源深度约2.8 km。 由于长清地震序列的地震数量较少, 为了更准确地了解长清地震序列的发震构造、 探索该序列的发生和发展过程, 本文采用CAP方法反演了主震的震源机制解, 其中, 节面Ⅰ走向223°、 倾角42°、 滑动角-160°, 节面Ⅱ走向117.9°、 倾角76.8°、 滑动角-49.8°, 最佳拟合震源矩心深度约2.8 km, 矩震级M_W4.2。 结合区域构造特征分析认为, 长清M_S4.1地震的发震断裂为孝里铺断裂和东阿断裂之间发育的一条浅层次生断裂。 在ENE向区域应力场作用下, 发震断裂产生高角度正断滑动, 并伴有左旋走滑分量, 从而引发长清地震序列。     

2020年7月12日唐山古冶MS 5.1地震震源参数

利用河北及邻区地震台网提供的震相观测报告,使用双差定位方法,对2020年7月12日唐山古冶M_S 5.1地震序列进行重新定位,并基于部分地震台站记录的波形资料,采用近震全波形方法,反演得到主震震源机制解。精定位结果显示,此次地震序列展布长度约8 km,余震自主震位置向SW扩展,震源深度优势分布范围在10—18 km,发震断层面较陡,倾向SE。震源机制解显示,主震为一次走滑型事件,结合序列展布形态和地震活动背景,SW—NE向近似直立节面为可能发震断层面,与1976年唐山M_S 7.8主震断层特征基本一致。综上所述,古冶M_S 5.1地震为唐山M_S 7.8地震的又一次余震,属唐山老震区正常地震活动。     

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

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

2021年9月4日新疆皮山MS 5.1和9月5日叶城MS 5.0地震总结

系统梳理2021年9月4日新疆皮山M_S 5.1和9月5日叶城M_S 5.0地震发生前地震活动异常特征、地球物理观测异常以及区域构造情况,结果如下：①地震活动异常：皮山M_S 5.1地震前出现前兆震群活动、阿图什至皮山4级地震平静、库车—沙雅4级地震窗口以及地震发生率指数异常;②地球物理观测异常：2项异常,其中形变1项,电磁1项,但与此次地震活动相关性较低;③震源机制：泽普断裂距皮山M_S 5.1和叶城M_S 5.0地震最近,间距分别约15 km和3 km,震源机制解显示2次地震均为逆冲型破裂;④序列分析：此次地震活动序列类型为震群型,主震前存在明显的前震活动,余震较为丰富,余震活动呈阶段性衰减特征。综合分析认为,2021年9月皮山M_S 5.1和叶城M_S 5.0地震组成震群型地震序列,震前地球物理观测异常较少,地震活动异常主要对发震地点具有预测意义。     

2016年射阳MS 4.4地震前后重力场变化

利用江苏重力测网2014-2017年重力场观测资料,采用绝对重力控制与相对重力联测相结合的平差方法,获取2016年射阳M_S 4.4地震前后重力场变化图。根据射阳M_S 4.4地震前后射阳地区各测线重力段差变化特点,绘制重力场等值线并进行对比分析,结合相关机理,探讨重力场变化与该地震的内在联系。同震观测数据显示：射阳M_S 4.4地震发生在重力异常值高梯度带附近,发震时震中地区位于NS挤压正异常、EW张拉负异常状态,震后区域重力梯度变化量开始减小,是一种典型的重力异常调整现象。     

2015年7月3日皮山6.5级地震发震构造初步研究

基于新疆区域数字地震台网记录,采用CAP（Cut and Paste）方法反演了2015年7月3日皮山6.5级主震和部分M_S3.6以上余震的震源机制解和震源深度;采用HypoDD方法重新定位了序列中M_L2.5以上地震序列的震源位置,并利用小震分布和区域应力场拟合了可能存在的发震断层面参数.基于上述研究,综合分析了皮山6.5级地震序列的震源深度、震源机制和震源破裂面特征,探讨可能的发震构造.结果显示,利用CAP方法得到的最佳双力偶机制解节面I：走向280°/倾角60°/滑动角90°;节面Ⅱ：走向100°/倾角30°/滑动角90°,矩心深度19 km,表明该地震为一次逆冲型地震事件.大部分M_S3.6以上余震震源机制与主震具有一定的相似性.双差定位结果显示,M_L2.5以上的余震序列主要分布在主震的西南方向,深度主要分布在0～15 km范围内,余震分布显示出与发震构造泽普隐伏断裂一致的倾向南西的特征.利用小震分布和区域应力场拟合得到发震断层参数为走向104°/倾角34°/滑动角94°,该结果与主震震源机制解中节面Ⅱ的滑动角较为接近,绝大多数余震发生在断层面附近10 km左右的区域.根据本研究得到的震源机制、精定位结果以及利用小震分布和区域应力场拟合得到的断层面的参数,结合震源区地质构造情况,初步给出了此次皮山6.5级地震的发震模式.     

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;同震水平位移场模拟结果显示贡嘎错断裂、康西瓦断裂和普鲁断裂等不同位置主应变特征具有差异性,这种差异特征是否影响断裂带以及周围区域的应力构造特征,值得关注.     

2018年5月6日称多5.3级地震序列重定位及发震构造分析

采用双差定位法对2018年5月6日称多5.3级地震及其余震序列进行重新定位,至2018年7月15日共获取129个地震的重新定位结果。结果显示,称多5.3级地震序列主要呈NWW或NNW向分布,其中长轴沿NWW向展布,长约11 km,震源深度主要分布在6-12 km范围内,优势分布为8-11 km。此次地震的震源机制解为走滑型,最佳波形拟合深度为10.1 km。结合精定位、震源机制等综合分析,认为主破裂面走向呈NNW向,发震构造应为巴颜喀拉山主峰断裂。     

九寨沟M_s7.0地震的InSAR观测及发震构造分析

2017年8月8日四川省九寨沟县发生M_s7.0地震.本文基于Sentinel-1 SAR影像,利用InSAR技术获取了此次地震的同震形变场,反演获得同震滑动分布,计算了同震位错对余震分布和周边断层的静态库仑应力变化,并对发震构造进行了分析讨论.结果表明:①InSAR同震形变场显示,九寨沟地震造成地表形变最大量级约为20 cm(雷达视线方向),同震形变存在非对称性分布特征.②同震位错以左旋走滑为主,主要发生在4~16 km深度,最大滑动量约为77 cm,位于9 km深处.反演得到的矩震级为Mw6.46.同震错动未破裂到地表.③大部分余震发生在库仑应力增加区.此次地震增加了震中周边地区一些断裂的库仑应力,如东昆仑断裂带东段、龙日坝断裂、虎牙断裂等.④东昆仑断裂东段的未来地震危险性值得关注.⑤九寨沟地震的发震断层为树正断裂,可能是虎牙断裂的北西延伸隐伏部分,此次地震是巴颜喀拉块体南东向运动受到华南块体的强烈阻挡过程中发生的一次典型构造事件.     

1966年邢台强震群震前、 同震和震后断层运动特征分析

在1966年3月邢台强震群水准观测基础上, 本文首先对同震位移和数值模拟结果进行了对比, 认为邢台强震群发震断层浅层是牛家桥阳断层和永福庄断层, 深部发震断层是倾向北西的北东向断层; 然后以区域水准观测为约束, 通过反演给出了邢台强震群震前断层无震滑移、 同震断层位错和震后断层余滑空间分布。 结果表明, 震前显著无震滑移主要发生在3月8日6.8级前震的发震断层上, 而3月20日宁晋7.2级和6.7级地震发震断层深部同震位错相对弱的断层段, 震前也存在较为明显的无震滑移; 震后断层余滑主要分布在束鹿凹陷西侧的牛家桥阳断层。 从邢台强震群过程中断层运动时空演化特征可以看出, 强震破裂成核过程中震前断层无震滑移既可能发生于前震发震断层, 也可能发生于主震发震断层, 震后断层位移空间分布与同震位错具有互补性; 而震前断层无震滑移是否为有关震前前兆的真正原因, 是值得进一步深入研究的问题。     

基于GNSS与InSAR约束的九寨沟MS7.0地震滑动模型及其库仑应力研究

2017年8月8日的九寨沟M_S7.0地震发生在岷江断裂、塔藏断裂及虎牙断裂交汇地区,地处青藏高原东北部的川甘交界地区,位于巴颜喀拉地块的东缘,地质构造复杂,对于九寨沟地震震中位置和发震断层的确定,存在不同意见.本文利用GNSS及升降轨InSAR观测,在获取九寨沟地震同震形变场的基础上,基于均匀弹性半无限位错模型,联合反演了发震断层的滑动分布模型,并计算了同震库仑应力变化.InSAR同震形变场显示,视线向最大沉降量和抬升量分别为0.21 m和0.16 m,形变场长轴为NW向,形变主要集中在断层西侧.距震中40 km和65 km的九寨和松潘两县,水平向的GNSS同震位移分别达14.31 mm和8.22 mm.联合GNSS和InSAR同震形变场反演得到的滑动分布主要集中在沿走向5~33 km,倾向2~20 km的范围内,平均滑动量为0.18 m,最大滑动量为0.91 m.发震断层长40 km,宽30 km,走向155°,倾角81°,滑动角-9.56°.同震位移场及滑移分布模型表明此次地震为一次左旋走滑为主的地震事件,地震破裂并未完全到达地表,与虎牙断裂北段的几何产状和运动学性质更为接近,结合精定位余震的分布,我们确定虎牙断裂北段为此次地震的发震断层,震中位于北纬33.25°,东经103.82°,震源深度10.86 km,矩震量为7.754×10¹⁸ Nm,相应的矩震级为M_W6.5,与美国地调局和哈佛大学给出的震源机制解基本一致.同震库仑应力导致了虎牙断裂北段延长线的东北和西南两端应力增强,其中塔藏断裂的罗叉段和马磨段未来强震的危险性值得关注.