2021年5月21日云南漾濞MS6.4地震序列地震矩张量及发震构造

2021年5月21日云南大理州漾濞县发生了MS6.4 地震.我们利用区域地震台网记录的地震波形资料,首先采用多点源地震矩张量反演方法确定了漾濞地震序列中 3 次MS≥5.0地震的矩心矩张量解,分析研究了地震矩释放的最佳模型;然后对序列中较大地震进行了绝对定位,结合余震序列重新定位结果研究了地震矩心在断层面上的位置.结果显示MS5.6 前震可用2 点源模型模拟,矩震级分别为MW 5.3、MW5.1,矩心时间相隔约30 s,矩心位置相距约 2 km.MS 6.4 主震可用单点源模型模拟,矩心与起始破裂点平面距离约 5 km.前震和主震的矩心均位于地表以下 6 km处,矩心与起始破裂点的位置关系显示两地震向南东方向单侧破裂,断层以"前震-主震"型地震序列典型的"撤退式"方式破裂,MS5.6 前震的发生降低了断层面的抗剪强度,从而发生了更大的MS6.4 主震.MS5.2余震可用单点源模型模拟,起始破裂点与矩心空间位置相近,在地表以下约 10 km处.余震区构造应力场反演结果显示漾濞 6.4级地震序列属于区域应力场触发的地震活动,地震序列震源机制解符合走滑断裂伴生的负花状构造系统内部断裂的运动特征,余震的空间分布图像显示花状构造系统内部的两条断裂发生了地震活动.     

1989—1999大同地震序列的隐伏断层研究:库仑应力分析和余震JHD重定位

1989年到1999年,大同—阳高地区发生了一系列MS≥5的中强地震.本文基于前人对1989年三次MS≥5地震的震源机制反演的结果,通过建立不同断层模型,利用库仑应力方法,计算前震对于主震,以及前震和主震对于余震的库仑应力触发关系,提出了一种可能的破裂模型,即1989年前震沿北西西方向发生左旋破裂,之后主震和余震沿北北东方向发生右旋破裂.根据这种破裂模式计算得出,前震发生后,主震震源处的库仑应力增加了约2×105 Pa,余震震源处的库仑应力出现下降;主震发生后,余震处的库仑应力出现回升,最后余震处的库仑应力几乎没有变化.基于大同地震台网的近场观测数据,用JHD(Joint Hypocentral determination)定位方法,对1999年11月1日MS=5.6地震后一个月的余震进行重定位,得到一条走向118°,倾角85°的左旋走滑断层,余震的深度分布在5km至20km范围内,显示该断层是隐伏断层.另外提出对主震震中位置约10km的修正.本文对1989年三次MS≥5地震序列和1999年MS=5.6地震余震空间分布的研究揭示该地区存在两条活跃的共轭隐伏走滑断层(1989年主震的北北东方向和1999年地震的北西西方向),并且推断已知的大王村断裂和团堡断裂是地下这两条共轭的隐伏走滑断层构造/地震活动在地表的响应.     

2014年云南鲁甸M_W6.1地震:一次共轭破裂地震

烈度与余震分布显示2014年云南鲁甸MW6.1(MS6.5)地震的发震构造较复杂.为深入了解鲁甸地震的发震断层与破裂特征,本文考虑了单一断层破裂和共轭断层破裂的情况,对震中距250km范围内的近震资料(宽频带资料和强震资料)和远震体波资料进行了反演,得到了鲁甸地震的破裂过程,探讨了滑动分布与余震分布之间的关系.根据反演得到的滑动分布、震源时间函数和波形拟合,认为鲁甸地震是一次在北西向主压应力与北东向主张应力的统一应力场下发生的两条共轭断层先后破裂的一次复杂地震事件.在破裂开始后0~2s,破裂主要发生在ENE—WSW向(近东西向)的断层上,随后NNW—SSE向(近南北向)断层开始破裂,释放了大部分的地震矩.由于近南北向断层南段(即震中以南)的破裂规模较大,且以左旋走滑为主,对近东西向断层的西段起到了一定程度的解锁作用,可能是震中以西无明显主震破裂但存在密集余震分布的主要原因.     

