2013年甘肃岷县-漳县MS6.6地震余震 序列目录完备性研究 基于对单台记录地震事件震中与震级的估计

中强地震余震序列地震目录编目是否完备、 震源参数是否准确,直接影响余震序列特征分析、 震后趋势快速判断和强余震预测等研究结果的科学性和可靠性. 2013年7月22日甘肃岷县-漳县M_S6.6地震余震序列目录中存在较多单台记录地震事件,地震观测报告仅给出其震级,而未给出震中位置. 由于余震波形间的相互交叠干扰,使得余震最大振幅的测量误差较大,造成地震观测报告给出的单台事件震级误差较大. 精确估计单台记录地震事件的震中和震级,能够补充完善现有地震目录,提高地震目录的完备性. 本文对单台记录地震事件震中和震级的估计不仅限于单个台站,而是通过分析区域台网中多个台站的波形记录实现. 首先以余震序列中震级较大、 波形记录信噪比较高的地震波形作为模板,使用波形互相关震相检测技术,检测单台记录的地震事件在多个台站的震相到时. 如果能在4个以上台站检测到震相,则利用测震台网常用的HYPOSAT方法估计其震中位置,并利用多个台站记录波形与模板地震的振幅比估计其震级. 之后计算主震发生后不同时间的最小完备性震级,并通过线性拟合得到最小完备性震级随时间变化的表达式,以分析此地震余震序列的目录完备性. 经过计算共得到253个单台记录地震事件的震级和其中177个事件的震中位置,其震中空间分布范围与余震序列中其它地震分布范围基本一致. 震级复测以及与人工拾取震相到时误差对比表明,该方法所得震相检测和震级估计结果具有较好的可靠性. 主震及最大余震发生后的短时间内,有较多数量单台事件的目录所给出的震级偏低,分析认为可能受主震与较大余震后续震相以及余震间相互干扰所致. 主震发生0.02—0.3天内,其余震序列最小完备性震级随时间的对数呈线性下降,在0.3天后最小完备性震级稳定在M_L1.1左右.      

九寨沟M7.0级地震余震目录完备性研究

2017年四川省九寨沟县M 7.0级地震发生后,截止到8月24日共16天内时间内地震台网中心记录到5994个九寨沟地震余震事件,为完备该地震的余震目录,本文利用模板匹配方法对本次地震进行了遗漏地震检测研究,选取震中附近台站记录波形为计算数据,重新检测主震后16天时间内的遗漏地震事件,共得到台网目录外地震目录1053个,与台网给出目录相比增加了0.18倍.为比较本次地震检测前后的最小完整性震级以及活动性b值,利用包络差峰值振幅对检测地震的震级进行了估算,根据检测后的余震目录得到九寨沟地震余震序列的最小完整性震级为M_L 0.9,地震活动性b值为0.80左右,较原目录的M_L 1.2和0.85均有降低.     

青海祁连MS5.2地震微震检测与目录完备性研究

利用模板匹配方法对2015年11月23日青海省祁连县M_S5.2地震进行遗漏地震检测研究,由于主震后短时间内目录中遗漏事件较多,故对主震后1天的连续波形进行检测。主震后1天内青海测震台网记录到的余震个数(包括单台)共62个,选取主震后M_L1.0以上余震30个作为模板事件,通过匹配滤波的方式扫描出遗漏地震31个,约为台网目录给出的0.5倍。基于包络差峰值振幅与震级的线性关系估测检测事件的震级参数,最后将检测后的余震目录与台网余震目录在主震后1天内的最小完备震级进行对比分析,结果发现检测后最小完备震级从M_L1.2降到了M_L0.7,得到青海测震台网在祁连地区最小完整性震级为M_L0.7。     

2017年九寨沟MS7.0地震序列活动特征分析

中强地震序列的主震发生后,短时间内受台站距震中较远、尾波干扰和波形重叠等因素的影响,往往会遗漏大量的地震,而地震目录的完整性会直接影响到震后趋势判定和余震序列特征分析的科学性和可靠性。本文利用基于GPU加速的模板匹配方法对2017年8月1～12日的连续波形进行扫描计算,检测九寨沟M_S7.0地震前后遗漏的地震事件,选取台网目录中信噪比较高的1033个地震事件作为模板,在主震前7天至震后5天期间识别出4854个检测地震事件,为台网可定位目录的3.3倍,除去对台网单台地震事件的修正外,还检测到1797个遗漏地震事件,将完备震级从1.6级降低到1.4级。基于补充了遗漏地震的完整地震目录,对2017年8月8日九寨沟M_S7.0地震序列活动特征进行分析。结果表明,前震序列在主震前短时间内出现了地震活动的密集增强,b值也显示为低值状态,可能是深部断层发生破裂之前的加速蠕动的结果。随着时间的推移,余震序列的完备震级逐渐下降并趋于稳定,b值存在缓慢升高的趋势,未来较长时期内余震序列仍将处于持续衰减的状态。     

