2025-11-25T20:55:18.544254

Earthquake Forecasting with ETAS.inlabru

Zhong
The ETAS models are currently the most popular in the field of earthquake forecasting. The MCMC method is time-consuming and limited by parameter correlation while bringing parameter uncertainty. The INLA-based method "inlabru" solves these problems and performs better at Bayesian inference. The report introduces the composition of the ETAS model, then provides the model's log-likelihood and approximates it using Taylor expansion and binning strategies. We also present the general procedure of Bayesian inference in inlabru. The report follows three experiments. The first one explores the effect of fixing one parameter at its actual or wrong values on the posterior distribution of other parameters. We found that $α$ and $K$ have an apparent mutual influence relationship. At the same time, fixing $α$ or $K$ to its actual value can reduce the model fitting time by more than half. The second experiment compares normalised inter-event-time distribution on real data and synthetic catalogues. The distributions of normalised inter-event-time of real data and synthetic catalogues are consistent. Compared with Exp(1), they have more short and long inter-event-time, indicating the existence of clustering. Change on $μ$ and $p$ will influence the inter-event-time distribution. In the last one, we use events before the mainshock to predict events ten weeks after the mainshock. We use the number test and Continuous Ranked Probability Score (CRPS) to measure the accuracy and precision of the predictions. We found that we need at least one mainshock and corresponding offspring to make reliable forecasting. And when we have more mainshocks in our data, our forecasting will be better. Besides, we also figure out what is needed to obtain a good posterior distribution for each parameter.
academic

ETAS.inlabru による地震予測

基本情報

  • 論文ID: 2510.13930
  • タイトル: Earthquake Forecasting with ETAS.inlabru
  • 著者: Ziwen Zhong
  • 分類: stat.AP（統計学 - 応用）
  • 発表日: 2025年10月15日（arXiv プレプリント）
  • 論文リンク: https://arxiv.org/abs/2510.13930v1

要旨

本論文は、ETAS モデルを用いた地震予測を研究し、従来の MCMC 法の計算時間が長く、パラメータ間の相関性が強いという問題に対して、INLA ベースの「inlabru」法を用いたベイズ推定を提案している。論文は 3 つの実験を通じて方法の有効性を検証している：(1) 固定パラメータが他のパラメータの事後分布に与える影響を探索し、α と K パラメータ間に明らかな相互影響関係があることを発見、(2) 実データと合成カタログの標準化事象間隔時間分布を比較、(3) 主震前の事象を用いて主震後 10 週間以内の地震事象を予測し、数量テストと連続ランク確率スコア (CRPS) を用いて予測精度を評価している。

研究背景と動機

問題の定義

地震は人命と財産に対する脅威となる自然災害の一つであり、突発性が高く、迅速な回避が困難である。現在のところ、主震予測という研究の最高理想の実現はまだ達成されていないが、余震予測の分野では良好な進展が見られている。信頼性の高い余震予測は、政府と救援関係者がより効果的に震後の救援活動を計画するのに役立つ。

重要性

余震予測の精度は以下に直結している：

  1. 被災者により多くの生存機会をもたらす
  2. 救援関係者の生命安全をより良く保護する
  3. 震後の救援資源配置を最適化する

既存手法の限界

  1. 最尤推定法 (MLE)：方法は直接的で単純だが、パラメータ推定の不確実性を得ることができない
  2. マルコフ連鎖モンテカルロ法 (MCMC)：ベイズ法として事後分布を得ることができるが、以下の問題がある：
    • 計算時間が長い
    • モデルパラメータ間の高度な相関性
    • 良好なパラメータ事後分布の取得を保証することが難しい

研究動機

「inlabru」法は積分入れ子ラプラス近似 (INLA) に基づいており、以下が可能である：

  • 事後分布の近似を通じてモデル適合時間を大幅に削減
  • 共分散パラメータと潜在変数間の依存関係の近似が線形に近い場合に優れた性能を発揮
  • 内部的にガウス近似を推定し、ベイズ推定の効果を向上させる

核心的貢献

  1. 方法論的貢献：INLA ベースの inlabru 法を ETAS モデルの地震予測に適用し、従来の MCMC 法の計算効率の問題を解決
  2. パラメータ相互関係の発見：系統的な実験を通じて ETAS モデル内のパラメータ間の相互影響関係を明らかにし、特に α と K パラメータの強い相関性を発見
  3. 計算効率の向上：α または K パラメータを真の値に固定することで、モデル適合時間を 50% 以上削減できることを実証
  4. 予測能力の検証：実際の地震データに基づく予測評価フレームワークを構築し、信頼性の高い予測に必要な最小データ要件を特定
  5. 理論的洞察：良好なパラメータ事後分布を得るために必要な条件を明確にし、実際の応用に指針を提供

方法の詳細

タスク定義