基于区域地震波形的2017年新疆精河MS6.6地震破裂方向性及发震构造研究

2017年8月9日的新疆精河MS6.6地震是近年来天山北缘发生的最大地震,震中位于由多条逆冲断层组成的库松木契克断裂带内.由于震源较深、构造形变复杂、区域地震台站相对稀疏,仅根据震源机制解、余震分布和InSAR观测结果等难以直接判定发震构造.本文针对倾滑型地震发展了一种基于区域地震波形的破裂方向性测定方法,利用余震作为参考地震进行路径校正,根据主震和参考地震的波形时移差和Pn-Pg到时差分别确定主震在水平方向和深度方向的破裂尺度,进而推断同震破裂的延展方向和延伸尺度.本文在反演了主震的点源参数后,应用新发展的方法测定了地震的破裂方向性.点源反演结果显示,精河地震是一个发生在中地壳的高角度逆冲地震,矩震级约6.2,质心深度21km,震源持续时间5.5s,两个双力偶节面分别为102°/45°/106°(NP1)和259°/47°/74°(NP2).破裂方向性分析结果显示,地震的破裂面为南倾的NP1节面,地震沿着破裂起始点向西南方向、向下破裂,总破裂长度约11.5km,其中,沿深度的破裂范围约7km,沿水平的破裂范围约9km,平均破裂速度约2.1km·s-1.综合区域地质资料、卫星影像等判定本次地震的发震断层为精河南断层,地震可能只破裂了断层的下段(17~25km),并未破出地表.     

2021年5月21日云南漾濞6.4级地震序列重定位及震源机制研究

2021年5月21日 21 时,云南省大理州漾濞县先后发生MS5.6 和MS6.4 地震,两次地震震中位置相距约7 km,均位于滇西地区,该地区地处青藏高原东南缘、南北地震带南段,地质构造复杂.地震序列跟踪结果显示漾濞MS6.4地震类型为前震-主震-余震型,MS5.6 地震为MS6.4地震的前震.本文基于云南地震台网的震相报告,采用双差定位方法对漾濞MS6.4地震早期序列(2021年5月18日至 25日,ML0.0 以上)进行重定位,同时,利用 Cut-And-Paste(CAP)震源机制波形反演方法,获得了序列中截止至 5月25日 31 次MS≥3.0 地震的震源机制解和矩心深度,对该序列的发震构造进行了初步分析.结果显示:(1)重定位后获得的 2159 个地震事件呈 NW-SE向展布,长约 25 km,宽约 5～10 km.MS 6.4主震的震源深度为 8.9 km,序列震源深度主要集中在 4～10 km,深度均值约7.5 km.(2)前震序列具有从中间开始破裂,然后向北西向破裂,继而向东南向破裂的特征,漾濞MS 6.4 主震位于余震区的北西端,最大余震MS 5.2 地震位于余震区的东南端.(3)CAP波形反演获得的 31 次MS≥3.0地震的震源矩心深度在4～11 km范围,深度均值约 6.5 km,与重定位结果接近,仅相差 1.0 km,说明重定位的震源深度分布是合理可靠的.(4)震源机制多为右旋走滑型,部分地震具有较为明显的正断分量,反演的区域应力场与目前已知的构造应力场水平主压应力方向一致,反映区域内构造活动主要受区域构造应力场控制.根据重定位后的序列空间分布、震源机制解及震源区构造特征,综合分析认为,此次漾濞地震序列的发震构造可能是维西—乔后—巍山断裂中段附近的次级断层.     

2016年日本熊本地震破裂时空过程联合反演

