A Variance-Reduction Framework for Practical Event-Based Monte Carlo PSHA

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

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

لعقود من الزمن، كانت الطريقة القياسية لتقدير مخاطر الزلازل هي محاكة خط زمني طويل للتاريخ، مليء بآلاف الزلازل العشوائية. في هذه المحاكاة، وفي كل مرة يتم فيها توليد زلزال، يقوم الكمبيوتر أيضًا باختيار رقم عشوائي لتمثيل مقدار اهتزاز الأرض. يمثل هذا الرقم العشوائي حقيقة أن زلزالين من نفس الحجم والمسافة يمكن أن ينتجا مستويات مختلفة من الاهتزاز. يتحقق الكمبيوتر بعد ذلك مما إذا كان ذلك الاهتزاز المولد عشوائيًا محددًا قويًا بما يكفي لتجاوز عتبة السلامة. إذا تجاوزها، يُحتسب كـ "إصابة"؛ وإذا لم يتجاوزها، يُحتسب كـ "إخفاق". ومن خلال تكرار هذه العملية ملايين المرات، يبني الكمبيوتر صورة عن مدى تكرار تجاوز العتبة. تعمل هذه الطريقة بشكل جيد عندما يكون الخطر مرتفعًا وتتجاوز العديد من الزلازل الحد، لكنها تعاني عند النظر في الأحداث الخطيرة والنادرة جدًا. في تلك الحالات، قد تؤدي الطبيعة العشوائية للاهتزاز إلى أن تفشل المحاكاة بالصدفة في رصد حدث خطير، أو تحسب حدثًا غير مرجح، مما يؤدي إلى نتائج متعرجة وغير موثوقة تبدو مثل الدرج بدلاً من منحنى سلس.

حدد عمل روباكييتي أن هذه العشوائية هي المسبب. فقد أدرك أنه بالنسبة للهدف المحدد المتمثل في حساب الاحتمالية الإجمالية للاهتزاز، فإن القيمة العشوائية الدقيقة لكل زلزال ليست ضرورية في الواقع. فبدلاً من اختيار رقم عشوائي ومعرفة ما إذا كان سيتخطى الخط، يمكن للكمبيوتر ببساطة حساب الاحتمالية الدقيقة بأن يتخطى الاهتزاز ذلك الخط لهذا السيناريو الزلزالي المحدد. هذا تغيير طفيف ولكنه قوي. الأمر يشبه محاولة تخمين متوسط طول حشد من الناس؛ إحدى الطرق هي قياس كل شخص وتوسيط الأرقام، أما الطريقة الأخرى، وهي ما يفعله هذا النهج الجديد، فهي النظر في الخصائص المعروفة للمجموعة وحساب المتوسط مباشرة دون قياس كل فرد. ومن خلال إجراء هذا الحساب لكل زلزال محاكى، يزيل الأسلوب "الضجيج" الناتج عن الاختيار العشوائي لقيم الاهتزاز. وتوضح الورقة أن هذا النهج، المعروف باسم "راو-بلاكليزيشن" (Rao-Blackwellization)، ينتج منحنى مخاطر أكثر سلاسة ودقة، خاصة لفترات العودة الطويلة التي يستخدمها المهندسون لتصميم البنية التحتية الحيوية.

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

للتأكد من أن هذا لم يكن مجرد تمرين نظري، طبق البحث الطريقة على نموذج واقعي لمصادر الزلازل في جنوب آيسلندا. هذه المنطقة معقدة، حيث توجد خطوط صدع متعددة وأنواع مختلفة من النشاط التكتوني. قام الباحثون بمحاكاة المخاطر لموقع صخري محدد في المنطقة، وقارنوا بين نهج أخذ العينات العشوائية التقليدي ونهج الاحتمال الشرطي الجديد. وقد عكست النتائج الاختبارات المحكومة؛ حيث أنتجت المحاكاة التقليدية منحنيات مخاطر تفاوتت بشكل كبير من جولة إلى أخرى، مما جعل من الصعب الوثوق بالأرقام للأحداث النادرة. بينما أنتج الأسلوب الجديد منحنيات ثابتة ومتسقة تطابقت مع حل مرجعي عالي الدقة. علاوة على ذلك، نظرت الدراسة في كيفية أداء هذه الأساليب عند حساب "طيف الخطر الموحد"، الذي يصف شدة الاهتزاز عبر ترددات اهتزاز مختلفة. غالبًا ما أنتجت الطريقة التقليدية أطيافًا متعرجة ومقطعة يصعب استخدامها هندسيًا، بينما حافظت الطريقة الجديدة على شكل سلس ومنطقي عبر جميع الترددات.

استكشفت الورقة أيضًا كيف يغير هذا التفكير الجديد الطريقة التي يحلل بها العلماء مصادر الخطر، وهي عملية تسمى "تفكيك التجميع" (deaggregation). في الطريقة التقليدية، ينظر العلماء إلى أي زلازل محددة تسببت فعليًا في تجاوز الاهتمام للحد في محاكاتهم العشوائية. ويمكن أن يكون هذا مضللاً إذا فاتت المحاكاة بعض الأحداث الرئيسية بالصدفة. أما الطالطريقة الجديدة، فهي تخصص وزنًا كسريًا لكل زلزال محاكى بناءً على احتمالية تسببه في خطر ما، بدلاً من مجرد عد الزلازل التي تجاوزت الخط بالفعل. وهذا يوفر صورة أوضح وأكثر استقرارًا حول سيناريوهات الزلازل التي تقود الخطر فعليًا. وفي مثال جنوب آيسلندا، قلل هذا النهج من عدم اليقين في تحديد أقوى أحجام المسافات الزلزالية بهامش كبير، مما أعطى المهندسين أساسًا أكثر موثوقية للتصميم.

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

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

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