2013年3月涿鹿微震群遗漏地震事件检测和发震构造分析

震群活动时,短时间发生大量地震,不同地震事件的记录波形相互交叠影响,易造成地震目录的遗漏,对震群发震构造分析等研究带来不利的影响.本文针对2013年3月3日至5日在河北涿鹿发生的微震震群,利用匹配滤波技术,以地震台网观测目录所记录地震事件的波形为模板,在连续波形记录中搜索与模板相似的信号,从而检测台网目录遗漏的地震.利用波形互相关标定新检测到地震事件的P波和S波到时,进而对其震中位置和震级做出估计.计算结果显示,通过互相关扫描检测到52个地震台网常规分析遗漏的地震,约为地震目录给出的45个事件的1.16倍.检测到的遗漏地震震级估算为M_L0.1~0.9,通过震级-频次统计分析,加入遗漏地震后地震目录的完整性在M_L0.3~0.8范围内有较明显的改善.根据地震事件精定位结果,推测此次震群的发震构造为北西走向倾角较大的断层,施庄断裂为发震构造的可能性较大.     

利用匹配滤波技术研究2013年灯塔MS5.1地震序列迁移

1 研究背景 地震序列活动的增强,尤其是震群活动或主震后的余震活动,短时间内区域地震台网会记录到大量波形资料,距震中稍远的台站常出现波形耦合现象,震级较小的地震波形常淹没在震级较大地震的尾波里,导致在地震编目中容易遗漏一定数量的地震,给利用余震序列研究发震构造带来不利影响,而且编目的完整性和震源参数的准确性对后续震情跟踪和相关研究结果的可靠性有重要影响(Woessner and Wiemer,2005;冯建刚等,2012).因此,检测遗漏地震,并测定震中参数是一个值得研究的问题.     

基于波形相似性的单台地震目录校正

本文基于波形相似性提出了一种适用于地震序列中单台事件目录校正的方法。以福建仙游地震序列为研究对象,应用双差地震定位（HypoDD）算法对多台事件进行重定位,依据记录波形相似性将仙游地震序列事件划分为161个事件丛集,以多台事件重定位目录为参考,对单台事件目录进行校正。结果表明,本方法能够有效地实现单台事件目录的合理校正,有助于提升地震序列目录的完整性和地震目录的实用性。以HypoDD重定位目录为参考的校正目录可用于地震参数计算。分析校正后的仙游地震序列目录认为,序列中的单台事件主要由金钟水库库区局部应力调整诱发。      

2016年唐山开平ML4.7地震活动性分析

2016年8月至2016年10月唐山余震区小震活动增强,尤其是8月21日至9月22日期间,其中9月10日发生震级最大地震为唐山开平ML4.7地震。为了解中强地震发震特征及地震序列变化,基于模板匹配定位识别技术,利用2016年8月1日至2016年10月31日连续波形进行唐山开平地震前后微震检测和定位。从研究区域内中国地震台网目录中选择283个地震事件作为模板事件,基于连续波形共检测识别出878个地震事件,其数量约为模板事件的3倍,约为中国地震台网目录的1.4倍,检测后的地震目录完备震级由原来的ML0.2提升至ML-0.1。利用双差定位方法对847个地震事件进行重定位,结果显示震中附近小震活动频繁,且在唐山开平地震震中周边分布较集中。分析震群分布特征认为唐山开平地震属于唐山老震区正常地震活动,发震断层为唐山-古冶断裂的可能性较大。     

基于Gutenberg-Richter定律快速估算最大余震震级:以2017年九寨沟MS7.0地震为例