为了深入认识2016年4月15日日本熊本地震破裂的复杂性,利用远场体波资料和同震InSAR资料联合反演了此次地震的震源破裂时空过程. 联合反演结果表明:熊本地震的震源破裂持续时间约为25 s,整个破裂过程释放的总标量矩为6.03×1019 N·m,对应于矩震级M_W7.1;同震滑动主要集中分布于浅部,破裂以右旋走滑为主,但在沿倾向0—5 km范围内,破裂呈较强的正断特征;此次地震破裂的最大同震滑动量约为4.9 m,且最大同震位错区位于背离断层走向上、距离起始破裂点约5—10 km的区域;破裂前期（0—7 s）,在倾向上向浅表发生破裂,在走向上向东北和西南两侧扩展;大约7 s后,破裂背离断层走向主要向东北方向扩展. 根据有限断层联合反演结果推测,此次熊本地震破裂可能出露至地表.      

玉树MS7.1级地震部分余震重新定位及发震构造分析

综合利用玉树震区应急流动台站观测数据和青海地震台网固定台站观测数据,依据最新的人工地震宽角反射/折射剖面的速度模型,采用Hypo2000地震定位法,对2010年4月18日至4月29日期间玉树震区发生的部分余震进行了重新定位.重新定位后,震源位置的水平和垂直方向平均误差分别为1.35 km和4.68 km,走时残差为0.49 s.震源深度分布范围为1.48～19.85 km,平均震源深度为10.28 km.定位研究结果表明:玉树地震余震沿北西-南东向的甘孜-玉树断裂带的北支,即玉树-隆宝断裂分布,长约97 km.余震分布特征在主震(微观震中)两侧存在差异,可能反映了两侧构造特征存在差异.截止到4月29日,主震东南仍是应力的主要释放区域,余震强度大且活动密集的区域位于主震东南距主震约5 km、横向范围约20 km.主震破裂区的大部分应力在主震过程中得以释放,主震时应力未释放的区域成为主要的余震分布区.余震的连续发生可能已造成主震破裂区相互连通,且破裂范围向西北方向扩展.玉树主震及余震的发震构造为甘孜-玉树断裂的北支,即玉树-隆宝断裂段,断层性质为北东倾向的高角度左旋走滑断层.发震断层的倾角和宽度在帮洞两侧有所不同,帮洞以东发震断层宽度约为12 km,倾角约为83°;而帮洞以西发震断层宽度约为6.5 km,断层倾角约减缓为63°.     

1997年中国西藏玛尼Mem>Ssub>7.9地震的时空破裂过程l

用中国数字地震台网（CDSN）的长周期波形资料反演了1997年11月8日中国西藏玛尼地区MS7.9地震的地震矩张量,用频率域里反褶积方法从P波和S波震相中分别提取了震源时间函数,并经反演依赖于方位的震源时间函数获取了断层面上破裂随时空变化的图象．矩张量反演结果表明:玛尼地震发震应力场的P轴和T轴均接近于水平,P轴在NNE方向（方位角29,倾角7）,T轴在SEE方向(方位角122,倾角23),断层错动以走滑为主;标量地震矩为3.41020 Nm,矩震级MW=7.6．由矩张量反演得到的震源时间函数显示,这次地震是由一次较小事件和较大事件组成的,较小事件大约持续5 s,较大事件持续约10 s．由余震分布可推断出玛尼地震的发震断层是走向为250、以走滑为主的左旋-逆断层,断层面的倾角比较陡,约88．根据反演结果计算了理论格林函数,然后用反褶积方法提取了震源时间函数．从不同台站的P波和S波中分别提取的震源时间函数一致表明这次地震破裂的时间历史比较简单,可用一宽度约10 s的正弦形的函数近似表示．进一步反演从不同台站上得到的、依赖于方位的P波和S波震源时间函数,获得了断层面上滑动的时空分布图象．从破裂的记忆式快照看,破裂开始于断层的西端,然后向东向下发展,总体上具有单侧破裂的特征．破裂面由3个破裂子区构成．一个在断层西端,深度约10 km(西区);另一个距断层西端约55 km,深度约35 km(东区);第3个距断层西端约30 km,深度约40 km(中区)．３个破裂子区构成约长70 km,宽60 km的破裂面．从破裂的遗忘式快照看,这次地震的破裂过程是相当复杂的,在不同时刻断层面上发生错动的地点并不相同,显示出这次地震的破裂过程具有愈合脉冲的特征,而且在断层面上的某些部位发生了多次错动;另一特征是最先和最后破裂的部位都不是主要的破裂区．根据标量地震矩计算了断层面上静态位错的分布,位错最高的3处（西区、东区和中区）的位错值分别为956 cm,743 cm和1 060 cm．由断层面上位错的分布推知,破裂主要集中在震中以东长约70 km的断层上;从余震的分布看,震中以西余震稀疏而震中以东余震密集．这些都表明这次玛尼MS7.9地震是北东东-南西西向至近东-西向断层向东扩展的结果．      

