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 (सांख्यिकी - अनुप्रयोग)
  • प्रकाशन समय: 25 अक्टूबर, 2025 (arXiv प्रीप्रिंट)
  • पेपर लिंक: https://arxiv.org/abs/2510.13930v1

सारांश

यह पेपर ETAS मॉडल का उपयोग करके भूकंप पूर्वानुमान के लिए अनुसंधान करता है। पारंपरिक MCMC विधि की समय-खपत और मजबूत पैरामीटर सहसंबंध की समस्याओं को संबोधित करते हुए, यह INLA-आधारित "inlabru" विधि का उपयोग करके बेयेसियन अनुमान प्रस्तावित करता है। पेपर तीन प्रयोगों के माध्यम से विधि की प्रभावशीलता को सत्यापित करता है: (1) निश्चित पैरामीटर के अन्य पैरामीटर पश्च वितरण पर प्रभाव की खोज, α और K पैरामीटर में स्पष्ट पारस्परिक प्रभाव संबंध पाया गया; (2) वास्तविक डेटा और सिंथेटिक कैटलॉग के मानकीकृत घटना अंतराल समय वितरण की तुलना; (3) मुख्य भूकंप से पहले की घटनाओं का उपयोग करके मुख्य भूकंप के बाद दस सप्ताह के भीतर भूकंप की घटनाओं की भविष्यवाणी करना, संख्या परीक्षण और निरंतर रैंकिंग संभाव्यता स्कोर (CRPS) का उपयोग करके पूर्वानुमान सटीकता का मूल्यांकन करना।

अनुसंधान पृष्ठभूमि और प्रेरणा

समस्या परिभाषा

भूकंप मानव जीवन और संपत्ति को खतरे में डालने वाली प्राकृतिक आपदाओं में से एक है, जिसमें अचानकपन की मजबूत विशेषता और समय पर बचने की कठिनाई है। हालांकि वर्तमान में मुख्य भूकंप की भविष्यवाणी करना अभी भी संभव नहीं है, लेकिन आफ्टरशॉक पूर्वानुमान में अच्छी प्रगति हुई है। विश्वसनीय आफ्टरशॉक पूर्वानुमान सरकार और बचाव कर्मियों को भूकंप के बाद बेहतर तरीके से बचाव कार्य की योजना बनाने में मदद कर सकता है।

महत्व

आफ्टरशॉक पूर्वानुमान की सटीकता सीधे संबंधित है:

  1. प्रभावित लोगों के लिए अधिक जीवित रहने का मौका सुनिश्चित करना
  2. बचाव कर्मियों की जान की बेहतर सुरक्षा
  3. भूकंप के बाद बचाव संसाधनों का अनुकूलन

मौजूदा विधि की सीमाएं

  1. अधिकतम संभावना अनुमान (MLE): विधि सीधी और सरल है, लेकिन पैरामीटर अनुमान की अनिश्चितता प्राप्त नहीं कर सकती
  2. मार्कोव चेन मोंटे कार्लो (MCMC): बेयेसियन विधि के रूप में पैरामीटर पश्च वितरण प्राप्त कर सकती है, लेकिन निम्नलिखित समस्याएं हैं:
    • लंबी गणना समय
    • मॉडल पैरामीटर के बीच उच्च सहसंबंध
    • अच्छे पैरामीटर पश्च वितरण प्राप्त करना कठिन

अनुसंधान प्रेरणा

'inlabru' विधि एकीकृत नेस्टेड लाप्लास सन्निकटन (INLA) पर आधारित है, जो:

  • पश्च वितरण के सन्निकटन के माध्यम से मॉडल फिटिंग समय में काफी कमी ला सकती है
  • सहप्रसरण पैरामीटर और अव्यक्त चर के बीच निर्भरता संबंध लगभग रैखिक होने पर उत्कृष्ट प्रदर्शन करती है
  • आंतरिक गॉसियन सन्निकटन का अनुमान लगाती है, बेयेसियन अनुमान प्रभाव में सुधार करती है

मुख्य योगदान

  1. पद्धति संबंधी योगदान: ETAS मॉडल के भूकंप पूर्वानुमान में INLA-आधारित inlabru विधि को लागू करना, पारंपरिक 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. सिंथेटिक कैटलॉग: 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 मॉडल बनाम एकल पैरामीटर निश्चित मॉडल
  • विभिन्न ऐतिहासिक मुख्य भूकंप संख्या के साथ मॉडल तुलना
  • वास्तविक डेटा बनाम सिंथेटिक डेटा वितरण तुलना

कार्यान्वयन विवरण

पूर्व वितरण:

  • μ ~ 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) के करीब बनाती है (समान Poisson प्रक्रिया)
  • p में वृद्धि छोटे अंतराल घटनाओं की संख्या में काफी वृद्धि करती है

प्रयोग 3: पूर्वानुमान प्रदर्शन मूल्यांकन