针对九寨沟M_S7.0地震之后不同时间段的余震序列目录,利用推定最大余震震级,给出了实际最大余震震级的估计值。结果表明,推定最大余震震级随主震后时间尺度的延长而趋于稳定,且该值与实际发生的最大余震的震级一致。需要强调的是,就九寨沟地震序列而言,当余震数据较为完备时,采用主震后较短时间段内（1～2天）的余震目录就可以较准确地估算出主震区域内可能发生的最大余震震级。实际上,主震后12h（0.5天）的余震数据已完全可以给出最大余震震级的有效下限。此外,计算中我们采用了里氏震级M_L和面波震级M_S的余震目录,结果显示,2种震级类型目录的估算结果完全一致,表明利用推定最大余震震级估算实际最大余震震级的方法不受震级类型的影响。据此,该最大余震震级快速评估方法可进一步推广应用于我国大陆地区中强震后强余震灾害分析评估中。目前的拟合技术也显示出随着测震技术的不断进步以及余震识别能力的提高,快速评估方法可以在主震后短时间（<1天）内准确地预测可能发生的最大余震震级。     

2014年云南鲁甸M_S6.5地震早期余震序列遗漏地震检测与衰减特征分析

本文对2014年8月3日云南鲁甸6.5级地震早期余震序列中因波形叠加而遗漏的地震进行了检测,并基于修正的大森公式对序列时间域衰减特征进行了分析.首先利用巧家台连续地震记录波形的高频包络差检测主震发震至当日24点约7.5小时内目录遗漏的地震事件.检测到目录遗漏的ML2.0以上地震159个,最大震级为ML4.2,其中ML4.0以上遗漏地震1个,ML3.0以上遗漏地震3个.主震后1000秒内检测到遗漏地震37个,约为目录给出余震数量的1.76倍.使用补充遗漏地震的目录,基于修正的大森公式拟合主震后10天余震频度随时间的变化.增补遗漏地震后拟合p值由0.85上升至0.99,c值和K值无明显变化,表明鲁甸地震早期余震序列衰减速率高于仅依靠目录得到的结果.遗漏地震检测对于估计早期余震序列衰减速率估计影响较大,早期余震的缺失易造成序列衰减速度的低估.     

灌云地震台面波震级偏差分析

选取2019年1月—2021年12月灌云地震台远震事件记录,将灌云地震台测算的面波震级与中国地震台网中心统一正式报目录震级进行对比,计算二者面波震级偏差平均值,并分析震级偏差与震级、震中距、震中方位之间的关系。结果表明：灌云地震台面波震级与中国地震台网中心统一正式报目录震级平均偏差为-0.08,总体偏差较小,表明该台地震记录精度较高。     

邢台强余震前的异常波形

本文通过对(1970年～1985年)红山地震台所记录的台站周围地震资料的分析,发现较大地震前的小震大多呈现一种异常波形。我们对这种波形的特征进行了分析,并从统计的角度探索了它作为强余震前兆信息的可能性。结果表明,异常小震的震中与未来较大余震的震中相近;异常地震与强余震之间的时间间隔在数分钟至几天之间,异常小震的震级与未来较大地震的震级无明显的对应关系。最后,对产生这种波形的可能机制从震源区介质。应力场方面作了初浅的讨论。     

2020年新疆伽师MS6.4地震震源区 微震检测与目录完备性分析

中强地震发生后,地震检测因受到尾波的干扰可能会遗漏部分微震事件,影响地震目录的完备性。文章利用波形模板匹配方法对2020年新疆伽师M_S6.4地震序列开展微震检测,相比原始的中国地震台网中心统一地震目录,新检测出1 756个微震事件,地震数量增加了1.3倍。基于检测后的余震目录计算最小完备震级为M_L1.2,地震活动性b值为0.76,较原始目录的M_L1.6和0.77均有所降低。通过伽师震源区地震序列活动特征分析,结果表明前震序列在主震前短时间内(前36小时)出现地震活动的密集增强,相应的b值显示为低值;主震发生后地震序列完备震级较高,随着时间的推移,完备震级缓慢降低并趋于稳定,并且呈周期性的波动。本研究提高了伽师震源区地震目录的完备性,为精细化描述该地区地震序列时空演化特征提供了关键数据基础。     

2017年9月临城震群遗漏地震检测及发震构造初探

以河北测震台网记录的连续波形数据为基础,采用模板匹配滤波方法,检测2017年9月4日—6日河北临城震群遗漏地震。利用波形互相关扫描,检测到19次被测震台网常规分析遗漏的地震,震级范围为M_(L )0.0—1.2,并标定新检测地震事件的P波和S波到时,估计其震中位置及震级。通过地震精定位和最大地震震源机制解,认为此次临城地震序列的主要发震构造应为近EW向活动断层。     