2019年6月17日四川长宁地震重定位及震源机制研究

采用双差重定位和W震相波形反演方法分析 “地震编目系统” 和中国地震台网中心提供的地震观测报告及区域地震波形数据,对2019年四川长宁地震序列进行了重定位,反演获取了M＞4.5地震的震源机制解。地震序列重定位结果显示,长宁地震序列沿NW优势方向呈条带状分布,集中分布于5—10 km深度范围,且发震断层面呈高倾角。震源机制反演结果表明,2019年6月17日四川长宁M_S6.0主震的两个可能发震断层面参数分别为:节面Ⅰ走向12°,倾角50°,滑动角139°;节面Ⅱ走向131°,倾角59°,滑动角48°,最优矩心深度为7.5 km,矩震级M_W5.74。此外几个M＞4.5余震的震源机制也基本与主震类似,均为以逆断为主外加少量走滑的地震破裂事件。综合分析长宁地震序列的重定位、震源机制反演结果以及震中和附近区域的地质构造背景信息推断,本次长宁主震的发震破裂面呈NW?SE走向,发震断层为长宁—双河背斜东北翼发育的逆冲断层。      

芦山7.0级地震序列的震源位置与震源机制解特征

基于中国国家和四川区域数字地震台网记录,采用HypoDD方法精确定位了四川芦山M_L2.0级以上地震序列的震源位置,采用CAP方法反演了36次M_L4.0级以上地震的最佳双力偶震源机制解,并利用小震分布和区域应力场拟合了可能存在的发震断层面参数,从而综合分析了芦山地震序列的震源深度、震源机制和震源破裂面特征,探讨可能的发震构造.结果显示,7.0级主震的震源位置为30.30°N、102.97°E,初始破裂深度为15 km左右,震源矩心深度为14 km左右,最佳双力偶震源机制解的两组节面分别为走向209°/倾角46°/滑动角94°和走向23°/倾角44°/滑动角86°,可视为纯逆冲型地震破裂,绝大多数M_L4.0级以上余震的震源机制也表现出与主震类似的逆冲破裂特征.M_L2.0级以上余震序列发生在主震两侧,集中分布的长轴为30 km左右,震源深度主要集中在5～27 km,M_L3.5级以上较大余震则集中分布在9～25 km的深度上,并揭示出发震断层倾向北西的特征.利用小震分布和区域应力场拟合得到发震断层参数为走向207°/倾角50°/滑动角92°,绝大多数余震发生在断层面附近10 km左右的区域.综合地震序列分布特征、主震震源深度和已有破裂过程研究结果,可以推测主震破裂过程自初始点沿断层的两侧扩展破裂,南侧破裂比北侧稍长,滑动量主要集中在初始破裂点附近,可能没有破裂到地表.综合本文研究成果、地震烈度分布和现有的科学考察结果,初步推测发震构造为龙门山山前断裂,也不排除主震震中东侧还存在一条未知的基底断裂发震的可能性.     

日本熊本MW7.0地震同震位移场和震源滑动模型反演

