2017年漾濞MS4.8和MS5.1地震序列的微震检测及重定位

基于云南地震台网记录到的2017年漾濞M_S4.8和M_S5.1地震序列的地震目录和波形数据,采用结合波形互相关的双差定位方法得到了M_S5.1主震前7天至主震后30天共258次地震事件更加准确的震源位置,然后利用匹配定位方法获得了更加丰富的地震目录,并对其余震活动的时空演化特征进行了初步分析。双差重定位获得了漾濞地震序列中193次地震的精确位置,水平方向和垂直方向的定位误差均值分别为1.3 km和1.9 km,重定位及震源机制解反演结果显示维西—乔后断裂中南段应为此次地震的发震断层。使用匹配定位方法检测得到了1 625次中小地震事件,事件数量是云南地震台网目录中的6倍,被检测事件主要为主震后早期余震及M_L≤0微震事件。对比2016年云龙M_S5.0地震序列的研究结果后初步分析认为,“震群（多震）型”和“主余型”地震序列的一个主要区别可能是“震群（多震）型”地震序列的余震数量相对更少、余震活动衰减得更快。      

利用剪切波分裂研究四川九寨沟MS7.0地震震源区的地壳应力场

本文讨论了2017年8月8日九寨沟M_S7.0 （M_W6.5）地震之后震源区的地震活动性,并针对震后九寨沟地区架设的4个流动地震台站记录到的地震波形数据进行了剪切波分裂参数的计算,分析了震源区地壳应力场的特征性变化。九寨沟M_S7.0地震的地震活动性表现为:8月8日主震后震源区地震活动明显增强,直至8月22日之后才有所下降,余震震级整体也呈降低趋势,震级多集中在M_S2.0以下。震源区4个流动台站的剪切波分裂结果表明:快波偏振方向近似为东西向,揭示了九寨沟地区的区域应力场方向为近似东西向;区域内慢波延迟时间的平均值为5.04 ms/km,慢波延迟时间在九寨沟M_S7.0地震之后降低,反映了随着余震的发生,地壳应力不断调整,应力持续释放,直至8月底地壳应力场逐渐趋于稳定状态。      

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

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

2022年1月8日青海门源MS6.9地震深部构造背景浅析

2022年1月8日青海境内的托莱山—冷龙岭断裂附近发生了门源M_S6.9地震。结合地壳厚度、速度结构及各向异性等资料探讨了门源地震的深部构造特征,揭示了门源地震的发震位置与地壳结构变化的密切关联。结果显示:门源M_S6.9地震发生在地壳厚度和v_P/v_S值都出现快速空间变化的区域;大约在10—20 km深度范围内,震源位于P波速度从浅到深由高速变低速的垂向过渡区,同时也是S波速度和泊松比分布呈现明显横向变化的过渡区域,震源下方存在明显的低速区;冷龙岭断裂两侧相速度的方位各向异性变化比较明显。1月12日的M_S5.2余震震中紧邻2016年M_S6.4地震震中,揭示出2022年门源M_S6.9地震及其余震活动导致了冷龙岭断裂比较充分的破裂,两次门源地震主震之间及邻区短时间内难以积累更大能量,因而短时间内发生更大地震的可能性不大。青藏高原东北缘的持续向北扩展所导致的地表隆升和地壳增厚是该地区强震频发的主要构造成因。      

九江—瑞昌序列震源位置及震源区速度结构的联合反演

采用震源位置和速度结构的联合反演方法确定2005年九江一瑞昌5.7级地震序列的分布和震源区的速度结构.结果显示九江一瑞昌地震序列存在NW330.方向和NE550方向的两组共轭分布,其中NW330.方向的优势分布更为明显;震源深度分布在3～18km的范围.地震主体破裂发生在上地壳-中地壳过渡层的P波速度低速区内,主震发生后余震向东南与西北方向两头扩展,其中往西北方向破裂时可能遇到了障碍体,积累的能量转向西南-北东方向继续破裂,引发了4.8级强余震.     

重庆垫江MS4.4地震:一个可能发生在浅地表的地震

为更好地认识2016年8月11日重庆垫江M_S4.4地震的发震成因,本文首先分析了近场地震台站记录到的该地震产生的短周期面波（Rg）强频散现象,该现象表明本次地震具有明显的浅震特征。然后采用CAP方法反演了该地震的震源机制解和震源深度,结果显示本次地震是走向为近EW向的逆冲型地震,矩震级为M_W4.1,最佳震源深度为1—2 km。在此基础上,结合现场震害、烈度等震线长轴方位和地质剖面图等资料,分析认为在区域应力集中加载作用下,震源区沉积盖层内部的软弱地层产生挤压变形,导致地层内斜交的裂隙发生切层失稳滑动,从而诱发了垫江M_S4.4地震,同时排除了震源区NE走向的黄钦垭口地表断层为发震构造的可能。      