点过程模型统计方法在地震学中的应用

引言地震目录大部分由多震国家地震部门定期出版。目录中给出每一个地震的震中、震源深度、发震时刻和震级大小。这种来自震源区内部的记录反映出地震时间序列是极其复杂的。地震月频次的大幅波动变化、     

2013年甘肃岷县漳县M_S6.6地震余震序列时域衰减特征分析

本文基于补充遗漏地震事件后的地震目录和修正的大森公式,对2013年7月22日甘肃岷县漳县6.6级地震余震序列时间域衰减特征进行了分析.鉴于主震后短时间内目录遗漏的余震较多,首先利用岷县台连续地震记录波形的高频包络差,检测主震后3h内目录遗漏的地震.经检测在主震后3h内共发现目录遗漏的ML1.0以上地震139个,最大震级为ML3.6.主震后1000s内检测到遗漏地震69个,约为目录给出余震数量的6倍.而后使用补充遗漏地震的目录,基于修正的大森公式分别拟合余震频度和余震地震矩随时间的变化.结果显示拟合p值约为1.07,表明岷县漳县地震余震序列衰减速率与全球平均水平接近,而未补充遗漏地震的频度拟合会造成余震序列衰减速率的低估.利用高频包络差直接计算地震频度曲线,通过三种衰减模式对地震频度曲线拟合参数比较,未观察到岷县漳县地震主震后存在早期余震缺失现象.分析认为,加入遗漏地震可以提高余震频度拟合估计衰减速率结果的准确度和精度,若缺少遗漏地震检测结果,则使用地震矩拟合所得衰减速率结果准确度较优,但需充分考虑其精度上的误差.在分析余震序列衰减特征的实际研究工作中,需根据地震目录完整性选择适当的拟合方法.     

川滇地区强余震的震级分布和距离分布特征

研究了川滇地区1966年以来12次主震震级MSge;6.5的地震序列中,主震与强余震（本文定义为所有MSge;5的余震）震级差分布特征和强余震与主震距离的分布特征. 结果表明,强余震与主震震级差服从截断的指数分布,据此推导出了强余震与主震震级差的概率密度函数; 强余震距离分布的优势范围是距主震10——39 km,且强余震与主震震中的距离服从正态分布.      

全球大震和中国及邻区中强震地震活动(2013年7-8月)

本目录中的地震参数来自"中国地震台站观测报告"(简称"月报").其中,国内及邻区给出M≥4.7的事件,全球给出M≥6.0的事件."月报"由中国地震台网中心按月做出.本目录中的发震时刻采用协调世界时(UTC);为了方便中国读者,也给出北京时(BTC),震中位置除给出经纬度外,还给出参考地区名,它仅用作查阅参考,不包含任何政治意义;还给出测定震源位置的台数(n)和标准偏差(SD).面波震级M_s是对中周期宽频带SK地震仪记录,采用北京台1965年面波震级公式M_s=lg(A_H/T)+     

全球大震和中国及邻区中强震地震活动(2015年1—2月)

<正>本目录中的地震参数来自"中国地震台站观测报告"(简称"月报").其中,国内及邻区给出M≥4.7的事件,全球给出M≥6.0的事件."月报"由中国地震台网中心按月做出.本目录中的发震时刻采用协调世界时(UTC);为了方便中国读者,也给出北京时(BTC).震中位置除给出经纬度外,还给出参考地区名,它仅用作查阅参考,不包含任何政治意义;还给出测定震源位置的台数(n)和标准偏差(SD).面波震级MS是用中周期宽频带SK地震仪记录,采用北京台1965年面波震级公式MS=lg(AH/T)+     

全球大震和中国及邻区中强震地震活动(2014年9—10月)

<正>本目录中的地震参数来自"中国地震台站观测报告"(简称"月报").其中,国内及邻区给出M≥4.7的事件,全球给出M≥6.0的事件."月报"由中国地震台网中心按月做出.本目录中的发震时刻采用协调世界时(UTC);为了方便中国读者,也给出北京时(BTC).震中位置除给出经纬度外,还给出参考地区名,它仅用作查阅参考,不包含任何政治意义;还给出测定震源位置的台数(n)和标准偏差(SD).面波震级Ms是用中周期宽频带SK地震仪记录,采用北京台1965年面波震级公式Ms=lg(A_H/T)+