利用2016年4月16日日本熊本M_W7.0地震震中周围94个近场强震动台的观测资料和新近改进的强震经验基线校正方法SMBLOC,尝试解算并绘制了一个内陆M7左右走滑型地震的同震位移场全貌,并反演了其震源滑动模型.与日本国土地理院(GSI)公布的该地震57个GPS同震位移结果的比较显示,两种完全不同资料、不同解算方法给出的水平同震位移场的最大幅值均为100 cm左右,均呈右旋走滑为主兼具部分正断分量的震源机制.强震最大水平和垂直永久位移分别为104.5 cm和58.0 cm, 分别出现在震中东北侧的KMMH162台和KMM005台.两种资料单独以及联合反演的震源滑动模型均表明,此次地震为北东侧破裂为主并呈双事件特征,且主要滑动均不在初始破裂点附近, 而是集中于第二次事件周围,即距离初始破裂点东北侧约20 km处的走向长约40 km、倾向宽约20 km的范围内.基于强震和GPS模型所得的最大滑动量分别为5.10 m和5.87 m,量级一致,反演矩震级均为M_W7.1左右;主破裂区近地表滑动量比野外调查结果略微偏大,可能与数值效应有关.此外,还利用不同方法得到的解算结果比较了熊本地震特有的12组台间距在3 km以内的GPS-强震台站对各自的三分量同震位移,其结果表明对于M7左右的地震而言,SMBLOC方法解算同震位移时方向和幅值的可靠性下限约为2 cm.      

2010年4月14日青海玉树地震破裂过程快速反演

基于远震资料的破裂过程反演方法,通过反演2010年4月14日玉树地震的全球宽频带地震垂直向P波波形记录,先后4次反演得到了玉树地震的破裂过程,并比较和讨论了这4次结果.结果表明,玉树地震的破裂过程具有如下基本特征:①地震主要由两次子事件组成,分别对应于震中附近以及震中东南方向上的两块滑动量集中的破裂区域,其中与第2次子事件对应的震中东南方向上滑动集中的区域破裂贯穿至地表;②最大滑动量和最大滑动速率分别为2.1m和1.1m/s,断层滑动速率较大;③玉树地震总体上是一次单侧破裂事件,破裂从初始破裂点(即仪器测定的震源位置)开始,主要向震中东南方向扩展,由"地震多普勒效应"导致在震中东南方向上产生强烈能量聚焦,是玉树城区遭受严重破坏在震源方面的主要原因.     

基于深度学习震相拾取和密集台阵数据构建青海玛多Ms7.4地震震源区高分辨率地震目录

高分辨率地震目录有助于描绘断层的精细结构和认识发震断裂的构造形态及发震机制.基于玛多地震科考布设的短周期台阵数据,本文利用深度学习自动拾取P/S波震相、震相关联、绝对定位、相对定位等定位流程,构建了玛多Ms7.4地震后第14天至第43天的高分辨率地震目录,揭示了玛多震源区主震西侧以及主震向东20 km区域范围内的地震序...     

2016年11月25日新疆阿克陶6.7级地震序列震源机制特征分析

本研究利用新疆区域数字地震台网的波形资料,采用CAP方法反演了2016年11月25日阿克陶6.7级地震的前震、主震及11次M_S ≥ 3.6余震序列的最佳双力偶震源机制解,得到阿克陶6.7级地震最佳双力偶机制解：节面Ⅰ走向20°/倾角69°/滑动角-10°;节面Ⅱ走向114°/倾角81°/滑动角-159°,表明此次阿克陶6.7级地震为一次走滑型地震事件,结合震源区的地震地质构造及余震序列空间分布等已有研究成果,判定节面Ⅱ代表了主震的发震断层面。主震最大主压力轴方位为339°,与震源区附近历史中强震P轴近NW向的优势方位基本一致。其4.8级前震的震源机制解为走滑型,与主震震源机制解具有较高的一致性。11次余震中有6次为走滑型地震,3次为逆断型地震,1次正断型地震,1次混合型地震,且多数地震具有近NW向的P轴方位。此次6.7级地震序列的震源深度分布于6~16km之间,而大部分地震为9~13km,与本文计算得到的主震的震源深度10km相差不大。此外,初步分析了兴都库什-帕米尔地区强震活动与此次阿克陶6.7级地震的关系。     

2014 年云南盈江Ms 6.1 地震震源机制研究