入力：発生時刻 t、位置 s、震度 m を含む過去の地震事象系列 出力：将来の時間窓内の地震事象の確率分布予測 制約：自己励起点過程の ETAS モデルフレームワークに基づく

モデルアーキテクチャ

ETAS モデルの基礎

ETAS モデルは Hawkes 過程に基づき、その強度関数は以下の通りである：

λ_ETAS(t,m|H_t) = (μ + Σ K·e^(α(m_h-M_0))·(t-t_h/c + 1)^(-p))·π(m)

ここで：

  • μ: 背景発生率
  • K: 生産率、トリガー関数の強度を制御
  • α: 震度スケーリング、過去の事象の震度影響を制御
  • c: 時間オフセット
  • p: 余震減衰パラメータ
  • π(m): Gutenberg-Richter 則に従う震度分布

対数尤度の近似

Taylor 展開とビニング戦略を用いて対数尤度を近似：

L(θ|H) = -Λ_0(T_1,T_2) - Σ Σ Λ_i(t_j^(bi), t_{j+1}^(bi)) + Σ log λ(t_i|H_{t_i})

ベイズ推定フロー

  1. 初期パラメータ θ_0 を用いて Taylor 級数近似対数尤度関数を取得
  2. 対数事前分布と加算し、事後分布の最頻値 θ_1* を計算
  3. 線探索によるパラメータ更新：θ* = αθ_0 + (1-α)θ_1*
  4. 収束まで反復（差異が 1% 未満）

技術的革新点

  1. ビニング戦略の最適化：適応的ビニングを採用し、左側の狭いビンで精度を向上させ、右側の広いビンで計算量を削減
  2. 逆確率積分変換：ガウス分布サンプルを目標事前分布サンプルに変換
  3. 線形近似：Taylor 展開を通じた効率的な尤度関数近似の実現

実験設定

データセット

  1. Aquila 地震：合成カタログ生成とパラメータ分析に使用
  2. Amatrice 地震（2016-2017 年イタリア）：3 つの主震を含む実地震系列
  3. 合成カタログ：Aquila 地震パラメータに基づいて生成、制御実験に使用

評価指標

  1. 数量テスト (N-test)：予測事象数の精度を評価 
    δ_2 = |{N_j|N_j ≤ N_obs, j=1,...,m}|/m
  2. 連続ランク確率スコア (CRPS)：精度と精密性を同時に評価 
    S(F,N_true) = Σ(F(k) - I(N_true ≤ k))²

比較手法

  • 完全 ETAS モデル vs 単一パラメータ固定モデル
  • 異なる過去主震数のモデル比較
  • 実データ vs 合成データの分布比較

実装の詳細

事前分布：

  • μ ~ Gamma(0.3, 0.6)
  • K, α, c ~ U(0, 10)
  • p ~ U(1, 10)

初期値：μ₀=0.5, K₀=0.1, α₀=1, c₀=0.1, p₀=1.1

実験結果

主要な結果

実験 1：パラメータ固定効果

主要な発見：

  • α と K に明らかな相互影響関係が存在：一方を真の値に固定すると、他方の事後分布が改善される
  • c と p 間にも相互影響関係が存在
  • μ は相対的に独立しており、他のパラメータを固定しても μ の事後分布に影響しない

計算効率の向上：

モデル構成カタログ 1カタログ 2カタログ 3
完全モデル41.73s46.51s44.92s
K 固定26.00s20.95s21.47s
α 固定18.25s7.75s15.15s

実験 2：事象間隔時間分布

  • 実データと合成カタログの標準化事象間隔時間分布は基本的に一致
  • 指数分布 Exp(1) と比較して、より多くの短間隔と長間隔事象が存在し、クラスタリングの存在を示唆
  • μ の増加により分布は Exp(1)（均一ポアソン過程）に近づく
  • p の増加は短間隔事象数を著しく増加させる

実験 3：予測性能評価

Amatrice 地震予測結果：

  • 第 2 次主震を使用：全 10 週間で過予測が存在するが、真の値は 95% 信頼区間内
  • 第 3 次主震を使用：δ₂ 値は約 0.5 に近く、予測がより正確
  • CRPS は第 3 次主震予測が第 2 次主震予測を大幅に上回ることを示す

アブレーション実験

過去主震数の影響：

  1. 主震なし：μ 以外のすべてのパラメータの事後分布が非常に悪い
  2. 1 つの主震：基本的に使用可能な予測が得られるが、精度は限定的
  3. 2 つの主震：予測性能が著しく改善され、特に α パラメータの事後分布が改善

ケース分析

パラメータ推定の難易度ランキング：

  1. μ（最も容易）：少量の事象で推定可能
  2. c, p（中程度）：より多くの地震事象が必要であり、静穏なカタログでも推定可能
  3. α, K（最も困難）：少なくとも 1 つの主震と対応する余震が必要；より多くの主震は α の事後分布を著しく改善できる

