In dieser Studie wird ein eigenständiges, durch ein bifaziales Photovoltaik- und Batteriespeichersystem (BESS) gespeistes Elektrofahrzeug-Ladesystem entwickelt, wie in Abbildung 1 dargestellt, unter Verwendung einer durch ein künstliches neuronales Netz (ANN) gesteuerten MPPT-Methode. Zur Modellierung des bifazialen Photovoltaik-Arrays wird eine Einzel-Dioden-Äquivalentschaltung verwendet, die Temperatur, Serienwiderstand, Shunt-Widerstand, Sonneneinstrahlung und Bestrahlung der Rückseite berücksichtigt. Die Ausgangsleistung der Photovoltaikanlage wird über einen Gleichstrom-Gleichstrom-Boost-Wandler (DC–DC) zur Elektrofahrzeug-Ladestation geleitet, während das BESS die Regelung der Gleichstromzwischenkreisspannung (DC-Bus) unterstützt. Der Secretary-Bird-Optimierungsalgorithmus (SBOA) wird zur Optimierung der PI-Reglerverstärkungen und der Filterparameter verwendet. Die Modellierung der Systemkomponenten sowie die vorgeschlagene Methodik werden in den folgenden Unterkapiteln beschrieben.
Modellierung von Komponenten
1. Bifaziales PV-System
Die PV-Module tragen dazu bei, die erforderliche Strom- und Spannungsausgabe zu erreichen. In der vorliegenden Studie wurden zwei monofaziale PV-Module verwendet, um die Energieerzeugungscharakteristika bifazialer Anlagen unter verschiedenen Oberflächen-Albedo-Bedingungen zu simulieren9. Die bifaziale Ertragssteigerung der Bestrahlungsstärke BGg wird durch Gleichung (1) bestimmt
(1)
Hier ist die Frontbestrahlung Gf und die Hinterseitenbestrahlungsstärke Gr. Der maximale bifaziale Photovoltaikstrom beträgt Imp, die Spannung Vmp. Schließlich ergibt sich die Leistung des bifazialen Systems mit der gewählten Anzahl an in Reihe Ns und parallel Np geschalteten Modulen gemäß Gleichung (2). Die Regelung der Photovoltaik ist in Abbildung 2 dargestellt.
(2)
2. Batteriespeichersystem (BSS)
Batterien sind eine entscheidende Möglichkeit zur Speicherung elektrischer Energie18. In der vorliegenden Arbeit werden Lithium-Ionen-Batterien als Energiespeichersystem verwendet. Unter den verfügbaren wiederaufladbaren Batterietechnologien liegt der Fokus dieser Arbeit auf Lithium-Ionen-Batterien, die aufgrund ihrer Effizienz und Energiedichte als optimal für Elektrofahrzeuge gelten. Die Gleichungen (3) und (4) beschreiben die Lithium-Ionen-Batterie mit einschränkenden Bedingungen.
(3)
Die Spezifikationen der PV-Anlage, der Speicherbatterie und des Umrichters sind in Tabelle 2 zusammengefasst, während Tabelle 3 die Verteilung des Leistungsflusses zwischen der PV-Generatoranlage, der Energiespeichereinheit und den angeschlossenen Lasten darstellt.
(4)
3. Gleichstrom–Wechselstrom-Wandler
Der Neutralpunkt-gekoppelte (NPC-)Wechselrichter wird in modernen Stromumwandlungssystemen für erneuerbare Energien weit verbreitet eingesetzt. Im Vergleich zur herkömmlichen Zweistufenstruktur bietet die NPC-Konfiguration eine geringere Gesamtharmonische Verzerrung (THD), geringere elektromagnetische Störungen und eine niedrigere Spannungsbelastung der Halbleiterbauelemente. Der Gleichstromzwischenkreis dieser Topologie wird durch gepaarte Kondensatoren und Klemmdioden, die den Neutralpunkt bilden, in drei Potentialstufen +Vdc/2, 0 und -Vdc/2 unterteilt. Durch geeignete Steuerung der Schaltvorgänge erzeugt der Wechselrichter drei diskrete Spannungszustände an den Ausgangsklemmen, wodurch die Wellenformqualität und die Effizienz der Leistungsumwandlung verbessert werden. Der Neutralpunkt des Gleichstrombusses in einem NPC-Wechselrichter spielt eine entscheidende Rolle bei der Ausgleichung der Kondensatorspannungen und der Sicherstellung eines ordnungsgemäßen Betriebs des Wechselrichters. Das Regelkonzept des Wechselrichters ist in Abbildung 3 dargestellt.
4. Gleichspannungswandler mit Spannungsanhebung
DC-DC-Wandler erhöhen die Eingangs-Gleichspannung auf eine höhere Ausgangsspannung, indem sie den Tastgrad eines hochfrequent schaltenden Bauelements steuern, wobei sie gleichzeitig die Effizienz aufrechterhalten. Er verfügt über zwei Betriebsarten. Im eingeschalteten Zustand (ON-Zustand) ist die Diode sperrvorgespannt, und die Eingangsquelle liefert Energie an die Spule, die diese als magnetisches Feld speichert, während der Ausgangskondensator Energie an die Last abgibt. Im ausgeschalteten Zustand (OFF-Zustand) wird die Diode in Durchlassrichtung gepolt. Hier wird die in der Spule gespeicherte Energie über die Diode an den am Ausgang angeschlossenen Kondensator und die Last abgegeben und kombiniert sich mit der Eingangsquelle, um eine Ausgangsspannung zu erzeugen, die höher ist als die Eingangsspannung. Durch Variation des PWM-Tastgrades ändert sich die Ausgangsspannung des Wandlers gemäß der Beziehung
, wobei D den Tastgrad bezeichnet. Er wird in Solar-, Brennstoffzellen-, Elektrofahrzeug- und Speichersystemen eingesetzt.
ANN-Regelungsschema für MPPT
In dieser Arbeit wurde ein auf einem künstlichen neuronalen Netzwerk (ANN) basierendes MPPT-Verfahren gewählt, um die maximale Leistung des Solarsystems zu erfassen17. Die für das Training des ANN verwendeten Datensätze zu Bestrahlungsstärke und Temperatur wurden in MATLAB/Simulink unter variierenden Umweltbedingungen erzeugt, die Solarstrahlungsintensitäten von 800–1000 W/m2 und Temperaturbereichen von 20–25 °C entsprechen, jeweils mit einem Tastverhältnis (D). Das ANN wird trainiert, den mittleren quadratischen Fehler (MSE) zwischen den erzielten und erforderlichen Ausgabewerten (Op,
) zu minimieren, um die maximale Leistung aus dem bifazialen Photovoltaiksystem zu gewinnen. Hierbei werden Solarbestrahlung und Temperatur als Eingangsgrößen für das ANN verwendet, während das Tastverhältnis die Ausgangsgröße zur Steuerung des Gleichspannungswandlers darstellt. Die Struktur des entwickelten ANN-Modells für MPPT, wie sie in Abbildung 4 dargestellt ist, wird in der vorliegenden Arbeit berücksichtigt. Ergänzende Tabelle 1 zeigt die Vor- und Nachteile der vorgeschlagenen Methode im Vergleich zu anderen Standardverfahren auf.
Ein künstliches neuronales Netz (ANN) besteht aus drei Hauptebenen: einer Ausgabeebene (OPL), einer Eingabeebene (IPL) und einer versteckten Ebene (HIL). Dies unterstützt den Datentransfer zwischen der IPL und der HIL. Anschließend werden die Daten durch die Gewichtungen in den Verbindungen zwischen der IPL und der HIL erzeugt. In diesem Zusammenhang werden Berechnungen mit einer bestimmten Vorspannung (Bias) an der HIL-Variable durchgeführt, und die resultierenden Ergebnisse werden in der OPL-Variable gesammelt. Hier wird ein ANN vom Typ LMBP [17] gewählt. Die Verbindungsgewichtungen werden während des Trainingsprozesses angepasst, indem der Fehler zur Erreichung des gewünschten Ausgabewerts gemessen wird. Dabei wird LMBP verwendet, um ANNs mit dem mittleren quadratischen Fehler (MSE) als Leistungsfunktion zu trainieren. Der LMBP-Algorithmus nutzt die abgeleiteten Gradienten, um die Gewichtungen zu aktualisieren, wodurch sich effizientes Lernen und eine beschleunigte Konvergenz zeigen.
Jedes Neuron in einem mehrschichtigen Perzeptron-Netzwerk verfügt über eine Summations- und eine Aktivierungsfunktion. Dennoch gibt es bestimmte numerische Gewichtungen (wpk), die diese Neuronen über die Schichten hinweg verbinden. Wenn Eingaben mit den Gewichtungen multipliziert werden, wie in Gleichung (5) angegeben. Die nichtlineare Sigmoidfunktion wird unter Berücksichtigung der in Gleichung (6) angegebenen MSE-Minimierung betrachtet.
(5)
(6)
SBOA optimierte die Filter- und Regelparameter des vorgeschlagenen Systems
Heutzutage spielen metaheuristische Algorithmen eine entscheidende Rolle bei der Lösung technischer Probleme. Abbildung 5 zeigt eine Klassifizierung von Algorithmen. Der SBOA-Algorithmus23 ahmt das Verhalten des Wappenvogels (SB) nach, um in seinem natürlichen Lebensraum zu überleben. Die Phasen der Erkundung und Nutzung des Jagdverhaltens des Wappenvogels werden im SBOA abgebildet. Die Erkundungsphase des Algorithmus modelliert das Verhalten des Wappenvogels beim Fang von Schlangen, während die Nutzungsphase sein Verhalten beim Ausweichen vor Raubtieren wie Adlern nachahmt. Der Lösungsprozess beginnt damit, zufällige Werte innerhalb der jeweiligen Grenzen für jeden Wappenvogel in der Population zu erzeugen, und berechnet den Wert der Zielfunktion durch Ausführen des Simulink-Modells, wobei jeder Wert des Wappenvogels als Entwurfsparameter behandelt wird. Basierend auf dem Zielfunktionswert ändert jeder Wappenvogel seine Position, um eine bessere Lösung darzustellen. Die mathematische Modellierung des SBOA wird im Folgenden erläutert:
Anfangsphase der Vorbereitung
Zunächst ist es notwendig, die anfänglichen Lösungen zu identifizieren, die verwendet wurden, um die Suche nach einer typischen Minimierung der Zielfunktion F(Y) zu beginnen. In diesem Fall wird die anfängliche zufällige Population Y = [Y1, Y2, ......YN] von SBs durch N anfängliche Lösungen gebildet. Gleichung (7) initialisiert die Population, wobei Yi eine Lösung darstellt.
(7)
Dabei werden die unteren und oberen Grenzen der Entscheidungsvariablen mit lb und ub bezeichnet. In [0, 1] ist r eine Zufallszahl. N ist die Dimension des Problems. Außerdem wird der Fitnesswert der Lösung Fi = F(Yi) verwendet, um die Qualität von Yi zu quantifizieren.
Jagdstrategie von SBs
Das Auffinden, Fressen und Angreifen der Beute sind die drei Hauptphasen der SB-Jagd. Der gesamte Jagdprozess wurde in drei gleich lange Zeitintervalle unterteilt: t < 1/3T, 1/3T < t < 2/3T und 2/3T < t < T. Dabei ist t die aktuelle Iteration und T die maximale Anzahl an Iterationen. Diese Intervalle entsprechen den drei Phasen der Räuberjagd des SB: der Suche, der Nahrungsaufnahme und dem Angriff auf die Beute. Diese Einteilung basiert auf biologischen Statistiken der Jagdphasen des Vogels sowie auf den zeitlichen Dauern jeder Phase. Daher wird jede Phase des SBOA wie folgt modelliert:
Die Suche nach Beute
Der SB muss in dieser Phase Beute aus sicherer Entfernung jagen. Um ausreichende Daten über die gesamte Suchregion zu sammeln, ist im ersten Schritt von Optimierungsalgorithmen eine stärkere Erkundung erforderlich. Mithilfe der Positionen der beiden anderen SBs als Leitfaden kann der SB zusätzliche mögliche Regionen untersuchen. Um die Diversität des Algorithmus weiter zu erhöhen, werden daher differenzielle Mutationsprozesse eingeführt. Wenn t kleiner als ein Drittel von T ist, wird Gleichung (8) verwendet, um die Position jedes Individuums zu aktualisieren. Er wird Beute suchen, konsumieren und angreifen. Dementsprechend wird jede Phase des SBOA wie folgt modelliert:
(8)
Dabei ist Yi(t) die aktuelle Position des iten SB, Ynewi(t) liefert die aktualisierte Lösung; zwei Individuen Yi(t), Yr2(t) wurden zufällig aus der vorhandenen Population ausgewählt. Ein Zufallsvektor R1 für die Parameterauswahlstrategie mit 1 × N Elementen wurde zufällig aus [0, 1] gewählt, um ein wirksames Gleichgewicht zwischen den Phasen der Nutzung und Erkundung zu erreichen. Die aktualisierte Lösung wird nur berücksichtigt, wenn sie einen verbesserten Zielfunktionswert im Vergleich zur zuvor erhaltenen Lösung ergibt.
Beute verzehren
SB nutzt bewegliche Agilität und bewegt sich, um nach dem Erkennen möglicher Beute um die Schlange zu schweben. Die Geduld der Beute wird genutzt, um deren Verteidigung zu senken, indem Gegner beim Umkreisen beobachtet und gereizt werden. Andere SB passen ihre Haltung in der Nähe der Beute an. Der Jagderfolg wird mit diesem Ansatz deutlich gesteigert. Gleichung (9) veranschaulicht, wie die Brownsche Bewegung (B) verwendet wird, um die zufällige Bewegung des SB zu simulieren, wenn 1/3 T < t < 2/3 T.
(9)
Dabei wird B in MATLAB unter Verwendung der Standardnormalverteilung durch den zufällig generierten Vektor r(1,M) implementiert. Anschließend verwenden die SBs Gleichung (10), um ihre Positionen zu aktualisieren.
(10)
Dabei ist Ybest(t) die beste Lösung, und der exponentielle Skalierungsfaktor exp((t/T)4) wird gewählt, um die Exploitationskapazität in jeder Iteration zu erhöhen.
Beute angreifen
Nach anhaltendem Verbrauch wird die Beute ihren Energievorrat erschöpfen. Die Sekretärsvögel sollten nun den Angriff starten. Hier werden unterschiedliche Angriffsmethoden, wie andauernde Schritte oder sporadische weite Sprünge über einen kurzen Zeitraum, mithilfe der Lévy-Flug-Technik simuliert. Die Eigenschaften werden durch die Gleichungen (11) und (12) beschrieben. Die Kandidatenoptionen liegen momentan nahe an der optimalen Lösung, da sich die SBs rasch der Beute nähern. Dieser Plan wird umgesetzt, wenn t > 2/3T ist.
(11)
(12)
Dabei steht LF für die Levy-Flug-Strategie und
, wobei s = 0,01 und φ = 1,5 feste Werte sind. μ und υ sind beliebige Werte im Intervall [0, 1]. Hierbei bezeichnet τ Gamma und η = 0,5; der Gewichtungsfaktor 0,5 wurde gemäß der ursprünglichen SBOA-Formulierung übernommen, um die Schrittweite zu regulieren und ein Gleichgewicht zwischen großen explorativen Sprüngen und lokaler Feinabstimmung zu gewährleisten. Tabelle 4 zeigt die im entwickelten System verwendeten Werte für das KNN und den SB-Algorithmus.
Fluchtstrategie für SB
Wenn sie in freier Wildbahn andere Beutetiere jagen, laufen SBs Gefahr, selbst gejagt zu werden. Adler, Habichte, Füchse und Schakale sind die Hauptfeinde, mit denen sie sich auseinandersetzen müssen. Sie müssen eine Vielzahl von Ausweichtechniken anwenden, um sich selbst oder ihre Beute zu schützen, sobald sie Gefahr wahrnehmen. Dieser Algorithmus simuliert die Fluchtmethoden, indem er Laufmodi (D2) und Tarnung (D1) modelliert. Wenn sie Feinden gegenüberstehen, verschmelzen Sekretärsvögel zunächst mit ihrer Umgebung, um sicher zu bleiben. Hier bemühen sich Sekretärsvögel, lokale Optima im Algorithmus zu vermeiden, indem sie ihre Positionen um die Beute (das beste Individuum) herum aktualisieren. Falls sie dem Feind nicht ausweichen können, wenden sie Flug- oder Schnelllauf-Strategien an, um sicher zu bleiben. Zu Referenzzwecken wird ein zufälliges Individuum Yrand als Leiter gewählt, um zu verhindern, dass man in einem lokalen Optimum feststeckt. Sekretärsvögel nutzen die Gleichungen (13) und (14), um ihre Positionen zu aktualisieren. Das Ablaufdiagramm des SBOA ist in Abbildung 6 dargestellt.
(13)
(14)
Darstellung der Designvariablen
Zu den Problemvariablen dieser Studie gehören die Entwurfsparameter, wie die PI-Reglerverstärkungen von Gleichspannungswandlern, sowie die Widerstände und Induktivitäten von Filtern und Wandlern. Gleichung (15) zeigt die Darstellung der Steuervariablen. Die Grenzen der Steuerungsparameter sind in Ergänzungstabelle 2 aufgeführt.
(15)
Fitness-Funktion (FF)
Die Minimierung der THD wird in dieser Studie als Zielfunktion (Obj) gewählt. Der Maximalwert von FF entspricht dem Minimalwert von Obj, wie in den Gleichungen (16) und (17) angegeben.
(16)
(17)