2017年8月8日九寨沟7.0级地震长波辐射异常特征

通过距平方法,研究2017年8月8日九寨沟7.0级地震前震中所在区域（95.00°—110.00°E,25.00°—45.00°N）长波辐射时空演化特征,研究结果表明:①2017年7月,去除背景之后的长波辐射场在震区附近出现显著增强现象,增强区域基本走向与地质构造走向一致,其主体区域沿着巴颜喀拉块体南缘边界带,重要分支横跨巴颜喀拉块体,直接延伸至九寨沟7.0级地震震中;②紧邻九寨沟7.0级地震震中的4个格点在去除背景变化后的长波辐射时序曲线变化特征基本一致,即在2017年7月出现显著大于其他月份的现象。     

2017年四川九寨沟MS7.0地震及余震精定位研究

2017年8月8日21点19分四川省九寨沟县发生M_S7.0级强震,本文应用对相对位置和绝对位置都有较好约束的hypoRelocate精定位方法,对8月15日前M_L大于1.5的854个余震重新定位.重定位使用了6178个Pg走时、3363个Sg走时、158929个地震对的相对走时差以及74个M_L大于3.5的34个地震构造的SH尾波互相关走时.定位结果显示,余震区呈NW向分布,长度约35 km、宽度约12 km.以主震为中心,余震形成西北和东南两个丛集区,其中西北区的余震展布较东南区宽,发震断层近似垂直,与主震的震源机制解吻合.余震深度主要集中在5~15 km,主震深度为16 km.从沿断层走向的垂直剖面可见,主震周围一直延伸到地表的区域,余震相对稀少,出现明显空区,这种现象可能与主震发生时该区域的滑动位移较大,应力得到充分释放有关.另外,我们还发现震后60 h之后余震的深度明显变浅,且西北向的深度明显浅于东南向.通过对比分析流动台站布设前后余震深度的分布特点,我们认为这种时空演化差异主要是由流动台站陆续布设和使用造成的,对余震的主要分布特征影响不大.     

2017年九寨沟MS 7.0地震震源深度测定

为了更好地确定2017年8月8日九寨沟M_(S )7.0地震震源深度其发震机理,利用四川、甘肃和青海区域地震台网的观测波形数据,采用多种方法研究了此次地震的震源深度。首先,采用gCAP方法反演了九寨沟M_(S )7.0地震的震源机制解和矩心深度,结果显示,节面Ⅰ走向243°/倾角87°/滑动角-158°,节面Ⅱ走向151°/倾角68°/滑动角-3°,矩震级为M_(W )6.5,矩心深度为8 km;然后,采用ISOLA近震全波形方法反演了此次地震的震源机制解,反演结果与gCAP方法结果相差不大,矩心深度为7 km;最后,通过sPn震相与Pn震相之间的走时差测定此次地震初始破裂震源深度,结果显示深度约为12 km。研究表明,九寨沟M_(S )7.0地震的矩心深度为7—8 km,初始破裂深度约为12 km。     

2017年8月8日九寨沟MS7.0地震震源谱及震中区域品质因子

2017年8月8日四川省九寨沟县发生M_S7.0地震,国家强震动台网中心正式发布了此次地震中收集的66组强震动记录,采用其中震源距不超过150 km且台站场地反应已确定的8组记录,基于S波傅氏谱反演方法确定了此次地震的震源谱及震中区域的地壳介质品质因子.结果为：品质因子Q（f）=84.9f^0.71,说明该区域的S波非弹性衰减较为强烈;根据ω^-2震源谱模型,确定该地震的地震矩为9.42×10¹⁸Nm,其对数标准差为0.12,相应的矩震级M_W=6.616±0.079,拐角频率为0.131±0.011 Hz,应力降为3.854 MPa,该值明显低于全球板内地震的应力降平均水平（4.89 MPa）,但高于2013年芦山M_W6.6地震的应力降.

     

2017年8月8日四川九寨沟M 7.0地震震源深度测定