実験的発見

  1. パラメータ相関性：α-K と c-p に強い相関性が存在し、共同で考慮する必要がある
  2. データ要件：信頼性の高い予測には少なくとも 1 つの完全な主震-余震系列が必要
  3. 予測改善：過去主震数の増加は予測性能を著しく向上させることができる
  4. 計算最適化：事前情報の合理的な利用により計算効率を大幅に向上させることができる

関連研究

主要な研究方向

  1. ETAS モデルの発展：Ogata (1988) の原始モデルから各種拡張版へ
  2. パラメータ推定法：MLE から MCMC を経て INLA ベース法への進化
  3. 地震予測評価：数量テスト、CRPS などの評価指標の発展

本論文と関連研究の関係

  • Hawkes 過程理論の基礎の上に構築
  • ETAS モデルの核心的思想を継承
  • Serafini et al. (2023) と Naylor et al. (2023) の inlabru 応用の基礎上での深化研究

相対的優位性

  1. 従来の MCMC 法と比べて計算効率が著しく向上
  2. パラメータ間の相互関係を系統的に分析
  3. 実用的な予測評価フレームワークを提供

結論と考察

主要な結論

  1. 方法の有効性：inlabru 法は ETAS モデルで優れた性能を発揮し、計算効率を著しく向上させる
  2. パラメータ関係：α と K の相互依存関係はモデル性能に影響する主要因である
  3. 予測能力：モデルは信頼性の高い余震予測能力を有し、予測品質は過去データの豊富さに伴い向上
  4. データ要件：信頼性の高い予測には少なくとも 1 つの主震-余震系列が必要

限界

  1. データ制限：合成カタログは単一主震のみを含み、結論の一般性を制限
  2. 事象間隔の定義：隣接事象間隔のみを考慮し、主震クラスタ内部構造を考慮していない
  3. 計算上の課題：主震数の増加に伴いモデル適合時間が著しく増加し、収束が困難
  4. 特性分析の不足：主震数以外の地震特性が事後分布に与える影響はまだ十分に探索されていない

今後の方向

  1. パラメータ更新戦略：増分学習法を開発し、前の主震の事後分布を次の事前分布として使用
  2. 複数主震実験：より多くの主震を含む合成カタログで実験を繰り返す
  3. 事象クラスタリング分析：主震クラスタ内部および間の事象間隔分布を研究
  4. 特性影響研究：各種地震特性がモデル性能に与える影響を系統的に分析

深度評価

利点

  1. 方法の革新性：INLA 法を ETAS モデルに成功裏に適用し、実際の計算問題を解決
  2. 実験設計の完全性：3 つの実験が異なる角度から方法性能を包括的に評価
  3. パラメータ分析の深さ：ETAS モデルパラメータ間の相互関係を系統的に明らかにする
  4. 実用価値の高さ：実際の地震予測応用に明確なデータ要件指針を提供
  5. 結果の信頼性：実地震データを用いた検証により、結論は実際的意義を有する

不足

  1. 理論分析の不足：パラメータ相関性に対する理論的説明が不足
  2. 実験範囲の限定：主にイタリア地震データに基づき、地域代表性が不足
  3. 方法比較の不十分さ：他の現代的ベイズ法との詳細な比較がない
  4. 不確実性の定量化：予測不確実性の分析がまだ十分に深くない

影響力

  1. 学術的貢献：地震学と統計学の交差領域に新しい方法論を提供
  2. 実用価値：実際の地震予測とリスク評価に直接応用可能
  3. 再現性：オープンソース R パッケージ inlabru に基づき、良好な再現性を有する
  4. 推進可能性：方法は他の自己励起点過程応用領域に拡張可能

適用シーン

  1. 余震予測：震後短期の余震活動予測
  2. 地震リスク評価：過去データに基づく地震リスク評価
  3. 応急対応計画：震後救援に科学的意思決定支援を提供
  4. 地震学研究：地震活動の時空分布規則の理解に使用

参考文献

主要な参考文献には以下が含まれる：

  • Ogata (1988): ETAS モデルの原始提案
  • Hawkes (1971): Hawkes 過程の理論基礎
  • Rue et al. (2017): INLA 法の理論フレームワーク
  • Serafini et al. (2023): Hawkes 過程における inlabru の応用
  • Naylor et al. (2023): ETAS モデルのベイズモデリング

総合評価：これは地震予測領域において重要な実用価値を有する論文であり、ETAS モデルの計算効率問題を成功裏に解決し、系統的な実験を通じて価値あるパラメータ分析と予測評価を提供している。理論的深さと実験の広さにおいてなお改善の余地があるが、その方法の実用性と結論の操作可能性により、本領域の重要な貢献となっている。