GEO-HYBRID: A Physics-Informed Python Framework for Liquefaction Potential Assessment and Integrated Geotechnical Interpretation
تقدم هذه الورقة GEO-HYBRID، وهو إطار عمل بلغة بايثون يتسم بالشفافية ويعزز تقييم مخاطر التميع والتفسير الجيوتقني من خلال دمج مبادئ ميكانيكا التربة الأساسية مع معايرة المعلمات القائمة على البيانات، وبذلك يتفوق على نماذج تعلم الآلة البحتة في دقة التنبؤ والموثوقية.
إن الأرض تحت أقدامنا ليست كتلة صخرية صلبة ومتجانسة، بل هي مزيج فوضوي من الحبيبات المعدنية الدقيقة والماء والهواء. هذه المادة، المعروفة باسم التربة، تتصرف بطرق غالباً ما تتحدى التنبؤ البسيط. فالمنحدر الذي يصمد بثبات لسنوات قد ينهار فجأة بعد هطول أمطار غزيرة واحدة، أو قد تهبط قاعدة بناء بشكل غير متوقع عندما تهتز الأرض أثناء زلزال. وتكمن الصعوبة الجوهرية في حقيقة أن خصائص التربة تتغير بشكل جذري من مكان إلى آخر، وأحياناً حتى ضمن متر واحد. إذ يمكن لتغير طفيف في كمية المياه التي تحتفظ بها التربة أن يغير قوتها بهامش كبير، محولاً تلة مستقرة إلى كتلة منزلقة. ولأن الأرض غير قابلة للتنبؤ، فقد اعتمد المهندسون لفترة طويلة على الاختبارات الفيزيائية التي تُجرى مباشرة في الموقع لفهم ما يكمن تحت السطح، بدلاً من الثقة في الحسابات النظرية وحدها.
لعقود من الزمن، اعتمد النهج المعياري للتنبؤ بفشل التربة على صيغ رياضية تعامل التربة كما لو كانت مادة متسقة وقابلة للتنبؤ. وبينما تعمل هذه الصيغ جيداً من الناحية النظرية، إلا أنها غالباً ما تفشل في الممارسة العملية لأن العالم الحقيقي فوضوي. وعندما يحاول المهندسون التنبؤ بما إذا كانت التربة ستتحول إلى حالة تشبه السائل أثناء الزلزال — وهي ظاهرة خطيرة تسمى "تسيّل التربة" (liquefaction) — فإن الطرق التقليدية تخطئ في تقديرها كثيراً، وأحياناً بفارق شاسع. وفي السنوات الأخيرة، اتجه الكثيرون إلى الذكاء الاصطناعي، آملين أن تتمكن النماذج الحاسوبية من تعلم أنماط سلوك التربة من البيانات الماضية. ومع ذلك، غالباً ما تواجه هذه النماذج الحاسوبية صعوبة عند مواجهة مواقف جديدة، حيث تقوم أساساً بحفظ الأمثلة القديمة دون فهم حقيقي للفيزياء الكامنة وراءها.
تقدم دراسة جديدة مساراً مختلفاً للمضي قدماً، مساراً يجسّر الفجوة بين القوانين الفيزيائية الراسخة وبيانات العالم الحقيقي دون الاعتماد على حيل حاسوبية غامضة. فقد طور الباحثون إطار عمل يسمى "GEO-HYBRID"، والذي يأخذ المعادلات الموثوقة التي حكمت الهندسة الجيوتقنية لأجيال ويقوم بمعايرتها باستخدام سجلات حالات فعلية من الميدان. وبدلاً من ترك الكمبيوتر يبتكر قواعد خفية لتناسب البيانات، يقوم هذا النهج بتعديل الأرقام المحددة ذات المعنى فقط داخل القوانين الفيزيائية المعروفة. وقد اختبر الفريق هذه الطريقة على أربع مجموعات متميزة من بيانات العالم الحقيقي، بما في ذلك مئات المواقع التاريخية للزلازل حيث حدث تسيّل للتربة، وآلاف قياسات العمق من اختبارات الاختراق المخروطية، ونتائج مختبرية لعينات من الطين والتربة من فنلندا والنمسا.
تظهر النتائج أنه عندما يتم تغذية المعادلات الفيزيائية بالبيانات الصحيحة والمحققة، فإنها تؤدي عملها بدقة ملحوظة. وفي اختبار رئيسي شمل 251 موقعاً تاريخياً للزلازل، اضطر الباحثون أولاً إلى تصحيح خطأ طويل الأمد في كيفية قراءة البيانات؛ حيث اكتشفوا أن عموداً رئيسياً في قاعدة البيانات، كان يُعتقد أنه يمثل نسبة الإجهاد، كان في الواقع يمثل عمق مستوى المياه الجوفية. وبمجرد إصلاح ذلك، قفزت دقة النموذج بشكل كبير. ومن خلال إضافة ثلاث قواعد محددة ومبررة فيزيائياً إلى النموذج — مثل إدراك أن الرمال شديدة الكثافة تحت الاهتزاز الشديد تسلك سلوكاً مختلفاً، أو تصفية الحالات التي كانت فيها أدلة تسيّل التربة غامضة — نجح النظام في التنبؤ بشكل صحيح بما إذا كانت التربة ستتسيّل في 90.1% من الحالات. ويمثل هذا تحسناً كبيراً مقارنة بالطرق السابقة.
في المقابل، عندما قام الباحثون بتدريب نموذج تعلم آلي قياسي على نفس البيانات تماماً، بدا وكأنه يتعلم بشكل مثالي، محققاً نسبة نجاح بلغت 92.8% في مجموعة بيانات التدريب. ومع ذلك، عندما تم اختباره على بيانات جديدة وغير مرئية، انهارت دقة أدائه لتصل إلى 76.3%، مما أثبت أنه كان يحفظ الأمثلة القديمة فحسب بدلاً من تعلم قواعد الأرض. وبالمقارنة، حافظ النموذج القائم على الفيزياء على دقته العالية في البيانات الجديدة. ويشير هذا إلى أن المعادلات نفسها ليست هي المشكلة، بل كانت المشكلة دائماً في المعايير (parameters) التي يتم تغذيتها بها. كما أظهرت الدراسة أن التغيرات الصغيرة في محتوى الرطوبة يمكن أن تقلل بشكل كبير من سلامة المنحدر، مما يحدد كمياً خطراً كان المهندسون يعرفونه نوعياً ولكن لم يكن بإمكانهم حسابه بسهولة.
كما تسلط الدراسة الضوء على أهمية معرفة ما يجب عدم تضمينه. فقد اختبر الباحثون خاصية تتعلق بمدى سرعة تصريف المياه عبر التربة، متوقعين أن تؤدي إلى تحسين التنبؤات، ولكن بدلاً من ذلك، أدى إضافة هذه الخاصية إلى جعل النموذج أسوأ، لأن البيانات الموجودة تحتوي بالفعل وبشكل ضمني على معلومات حول التصريف. ومن خلال توثيق هذا الفشل، تمنع الدراسة الباحثين المستقبليين من إضاعة الوقت في عوامل مكررة. ويخلص العمل إلى أن الطريقة الأكثر موثوقية للتنبؤ بسلوك التربة ليست استبدال الفهم البشري بخوارزميات "الصندوق الأسود"، بل صقل فهمنا للثوابت الفيزيائية باستخدام أدلة من العالم الحقيقي. يوفر إطار العمل طريقة شفافة وقابلة للتكرار لتحويل بيانات التربة غير المؤكدة إلى قرارات هندسية موثوقة، مما يضمن أننا عندما نبني على الأرض، فإننا نبني على أساس من الحقيقة المتحقق منها وليس على مجرد التخمين.
ملخص تقني: إطار عمل GEO-HYBRID للتنبؤ بالهندسة الجيوتقنية
بيان المشكلة
تواجه الهندسة الجيوتقنية تحديًا جوهريًا: التربة والصخور هي مواد متعددة الأطوار، منفصلة، وتتسم بتباين مكاني عالٍ واستجابات ميكانيكية غير متوقعة. غالبًا ما تفشل النماذج التأسيسية الحتمية الكلاسيكية (مثل نموذج مور-كولوم، وكام-كلاي) في التنبؤ بظواهر مثل تميع التربة، والتعرية الداخلية، والانزلاقات الأرضية الناجمة عن الأمطار بدقة مقبولة. إذ تظهر التنبؤات النظرية المعزولة أخطاءً تتجاوز 50% بشكل روتيني، ويمكن لتقلبات طفيفة في الظروف البيئية (مثل تغير بنسبة 2% في المحتوى المائي) أن تغير عوامل أمان استقرار المنحدرات بنسبة تتراوح بين 10 إلى 30%.
وبينما يعتمد المجال بشدة على الاختبارات الموقعية (SPT، CPT) والطرق شبه التجريبية (مثل منحنيات بولانجر-إدريس)، لا تزال هذه النهج تسجل معدلات خطأ في التصنيف المتبقي تتراوح بين 10-15%. وفي المقابل، فإن نماذج تعلم الآلة (ML)، رغم قدرتها على تحقيق دقة تدريب عالية، تعاني غالبًا من عدم الشفافية، وتتطلب مجموعات بيانات ضخمة نادرة التوفر في الجيوتقنية، وتميل إلى الإفراط في التخصيص (overfitting) لمجموعات الحالات التاريخية الصغيرة وغير المتوازنة، مما يؤدي إلى ضعف التعميم على البيانات غير المرئية. إن المشكلة الجوهرية المحددة ليست في المعادلات الفيزيائية الحاكمة، بل في تقدير المعاملات والافتقار إلى المعايرة الصارمة والقابلة للتفسير فيزيائيًا مقابل البيانات الواقعية.
المنهجية
يقدم البحث إطار عمل GEO-HYBRID، وهو إطار عمل بلغة بايثون (Python) قائم على الفيزياء ومعتمد على المعايرة، صُمم لربط المعادلات النظرية الموثقة بالبيانات الواقعية دون الاعتماد على الأوزان الضربيه الغامضة. تتبع المنهجية بروتوكولًا مكونًا من ست مراحل:
تفكيك النطاق: يتم رسم خرائط للتنبؤات المستهدفة (عامل أمان التميع، قدرة التحمل، الانضغاط، زاوية الاحتكاك) إلى المعادلات الحاكمة الموثقة (مثل بولانجر-إدريس، تيرزاجي-فيسيك، تايلور، فريدلوند).
حصر المعاملات: تُصنف المدخلات كمتغيرات قابلة للتعديل، أو ثابتة، أو مشتقة، مع فحص صارم للوحدات.
بناء المعادلات: يتم تنفيذ المعادلات الفرعية كدوال مستقلة ومحققة الوحدات، ثم تُسلسل لتشكيل تدفقات عمل كاملة (على سبيل المثال: ملف تعريف CPT ← الإجهادات ←qc1Ncs← CRR ← FS).
التحقق من الصحة: يتم التحقق من المعادلات الفرعية مقابل الجداول والمنحنيات المنشورة، محققة أخطاء تتراوح بين 0-2%.
المعايرة: يتم ضبط الثوابت ذات المعنى الفيزيائي (الميول، التقاطعات، المئويات) مقابل بيانات حقيقية باستخدام بروتوكولات صارمة خارج عينة التدريب، بما في ذلك استبعاد 70/30 الطبقي، والتحقق من صحة "ترك حدث واحد خارجًا" (LOEO) لحالات التميع.
اكتشاف الميزات: تُضاف الميزات الإضافية فقط إذا كانت تمثل آليات فيزيائية موثقة. ويتم الإبلاغ صراحة عن النتائج السلبية (الميزات المرفوضة).
يُطبق مبدأ مركزي هنا: لا توجد أوزان غامضة. كل معامل هو إما ثابت أدبي ذو معنى فيزيائي أو ثابت تمت معايرته مع تفسير فيزيائي صريح.
المساهمات الرئيسية
يقدم البحث سبع مساهمات رئيسية:
إطار عمل بايثون قابل لإعادة الاستخدام: نظام معياري يربط المعادلات الموثقة، ويعاير الثوابت مقابل البيانات الواقعية، ويتحقق من الصحة عبر بروتوكولات صارمة خارج عينة التدريب.
إعادة بناء قاعدة البيانات: أول إعادة بناء عامة وقابلة للقراءة آليًا لقاعدة بيانات بولانجر-إدريس (2014) الخاصة بـ CPT للتميع (253 حالة)، مع التحقق من هويات الأعمدة مقابل التقرير الرسمي (UCD/CGM-14/01).
تصحيح حساب CSR: تحديد وتصحيح خطأ في معالجة البيانات حيث خزن الجدول المنشور عمق المياه الجوفية والإجهاد الفعال بدلًا من نسبة الإجهاد. هذا التصحيح وحده حسن الدقة بمقدار +5.6 نقطة مئوية.
ميزات الحدود الفيزيائية الجديدة: إدخال ثلاث ميزات مبررة فيزيائيًا:
غطاء رملي شديد الكثافة (لـ qc1Ncs>130، FC<5%، amax/g>0.5).
غطاء شديد القوة (لـ Mw>8.5 مع هزات منخفضة).
مرشح مظاهر لوما بريتا (استبعاد الحالات التي تفتقر إلى أدلة سطحية موثوقة).
توثيق النتائج السلبية: الرفض الصريح لميزة تصريف قائمة على النفاذية، مما يثبت أن قاعدة البيانات تشفر بالفعل تأثيرات التصريف.
مقارنة تعلم الآلة بالفيزياء: عرض موثق يوضح أن نموذج XGBoost الذي تم تدريبه على نفس البيانات قام بحفظ مجموعة التدريب (دقة 92.8%) ولكنه انهار عند مواجهة بيانات غير مرئية (76.3%)، مما جعله أقل أداءً من النموذج الفيزيائي المعاير.
الارتباطات الإقليمية: ارتباطات معايرة وذات معنى فيزيائي لكل من معامل الانضغاط (طين فنلندا) وزاوية الاحتكاك (تربة النمسا)، بما في ذلك ميزة حجم الحبيبات (d50) التي تؤكد فيزياء الاحتكاك بين الجزيئات.
النتائج
تم اختبار إطار العمل على أربع قواعد بيانات مستقلة:
التميع (251 حالة): حقق النموذج الرياضي الصرف دقة أولية بلغت 79.7%. بعد تصحيح حساب CSR، وإضافة ميزات الحدود الثلاث، وتطبيق مرشح لوما بريتا، ارتفعت الدقة إلى 90.1% (89.3% مع LOEO). في المقابل، حقق منحنى روبرتسون-وايد (1998) دقة 72.1%، بينما انخفض نموذج XGBoost إلى 76.3% في مجموعة الاختبار.
تنفيذ CPT (2,765 عمقًا): طابق التنفيذ الداخلي لإجراء بولانجر-إدريس تنفيذ liquepy/CLiq بمتوسط خطأ نسبي قدره 0.04%، مما يؤكد أن المعادلات نفسها ليست مصدر الخطأ عند معرفة المعاملات.
طين فنلندا (114 عينة): حقق نموذج معامل الانضغاط المعاير (Cc مقابل المحتوى المائي) متوسط خطأ مطلق (MAE) قدره 0.514، متفوقًا بشكل كبير على ارتباط تيرزاجي-بيك (MAE 1.055). وأشار الميل المعاير إلى أن طين فنلندا أكثر انضغاطًا بنحو 2.7 مرة لكل وحدة محتوى مائي مما يشير إليه الارتباط الكلاسيكي.
تربة النمسا (39 زوجًا): قلل نموذج زاوية الاحتكاك الذي يتضمن اللدونة وحجم الحبيبات (d50) من متوسط الخطأ المطلق (MAE) في مجموعة الاختبار إلى 2.52°، مقارنة بـ 13.03° للعلاقات القياسية من نوع كيني. وأكدت إشارة معامل حجم الحبيبات أن الحبيبات الأكبر تزيد من الاحتكاك.
الأهمية والادعاءات
يزعم البحث أن المصدر المهيمن للخطأ في التنبؤات الجيوتقنية هو تقدير المعاملات، وليس المعادلات الحاكمة. وتكمن أهمية إطار عمل GEO-HYBRID في قدرته على:
استعادة التفسير الفيزيائي: من خلال معايرة الثوابت بمعانٍ فيزيائية بدلاً من الأوزان الغامضة، يظل النموذج قابلًا للتفسير والنقل عبر المواقع المختلفة.
التفوق على تعلم الآلة في البيانات الصغيرة: يتفوق النموذج الفيزيائي المعاير في التعميم على نماذج تعلم الآلة عالية القدرة في مجموعات البيانات الجيوتقنية الصغيرة وغير المتوازنة، متجنبًا فخ "الحفظ".
تحديد كمية عدم اليقين: يدمج إطار العمل حساسية الرطوبة وتحليل العزم الأول من الدرجة الثانية (FOSM) لتحويل التنبؤات النقطية إلى مؤشرات موثوقية (β) واحتمالات الفشل (Pf)، مع الإقرار بالعشوائية غير القابلة للاختزال في سلوك التربة.
إرساء القابلية للتكرار: يتم توفير جميع البيانات، والمحللات (parsers)، والشيفرات البرمجية كموارد مفتوحة بلغة بايثون مع بذور (seeds) ثابتة، مما يضمن إمكانية إعادة توليد كل رقم تم تسجيله.
يخلص المؤلف إلى أنه بينما تظل بعض الأخطاء المتبقية بسبب التشتت الفيزيائي غير القابل للاختزال للتربة وتحيزات تفسير قواعد البيانات (مثل حالات لوما بريتا)، فإن الانضباط المقترح المتمثل في "معادلات موثقة + ثوابت معايرة + تحقق صارم" يوفر مسارًا قويًا للتنبؤ الجيوتقني.