2017年8月8日四川阿坝州九寨沟县发生M 7.0地震,随后发生一系列余震。该地震发生在川甘交界处四川一侧。国内外科研机构对主震的震源深度给出不同结果。本文采用PTD方法,采用四川省地震台网测定的震相数据,基于四川2015地壳速度模型,计算得到M 7.0主震震源深度为13 km,利用甘肃省地震台网震相数据,基于甘肃2015地壳速度模型,计算得到主震震源深度为12 km。震后流动台架设后,由近台证明PTD测定深度结果的可靠性。另外,计算M 3.0以上余震的震源深度。所有震相数据来自中国地震台网中心编目数据库。     

采用IDS方法反演强震数据确定2017年8月8日九寨沟地震的破裂过程

近场强震数据是地震观测的主要数据类型之一,但尚未广泛用于震源参数快速和常规测定.作者利用2017年8月8日九寨沟地震的近场强震数据,分别基于中国地震台网中心和美国地质调查局的定位结果,采用迭代反褶积和叠加（Iterative Deconvolution and Stacking,IDS）方法反演得到了两个破裂模型.反演结果显示,两个破裂模型除在时间上存在大约2.5 s的偏移外,在其他破裂特征上都极其相似,表明强震数据反演并不严重依赖地震定位结果.破裂过程大约持续8~10 s,释放的地震矩在6.9×10¹⁸ Nm左右,对应的矩震级约为M_W6.5.走向方向上,破裂覆盖了震中西北15 km至震中东南10 km的区域,整体表现为不对称的双侧破裂模式,其中西北方向的破裂略占优势.深度方向上,主要破裂集中在0~10 km的较浅区域.这些特征与目前的地震学和大地测量学结果基本一致,表明在当前的数据条件下,即使定位结果存在不确定性,通过反演强震数据也能够独立可靠地确定中强地震的破裂过程.特别地,未来强震数据的数据质量和获取效率得到进一步改进和提升后,这一工作有望成为震后快速确定破裂过程的主要方式之一,在震后应急工作中发挥积极作用.

     

近震全波形反演2017年九寨沟M7.0地震序列震源机制解

搜集了四川地震台网的波形资料,采用全波形反演2017年8月8日九寨沟M7.0地震序列震源机制解.反演结果显示,九寨沟主震矩震级为M_W6.36,震源深度为22 km,节面Ⅰ走向为150°,倾角为80°,滑动角为-20°;节面Ⅱ走向为244°,倾角为70°,滑动角为-169°.余震主要分布在14~22 km深度范围内,震源机制以走滑型为主,其中正断型地震2个,逆冲型地震2个,走滑型地震24个,混合型地震8个.断层面优势方向为SSE向,与塔藏断裂和虎牙断裂走向基本一致,但与塔藏断裂最南段存在明显差异.倾角变化集中在60°~80°,滑动角主要分布在0°附近,表明九寨沟地震序列主要受SSE走向、近似直立的左旋走滑断层控制.P轴优势方位为SEE向,仰角主要分布在30°以内,与区域应力场基本一致.震源区的机制类型和应力状态均存在空间分段差异.本文推测此次九寨沟M7.0地震序列可能发生在虎牙断裂向北延伸的隐伏断裂上,但不排除地震引起了塔藏断裂南段和虎牙断裂以北隐伏断裂同时破裂的可能.

     

2017年8月8日四川九寨沟M7.0地震及其余震序列的震源参数

2017年8月8日21时19分46秒,四川北部阿坝州九寨沟县发生M7. 0地震。该地震发生在东昆仑断裂带东段南端的一条分支断裂上,地处玛曲—玛沁地震空区和松潘—平武地震破裂区的过渡部位。为了解这次地震的震源参数特征以及震区构造环境,文中基于四川数字波形资料及中国地震台网中心的震相数据,采用Loc3D定位方法获取初始地震目录参数,并利用双差定位方法重新定位地震序列中的较小地震,选用CAP(Cut-and-Paste)方法分析29次中强地震的震源机制解与震源矩心深度,同时对186次地震(2. 0≤ML≤5. 5)进行了震源谱恢复,绘制出余震序列的震源应力降空间分布图像。重新定位结果显示,该地震序列以M7. 0主震震中为界表现出明显的空间分段差异:余震区的NW段地震数量较少、震源应力降较大、震源机制解类型复杂; SE段分布着多数4. 0级以上中强地震,且几乎全部为纯走滑型地震,余震序列的应力降随时间逐渐减小。结合震源机制解的研究结果分析,发震构造可能为NW-SE向、断续的隐伏断裂或虎牙断裂的北延推测段,断层面在主震南、北两侧发生明显转折,同时由挤压逆冲为主逐渐转换为以走滑运动为主。     

