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

المعلومات الأساسية

  • معرّف الورقة: 2510.13930
  • العنوان: التنبؤ بالزلازل باستخدام ETAS.inlabru
  • المؤلف: Ziwen Zhong
  • التصنيف: stat.AP (الإحصاء - التطبيقات)
  • تاريخ النشر: 15 أكتوبر 2025 (نسخة arXiv المسبقة)
  • رابط الورقة: https://arxiv.org/abs/2510.13930v1

الملخص

تدرس هذه الورقة استخدام نموذج ETAS للتنبؤ بالزلازل، وتقترح استخدام طريقة "inlabru" القائمة على INLA للاستدلال البايزي لحل مشاكل الطرق التقليدية MCMC المتعلقة بالحسابات المستهلكة للوقت والارتباط القوي بين المعاملات. تتحقق الورقة من فعالية الطريقة من خلال ثلاث تجارب: (1) استكشاف تأثير تثبيت المعاملات على التوزيعات اللاحقة للمعاملات الأخرى، مع اكتشاف علاقة تأثير متبادل واضحة بين معاملات α و K؛ (2) مقارنة توزيع الفترات الزمنية بين الأحداث المعيارية للبيانات الحقيقية والفهارس الاصطناعية؛ (3) استخدام الأحداث السابقة للهزة الرئيسية للتنبؤ بأحداث الزلازل خلال عشرة أسابيع بعد الهزة الرئيسية، مع استخدام اختبار الكمية ودرجة الاحتمالية المرتبة المستمرة (CRPS) لتقييم دقة التنبؤ.

الخلفية البحثية والدافع

تعريف المشكلة

تعتبر الزلازل من الكوارث الطبيعية التي تهدد حياة الإنسان وممتلكاته، وتتسم بطابع مفاجئ وصعوبة الهروب في الوقت المناسب. على الرغم من عدم إمكانية تحقيق التنبؤ بالهزات الرئيسية - وهو الهدف الأسمى للبحث - إلا أنه تم إحراز تقدم جيد في التنبؤ بالهزات الارتدادية. يمكن للتنبؤ الموثوق بالهزات الارتدادية أن يساعد الحكومات وفرق الإنقاذ على التخطيط بشكل أفضل لعمليات الإنقاذ بعد الزلزال.

الأهمية

ترتبط دقة التنبؤ بالهزات الارتدادية بشكل مباشر بـ:

  1. توفير فرص بقاء أكثر للأشخاص المتضررين
  2. حماية أفضل لحياة فرق الإنقاذ
  3. تحسين توزيع موارد الإنقاذ بعد الزلزال

قيود الطرق الحالية

  1. تقدير الحد الأقصى للاحتمالية (MLE): الطريقة مباشرة وبسيطة، لكنها لا تستطيع الحصول على عدم اليقين في تقدير المعاملات
  2. سلسلة ماركوف مونت كارلو (MCMC): كطريقة بايزية يمكنها الحصول على التوزيع اللاحق للمعاملات، لكنها تواجه المشاكل التالية:
    • استهلاك حسابي طويل
    • ارتباط عالي بين معاملات النموذج
    • صعوبة ضمان الحصول على توزيع لاحق جيد للمعاملات

الدافع البحثي

يمكن لطريقة 'inlabru' القائمة على التقريب اللابلاسي المتكامل المتداخل (INLA) أن:

  • تقلل بشكل كبير من وقت ملاءمة النموذج من خلال تقريب التوزيع اللاحق
  • تظهر أداء ممتازة عندما تكون العلاقة بين معاملات التغاير والمتغيرات الكامنة تقريباً خطية
  • تقدر التقريب الغاوسي داخلياً، مما يحسن الاستدلال البايزي

المساهمات الأساسية

  1. المساهمة المنهجية: تطبيق طريقة inlabru القائمة على INLA على نموذج ETAS للتنبؤ بالزلازل، مما يحل مشكلة الكفاءة الحسابية للطرق التقليدية MCMC
  2. اكتشاف العلاقات بين المعاملات: الكشف المنهجي عن العلاقات المتبادلة بين المعاملات في نموذج ETAS من خلال التجارب، خاصة الارتباط القوي بين معاملات α و K
  3. تحسين الكفاءة الحسابية: إثبات أن تثبيت معاملات α أو K إلى القيم الحقيقية يمكن أن يقلل وقت ملاءمة النموذج بأكثر من النصف
  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

تقريب لوغاريتم الاحتمالية

استخدام توسيع تايلور واستراتيجية التجميع لتقريب لوغاريتم الاحتمالية:

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

عملية الاستدلال البايزي

  1. استخدام المعاملات الأولية θ_0 للحصول على تقريب سلسلة تايلور لدالة لوغاريتم الاحتمالية
  2. إضافتها إلى لوغاريتم السابق، وحساب نمط التوزيع اللاحق θ_1*
  3. تحديث المعاملات من خلال البحث الخطي: θ* = αθ_0 + (1-α)θ_1*
  4. التكرار حتى التقارب (الفرق أقل من 1%)

نقاط الابتكار التقني

  1. تحسين استراتيجية التجميع: استخدام تجميع تكيفي، مع تجميع ضيق على الجانب الأيسر لزيادة الدقة وتجميع واسع على الجانب الأيمن لتقليل الحسابات
  2. تحويل الاحتمالية العكسية المتكاملة: تحويل عينات التوزيع الغاوسي إلى عينات التوزيع السابق المستهدف
  3. التقريب الخطي: تحقيق تقريب فعال لدالة الاحتمالية من خلال توسيع تايلور

إعداد التجارب

مجموعات البيانات

  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) (عملية بواسون موحدة)
  • زيادة 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
  • تجري بحثاً متعمقاً على أساس تطبيق inlabru من قبل Serafini et al. (2023) و Naylor et al. (2023)

المزايا النسبية

  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): تطبيق inlabru في عملية Hawkes
  • Naylor et al. (2023): النمذجة البايزية لنموذج ETAS

التقييم الشامل: هذه ورقة ذات قيمة عملية مهمة في مجال التنبؤ بالزلازل، وقد حلت بنجاح مشكلة الكفاءة الحسابية لنموذج ETAS، وقدمت من خلال التجارب المنهجية تحليلاً قيماً للمعاملات وتقييماً للتنبؤ. على الرغم من وجود مجال للتحسن في العمق النظري واتساع التجارب، فإن العملية العملية للطريقة وقابلية تشغيل الاستنتاجات تجعلها مساهمة مهمة في هذا المجال.