四川长宁MS6.0地震震源干涉成像定位

文中利用干涉成像定位方法对2019年6月17日四川长宁M_S6.0主震及部分余震进行定位。天然地震波振幅的极性和大小随震源机制和辐射花样的变化而变化,利用原始波形的特征函数可消除震源在不同方位上远场P波的初动极性和大小的不一致性。文中将干涉成像技术应用于天然地震定位,通过对干涉波形进行偏移和叠加处理,分别对震源的水平位置和深度进行偏移成像,确定了主震及较大余震(M_S≥4.0)的震源位置参数,其中主震的位置为(28.38°N,104.88°E),震源深度为8.0km。此外,还测试了4种不同的速度模型对定位的影响,结果显示利用震源干涉成像定位方法获得的结果较为稳定。通过计算台阵、台网的响应函数,评估了台站分布及特征函数周期长短对定位结果的影响。计算得到的特征函数优势周期约为4s,该周期的台阵、台网响应函数在水平和垂直方向均显示了较好的收敛性和稳定性。     

2017年8月8日九寨沟M7.0地震及余震震源机制解与发震构造分析

2017年8月8日在青藏高原东缘四川省九寨沟县发生M7.0级强烈地震,极震区烈度达Ⅸ度,但无明显地表破裂,一定程度上限制了发震构造的确定和后续地震危险性判定.本文基于截止至2017年8月14日的地震资料,采用多阶段定位方法,对主震及余震进行了重新定位,同时,利用CAP波形反演方法,获得了M7.0主震与13次M_L ≥ 4.0级余震的震源机制解和震源矩心深度,进而初步分析了本次地震的发震构造.结果显示,九寨沟M7.0地震的矩震级M_W6.4,震源矩心深度5 km,表明主震发生在上地壳浅部,与2003年伊朗巴姆（Bam）M_W6.5地震特征极为相似;12次M_L ≥ 4.0级余震的震源矩心深度6~12 km,显示这些余震发生在主震下部,仅1次例外.重新定位后的余震震中呈NW-SE向窄带展布,位于近NS向的岷江断裂与近EW向的东昆仑断裂带东端分支塔藏断裂所夹持的区域,余震带长轴长约38 km,主震位于余震带中部.根据余震震中分布、主震及余震震源机制解等,推测本次九寨沟M7.0地震及其余震的主发震构造为位于岷江断裂与塔藏断裂之间的树正断裂.震源机制解揭示,树正断裂呈左旋走滑,走向约152°,近SE,倾向SW,倾角约70°,该断裂应属于东昆仑断裂东端的分支断裂之一,或与东南侧的虎牙断裂构成统一断裂系.     

2017年九寨沟7.0级地震断层破裂尺度初步研究

本文利用2017年九寨沟7.0级地震的余震序列数据,结合震源机制解确定的地震断层面,研究了地震断层的破裂尺度。基于余震序列在断层面上的投影,分析直接余震和间接余震以及不同起始震级条件下余震区的长度和震源深度分布,估算了地震断层破裂面的破裂长度和破裂深度,并探讨了地震断层破裂对地震烈度分布的影响。研究结果显示,九寨沟7.0级地震断层的破裂长度约33—35km,破裂深度约23—26km;地震断层破裂对Ⅷ度区分布有明显的控制作用,地震断层破裂长度接近Ⅷ度区的长度。     

2017年九寨沟MS7.0地震震源区速度结构与余震分布

采用九寨沟M_S7.0 （M_W6.5）地震的余震直达P波、S波走时数据,通过体波走时层析成像方法,获得了震源区及其邻区的P波和S波速度结构,并利用成像结果对余震进行了重定位。结果显示:余震主要集中分布于高、低速异常交界处偏低速异常一侧,呈走向NNW,倾向SW,倾角较高的分布特征;余震序列两侧的P波、S波速度结构揭示了发震断层两侧介质性质的差异,即上盘为刚性较强的高地震波速度区,下盘为刚性较弱的低地震波速度区。由余震分布特征和地震波速度结构推断:九寨沟地震发生在上地壳底部,发震断层具有上盘地震波速度高、下盘地震波速度低的特征;主震引起的后续破裂在上地壳内部的剧烈形变区内传播,破裂能量终止于25 km深度附近。      

2017年九寨沟MS7.0级地震震源过程反演与烈度估计

2017年8月8日,在我国四川省九寨沟县发生一次M_S7.0级地震.在快速响应的基础上,重新筛选远场地震波形资料并收集覆盖震中区的InSAR资料对主震的震源破裂过程重新进行了反演分析,收集震后约13 h的余震震相数据对余震进行了双差定位,并基于此对发震断层的复杂性进行了讨论,提出了有待进一步研究的问题.最后,利用反演得到的有限动态破裂模型对地震烈度进行了估计.

     