对2017年九寨沟MS7.0地震的预测分析

2017年8月8日四川九寨沟发生M_S7.0地震。震前我们实验室主要做出以下预测：（1）根据我国强震的空间分布有序性规律,在2015年曾做出预测：甘肃省东南部、四川交界处将来可能会发生8级地震;（2）根据44年周期和7年周期,2016年预测四川区域在2017年会有7.5级地震发生;（3）依据地磁仪器观测到的数据异常,2017年7月28日认为近期将有6~7级地震发生。九寨沟地震发生后,经综合分析我们的三次预测,以及6—7月在四川茂县出现的多次滑坡事件,认为若将这些因素结合起来,将会是一个很好的短临预测方式。     

2017年8月8日四川九寨沟7.0级地震强震记录及特征分析

2017年8月8日四川九寨沟县发生7.0级地震,中国数字强震动台网布设在四川、甘肃、陕西、宁夏的66个强震台获得主震加速度记录。本文首先对198条三分向强震记录进行常规处理,计算出近场强震记录的加速度峰值随震中距的分布情况;再根据2个典型台站的加速度时程记录,通过计算其加速度反应谱并与设计反应谱比较,分析本次地震的基本特征;然后将实际观测数据与意大利新一代地震动衰减公式对比,分析峰值加速度（PGA）及谱加速度的衰减关系;最后结合已有的工程场地钻探资料,采用H/V谱比法对2个不同类别的典型台站进行地场地效应分析,发现该方法能很好地反映实际台站场地的反应特征。     

应用地震干涉法定位四川九寨沟7.0级地震震源位置

利用地震干涉法对2017年8月8日九寨沟M7.0主震及部分余震进行定位.地震干涉法在可控源地震勘探领域有着广泛应用,但由于震源机制和震源辐射花样的复杂性,在天然地震研究中,地震干涉震源成像应用较少.本研究将地震干涉震源成像技术应用到了天然地震定位中,通过计算原始地震波形记录的特征函数,消除了由于震源辐射花样不同而引起的初动方向不一致性.通过对干涉波形进行偏移叠加,选择互相关型和卷积型两种偏移核函数,分别对震源水平位置和深度进行偏移成像,确定了九寨沟M7.0主震及11次M > 3.5余震的震源位置参数.结果显示,这次地震发生在一个之前未被关注的无名断裂（现已命名为树正断裂）上,研究分析推测,树正断裂与虎牙断裂存在连接的可能性.

     

2008年汶川地震同震-震后应力演化及其对2017年九寨沟M_S 7.0地震的影响

2017年8月8日四川九寨沟发生了MS7.0地震,距离汶川地震震中约为140km.地震的发生会调整区域应力分布和影响区域地震活动性变化.那么,汶川地震的发生对此次九寨沟地震有何影响?区域应力如何演化?为回答这些问题,本研究基于高性能并行有限元方法和网格自适应加密技术,采用含地形起伏的黏弹非均匀椭球型地球模型,分别计算了汶川地震引起的同震-震后应力时空动态演化及其对此次九寨沟地震的影响.同时,考虑同震静态应力和震后黏弹性应力调整,计算了汶川、芦山和九寨沟地震的发生对周边断层应力积累的影响.计算结果显示汶川地震对九寨沟发震断层的同震库仑应力加载约为0.008 MPa,震后9年的黏弹性应力加载约为0.012~0.016 MPa,可能会使九寨沟地震提前发生.此外,汶川、芦山和九寨沟地震的发生引起龙门山断裂映秀以南段、东昆仑断裂、龙日坝断裂东段和岷江断裂北段及鲜水河断裂康定段应力加载大于0.01 MPa,增加了这些断裂带发生地震的危险性.     

