Sensitivity Analysis of Shield Tunnel-Induced Ground Settlement to Stratigraphic Spatial Variability: From Random Field Simulation to Causal Decoupling
تقدم هذه الورقة إطاراً متكاملاً يجمع بين محاكاة المجال العشوائي، والتجارب العددية المنضبطة، والانحدار الرمزي لتحديد كيفية تأثير التباين المكاني الجيولوجي في الجرانيت المتجوّي على الهبوط الناجم عن أنفاق الدروع، مع تحديد معامل التشوه للطبقات المحفورة والطبقات المغطية كعوامل مهيمنة، واستخلاص نموذج تنبؤي قابل للتفسير فيزيائياً.
تخيل الأرض تحت أقدامنا ليس ككتلة صخرية صلبة وثابتة، بل ككعكة طبقات ضخمة وفوضوية. أحياناً تكون طبقة الكريمة سميكة، وأحياناً رقيقة؛ وأحياناً تكون الإسفنجة ناعمة ولينة، وفي أحيان أخرى تكون صلبة كالقالب. هذا هو واقع الجيولوجيا، خاصة عندما نحاول حفر أنفاق للمترو من خلالها. لعقود من الزمن، حاول المهندسون التنبؤ بمدى هبوط الأرض (الاستقرار) عندما تلتهم آلة حفر الأنفاق العملاقة (TBM) طريقها عبرها. كانوا في السابق يخمنون بناءً على معادلات بسيطة أو ينظرون إلى الأرض كما لو كانت متجانسة تماماً، مثل نكهة واحدة من الآيس كريم. لكننا نعلم الآن أن الأرض مليئة بالمفاجآت؛ فهي تمتلك "تباينًا مكانيًا"، وهو تعبير معقد يعني أن خصائص التربة تتغير من مكان إلى آخر، حتى عبر مسافات قصيرة. إذا حفرت في بقعة ما، قد تكون الأرض صلبة؛ وإذا تحركت بضعة أمتار، قد تصبح لينة. هذا عدم القدرة على التنبؤ يمثل صداعاً كبيراً للمهندسين، لأنه إذا أخطأوا في التقدير، فقد تهبط الأرض بشكل كبير، مما يؤدي إلى تشقق الأرصفة والمباني فوقها. لذا، السؤال الكبير هو: كيف نتنبأ بهذا الهبوط عندما تكون الأرض فوضوية ومتغيرة إلى هذا الحد؟
تغوص هذه الورقة البحثية في ذلك الواقع الفوضوي، وتنظر تحديداً في الأنفاق الدرعية (من النوع الذي يستخدم رأس قطع دوار) التي تحفر عبر الجرانيت المتجوّي في شنتشن، الصين. أراد الباحثون معرفة مدى تأثير "الفوضوية" الطبيعية للأرض على الهبوط بدقة، والأهم من ذلك، تحديد أي جزء من الأرض هو الأكثر أهمية. لم يكتفوا بمجرد التخمين؛ بل بنوا توأماً رقمياً متطوراً للنفق باستخدام برنامج حاسوبي قوي يسمى PLAXIS 3D. لقد قاموا بمحاكاة الأرض كـ "حقل عشوائي"، وهو ما يشبه رمي النرد لإنشاء آلاف النسخ المختلفة والواقعية للطبقات الموجودة تحت الأرض، حيث تمتلك كل نسخة قوة وصلابة مختلفة قليلاً. ووجدوا أنه بمجرد تغيير التوزيع الطبيعي للصخور والتربة، يمكن أن يتفاوت مقدار هبوط الأرض بنحو 3.38 ملم. قد يبدو هذا الرقم ضئيلاً، لكنه في عالم حفر الأنفاق يعد فرقاً هائلاً قد يعني الفرق بين نفق آمن ومبنى متصدع.
الجزء الأكثر إثارة في اكتشافهم هو كيف عرفوا لماذا تهبط الأرض. لقد عاملوا الأرض كفريق مكون من جزئين: "الطبقة المحفورة" (الطبقة التي يقطعها النفق بالفعل) و"الغطاء الصخري" (طبقات التراب والصخر التي تقع فوق النفق). أجروا 37 تجربة حاسوبية مختلفة، مع تغيير شيء واحد في كل مرة، تماماً مثل عالم يضبط مقادير وصفة ما. ووجدوا أن الصلابة (مدى قساوة أو ليونة المادة) للطبقة التي يقطعها النفق هي "الرئيس المطلق". فإذا كانت تلك الطبقة أكثر صلابة، يهبط سطح الأرض بشكل أقل. كما أن صلابة الطبقات التي تعلوها مهمة أيضاً، لكنها تأتي في المرتبة الثانية من حيث الأهمية. ومن المثير للدهشة، وجدوا أن عوامل أخرى، مثل مدى "تماسك" التربة (التماسك) أو زاوية انزلاقها (زاوية الاحتكاك)، لا تشكل فرقاً يُذكر على الإطلاق. حتى العدد الدقيق للطبقات أو كيفية تكدسها لم يكن مهماً طالما أن الصلابة الإجمالية هي نفسها. الأمر يشبه إدراك أنك عندما تدفع صندوقاً ثقيلاً، فإن وزن الصندوق يهم أكثر بكثير من لون السجادة التي يستقر عليها.
ولجعل هذا مفيداً للمهندسين في العالم الحقيقي، استخدم الفريق خدعة رياضية ذكية تسمى "الانحدار الرمزي". فكر في الأمر كإنسان آلي ينظر إلى جميع بياناتهم الحاسوبية ويحاول كتابة جملة بسيطة ومفهومة (معادلة) تشرح النتائج، بدلاً من مجرد إعطاء إجابة مبهمة (صندوق أسود). لقد توصلوا إلى صيغة تستخدم ثلاثة أرقام فقط — صلابة طبقة النفق، وصلابة الطبقات العلوية، وتماسك طبقة النفق — للتنبؤ بمدى هبوط الأرض. كانت هذه الصيغة دقيقة بشكل مذهل، حيث أعطت تنبؤات صحيحة بفارق يقل عن نصف مليمتر في المتوسط. كما استخدموا أداة تسمى SHAP للتحقق من عملهم، والتي أكدت بالفعل أن صلابة طبقة النفق هي بالفعل "نجم العرض" الحقيقي.
تجادل الورقة البحثية أساساً ضد الفكرة القديمة التي تقول إننا بحاجة لمعرفة كل تفصيل صغير ودقيق عن الأرض لتقديم تنبؤ جيد. بدلاً من ذلك، يقترحون أنه إذا ركزنا على الحصول على "الصلابة" الصحيحة للطبقات المناسبة، فيمكننا تقديم تنبؤات أفضل بكما دون الضياع في التفاصيل. لقد استبعدوا صراحةً فكرة أن زاوية الاحتكاك أو العدد الدقيق للطبقات هي عوامل رئيسية في هذا النوع المحدد من أرض الجرانيت. وبينما تعتمد نتائجهم على عمليات المحاكاة الحاسوبية وليس الحفر الفعلي في العالم الحقيقي، فإن الاتساق بين تجاربهم المنضبطة وصيغتهم الجديدة يمنحهم ثقة عالية في نتائجهم. هم لا يدعون أنهم حلوا كل مشاكل الأنفاق في الكون، لكنهم قدموا طريقة أكثر وضوحاً وبساطة ودقة لفهم كيف يؤثر التباين الطبيعي للأرض على سلامة الأنفاق، مما يساعد المهندسين على بناء مترو أكثر أماناً في المستقبل.
ملخص تقني: تحليل الحساسية لهبوط الأرض الناجم عن أنفاق الدروع تجاه التباين المكاني للطبقات الجيولوجية
بيان المشكلة تؤدي عمليات حفر الأنفاق باستخدام درع الحماية (Shield Tunnelling) في البيئات الجيولوجية المعقدة، ولا سيما في تكوينات الجرانيت المتجوية التي تتميز بواجهات طبقية غير منتظمة وتباينات حادة في الصلابة، إلى إدخال عدم يقين كبير في تنبؤات هبوط الأرض. غالبًا ما تفترض التحليلات الحتمية التقليدية وجود تربة متجانسة أو تعتمد على القيم المتوسطة، مما يؤدي إلى الفشل في رصد الطبيعة المزدوجة لعدم اليقين الجيولوجي: وهي التباين المنفصل في التتابعات الطبقية، والتباين المكاني المستمر للمعاملات الجيوتقنية داخل الطبقات. علاوة على ذلك، فإن تحليلات الحساسية القائمة غالبًا ما تعامل الملف التربي بأكمله كنظام موحد، متجاهلة التمييز الميكانيكي الجوهري بين طبقة الحفر (مصدر التشوه) والطبقات التي تعلوها (مسار انتقال التشوه). تعالج هذه الدراسة الحاجة إلى قياس وتفسير والتنبؤ بعدم اليقين في الهبوط الناتج عن هذه التباينات المكانية، وتحديدًا في سياق سكة حديد شنتشن البينية.
المنهجية طوّر المؤلفون إطار عمل متكامل يجمع بين محاكاة الحقل العشوائي، والتجارب العددية المنضبطة، والانحدار الرمزي:
النمذجة العددية والتحقق من الصحة: تم إنشاء نموذج ثلاثي الأبعاد باستخدام برنامج PLAXIS 3D، يحاكي نفق درع يعمل بضغط التربة المتوازن (EPB) بقطر 8.0 م وعمق دفن متوسط يبلغ 20 م. استخدم النموذج نموذج "موهر-كولوم" (Mohr–Coulomb) لخصائص التربة ونماذج مرنة خطية للمكونات الإنشائية. وقد تم التحقق من صحة النموذج مقابل بيانات الرصد الميداني وطريقة "بيك" (Peck) التجريبية، مما أظهر دقة عالية في إعادة إنتاج أخاديد الهبوط السطحي في الجرانيت المتجوي.
إطار عمل اقتران عدم اليقين المزدوج: لمحاكاة التباين المكاني الجيولوجي، استخدمت الدراسة نهجًا يتكون من وحدتين:
محاكاة الطبقات: استُخدمت سلاسل ماركوف المقترنة (CMC) لتوليد 11 واقعة جيولوجية مستقلة، مما ساهم في رصد عدم اليقين المنفصل في التتابعات الطبقية (مثل تذبذب واجهات التجوية).
تعيين المعاملات: طُبقت الحقول العشوائية الشرطية (CRF) مع دوال الارتباط الذاتي الأسية ضمن كل واقعة لتخصيص قيم متغيرة مكانيًا لمعامل التشوه (E)، والتماسك (c)، وزاوية الاحتكاك (ϕ)، مما رصد التباين البارامتري المستمر.
التجارب البارامترية المنضبطة: أُجريت سلسلة من 37 تجربة عددية بنهج "تغيير متغير واحد في كل مرة" (OVAT) لفصل آثار المعاملات الجيولوجية. درست الدراسة بشكل مستقل تأثيرات E و c و ϕ لكل من طبقة الحفر والطبقات التي تعلوها، بالإضافة إلى عدد الطبقات، لعزل مساهماتها السببية في أقصى هبوط سطحي (Smax).
الانحدار الرمزي والتفسير: تم توسيع مجموعة البيانات باستخدام أخذ العينات بالهيبيركوب اللاتيني (LHS) لتغطية فضاء المعاملات بشكل أكثر تجانسًا. ثم استُخدم الانحدار الرمزي لاشتقاق معادلة تنبؤ صريحة وقابلة للتفسير. وأخيرًا، طُبق تحليل SHAP (تفسيرات شابلي المضافة) على النموذج الناتج للتحقق من ترتيب الحساسية وضمان التفسير الفيزيائي.
النتائج الرئيسية
تأثير التباين المكاني: أنتجت الـ 11 واقعة للحقول العشوائية أقصى هبوط سطحي تراوح بين $-11.25ملمو-14.63ملم،بتباينقدره3.38$ ملم. وهذا يؤكد أن التباين المكاني الجيولوجي وحده يمكن أن يسبب عدم يقين في الهبوط ذا أهمية عملية، وهو ما قد تقلل النماذج الحتمية القائمة على القيم المتوسطة من تقديره.
ترتيب الحساسية: أرست التجارب المنضبطة تسلسلًا هرميًا واضحًا للحساسية: Eexc (معامل تشوه طبقة الحفر) >Eupper (معامل تشوه الطبقات العليا) >cexc (تماسك طبقة الحفر).
تسببت التغيرات في Eexc (35–75 ميجا باسكال) في تغيير الهبوط بمقدار $3.72$ ملم.
تسببت التغيرات في Eupper (5–30 ميجا باسكال) في تغيير قدره $2.57$ ملم.
تسببت التغيرات في cexc (25–48 كيلوباسكال) في تغيير قدره $1.09$ ملم.
كان لزاوية الاحتكاك (ϕ) وتعقيد الطبقات تأثير ضئيل، حيث كانت تغيرات الهبوط أقل من $0.18$ ملم.
الأدوار الميكانيكية: حددت الدراسة أن طبقة الحفر تعمل كمصدر للتشوه، حيث تتحكم كل من الصلابة والقوة (إلى حد أقل) في توليد التشوه. وفي المقابل، تعمل الطبقات العليا أساسًا كمسار لانتقال التشوه، حيث تعد الصلابة هي العامل الحاكم الوحيد؛ ولم تظهر معاملات القوة (c و ϕ) في الطبقات العليا أي تأثير ملموس على الهبوط السطحي.
النموذج التنبئي: حققت معادلة الانحدار الرمزي المستمدة، والتي تتضمن Eexc و Eupper و cexc، معامل تحديد R2 قدره $0.806،ومتوسطخطأمطلق(MAE)قدره0 .265ملم،وجذرمتوسطمربعالخطأ(RMSE)قدره0.454ملم.وقدأعادتحليلSHAPعلىهذاالنموذجإنتاجترتيبالحساسية(E_{exc} \gg E_{upper} > c_{exc}$)، مما أكد اتساقه الفيزيائي.
الأهمية والادعاءات يدعي البحث أنه يطور تحليل هبوط أنفاق الدروع من مجرد توصيف بسيط لعدم اليقين إلى نهج قائم على الآلية وفصل الأسباب. وتتمثل مساهماته الرئيسية في:
الابتكار المنهجي: تطوير إطار عمل لعدم اليقين المزدوج الذي ينمذج في آن واحد عدم اليقين في التتابع الطبقي (عبر CMC) والتباين البارامتري داخل الطبقة (عبر CRF)، مما يقدم تمثيلًا أكثر واقعية لتكوينات الجرانيت المتجوي مقارنة بالنماذج الحتمية أو النماذج العشوائية ذات المصدر الواحد.
الرؤية الميكانيكية: كشف الفصل الصريح بين طبقة الحفر والطبقات العليا أن لكل منهما أدوارًا ميكانيكية مختلفة جوهريًا. وهذا يتحدى تحليلات الحساسية العالمية التقليدية التي تعامل الملف التربي كوحدة واحدة، مما يوفر أساسًا أكثر دقة من الناحية الفيزيائية لتقييم المعاملات.
المنفعة العملية: تثبت الدراسة أن معامل التشوه هو المعامل المهيمن الذي يتحكم في الهبوط، وهو أكثر تأثيرًا بكثير من التماسك أو زاوية الاحتكاك. كما تقترح معادلة انحدار رمزي صريحة وقابلة للتفسير تعمل كبديل شفاف لنماذج التعلم الآلي "الصندوق الأسود". وتسهل هذه الأداة التنبؤ الأولي السريع بالهبوط، وفحص المعاملات الجيولوجية، وتقييم المخاطر دون التضحية بالتفسير الفيزيائي.
يضع المؤلفون هذا العمل كخطوة نحو النمذجة التنبؤية القائمة على الفيزياء في الهندسة الجيوتقنية، مما يسد الفجوة بين المحاكاة العشوائية المعقدة واتخاذ القرارات الهندسية العملية.