2013年4月20日四川芦山7.0级地震震源破裂特征

2013年4月20日四川芦山发生7.0级地震.本文运用反投影远震P波的方法研究了中心频率为1 Hz的芦山地震震源破裂特征.研究结果显示2013年芦山7.0级地震破裂长度约为20 km,震源破裂时间约为26 s.本文认为在芦山地震的开始阶段(0~4 s)震源的破裂是向震中位置两侧进行的.而芦山地震破裂的第二阶段(5~26 s)是单侧破裂.芦山地震最大能量释放区域位于震中以北.本文对比了运用相同方法研究的发生在同一断裂带上的2008年汶川地震震源破裂特征.发现2013年芦山地震和2008年汶川地震有三点相似之处,即,破裂主要沿北东走向的龙门山断裂带发展;最大能量释放区域没有位于震中;能量都是通过多次子事件来释放的,且第二次能量释放是最大能量释放.对比两次地震破裂区域,可以看出芦山地震的破裂区域是在2008年汶川地震破裂区域的西南端发展的.两次地震的破裂区域占了整个龙门山断裂带的三分之二.     

利用多台阵压缩传感方法反演尼泊尔M_W7.9地震破裂过程

2015年4月25日,尼泊尔地区发生MW7.9地震,震中位于28.1°N,84.7°E.为了详细地研究此次破坏性极强的地震的破裂过程,本文利用多台阵压缩传感方法,使用了阿拉斯加、欧洲和澳大利亚三个台网的共计179个台站的远场P波垂直分量的数据来反演,结果表明本次地震的破裂过程是一个清晰的南东东方向的单侧破裂,破裂尺度约为105km,整体持续时间约为58s.在破裂初始的前15s,能量辐射基本围绕在震源附近,16s后破裂开始向南东东方向以1.9km·s-1的速度破裂.释放能量最大的时间为第38s,位于距震中70km处.该位置从第29秒开始破裂,并持续释放能量长达30s之久.     

2022年9月5日四川泸定MW6.6地震破裂过程研究

2022年9月5日四川省甘孜州泸定县发生M_W6.6地震,利用国家烈度速报与预警工程项目建成的基本站强震动数据,使用迭代反褶积和叠加法(IDS)进行破裂过程反演。反演所得破裂模型显示,破裂面呈NNW—SSE走向,破裂持续时间为15 s,分为4个阶段：首个阶段发生在震后3 s,破裂朝着断层面上倾方向以及SE侧传播;第二阶段为震后6—9 s,破裂继续向SE侧传播并在震中SE侧10 km处迅速加剧,此时破裂滑动速率达到峰值;第三阶段在震后9—12 s,破裂能量继续在SE侧释放,破裂滑动速率逐渐减小,破裂静态滑动累积量达到峰值并趋于稳定;第四阶段在震后12—15 s,破裂能量基本释放完毕,破裂结束。整个破裂由震中向SE方向延伸,由深部向浅部扩展。最大破裂点位于震中SE向10 km附近地下5 km处,最大滑动量为0.8 m,破裂可能出露地表。     

利用阿拉斯加台阵资料分析2008年汶川大地震的破裂过程

2008年5月12日,四川省汶川县发生M_S8.0地震.我们利用美国阿拉斯加区域台网的部分宽频带地震台构成广义台阵, 应用非平面波台阵技术——迁移叠加方法,获得了这次地震的高频（>0.1 Hz）能量辐射源随时间和空间变化的图像.图像表明,这次地震破裂从震中开始向北东方向扩展约300 km,震源过程至少长达90 s,平均破裂速度为3.4 km/s;整个过程可分为两大阶段,前段持续时间50 s,破裂长度约110 km,破裂传播的平均速度为2.2 km/s,后段持续时间约40 s,长度约190 km,平均破裂速度为4.8 km/s.这意味着地震过程的后期似乎发生了超S波破裂,而且后段很可能为前段动态触发所致.     

联合高频GNSS和强震数据快速反演四川泸定6.8级地震断层破裂特征

