W niniejszej pracy opracowano autonomiczny system ładowania pojazdów elektrycznych (EV) zasilany dwustronnymi modułami fotowoltaicznymi (PV) oraz systemem magazynowania energii z baterii (BESS), jak pokazano na Rysunku 1, z wykorzystaniem metody MPPT sterowanej przez sztuczną sieć neuronową (ANN). Do modelowania dwustronnego układu PV zastosowano równoważny obwód jednodiodowy, uwzględniający temperaturę, rezystancję szeregową, rezystancję równoległą, natężenie promieniowania słonecznego oraz natężenie promieniowania na tylnej stronie. Moc z układu PV jest dostarczana do stacji ładowania EV za pomocą przetwornicy DC–DC typu boost, podczas gdy BESS wspomaga regulację napięcia szyny DC. Do optymalizacji wzmocnień regulatora PI oraz parametrów filtrów wykorzystano algorytm optymalizacji Secretary Bird (SBOA). Modelowanie komponentów systemu oraz proponowana metodologia zostały opisane w następnych podrozdziałach.
Modelowanie komponentów
1. Dwustronny system fotowoltaiczny (PV)
Moduły PV pomagają w osiągnięciu wymaganego natężenia prądu i napięcia wyjściowego. W niniejszym badaniu wykorzystano dwa jednostronne moduły PV do emulacji charakterystyki generowania energii przez macierze dwustronne w różnych warunkach albedo powierzchni9. Zysk dwustronny z promieniowania BGg jest określany przez równanie (1)
(1)
W tym przypadku natężenie promieniowania przedniego to Gf, a natężenie promieniowania tylnego to Gr. Maksymalny prąd dwustronnego modułu PV to Imp, a napięcie to Vmp. Ostatecznie moc wyjściowa generowana przez system dwustronny z wybraną liczbą modułów połączonych szeregowo Ns oraz równolegle Np jest określona przez równanie (2). Układ sterowania PV przedstawiono na Rysunku 2
(2)
2. System magazynowania energii w akumulatorach (BSS)
Baterie są kluczowym sposobem przechowywania energii elektrycznej18. W niniejszej pracy jako system magazynowania energii wykorzystano baterie Li-ion. Spośród dostępnych opcji akumulatorów wielokrotnego ładowania, praca ta koncentruje się na bateriach Li-ion, które ze względu na swoją wydajność i gęstość mocy uznaje się za optymalne dla pojazdów elektrycznych (EV). Równania (3) i (4) opisują baterię Li-ion z uwzględnieniem ograniczeń.
(3)
Specyfikacje PV, akumulatora i konwertera zostały podsumowane w Tabeli 2, natomiast Tabela 3 przedstawia alokację przepływu mocy pomiędzy tablicą PV, jednostką magazynowania energii a podłączonymi obciążeniami.
(4)
3. Przetwornica DC-AC
Inwerter z zaciskaniem punktu neutralnego (NPC) jest powszechnie stosowany w nowoczesnych systemach konwersji energii dla odnawialnych źródeł energii. W porównaniu z konwencjonalną strukturą dwupoziomową, konfiguracja NPC oferuje zredukowane całkowite harmoniczne zniekształcenia (THD), niższe zakłócenia elektromagnetyczne oraz mniejsze naprężenia napięciowe w urządzeniach półprzewodnikowych. Szyna DC tej topologii jest podzielona na trzy poziomy potencjału +Vdc/2, 0 i -Vdc/2 za pomocą sparowanych kondensatorów i diod zaciskowych tworzących punkt neutralny. Poprzez odpowiednie sterowanie przełączaniem, inwerter generuje trzy dyskretne stany napięcia na zaciskach wyjściowych, co poprawia jakość przebiegu i sprawność konwersji mocy. Punkt neutralny szyny DC w inwerterze NPC odgrywa kluczową rolę w równoważeniu napięć kondensatorów i zapewnianiu prawidłowego działania inwertera. Schemat sterowania inwertera przedstawiono na Rysunku 3.
4. Przetwornica step-up DC-DC
Przetwornice DC-DC podwyższają wejściowe napięcie stałego prądu do wyższego napięcia wyjściowego poprzez sterowanie cyklem pracy urządzenia przełączającego o wysokiej częstotliwości, przy jednoczesnym zachowaniu sprawności. Posiada ona dwa tryby pracy. W stanie włączonym (ON) dioda jest spolaryzowana w kierunku zaporowym, a źródło wejściowe dostarcza energię do cewki, która magazynuje ją w postaci pola magnetycznego, podczas gdy kondensator wyjściowy zasila odbiornik. W stanie wyłączonym (OFF) dioda zostaje spolaryzowana w kierunku przewodzenia. Wówczas energia zmagazynowana w cewce jest uwalniana poprzez diodę do kondensatora podłączonego do wyjścia i odbiornika, sumując się z napięciem źródła wejściowego, co pozwala uzyskać napięcie wyjściowe wyższe od wejściowego. Poprzez zmianę wypełnienia sygnału PWM napięcie wyjściowe przetwornicy zmienia się zgodnie z zależnością
, gdzie D oznacza cykl pracy. Przetwornice te są stosowane w systemach solarnych, ogniwach paliwowych, pojazdach elektrycznych (EV) oraz systemach magazynowania energii.
Schemat sterowania ANN dla MPPT
W niniejszej pracy wybrano algorytm MPPT oparty na sztucznej sieci neuronowej (ANN) w celu pozyskania maksymalnej mocy generowanej przez system solarny17Zbiory danych dotyczących promieniowania i temperatury wykorzystane do trenowania sieci ANN zostały wygenerowane w środowisku MATLAB/Simulink przy zmiennych warunkach pracy środowiskowej, odpowiadających poziomom nasłonecznienia od 800 do 1000 W/m²2 oraz zakresy temperatur od 20–25 °Cz cyklem pracy (D). Sieć ANN jest trenowana w celu zminimalizowania błędu MSE pomiędzy uzyskanymi a wymaganymi wartościami wyjściowymi (Op,
) w celu wydobycia maksymalnej mocy z bifacjalnego systemu fotowoltaicznego. W tym przypadku natężenie promieniowania słonecznego oraz temperatura służą jako dane wejściowe dla sztucznej sieci neuronowej (ANN), a cykl pracy jest sygnałem wyjściowym służącym do sterowania przetwornicą step-up (boost). Struktura opracowanego modelu ANN dla śledzenia punktu maksymalnej mocy (MPPT), przedstawiona w Rycina 4zostało uwzględnione w pracy. Tabela uzupełniająca 1 przedstawia zalety i wady proponowanej metody w porównaniu z innymi standardowymi metodami.
Sztuczna sieć neuronowa (ANN) składa się z trzech głównych warstw: warstwy wyjściowej (OPL), warstwy wejściowej (IPL) oraz warstwy ukrytej (HIL). Pozwala to na przekazywanie danych pomiędzy IPL a HIL. Następnie dane są przetwarzane przez wagi w połączeniach łączących IPL i HIL. W tym kontekście obliczenia są wykonywane z zastosowaniem określonego obciążenia (bias) zmiennej HL, a uzyskane wyniki są akumulowane w zmiennej OPL. W niniejszej pracy wybrano ANN typu LMBP [17]. Wagi połączeń są dostrajane podczas trenowania poprzez pomiar błędu w celu osiągnięcia docelowego wyniku wyjściowego. Tutaj LMBP jest wykorzystywany do trenowania sieci ANN z zastosowaniem MSE jako funkcji wydajności. Algorytm LMBP wykorzystuje wyznaczone pochodne do aktualizacji wag, co wykazuje właściwości efektywnego uczenia się i przyspieszonej zbieżności
Każdy neuron w sieci wielowarstwowego perceptrona posiada funkcję sumowania i funkcję aktywacji. Ponadto istnieją określone wagi liczbowe (wpk), które łączą te neurony pomiędzy poszczególnymi poziomami. Wejścia są mnożone przez wagi zgodnie z Równaniem (5). W celu minimalizacji MSE, zgodnie z Równaniem (6), zastosowano nieliniową funkcję sigmoidalną.
(5)
(6)
Algorytm SBOA zoptymalizował parametry filtra i sterowania proponowanego systemu
Obecnie algorytmy metaheurystyczne odgrywają kluczową rolę w rozwiązywaniu problemów inżynieryjnych. Rycina 5 przedstawia klasyfikację algorytmów. Algorytm SBOA23 naśladuje zachowanie ptaka sekretarza w celu przetrwania w jego naturalnym środowisku. W SBOA zamodelowano etapy eksploracji i eksploatacji w zachowaniach łowieckich ptaka sekretarza. Faza eksploracji algorytmu naśladuje zachowanie ptaka sekretarza podczas łapania węży, natomiast faza eksploatacji naśladuje jego zachowanie podczas unikania drapieżników, takich jak orły. Proces poszukiwania rozwiązania rozpoczyna się od generowania wartości losowych w odpowiednich granicach dla każdego osobnika w populacji oraz od obliczenia wartości funkcji celu poprzez uruchomienie modelu Simulink, traktując wartość każdego osobnika jako parametr projektowy. Na podstawie wartości celu każdy osobnik zmienia swoją pozycję, aby reprezentować lepsze rozwiązanie. Modelowanie matematyczne SBOA omówiono poniżej:
Wstępna faza przygotowawcza
W pierwszej kolejności należy zidentyfikować rozwiązania początkowe, które zostały wykorzystane do rozpoczęcia poszukiwań typowego minimum funkcji celu F(Y). W tym przypadku początkowa losowa populacja Y = [Y1, Y2, ......YN] SB jest tworzona poprzez N rozwiązań początkowych. Równanie (7) inicjuje populację, a Yi reprezentuje poszczególne rozwiązanie.
(7)
Gdzie dolne i górne granice zmiennych decyzyjnych są oznaczone jako lb oraz ub. W przedziale [0, 1], r jest liczbą losową. N to wymiar problemu. Ponadto wartość dopasowania rozwiązania Fi = F(Yi) jest wykorzystywana do określenia jakości Yi.
Strategia łowiecka SB
Znalezienie, zjedzenie i zaatakowanie ofiary to trzy główne fazy polowania SB. Cały proces polowania został podzielony na trzy równe odstępy czasu: t < 1/3T, 1/3T < t < 2/3T oraz 2/3T < t < T. Tutaj t oznacza bieżącą iterację, a T maksymalną liczbę iteracji. Przedziały te odpowiadają trzem fazom drapieżnictwa SB: poszukiwaniu, konsumpcji i atakowaniu ofiary. Podział ten opiera się na statystykach biologicznych faz polowania ptaka oraz czasie trwania każdej z tych faz. W związku z tym, każda faza SBOA jest modelowana w następujący sposób:
Poszukiwanie ofiary
W tej fazie SB musi polować na ofiarę z bezpiecznej odległości. Aby zebrać wystarczającą ilość danych z całego obszaru poszukiwań, pierwszy krok w algorytmach optymalizacji wymaga silniejszej eksploracji. Wykorzystując lokalizacje dwóch pozostałych SB jako wskazówkę, SB może zbadać dodatkowe potencjalne obszary. W celu dalszego zwiększenia różnorodności algorytmu wprowadzono zatem procesy mutacji różnicowej. Gdy t jest mniejsze niż 1/3 T, do aktualizacji pozycji każdego osobnika stosuje się równanie (8). Będzie on poszukiwać, konsumować i atakować ofiarę. W konsekwencji, poniżej przedstawiono modelowanie każdej fazy SBOA:
(8)
Gdzie Yi(t) to aktualna pozycja itej SB, Ynewi(t) przedstawia zaktualizowane rozwiązanie, a dwie jednostki Yi(t) i Yr2(t) zostały randomly wybrane z istniejącej populacji. Losowy wektor dla strategii wyboru parametrów, R1 o wymiarach 1 × N, został wybrany losowo z przedziału [0, 1], aby zapewnić efektywną równowagę między fazami eksploatacji i eksploracji. Zaktualizowane rozwiązanie jest brane pod uwagę tylko wtedy, gdy prowadzi do poprawy wartości funkcji celu w porównaniu z wcześniej uzyskanym rozwiązaniem.
Konsumpcja ofiary
SB wykorzystuje zwinną zwinność i porusza się, by krążyć wokół węża po dostrzeżeniu potencjalnej ofiary. Cierpliwość ofiary zostanie wykorzystana do obniżenia jej zapór obronnych poprzez obserwację i wabić przeciwników podczas ich krążenia. Inne SB dostosowują swoją postawę w pobliżu ofiary. Dzięki takiemu podejściu sukces w polowaniu zostanie znacznie zwiększony. Równanie (9) ilustruje, w jaki sposób ruchy Browna (B) są wykorzystywane do symulowania przypadkowego ruchu SB, gdy 1/3 T < t < 2/3 T.
(9)
Wartość B jest implementowana w programie MATLAB przy użyciu standardowego rozkładu normalnego za pomocą losowo generowanego wektora r(1,M). Następnie SB wykorzystują równanie (10) do aktualizacji swoich pozycji.
(10)
Gdzie Ybest(t) jest najlepszym rozwiązaniem, a wykładniczy czynnik skalujący exp((t/T)4) został wybrany w celu zwiększenia zdolności eksploatacji w każdej iteracji.
Atakowanie ofiary
Po ciągłym zużyciu energii ofiara zostanie z niej całkowicie pozbawiona. W tym momencie ptaki sekretarze powinny rozpocząć atak. Tutaj, z wykorzystaniem techniki lotu Lévy’ego, symulowane są różne metody ataku, takie jak nieustanne kroki lub sporadyczne długie skoki w krótkim przedziale czasowym. Charakterystyka tych działań jest opisana w równaniach (11) i (12). Obecnie potencjalne opcje znajdują się blisko rozwiązania optymalnego, ponieważ ptaki sekretarze szybko zbliżają się do ofiary. Plan ten zostanie zaimplementowany, gdy t > 2/3T.
(11)
(12)
Gdzie LF reprezentuje strategię lotu Levy'ego i
, s= 0.01 oraz φ = 1.5 są wartościami stałymi. μ i υ to wartości dowolne z przedziału [0, 1]. Gdzie τ oznacza gamma, a η = 0.5; przyjęto współczynnik wagowy 0.5, zgodnie z propozycją w oryginalnej formulacji SBOA, aby regulować wielkość kroku i utrzymać równowagę między dużymi skokami eksploracyjnymi a lokalnym doprecyzowaniem. Tabela 4 przedstawia wartości wykorzystane dla ANN i algorytmu SB w opracowanym systemie.
Strategia ucieczki dla SB
Podczas polowania na inne ofiary w naturze, ptaki sekretarze ryzykują stanie się zwierzyną łowną. Głównymi przeciwnikami, z którymi muszą się zmagać, są orły, jastrzębie, lisy i szakale. Aby chronić siebie lub swoją ofiarę w obliczu niebezpieczeństwa, muszą stosować różnorodne techniki unikania. Niniejszy algorytm symuluje metody ucieczki poprzez modelowanie trybów biegu (D2) i kamuflażu (D1). W starciu z przeciwnikami ptaki sekretarze początkowo wtapiają się w otoczenie, aby zapewnić sobie bezpieczeństwo. W tym procesie ptaki sekretarze dążą do uniknięcia optimów lokalnych w algorytmach poprzez aktualizację swoich pozycji wokół ofiary (najlepszego osobnika). Jeśli nie są w stanie uniknąć wroga, stosują taktykę lotu lub szybkiego biegu, aby przetrwać. Dla odniesienia, jako lider wybierany jest losowy osobnik Yrand, co zapobiega utknięciu w optimum lokalnym. Ptaki sekretarze wykorzystują równania (13) i (14) do aktualizacji swoich pozycji. Schemat blokowy SBOA przedstawiono na Rysunku 6.
(13)
(14)
Reprezentacja zmiennych projektowych
Zmienne problemowe w niniejszym badaniu obejmują parametry projektowe, takie jak wzmocnienia regulatorów PI przetwornic DC-DC, a także rezystancje i indukcyjności filtrów i przetwornic. Równanie (15) przedstawia reprezentację zmiennych sterowania. Zakresy parametrów sterowania wykazano w Uzupełniającej Tabeli 2.
(15)
Funkcja przystosowania (FF)
W niniejszym badaniu jako funkcję celu (Obj) wybrano minimalizację THD. Maksymalna wartość FF odpowiada minimalnej wartości Obj, zgodnie z równaniami (16) i (17).
(16)
(17)