2019年6月17日四川长宁MS6.0地震早期序列参数分析及强余震概率预测

2019年6月17日四川长宁发生M_S6.0地震, 之后发生了一系列的强余震, 为更好地分析此次地震的序列特征以及强余震的可预测属性, 采用国际上对复杂序列拟合相对较好的“传染型余震序列”(ETAS)模型以及基于Reseanberg-Jones(R-J)模型发展的Omi-R-J模型, 通过连续滑动、 拟合和余震发生率预测, 对地震序列的模型参数稳定性、 预测结果进行了比较研究, 并利用N-test、 T-test检验方法对预测结果进行了效能评估。 结果表明, 相比于其他中强震序列参数, 此次长宁M_S6.0地震序列参数中反映激发能力的αETAS较其他序列明显偏小, 而反映衰减能力的pORJ值和应力累积水平的bORJ值相对较小, 与此次余震序列丰富、 持续时间相对较长相吻合; ETAS和Omi-R-J模型对于复杂序列在[3.0, 3.5, 4.0]三个震级档的强余震仍具有一定的预测能力; 总体的“每个地震的信息增益”(IGPE)计算结果显示, ETAS模型略优于Omi-R-J模型, 前者或更适合复杂地震序列的余震预测。     

利用早期余震预测的Omi-R-J方法对2017年四川九寨沟MS7.0地震的应用研究

在地震序列的早期阶段获得科学可靠的余震预测结果,始终是地震预测研究的前沿课题.针对新近发展的Omi-R-J方法的适用性问题,以及相比于传统的Reseanberg-Jones （R-J）模型在地震序列参数拟合、余震预测效能上的差异等问题,本文以2017年四川九寨沟M_S7.0地震序列为例,利用多时间窗的连续滑动拟合、预测和检验,以及通过构建重采样随机地震序列目录的系统检验等方式,开展了应用研究和比较研究.结果表明：相比于R-J模型,Omi-R-J方法可在余震记录较不完整的地震序列早期阶段获得稳定、可靠的序列参数,p值、c值和b值的标准差相比于R-J模型参数明显减小.N-test方法检验结果表明,在R-J模型等传统方法无法获得预测结果的地震序列早期阶段,Omi-R-J方法预测结果可以较高的比例通过N-test检验,在后续时段的预测效果也明显优于R-J模型.利用随机地震序列目录的测试结果表明,余震记录完整性下降对Omi-R-J方法预测结果的影响相对较小,在全部时段的完整性下降的影响高于在部分时段的完整性下降.上述结果对进一步地将Omi-R-J方法应用在震后早期余震预测中具有一定科学借鉴意义.     

地震序列早期参数估计和余震概率预测中截止震级Mc的影响：以2013年甘肃岷县—漳县6.6级地震为例

针对“时空传染型余震序列”（英文简称ETAS）模型在地震序列参数的早期特征和余震短期概率预测研究中的应用问题,重点考察了不同截止震级M_c选取对结果的影响.以甘肃岷县—漳县6.6级地震序列的震后1.677天内的早期阶段为例,考察了ETAS模型和修正的Omori-Utsu公式的适用性问题,发现ETAS模型AIC值在各截止震级M_c下均小于修正的Omori-Utsu公式的结果,表明其适用效果更好.设定M_c=M_L1.0,1.1,…,2.5,分别考察了ETAS模型中α值和p值的稳定性,并与2013年芦山7.0级地震序列进行了对比.结果表明,M_c对α值的影响相对较小,p值影响较大.此外,对基于ETAS模型和“瘦化算法”的余震短期概率预测结果进行了N-test检验,结果表明M_c的设定对余震短期概率预测影响较大,对甘肃岷县—漳县6.6级地震,仅当M_c=M_L1.0或M_L1.1时可获得较好的预测结果.由此,在真正的“向前”的预测实践中,需要首先考察不同的M_c下的余震预测效果.     

威远MS 5.4地震序列活动特征与成因分析