北京时间2022年9月5日12时52分,四川甘孜藏族自治州泸定县发生6.8级地震.利用震中附近1 Hz高频GNSS观测数据获取了同震速度和位移波形,并快速测定了泸定地震的震中和震级.实验结果表明：高频GNSS反演的震中与美国地质调查局(USGS)发布的震中相差32 km,与中国地震台网中心发布值相差16 km;高频GNSS反演的震级,与两个机构均仅差0.1个震级单位.针对地震预警、震后快速响应等时效性应用,提出了一种联合高频GNSS和强震数据的线源破裂特征快速反演方法.泸定地震实验结果表明：在震后20 s时可获得稳定的线源模型,破裂长度、方向和破裂模式值分别为33.3 km、151°和0.6,破裂方向与USGS震源机制解断层走向相差14°,反演的断层破裂模式为双侧破裂.提出的地震断层破裂特征快速反演方法可用于地震预警、震后灾害快速评估以及紧急响应,同时可为今后联合高频GNSS和强震数据快速测定地震破裂特征提供参考.     

1976年唐山地震震源动力过程的数值模拟

用新LDDA(Lagrangian DiscontinuousDeformation Analysis)方法模拟了唐山地震断层的破裂、错动和应力释放的整个动力过程．模拟结果表明,唐山地震的震源滑动过程在发震断层上各处不一样．近场位移受断层的曲率影响,断层凹侧的位移大于断层凸侧的位移．滑动过冲现象在震中处最大,并向断层两端衰减．我们发现,唐山地震断层的破裂速度和应力降与断层上的初始剪应力大小有关．唐山发震断层的最大动态、准静态位错量和剪应力降均发生在中间部位,分别是7.1 m、6.2 m和8.1 MPa、5.4 MPa,发震断层的平均准静态位错量和剪应力降分别为4.5 m和3.3 MPa,断层破裂的传播速度从震中向东南和西北方向分别为3.08 km/s和1.18 km/s．      

基于GIS研究南黄海活动断层与地震的关系

基于GIS缓冲区分析和叠加分析的空间分析方法,定量研究南黄海活动断层与地震的关系.通过对震中与10 km和20 km断层缓冲区的叠加分析得出断层缓冲区内地震的数量,从中看出:走向NE～NEE向断裂是控制地震发生的主要活动断裂;通过对断层和震中缓冲区的叠加分析得出对地震发生有影响的断层的数量.     

从全球长周期波形资料反演2001年11月14日昆仑山口地震时空破裂过程

利用从全球数字地震台网记录的资料中选择出的震中距小于90°且震相清晰的20个台站的长周期垂直分量P波震相, 通过反演得到了2001年11月14日昆仑山口地震的震源时空破裂过程. 结果表明, 这次地震由三次子事件构成. 第一次子事件的破裂从震中位置(35.97°N, 90.59°E)开始并向东西两侧扩展, 向西以4.0 km/s的破裂速度扩展了140 km, 向东以 2.2 km/s的破裂速度扩展了80 km, 表现为以自东向西为主的不对称双侧破裂, 形成了约220 km长的断层. 在第一次子事件的破裂开始后大约52 s, 在震中以西约220 km的地方, 第二次子事件的破裂开始. 此时, 第一次事件没有结束, 但已进入愈合阶段. 第二次子事件的破裂向东西两侧扩展, 向西以2.2 km/s的破裂速度扩展了50 km, 向东以5.8 km/s的破裂速度扩展了70 km, 表现为以自西向东为主的不对称双侧破裂, 形成了约120 km长的断层. 在第二次子事件开始后大约 12 s, 第二次子事件的破裂与第一次子事件的破裂在震中以西约140 km处发生了聚合. 在第一次子事件的破裂开始后大约56 s, 在震中以东约220 km的地方, 第三次子事件开始. 此时, 第一次事件仍未结束, 但已进入愈合阶段的尾声. 第三次子事件的破裂向东西两侧扩展, 向西以4.0 km/s的破裂速度扩展了140 km, 向东以3.7 km/s的破裂速度扩展了130 km, 基本上是一次不对称双侧破裂, 形成了约270 km长的断层. 在第三次子事件开始后大约36 s, 第三次子事件的破裂与第一次子事件的破裂在震中以东约80 km处发生了聚合. 此后, 震源过程主要是第一次子事件与第三次子事件聚合后的断层运动过程.     

用频率域台阵技术推测2001年昆仑山口西大地震的破裂时间与几何特征

