A tri-level optimization framework for the resilient design of water distribution networks under pipe failures
تقترح هذه الورقة إطار عمل للتحسين ثلاثي المستويات يتسم بالكفاءة الحسابية، يجمع بين الاسترخاءات الخطية وخوارزمية توليد الأعمدة والقيود المتداخلة لتصميم شبكات توزيع مياه مرنة قادرة على الحفاظ على الإمداد في ظل أسوأ حالات فشل الأنابيب.
تُعد المياه شريان الحياة للمدن الحديثة، حيث تتدفق عبر متاهة واسعة وخفية من الأنابيب التي توصل مياه الشرب الآمنة إلى الصنابير، وأدوات الاستحمام، والنوافير. ولعقود من الزمن، صمم المهندسون هذه الشبكات للتعامل مع ضغط الطلب اليومي، لضمان بقاء الضغط مستقرًا والإمداد مستمرًا. ومع ذلك، فإن واقع البنية التحتية المتقادمة وعدم القدرة على التنبؤ المتزايدة بالطقس القاسي يعني أن الأنابيب ستتعطل حتمًا. وعندما ينفجر أنبوب رئيسي، لا يكمن التحدي المباشر في إصلاح الكسر فحسب، بل في الحفاظ على تدفق المياه إلى أكبر عدد ممكن من الناس مع احتواء الضرر. وتُعرف قدرة النظام على مواصلة العمل رغم الصدمة باسم "المرونة". وبينما استخدم المهندسون منذ فترة طويلة عمليات المحاكاة الحاسوبية للتحقق من كيفية تعامل الشبكة مع كسر أنبوب واحد، إلا أن هذه الفحوصات غالبًا ما تُجرى بعد الانتهاء من التصميم، أو تعتمد على تخمينات مبسطة حول سلوك المياه. ويبقى السؤال: كيف يمكننا تصميم شبكة منذ الب البداية تكون مضمونة رياضيًا لتقديم أفضل أداء لها في أسوأ سيناريو ممكن؟
لقد تصدى فريق من الباحثين في "إمبريال كوليدج لندن" لهذا التحدي من خلال ابتكار طريقة جديدة لتخطيط شبكات المياه تتعامل مع الفشل كجزء لا يتجزأ من عملية التصميم. فبدلاً من الأمل في الأفضل، بنوا إطار عمل يتخيل أسوأ سيناريو ثم يصمم نظامًا للنجاة منه. لقد تناولوا المشكلة كـ "لعبة استراتيجية" تتكون من ثلاث حركات. أولاً، يقرر "المدافع"، الذي يمثل مخططو المدن، أين يتم تركيب أنابيب أو صمامات جديدة لتعزيز الشبكة. ثانيًا، يختار "المهاجم"، الذي يمثل قوى الطبيعة، الأنبوب الأكثر ضررًا ليقوم بكسره. وأخيرًا، يحصل المدافع على حركة ثانية: يجب عليه فورًا ضبط الصمامات وأدوات التحكم لإعادة توجيه المياه والحفاظ على إمداد أكبر عدد ممكن من المستهلكين. والهدف هو إيجال التصميم الأولي الذي يسمح للنظام بتوفير أكبر كمية من المياه، حتى بعد أن يختار المهاجم الفشل الأكثر تدميرًا وبعد أن يبذل المشغلون قصارى جهدهم للاستجابة.
تكمن الصعوبة في فيزياء المياه؛ فالمياه لا تتدفق مثل الكهرباء في سلك بسيط، بل تتبع قواعد غير خطية معقدة حيث تتغير الضغوط والسرعات بطرق لا يمكن التنبؤ بها، خاصة عند إغلاق الأنابيب أو فتحها. وتواجه الطرق الحاسوبية التقليدية صعوبة في حل هذه اللعبة ثلاثية الخطوات لأن الرياضيات تصبح ثقيلة جدًا للحساب في شبكة مياه لمدينة حقيقية. وقد طور الباحثون حلاً ذكيًا؛ حيث ابتكروا نسخة خطية مبسطة من فيزياء المياه يسهل على الحواسيب حلها بسرعة، مما يعطيهم تقديرًا جيدًا جدًا للنتيجة المثلى. ثم استخدموا خوارزمية متخصصة تضيف بشكل تكراري سيناريوهات الفشل الأكثر حرجًا إلى المزيج، لتقوم بتنقيح التصميم حتى يصبح مثاليًا تقريبًا. ولضمان عدم تضليلهم من قبل رياضياتهم المبسطة، أجروا فحصًا نهائيًا مفصلاً لتصاميمهم المقترحة باستخدام معادلات الفيزياء الكاملة والمعقدة لمعرفة كمية المياه التي سيتم تسليمها بالفعل.
وعندما طبقوا هذه الطريقة على شبكة مياه واسعة وحقيقية في شرق إنجلترا، كانت النتائج مذهلة. تخدم الشبكة التي اختبرها الباحثون أكثر من ألف نقطة ربط، وهي مقسمة حاليًا إلى مناطق منفصلة للمساعدة في مراقبة التسربات. وسأل الباحثون عن عدد الصمامات الجديدة التي يمكن التحكم فيها عن بُعد والتي يجب إضافتها لتحسين قدرة النظام على التعامل مع انفجار أنبوب. ووجدت طريقتهم العدد الأمثل وأماكن وضع هذه الصمامات في غضاء ساعات قليلة. وفي المقابل، فإن محاولة اختبار كل مجموعة ممكنة من الصم valves باستخدام الطريقة القديمة الشاملة كانت ستستغرق سنوات من وقت الحاسوب. فعلى سبيل المثال، كان اختبار وضع خمسة صمامات جديدة فقط باستخدام النهج القديم سيتطلب أكثر من 1.5 مليار ثانية من وقت الحوسبة، وهي مهمة مستحيلة عمليًا. بينما حلت الطريقة الجديدة نفس المشكلة في أقل من خمس ساعات.
وقد حسنت التصاميم الناتجة عن هذا الإطار الجديد بشكل كبير من قدرة الشبكة على الحفاظ على تدفق المياه. ففي أسوأ سيناريو حيث يفشل أنبوب واحد، لم تتمكن الشبكة الحالية من تلبية سوى حوالي نصف إجمالي طلب المستهلكين. ومن خلال إضافة عدد قليل فقط من الصمامات الجديدة في الأماكن الصحيحة، أظهر الباحثون أن الشبكة يمكن أن تلبي ما يصل إلى 77% من الطلب أثناء الأزمة. وحتى في أكثر حالات الاختبار تعقيدًا مع خمسة صمامات جديدة، وُجد أن التصميم أقل كفاءة بنسبة 8% كحد أقصى من الحل المثالي النظري، وهو هامش خطأ ضئيل للغاية بالنسبة لمثل هذا النظام المعقد. وتؤكد الدراسة أنه من خلال استخدام نهج التخطيط الاستراتيجي ثلاثي الخطوات هذا، يمكن للمدن بناء شبكات مياه أكثر متانة ضد الفشل، مما يضمن أنه عندما ينفجر أنبوب، تستمر المياه في التدفق إلى الناس الذين يحتاجون إليها.
ملخص تقني: إطار تحسين ثلاثي المستويات للتصميم المرن لشبكات توزيع المياه تحت ظروف فشل الأنابيب
صياغة المشكلة تتناول هذه الورقة التحدي المتمثل في تصميم شبكات توزيع المياه (WDNs) المرنة القادرة على الحفاظ على استمرارية الإمداد للعملاء في ظل أسوأ سيناريوهات فشل الأناب الفردية. صاغ المؤلفون ذلك في شكل مسألة تحسين ثلاثية المستويات من نوع (المدافع-المهاجم-المدافع) (DAD).
المستوى العلوي (المدافع): يختار مجموعة من الروابط الجديدة (أنابيب أو صمامات) لتركيبها، مع مراعاة قيود الميزانية أو التصميم، لتعظيم مرونة الشبكة.
المستوى الأوسط (المهاجم): يختار سيناريو فشل أنبوب واحد هو الأسوأ (عزل أنبوب موجود) لتقليل أداء الشبكة.
المستوى الأدنى (المدافع): ينفذ استراتيجية استجابة للحوادث مثالية ومتوازنة هيدروليكيًا (ضبط حالات الصمامات وإدارة التدفقات) لتعظيم إجمالي الطلب المورد بعد الفشل.
تعد مسألة المستوى الأدنى مسألة برمجة خطية مختلطة بأعداد صحيحة غير خطية (MINLP) بسبب تضمين متغيرات التحكم الثنائية في الصمامات والقيود الهيدروليكية غير المحدبة (فقدان الضغط بالاحتكاك والطلبات المعتمدة على الضغط). ويشير المؤلفون إلى أن حل هذه المسألة (MINLP) الثلاثية الدقيقة للشبكات التشغيلية واسعة النطاق أمر مستعصٍ حاسوبيًا باستخدام الطرق الدقيقة الموجودة.
المنهجية للتغلب على الاستعصاء الحسابي، يقترح المؤلفون إطارًا يجمع بين التخفيفات الخطية وخوارزمية توليد الأعمدة والقيود (CCG) المتداخلة المصممة خصيصًا لهذا الغرض.
التخفيف الخطي: تم تخفيف القيود الهيدروليكية غير المحدبة في مسألة المستوى الأدنى باستخدام تقريبات خارجية متعددة الأوجه. وتحديدًا:
تم استبدال دوال الطلب المعتمدة على الضغط بمتعالية خطية مجزأة.
تم استبدال دوال فقدان الاحتكاك بمتعالية ومرجعية خطية مجزأة (أعلى وأسفل). هذا يحول مسألة المستوى الأدنى إلى مسألة برمجة خطية مختلطة بأعداد صحيحة (MILP)، مما يسمح بحل المسألة الثلاثية كنسخة مخففة (L-DAD).
خوارزمية CCG المتداخلة: يتم حل مسألة الـ (DAD) المخففة باستخدام نهج CCG متداخل:
الحد العلوي: تقوم حلقة CCG خارجية بإضافة سيناريوهات فشل (أعمدة) وقيود بشكل تكراري إلى مسألة رئيسية، مما يوفر حدًا علويًا للحل الأمثل للمسألة المخففة. ولضمان القدرة على المعالجة، تستخدم الخوارزمية تخفيفًا أقل صرامة للمسألة الرئيسية لتجنب الاعتماد المعقد على متغيرات التصميم أثناء تشديد الحدود.
الحد الأدنى: لكل تصميم مرشح يتم توليده في الحلقة العليا، تقوم خطوة الحد الأدنى بحل المسألة ثنائية المستوى (المهاجم مقابل المستجيب) لإيجاد أسوأ فشل واستجابة مثالية. تستخدم هذه الخطوة حلقة CCG داخلية لأن مسألة المستوى الأدنى تحتوي على متغيرات ثنائية.
التقارب: تتكرر الخوارمازمية حتى يصبح الفرق بين الحدين العلوي والسفلي ضمن نطاق سماحية محدد أو حتى استنفاد مجموعة التصاميم الممكنة.
التحقق ما بعد التحسين (تحليل الرابط الحرج): بما أن الحل يعتمد على قيود مخففة، يقوم المؤلفون بإجراء تحليل الرابط الحرج (CLA) "بعدي" (a posteriori) على التصميم الأمثل المحسوب. يتضمن ذلك حل مسألة (MINLP) غير الخطية الأصلية لسيناريو الفشل الأسوأ المحدد للحصول على حل ممكن وحد أدنى حقيقي لمسألة (DAD) الأصلية.
المساهمات الرئيسية
صياغة DAD لشبكات توزيع المياه: قامت الورقة بتكييف إطار عمل DAD، الذي كان شائعًا في أنظمة الطاقة، لنمذجة الهيدروليكا، والتحكم في الصمامات، وأحداث الفشل بشكل صريح في تخطيط شبكات المياه.
خوارزمية حل مخصصة: طور المؤلفون خوارزمية CCG متداخلة قادرة على التعامل مع المتغيرات الثنائية في كل من المستويات العليا/الوسطى والمستوى الأدنى، مقترنة بالتخفيفات الخطية للقيود الهيدروليكية غير المحدبة. يوفر هذا حدودًا عالمية للمسألة المخففة.
القابلية للتوسع: تم إثبات قدرة الطريقة على حل مشكلات التوسع الأمثل للشبكات التشغيلية واسعة النطاق بشكل أسرع بكثير من طرق البحث الشامل.
النتائج طُبق الإطار على مختبر ميدان مياه أنجليان (AWFL)، وهو شبكة تشغيلية واسعة النطاق في المملكة المتحدة تضم 1,167 نقطة التقاء و1,315 أنبوبًا. قامت الدراسة بتحسين وضع ما يصل إلى 5 صمامات جديدة للتحكم في الضغط ومنع الارتداد عن بُعد (PRVs).
الكفاءة الحسابية: قلل نهج CCG المتداخل من وقت الحوسبة بعدة مراتب مقارنة بالبحث الشامل القائم على CLA. على سبيل المثال، كان البحث الشامل لـ 5 صمامات جديدة سيستغرق أكثر من 1.5 مليار ثانية، بينما وصل نهج CCG المتداخل إلى الحد الزمني البالغ 5 ساعات (18,000 ثانية)، بينما استغرق التحقق اللاحق عبر تحليل الرابط الحرج (CLA) حوالي 23,567 ثانية.
فجوات المثالية: وُجد أن التصاميم المحسوبة كانت أقل مثالية بما لا يتجاوز 8% من الأداء الأسوأ الحقيقي (الفجوة بين الحد العلوي المخفف والحد الأدنى المتحقق منه).
مكاسب المرونة: حسنت التصاميم المحسنة أداء الشبكة في أسوأ الحالات بشكل كبير. أدى التصميم الأساسي المقسم (صفر صمامات جديدة) إلى سيناريو أسوأ أدى لقطع خمس مناطق إدارة توزيع (DMAs) متتالية تمثل 50% من إجمالي طلب الشبكة. أدى تركيب 3 صمامات جديدة إلى زيادة إجمالي الطلب المورد في ظل ظروف ذروة التحميل بنسبة 9% مقارنة بالتصميم الذي يحتوي على صمام واحد، كما استعاد الإمداد للمناطق التي كانت منفصلة في السيناريو الأساسي.
دقة التخفيف: بينما مالت التخفيفات الخطية إلى المبالغة في تقدير الرؤوس الهيدروليكية والطلب المورد، إلا أن الخطأ ظل ضمن نطاق يمكن إدارته. وأشارت النتائج إلى أن التخفيفات كانت أكثر إحكامًا لقيم التدفق الكبيرة (في الأنابيب الرئيسية)، والتي تهيمن على فقدان الضغط في الشبكة، مما أبقى عدم المثالية الإجمالية منخفضًا.
الأهمية والادعاءات تدعي الورقة أن هذا الإطار يقدم بديلًا فعالًا حاسوبيًا لمقاييس المرونة القائمة على النماذج البديلة أو عمليات البحث الشامل. ومن خلال نمذجة التفاعل بين قرارات التصميم، وأحداث الفشل، والاستجابات الهيدروليكية بشكل صريح، توفر الطريقة تصاميم تحسن بشكل كبير أداء الشبكة في أسوأ الحالات. ويؤكد المؤلفون أنه بينما يعتمد النهج على التخفيفات الخطية، فإن التصاميم الناتجة قوية، مع فجوات مثالية مثبتة تقل عن 8% للسيناريوهات المختبرة. ويخلصون إلى أن الطريقة مناسبة لقرارات التخطيط والتوسع في الشبكات التشغيلية واسعة النط scales، رغم إقرارهم بأن العبء الحسابي وفجوات المثالية يزدادان مع عدد التركيبات المحتملة، مما يشير إلى عمل مستقبلي على تقليم الاستدلال (heuristic pruning) والصياغات غير الخطية المباشرة.