2018年5月28日吉林松原MS5.7地震发震构造分析

利用黑龙江、吉林、内蒙古三省数字地震台网的三分量宽频带波形资料对2018年5月28日吉林省松原市宁江区发生的M_S5.7地震进行了全波形矩张量反演,获取了此次地震的震源机制解和矩心位置;并使用震源-矩心方法讨论了该地震的发震断层。研究结果显示:松原M_S5.7地震的矩震级为M_W5.2,矩心位置为（45.225°N,124.685°E）,矩心深度为7 km。震源机制解参数显示:该地震为走滑型;节面Ⅰ的走向、倾角和滑动角分别为217°,82°和164°;节面Ⅱ的走向、倾角和滑动角分别为309°,74°和8°;双力偶成分占96.4%,方差减少为93%。震源?矩心图显示震源更接近节面Ⅰ,与北东向的扶余—肇东断裂走向及倾角一致,因此,推测扶余—肇东断裂为发震断层。      

新丰江水库中上地壳P波三维速度结构特征及库水的渗透影响

利用震源位置和速度结构的联合反演得到2007年6月—2014年7月新丰江水库地区地震序列的震源位置及中上地壳P波三维速度结构模型,并进一步研究库区序列分布及速度结构特征.结果显示:库区中上地壳不同深度P波速度存在显著横向不均匀性,浅部库区速度高于周缘,在5~10 km上地壳从库坝下游白田至库尾锡场NW方向存在高速异常体以及2个低速间断区域,低速间断区分别位于人字石断裂与南山—坳头断裂交汇处以及1962年6.1级地震震源区,库水可能沿低速间断区的人字石断裂、石角—新港—白田断裂下渗至13~14 km的地壳.在10~14 km地壳以NE走向的大坪—岩前断裂为界,NW侧为最高速度6.2 km·s~(-1)的高速区域,SE侧从库区中部回龙至库坝下游白田为显著低速异常区域,是可能的库水渗透影响区域,亦是库区中强地震集中区.库区地震多发生在高速体内部、高低速过渡带或低速的渗水通道两侧.     

由中小地震的震级ML和MS估算1970年通海地震前后其周围地区的剪应力值

本文研究了震源破裂过程及体波和面波的传播特点,导出了三个用中小地震的震级M_L和M_S估算剪应力值τ₀的公式和一个估算震源尺度2a的公式: logτ₀=1.5M_L—M_S—0.68 M_L=2logτ₀+2log（2a）—0.06 M_S=2logτ₀+3log（2a）—0.77 log（2a）=M_S—M_L+0.71用这些公式可以比较容易地估算出τ₀值和2a值。 通过上述公式估算的τ₀值,研究了通海大震前后应力场的演变过程:大震前,通海附近的峨山显示出很高的应力值,达282巴。而其他地区0值都在100巴以下。大震发生后,破裂带中部0值都下降到100巴左右。在两端还有少数高应力值点,当几个5级多强余震发生后,应力值都降到100巴左右。      

四川芦山MS7.0地震震源区及其周边区域P波三维速度结构研究

利用四川数字地震台网和流动地震台站在芦山M_S7.0地震震后(2013年4月20日—6月23日)记录到的2026次区域地震事件的28188条P波到时资料,采用地震层析成像方法反演得到了芦山地震震源区及其周边区域中上地壳P波三维速度结构. 结果表明,浅部地壳的P波速度异常分布特征与地表地质构造、 地形和岩性密切相关,即成都断陷盆地表现出与第四纪沉积有关的低速异常区;犍为、 乐山一带的川中微升区和川青块体龙门山以西的邻近地带均表现为与构造抬升有关的高速异常;宝兴、 康定附近分布的基性火山岩及火山碎屑岩均呈局部高速异常分布. 芦山地震震源位于高低速异常分界线附近且偏向高速体一侧,其下方存在明显的低速异常分布,可能与流体的存在有关. 流体的作用导致中上地壳内部发震层的弱化,使孕震断层易于破裂,可能对芦山地震起到了触发作用. 芦山地震与汶川地震两次地震的余震密集区相距50 km,这50 km地震空区震源体的深度范围附近目前正处于高速异常区内,加之龙门山断裂带西南段又具有比较典型的断错地貌发育,使得该段地震空区(大邑—邛崃活动断裂破裂空段)现在所处的深浅部构造环境变得复杂,其潜在的地震危险性仍值得进一步关注.      

利用远震P波层析成像研究巴颜喀拉地块东缘及邻区的深部结构

