Simulation water movement through the GCL cover using a physical model and HYDRUS-1D software
تقيم هذه الدراسة أداء برنامج HYDRUS-1D في محاكاة حركة المياه عبر تكوينات مختلفة لمبطنات الطين الجيوسنثيتية (GCL)، حيث وجدت أنه بينما يلتقط النموذج اتجاهات التدفق بدقة وكفاءة عالية، فإنه يقلل بشكل منهجي من أحجام المياه المطلقة بسبب فشله في مراعاة خاصية كره الماء لطبقات المنسوجات الأرضية.
تخيل الأرض تحت أقدامنا كأنها إسفنجة عملاقة غير مرئية. لعقود من الزمن، اعتقد الناس أن هذه الإسفنجة بارعة جداً في تنظيف الانسكابات الفوضوية، لدرجة أنه إذا دفنا النفايات في مكب للنفايات، فإن أي سائل يتسرب منها (يسمى الرشاحة) سيتم تنقيته ببساطة أثناء تغلغله في التربة. لكن العلماء أدركوا في النهاية أن هذه "الإسفنجة" ليست سحرية؛ فإذا كانت النفايات شديدة السمية أو كان التسرب سريعاً جداً، فإن الإسفنجة ستتعرض للضغط وتفشل، ويمكن للمياه الملوثة أن تسمم المياه الجوفية التي نشربها. ولإيقاف ذلك، اخترع المهندسون "معطف مطر" خاصاً لمكبات النفايات يسمى المبطن الطيني الاصطناعي (GCL). فكر في الـ GCL كأنه شطيرة عالية التقنية: طبقتان من القماش (مثل شبكة متينة) تحصر بينهما طبقة سميكة من الطين فائق الامتصاص. عندما يبتل هذا الطين، ينتفخ مثل الإسفنجة، مما يؤدي إلى إحكام إغلاقه لمنع مرور الماء. ولكن الجزء الصعب هنا هو: كيف نعرف بالضبط كمية الماء التي ستمر عبر هذه الشطيرة قبل أن نبني مكباً ضخماً للنفايات؟ لا يمكننا الانتظار لنرى النتيجة، لأن الضرر قد يكون قد وقع بالفعل. وهنا يأتي دور النماذج الحاسوبية؛ حيث يستخدم العلماء البرمجيات لإنشاء "توأم رقمي" للمكب، يحاكي كيفية تحرك الماء عبر هذه المبطنات للتنبؤ بما إذا كانت ستصمد أم لا.
في هذه الدراسة، قرر فريق من الباحثين وضع هذا التوأم الرقمي قيد الاختبار. فقد قاموا ببناء تجربة فيزيائية تتضمن اثني عشر نسخة مختلفة من هذه الشطائر الطينية، مع تغيير المكونات: استخدموا نوعين من الطين (النوع A والنوع I)، ونوعين من القماش (PP و PET)، وثلاث كميات مختلفة من وزن الطين (4، 4.5، و5 كجم لكل متر مربع). قاموا بتمرير الماء عبر هذه العينات في المختبر، وقاسوا بدقة كمية الماء التي تسربت من خلالها. ثم أدخلوا هذه البيانات في برنامج حاسوبي شهير يسمى HYDRUS-1D ليروا ما إذا كان البرنامج يمكنه التنبؤ بالنتائج ببراعة تضاهي البشر.
جاءت النتائج مزيجاً من "عمل رائع" و"شبه مثالي". كان النموذج الحاسوبي بارعاً للغاية في تخمين نمط حركة الماء. فعندما قارن الباحثون تنبؤات الكمبيوتر بالقياسات المخبرية الحقيقية، كان التطابق قوياً، حيث أصاب النموذج حوالي 97% من اتجاهات الترشيح (تغلغل الماء) و81% من اتجاهات التصريف. ومع ذلك، لم يكن الكمبيوتر مثالياً في تخمين الكمية الدقيقة؛ فقد توقع باستمرار أن كمية أقل قليلاً من الماء ستمر مما حدث بالفعل في المختبر. وجد الباحثون أن النموذج قلل من تقدير حجم الماء بهامش صغير ولكنه ملحوظ. ويعتقدون أن هذا حدث لأن الكمبيوتر يعامل طبقات القماش كتربة عادية محبة للماء (hydrophilic)، بينما في الواقع، الأقمشة الاصطناعية في الـ GCL طاردة للماء (hydrophobic). الأمر يشبه أن الكمبيوتر يعتقد أن القماش عبارة عن منشفة عطشى تمتص الماء فوراً، بينما القماش الحقيقي يشبه سترة واقية من المطر تقاوم الماء للحظة قبل أن تسمح له بالمرور.
على الرغم من هذا "الخلل" الصغير في أرقام الحجم، أكدت الدراسة أن النموذج الحاسوبي أداة قوية للمهندسين. فقد أثبت أن نوع الطين مهم جداً: المبطن المصنوع من الطين النوع A (الذي يحتوي على المزيد من معدن المونتموريلونيت) كان أفضل بكثير في حجب الماء من النوع I. ومن المثير للدهشة أن نوع القماش المستخدم (PP مقابل PET) لم يشكل فرقاً كبيراً في مدى قدرة المبطن على حجب الماء. لذا، بينما يحتاج النموذج الحاسوبي إلى تعديل بسيط ليأخذ في الاعتبار الطبيعة "المقاومة للماء" للقماش، فقد نجح في إظهار أن هذه المبطنات الطينية هي حواجز فعالة، وأنه يمكننا الاعتماد على المحاكاة للمساعدة في تصميم مكبات نفايات أكثر أماناً دون الحاجة لبناء عشرات النماذج الأولية الفيزيائية لكل مشروع على حدة.
ملخص تقني: محاكاة حركة المياه عبر أغطية البطانات الطينية الجيوسنثيتيكية (GCL) باستخدام نموذج فيزيائي وبرنامج HYDRUS-1D
بيان المشكلة لا تزال تلوث المياه الجوفية الناتج عن تسرب الرشاحة في مواقع طمر النفايات يمثل تحدياً بيئياً حرجاً. وبينما كانت الآراء التقليدية تشير سابقاً إلى أن طبقات التربة يمكنها معالجة الرشاحة طبيعياً، فقد أظهرت الأبحوث اللاحقة أن مدافن النفايات غالباً ما تشكل مخاطر جسيمة على موارد المياه الجوفية. وللتخفيف من ذلك، تم اعتماد البطانات الطينية الجيوسنثيتيكية (GCL) على نطاق واسع كحواجز هيدروليكية فعالة. ومع ذلك، فإن أداء هذه البطانات يعتمد على تفاعلات معقدة بين أنواع البنتونيت، وطبقات الجيوتكستايل، وتوزيع الكتلة. وعلاوة على ذلك، بينما تُستخدم النماذج العددية مثل HYDRUS-1D بشكل شائع لمحاكاة تدفق المياه، فإن قدرتها على التنبؤ بدقة بحركة المياه عبر الطبقات الهيدروفوبية (الكارهة للماء) وغير المتجانسة الخاصة بـ GCL تتطلب تحققاً صارماً مقابل البيانات التجريبية الفيزيائية.
المنهجية اعتمدت الدراسة نهجاً مزدوجاً يجمع بين التجريب الفيزيائي والمحاكاة العددية:
1. النمذجة الفيزيائية:
العينات: تم اختبار اثنتي عشرة تهيئة مختلفة لـ GCL، تتكون من نوعين من صوديوم بنتونيت (النوع A والنوع I)، ونوعين من الجيوتكستايل (البولي بروبلين [PP] والبوليستر [PET])، وثلاث كتلات بنتونيت لكل وحدة مساحة (4، 4.5، و5 كجم/م²).
التوصيف: تم تحليل عينات البنتونيت لتحديد التركيب الكيميائي (XRD، XRF، قياس اللهب، المعايرة) ومؤشر الانتفاخ الحر.
اختبار النفاذية: استُخدم جهاز مخصص لاختبار النفاذية لقياس الموصلية الهيدروليكية تحت ضغط ثابت قدره 30 كيلو باسكال (محاكاة لظروف قاع مدفن النفايات) لمدة 6 ساعات. أُجريت الاختبارات في ثلاث مكررات لكل تهيئة.
جمع البيانات: تم تسجيل أحجام الارتشاح التراكمي ومياه التصريف عند فترات زمنية مختلفة.
2. المحاكاة العددية (HYDRUS-1D):
إعداد النموذج: تم إنشاء ملف تربة أحادي البعد (1D) ليتطابق تماماً مع سمك عينات GCL.
المعايرة: استُخدمت خوارزمية تحسين "ليفنبرج-ماركوارت" لمعايرة معاملات فان جينين-موليم الهيدروليكية (α,n,Ks,θr,θs) باستخدام بيانات المكرر التجريبي الأول. تم تثبيت معامل اتصال المسام (l) عند 0.5.
التحقق: طُبقت المعلمات المعايرة على نموذج مباشر (أمامي) لمحاكاة التدفق للمكرر الثاني والثالث دون أي تعديل إضافي.
مقاييس الأداء: تم تقييم دقة النموذج باستخدام متوسط الخطأ المطلق (MAE)، وجذر متوسط مربع الخطأ (RMSE)، ومعامل التحديد (R2)، وكفاءة النمذجة (EF)، ومعامل كتلة البواقي (CRM).
النتائج الرئيسية
الموصلية الهيدروليكية: كشف التحليل الإحصائي (ANOVA) أن نوع البنتونيت وكتلته لكل وحدة مساحة يؤثران بشكل كبير على الموصلية الهيدروليكية المشبعة (p<0.01)، بينما لم يكن لنوع الجيوتكستايل (PP مقابل PET) أي تأثير معنوي.
النوع A مقابل النوع I: خلافاً للتوقعات الأولية بناءً على المحتوى المعدني، أظهرت العينات ذات بنتونيت النوع A (التي تحتوي على نسبة أعلى من المونت موريلونيت ومؤشر انتفاخ حر أعلى) موصلية هيدروليكية أقل (أداء أفضل) من عينات النوع I. على سبيل المثال، أظهرت عينة PA5 (النوع A، 5 كجم/م²) أدنى موصلية (0.05 مم/دقيقة)، بينما أظهرت عينة TI4 (النوع I، 4 كجم/م²) أعلى موصلية (0.18 مم/دقيقة).
أداء النموذج:
الارتباط: كان هناك ارتباط قوي بين الاتجاهات المحاكات والمقاسة. تراوحت قيم R2 للارتشاح بين 0.98 و0.99، وللتصريف بين 0.9973 و0.9998.
الكفاءة: بلغت كفاءات النمذجة (EF) نحو 0.81 للتصريف و0.97 للارتشاح.
تحليل الخطأ: أظهر النموذج نقصاً منهجياً في تقدير أحجام المياه المطلقة. بلغ جذر متوسط مربع الخطأ (RMSE) 1.45 للتصريف و4.71 للارتشاح. وكان معامل كتلة البواقي (CRM) سالباً (-0.06 للتصريف، -0.26 للارتشاح)، مما يؤكد هذا الميل للتقليل من التقدير.
الأهمية والادعاءات تؤكد الورقة البحثية أنه بينما يوفر برنامج HYDRUS-1D تقديراً قوياً لديناميكيات التدفق الإجمالية وهو فعال للغاية في محاكاة حركة المياه في GCL، إلا أن لديه قيوداً محددة فيما يتعلق بالخصائص الفيزيائية للمواد الاصطناعية.
قيد الهيدروفوبية (كراهية الماء): يعزو المؤلفون التقليل المنهجي في تقدير أحجام المياه إلى اعتماد النموذج على دوال هيدروليكية تقليدية للتربة (Mualem-van Genuchten)، والتي تفترض وجود وسط مسامي محب للماء. فالنموذج لا يأخذ في الحسبان بشكل كامل الطبيعة الهيدروفوبية لطبقات الجيوتكستايل (PP و PET)، والتي تقاوم البلل الأولي حتى يتم الوصول إلى ضغط الاختراق.
الآثار العملية:
اختيار المواد: خلصت الدراسة إلى أن الخصائص المعدنية للبنتونيت (تحديداً النوع A) هي العامل المهيمن على فعالية الحاجز، متفوقة على الاختلافات الطفيفة في الكتلة لكل وحدة مساحة.
استقلالية الجيوتكستايل: بما أن نوع الجيوتكستايل لم يظهر أي تأثير معنوي على الموصلية الهيدروليكية، يمكن للمهندسين اختيار الجيوتكستايل بناءً على القدرة على البقاء الميكانيكي، أو المقاومة الكيميائية، أو التكلفة دون المساس بالسلامة الهيدروليكية.
أداة فحص: رغم التباينات الطفيفة الناتجة عن الهيدروفوبية، فإن الكفاءة الإحصائية العالية تثبت أن HYDRUS-1D يعد أداة فحص أولية قوية؛ حيث يمكنه تقليل الوقت والتكاليف المالية المرتبطة بالاختبارات الفيزيائية الشاملة لأجهزة النفاذية في المشاريع الجيوتقنية واسعة النطاق بشكل كبير.
الخلاصة نجحت الدراسة في التحقق من استخدام HYDRUS-1D لمحاكاة GCL مع التأكيد على ضرورة تفسير النتائج مع فهم الخصائص الهيدروفوبية للمواد. وقد أثبتت أن النوع A من البنتونيت هو الخيار الأفضل لبطانات مدافن النفايات من بين التكوينات المختبرة، وأكدت إمكانية فصل اختيار الجيوتكستايل عن مخاوف الأداء الهيدروليكي.