Active learning-driven experimental design for torrefaction process optimization using multi-objective Gaussian process surrogates
تقدم هذه الورقة إطار عمل للتعلم النشط يستخدم بدائل العمليات الغاوسية متعددة الأهداف لتحسين عملية تحلل الكتلة الحيوية حرارياً بكفاءة، وذلك من خلال التوصية بتجارب عالية القيمة توازن بين مردود الكتلة، والقيمة الحرارية العليا، ومردود الطاقة، مع تقليل التكاليف التجريبية بشكل كبير.
المؤلفون الأصليون:Mohammed R. AlAttar¹, Amr Y. Elbanhawy¹, Zakaria Ghoneim²
تُعد الكتلة الحيوية، وهي المواد العضوية المستمدة من النباتات والأشجار، واعدة للغاية كمصدر للطاقة المتجددة، ولكن في حالتها الخام، غالبًا ما يكون من الصعب التعامل معها ونقلها وحرقها بكفاءة. فهي مليئة بالرطوبة، وتتفتت بسهء، وتفتقر إلى الكثافة الطاقية التي يتمتع بها الفحم. ولحل هذه المشكلة، يستخدم العلماء عملية تسمى "التوريفاكسيون" (torrefaction)، وهي تسخين الكتلة الحيوية بلطف في بيئة مغلقة وخالية من الأكسجين. هذا الطهي الخفيف، الذي يحدث عادةً بين 200 و320 درجة مئوية، يحول المادة الخام إلى وقود صلب وجاف يشبه الفحم، مما يسهل طحنه ويجعله يحترق بشكل أنظف. ومع ذلك، فإن العثور على ظروف الطهي المثالية يمثل لغزًا معقدًا؛ إذ تعتمد جودة الوقود النهائي بشكل كبير على عاملين رئيسيين: مدى سخونة الفرن ومدة بقاء الكتلة الحيوية بداخله. فإذا كانت درجة الحرارة منخفضة جدًا أو الوقت قصيرًا جدًا، سيظل الوقود ضعيفًا؛ وإذا كانت مرتفعة جدًا أو الوقت طويلًا جدًا، فسيُفقد الكثير من المادة القيمة على شكل دخان وغاز. وبما أن الكتلة الحيو تختلف من دفعة إلى أخرى، وأن عدد مجموعات درجات الحرارة والوقت الممكنة هائل، فإن العثور على الإعدادات المثالية من خلال تجارب المحاولة والخطأ التقليدية هو أمر بطيء ومكلف وغير فعال في كثير من الأحيان.
لقد طور باحثون في جامعة عين شمس في مصر طريقة جديدة لحل هذا اللغز، تعمل بمثابة دليل ذكي للمختبر. فبدلاً من اختبار كل مجموعة ممكنة من الحرارة والوقت بشكل عشوائي، قاموا ببناء نموذج حاسوبي يتعلم من البيانات الموجودة للتنبؤ بنتيجة التجارب الجديدة قبل إجرائها حتى. هذا النموذج، المعروف باسم "العملية الغاوسية" (Gaussian process)، مميز لأنه لا يكتفي بتخمين النتيجة فحًا، بل يحسب أيضًا مدى عدم تأكده من ذلك التخمين. فكر في الأمر كأنه توقع جوي لا يخبرك فقط بأنها ستمطر، بل يخبرك أيضًا بمدى ثقة خبير الأرصاد الجوية في هذا التوقع. ومن خلال استخدام مستوى الثقة هذا، أنشأ الباحثون نظام "تعلم نشط" (active learning). يعمل هذا النظام كالمستكشف الفضولي الذي يعرف بالضبط أين يبحث تاليًا؛ فهو يوازن بين هدفين: اختبار الظروف التي يكون فيها النموذج جيدًا بالفعل في التنبؤ بالنتائج لصقل الإجابة، واختبار الظروف التي يكون فيها النموذج غير متأكد، ليتعلم شيئًا جديدًا. كما يضمن النظام أن تكون التجارب التي يقترحها موزعة عبر النطاق الكامل من الاحتمالات، مما يمنع الباحثين من إضاعة الوقت في اختبار نفس الظروف مرتين.
قام الفريق بتغذية نموذجهم الحاسوبي ببيانات من 486 تجربة توريفاكسيون سابقة، تضمنت تفاصيل حول درجة حرارة تسخين الكتلة الحيوية، والوقت الذي سُخنت فيه، والتركيب الكيميائي للمادة النباتية. كُلِّف النموذج بمهمة التنبؤ بثلاث نتائج رئيسية: مقدار الوقود الصلب المتبقي بعد التسخين، وكمية الطاقة التي يحتوي عليها هذا الوقود، وإجمالي عائد الطاقة. وعند اختباره مقابل بيانات لم يسبق له رؤيتها، أثبت النموذج دقة ملحوظة، حيث تنبأ بشكل صحيح بمحتوى الطاقة للوقود في ما يقرب من 97 في المائة من الحالات، وبكمية الوقود المتبقية في حوالي 85 في المائة من الحالات. وقد تفوق أداؤه بشكل كبير على تقنيات النمذجة الحاسوبية الشائعة الأخرى، مثل "الغابات العشوائية" (random forests) أو "الشبكات العصبية" (neural networks)، والتي تُستخدم غالبًا لمهام مماثلة. والأهم من ذلك، تحقق الباحثون من أن تقديرات عدم اليقين في النموذج موثوقة؛ فقد وجدوا أنه عندما يقول النموذج إنه واثق بنسبة 95 في المائة من تنبؤ ما، فإن النتيجة الفعلية تقع ضمن هذا النطاق في 93 في المائة من الحالات بالنسبة لعائد الوقود، وفي 95 في المائة من الحالات بالنسبة لإجمالي الطاقة. وهذا الموثوقية ضرورية، لأنها تعني أن اقتراحات النظام حول مكان إجراء التجربة التالية تستند إلى أرضية صلبة.
باستخدام هذا النموذج الموثوق، طلب الباحثون من الكمبيوتر التوصية بالعشر تجارب القادمة التي ستعلمهم أكثر قدر ممكن. لم يكتفِ النظام باختيار أفضل ظرف واحد، بل اختار مجموعة متنوعة من عشرة إعدادات مختلفة تغطي الطيف الكامل لعملية التوريفاكسيون. وقد شملت التجارب الموصى بها درجات حرارة تتراوح من 200 إلى 315 درجة مئوية وأوقات تسخين تتراوح من 30 ثانية إلى ما يقرب من 1,100 ثانية. وأظهرت النتائج وجود مقايضة واضحة: فالظروف الأخف، حوالي 200 إلى 215 درجة، أنتجت كمية أكبر من الوقود الصلب ولكن بكثافة طاقية أقل قليلاً، بينما أنتجت الظروف الأكثر قسوة، بالقرب من 315 درجة، وقوداً أقل ولكن بمحتوى طاقي أعلى بكثير. وقد وقعت اقتراحات الكمبيوتر عند حافة هذه المقايضة تمامًا، مما قدم مجموعة متوازنة من التجارب التي يمكن أن تساعد العلماء في تحديد ما إذا كانوا يريدون تعظيم كمية الوقود المنتجة أو الطاقة المعبأة في كل كيلوغرام. كما كشف التحليل أن درجة الحرارة والوقت كانا العاملين الأكثر أهمية بلا منازع، بينما كانت نوع المادة النباتية أقل أهمية مما كان متوقعًا، على الأقل ضمن نطاق البيانات التي درسوها.
يقدم هذا النهج خارطة طريق عملية للعلماء والمهندسين العاملين في تحسين أنواع الوقود المتجددة. فمن خلال استخدام نظام يتعلم من البيانات السابقة ويختار التجارب المستقبلية بذكاء، يمكن للمختبرات تجاوز التخمين والتركيز على الاختبارات الأكثر إفادة. لقد أثبت الباحثون أن طريقتهم قوية، بمعنى أن التغييرات الصغيرة في كيفية وزن الكمبيوتر لخياراته لا تغير التوصيات النهائية بشكل جذري. وبينما اعتمدت الدراسة الحالية على نوع ثابت من المواد النباتية وجولة واحدة من التوصيات، فإن الإطار مصمم ليتم تحديثه. وفي المستقبل، يمكن استخدام هذا النظام نفسه في حلقة مستمرة، حيث يتم إدخال نتائج التجارب الجديدة مرة أخرى في النموذج لصقل تنبؤاته وتوجيه مجموعة التجارب التالية. أما في الوقت الحالي، فيوفر هذا العمل مسارًا واضحًا مدفوعًا بالبيانات، موضحًا أنه مع الأدوات المناسبة، يمكن تحسين العملية المعقدة لتحويل النباتات إلى وقود عالي الجودة مع عدد أقل بكثير من التجارب وفهم أوضح لما يعمل بشكل أفضل.
ملخص تقني لورقة بحثية: التصميم التجريبي المدفوع بالتعلم النشط لتحسين عملية التوريفيكشن باستخدام بدائل العمليات الغاوسية متعددة الأهداف
بيان المشكلة
يواجه تحسين عملية "التوريفيكشن" (Torrefaction) للكتلة الحيوية — وهي عملية تفحم جزئي تُستخدم لرفع جودة وقود الكتلة الحيوية الخام وتحويله إلى وقود صلب شبيه بالفحم — عائقاً يتمثل في التكلفة العالية والوقت الطمر المطلوب للاستكشاف التجريبي الواسع. وتتسم هذه العملية بحساسية عالية لظروف التشغيل (درجة الحرارة ووقت المكوث) وتكوين المادة الخام، مما يخلق فضاءً بارامترياً عالي الأبعاد، تجعل من أساليب "تغيير عامل واحد في كل مرة" التقليدية غير فعالة، ومن التصاميم التجريبية كاملة العوامل (Full Factorial Designs) باهظة التكلفة. وبينما ظهرت النمذجة البديلة القائمة على البيانات كبديل، ركزت الدراسات السابقة في مجال التوريفيكشن بشكل كبير على دقة التنبؤ وحدها، وفشلت في استغلال عدم اليقين التنبؤي لتصميم تجارب جديدة وأكثر إفادة. علاوة على ذلك، فإن التطبيقات الحالية لـ "التحسين البايزي" (Bayesian Optimization) في المجالات ذات الصلة تفتقر غالباً إلى اختيار الدفعات (Batches) المتنوعة، والمعايرة الصارمة لعدم اليقين، والتحليل المنهجي لأهمية الميزات المستمدة من النموذج نفسه.
المنهجية
طوّر المؤلفون إطار عمل شاملاً للتعلم النشط مصمماً للتوصية بأكثر التجارب إفادة للخطوات القادمة بهدف تحسين إنتاج الكتلة (Mass Yield)، والقيمة الحرارية العليا (HHV)، وإنتاج الطاقة (Energy Yield). تتكون المنهجية من العناصر التالية:
مجموعة البيانات والمعالجة المسبقة: استخدمت الدراسة مجموعة بيانات منقحة مكونة من 486 ملاحظة مستمدة من تجارب داخلية ومن الأدبيات العلمية. شملت الميزات المدخلة درجة الحرارة، ووقت المكوث، والتحليل التقريبي (الرماد، المواد المتطايرة، الكربون الثابت). تم تعويض القيم المفقودة باستخدام المتوسطات لكل عمود، وتم توحيد البيانات (Standardization).
النمذجة البديلة (Surrogate Modeling): تم تدريب نماذج انحدار العمليات الغاوسية (GPR) مع نوى "تحديد الأهمية التلقائية" (ARD) من النوع الأسي المربع لكل متغير مستهدف. تم تحسين المعلمات الفائقة (Hyperparameters) عبر تعظيم الاحتمالية الهامشية. ولضمان قابلية إعادة الإنتاج، تم إعادة تدريب النماذج بشكل حتمي باستخدام المعلمات الفائقة المستخرجة.
المقارنة المرجعية الأساسية: تمت مقارنة نماذج GPR مع نماذج الغابة العشوائية (RF)، وانحدار متجه الدعم (SVR)، والشبكة العصبية الاصطناعية الضحلة (ANN) باستخدام التحقق المتقاطع خماسي الطيات (5-fold cross-validation).
فضاء المرشحين: تم إنشاء فضاء للمرشحين فوق نافذة التوريفيكشن التقليدية (200-320 درجة مئوية) مع هامش توسع بنسبة 10%. تضمن ذلك شبكة (10×10) لدرجة الحرارة ووقت المكوث، بالإضافة إلى 500 نقطة تم توليدها عبر أخذ عينات "لاتين فائق التوزيع" (LHS) لضمان تغطية واسعة.
دالة الاستحواذ متعددة الأهداف: تم بناء دالة استحواذ موزونة توازن بين ثلاثة مكونات معيارية:
عدم اليقين (40%): متوسط التباين التنبؤي عبر جميع الأهداف الثلاثة لتعزيز الاستكشاف.
التحسين المتوقع لإنتاج الكتلة (30%): لتعظيم الإنتاج.
التحسين المتوقع للقيمة الحرارية العليا (30%): لتعظيم كثافة الطاقة. ملاحظة: لم يتم تخصيص حد "تحسين متوقع" مستقل لإنتاج الطاقة، كونه كمية مشتقة؛ حيث يتم التقاط تأثيره بشكل غير مباشر من خلال مكون عدم اليقين.
اختيار الدفعة الواعي بالتنوع: اختارت خوارزمية جشعة متسلسلة دفعة مكونة من 10 تجارب. قامت هذه الخوارزمية بمعاقبة المرشحين القريبين من النقاط المختارة بالفعل باستخدام عامل مكافأة التنوع (λ=0.3) لضمان التغطية المكانية في فضاء التصميم.
التحقق والتحليل: تضمنت الدراسة التحقق المتقاطع لمعايرة عدم اليقين (للتحقق مما إذا كانت البواقي المعيارية تتبع التوزيع الطبيعي القياسي)، واستخراج أهمية الميزات بناءً على ARD، ودراسات الاستئصال (Ablation studies) لاختبار مساهمات المكونات، وتحليلات الحساسية تجاه تغيرات المعلمات الفائقة.
النتائج الرئيسية
أداء النموذج: تفوقت نماذج GPR على جميع النماذج المرجعية (RF, SVR, ANN) في جميع الأهداف. بلغت درجات R2 في التحقق المتقاطع 0.854 لإنتاج الكتلة، و 0.971 للقيمة الحرارية العليا، و 0.869 لإنتاج الطاقة. حقق نموذج GPR للقيمة الحرارية العليا أعلى دقة، بينما أظهر إنتاج الكتلة تبايناً ذاتياً أعلى.
معايرة عدم اليقين: سجلت الدراسة أول معايرة لعدم اليقين عبر التحقق المتقاطع لنماذج GPR في مجال التوريفيكشن؛ حيث كانت نسبة التغطية ضمن الفترة الموثوقة 95% هي 93.2% لإنتاج الكتلة و 95.1% لإنتاج الطاقة، مما يشير إلى معايرة مثالية تقريباً. أظهر نموذج القيمة الحرارية العليا معايرة متحفظة (تغطية 99.8%)، ويرجع ذلك على الأرجح إلى دقتها العالية التي تؤدي إلى فترات تنبؤية أوسع.
أهمية الميزات: حددت أطوال النطاق (Length scales) في نظام ARD أن درجة الحرارة هي الميزة المهيمنة لجميع الأهداف (مسؤولة عن ما يقرب من 80% من الأهمية للقيمة الحرارية العليا)، تليها زمن المكوث. أما متغيرات التحليل التقريبي (الرماد، المواد المتطايرة، الكربون الثابت) فقد أظهرت أهمية منخفضة باستمرار، مما يبرر التركيز على درجة الحرارة والزمن في فضاء المرشحين.
التوصيات التجريبية: أوصى إطار العمل بدفعة من 10 تجارب تمتد بين 200–315 درجة مئوية و 30–1089 ثانية. انقسمت الدفعة طبيعياً إلى مجموعتين:
التوريفيكشن الخفيف (200–215 درجة مئوية): إنتاج كتلة أعلى (65–75%) ولكن قيمة حرارية عليا أقل.
التوريفيكشن الشديد (302–315 درجة مئوية): إنتاج كتلة أقل (45–55%) ولكن قيمة حرارية عليا أعلى بكثير (26–28 ميجاجول/كجم).
يعكس هذا التوزيع المقايضة بين الحفاظ على الإنتاج وتكثيف الطاقة.
الاستئصال والحساسية: أكدت دراسات الاستئصال أن الجمع بين عدم اليقين، والتحسين المتوقع، وقيود التنوع كان ضرورياً لتحقيق انتشار مكاني واسع (4.1 في الفضاء المعياري) مقارنة باستخدام عدم اليقين أو التحسين المتوقع بمفردهما. كما أظهر تحليل الحساسية أن المجموعة الموصى بها قوية تجاه التغيرات المتوسطة في أوزان الاستحواذ وعوامل التنوع.
جبهة باريتو (Pareto Front): وُجد أن التجارب الموصى بها تقع بالقرب من جبهة باريتو المثلى لإنتاج الكتلة مقابل القيمة الحرارية العليا، مما يؤكد قدرة دالة الاستحواذ على تحديد النقاط التي تدفع حدود المقايضة.
الأهمية والادعاءات
تدعي الورقة أنها الدراسة الأولى في أدبيات التوريفيكشن التي تدمج في آن واحد تنوع الدفعات، ومعايرة عدم اليقين عبر التحقق المتقاطع، وأهمية الميزات المستمدة من ARD ضمن إطار تعلم نشط متعدد الأهداف.
يضع المؤلفون عملهم كمنهجية منظمة، قابلة للتفسير، وواعية بعدم اليقين لتسريع تطوير العمليات. وتشمل المساهمات الرئيسية ما يلي:
النمذجة البديلة الصارمة: إثبات أن نماذج GPR مع نوى ARD توفر دقة فائقة وتقديرات عدم يقي ضرورية مقارنة بالنماذج المرجعية الأخرى للتعلم الآلي.
عدم يقين موثوق: التحقق من أن تقديرات عدم اليقين في نماذج GPR موثوقة (عبر المعايرة) لتوجيه المقايضة بين الاستكشاف والاستغلال.
تصميم قابل للتنفيذ: تقديم دفعة محددة وواعية بالتنوع مكونة من 10 تجارب تغطي كلاً من أنظمة المعالجة الخفيفة والشديدة ضمن نافذة التوريفيكشن التقليدية، مما يجعلها قابلة للتطبيق الفوري في المختبر.
النمطية: يُقدم إطار العمل كمنهج فعال حاسوبياً وقابل للنقل إلى مشاكل تحسين العمليات الكيميائية متعددة الأهداف الأخرى.
يصرح المؤلفون صراحةً أن إطار العمل يهدف إلى تحديد القياسات المستقبلية الأكثر إفادة في ظل الموارد المحدودة بدلاً من استبدال التجريب. كما يقرون بوجود قيود، مثل افتراض تكوين ثابت للمادة الخام وغياب التغذية الراجعة التجريبية التكرارية في الدراسة الحالية، مشيرين إلى أنه سيتم معالجة هذه الأمور في أعمال مستقبلية.