KOOPMAN-Luenberger Observer Design for Nonlinear Systems with Application to the Monitoring of a Latent Thermal Energy Storage
تقترح هذه الورقة إطار عمل لمراقب حالة خطي قائم على "كوبمان" (KOOPMAN-LSO) يقوم برفع الديناميكيات غير الخطية إلى فضاء ذي أبعاد أعلى لتمكين تصميم مراقب خطي فعال، مما يثبت كفاءته في تقدير الحالات الداخلية غير المقاسة لنظام تخزين طاقة حرارية كامن بدقة عالية في ظل ظروف استشعار محدودة.
تخيل أنك تحاول تخمين حالة الطقس داخل جرة سحرية ضخمة ومغلقة بإحكام. لا يمكنك رؤية ما بداخلها، ولا يمكنك إدخال ميزان حرارة فيها دون كسر الجرة. كل ما تملكه هو نافذة صغيرة في الخارج لتشعر بحرارة الهواء المتدفق منها. إذا كانت الجرة بسيطة، مثل دلو من الماء، فمن المحتمل أن تتمكن من تخمين ما يحدث بالداخل بمجرد الشعور بهذا النسيم. ولكن ماذا لو كانت الجرة تحتوي على مادة غامضة متغيرة الشكل، تتجمد وتذوب، وتغير شخصيتها وسلوكها اعتماداً على مدى حرارتها أو برودتها؟ هذا هو اللغز الذي يواجه العلماء مع أنظمة "تخزين الطاقة الحرارية الكامنة" — وهي بطاريات ضخمة تخزن الحرارة عن طريق صهر مواد مثل الملح أو الشمع. هذه الأنظمة فعالة للغاية ولكنها صعبة المراقبة بشكل ملحوظ لأن داخلها مخفي، كما أن سلوكها غير خطي بشكل جامح، مما يعني أن تغيراً طفيفاً في المدخلات يمكن أن يسبب تحولاً هائلاً وغير متوقع في درجة الحرارة.
لحل هذه المعضلة، يحاول المهندسون عادةً بناء نماذج رياضية معقدة تحاكي فيزياء عملية الصهر، أو يعتمدون على شبكات مكلفة وفوضوية من المستشعرات. ولكن هناك طريقة أذكى وأكثر حداثة: استخدام البيانات لإيجاد الأنماط الخفية. فكر في الأمر كأنك تأخذ صورة ضبابية وفوضوية لعاصفة ثم تستخدم مرشحاً خاصاً لتحويلها إلى رسم بياني واضح وخط مستقيم يتنبأ بهبة الرياح التالية. تقدم هذه الورقة البحثية "مرشحاً" ذكياً جديداً يسمى "مراقب الحالة الخطي القائم على كوبمان" (KOOPMAN-based Linear State Observer). إنها طريقة تأخذ نظاماً فوضوياً وغير خطي وتعمل على "رفعه" إلى أبعاد أعلى حيث يبدو التخبط كخطوط مستقيرة بسيطة. وبمجرد تسوية هذا التخبط، يمكن للفريق استخدام رياضيات بسيطة وموثوقة (من النوع المستخدم في الهندسة الأساسية) لتخمين درجات الحرارة المخفية داخل الجرة بمجرد النظر إلى الهواء المتدفق خارجها. الأمر يشبه تحويل كرة متشابكة من الخيوط إلى خط مستقيم حتى تتمكن من سحب طرفها بسهءولة.
طبق الباحثون، بقيادة مصطفى حبيب وفريقه، هذه الفكرة على وحدة تخزين حراري حقيقية مملوءة بمادة خاصة تعتمد على الملح تسمى SP31. تعمل هذه المادة مثل إسفنجة حرارية، تمتص الحرارة أثناء ذوبانها وتطلقها أثناء تجمدها. كان التحدي يكمل في امتلاكهم مستشعراً واحداً فقط يقيس درجة حرارة الماء الخارج من النظام، ومع ذلك كانوا بحاجة لمعرفة درجة حرارة المادة في عمق النظام. قاموا ببناء "توأم رقمي" للنظام باستخدام تقنية eDMDc، التي تتعلم كيفية سلوك النظام من خلال مراقبة دورات الشحن والتفريغ. وبدلاً من محاولة نمذجة الفيزياء المعقدة للصهر مباشرة، استخدموا حيلًا رياضية "مستوحاة من الفيزياء" — وهي دوال تشبه الطريقة التي يحدث بها الصهر — لرفع البيانات إلى ذلك الفضاء عالي الأبعاد حيث تصبح القواعد بسيطة وخطية.
بمجرد "رفع" النظام، قاموا بإرفاق "مراقب لونبرجر" (Luenberger Observer)، وهو في الأساس "مخمن ذكي" يصحح توقعاته باستمرار. تخيل طالباً يؤدي اختباراً وهو يعلم أن نموذج الإجابة يختلف قليلاً عن تخمينه؛ في كل مرة يحصل فيها على دليل جديد (درجة حرارة المخرج)، يقوم بتعديل إجابته ليقترب من الحقيقة. قام الفريق بضبط هذا المخمن باستخدام طريقة تسمى LQR، والتي تعمل مثل "مقبض التحكم في الصوت"، حيث توازن بين سرعة تعلم المخمن وبين مدى ثقته في بيانات المستشعر المشوشة. اختبروا هذا الإعداد في مختبر باستخدام خزان سعة 200 لتر، مع تشغيله عبر سيناريوهات مختلفة بمعدلات تدفق مياه مختلفة (2، 4، و8 لتر في الدقيقة) ودرجات حرارة متفاوتة.
كانت النتائج مبهرة. عندما كان النظام يعمل بوتيرة ثابتة، استطاع المراقب التنبؤ بدرجة حرارة الماء الخارج من الخزان بدقة مذهلة، حيث لم يتجاوز الخطأ 0.0819 درجة مئوية. والأهم من ذلك، نجح في إعادة بناء درجات الحرارة المخفية للمادة التي تذوب بالداخل، محافظاً على متوسط خطأ أقل من 1.0 درجة مئوية لمعظم المستشعرات. ومع ذلك، كانت الورقة صريحة بشأن حدودها؛ فعندما تم رفع تدفق المياه إلى 8 لتر في الدقيقة، أو عندما تحول النظام فجأة من التسخين إلى التبريد، استغرق "التخمين" وقتاً أطول للحاق بالواقع، وزادت الأخطاء قليلاً لتصل إلى 2.3072 درجة مئوية في بعض النقاط المحددة داخل الخزان. ويوضح المؤلفون أن هذا يحدث لأن بعض أجزاء الخزان بعيدة جداً عن المستشعر الوحيد بحيث لا يمكن رؤيتها بوضوح، تماماً مثل محاولة سماع همس من خلف غرفة صاخبة.
من الأهمية بمكان أن الورقة تجادل ضد فكرة أنك بحاجة إلى نموذج مثالي كثيف الفيزياء أو مستشعر في كل زاوية لمراقبة هذه الأنظمة. كما تقترح أنه بينما توجد طرق أخرى مثل "مرشح كالمان الموسع" (Extended Kalman Filter)، إلا أنها قد تكون ثقيلة حسابياً ويصعب ضبطها. يقدم هذا النهج الجديد "حلاً وسطاً": فهو ليس كرة بلورية مثالية ترى كل شيء فوراً، ولكنه أداة قوية وفعالة من الناحية الحسابية تحول مشكلة غير خطية وفوضوية إلى مشكلة خطية يمكن إدارتها. يشير المؤلفون إلى أنه بينما يعمل هذا الأسلوب جيداً في الظروف المستقرة، فإن العمل المستقبلي قد يجمع بينه وبين تقنيات أخرى تعتمد على البيانات للتعامل مع نطاقات أوسع من ظروف التشغيل. وفي الوقت الحالي، يظل هذا النهج طريقة واعدة لمراقبة القلب الخفي للبطاريات الحرارية دون الحاجة إلى حفر ثقوب فيها.
ملخص تقني: تصميم مراقب كوپمان-لونبرجر للأنظمة غير الخطية مع التطبيق على تخزين الطاقة الحرارية الكامنة
بيان المشكلة لا يزال التقدير الدقيق لحالة الأنظمة الديناميكية غير الخطية، وخاصة أنظمة تخزين الطاقة الحرارية الكامنة (LTES) القائمة على المواد متغيرة الطور (PCMs)، يمثل تحدياً كبيراً. تظهر هذه الأنظمة عدم خطية قوية بسبب تحولات الطور، والتباطؤ (hysteresis)، وتدرجات درجات الحرارة الموزعة مكانياً. وفي التطبيقات العملية، تكون الحالات الداخلية (مثل درجات حرارة مادة الـ PCM المحددة) غير قابلة للوصول إلى حد كبير، حيث يقتصر الرصد غالباً على قياسات حدودية متفرقة مثل درجات حرارة دخول وخروج السائل. تواجه الحلول الحالية مقايضات: فالنماذج عالية الدقة القائمة على الفيزياء (مثل ديناميكا السوائل الحسابية CFD) مكلفة حاسوبياً للاستخدام في الوقت الفعلي؛ والنماذج منخفضة الرتبة تعاني غالباً من الحساسية تجاه المعلمات؛ أما النهج القائم كلياً على البيانات (مثل الشبكات العصبية) فيفتقر إلى التفسير والقدرة على التكيف خارج نطاق تدريبها. علاوة على ذلك، فإن المراقبين غير الخطيين التقليديين (مثل مرشح كالمان الموسع EKF) يمكن أن يكون من الصعب ضبطهم ويفتقرون إلى ضمانات الاستقرار العالمي.
المنهجية يقترح المؤلفون إطار عمل مراقب الحالة الخطي القائم على كوپمان (KOOPMAN-LSO). يستفيد هذا النهج من نظرية مؤثر كوپمان لرفع (lifting) ديناميكيات النظام غير الخطي إلى فضاء أعلى من الأبعاد للمتغيرات المرصودة حيث يمكن تقريب الديناميكيات كأنها خطية. وتتقدم المنهجية عبر ثلاث مراحل رئيسية:
الرفع المستند إلى الفيزياء وتحديد eDMDc:
يتم رفع النظام غير الخطي إلى فضاء متغيرات مرصودة أعلى الأبعاد باستخدام دوال أساس مستندة إلى الفيزياء. وتحديداً، بالنسبة لديناميكيات الـ PCM، يدمج المؤلفون دالة سيجمويد (Sigmoid) لنمذجة تحول الطور (صلب-سائل) ودالة غاوس (Gaussian) لنمذجة تغير السعة الحرارية النوعية حول نقطة الانصهار.
بخلاف عملية الرفع القائمة كلياً على البيانات (مثل المشفرات التلقائية)، فإن هذا الاختيار يدمج معرفة فيزيائية مسبقة لتحسين التعميم والقدرة على التفسير.
يتم تحديد متنبئ خطي مع مدخلات التحكم باستخدام طريقة تحليل النمط الديناميكي الموسع مع التحكم (eDMDc). ينتج عن ذلك نموذج حالة في الزمن المنفصل (zk+1=Azk+Buk) يقرب الديناميكيات غير الخطية الأصلية.
تخليق المراقب:
يتم بناء مراقب لونبرجر (Luenberger observer) في الزمن المنفصل ضمن فضاء الإحداثيات المرفوع.
يتم تحديد كسب المراقب (L) باستخدام صيغة التحكم التربيعي الخطي (LQR) المزدوج. هذا النهج الحتمي يقلل من دالة الهدف التي توازن بين تقارب خطأ التقدير والتنبؤ بالقياس، متجنباً افتراضات الضوضاء العشوائية المطلوبة في ترشيح كالمان.
يتم رسم خرائط الحالات المرفوعة المقدرة عائدة إلى الحالات الفيزيائية (درجات حرارة الـ PCM) باستخدام مصفوفة اختيار.
التحقق التجريبي:
تم اختبار إطار العمل على نظام LTES مخبري يستخدم مادة PCM غير عضوية من نوع SP31.
تعرض النظام لدورات شحن وتفريغ مع معدلات تدفق مختلفة لسائل نقل الحرارة (HTF) وهي (2، 4، و8 لتر/دقيقة).
صُمم المراقب لإعادة بناء تسعة حالات داخلية لدرجة حرارة الـ PCM باستخدام قياس واحد فقط: درجة حرارة مخرج سائل الـ HTF.
المساهمات الرئيسية
تطوير إطار العمل: إنشاء إطار KOOPMAN-LSO الذي يتيح تخليق مراقب خطي للأنظمة الحرارية غير الخطية مع مدخلات التحكم.
الرفع المستند إلى الفيزياء: تقديم دوال أساس محددة للمجال (السيجمويد والغاوس) لديناميكيات الـ PCM، مما يعزز متانة النموذج وقابليته للتفسير مقارنة بالطرق المعتمدة على البيانات "الصندوق الأسود".
تصميم الكسب الحتمي: صياغة مراقب لونبرجر مع ضبط الكسب بناءً على LQR في الفضاء المرفوع، مما يوفر تقديراً مستقراً وقابلاً للضبط دون الاعتماد على افتراضات الضوضاء العشوائية.
التحقق التجريبي: التحقق الشامل الذي أظهر القدرة على إعادة بناء الحالات الحرارية الداخلية من قياسات حدودية متفرقة تحت ظروف تشغيل واقعية ومتغيرة.
النتائج تقريرت نتائج الدراسة وفق مقايطات الأداء التالية:
درجة حرارة المخرج: حقق المراقب إعادة بناء عالية الدقة لدرجة حرارة مخرج الـ HTF، مع جذر متوسط مربع الخطأ (RMSE) يصل إلى 0.0819 درجة مئوية.
درجات حرارة الـ PCM الداخلية: في ظروف التدفق القياسية (2 لتر/دقيقة)، كان RMSE لدرجات حرارة الـ PCM الداخلية عموماً أقل من 1.0 درجة مئوية.
الحساسية لمعدل التدفق: انخفضت دقة التقدير مع زيادة معدلات التدفق. عند 8 لتر/دقيقة، وصلت قيم RMSE لبعض المستشعرات إلى 2.3072 درجة مئوية، لا سيما أثناء التحولات السريعة بين الشحن والتفريغ.
حدود القابلية للملاحظة: كشف تحليل القابلية للملاحظة عن رتبة 8/10، مما يشير إلى أن اتجاهين من درجات الحرارة الداخلية غير قابلين للملاحظة هيكلياً من قياس المخرج الواحد. وبناءً على ذلك، يلتقط المراقب المحتوى الطاقي العام والاتجاهات المتوسطة، لكنه لا يستطيع تمييز الميزات المكانية الدقيقة مثل انتشار جبهة الطور بشكل فريد.
السلوك العابر: أظهر المراقب انحرافات مؤقتة أثناء التغيرات الحادة في درجة حرارة المدخل (التبديل بين التسخين والتبريد)، مما يتطلب عدداً محدداً من التكرارات لاستقرار ديناميكيات الخطأ.
الأهمية والادعاءات يدعي البحث أن KOOPMAN-LSO يمثل "منطقة وسطى" عملية بين النمذجة الفيزيائية المعقدة (الصندوق الأبيض) والنهج المعتمد كلياً على البيانات (الصندوق الأسود). تكمن أهميته الأساسية في:
البساطة التحليلية: من خلال تحويل الديناميكيات غير الخطية إلى تمثيل خطي، يسمح إطار العمل باستخدام أدوات التحكم والتقدير الخطية الناضجة والبسيطة تحليلياً (مراقبي لونبرجر، LQR) والتي توفر ضمانات استقرار أقوى من البدائل غير الخطية مثل EKF.
الكفاءة الحسابية: هذا النهج خفيف حسابياً، مما يجعله مناسباً لتطبيقات المراقبة والتحكم في الوقت الفعلي حيث يكون استخدام CFD عالي الدقة أو التحسين غير الخطي المعقد (مثل تقدير الأفق المتحرك MHE) غير ممكن.
المتانة تجاه الاستشعار المتفرق: يثبت أنه يمكن تحقيق تقدير دقيق للحالات الحرارية العالمية (حالة الشحن) حتى مع نشر ضئيل للمستشعرات، بشرভাবে استخدام متغيرات مرصودة ذات معنى فيزيائي للرفع.
يقر المؤلفون بتواضع بأن الطريقة ليست تهدف إلى تحقيق دقة مطلقة لكل حالة غير مقاسة، خاصة تلك غير القابلة للملاحظة هيكلياً. كما يشيرون إلى أن الأداء محكوم بقوة التبعيات الوظيفية بين القياسات المتفرقة ومدخلات النظام. وقد تم تحديد العمل المستقبلي في استقصاء دمج "تحديد الديناميكيات غير الخطية المتفرقة" (SINDy) لتعزيز التفسير وزيادة نطاق صلاحية النماذج عبر نوافذ تشغيل أوسع.