利用中国数字地震台网记录的区域宽频带波形,通过频率域和时间域多步反演,研究了2014年云南盈江Ms6.1地震基于点源模型的震源机制解和有限断层模型.考虑到使用不同的波形资料类型和简化的一维速度模型等因素对震源参数反演结果的影响,进行了大量的测试比较.结果表明,使用近震波形和本区域简化一维速度模型M1,波形拟合误差最小.基于点源模型的震源机制解显示此次地震发震断层面参数分别为:走向176°/倾角84°/滑动角-173°,表现为一次右旋走滑错动为主的事件.矩心在水平方向上位于震中(24.99°N,97.84.E)北东向约7km,最佳波形拟合矩心深度7km.平均总标量地震矩M0为7.56×1017N·m,计算成矩震级为Mw5.8.进一步模拟高达0.5Hz的高频波形,获得了盈江地震的有限断层模型,结果显示此次地震未表现出明显的破裂方向性.破裂半径约10km,整个破裂面积为267.2km2,平均滑动量约0.05m,破裂在5 s内释放了大多数能量.震后0～2s内,破裂以孕震点为中心向四周同时扩展,在深度7～ 17km内释放了部分能量.2s后,破裂朝断层面顶部和沿走向两侧进一步延伸,约5s后破裂基本停止.     

2011年日本MW9.0地震相关研究综述

2011年3月11日, 日本海沟发生的9级地震造成重大人员伤亡, 受到社会普遍关注, 本文基于此次日本9级地震相关研究结果, 尝试从不同侧面分析此次地震的观测、 现象和认识, 主要包括如下几点: ① 此次地震发生在太平洋板块西北边界上日本海沟俯冲带上, 同震破裂可能存在深浅两个位错集中区, 较深的位错集中区位错量相对较小, 但历史上7级地震多发; 而较浅的位错集中区位错较大, 但历史上强震活动相对较弱; ② 基于GPS观测资料为约束的相关断层运动研究结果表明, 日本海沟断层运动背景以大范围稳定闭锁为主(闭锁区空间尺度与同震破裂尺度相当), 自2003年日本北海道8级地震后日本海沟地区断层运动开始出现扰动, 2008年以后有几次7级左右地震震后余滑分布明显比主震位错量要大, 之后分别于2008年和2011年观测到显著慢滑移事件, 最后分别于2011年3月9日和3月11日发生7级前震和9级主震, 震前日本海沟俯冲带断层运动变化过程比较清楚; ③ 可能是由于监测的原因, 传统上的前兆观测并未出现显著异常, 其震前异常主要为: 部分地震活动参数表明强震震源区震前应力状态相对较高、 区域地表运动速率的短期异常等; ④ 对于震源区物理性质的分析引起了更多的科学问题, 例如, 震源区介质物性是否与周边存在显著差异、 断层摩擦性质是否决定了发震能力和破裂过程、 震前断层运动是否存在预滑、 震前深部流体是否影响到震源区断层运动等。 他山之石可以攻玉, 希望本文对地震预测预报基础研究工作能起到抛砖引玉的作用。     

2008年新疆乌恰Mw6.7地震震源机制与形变特征的InSAR研究

2008年10月5日新疆乌恰M_w6.7级地震发生在南天山、帕米尔高原及塔里木盆地交汇地带,基于地震波反演的震源机制解确定的震源深度存在较大差异.本文利用日本ALOS卫星的PALSAR图像,获得了本次地震的同震形变场,基于卫星视线向（LOS）和方位向（Azimuth）的形变,采用均匀弹性半无限位错模型和有界最小二乘（BVLS）算法,以网格矩形位错元法对发震断层的几何产状、滑移及分布进行了估算,结果表明本次地震以逆断破裂为主,断层面上最大位错量接近3.4 m,形变中心位于73.8040°E,39.5335°N,深度约5 km,震级估算为M_w6.6;地震发生在走向46°,倾角48°的断层上,发震断层长30 km,宽14 km,闭锁深度9 km,符合该地区浅源地震多发的构造特点,发震断层为乌合沙鲁断裂带.InSAR反演的滑移形变主要集中于地下2~7 km,表明乌恰地震为浅源地震,可能与该断层附近历史地震未完全释放的残余应力积累有关.同时,InSAR反演的断层位错分布呈现双破裂特征,震级分别为M_w6.5和M_w6.1,可能与本次地震的主震和余震相对应,也可能是由主震激发而产生的两组破裂.     