通过对威远M_S 5.4地震背景与活动特征研究,表明自贡地震监测中心加强流动监测提高了威远M_S 5.4地震序列地震参数的精度和地震信息报送的准确性。威远M_S 5.4地震前经历了区域地震平静至活跃、再发震的过程,威远M_S 5.4地震为一般主—余型地震,其ρ、b值异常可能反映了其后将发生资中M_S 5.2地震。威远M_S 5.4地震序列及其附近的华蓥山断裂地震带和马边地震带震后时间扫描地震学参数均正常,威远M_S 5.4地震区近期再次发生更大地震的可能性较小;并认为威远M_S 5.4发震断裂为荣县—威远基底断裂,华蓥山断裂带及其附近荣县—威远—资中地区2018年开始的4—6级地震活跃与川西地区东昆仑断裂带2017年8月8日九寨沟M_S 7.0地震的发生存在呼应关系。     

能量震级及其测定

对于地震灾害与风险评估,人们更关注的是地震辐射能量E_S和能量震级M_e的大小,能量震级M_e反映震源动态特征,适合描述地震的潜在破坏性。本文介绍地震波能量E_S和地震矩M₀的物理意义及能量震级M_e的定义和测定方法,并测定得到2017年8月8日四川九寨沟M_S 7.0地震的能量震级M_e为6.3。     

2017年九寨沟MS7.0地震余震危险性概率分析

本文介绍了余震危险性概率分析的概念,完整列出其相关公式的解析表达式并阐述了余震危险性概率分析的具体方法。以2017年九寨沟M_S7.0地震为例,首先,求得该地震余震序列的相关参数,结果显示:九寨沟M_S7.0地震余震序列理论最大余震震级约为M_L5.3;b值约为0.784 1,明显低于中国西南地区同类型的其它地震,表明九寨沟地震余震区应力水平相对较高;p值约为1.109 7,明显高于中国西南地区同类型的其它地震,表明九寨沟地震余震序列随时间衰减较快。其次,计算此次地震余震与主震释放能量的比例关系,结果显示:在九寨沟M_S7.0地震事件中,截止到10月22日约99.69%的能量为主震所释放,0.31%的能量为余震所释放。最后,利用九寨沟M_S7.0余震序列参数结果并结合衰减关系,分别计算了主震后0—1天、1—10天、10—30天和90—100天内不同断层距内不同水平的峰值加速度和峰值速度的超越概率,结果显示:随着断层距增加,在相同时间间隔内超越概率的值呈下降趋势;在震后第一天内余震危险性最高,随着震后时间增加,相应的超越概率值呈现出明显下降趋势,表明九寨沟M_S7.0地震余震危险性主要来自早期余震。本文的相关成果可以为震后地震危险性分析提供参考,并为短时间内应急救援及灾后重建提供辅助决策意见。      

华北地区地震序列参数的分布特征

为系统地考察华北地区地震序列参数的分布特征,以期构建适合区域地震活动特征的短期概率预测模型和评估地震危险性,本文利用当前较为前沿的ETAS模型和R-J模型,采用连续滑动、多时段拟合的方式,对华北地区1970年以来的16个地震序列进行了参数拟合,并对参数的整体情况、参数之间的关系、参数与大地热流之间的关系等特征进行了分析研究.研究结果表明,主震发生后的早期阶段,地震序列参数变化相对较为剧烈,误差也较大;两种地震预测模型的序列参数呈现出一种优势分布特征,主要参数的平均值分别为α_ETAS=1.7404±0.3420,p_ETAS=0.9769±0.1396,a_OML=-1.6638±0.5284,b_OML=0.8312±0.1658,p_OML=0.9053±0.1527,这与国际上其他区域的研究具有很强的一致性;大多数情况下,ETAS模型参数p_ETAS高于R-J模型参数p_OML,平均偏高0.0716;序列参数α_ETAS与大地热流值整体上呈现一种负相关关系,激发次级余震的能力与大地热流具有一定的相关性.地震序列参数的分布特征以及与大地热流之间的关系对地震序列类型的判断,及基于地震序列参数构建定量化的短期概率预测模型,均有重要的应用价值.     

芦山两次强震序列活动特征及余震预测效能对比分析