为了揭示巴颜喀拉地块东缘及邻区的壳幔速度结构差异,获取2017年九寨沟M_S7.0地震的深部构造背景,本文收集了2009年5月至2016年8月期间四川及邻区数字测震台网的203个地震台站所记录到的远震P波走时数据,应用有限频体波走时层析成像方法,反演得到了巴颜喀拉地块东缘及邻区50—600 km深度范围内的三维壳幔P波速度结构。反演结果表明:巴颜喀拉地块东缘及邻区的壳幔速度结构具有明显的横向不均匀性和分区特征,松潘—甘孜地槽褶皱系、西秦岭和祁连山褶皱系的整体速度异常较低,研究区东部具有克拉通性质的四川盆地西北缘和鄂尔多斯地块南缘则呈明显的高速异常。上地幔P波速度结构特征差异表明松潘—甘孜地块的抬升可能与地幔上涌有关,巴颜喀拉地块东缘九寨沟震区及周边50—250 km深度范围内的上地幔存在低速异常,在400—600 km地幔过渡带深度范围内表现为明显的高速异常特征。巴颜喀拉地块向东南方向运移受到东部高速、高强度的扬子克拉通地块对青藏高原物质东向挤出的强烈阻挡,而九寨沟震区处于松潘—甘孜地块重要的北东边界断裂交会处附近,应力容易在此集中,这些因素均可能是东昆仑断裂塔藏段与岷江断裂北段交会处附近发生九寨沟M_S7.0地震的深部动力学背景。      

2014年新疆于田MS7.3地震序列重定位及其发震构造初步研究

2014年2月12日在新疆于田县发生了M_S7.3地震,主震前一天在震区发生了M_S5.4前震,震后余震活动频繁,由于震区台站十分稀疏和不均匀、地壳速度结构复杂,台网常规定位结果精度有限,很难从中获得序列的空间分布特征和活动趋势的正确认识.本文首先利用位于震区附近的于田地震台5年记录的远震波形数据,采用接收函数方法研究了震区附近的地壳结构,建立了震源区的地壳速度模型.在此基础上,联合震相到时和方位角对2014年于田M_S7.3地震序列(从2014年02月11日-2014年04月30日,共计577次地震)进行了重新绝对定位.结果显示,(1) 重定位后的前震和主震震中位置明显向地表破裂带及其附近的阿尔金分支断裂(南肖尔库勒断裂和阿什库勒-肖尔库勒断裂)靠近,两者相距5.4 km,主震位置为36.076°N、82.576°E,震源深度为22 km, 前震位置为36.055°N、82.522°E,震源深度为19 km;(2) 本文重定位结果显示,余震序列沿NEE-SWW展布,优势分布长度约73 km、宽度约16 km,平均震源深度为14.8 km,其中77%的余震分布在地表破裂带的西南端,这部分余震中少数沿阿什库勒-肖尔库勒断裂分布,绝大多数沿北东东向的南肖尔库勒断裂分布,位于地表破裂带东北端的余震沿阿什库勒-肖尔库勒断裂分布,但发生在地表破裂带的余震极少;重定位后,位于地表破裂带西南侧的震中分布由台网目录的近南北向变为北东向,与地表破裂带、南肖尔库勒断裂和阿什库勒-肖尔库勒断裂走向一致;(3) 沿重定位剖面的地震分布,可推断位于地表破裂带西南段的南肖尔库勒断裂与位于北东段的阿什库勒-肖尔库勒断裂倾向反向,南肖尔库勒断裂的倾向为SE,阿什库勒-肖尔库勒断裂的倾向为NW,这与本次地震野外考察得到的断裂性质一致.综合重定位结果、地表破裂带分布、震源机制解、南肖尔库勒断裂和阿什库勒-肖尔库勒断裂的性质认为,2014年于田M_S7.3地震的发震构造为阿尔金断裂西南尾段的两条分支断裂——南肖尔库勒断裂和阿什库勒-肖尔库勒断裂.     

基于GRACE数据的尼泊尔MS8.1地震北向重力梯度变化

利用德克萨斯大学空间研究中心发布的GRACE RL05月重力场模型数据,采用300 km扇形滤波,得到2015年尼泊尔M_S8.1地震北向重力梯度的时空分布,之后利用最小二乘拟合方法分析尼泊尔及其邻近区域北向重力梯度的长期变化趋势及研究区内6个特征点的北向重力梯度月变化时间序列,并结合黏弹性分层位错理论对GRACE检测尼泊尔M_S8.1地震北向重力梯度同震变化的可能性进行分析. 研究结果表明:在尼泊尔M_S8.1地震发生前北向重力梯度表现出比较明显的正负异常变化,而该地震即发生在北向重力梯度正负变化的零值线附近;研究区北向重力梯度年变率在印度板块边界及其垂直方向所形成的四象限呈正负相间分布,6个特征点的北向重力梯度在2012年4月出现比较明显的跳变. 由此推断,北向重力梯度的动态变化过程反映了震前区域物质迁移和震后壳幔物质黏滞性调整等问题.      

