2017年四川九寨沟MS7.0强震的余震重定位及主震震源机制反演

本文采用双差定位法对2017年8月8日至10月31日期间四川九寨沟M_S7.0主震及5200个余震序列进行相对定位,得到4036个重定位地震事件.采用中国区域地震台网观测到的宽频带垂直分向波形数据和W震相反演方法,得到了主震震源机制解.重定位结果显示,余震序列分别沿NNW和SSE两个方向扩展,展布长度约58 km,且这些余震主要集中在22 km深度之上.余震分布的另一个重要特点是具有分区特性,即在主震NNW方向约5 km处存在明显的西北和东南两区余震活动分界线;西北区的余震由深至浅具有较好连续性,而东南区却在约10 km深度处存在不连续性.余震分布的这种分区特征,说明九寨沟地震震源区的地壳结构存在强烈的不均匀性.余震分布与主震破裂特征的一致性,证实了我们定位结果的可靠性.主震的震源机制解展示出节面Ⅰ的走向/倾角/滑动角分别为246°/83.7°/-177°,而节面Ⅱ的走向/倾角/滑动角为155.7°/87.1°/-6.3°,最佳质心深度为15.5 km,矩震级M_W为6.5.根据余震分布较为垂直和主震震源机制解两节面的倾角均在80°以上,并结合野外地质调查结果,推测此次九寨沟地震为与节面Ⅱ参数相近的一次高角度的左旋走滑型事件.     

历史大地震断层滑动模型的建立及其对同震数值计算的影响——以1920年宁夏海原MS8.5大地震为例

大地震的发生会引起区域位移场和应力场发生变化,进而改变区域内及临近断层的应力状态和地震活动性.目前,研究学者可据已有的断层滑动模型来计算分析大地震同震应力变化,同时采用库仑应力触发理论来进一步分析震后余震分布和断层危险性.然而,历史上曾经发生过不少大地震,例如,1920年的海原M_S8.5大地震,是全球范围内少见的特大地震之一.局限于无确切的地震台站地震波等资料,前人在研究历史地震的影响时往往给出一些简单的断层滑动模型,将断层面上错动量视为均匀分布.为更准确地了解历史地震对后续地震的影响,基于前人研究和一般地震滑动形态分布规律及地震反射剖面等资料,以海原M_S8.5大地震为例,探讨了如何建立海原大地震断层滑动模型,并分别搭建了简单断层滑动模型和复杂断层滑动模型的全球同震横向不均匀并行椭球型地球模型.通过对海原M_S8.5地震的同震位移场和应力场的计算,发现采用复杂断层滑动模型比简单断层滑动模型地表位错分布更切合实际.同时,进一步计算和分析了此次大地震对青藏高原东北缘近100年历史地震和周围断层的应力触发作用,得出断层滑动模型对同震计算结果的影响集中在发震断层附近而对远场影响较小.     

2001年昆仑山口西MS8.1地震对周围断层的应力影响数值分析

大地震的发生往往会引起周围区域形变场和应力场变化,且对临近断层上的应力状态也有影响.2001年11月4日,昆仑山口西发生了半个世纪以来中国最大的M_S8.1级地震.本文基于已有的滑动模型,建立了三维含地形高程的横向不均匀性椭球型地球有限元模型,采用等效体力方法,分析了此次M_S8.1地震产生的全球同震位移和应力场变化.与解析方法相比,该模型考虑了地形、Moho面起伏和地球介质横向不均匀性;与一般的有限元数值模拟相比,该模型考虑了地球曲率和椭率,合理地规避了有限块体模型假定边界位移为零所引入的误差.计算得出同震位移与GPS观测数据可以很好地吻合.据库仑破裂应力准则和震源参数,计算得出昆仑山口西M_S8.1地震的发生造成了汶川、芦山、改则和当雄地震的发震断层上库仑应力增加,对这些地震的发生起促进作用;而造成玉树和德令哈地震发震断层上的库仑应力变化为负值,在一定程度上抑制了这些断层的地震活动性.此外,计算结果显示地球地形高程、介质非均匀性和椭率对昆仑山口西M_S8.1地震同震变化计算有一定的影响,其中地形和椭率造成的同震位移场相对误差约10%.     

深度神经网络拾取地震P和S波到时

从地震波形数据中快速准确地提取各个震相的到时是地震学中的基础问题.本文针对上述问题提出了利用深度神经网络拾取到时的新方法,建立了用于地震到时提取的17层Inception深度网络模型,在对原始三分量数据进行高通滤波和归一化处理后输入网络直接输出到时信息.整个过程基于神经网络自适应提取波形特征,自动输出结果.通过对100组加了不同强度的噪声数据进行了可靠性检验,相比于其他方法神经网络方法对于噪声具有较高的容忍度以及稳定性,并且与地震目录数据有较高的相似性.相比于AR-AIC+STA/LTA,深度神经网络虽然运算速度稍慢,但整个过程不需设定时窗与阈值,同时具有更高的可用性,并且可以迭代升级以提高精度.此方法作为人工智能方法,为波形到时拾取提供了新思路.     