利用可充分考虑小震信息的Omi-R-J模型,对2013年芦山M_S7.0地震和2022年芦山M_S6.1地震序列活动特征进行对比分析,发现芦山M_S7.0地震序列较芦山M_S6.1地震序列完备性随时间变化更加明显,震后早期检测能力较低,两次地震序列稳定时段检测得到的完整性震级均为1.8级左右。芦山两次地震序列p值差异较小,展示了相对正常的衰减过程;芦山M_S7.0地震序列k值明显大于芦山M_S6.1地震序列,这或许与芦山M_S7.0地震较为发育的余震或强余震有关;芦山M_S7.0地震序列b值小于芦山M_S6.1地震序列,表明M_S7.0地震之后芦山地区仍处于较高的应力状态,而M_S6.1地震之后应力处于相对较低的水平。利用N-test方法开展效能评估,结果显示,两次地震序列初期阶段余震发生率展示了较好的预测效能,基于震后初期1天内的数据开展的未来1天的余震预测中,仅有1...     

基于ETAS模型对青海地区余震预测效能评估的研究

为系统评估青海地区余震短期发生率的预测效能,以及构建适合地震活动特点的余震早期预测策略和预测指标体系,利用国际上当前较为前沿的时间序列ETAS模型和"瘦化算法"对青海地区2009年以来的8个地震的早期余震序列参数进行拟合,并利用N-test检验方法对预测结果进行回溯性的效能评估。研究表明:ETAS模型和"瘦化算法"对青海地区的余震发生率具有很好的预测能力,建议采用3天的预测时间窗,且对序列早期阶段进行应用,或可取得"最佳"的预测效果。     

岫岩地震序列视应力变化及其预测意义

选取辽宁省地震局沈阳数字化地震台网中的岫岩、 营口、 本溪、 北镇、 新民5个台的波形记录, 计算了1999年11月29日岫岩M_S5.4地震序列视应力值。 5个台资料的计算结果具有很好的一致性, 即岫岩序列的视应力变化特征具有显著的共性特征: 震前低视应力阶段、 震时高视应力阶段和震后逐渐降低的调整至稳定阶段, M_S5.4主震和5.1级最大强余震前视应力升高, 表明视应力可作为中强震及震后地震趋势预测的判定依据。     

余震的序列参数稳定性和余震短期发生率预测效能的连续评估——以2014年云南鲁甸6.5级地震为例

为考察目前国际上广泛使用、对真实地震序列描述最好的"传染型余震序列模型"(ETAS)在主震后的序列参数拟合、余震短期发生率预测的效能,本研究以2014年云南鲁甸MS6.5地震序列为例,采用滑动连续拟合与预测的方式,考察了ETAS模型参数的动态变化和余震短期发生率预测的实际效能.连续滑动拟合结果表明,在主震发生后的早期阶段,α值有明显的不稳定变化,在震后5.10天稳定在1.6~2.0;p值在震后25.00天内由1.07逐渐下降至0.78左右,其后稳定在0.72~0.85;b值在震后35.00天内逐渐由0.80增加至0.95,其后稳定在0.93~0.97.对连续滑动预测结果的N-test检验表明,余震发生率预测会出现部分失效现象,1天预测时间窗失效比例约为12%、3天预测时间窗失效比例为6%.建议可在震后早期采用1天的较短预测时间窗,而在序列参数较为稳定时段采用较长的3天预测时间窗.     

2020年4月1日四川石渠5.6级地震总结

系统梳理2020年4月1日四川石渠5.6级地震震前地震活动异常特征、地球物理观测异常以及区域构造情况,结果如下:①地震活动:石渠5.6级地震震中区域属于历史弱震区,震前出现前兆震群,也是其前震活动;②地球物理观测:出现8项异常,其中形变5项,电磁3项,且电磁异常对于该地震具有预测意义;③综合方法:地震发生前,震中附近存在W_q值异常。石渠5.6级地震发生在五道梁—长沙贡玛断裂带上,震源机制解显示为走滑型破裂。该序列类型为前震—主震—余震型,主震前存在明显的前震活动,余震较为丰富,序列活动呈阶段性衰减特征。综合分析认为,对于历史弱震区,可以通过分析地球物理观测异常和地震活跃时段的b值,为区域地震危险性评估提供依据,同时可利用地震序列参数h值和b值,对震后余震水平进行有效评估。     

2020年5月18日云南巧家MS 5.0地震总结

