Flow rate and leakage-channel geometry govern thermal anomalies in tunnel lining blocks
تُظهر هذه الدراسة، من خلال التجارب المختبرية والنمذجة ثلاثية الأبعاد، أن معدل تدفق المياه وهندسة قناة التسرب يحكمان بشكل حاسم شدة الشذوذ الحراري وتوزيعه المكاني وديناميكياته الزمنية في كتل بطانة النفق، مما يستوجب تفسيراً لبيانات درجة الحرارة يعتمد على الموقع المحدد.
تُعد الأنفاق الأرضية من عجائب الهندسة، لكنها تواجه عدواً مستمراً: الماء. فعندما يجد الماء الجوفي طريقاً عبر الجدران الخرسانية التي تحجز التربة، فإنه لا يكتفي بمجرد إحداث بقعة رطبة؛ إذ تعمل المياه المتحركة على سحب الحرارة بعيداً عن الصخور والخرسانة، مما يخلق بصمة باردة يمكن الشعور بها على السطح. وبالنسبة للمهندسين، فإن رصد انخفاضات درجات الحرارة هذه يعد وسيلة واعدة لاكتشاف التسريبات الخفية قبل أن تسبب أضراراً جسيمة. ومع ذلك، فإن قراءة هذه الإشارات الحرارية أمر معقد؛ إذ يعتمد شكل وحجم البقعة الباردة على سرعة حركة الماء والمسار الدقيق الذي يسلكه عبر الشقوق والفجوات خلف الجدار. فإذا تدفق الماء عبر شق مستقيم أو قناة متعرجة ومتحولة، فإن البقعة البردة الناتجة ستبدو مختلفة، مما يجعل من الصعب معرفة ما يحدث بالضبط تحت الأرض بمجرد النظر إلى خريطة حرارية.
ولحل هذه المعضلة، وضع فريق من الباحثين في جامعة شاندونغ نصب أعينهم فهم كيف تتحكم سرعة تسرب الماء وهندسة مسار التسريب في أنماط درجات الحرارة على جدار النفق. قاموا ببناء نموذج فيزيائي يحاكي مقطعاً من تبطين نفق، باستخدام كتلتين كبيرتين من الخرسانة موضوعتين فوق بعضهما البعض لتمثيل وصلة شائعة حيث تحدث التسريبات غالباً. وداخل هذا الإعداد، أنشأوا قنوات ضيقة مملوءة بالماء لمحاكاة المسارات المخفية التي قد يتسرب من خلالها الماء الجوفي. ومن خلال ضخ الماء بسرعات مختلفة وتعديل موقع وزاوية هذه القنوات، تمكنوا من مراقبة كيفية انتشار الماء البارد عبر الخرسانة. وقد قرنوا هذه التجارب الفيزيائية بمحاكاة حاسوبية دقيقة تتبع حركة الحرارة والماء، مما سمح لهم باختبار سيناريوهات أكثر بكثير مما قد يكون ممكناً في المختبر وحده.
وجد الباحثون أن سرعة الماء هي عامل مهيمن في مدى برودة جدار النفق. فعندما يتدفق الماء ببطء، يظل تأثير التبريد مركزاً بالقرب من الأعلى حيث يدخل الماء. ولكن مع زيادة معدل التدفق، ينتقل الماء البارد لمسافة أبعد في القناة وينشر برودته بشكل أعمق في الجزء السفلي من الهيكل الخرساني. وهذا يعني أن التسرب الأسرع يخلق منطقة باردة أكبر وأكثر انتشاراً. ومع ذلك، هناك حد لهذا التأثير؛ فبمجرد أن يتدفق الماء بسرعة كافية، تتوقف درجة الحرارة بالقرب من المدخل عن الانخفاض بشكل ملحوظ لأن الخرسانة هناك تصبح بالفعل ببرودة الماء نفسه. غير أن الاكتشاف الأكثر لفتاً للنظر جاء من مراقبة شكل مسار التسرب؛ فعندما قام الباحثون بإزاحة القناة السفلية جانبياً بالنسبة للقناة العلوية، لم تتحرك المنطقة الباردة في خط مستقيم ببساطة، بل أعاد التأثير التبريدي توزيع نفسه بطرق معقدة. فقد أصبحت بعض النقاط على الجدار أكثر برودة بكove مع تحرك القناة، بينما أصبحت نقاط أخرى أكثر دفئاً، مما يظهر أن موقع البقعة الباردة حساس للغاية للمحاذاة الدقيقة للتسرب.
وعلى النقيين من التغييرات الدراماتيكية الناتجة عن إزاحة القناة جانبياً، كان تدوير القناة قليماً له تأثير ضئيل جداً على نمط درجة الحرارة الملحوظ على السطح. إذ إن تدوير القناة السفلية بزاوية ثلاثين أو ستين درجة غير درجة الحرارة بجزء ضئيل جداً، وهو فرق صغير جداً يصعب اكتشافه دون أدوات دقيقة. وهذا يشير إلى أنه لأغراض اكتشاف التسريبات، فإن الموقع الأفقي لمسار التسرب يهم أكثر بكثير من ميله الزاوي الطفيف. كما كشفت الدراسة أن الوقت يلعب دوراً حاسماً في كيفية تطور هذه البقع الباردة؛ إذ بردت النقاط في الخرسانة الأقرب إلى قناة الماء بسرعة كبيرة، حيث انخفضت معظم درجات حرارتها في الساعات القليلة الأولى. أما النقاط الأبعد فقد بردت ببطء أكبر، واستمرت في انخفاض درجة حرارتها لفترة طويلة بعد التدفق الأولي. وهذا يعني أنه لا توجد لحظة واحدة يستقر فيها نمط درجة الحرارة؛ فالجدار في تغير مستمر مع انتشار البرودة.
تساعد هذه النتائج المهندسين على فهم أن الصورة الحرارية للنفق ليست مجرد لقطة بسيطة لتسرب واحد؛ فالبقعة الباردة المشاهدة على السطح هي نتيجة تفاعل ديناميكي بين سرعة تدفق الماء والمسار المحدد، والذي غالباً ما يكون غير منتظم، الذي يسلكه عبر الصخور والخرسانة. وبينما قدمت النماذج الحاسوبية والاختبارات المختبرية صورة واضحة لهذه الآليات، يشير الباحثون إلى أن تطبيق هذه المعرفة على الأنفاق في العالم الحقيقي سيتطلب مزيداً من التحقق. فقد افترضت النماذج وجود قناة مستقيمة مثالية وخصائص مواد ثابتة، في حين أن الأنفاق الحقيقية تحتوي على شقوق معقدة وفوضوية وظروف متغيرة. ومع ذلك، فإن هذا العمل يضع أساساً حيوياً: فمن أجل تشخيص التسرب بدقة، يجب على المرء النظر في كل من شدة التبريد والتوقيت المحدد لكيفية انتشار تلك البرودة، مع الإدراك بأن شكل البقعة الباردة يروي قصة عن الهندسة الخفية للتسرب نفسه.
ملخص تقني: معدل التدفق وهندسة قناة التسرب يحكمان الشذوذ الحراري في كتل تبطين الأنفاق
بيان المشكلة لا يزال تسرب المياه يمثل مشكلة مستمرة تتعلق بالصلاحية التشغيلية والمتانة في الأنفاق تحت الأرض، لا سيما في الأنفاق المحفورة بالمتفجرات (drill-and-blast)، حيث تكون مسارات التسرب غير متجانسة، وتتضمن شقوقاً منفصلة، وفواصل إنشائية، ومناطق تلامس غير مثالية. وبينما يتم اكتشاف التسرب غالباً عبر علامات الرطوبة السطحية، فإن هذه المؤشرات البصرية لا تكشف عن مورفولوجيا قناة التسرب أو شدة الرشح. وعلى الرغم من أن درجة الحرارة تعد مؤشراً واعداً لأن مياه التسرب تنقل الحرارة عبر الحمل الحراري (advection) بينما تستجيب البطانة عبر التوصيل (conduction)، إلا أن طرق الكشف الحالية تعتمد بشكل كبير على الصور؛ فهي تحدد الشذوذ الحراري ولكنها غالباً ما تفشل في تفسير كيف تتحكم بارامترات التسرب المحددة — مثل معدل التدفق، وإزاحة القناة، واتجاه القناة — في شكل وسعة والتطور الزمني لهذه البصمات الحرارية. هناك حاجة إلى تفسير قائم على الآليات يربط هيدروليكا التسرب وهندسة القناة بالتطور الزمني والمكاني لمجال درجة الحرارة.
المنهجية استخدمت هذه الدراسة نهجاً تجريبياً ورقمياً مشتركاً لاستقصاء الشذوذات الحرارية الناتجة عن التسرب في نظام تبطين نفق مبسط مكون من كتلتين.
الإعداد التجريبي: أُجريت اختبارات مخبرية باستخدام كتلتين من الخرسانة من فئة C50 (بأبعاد 1.0 م × 0.5 م × 0.5 م لكل منهما) تم رصهما فوق بعضهما لمحاكاة منطقة تبطين محلية. وقام صندوق مبطن بالفولاذ بمحاكاة بيئة مياه جوفية مضغوطة. تم اختبار ثلاثة نماذج مرجعية (النماذج A–C) مع تغيير معدلات التدفق (100–150 مل/دقيقة)، وإزاحات القناة (100–200 ملم)، وزوايا الدوران (0° و30°). وقامت خمسة وعشرون مستشعراً من نوع PT1000 مرتبة في مصفوفة 5×5 بمراقبة درجات الحرارة على مدار 6 ساعات بفواصل زمنية قدرها 10 ثوانٍ.
النمذجة الرقمية: تم تطوير نموذج ثلاثي الأبعاد لتدفق الكسور وانتقال الحرارة في برنامج COMSOL Multiphysics. ربط النموذج واجهة تدفق الكسور مع انتقال الحرارة في المواد الصلبة والسائلة. واستخدم صياغة "دارسي" للكسور بفتحة قدرها 0.5 ملم، وافترض خصائص مواد متجانسة. وقد حاكت النطاق الرقمي هندسة الكتلتين مع تكثيف مصفوفة المراقبة إلى 45 نقطة لالتقاط التدرجات المكانية الأكثر دقة.
التحقق والتحليل: تم التحقق من صحة النموذج الرقمي مقابل البيانات التجريبية لمدة 6 ساعات باستخدام جذر متوسط مربع الخطأ (RMSE)، ومتوسط الخطأ المطلق بالنسبة المئوية (MAPE)، ومعامل ارتباط بيرسون المربع (R2). عقب التحقق، أُجريت تحليلات رقمية أحادية العامل على مدار 24 ساعة لعزل تأثيرات معدل التدفق (10–200 مل/دقيقة)، وإزاحة القناة (0–400 ملم)، ودوران القناة (0°–60°) على كل من التوزيعات المكانية لدرجة الحرارة وتاريخ التبريد العابر.
النتائج الرئيسية كممت الدراسة كيف تتحكم شدة التسرب وهندسة القناة في الاستجابات الحرارية:
التحقق من النموذج: نجح النموذج الرقمي في التقاط التباينات المكانية الرئيسية لدرجات الحرارة واستجابات التبريد العابرة الملاحظة في التجارب. وعبر 30 سجلاً تاريخياً للمراقبة، تراوحت قيم RMSE بين 0.226 و0.952 درجة مئوية، وتراوحت قيم MAPE بين 1.126% و4.508%. ومع ذلك، تباين الارتباط الزمني (R2) بشكل كبير حسب الموقع، حيث تراوح بين أقل من 0.001 و0.994، مما يشير إلى أنه بينما كانت أخطاء المقدار ضئيلة، فإن توقيت الاستجابات كان يعتمد على الموقع.
تأثيرات معدل التدفق: أدى رفع معدل التدفق من 10 إلى 200 مل/دقيقة إلى تكثيف التبريد وتوسيع المدى المكاني لمنطقة درجة الحرارة المنخفضة على طول قناة التسرب. ومع ذلك، بالقرب من المدخل، أظهر انخفاض درجة الحرارة عوائد متناقصة مع اقتراب درجات الحرارة المحلية من درجة حرارة مياه المدخل (10 درجات مئوية).
تأثيرات إزاحة القناة: أدت زيادة الإزاحة الأفقية بين القناتين العلوية والسفلية إلى إعادة توزيع التبريد جانبياً. أنتج هذا استجابات حرارية تعتمد على الموقع وغير رتيبة في بعض الأحيان. فعلى سبيل المثال، عند نقاط مراقبة معينة، أدت زيادة الإزاحة في البداية إلى خفض درجات الحرارة قبل أن ترفعها مع ابتعاد منطقة التبريد. كان هذا التأثير أقوى بكثير من تأثير دوران القناة.
تأثيرات دوران القناة: كان تدوير القناة السفلية حول المحور الرأسي (0° إلى 60°) له تأثير ضعيف نسبياً على قسم المراقبة المركزي. بلغت أقصى فروق درجات الحرارة عبر زوايا الدوران حوالي 0.15 درجة مئوية، مما يشير إلى أن الدوران يعد مميزاً حرارياً أقل حساسية من معدل التدفق أو الإزاحة بالنسبة للهندسة المحددة وقسم المراقبة المدروس.
التطور العابر: لم تكن معدلات التبريد موحدة عبر قسم المراقبة. فالنقاط القريبة من القناة ومنطقة الاتصال بالواجهة استجابت بسرعة (مما شكل حوالي 95% من إجمالي تبريد الـ 24 ساعة خلال أول 6 ساعات)، بينما أظهرت النقاط الأكثر بعداً استجابات متأخرة. يشير هذا إلى عدم وجود زمن استقرار واحد لقسم المراقبة؛ حيث يتواجد التبريد السريع بالقرب من القناة جنباً إلى جنب مع التبريد المتأخر المدفوع بالتوصيل في أماكن أخرى.
الأهمية والقيود تزعم الورقة أن نتائجها تدعم التفسير المعتمد على الموقع للاستجابات الحرارية الداخلية، منتقلة من مجرد التعرف الظاهري نحو الفهم القائم على الآليات. ومن خلال إثبات أن معدل التدفق وإزاحة القناة يغيران بشكل كبير البصمة المكانية والتطور الزمني للشذوذات الحرارية، توفر الدراسة أساساً لربط التصوير الحراري السطحي بمعايير التسرب الخفية.
ومع ذلك، حافظ المؤلفون على موقف متواضع بشأن التطبيق المباشر لهذه النتائج:
تركز الدراسة على قسم مراقبة محلي لكتلتين من الخرسانة ولا تفترض أن نمط درجة الحرارة يحدد بشكل فريد حالة التسرب الكامنة.
غطت عملية التحقق التجريبي أول 6 ساعات فقط؛ والنتائج الخاصة بـ 24 ساعة هي تنبؤات لما بعد النافذة التجريبية ولا تثبت حالة التوازن التي تم التحقق منها تجريبياً.
تتطلب فروق درجات الحرارة الصغيرة الناتجة عن دوران القناة مزيداً من التدقيق الرقمي والقدرة على التحليل التجريبي لاعتبارها بصمات كمية قوية.
تتطلب قابلية نقل هذه النتائج إلى أنظمة الأنفاق كاملة النطاق ذات سماكات التبطين المتغيرة، ودرجات حرارة المياه الجوفية، والتهوية، واتصالات القنوات المعقدة مزيداً من التحقق.
في النهاية، تضع الدراسة إطاراً تم التحقق منه لمحاكاة الشذوذات الحرارية الناتجة عن التسرب، لكنها تؤكد أن التصوير الحراري السطحي والتقدير الكمي للتسرب بناءً على هذه الآليات يتطلبان مزيداً من التحقق.