基于余震重定位和震源机制解研究青海玛多MS7.4地震序列的发震构造和断裂形态

2021年5月22日02时04分(北京时间),青海果洛州玛多县发生M_S7.4地震,震后余震不断.地质调查和卫星观测对地表断裂痕迹有较好的约束.然而,对于理解区域应力场、地震的产生、传播和终止具有重要意义的地下断层几何结构的约束精度略显不足.利用国家地震台网的连续波形记录,本研究首先基于双差定位法对玛多地震震后25天的余震序列进行重定位,结果显示余震序列大致沿NWW向的江错断裂呈线性分布,位于主震震中两侧,延伸总长~170 km.主震东南侧存在一余震稀疏区,在断裂带东西两端余震分布转向且出现分叉现象,反映出发震断层的复杂几何形态,这与前人研究结果基本一致.进一步采用波形反演方法和P波初动极性反演方法,获得了玛多震源区132个中小余震的震源机制解与震源矩心深度,并基于此对该主余型地震的发震构造与断裂形态进行了初步分析.震源机制解结果表明,玛多M_S7.4主震的发震断裂主要为左旋走滑性质,余震与主震性质整体相同,在断裂带东段存在部分逆冲型余震.震源机制解约束的区域主应力方向约N60°E,与区内整体走滑断裂作用相一致.余震震源深度略微起伏,主要集中在10~12 km,且浅部余震较少,表明浅部应力可能主要通过主震释放,余震深度分布可能限定了主震同震破裂的下边界.玛多主震破裂起始于断裂带走向和倾向发生明显变化的位置,表明断裂带的复杂几何结构可能是此次玛多M_S7.4地震初始破裂空间分布的决定因素.主震破裂结束的两端都有"马尾状"构造(或次级断层),表明这种分叉断层复杂的几何形态可能控制着主震破裂的最终位置.     

2021年5月21日云南漾濞Ms6.4地震序列重新定位与发震构造分析

本文利用云南区域地震台网提供的震相报告,采用"多阶段定位"方法,对2021年5月21日云南漾濞Ms6.4地震序列的震源位置进行了重新定位.选取序列中13次具有高信噪比宽频带波形记录的Ms≥4.0地震事件,应用CAP波形反演方法,获取了这些事件的震源机制解和震源矩心深度.结果显示,漾濞地震序列震中距离NE侧的维西—乔后断裂约3～10 km,余震区长轴呈NW-SE向展布,长约22 km,震源优势分布深度为4～8 km,反映本次地震发生在上地壳浅部.地震分布北西窄、南东宽的分段特征明显,北西段地震更浅、更密集,南东段事件更深、更分散且存在向S偏转的现象.余震区不同位置的短轴剖面分析表明,发震构造倾向SW,南东段倾角缓于北西段.北西段还存在规模较小的NE向共轭断裂,长度约5 km.5月18日开始的前震自余震区中段向NW逐级迁移扩展,前震序列密集且存在相互触发,符合前震发生的级联模式.主震位于余震密集区的NW端,表明其具备SE向单侧破裂特征.多数4级以上地震的震源机制解都具有与序列长轴走向一致的NW-SE走向、高倾角SW倾向节面;基于震源机制解反演获得的区域构造应力场显示,发震构造受NNW-SSE向近水平主压应力作用发生右旋走滑运动,揭示主发震断层产状和错动类型与维西—乔后断裂基本一致.因此,可初步判定发震构造为维西—乔后断裂的平行伴生断裂,这一构造的形成可能与川滇块体SE向滑移和滇西南块体的顺时针旋转有关.     

基于视震源时间函数的震源过程复杂性分析新方法研究