概述频率域远场台阵技术应用于大地震能量辐射源的追踪原理,介绍利用这一技术追踪能量辐射源的一般步骤,并以2001年昆仑山口西大地震为例,分析资料处理中主要参数设置所需要考虑的关键因素. 通过分析埃塞俄比亚/肯尼亚地震台阵（EK台阵）的宽频带波形资料,获得了2001年昆仑山口西大地震的起始破裂点的位置,能量辐射源的空间分布特征,以及破裂持续时间. 根据EK台阵资料的分析结果,2001年昆仑山口西大地震的起始破裂点位于布喀达坂峰东侧（35.92deg;N,91.70deg;E）,破裂持续时间不超过160 s,破裂长度约520 km. 其中,初始破裂点以西为180 km,以东为340 km. 断层西段在布喀达坂峰附近向西南方向发生弯曲,与地表破裂吻合. 断层东段似乎在西大滩断裂附近向东北方向发生弯曲,其方向与西大滩断层的走向一致,但二者相距约30 km. 另外,分析结果似乎表明, 地震断层的西段近乎直立,而东段逐渐南倾.      

基于区域地震波形的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),并未破出地表.     

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后,破裂背离断层走向主要向东北方向扩展. 根据有限断层联合反演结果推测,此次熊本地震破裂可能出露至地表.      

芦山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左右的区域.综合地震序列分布特征、主震震源深度和已有破裂过程研究结果,可以推测主震破裂过程自初始点沿断层的两侧扩展破裂,南侧破裂比北侧稍长,滑动量主要集中在初始破裂点附近,可能没有破裂到地表.综合本文研究成果、地震烈度分布和现有的科学考察结果,初步推测发震构造为龙门山山前断裂,也不排除主震震中东侧还存在一条未知的基底断裂发震的可能性.     

2016年印尼苏门答腊岛海域MW7.8地震震源运动学特征

根据中国和全球地震台网记录的波形记录,采用W震相矩张量反演、反投影分析及有限断层模型反演方法,研究了2016年3月2日印尼7.8级地震破裂过程,分析讨论印尼地震震源运动学特征.结果表明：此地震为一次对称的双侧破裂走滑型事件,北北东─南南西向的断层节面（走向5°/倾角85°）为发震断层面.标量地震矩约6.19×10²⁰ Nm,矩震级为7.79,最大的滑动量约11 m,位于破裂起始点北东,沿着断层走向约30 km处.破裂平均速度2.0~2.2 km·s^-1,破裂持续时间35 s,破裂在5~25 s内释放的能量,约占总能量的97%.最终形成了总长度90 km左右的断层.印尼地震具有破裂持续时间短、破裂速度慢、高滑动能量带相对集中等显著特点.本研究对进一步增进海洋岩石圈地震的震源特性认识有重要参考意义.     

龙门山断裂带北段深部结构与反射地震特征

2008年5月12日汶川M_W7.9特大地震发生在龙门山断裂带,龙门山断裂带深部结构的复杂性制约了地震的破裂过程.通过对研究区区域地质、汶川地震前后采集的地震反射剖面等研究,在对龙门山北段汶川地震断裂带的深部结构和反射地震特征进行了分析的基础上,探讨了它对地表破裂过程的制约.研究结果表明,在地震剖面上,断裂带表现为能量破碎、联系性差;频率剖面上显示整体剖面频率在5~45 Hz,断裂带呈现频率低（15~26 Hz）等特征.龙门山北段映秀-北川断裂在10 km以上是一条倾向北西的高角度走滑兼逆冲性质的断裂,倾角50°~70°.它分割了西侧的轿子顶杂岩和东侧的唐王寨推覆体,错断了早期形成的逆冲岩片,从南到北总位移量由大变小.它高角度的几何形态约束了断裂以走滑为主兼逆冲分量的运动性质,降低了地表滑移量,影响了地震破裂过程以及余震沿断裂带两侧分布的特性.     

远震体波有限破裂过程反演敏感性测试及其在汶川地震上的应用

本文利用远震P和SH波反演得到2008年5月12日汶川大地震(M_W=7.9)的一系列有限破裂模型。使用一种基于小波变换的模拟退火非线性反演方法, 我们将主断层划分成若干个子断层, 在反演时同时确定每个子断层上的滑移量、 滑动角、 上升时间(rise time)以及平均破裂速度。我们首先根据一个假定的破裂模型生成理论地震图, 将该理论地震数据作为输入进行反演, 对该有限破裂反演方法进行了一系列测试, 以验证反演对断层倾角、 平均破裂速度、 最大破裂深度等参数的敏感性。然后我们采用4个不同倾角的断层面来对汶川地震远震体波记录进行反演。结果表明, 若对只在一个断层面上模拟该地震, 30°倾角是个较为合适的值。反演的结果还表明, 此次地震有两个主要的能量释放区域, 并且主断层面存在倾角变化的可能性。在将来的研究中, 可以结合GPS, InSAR测地学以及强震等数据, 来对强震的破裂过程做更细致的研究。     

