Bu araştırmada, Şekil 1'de gösterildiği gibi, ANN kontrollü bir MPPT yöntemi kullanılarak bağımsız, bifasiyal PV ve batarya enerji depolama sistemli (BESS) beslemeli bir EV şarj sistemi geliştirilmiştir. Bifasiyal PV dizisini modellemek için sıcaklık, seri direnç, paralel direnç, güneş ışınımı ve arka yüzey ışınımını hesaba katan tek diyotlu eşdeğer devre kullanılmıştır. PV çıkışı, bir DC–DC yükseltici konvertör aracılığıyla EV şarj istasyonuna sağlanırken, BESS ise DC-bara voltaj regülasyonunu desteklemektedir. PI kontrolcü kazançlarını ve filtre parametrelerini optimize etmek için Sekreter Kuşu Optimizasyon Algoritması (SBOA) kullanılmıştır. Sistem bileşenlerinin modellenmesi ve önerilen metodoloji sonraki alt bölümlerde açıklanmıştır.
Bileşenlerin Modellenmesi
1. Bifacial PV sistemi
PV modülleri, gerekli akım ve voltaj çıkışının elde edilmesine yardımcı olur. Mevcut çalışmada, farklı yüzey albedosu koşulları altında bifasiyal dizilerin enerji üretim karakteristiklerini taklit etmek için iki adet monofasiyal PV modülü kullanılmıştır9. Işınım bifasiyal kazancı BGg, Denklem (1) ile belirlenir
(1)
Burada, ön ışınım Gf ve arka ışınım Gr olarak tanımlanmıştır. Maksimum çift yüzeyli PV akımı Imp ve voltaj Vmp'dir. Son olarak, seçilen seri modül sayısı Ns ve paralel modül sayısı Np ile çift yüzeyli sistem tarafından üretilen çıkış, Denklem (2) ile verilmiştir. PV kontrolü Şekil 2'de gösterilmektedir.
(2)
2. Batarya depolama sistemi (BSS)
Piller, elektrik enerjisinin depolanması için kritik bir araçtır18. Bu çalışmada enerji depolama sistemi olarak Li-ion piller kullanılmıştır. Mevcut şarj edilebilir pil seçenekleri arasında bu çalışma, verimlilikleri ve güç yoğunlukları nedeniyle elektrikli araçlar için optimal kabul edilen Li-ion pillere odaklanmaktadır. Denklem (3) ve (4), kısıtlamalarla birlikte Li-ion pili tanımlamaktadır.
(3)
PV, depolama bataryası ve dönüştürücünün özellikleri Tablo 2'de özetlenmiş; PV dizisi, enerji depolama ünitesi ve bağlı yükler arasındaki güç akışı tahsisi ise Tablo 3'te sunulmuştur.
(4)
3. DC –AC dönüştürücü
Nötr Noktası Kenetli (NPC) evirici, yenilenebilir enerji için modern güç dönüştürme sistemlerinde yaygın olarak benimsenmiştir. Geleneksel iki seviyeli yapı ile karşılaştırıldığında, NPC konfigürasyonu azaltılmış toplam harmonik distorsiyon (THD), daha düşük elektromanyetik girişim ve yarı iletken cihazlar üzerinde daha düşük voltaj stresi sunar. Bu topolojinin DC bağlantısı, nötr eklemini oluşturan eşleşmiş kapasitörler ve kenetleme diyotları aracılığıyla +Vdc/2, 0 ve -Vdc/2 olmak üzere üç potansiyel seviyeye bölünmüştür. Anahtarlama uygun şekilde kontrol edilerek, evirici çıkış terminallerinde üç ayrı voltaj durumu oluşturur ve böylece dalga formu kalitesini ve güç dönüştürme verimliliğini artırır. Bir NPC eviricideki DC veri yolunun nötr noktası, kapasitör voltajlarını dengelemede ve eviricinin düzgün çalışmasını sağlamada kilit bir rol oynar. Eviricinin kontrol şeması Şekil 3'te gösterilmiştir.
4. DC-DC Boost Dönüştürücü
DC-DC dönüştürücüler, verimliliği korurken yüksek frekanslı bir anahtarlama cihazının görev döngüsünü (duty cycle) kontrol ederek giriş DC gerilimini daha yüksek bir çıkış gerilimine yükseltirler. İki çalışma modu vardır. AÇIK (ON) durumu sırasında, diyot ters kutupludur ve giriş kaynağı enerjiyi bir manyetik alan olarak depolayan indüktöre sağlarken, çıkış kapasitörü yüke enerji sağlar. KAPALI (OFF) durumu sırasında, diyot iletim durumuna geçer. Burada, indüktörde depolanan enerji diyot üzerinden çıkışa bağlı kapasitöre ve yüke aktarılır; bu enerji giriş kaynağıyla birleşerek giriş geriliminden daha büyük bir çıkış üretir. PWM görev döngüsü değiştirilerek, dönüştürücünün çıkış gerilimi
bağıntısına göre değişir; burada D görev döngüsünü temsil eder. Bu sistemler güneş enerjisi, yakıt hücresi, EV ve depolama sistemlerinde kullanılır.
MPPT için YSA kontrol şeması
Bu çalışmada, güneş sisteminden üretilen maksimum gücü toplamak için YSA tabanlı bir MPPT seçilmiştir17. YSA eğitimi için kullanılan ışınım ve sıcaklık veri setleri; 800–1000 W/m2 güneş ışınımı seviyeleri ve 20–25 °C sıcaklık aralıklarına karşılık gelen çeşitli çevresel çalışma koşulları altında, bir görev döngüsü (D) ile MATLAB/Simulink'te oluşturulmuştur. Bifacial PV sisteminden maksimum gücü çıkarmak için YSA, elde edilen ve gereken çıktılar arasındaki MSE'yi (Op,
) minimize edecek şekilde eğitilmiştir. Burada, güneş ışınımı ve sıcaklık YSA'ya girdi olarak kullanılırken, boost konvertörü kontrol etmek için görev döngüsü çıktı olarak kullanılmıştır. Çalışmada, Şekil 4'te verildiği gibi MPPT için geliştirilen YSA modelinin yapısı ele alınmıştır. Ek Tablo 1, önerilen yöntemin diğer standart yöntemlerle karşılaştırıldığında avantajlarını ve dezavantajlarını sunmaktadır.
Bir YS Ağ (ANN), üç ana katmandan oluşur: bir çıktı katmanı (OPL), bir giriş katmanı (IPL) ve bir gizli katman (HIL). Bu yapı, verilerin IPL ile HIL arasında aktarılmasına yardımcı olur. Ardından, IPL ve HIL'yi birbirine bağlayan bağlantılardaki ağırlıklar tarafından işlem üretilir. Bu bağlamda, HL değişkenine belirli bir sapma (bias) uygulanarak hesaplamalar yürütülür ve elde edilen sonuçlar OPL değişkeninde biriktirilir. Burada, LMBP tipi YS Ağ [17] seçilmiştir. Hedef çıktıya ulaşmak için hata ölçülerek, eğitim boyunca bağlantı ağırlıkları optimize edilir. Burada, performans fonksiyonu olarak MSE kullanılarak YS Ağların eğitilmesi için LMBP kullanılır. LMBP algoritması, etkili öğrenme ve hızlandırılmış yakınsama özellikleri sergileyen ağırlıkları güncellemek için türetilmiş türevlerden yararlanır
Çok katmanlı bir algılayıcı ağdaki her nöronun bir toplama ve bir aktivasyon fonksiyonu vardır. Bununla birlikte, bu nöronları katmanlar arasında birbirine bağlayan belirli sayısal ağırlıklar (wpk) bulunur. Girişler, Denklem (5)'te belirtildiği gibi ağırlıklarla çarpıldığında; Denklem (6)'da verilen MSE minimizasyonu ile birlikte doğrusal olmayan sigmoid fonksiyonu değerlendirilir.
(5)
(6)
SBOA, önerilen sistemin filtre ve kontrol parametrelerini optimize etmiştir
Günümüzde meta-sezgisel algoritmalar, mühendislik problemlerinin çözümünde kilit bir rol oynamaktadır. Şekil 5, algoritmaların bir sınıflandırmasını sunmaktadır. SBOA algoritması23, sekreter kuşunun (SB) doğal habitatında hayatta kalma davranışını taklit eder. Sekreter kuşunun avlanma davranışının keşif ve sömürü aşamaları SBOA'da modellenmiştir. Algoritmanın keşif aşaması SB'nin yılan yakalama davranışını, sömürü aşaması ise kartal gibi yırtıcılardan kaçma davranışını taklit eder. Çözüm süreci, popülasyondaki her bir SB için kendi sınırları dahilinde rastgele değerler üretilmesiyle ve her SB'nin değerini bir tasarım parametresi olarak kabul ederek Simulink modelinin çalıştırılmasıyla amaç fonksiyonu değerinin hesaplanmasıyla başlar. Amaç değerine dayanarak, her SB daha iyi bir çözümü temsil etmek için konumunu değiştirir. SBOA'nın matematiksel modellemesi aşağıda tartışılmaktadır:
İlk hazırlık aşaması
Öncelikle, hedef fonksiyonun F(Y) tipik minimizasyonu için aramayı başlatmak amacıyla kullanılan başlangıç çözümlerinin belirlenmesi gerekir. Bu durumda, SB'lerden oluşan başlangıç rastgele popülasyonu Y = [Y1, Y2, ......YN], N adet başlangıç çözümü aracılığıyla oluşturulur. Denklem (7) popülasyonu başlatır ve Yi bir çözümü temsil eder.
(7)
Burada seçim değişkenlerinin alt ve üst sınırları lb ve ub ile gösterilmiştir. [0, 1] aralığında r rastgele bir sayıdır. N problemin boyutudur. Ayrıca, çözümün uygunluk değeri Fi = F(Yi), Yi'nin kalitesini belirlemek için kullanılır.
SB'lerin avlanma stratejisi
Avı bulma, yeme ve saldırı, SB avcılığının üç ana aşamasıdır. Tüm avlanma süreci üç eşit zaman aralığına bölünmüştür: t < 1/3T, 1/3T < t < 2/3T ve 2/3T < t < T. Burada t mevcut iterasyonu, T ise maksimum iterasyon sayısını temsil eder. Bu aralıklar, SB'nin avlanma sürecinin üç aşamasına karşılık gelir: avı arama, tüketme ve avl sayfası. Bu bölümler, kuşun avlanma aşamalarının biyolojik istatistiklerine ve her bir aşamanın sürelerine dayanmaktadır. Buna bağlı olarak, her bir SBOA aşaması şu şekilde modellenmiştir:
Av arama
SB'nin bu aşamada avını güvenli bir mesafeden takip etmesi gerekir. Tüm arama bölgesi boyunca yeterli veri toplamak için, optimizasyon algoritmalarındaki ilk adım daha güçlü bir keşif gerektirir. Diğer iki SB'nin konumlarını rehber olarak kullanan SB, ek olası bölgeleri araştırabilir. Algoritma çeşitliliğini daha da artırmak için diferansiyel mutasyon süreçleri tanıtılmıştır. t, T'nin 1/3'ünden küçük olduğunda, her bir bireyin konumunu güncellemek için Denklem (8) kullanılır. SB avı arayacak, tüketecek ve saldıracaktır. Sonuç olarak, her bir SBOA aşaması şu şekilde modellenmiştir:
(8)
Burada, Yi(t) i.ncil SB'nin mevcut konumunu, Ynewi(t) ise güncellenmiş çözümü temsil eder ve mevcut popülasyondan Yi(t) ile Yr2(t) şeklinde iki birey rastgele seçilmiştir. Parametre seçim stratejisi için kullanılan rastgele bir vektör olan R1, sömürü (exploitation) ve keşif (exploration) aşamaları arasında etkili bir denge sağlamak amacıyla [0, 1] aralığından rastgele seçilen 1 × N boyutundadır. Güncellenmiş çözüm, yalnızca daha önce elde edilen çözüme kıyasla iyileştirilmiş bir amaç fonksiyonu değeri sağladığı takdirde dikkate alınır.
Avın tüketilmesi
SB, olası avı tespit ettikten sonra çevik bir hareketlilik sergileyerek yılanın etrafında süzülür. Avın sabrı, rakipler etrafında dönerken onları gözlemleyip cezbederek savunmalarını düşürmek için kullanılır. Diğer SB'ler, avın yakınında duruşlarını ayarlar. Bu yaklaşımla avlanma başarısı önemli ölçüde artacaktır. Denklem (9), 1/3 T < t < 2/3 T olduğunda SB'nin rastgele hareketini simüle etmek için Brownian hareketinin (B) nasıl kullanıldığını göstermektedir.
(9)
Buradaki B, MATLAB'da rastgele oluşturulan r(1,M) vektörü kullanılarak standart normal dağılım ile uygulanır. SB'ler daha sonra konumlarını güncellemek için Denklem (10)'u kullanırlar.
(10)
Burada, Ybest(t) en iyi çözümü temsil eder ve üstel ölçekleme faktörü exp((t/T)4), her iterasyon için işletim kapasitesini artırmak amacıyla seçilmiştir.
Avına saldırma
Sürekli tüketimin ardından, kurbanın enerjisi tükenecektir. Sekreter kuşlar saldırıyı şimdi başlatmalıdır. Burada, Lévy uçuş tekniği kullanılarak kısa bir zaman dilimi boyunca aralıksız adımlar/sporadik uzun sıçramalar gibi farklı saldırı yöntemleri simüle edilir. Karakteristikler Denklem (11) ve (12) ile tanımlanır. SB'ler ava hızla yaklaştığı için aday seçenekler şu an optimal cevaba yakındır. Bu plan, t > 2/3T olduğunda uygulanacaktır.
(11)
(12)
Burada, LF Levy uçuş stratejisini ve
, s= 0.01 ve φ = 1.5 sabit değerlerini temsil eder. μ ve υ, [0, 1] aralığında yer alan rastgele değerlerdir. Burada, τ gammayı ve η = 0.5'i ifade eder; adım boyutunu düzenlemek ve geniş keşif sıçramaları ile yerel iyileştirme arasında bir denge sağlamak için orijinal SBOA formülasyonunda önerildiği üzere 0.5 ağırlık faktörü benimsenmiştir. Tablo 4 , geliştirilen sistemdeki ANN ve SB algoritması için kullanılan değerleri sunmaktadır.
SB için kaçış stratejisi
Yaban hayatında diğer avların peşindeyken, SBs'ler avlanma riskiyle karşı karşıyadırlar. Kartallar, şahinler, tilkiler ve çakallar, mücadele etmeleri gereken temel rakiplerdir. Tehlikeyi sezdiklerinde, kendilerini veya avlarını korumak için çeşitli kaçınma teknikleri kullanmak zorundadırlar. Bu algoritma, koşu modlarını (D2) ve kamuflajı (D1) modelleyerek kaçış yöntemlerini simüle eder. Rakiplerle karşılaştıklarında, sekreter kuşları güvende kalmak için başlangıçta çevrelerine uyum sağlarlar. Burada sekreter kuşları, konumlarını avın (en iyi birey) etrafında güncelleyerek algoritmalardaki yerel optimumlardan kaçınmaya çalışırlar. Eğer düşmandan kaçamazlarsa, güvende kalmak için uçuş veya hızlı koşu taktiklerini kullanacaklardır. Referans olarak, bir yerel optimumla sınırlı kalmamak için rastgele bir birey Yrand lider olarak seçilir. Sekreter kuşları konumlarını güncellemek için Denklem (13) ve (14)'ü kullanırlar. SBOA akış şeması Şekil 6'da gösterilmiştir.
(13)
(14)
Tasarım değişkenlerinin temsili
Bu çalışmadaki problem değişkenleri; DC-DC dönüştürücülerin PI kontrolör kazançları gibi tasarım parametrelerinin yanı sıra filtrelerin ve dönüştürücülerin direnç ve endüktanslarını içermektedir. Denklem (15), kontrol değişkenlerinin temsilini vermektedir. Kontrol parametresi sınırları Ek Tablo 2'de listelenmiştir.
(15)
Uygunluk fonksiyonu (FF)
THD'nin minimize edilmesi, bu çalışmada amaç fonksiyonu (Obj) olarak seçilmiştir. FF'nin maksimum değeri, Denklem (16) ve (17)'de verilen Obj'nin minimum değeridir.
(16)
(17)