在全面学习和了解震源破裂过程研究历史的基础上,从地震位移表示定理出发,梳理推导了地震震源时空破裂过程反演研究的基本原理、总结了地震破裂过程反演研究中惯用的约束条件、典型的视震源时间函数提取方法和已有的多种震源破裂过程反演方法,并在已有研究的基础上,尝试性地提出了两种新方法:等矩反投影方法和子事件成像方法。首先,为了提高区域(中强)地震震源过程的时空分辨能力,借助于2013年芦山M_W6.6地震和2014年康定M_W6.3以及M_W5.8地震等实际震例,我们探讨了利用经验格林函数技术从区域地震记录中提取视震源时间函数、进而利用不同视震源时间函数分析和反演中强地震震源过程的效果。其次,为了快速获取震源破裂过程信息,用于震后快速响应以及早期预警,在充分吸收了传统高频反投影方法的优点的基础上,提出了等矩反投影方法,用于对震源破裂过程激发的低频能量源的直接成像。通过多组数值实验测试了新发展的等矩反投影方法的可行性,并将该方法应用于2014年11月在我国四川省康定发生的两次(M_W6.3和M_W5.8)中强地震的震源破裂过程研究以检验新方法的实用性。结果表明,等矩反投影方法对于提取低频能量迁移的主要特征是有效的,是对高频能量反投影方法的必要补充。然后,基于等矩反投影方法的基本思想,并借鉴于传统的震源破裂过程反演方法,新发展了用于对震源破裂时间过程和同震滑动量分布进行直接成像的子事件成像方法。相较于震源破裂过程研究中传统的反演方法而言,新方法不需要事先经验性地给定最大破裂速度以及子断层滑动上升时间,故而在对震源破裂过程进行成像时,允许破裂速度发生变化,并且允许子断层发生多次破裂。同样,利用数值实验测试了新方法的可靠性,且为了检验该方法的实用性,将子事件成像方法用于2014年11月在我们四川省康定县发生的两次(M_W6.3和M_W5.8)中强地震的震源破裂过程研究。结果表明,子事件成像方法能够获取震源过程的主要破裂信息,但与视震源时间函数反演得到的结果相比,破裂尺度略小。最后,将多种方法应用于2015年4月发生的尼泊尔M_W7.8地震及其强余震,获得了这些地震的震源破裂过程图像,对比讨论了这些方法在揭示震源时空破裂过程复杂性信息方面各自的优缺点。具体讲,首先利用远场P波视震源时间函数分析了主震及其强余震的破裂方向,并利用P波波形资料和InSAR资料联合反演了这些事件的震源破裂过程;然后利用等矩反投影方法和子事件成像方法分析了这些事件的低频能量的时空迁移过程和震源破裂过程;最后对利用不同方法得到的结果进行了对比讨论。对于M_W7.8主震而言,联合反演结果表明其以单侧破裂为主,向南东东方向拓展约150km,并沿断层倾向方向向深部拓展约70km,这与利用等矩反投影方法和子事件成像方法直接成像结果得到的主要破裂特征较为一致,但破裂空间拓展尺度要略大,这可能是由于联合反演所给定的3km/s的最大破裂速度要大于直接成像获取的主震破裂速度所致;对于M_W6.7强余震而言,直接波形反演方法与等矩反投影方法和子事件成像方法得到的结果一致表明,该事件破裂主要沿断层倾向方向向深部拓展,并且利用直接波形反演方法搜索到的此次事件最大破裂速度为4km/s,与等矩反投影方法直接成像结果中获取的3.9km/s的破裂速度较为一致,均表明此次强余震可能是超剪破裂事件;对于M_W7.3强余震而言,由于融合了同震InSAR资料的缘故,联合反演结果表明,此次事件以朝向南—东南方向破裂为主,这与基于等矩反投影方法和子事件成像方法的直接成像结果中破裂以朝向南南西—南方向为主的主要破裂特征略有差别,并且直接成像结果表明此次事件破裂速度较慢。概而言之,本论文总结了震源破裂过程研究的发展历史,梳理推导了已有的方法;在此基础上提出了两种新方法,并通过必要的数值实验和对实际震例的应用,阐述了新方法的实用效果和优缺点。