2015年尼泊尔地震破裂过程的统一模型

模拟2015年尼泊尔地震（主震M_W7.8及最大余震M_W7.3） GPS/InSAR同震位移、远震体波、高频GPS位移波形和强震加速度记录,构建统一震源模型.统一模型分布特征主要由InSAR观测决定,地震矩释放过程则与P波模型相似,静态与高频GPS观测增加了对破裂时空特征的约束强度;各种比对表明,该模型对各基于单一类型反演模型具有很好的兼容性,棋盘测试展现其具有更优空间分辨率,最小可恢复20 km×20 km尺度的空间特征,压缩了非同震信号或误差导致的零散瑕疵,主、余震破裂具有更好的空间对应关系.主震展布范围为140 km×80 km;4 m以上破裂集中在加德满都以北30 km、深度15 km的狭长区域内,最大滑动量为7.4 m;破裂持续总时长为60 s,破裂速度为3.3 km·s^-1,子断层上升时间在10 s内.M_W7.3余震破裂区域位于主震东侧边缘,滑动量围绕震中扩散,扩展范围为30 km×20 km,最大滑动量约为4.4 m,总破裂持续时间为35 s.本次地震中静态和高频的GPS观测亦具备独立约束主震破裂扩展过程的能力.     

2021年吴忠—灵武ML3.6震群重新定位及震源机制研究

2021年7月18日—8月7日,宁夏吴忠—灵武地区发生M_L3.6显著震群活动。本文利用多阶段定位方法对该震群进行了重新定位,并根据gCAP方法反演了2021年7月20日灵武M_L3.6地震的震源机制及震源矩心深度,采用Snoke方法计算了震群中3次M_L3.0以上地震的震源机制,测定了同一地震多个震源机制的中心解。结果表明,该震群中最大的地震即7月20日02时40分M_L3.6地震的震源机制为节面Ⅰ走向289°,倾角72°,滑动角?22°,节面Ⅱ走向26°,倾角69°,滑动角?161°,震源矩心深度为12 km,初始破裂深度为12.5 km;7月20日03时15分M_L3.2地震的震源机制为节面Ⅰ走向290°,倾角82°,滑动角?2°,节面Ⅱ走向20°,倾角88°,滑动角?172°,初始破裂深度为11.9 km;7月21日04时55分M_L3.1地震的震源机制为节面Ⅰ走向285°,倾角53°,滑动角2°,节面Ⅱ走向194°,倾角88°,滑动角143°,初始破裂深度为11.6 km,这些地震震源机制的主压应力轴主要为NE向。该震群序列的震源深度主要相对集中在7—15 km之间,其中M_L3.0以上地震的震源深度主要介于11—13 km,震源深度剖面显示震群相对集中的区域由深到浅大体呈现近似于陡立的展布。本文进一步研究发现区域应力场在灵武M_L3.6地震震源机制NNE向节面产生的相对剪应力为0.393,而在NWW向节面产生的相对剪应力为0.945。结合地质构造和已有断层资料初步分析认为,若NNE向的崇兴隐伏断裂为灵武M_L3.6地震的发震断层,则表明崇兴断裂可能是一条断裂薄弱带,地震破裂方式主要为右旋走滑;若NWW向的未知隐伏断裂为发震断层,则表明NWW向断裂可能为该地震在区域应力场下的剪应力相对最大释放节面,其破裂方式为左旋走滑。      

反投影全球子台网P波记录研究2010年4月14日玉树地震破裂过程

我们运用反投影远震P波研究了2010年4月14日玉树地震的破裂过程.通过台站响应函数对从IRIS下载的经过波形筛选的206个震中距在30°～95°间的台站进行选取.29个能够得到最好台站响应函数且波形相似度高的台站的记录被用于最终的反投影研究.在研究中我们运用4阶方根叠加法和2 s窗长、1 s滑动的滑动窗对总时长为10...