基于手机位置数据的四川九寨沟7.0级地震人流分析

2017年8月8日四川九寨沟7.0级震后,浙江省地震局利用大数据采集了1031.99万条手机位置记录和4.46万个空间格网位置。本文结合24hr连续定位的手机数据,使用手机位置数据分别对九寨沟灾区人口从时间和地理维度上进行了量化分析,估计了多维人口分布的偏差,同时,探讨了剔除微观误差数据用户来估算灾区通讯基站退服分布的方法。该项工作为震后快速获取灾区人口实时动态分布提供了有效途径,同时也为地震灾害评估提供了较准确的依据。     

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之后余震的深度明显变浅,且西北向的深度明显浅于东南向.通过对比分析流动台站布设前后余震深度的分布特点,我们认为这种时空演化差异主要是由流动台站陆续布设和使用造成的,对余震的主要分布特征影响不大.     

四川九寨沟7.0级地震微博灾情信息特征分析

通过调用新浪微博官方API的方式,获取了四川九寨沟7.0级地震震后微博信息,对获取到的数据进行了分类统计、微博词频统计与时空特征分析。统计分析结果表明:在分类统计方面,地震发生之后24h内,与地震相关微博中人的反应占比达到73%,救援行动占比达到11%,这时由于微博本身特点以及民众对救援的期望造成的。在词频统计方面,高频名词表明震后24h内微博上的热点事件,而高频动词与形容词表明民众会在震后变得焦虑和不安并通过在网上互相激励来缓解震后的不安和焦虑。时间特征表明震后0~4h内,有大量和地震相关的灾情信息会通过微博博文内容的方式发布,而随着应急救援行动的进行,在震后12h以后,民众的情绪会逐渐平复,社会及民众的关注度也随之降低。空间特征表明,由于震后通信中断与网络堵塞,微博灾情基本成点状分布,随后随着通信与网络的恢复、应急救援的进行微博灾情逐渐变为带状分布并进而趋向于形成一个面。     

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

采用特定震源台站校正定位方法对九寨沟M_S7.0地震及其余震进行了精定位研究,结果显示九寨沟地震主震深度为19km左右,余震主要分布于5~15km深度范围;与单事件定位方法相比,其余震展布和断层的线性特征更吻合,更具有平面特征;与双差定位法相比,精定位事件覆盖范围更广。精定位结果显示,断层整体为虎牙断裂西北向的延伸,根据余震分布和特定震源区的校正项分布,认为此发震断层应分南、北两段;余震区长度为40km左右,西北段宽度为10km左右,东南段宽度为6.5km左右。     

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,与地震波反演结果相当.结合同震滑动分布、同震主应变分布、余震分布和震源机制解等特征,推测此次地震破裂极值区累积的能量得到较充分释放,进一步分析得出此次地震在塔藏断裂、岷江断裂和虎牙断裂处产生了一定的应力变化,值得持续关注.

     

四川九寨沟7.0级地震前震情跟踪概述及震后总结

对于强震可否预测与如何科学预测的问题,通过九寨沟地震前的震情跟踪及地球物理场的时空演化总结,科学地认识地震预测工作.四川九寨沟7.0级地震前距离震中400 km范围内,在1～2年内出现形变、地电阻率等地球物理场的异常.异常时间主要从2016年开始,空间上主要沿北东方向分布,在200 km内有流动水准4项异常、GNSS观测1项;200～400 km流动水准1项异常、定点形变2项,GNSS观测1项,地电阻率3项.这些异常变化出现后,地震部门组织了多次现场核实和跟踪;最为典型的核实是2017年6月27日至29日中国地震局台网中心、宁夏地震局、甘肃省地震局组成异常核实小组,对甘宁川交界地区多台地电阻率变化进行了分析研判,2016年7月平凉台井下观测的地电阻率有1年的趋势性破年变异常,NS和EW向地电阻率异常变化率比背景值(背景值年变化为0.6%、0.4%)的明显增大,分别为1.1%、0.6%;天水深井地电阻率从2月份三个测道地电阻率值偏离背景值幅度达到0.4%,宁夏固原EW测项异常下降幅度约为0.4%,计算分析这组异常是由地层电性变化引起;针对上述变化再结合2015年全国7级地震(M 7)大形...     

利用早期余震预测的Omi-R-J方法对2017年四川九寨沟MS7.0地震的应用研究

