في هذا البحث، تم تطوير نظام شحن للمركبات الكهربائية (EV) يعمل بنظام مستقل يعتمد على الخلايا الكهروضوئية ثنائية الوجه (bifacial PV) ونظام تخزين طاقة البطاريات (BESS)، كما هو موضح في الشكل 1، باستخدام طريقة تتبع نقطة الطاقة القصوى (MPPT) التي يتم التحكم فيها بواسطة الشبكة العصبية الاصطناعية (ANN). تم استخدام دائرة مكافئة أحادية الثنائيات لنمذجة المصفوفة الكهروضوئية ثنائية الوجه، مع مراعاة درجة الحرارة، والمقاومة المتسلسلة، ومقاومة التوازي، والإشعاع الشمسي، والإشعاع من الجانب الخلفي. ويتم تزويد محطة شحن المركبات الكهربائية بمخرجات النظام الكهروضوئي من خلال محول رفع تيار مستمر-تيار مستمر (DC–DC boost converter)، بينما يدعم نظام BESS تنظيم جهد ناقل التيار المستمر (DC-bus). كما تم استخدام خوارزمية تحسين طائر السكرتير (SBOA) لتحسين مكاسب متحكم PI ومعاملات المرشح. وفيما يلي وصف لنمذجة مكونات النظام والمنهجية المقترحة في الأقسام الفرعية التالية.
نمذجة المكونات
1. نظام الخلايا الكهروضوئية ثنائية الوجه
تساعد الوحدات الكهروضوئية (PV) في تحقيق ناتج التيار والجهد المطلوبين. في هذه الدراسة، تم استخدام وحدتين كهروضوئيتين أحاديتي الوجه لمحاكاة خصائص توليد الطاقة للمصفوفات ثنائية الوجه تحت ظروف مختلفة من وضاءة السطح (albedo)9. يتم تحديد الكسب ثنائي الوجه للإشعاع BGg بواسطة المعادلة (1)
(1)
هنا، يمثل Gf الإشعاع الأمامي، ويمثل Gr الإشعاع الخلفي. أما أقصى تيار للخلايا الكهروضوئية ثنائية الوجه فهو Imp، والجهد هو Vmp. وأخيرًا، يتم تحديد المخرجات الناتجة عن النظام ثنائي الوجه مع العدد المختار من الوحدات المتصلة على التوالي Ns وعلى التوازي Np بواسطة المعادلة (2). يظهر التحكم في الخلايا الكهروضوئية في الشكل 2
(2)
2. نظام تخزين البطارية (BSS)
تعد البطاريات وسيلة حاسمة لتخزين الطاقة الكهربائية18 . يستخدم هذا العمل بطاريات Li-ion كنظام لتخزين الطاقة. ومن بين خيارات البطاريات القابلة لإعادة الشحن المتاحة، يركز هذا العمل على بطاريات Li-ion، والتي تُعتبر مثالية للمركبات الكهربائية (EVs) نظرًا لكفاءتها وكثافة طاقتها. تصف المعادلتان (3) و(4) بطارية Li-ion مع القيود.
(3)
تم تلخيص مواصفات الخلايا الكهروضوئية (PV)، وبطارية التخزين، والمحول في الجدول 2، بينما يوضح الجدول 3 توزيع تدفق الطاقة بين مصفوفة الخلايا الكهروضوئية، ووحدة تخزين الطاقة، والأحمال المتصلة.
(4)
3. محول التيار المستمر إلى تيار متردد (DC –AC converter)
يتم اعتماد العاكس ذو نقطة التعادل الموحدة (NPC) على نطاق واسع في أنظمة تحويل الطاقة الحديثة للطاقة المتجددة. ومقارنة بالبنية التقليدية ذات المستويين، يوفر تكوين NPC تشويهاً توافقياً كلياً (THD) أقل، وتداخلاً كهرومغناطيسياً أقل، وإجهاد جهد أقل على أجهزة أشباه الموصلات. ينقسم رابط التيار المستمر (DC link) في هذه الطوبولوجيا إلى ثلاثة مستويات جهد هي +Vdc/2 و 0 و -Vdc/2، وذلك عن طريق مكثفات مزدوجة وصمامات ثنائية للتثبيت تشكل وصلة التعادل. ومن خلال التحكم المناسب في عملية التبديل، يولّد العاكس ثلاث حالات جهد منفصلة عند أطراف المخرجات، مما يؤدي إلى تحسين جودة شكل الموجة وكفاءة تحويل الطاقة. وتلعب نقطة التعادل لناقل التيار المستمر في عاكس NPC دوراً رئيسياً في موازنة جهود المكثفات وضمان التشغيل السليم للعاكس. يظهر مخطط التحكم في العاكس في الشكل 3.
4. محول رفع الجهد من تيار مستمر إلى تيار مستمر (DC-DC Boost Converter)
تقوم محولات التيار المستمر إلى تيار مستمر (DC-DC converters) برفع جهد التيار المستمر الداخل إلى جهد خرج أعلى من خلال التحكم في دورة التشغيل (duty cycle) لجهاز تبديل عالي التردد مع الحفاظ على الكفاءة. وللمحول وضعان للتشغيل؛ فخلال حالة التشغيل (ON state)، يكون الديود في حالة انحياز عكسي، ويقوم مصدر الدخل بتزويد المحث بالطاقة التي يخزنها على شكل مجال مغناطيسي، بينما يقوم مكثف الخرج بتزويد الحمل بالطاقة. أما خلال حالة الإيقاف (OFF state)، فيصبح الديود في حالة انحياز أمامي، وهنا تتحرر الطاقة المخزنة في المحث عبر الديود إلى المكثف المتصل بالخرج والحمل، وتتحد مع مصدر الدخل لإنتاج جهد خرج أكبر من جهد الدخل. ومن خلال تغيير دورة تشغيل تعديل عرض النبضة (PWM)، يتغير جهد خرج المحول وفقاً للعلاقة
، حيث تمثل D دورة التشغيل. وتُستخدم في الأنظمة الشمسية، وخلايا الوقود، والمركبات الكهربائية، وأنظمة التخزين.
مخطط تحكم بالشبكة العصبية الاصطناعية لتتبع أقصى نقطة قدرة
في هذا العمل، تم اختيار تتبع نقطة الطاقة القصوى (MPPT) القائم على الشبكة العصبونية الاصطناعية (ANN) لجمع أقصى قدر من الطاقة المولدة من النظام الشمسي.17تَمَّ توليد مجموعات بيانات التشعيع ودرجة الحرارة المستخدمة في تدريب الشبكة العصبية الاصطناعية (ANN) في برنامج MATLAB/Simulink تحت ظروف تشغيل بيئية متغيرة، بما يتوافق مع مستويات إشعاع شمسي تتراوح بين 80 و100 واط/م²2 ونطاقات درجات حرارة تتراوح بين 20-25 °C، مع دورة تشغيل (D). وتُدرَّب الشبكة العصبية الاصطناعية (ANN) لتقليل متوسط مربع الخطأ (MSE) بين المخرجات التي تم الحصول عليها والمخرجات المطلوبة (أوp,
لاستخراج أقصى قدرة من نظام الخلايا الكهروضوئية ثنائية الوجه. وهنا، تُستخدم شدة الإشعاع الشمسي ودرجة الحرارة كمدخلات للشبكة العصبية الاصطناعية (ANN)، بينما تكون دورة التشغيل (duty cycle) هي المخرج للتحكم في محول الرفع (boost converter). وبنية نموذج الشبكة العصبية الاصطناعية المطور لتتبع نقطة الطاقة القصوى (MPPT)، كما هو موضح في الشكل 4، يتم أخذها في الاعتبار في هذا العمل. الجدول التكميلي 1 يعرض المزايا والعيوب للطريقة المقترحة مقارنةً بالطرق القياسية الأخرى.
تكون الشبكة العصبية الاصطناعية (ANN) من ثلاث طبقات رئيسية: طبقة المخرجات (OPL)، وطبقة المدخلات (IPL)، وطبقة خفية (HIL). ويساعد ذلك في نقل البيانات بين طبقة المدخلات (IPL) والطبقة الخفية (HIL). وبعد ذلك، يتم إنتاج البيانات بواسطة الأوزان في الروابط التي تربط بين طبقة المدخلات (IPL) والطبقة الخفية (HIL). وفي هذا السياق، تُنفذ الحسابات مع تطبيق انحياز معين على متغير الطبقة الخفية (HL)، وتتراكم النتائج الناتجة في متغير طبقة المخرجات (OPL). وقد تم هنا اختيار الشبكة العصبية الاصطناعية من نوع LMBP [17]. ويتم ضبط أوزان الاتصال أثناء التدريب من خلال قياس الخطأ لتحقيق المخرج المستهدف. وهنا، يُستخدم خوارزم LMBP لتدريب الشبكات العصبية الاصطناعية (ANNs) باستخدام متوسط مربع الخطأ (MSE) كدالة للأداء. ويستخدم خوارزم LMBP المشتقات المستمدة لتحديث الأوزان، والتي تظهر خصائص التعلم الفعال والتقارب المتسارع
تحتوي كل عصبونة في شبكة البيرسيبترون متعددة الطبقات على دالة جمع ودالة تنشيط. ومع ذلك، توجد أوزان عددية معينة (wpk) تربط هذه العصبونات عبر المستويات. يتم ضرب المدخلات في الأوزان كما هو محدد في المعادلة (5). وتُستخدم دالة السينيج (sigmoid) غير الخطية مع تقليل متوسط مربع الخطأ (MSE) الموضح في المعادلة (6).
(5)
(6)
قام خوارزمية SBOA بتحسين معلمات المرشح والتحكم للنظام المقترح
تلعب الخوارزميات الميتا-هيو رستيكية (metaheuristic) في هذه الأيام دوراً رئيسياً في حل المشكلات الهندسية. يوضح الشكل 5 تصنيفاً للخوارزميات. تحاكي خوارزمية SBOA23 سلوك طائر السكرتير (SB) للبقاء على قيد الحياة في موطنه الطبيعي. وقد تم نمذجة مرحلتي الاستكشاف والاستغلال في سلوك الصيد لدى طائر السكرتير في خوارزمية SBOA؛ حيث تحاكي مرحلة الاستكشاف في الخوارزمية سلوك طائر السكرتير في صيد الثعابين، بينما تحاكي مرحلة الاستغلال سلوكها في تفادي المفترسات مثل العقبان. تبدأ عملية الحل بتوليد قيم عشوائية ضمن الحدود المخصصة لكل طائر سكرتير في المجتمع، ومن ثم حساب قيمة الدالة الهدف عن طريق تشغيل نموذج Simulink، مع اعتبار قيمة كل طائر سكرتير كمعلمة تصميم. وبناءً على القيمة الهدف، يغير كل طائر سكرتير موضعه لتمثيل حل أفضل. ويناقش فيما يلي النمذجة الرياضية لخوارزمية SBOA:
مرحلة التحضير الأولية
أولاً، من الضروري تحديد المحاليل الأولية التي استُخدمت لبدء البحث عن تقليل نموذجي لدالة الهدف F(Y). في هذه الحالة، يتم تكوين المجتمع العشوائي الأولي Y = [Y1, Y2, ......YN] من SB‘s عبر N من المحاليل الأولية. تقوم المعادلة (7) بتهيئة المجتمع، ويمثل Yi حلاً ما.
(7)
حيث يُرمز للحدود الدنيا والعليا لمتغيرات الاختيار بـ lb و ub. ويمثل r رقماً عشوائياً في النطاق [0, 1]. بينما تمثل N بُعد المسألة. علاوة على ذلك، تُستخدم قيمة لياقة الحل Fi = F(Yi) لقياس جودة Yi.
استراتيجية الصيد لـ SB’s
تعد عمليات البحث عن الفريسة وأكلها ومهاجمتها المراحل الثلاث الرئيسية لصيد طائر الـ SB. وقد قُسمت عملية الصيد بأكملها إلى ثلاثة فترات زمنية متساوية، t < 1/3T، و 1/3T < t < 2/3T، و 2/3T < t < T. حيث تمثل t التكرار الحالي، وتمثل T الحد الأقصى لعدد التكرارات. وتتوافق هذه الفترات مع المراحل الثلاث لافتراس الـ SB: البحث، والاستهلاك، ومهاجمة الفريسة. وتستند هذه التقسيمات إلى الإحصاءات البيولوجية لمراحل صيد الطائر والمدد الزمنية لكل مرحلة. وبناءً على ذلك، يتم نمذجة كل مرحلة من مراحل SBOA على النحو التالي:
البحث عن الفريسة
يجب على SB صيد الفريسة من مسافة آمنة خلال هذه المرحلة. ولجمع بيانات كافية عبر منطقة البحث بأكملها، تتطلب الخطوة الأولى في خوارزميات التحسين استكشافاً أقوى. ومن خلال استخدام مواقع اثنين من SB الآخرين كدليل، يمكن لـ SB استقصاء مناطق محتملة إضافية. وبالتالي، تم إدخال عمليات طفرة تفاضلية لزيادة تنوع الخوارزمية بشكل أكبر. عندما تكون t أقل من 1/3 من T، يتم استخدام المعادلة (8) لتحديث موضع كل فرد، حيث سيقوم بالبحث عن الفريسة واستهلاكها ومهاجمتها. وبناءً على ذلك، يتم نمذجة كل مرحلة من مراحل SBOA على النحو التالي:
(8)
حيث تمثل Yi(t) الموضع الحالي لـ SB رقم ith، وتعطي Ynewi(t) الحل المحدث، وقد تم اختيار فردين Yi(t) و Yr2(t) عشوائيًا من المجتمع الحالي. ومن أجل استراتيجية اختيار المعلمات، تم اختيار متجه عشوائي R1 بحجم 1 × N عشوائيًا من [0, 1] لتحقيق توازن فعال بين مرحلتي الاستغلال والاستكشاف. ولا يتم اعتبار الحل المحدث إلا إذا أدى إلى تحسين قيمة الهدف مقارنة بالحل الذي تم الحصول عليه سابقًا.
افتراس الفريسة
تستخدم SB رشاقة فائقة وتتحرك لتحوم حول الثعبان بعد رصد فريسة محتملة. وسيتم استغلال صبر الفريسة لخفض دفاعاتها من خلال مراقبة الخصوم وإغرائهم أثناء دورانهم. وتقوم SB الأخرى بتعديل وضعياتها بالقرب من الفريسة. وبفضل هذا النهج، ستزداد فرص نجاح الصيد بشكل ملحوظ. وتوضح المعادلة (9) كيفية استخدام الحركة البراونية (B) لمحاكاة الحركة العشوائية لـ SB عندما يكون 1/3 T < t < 2/3 T.
(9)
يتم تنفيذ B في برنامج MATLAB باستخدام التوزيع الطبيعي القياسي بواسطة المتجه r(1,M) الذي يتم توليده عشوائياً. بعد ذلك، تستخدم SB المعادلة (10) لتحديث مواقعها.
(10)
حيث يمثل Ybest(t) الحل الأمثل، ويتم اختيار عامل القياس الأسي exp((t/T)4) لزيادة قدرة الاستغلال في كل تكرار.
مهاجمة الفريسة
بعد الاستهلاك المستمر، ستنفد طاقة الضحية. ويجب على طيور السكرتير شن الهجوم الآن. هنا، يتم محاكاة طرق هجوم مختلفة، مثل الخطوات المتواصلة / القفزات الطويلة المتقطعة خلال فترة زمنية قصيرة، باستخدام تقنية رحلة ليفي (Lévy flight). وتوصف هذه الخصائص بواسطة المعادلتين (1) و (12). وتكون الخيارات المرشحة قريبة من الإجابة المثلى في هذه اللحظة لأن طيور السكرتير تقترب بسرعة من الفريسة. سيتم تنفيذ هذه الخطة عندما تكون t > 2/3T.
(11)
(12)
أين، تصفية السوائل (Liquid Filtration) تمثل استراتيجية طيران ليفي و
، وقيم ثابتة لـ s = 0.01 و φ = 1.5. μ و υ هي قيم اعتباطية تقع في النطاق [0، 1]. حيث، τ يرمز إلى غاما و η = 0.5، حيث اعتُمد معامل الوزن 0.5، كما هو مقترح في صياغة SBOA الأصلية، لتنظيم حجم الخطوة والحفاظ على التوازن بين القفزات الاستكشافية الكبيرة والتحسين المحلي. الجدول 4 تعرض القيم المستخدمة للشبكة العصبية الاصطناعية (ANN) وخوارزمية SB في النظام المطور.
استراتيجية الهروب لـ SB
عندما تطارد طيور السكرتير فرائس أخرى في البرية، فإنها تخاطر بأن تصبح هي نفسها مطاردة. وتعد النسور والصقور والثعالب وابن آوى الخصوم الرئيسيين الذين يجب أن تتعامل معهم. لذا، يتعين عليها استخدام مجموعة متنوعة من تقنيات التهرب لحماية نفسها أو فرائسها عند استشعار الخطر. تحاكي هذه الخوارزمية طرق الهروب من خلال نمذجة أوضاع الركض (D2) والتمويه (D1). فعند مواجهة الخصوم، تندمج طيور السكرتير في البداية مع محيطها للبقاء آمنة؛ وهنا تسعى طيور السكرتير إلى تجنب الوصول إلى القيم المثلى المحلية في الخوارزميات عن طريق تحديث مواقعها حول الفريسة (أفضل فرد). وإذا لم تتمكن من التهرب من العدو، فإنها تستخدم تكتيكات الطيران أو الركض السريع للبقاء آمنة. وللمرجعية، يتم اختيار فرد عشوائي Yrand كقائد لتجنب التقيد بقيمة مثلى محلية. وتستخدم طيور السكرتير المعادلتين (13) و(14) لتحديث مواقعها. يظهر المخطط الانسيابي لخوارزمية SBOA في الشكل 6.
(13)
(14)
تمثيل متغيرات التصميم
تضمن متغيرات المشكلة في هذه الدراسة معاملات التصميم، مثل مكاسب وحدة التحكم PI للمحولات DC-DC، بالإضافة إلى مقاومات ومحثات المرشحات والمحولات. وتوضح المعادلة (15) تمثيل متغيرات التحكم. أما حدود معاملات التحكم فهي مدرجة في الجدول التكميلي 2.
(15)
دالة اللياقة (FF)
تم اختيار تقليل THD كدالة هدف (Obj) في هذه الدراسة. وتكون القيمة القصوى لـ FF هي القيمة الدنيا لـ Obj، كما هو موضح في المعادلتين (16) و(17).
(16)
(17)