2017年四川九寨沟MS7.0地震前区域重力场变化

利用南北地震带2014-2017年期间的流动重力观测资料,系统分析了区域重力场变化及其与2017年8月8日四川九寨沟7.0级地震发生的关系.结果表明：① 区域重力场异常变化与北西西向塔藏断裂和南北向岷江断裂带在空间上关系密切,反映了沿控震断裂在2013-2017年期间发生了引起地表重力变化效应的地壳变形和构造活动.② 九寨沟地震前,测区内出现了大空间范围的区域性重力异常,而震源区附近产生了局部重力异常,沿塔藏断裂带形成了重力变化高梯度带,其中,甘肃玛曲、迭部、青海河南蒙古族自治县、四川若尔盖、九寨沟一带重力差异变化达100×10^-8m·s^-2以上;这些可能反映九寨沟地震前,区域及震源区附近均产生与该地震孕育、发生有关的构造运动或应力增强作用.③ 九寨沟地震震中位于重力差异运动剧烈的鞍部等值线附近,与断裂走向基本一致的重力变化高梯度带零值线上,这一观测事实进一步佐证了重力场动态变化图像对强震地点预测具有重要的指示意义.

     

东昆仑断裂带东端的构造转换与2017年九寨沟MS7.0地震孕震机制

2017年四川九寨沟M_S7.0地震是继2008年汶川M_S8.0地震和2013年芦山M_S7.0地震之后,青藏高原东缘在不到十年的时间内发生的第三个震级M_S7.0以上的强震.这次地震发生在东昆仑断裂带东端,作为青藏高原东北缘的一条大型左旋走滑断裂带,东昆仑断裂带与东端其它构造之间的转换关系仍不清楚,因区内地质构造和地形复杂,东昆仑断裂带东端的主要构造仍缺少深入的研究.本文在总结区域地震构造活动特征、历史地震和现代地震基础上,通过东昆仑断裂带东端已有的和最近开展的活动构造定量研究结果,并结合现今GPS变形场资料和2017年九寨沟M_S7.0地震灾害特征分析,发现东昆仑断裂带最东段塔藏断裂上的左旋走滑除了一小部分继续向东传播转移到文县断裂带上外,大部分转化为其南侧的龙日坝断裂带北段、岷江断裂和虎牙断裂上的近东西向地壳缩短,这可能是岷山隆起的构造机制,而2017年九寨沟M_S7.0地震正是左旋走滑的东昆仑断裂带在东端继续向东扩展的结果.

     

GNSS观测的九寨沟7.0级地震同震位移初步结果

2017年8月8日四川省九寨沟县发生了7.0级地震,中国大陆构造环境网络与北斗地基增强系统的GNSS连续观测共同监测到了此次地震的同震位移（坐标：东向为正,北向为正）,结果显示：3个站点记录到了明显的同震位移,距离震中43 km的九寨沟台站（SCJZ）在东西向的位移为-9.8±1.5 mm,在南北向的位移为3.3±0.7 mm;距离震中65 km的松潘站（SCSP）在东西向的位移为-1.8±0.7 mm,在南北向的位移为-7.7±0.6 mm;距离震中77 km的舟曲站（GSZQ）在东西向的位移为0.4±1.2 mm,在南北向的位移为3.6±0.8 mm.通过同震位移分布特征,可以推测此次地震为一次左旋走滑型事件,引起水平向同震位移大致不超过150 km范围,地震对东南侧的龙门山断裂带影响非常小,对北侧的塔藏断裂和西侧的岷江断裂处引起的同震位移为厘米级.同震位移的反演结果显示：断层面上滑动量主要集中在7 km深度,最大量值约为0.4 m,平均滑动角为-15°,利用滑动分布计算的相应矩震级为M_W6.4,与地震波反演结果相当.结合同震滑动分布、同震主应变分布、余震分布和震源机制解等特征,推测此次地震破裂极值区累积的能量得到较充分释放,进一步分析得出此次地震在塔藏断裂、岷江断裂和虎牙断裂处产生了一定的应力变化,值得持续关注.