Monte Carlo Simulation of Earthquake Damage to Mexico City:A Cadastral Exposure Model with Voronoi-Weighted Spatio-Temporal Point Processes
تقدم هذه الورقة نموذج محاكاة مونت كارلو يدمج بيانات العقارات الكادسترالية، وتكتل المصادر الزلزالية بنموذج (ETAS)، وتضخيم الحركة الأرضية الخاص بكل موقع لتقدير الخسائر المالية المستقبلية الناجمة عن الزلازل في مدينة مكسيكو، مما يكشف عن متوسط خسارة سنوية قدرها 2.26 مليار بيزو مكسيكي وخسارة تبلغ 40 مليار بيزو مكسيكي تحدث مرة كل 100 عام.
المؤلفون الأصليون:Geronimo Rojas Barragan, Nelson Muriel Torrero
تقع مدينة مكسيكو فوق فخ جيولوجي. فقد بُنيت فوق القاع المجفف لبحيرة قديمة، وتستقر المدينة فوق طبقات من الطين الناعم المشبع بالمياه. عندما تنتقل الموجات الزلزالية عبر الصخور الصلبة للجبال المحيطة، فإنها تتحرك بسرعة وبأثر ضئيل. ولكن عندما تصطدم تلك الموجات نفسها بالرواسب البحيرية الناعمة تحت المدينة، فإنها تتباطأ وتتضخم في سعتها، تماماً مثل قطعة هلام (جيلي) تتذبذب فوق طبق. هذه الظاهرة، المعروفة باسم "تضخيم الموقع"، تحول الزلازل البعيدة والمتوسطة إلى كوارث محلية مدمرة. لقد شعرت المدينة بهذه القوة مرتين في العقود الأخيرة: أولاً في عام 1985، عندما تسبب زلزال هائل على بعد مئات الكيلومترات في انهيارات واسعة النطاق، ومرة أخرى في عام 2017، عندما كرر حدث أصغر وأقرب الدمار نفسه. ولأن حركة الأرض تسلك سلوكاً مختلفاً تماماً هنا عما هي عليه في أي مكان آخر، فإن الطرق التقليدية للتنبؤ بالأضرار الزلزالية تفشل. ولفهم المخاطر المالية الحقيقية التي تواجه هذه المدينة الضخمة، احتاج الباحثون إلى نموذج يحترم الفيزياء الفريدة لقاع البحيرة والقيمة المحددة للمباني القائمة فوقها.
قام فريق من الباحثين في جامعة "إيبيرواميريكانا" ببناء مثل هذا النموذج، حيث أنشأوا محاكاة تفصيلية لما قد تكلفه الزلازل المستقبلية مدينة مكسيكو. وبدلاً من الاعتماد على التقديرات الاقتصادية العامة أو أعداد السكان، استخدموا السجل العقاري الرسمي للمدينة، المعروف بالسجل الكاداسترالي. توفر هذه القاعدة البيانية الحكومية القيمة المقدرة لكل قطعة أرض والمنشآت المقامة عليها، مما يوفر خريطة دقيقة لما هو معرض للخطر بالضبط. وقد دمج الباحثون هذه الخريطة المالية مع نموذج حاسوبي متطور لكيفية حدوث الزلازل. حيث قاموا بتحليل عقود من السجلات الزلزالية لتحديد نوعين متميزين من المصادر الزلزالية التي تهدد المنطقة: الأحداث الضحلة بالقرب من الساحل حيث تنزلق الصفائح التكتونية بجانب بعضها البعض، والأحداث الأعمق التي تحدث داخل صفيحة مندسة بعيداً تحت السطح. ومن خلال التعامل مع هذه المصادر كأحداث منفصلة، يمكن للنموذج توليد آلاف السيناريوهات الزلزالية الاصطناعية، لكل منها موقع وزمن وقوة محددة.
يتضمن جوهر الدراسة تشغيل هذه السيناريوهات عبر نسخة افتراضية من المدينة. فمقابل كل زلزال محاكى، يحسب النموذج مقدار اهتزاز الأرض في موقع المدينة، مع مراعاة المسافة من مركز الزلزال والتضخيم الهائل الناتج عن رواسب البحيرة. ثم يترجم ذلك الاهتزاز إلى أضرار مادية لكل مبنى في السجل، محولاً الضرر الإنشائي إلى خسارة مالية. ولضمان أن يعكس النموذج الواقع، قام الباحثون بمعايرته بناءً على زلزال عام 2017. لم يكتفوا بمجرد تعديل الأرقام لتتطابق مع الفاتورة النهائية؛ بل اختبروا النموذج مقابل خمس حقائق مستقلة من ذلك الحدث: إجمالي الأموال المفقودة، وعدد المباني التي انهارت، وعدد المنشآت المتضررة بشدة، والمساحة الجغرافية التي غطاها الضرر، وعمر المباني التي سقطت. وقد أكد هذا الفحص الصارم أن النموذج يعيد إنتاج الأنماط المعقدة للدمار الواقعي بدقة.
تكشف نتائج هذه المحاكاة الضخمة، التي استمرت عبر مائتي مليون سنة من الوقت الافتراضي، عن نطاق واسع من النتائج المحتملة. ويُقدر متوسط الخسارة السنوية بنحو 2.26 مليار بيزو مكسيكي، لكن هذا المتوسط يخفي الطبيعة المتطرفة للمخاطر. ففي عام نموذجي، قد لا تشهد المدينة أي أضرار على الإطلاق. ومع ذلك، يتوقع النموذج أن حدثاً يحدث مرة واحدة في القرن قد يتسبب في أضرار بقيمة 40 مليار بيزو. إن عدم اليقين في هذه الأرقام كبير، وهو مدفوع إلى حد كبير بمدى تضخيم الأرض الناعمة للاهتزاز. ووجد الباحثون أن الخسارة المحتملة قد تكون منخفضة بقدر 14 مليار بيزو أو مرتفعة بقدر 439 مليار بيزو لحدث واحد بمقياس القرن، اعتماداً على السلوك المحدد للتربة. وهذا الفارق الواسع ليس فشلاً في النموذج، بل هو انعكاس للنقص الحالي في البيانات الدقيقة حول طبقات التربة في المدينة.
كما كشفت الدراسة عن طبقة خفية من الواقع الاقتصادي فيما يتعلق بكيفية الإبلاغ عن الأضرار. فعندما تمت معايرة النموذج لحدث عام 2017، أظهر أن الإحصاء الرسمي الحكومي للخسائر البالغ 16 مليار بيزو لم يسجل سوى جزء ضئيل من الضرر المادي الفعلي. ويشير النموذج إلى أن التكلفة الاقتصادية الحقيقية لذلك الزلزال كانت أقرب إلى 44 مليار بيزو، حيث يتمثل الفرق في الأضرار الإنشائية الطفيفة التي قام المالكون بإصلاحها دون الإبلاغ عنها للسلطات. وتوحي هذه النتيجة بأن الاعتماد فقط على التقارير الرسمية يقلل من تقدير العبء المالي الحقيقي للمخاطر الزلزالية. ويؤكد الباحثون أنه بينما توفر محاكاتهم إطاراً قوياً لفهم ضعف المدينة، فإن الخطوة التالية الأكثر أهمية هي الحصول على خرائط تربة مفصلة وقيم إنشائية رسمية لتضييق النطاق الواسع للخسائر المحتملة. وحتى ذلك الحين، تظل المحاكاة هي التقدير الأكثر دقة المتاح، وهي ترسم صورة لمدينة تظل فيها الأرض التي تحتها هي المتغير الأكثر قوة في سلامتها المستقبلية.
ملخص تقني: محاكاة مونت كارلو للأضرار الناجمة عن الزلازل في مدينة مكسيكو
بيان المشكلة تواجه مدينة مكسيكو هشاشة زلزالية استثنائية بسبب موقعها فوق رواسب بحيرية مشبعة، والتي تعمل على تضخيم الموجات الزلزالية طويلة المدى الناتجة عن أحداث التصدع البعيدة. وبينما أظهرت الأحداث التاريخية، مثل زلزال ميتشواكانان (Mw 8.0) عام 1985 وزلزال بويبلا-موريلوس (Mw 7.1) عام 2017، هذه الديناميكية، فإن تحديد الأضرار المالية المستقبلية يتطلب دمج الزلزالية الاحتمالية، وتوهين الحركة الأرضية، وجرد دقيق للأصول المعرضة للخطر. تعتمد النماذج الحالية غالبًا على مجاميع اقتصادية كلية أو مؤشرات سكانية إقليمية، مما يفتقر إلى التفصيل اللازم لنمذجة الكوارث بدقة في مدينة ضخمة معقدة. علاوة على ذلك، غالبًا ما تفشل طرق المعايرة القياسية في مراعاة التفاوت بين الضرر المادي والخسائر المبلغ عنها رسميًا، أو تفترض مصادر زلزالية متجانسة لا تعكس الأنظمة التكتونية المتنوعة في المكسيك.
المنهجية يبني المؤلفون إطارًا موحدًا يتكون من ثلاث طبقات ملائمة، ترتكز على السجل العقاري لمدينة مكسيكو (مستخرج عام 2021) كقاعدة للتعرض.
نموذج التعرض: تستخدم الدراسة السجل العقاري البلدي، الذي يغطي 2.29 مليون سجل. وإدراكًا منها بأن المستخرج العام يغفل قيم البناء، استنتج المؤلفون التعرض الهيكلي عبر تطبيق تكلفة استبدال ثابتة قدرها 12,000 بيزو مكسيكي/متر مربع على المساحة المبنية المسجلة، مما أدى إلى إجمالي تعرض هيكلي قدره 5.12 تريليون بيزو مكسيكي. ويؤكد النموذج أن جمع القيم على مستوى السجل قد وزع حصص الأراضي بشكل صحيح دون تكرار في العد.
الزلزالية والعملية النقطية: قام المؤلفون بتكييف نموذج تسلسل الهزات الارتدادية الوبائي (ETAS).
التقسيم المكاني: تم استخدام تقسيم "فورونوي" (Voronoi tessellation) لحقل مركز الزلزال في المكسيك، مع وزنه بمقلوب مساحة الخلية لتقدير الكثافة المكانية.
تجميع المصادر: تم تقسيم فهرس الزلازل (1980–2025، M ≥ 5.0) إلى مصدرين زلزاليين متميزين باستخدام تجميع الوسائل الكيمية (k-means clustering) بناءً على القدر والعمق البؤري (بدون توحيد المعايير). وهذا يفصل بين أحداث "الواجهة الضحلة" (متوسط العمق ~21 كم) وأحداث "الصفيحة الداخلية العميقة" (متوسط العمق ~109 كم).
المحاكاة الزمنية: تم نمذجة كل مصدر كعملية بواسون مكانية-زمانية ذات بنية كثافة ضربية، مما أنتج ناقلات زلزالية اصطناعية (الموقع، الوقت، القدر) على مدار 200,000 قرن (أي 2 × 10⁷ سنة محاكاة).
الحركة الأرضية والهشاشة:
التوهين: يتم التنبؤ بالحركة الأرضية باستخدام علاقات توهين خاصة بكل نظام (منفصلة عند عمق 40 كم) مرتكزة إلى سجلات الأرض الصلبة المرصودة من حدثي 1985 و2017. ويتم تقريب تضخيم الموقع باستخدام بصمة "المنطقة صفر" (Zona Cero) من أطلس المخاطر لعام 2017، عبر تخصيص عوامل قدرها 4.8 (للترسبات البحيرية) و1.5 (للمناطق الانتقالية) لمراعاة تأثيرات الحوض.
الهشاشة: يقوم نموذج حالة الضرر المنفصل بتحويل الحركة الأرضية (PGA) إلى نسب خسارة مالية. وبخلاف المنحنيات اللوغاريتمية المستمرة، يستخدم هذا النموذج أربع حالات ضرر (طفيفة، متوسطة، شديدة، كاملة) مع نسب عواقب محددة.
المعايرة: تمت معايرة بارامترات الهشاشة عبر تقدير الاحتمال الأقصى المشترك لخمس ملاحظات مستقلة من زلزال 2017 (Mw 7.1): الخسارة المالية المبلغ عنها، عدد المباني المنهارة، عدد المباني المتضررة بشدة، مساحة البصمة المتضررة، وحصة الانهيارات في المباني التي بنيت قبل عام 1985. ومن الأهمية بمكان أن النموذج يتضمن بارامتر اكتمال الإبلاغ (ρ) للتوفيق بين الضرر المادي والإحصاءات الرسمية، مع الاعتراف بأن الأضرار الطفيفة غالبًا ما يتم عدم الإبلاغ عنها.
المساهمات الرئيسية
قاعدة تعرض دقيقة: أرست الورقة أول نموذج تعرض عالي الدقة لمدينة مكسيكو بناءً على السجل العقاري، متجاوزةً مؤشرات السكان إلى تقييم الأصول على مستوى المربعات السكنية.
تفكيك المصادر: من خلال تجميع المصادر الزلزالية بناءً على العمق والقدر دون توحيد المعايير، يلتقط النموذج السلوك الفيزيائي المتميز لأحداث الواجهة والصفيحة الداخلية، والتي تفرض أنظمة حركة أرضية مختلفة.
المعايرة المشتركة مع انحياز الإبلاغ: ترفض الدراسة ثوابت القياس البسيطة لصالح ملاءمة إحصائية مشتركة تنمذج صراحةً نقص الإبلاغ عن الأضرار الطفيفة إلى المتوسطة. وهذا يسمح للنموذج بإعادة إنتاج البصمة وعدد الانهيارات المرصودة مع تقدير إجمالي الخسائر الاقتصادية، والتي وُجد أنها أعلى بكماً من التقارير الرسمية.
الكفاءة الحسابية: يستخدم محرك المحاكاة أخذ العينات الدقيق لـ (inverse-CDF) لتوابع ETAS وتقليل إزاحة القدر للحركة الأرضية، مما يسمح بتنفيذ 200,000 واقعة طولها قرن لكل حد من حدود البارامتر في جزء من الملي ثانية لكل قرن.
النتائج تقدم محاكاة مونت كارلو مقاييس المخاطر التالية لمدينة مكسيكو:
متوسط الخسارة السنوية (AAL): التقدير المركزي هو 2.26 مليار بيزو مكسيكي (النطاق: 0.14–27.0 مليار بيزو مكسيكي)، وهو ما يمثل حوالي 0.044% من التعرض الهيكلي.
فترات العودة: يُقدر إجمالي الخسارة السنوية لفترة 1 في 100 عام بـ 40.25 مليار بيزو مكسيكي (النطاق: 1.4–439 مليار بيزو مكسيكي).
مخاطر الذيل: تصل القيمة المعرضة للخطر (VaR) عند مستوى 99.9% إلى 301 مليار بيزو مكسيكي (مركزياً)، بينما تصل قيمة الذيل المعرضة للخطر (TVaR) عند نفس المستوى إلى 560 مليار بيزو مكسيكي، مما يشير إلى ذيل ثقيل حيث تهيمن الأحداث المتطرفة على توزيع الخسائر.
اكتمال الإبلاغ: تشير المعايرة إلى أن 11.7% فقط من الأضرار الطفيفة و63.6% من الأضرار المتوسطة دخلت في الإحصاء الرسمي لعام 2017، مما يعني أن إجمالي الخسارة الاقتصادية المادية يبلغ حوالي 44.4 مليار بيزو مكسيكي (مقابل 16 مليار بيزو مكسيكي المبلغ عنها).
التوقعات العشرية: على مدى القرن القادم، تُقدر الخسارة التراكمية المتوسطة بـ 226 مليار بيزو مكسيكي (مركزياً)، رغم أن متوسط العقد الواحد يسجل خسائر أقل بكثير (حوالي 5 مليارات بيزو مكسيكي) بسبب هيمنة الأحداث الكبيرة والنادرة.
الأهمية والقيود تفترض الورقة أن نطاق عدم اليقين الواسع (الذي يمتد عبر عامل قدره ~77 بين الحدود الدنيا والعليا) ليس عيباً في العرض، بل هو انعكاس صادق لمحدودية البيانات، مدفوعاً أساساً بنقص بيانات التمنطق الدقيق (microzonation). ويؤكد المؤلفون أن "النطاق هو النتيجة"، جادلين بأن الاكتفاء بذكر القيمة المركزية فقط سيكون مضللاً.
تُقر الدراسة صراحةً بعدة قيود:
مقياس الشدة: يُنظر إلى استخدام تسارع الأرض الذروي (PGA) على أنه غير مثالي لمدينة مكسيكو، حيث ينتج الضرر عن الرنين الطيفي في المباني متوسطة الارتفاع؛ ومع ذلك، فإن بيانات دورة المباني اللازمة غير متوفرة في السجل العقاري.
التعرض الساكن: يفترض التوقع "مدينة مجمدة" بقيم أصول وخصائص بناء من عام 2021، متجاهلاً التوسع الحضري وتحسينات كود البناء المستقبلية، مما قد يميل بالتوقعات نحو التشاؤم.
تثبيت التوهين: تم تثبيت نموذج التوهين على ملاحظتين فقط بدلاً من إجراء تحليل انحدار كامل، مما يؤدي إلى نقص في التقدير عند المسافات الكبيرة جداً (مثل حدث 2017 Mw 8.2).
الاستقلال المكاني: يتم إعادة أخذ عينات مواقع الهزات الارتدادية بشكل مستقل عن حقل الهزة الرئيسية بدلاً من تجميعها مكانياً، وإن كان المؤلفون يشيرون إلى أن هذا له تأثير ضئيل على خسائر مدينة مكسيكو نظراً لبعد جميع المصادر.
يخلص المؤلفون إلى أن الخطوات التالية الأكثر قيمة تتضمن الحصول على بيانات التمنطق الدقيق (NTC-Sismo)، وترجمة جداول قيم البناء الرسمية لاستبدال نموذج تكلفة الاستبدال، وتأمين تقديرات مستقلة لاكتمال الإبلاغ لتثبيت بارامتر ρ خارجياً.