[论文解读] Dynamics of the Markov Time Scale of Seismic Activity May Provide a Short-Term Alert for Earthquakes

Reza Mansouri, Muhammad Sahimi+10 more|ArXiv.org|Oct 6, 2005

Complex Systems and Time Series Analysis参考文献 3被引用 3

一句话总结

该论文提出了一种基于地震数据中马尔可夫时标（$t_M$）和扩展自相似性（$T_1$）动力学的短期地震预测方法。当 $t_M$ 和 $T_1$ 同时超过临界阈值时，会触发可靠警报，该方法在伊朗地震测试中实现了近乎零的误报率。

ABSTRACT

We propose a novel method for analyzing precursory seismic data before an earthquake that treats them as a Markov process and distinguishes the background noise from real fluctuations due to an earthquake. A short time (on the order of several hours) before an earthquake the Markov time scale $t_M$ increases sharply, hence providing an alarm for an impending earthquake. To distinguish a false alarm from a reliable one, we compute a second quantity, $T_1$, based on the concept of extended self-similarity of the data. $T_1$ also changes strongly before an earthquake occurs. An alarm is accepted if {\it both} $t_M$ and $T_1$ indicate it {\it simultaneously}. Calibrating the method with the data for one region provides a tool for predicting an impending earthquake within that region. Our analysis of the data for a large number of earthquakes indicate an essentially zero rate of failure for the method.

研究动机与目标

• 基于地震时间序列的统计动力学，开发一种短期地震预测方法。
• 通过马尔可夫过程分析，区分真实前兆信号与背景噪声。
• 通过要求两个独立统计指标 $t_M$ 和 $T_1$ 同时发出信号，减少误报。
• 提供与地震震级相关的可量化预警时间，实现可操作的早期预警。
• 利用历史地震数据在区域范围内校准该方法，以实现实际应用中的实时监测。

提出的方法

• 该方法将马尔可夫时标 ($t_M$) 计算为地震数据可被建模为马尔可夫过程的时间段，通过检验条件概率的 Chapman-Kolmogorov 方程来确定。
• 通过比较不同时间延迟下直接测量的转移概率与通过 Chapman-Kolmogorov 方程计算出的转移概率，评估马尔可夫假设的有效性。
• 第二个指标 $T_1$ 基于扩展自相似性的概念，用于检测地震波动在大地震前的异常标度行为。
• 利用特定区域内先前发生地震的数据，校准临界阈值 $t_{cM}$ 和 $T_{c1}$。
• 仅当 $t_M > t_{cM}$ 且 $T_1 > T_{c1}$ 同时满足时才触发警报，从而最大限度减少误报。
• 通过对数关系估算预警时间 $t_a$：$ ext{log } t_a = -1.35 + 2.4 ext{ log } M$，其中 $t_a$ 以小时为单位，$M$ 为地震震级。

实验结果

研究问题

• RQ1 地震活动的马尔可夫时标 ($t_M$) 是否可作为大地震的可靠短期前兆？
• RQ2 扩展自相似性指标 $T_1$ 是否能通过减少误报来提高基于 $t_M$ 的警报可靠性？
• RQ3 是否存在地震震级与该方法提供的预警时间之间的定量关系？
• RQ4 该方法是否可区域化校准，以提供准确、实时的短期地震预警？
• RQ5 当 $t_M$ 和 $T_1$ 必须同时触发警报时，该方法的失败率是多少？

主要发现

• 当 $t_M$ 和 $T_1$ 同时发出警报时，该方法在两年分析期内仅出现两次失败，失败率几乎为零。
• 对于震级 $M = 4.5$ 的地震，根据关系式 $ ext{log } t_a = -1.35 + 2.4 ext{ log } M$，预测的预警时间为约 2 小时。
• 预警时间 $t_a$ 随震级增加而增加，表明较大震级的地震具有更长的预警时间，当 $M > 5.7$ 时，$t_a$ 可达 10 小时。
• 该方法成功为巴勒德地震提供了七小时的预警，尽管该地区无前震，证明其对前兆信号的高度敏感性。
• 位于深度超过 50 米且垂直于活动断层的测站，其数据相关性更强，警报更准确。
• 该方法通过 173 次伊朗地震（震级介于 3.2 至 6.3 之间）得到验证，证实其在广泛震级范围内的鲁棒性。