Amatrice भूकंप पूर्वानुमान परिणाम:

  • दूसरे मुख्य भूकंप का उपयोग: सभी दस सप्ताह में अधिक पूर्वानुमान मौजूद है, लेकिन वास्तविक मान 95% आत्मविश्वास अंतराल के भीतर है
  • तीसरे मुख्य भूकंप का उपयोग: δ₂ मान 0.5 के करीब है, पूर्वानुमान अधिक सटीक है
  • CRPS तीसरे मुख्य भूकंप पूर्वानुमान को दूसरे से काफी बेहतर दिखाता है

विलोपन प्रयोग

ऐतिहासिक मुख्य भूकंप संख्या का प्रभाव:

  1. कोई मुख्य भूकंप नहीं: μ को छोड़कर सभी पैरामीटर पश्च वितरण बहुत खराब हैं
  2. एक मुख्य भूकंप: मूलतः उपयोगी पूर्वानुमान प्राप्त कर सकते हैं, लेकिन सटीकता सीमित है
  3. दो मुख्य भूकंप: पूर्वानुमान प्रदर्शन में काफी सुधार, विशेष रूप से α पैरामीटर के पश्च वितरण में

केस विश्लेषण

पैरामीटर अनुमान कठिनाई क्रम:

  1. μ (सबसे आसान): कम संख्या में घटनाओं से अनुमान लगा सकते हैं
  2. c, p (मध्यम): अधिक भूकंप घटनाओं की आवश्यकता, शांत कैटलॉग में अनुमान लगा सकते हैं
  3. α, K (सबसे कठिन): कम से कम एक मुख्य भूकंप और संबंधित आफ्टरशॉक की आवश्यकता; अधिक मुख्य भूकंप α के पश्च वितरण में काफी सुधार कर सकते हैं

प्रयोग निष्कर्ष

  1. पैरामीटर सहसंबंध: α-K और c-p में मजबूत सहसंबंध मौजूद है, संयुक्त विचार की आवश्यकता है
  2. डेटा आवश्यकता: विश्वसनीय पूर्वानुमान के लिए कम से कम एक पूर्ण मुख्य भूकंप-आफ्टरशॉक अनुक्रम की आवश्यकता है
  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. डेटा सीमा: सिंथेटिक कैटलॉग में केवल एकल मुख्य भूकंप शामिल है, निष्कर्षों की सामान्यता को सीमित करता है
  2. घटना अंतराल परिभाषा: केवल आसन्न घटना अंतराल पर विचार करता है, मुख्य भूकंप समूह के भीतरी संरचना पर विचार नहीं करता
  3. गणना चुनौती: मुख्य भूकंप संख्या बढ़ने पर मॉडल फिटिंग समय में काफी वृद्धि, अभिसरण कठिन
  4. विशेषता विश्लेषण अपर्याप्त: मुख्य भूकंप संख्या के अलावा, अन्य भूकंप विशेषताओं के पश्च वितरण पर प्रभाव को पूरी तरह से अन्वेषित नहीं किया गया है

भविष्य की दिशाएं

  1. पैरामीटर अपडेट रणनीति: वृद्धिशील शिक्षण विधि विकसित करना, पिछले मुख्य भूकंप के पश्च वितरण को अगले के पूर्व के रूप में उपयोग करना
  2. बहु-मुख्य भूकंप प्रयोग: अधिक मुख्य भूकंपों वाले सिंथेटिक कैटलॉग पर प्रयोग दोहराना
  3. घटना समूहीकरण विश्लेषण: मुख्य भूकंप समूह के भीतर और बीच की घटना अंतराल वितरण का अध्ययन करना
  4. विशेषता प्रभाव अनुसंधान: विभिन्न भूकंप विशेषताओं के मॉडल प्रदर्शन पर प्रभाव का व्यवस्थित विश्लेषण

गहन मूल्यांकन

शक्तियां

  1. विधि नवाचार: INLA विधि को ETAS मॉडल में सफलतापूर्वक लागू करना, व्यावहारिक गणना समस्या को हल करना
  2. पूर्ण प्रयोग डिजाइन: तीन प्रयोग विभिन्न कोणों से विधि प्रदर्शन का व्यापक मूल्यांकन करते हैं
  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 मॉडल की गणनात्मक दक्षता समस्या को सफलतापूर्वक हल करता है, और व्यवस्थित प्रयोगों के माध्यम से मूल्यवान पैरामीटर विश्लेषण और पूर्वानुमान मूल्यांकन प्रदान करता है। हालांकि सैद्धांतिक गहराई और प्रयोग व्यापकता में सुधार की गुंजाइश है, लेकिन इसकी विधि की व्यावहारिकता और निष्कर्षों की कार्यान्वयनशीलता इसे इस क्षेत्र का एक महत्वपूर्ण योगदान बनाती है।