系统总结2020年5月18日云南巧家M_S 5.0地震发生前出现的地震活动和地球物理观测异常,其中：①地震活动异常：震前存在震群活动、Benioff应变比和地震发生率异常等中短期异常;②地球物理观测异常：震中200 km范围内分布重力、地磁、水位、水温、断层氢等测项,其中异常测项15项,大部分为趋势异常,异常台项比较低（5%）,短临异常不突出,巧家电磁波为地震的短期异常;③综合方法：多参数概率谱计算结果在震前出现异常。对巧家M_S 5.0地震序列进行跟踪研究,发现：主震与最大余震的震级差为1.6（M_L震级）,主震释放能量占序列总能量的99.64%,属主余型地震序列;主震震源机制显示为走滑型,根据余震分布,认为发震断层位于莲峰断裂和昭通-鲁甸断裂之间的NS向次级断裂。综合分析认为,在巧家M_S 5.0地震发生前,地震活动具有背景异常逐渐向中短期转变的过程,地球物理观测异常以趋势异常为主,短期异常不突出。     

2021年3月24日新疆拜城MS 5.4地震总结

系统梳理2021年3月24日新疆拜城M_S 5.4地震前地震活动异常特征、地球物理观测异常以及区域构造情况,结果如下：①地震活动：震前存在地震平静、高频、带状分布等异常;②地球物理观测：出现5项异常,其中形变3项,电磁2项,且形变异常对于该地震具有预测意义;③综合方法：地震发生前,震中附近存在多参数概率谱异常。拜城M_S 5.4地震发生在库车坳陷西南部,震源机制解显示为走滑型破裂。该序列类型为前震—主震—余震型,主震前存在前震活动,余震较少,序列活动呈持续衰减特征。综合分析认为,拜城M_S 5.4地震前异常较多,其中地震活动相关异常为后续中强地震的预测提供了有效依据。     

2021年3月19日西藏比如MS 6.1地震总结

2021年3月19日西藏比如M_(S )6.1地震发生在班戈—安多地区的安多盆地南缘断裂附近,震源机制解反演结果显示,本次地震为近东西走向断层产生正断兼走滑型错动的结果,与区域构造特征一致;截止至4月15日,地震序列跟踪发现,主震与余震震级差为1.4,主震释放能量占序列的98.96%,表明6.1级主震释放了此次地震序列的绝大部分能量。通过总结地震周边的地震活动异常发现,震前存在6级地震空区及3级地震条带交汇等异常,同时存在地震发生率、b值、D值、RTL值、地震调制比等地震学参数异常,地球物理观测统计显示,震中300 km范围内无地球物理观测测项。此外,比如地震打破了大陆地区长时间的5级以上地震平静,对后续6级地震的发生具有一定预测意义。     

2020年7月12日唐山古冶MS 5.1地震震源参数

利用河北及邻区地震台网提供的震相观测报告,使用双差定位方法,对2020年7月12日唐山古冶M_S 5.1地震序列进行重新定位,并基于部分地震台站记录的波形资料,采用近震全波形方法,反演得到主震震源机制解。精定位结果显示,此次地震序列展布长度约8 km,余震自主震位置向SW扩展,震源深度优势分布范围在10—18 km,发震断层面较陡,倾向SE。震源机制解显示,主震为一次走滑型事件,结合序列展布形态和地震活动背景,SW—NE向近似直立节面为可能发震断层面,与1976年唐山M_S 7.8主震断层特征基本一致。综上所述,古冶M_S 5.1地震为唐山M_S 7.8地震的又一次余震,属唐山老震区正常地震活动。     

2020年5月6日新疆乌恰5.0级和5月9日柯坪5.2级地震总结

2020年5月6日、5月9日,新疆地区南天山西段先后发生乌恰5.0级和柯坪5.2级地震,系统总结2次地震发生前出现的地震活动和地球物理观测异常,其中：①地震活动：震前存在调制地震集中、地震窗、5级以上地震成组等中短期异常;②地球物理观测：2次地震震中附近震前出现形变、电磁和流体观测异常,其中形变异常3项、电磁异常4项、流体异常1项,主要分布在柯坪5.2级地震震中附近。通过对2次地震序列进行跟踪,发现：乌恰5.0级地震余震较少,震后60天内共记录M_L 3.0以上余震4次,最大震级为M_L 4.5;柯坪5.2级地震后余震较丰富,震后60天内共记录M_L 3.0以上余震10次,最大震级为M_L 4.7,计算得到序列h值为1.6,b值为0.73。综合分析认为,2020年5月新疆地区2次5级以上地震前存在的地震活动异常较少,但区域地震活动水平较强,主要存在具有中短期指示意义的地球物理观测异常。