InSAR大气误差改正及其在活动断层形变监测中的应用

为了获取青藏高原东北缘老虎山断裂带精确的震间形变速率场,通过对短时间基线干涉图大气改正效果的评价,从3种外部大气数据(MERIS, ERA-Ⅰ, WRF)中确定出最优的大气改正方法,用于长时间基线干涉图中的大气信号改正;然后利用层叠法(stacking)累积平均经大气和轨道改正后的干涉图,获取了研究区的震间形变速率场.结果显示:海原断裂系统区域内,MERIS和ERA-Ⅰ的大气改正效果优于WRF;MERIS和ERA-Ⅰ的改正结果给出了相似的形变速率场,断层两盘相对形变速率为视线向2.5 mm/a,转换成平行于断层方向为6.5 mm/a,与GPS结果一致;在近断层5 km的范围内,出现了较大的形变梯度,揭示了浅层蠕滑的存在.      

利用初至P震相测定石嘴山ML4.4地震序列的震源深度

基于初至P震相较易辨识、读取误差小、有效震相丰富等优点,提出联合Pg与Pn震相到时数据测定震源深度的方法,并将其应用于石嘴山M_L4.4地震序列(11次)的震源深度计算,结果显示该地震序列的震源较浅,这也是此次有感地震震感强烈的原因之一．然后利用初至P震相到时方法分别计算了双层和4层地壳速度模型下的震源深度,并与双差定位所得深度进行对比,结果显示: 利用初至P震相到时方法较双差定位方法得到的震源深度整体上一致性较好,但前者得到的震源深度较之后者略深;对于双层与4层地壳速度模型,二者的震源深度计算结果具有较好的一致性．因此,基于可靠区域地壳速度模型下的初至P震相到时方法可以应用于震源深度的计算,并能够取得较好的结果．另外,石嘴山M_L4.4地震序列的主震深度为7—8 km,最大余震深度为6 km,根据银川盆地的速度成像结果可知,主震和最大余震发生于银川盆地基底下部;其余地震的震级偏小,多集中于3—5 km深度,主要发生于银川盆地基底顶部或基底覆盖层内．      

气汞的第一个同震效应记录

利用杭州超距科技有限公司最近研发的ATG-6138型测汞仪在云南省弥勒井气汞（Hg）观测中记录到尼泊尔8.1级地震的同震效应。这是我国气汞观测的第一个同震效应的纪录。这一事实表明了进一步提升地震地下流体化学量观测仪器的观测精度与采样率,有可能记录到固体潮与同震效应等更多的地壳动力作用信息,从而进一步提升其地震前兆监测能力。     

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甲仙地震也同该断层有关.     

九寨沟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.同震错动未破裂到地表.③大部分余震发生在库仑应力增加区.此次地震增加了震中周边地区一些断裂的库仑应力,如东昆仑断裂带东段、龙日坝断裂、虎牙断裂等.④东昆仑断裂东段的未来地震危险性值得关注.⑤九寨沟地震的发震断层为树正断裂,可能是虎牙断裂的北西延伸隐伏部分,此次地震是巴颜喀拉块体南东向运动受到华南块体的强烈阻挡过程中发生的一次典型构造事件.     

利用近震高频接收函数研究四川理县西山村滑坡体结构

远震接收函数已广泛用于反演台站下方的结构,然而由于地球的非弹性衰减作用,远震数据较难获得高频接收函数,对浅地表结构约束不足.为了克服这一问题,我们使用近震数据的高频接收函数来研究浅表速度结构,并应用于四川理县西山村滑坡体上3个宽频带地震仪记录到的近震事件.本文发展了接收函数V_P-k(V_P/V_S)叠加方法,结合接收函数H-k叠加和波形反演方法获得了台站下方滑坡体的厚度、S波速度和平均V_P/V_S比,并与钻孔得到的滑坡体厚度进行对比.结果表明,滑坡体具有小尺度的横向不均匀性,台站下方滑坡体的平均V_P/V_S比在2.4~3.1之间变化并且在底层存在78~143m·s~(-1)左右的S波低速层.本文观测到的高V_P/V_S比和底层低的S波速度结构,与电磁法获得的滑坡体底层低的电阻率和底部富水特征一致,表明滑坡体h1底界面的抗剪强度相对较弱,是潜在的滑坡危险区域.本文研究结果表明,利用近震接收函数能有效约束浅表的速度结构,进而能为滑坡灾害治理提供一定的地震学参考.