利用两次新疆于田MS7.3地震对 图像信息方法可预测性的检验

以青藏构造块体为研究区域, 采用以往利用图像信息方法预测M7.0地震得到的较好的计算参数(网格尺度为1平方度, 预测时间窗口为8年), 系统计算了青藏构造块体1993年以来的“地震热点”图像, 特别研究了两次于田M_S7.3地震前后的图像信息演化过程, 并通过ROC方法和R值评分方法检验, 对图像信息方法在该地区的预测效能进行评价. 结果表明: ① 2008年3月21日于田M_S7.3地震和2014年2月12日于田M_S7.3地震均发生在地震热点上, 且发生在平均发震概率随时间变化的峰值点年份. ② 在回溯性检验的时间段内青藏构造块体共发生9次M≥7.0地震, 其中6次地震前震中附近均有稳定的地震热点出现, 包括两次于田M_S7.3地震. 1996年喀喇昆仑山口M_S7.1地震及2008年汶川M_S8.0地震只在个别预测时间窗内出现地震热点, 并没有出现稳定的地震热点演化过程, 而1996年丽江地震发生前震中附近则没有出现地震热点. 该结果与前人研究结果存在差异, 值得进一步研究. ③ ROC检验和R值评分结果均显示图像信息方法在本文采用的模型参数下对大地震具有较好的中长期预测效果. 本文还根据目前存在的地震热点给出了青藏高原块体2014年以后M7.0地震发生概率高的3个可能地区.      

利用P波参数阈值实时估算地震预警潜在破坏区范围

由于传统的潜在破坏区范围估算方法只能在已获取到震中位置和地震事件结束后才能产出,且往往需要数分钟的耗时,其实时性已无法满足地震预警要求。因此,为了快速产出潜在破坏区范围估算结果并将其用于预警,本文采用了一种结合现地预警技术和区域预警技术、基于预警参数（位移幅值P_d和特征周期τ_c）阈值的实时潜在破坏区范围估算方法。首先利用国内地震事件（4.0≤M_S≤8.0）的记录数据和日本强震动观测事件（6.5≤M_J≤8.0）的数据拟合出特定的适应于我国的参数关系式,包括τ_c与震级M的相关性、P_d与峰值速度PGV的相关性以及P_d与τ_c和震源距R的相关性;其次,根据最小震级（M_S6.0）和仪器烈度（Ⅶ度）定义相应的参数阈值（P_d＝0.1 cm和τ_c＝1.1 s）;最后,利用已有的3次破坏性地震事件数据开展线下模拟,对该方法的适应性和时效性进行了验证。结果表明,对于2013年M_S7.0四川芦山和2014年M_S6.5云南鲁甸两次中强地震,震后约10 s即可获取到比较稳定的潜在破坏区范围估计结果;而对于2008年M_S8.0汶川特大地震,在其记录台站分布密度不高的情况下,震后40 s左右的估算结果始呈稳定状态。      

川滇地区地壳密度变化与强震孕育关系

为研究川滇地区地壳密度的时空演化特征与强震孕育的关系,本文基于2011—2014年川滇地区的重力复测资料,利用阻尼最小二乘反演算法,获得了川滇地区0—50 km深度范围内分辨率为55 km×55 km×10 km （长×宽×高）的三维动态密度变化模型。以所获取的动态密度变化为依据,分析了川滇地区三维密度变化特征与2013年四川芦山M_S7.0、2014年云南鲁甸M_S6.5和四川康定M_S6.3地震的关系,并由此对强震重点构造部位的深部地壳结构特征、孕震背景及区域动力学过程进行了深入分析。结果显示:川滇地区出现多个与主要活动断裂带展布方向基本一致的密度变化高梯度带,在三次地震的震中区及其附近观测到明显的区域性密度变化异常。15—35 km深度范围内的密度变化水平剖面显示:强震容易发生在上地壳密度变化正、负异常过渡的高梯度带和密度变化四象限分布的中心;中地壳深度密度变化低异常是强震孕育的主要介质条件;下地壳深度密度变化低异常或密度变化高梯度带均有可能是孕育地震的主要介质结构。0—50 km深度的垂直剖面上的密度变化结果显示,地震震中区及附近浅部、深部地壳呈现解耦变化。壳内垂向正负密度变化过渡带可能是强震孕育的又一个主要特征构造。