在地震序列的早期阶段获得科学可靠的余震预测结果,始终是地震预测研究的前沿课题.针对新近发展的Omi-R-J方法的适用性问题,以及相比于传统的Reseanberg-Jones（R-J）模型在地震序列参数拟合、余震预测效能上的差异等问题,本文以2017年四川九寨沟M_S7.0地震序列为例,利用多时间窗的连续滑动拟合、预测和检验,以及通过构建重采样随机地震序列目录的系统检验等方式,开展了应用研究和比较研究.结果表明：相比于R-J模型,Omi-R-J方法可在余震记录较不完整的地震序列早期阶段获得稳定、可靠的序列参数,p值、c值和b值的标准差相比于R-J模型参数明显减小.N-test方法检验结果表明,在R-J模型等传统方法无法获得预测结果的地震序列早期阶段,Omi-R-J方法预测结果可以较高的比例通过N-test检验,在后续时段的预测效果也明显优于R-J模型.利用随机地震序列目录的测试结果表明,余震记录完整性下降对Omi-R-J方法预测结果的影响相对较小,在全部时段的完整性下降的影响高于在部分时段的完整性下降.上述结果对进一步地将Omi-R-J方法应用在震后早期余震预测中具有一定科学借鉴意义.

     

基于Pn/Pg相对定位方法研究2017年8月8日九寨沟M7.0地震起始破裂深度

本文利用Pn/Pg相对定位方法,测定了2017年8月8日四川九寨沟M7.0主震及部分余震的起始破裂深度.地震的起始破裂深度是理解地震孕震机理的重要参数,而九寨沟地区台网稀疏,地壳速度结构复杂,基于传统的到时定位方法测定地震起始破裂深度误差较大.Pn/Pg相对定位方法首先基于流动观测近台记录对2017年8月10日M4.1和11月7日M4.5余震震源位置进行测定,选择其为参考事件,再利用Pg校正主震水平位置,Pn约束震源起始深度,基于参考事件可以有效降低速度模型对震源位置测定的影响.结果显示：九寨沟主震起始破裂深度约9 km,早期余震的震源深度分布在7~13 km,主要集中在主震起始破裂深度附近.     

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

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

     

九寨沟7.0级地震地表地震动位移及静态位移的模拟研究

北京时间2017年8月8日21时19分在四川省阿坝州九寨沟县（北纬33.2°、东经103.8°）发生7.0级地震,给九寨沟景区造成了滑坡坍塌等地质灾害.本文利用Crust1.0资料和反演的震源破裂过程建立了研究区三维介质模型及三维震源模型,采用谱元数值模拟方法模拟了九寨沟近断层地表地震动速度时程,并与强震台站实际速度记录的低频做了对比,两种结果对比的一致性,验证了介质模型及震源模型的有效性.在此基础上近一步模拟了地表位移时间序列,获得了位移峰值PGD及静态位移的空间分布,分析并给出了九寨沟地震地表位移的主要运动特征.

     

2017年九寨沟MS7.0地震对周围地区的静态应力影响

基于九寨沟M_S7.0地震的破裂模型及均匀弹性半空间模型,本文计算了该地震在周围主要活动断层上产生的库仑应力变化、在周围地区产生的应力场和位移场和同震库仑应力变化对余震的触发.结果表明：（1）九寨沟地震造成虎牙断裂中段库仑破裂应力有较大增加,已经超过0.01 MPa的阈值,虎牙断裂北段、塔藏断裂中段和岷江断裂北段北部的库仑破裂应力有较大降低,因此尤其要注意虎牙断裂中段的危险性.（2）水平面应力场在该地震震中东西两侧增加（拉张）,张应力起主要作用.在震中南北两侧降低（压缩）,压应力起主要作用.从水平主压和主张应力方向来看,均呈现出条形磁铁的磁场形态.从剖面上的应力场来看,在上盘的面膨胀区域内,大部分点的主张应力方向与地表是垂直的,在其他区域内,主张应力和主压应力均以震中为中心,向外呈辐射状.（3）从地表水平位移场来看,震中东西两侧物质朝震中位置汇聚,南北两侧物质向外流出,在震中处的最大水平位移量达43 mm.从地表垂直位移场来看,震中南北两侧出现明显的隆升,隆升最大值达56.8 mm.震中东西两侧出现明显的沉降,沉降最大值达74.5 mm.从剖面的位移场来看,九寨沟地震为左旋走滑地震,且有一定的正断成分.由分析可以推测该断层破裂在大致22~26 km的深度上就截止了.并推测下盘物质运动的动力来自震源北东东方向（四川块体）深度在6~30 km的下盘下层物质,上盘物质运动的动力来自震源北西西方向（巴颜喀拉块体）深度在0~6 km的上盘上层物质.（4）通过计算不同深度上主震对余震的触发作用可知,主震后的最大余震受到了主震的触发作用,多数其他余震也受到主震的触发作用.主震的发生促使了库仑应力增加地区余震的发生,抑制了一部分库仑应力减少地区余震的发生.

     

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地震的破裂空区完全破裂.考虑到塔藏断裂和虎牙断裂的运动性质,可初步判定发震断层为虎牙断裂北侧延伸分支.基于三种不同接收断层模型的同震库仑应力变化计算结果反映出该区域以应力释放为主,进一步触发较大走滑型余震的可能性不大.