基于数据驱动的时间序列 b 值计算新方法(TbDD): 以2021年云南漾濞 M S 6.4地震序列为例

时间序列的b值在天然地震和工业开采诱发地震的危险性分析中具有重要的应用潜力，但长期以来受到计算规则设置的人为主观性、计算结果的可靠性和时序对突变识别精度不高等问题影响，制约了不同结果的可比较性和共识性科学认识的提炼.本文借鉴基于数据驱动（data-driven）的地震活动参数计算思路，采用连续函数形式的OK1993模型、时间轴随机段落划分、贝叶斯信息准则模型选择等技术环节，构建了基于数据驱动的时间序列b值计算新方法TbDD.利用合成地震目录的理论测试，并分别与固定地震数目的窗长和步长、固定地震数目的步长和累积窗长等传统的固定窗口法进行了比较研究.结果表明，TbDD方法可较好地还原合成地震目录的b0值输入参数，在计算规则设置的客观性和对b值突变过程的准确识别上具有明显优势.此外，我们还对新近发生的2021年5月21日云南漾濞MS6.4地震序列进行了实际案例应用.结果显示，此次序列的b值在MS6.4主震前为0.7左右、震前20 h出现了约0.1幅度的下降，表明在序列发生前震区的差应力水平较高.而b值在MS6.4主震发生后起伏明显、逐渐增加至0.8左右，这一现象可能与震区在主震后早期较为剧烈的应力调整有关.进一步针对随机模型的数量以及时间轴的随机段落划分设置对TbDD方法b值计算结果的影响程度进行了测试，发现b值受随机模型数量影响较小、具备较强的稳定性，时间轴的随机段落划分设置可影响b值时序微观起伏变化的识别.本文发展的TbDD方法在对时间序列b值计算的准确性、余震趋势跟踪的高精度要求，以及工业开采诱发地震风险管控等领域有较好的应用潜力，所获得的2021年云南漾濞MS6.4地震序列的b值计算结果也对理解此次地震序列的孕育过程有参考价值.