Dans cette recherche, un système autonome de recharge de véhicules électriques (VE) alimenté par un système photovoltaïque bifacial et un stockage d'énergie par batterie (BESS) est développé, comme illustré dans la Figure 1, en utilisant une méthode MPPT commandée par un réseau de neurones artificiels (ANN). Un circuit équivalent à une seule diode est utilisé pour modéliser le champ photovoltaïque bifacial, en tenant compte de la température, de la résistance série, de la résistance shunt, de l'irradiance solaire et de l'irradiance côté arrière. La sortie du champ photovoltaïque est acheminée vers la station de recharge de VE via un convertisseur élévateur continu-continu (DC–DC), tandis que le BESS assure la régulation de la tension du bus continu. L'algorithme d'optimisation inspiré du comportement du secrétaire (Secretary Bird Optimization Algorithm, SBOA) est utilisé pour optimiser les gains du régulateur PI ainsi que les paramètres du filtre. La modélisation des composants du système et la méthodologie proposée sont décrites dans les sous-sections suivantes.
Modélisation des composants
1. Système photovoltaïque bifacial
Les modules photovoltaïques permettent d'atteindre le courant et la tension de sortie requis. Dans la présente étude, deux modules photovoltaïques monofaciaux ont été utilisés pour simuler les caractéristiques de production d'énergie de panneaux bifaciaux dans différentes conditions d'albédo de surface9. Le gain bifacial du rayonnement BGg est déterminé par l'équation (1)
(1)
Ici, l'irradiation frontale est Gf, et l'éclairement arrière est Gr. Le courant maximal photovoltaïque bifacial est Imp, et la tension est Vmp. Enfin, la production de sortie du système bifacial avec le nombre choisi de modules en série Ns et en parallèle Np est donnée par l'équation (2). La commande du photovoltaïque est indiquée dans la Figure 2
(2)
2. Système de stockage par batterie (SSB)
Les batteries constituent un moyen essentiel de stockage de l'énergie électrique18. Cette étude utilise des batteries au lithium-ion comme système de stockage d'énergie. Parmi les options de batteries rechargeables disponibles, ce travail se concentre sur les batteries au lithium-ion, considérées comme optimales pour les véhicules électriques en raison de leur efficacité et de leur densité énergétique. Les équations (3) et (4) décrivent la batterie au lithium-ion avec des contraintes.
(3)
Les caractéristiques du panneau photovoltaïque, de la batterie de stockage et du convertisseur sont résumées dans le Tableau 2, tandis que le Tableau 3 présente la répartition du flux de puissance entre le champ photovoltaïque, l'unité de stockage d'énergie et les charges connectées.
(4)
3. Convertisseur CC – CA
L'onduleur à point neutre calé (NPC) est largement utilisé dans les systèmes modernes de conversion d'énergie pour les énergies renouvelables. Par rapport à la structure classique à deux niveaux, la configuration NPC offre une distorsion harmonique totale (THD) réduite, une interférence électromagnétique plus faible et une contrainte de tension moindre sur les dispositifs semi-conducteurs. Le bus continu de cette topologie est divisé en trois niveaux de potentiel +Vdc/2, 0 et -Vdc/2, grâce à des condensateurs appariés et des diodes de calage formant le point neutre. En contrôlant de manière appropriée la commutation, l'onduleur génère trois états de tension discrets aux bornes de sortie, améliorant ainsi la qualité de l'onde et l'efficacité de la conversion d'énergie. Le point neutre du bus continu dans un onduleur NPC joue un rôle clé dans l'équilibrage des tensions des condensateurs et le bon fonctionnement de l'onduleur. Le schéma de commande de l'onduleur est illustré dans Figure 3.
4. Convertisseur élévateur CC-CC
Les convertisseurs DC-DC élèvent la tension d'entrée continue à une tension de sortie plus élevée en contrôlant le cycle de service d'un dispositif de commutation à haute fréquence tout en maintenant l'efficacité. Il existe deux modes de fonctionnement. Pendant l'état ON, la diode est polarisée en inverse, et la source d'entrée fournit de l'énergie à l'inductance, qui la stocke sous forme de champ magnétique, tandis que le condensateur de sortie alimente la charge. Pendant l'état OFF, la diode devient polarisée en direct. L'énergie stockée dans l'inductance est alors libérée à travers la diode vers le condensateur connecté en sortie et vers la charge, se combinant avec la source d'entrée pour produire une tension de sortie supérieure à la tension d'entrée. En faisant varier le cycle de service du PWM, la tension de sortie du convertisseur varie selon la relation
, où D représente le cycle de service. Ce dispositif est utilisé dans les systèmes solaires, les piles à combustible, les véhicules électriques et les systèmes de stockage.
Schéma de commande ANN pour le MPPT
Dans ce travail, un MPPT basé sur un réseau de neurones artificiels (ANN) est sélectionné afin de collecter la puissance maximale générée par le système solaire17. Les ensembles de données d'irradiation et de température utilisés pour l'apprentissage du réseau ANN ont été générés sous MATLAB/Simulink dans des conditions environnementales variables, correspondant à des niveaux d'irradiance solaire de 800 à 1000 W/m2 et à des plages de température de 20 à 25 °C, avec un rapport cyclique (D). Le réseau ANN est entraîné pour minimiser l'erreur quadratique moyenne (MSE) entre les sorties obtenues et celles requises (Op,
) afin d'extraire la puissance maximale du système photovoltaïque bifacial. Ici, l'irradiation solaire et la température sont utilisées comme entrées du réseau ANN, tandis que le rapport cyclique constitue la sortie permettant de commander le convertisseur élévateur. La structure du modèle ANN développé pour le MPPT, présentée dans la Figure 4, est prise en compte dans cette étude. Le Tableau supplémentaire 1 présente les avantages et inconvénients de la méthode proposée par rapport à d'autres méthodes classiques.
Un réseau de neurones artificiels (ANN) se compose de trois couches principales : une couche de sortie (OPL), une couche d'entrée (IPL) et une couche cachée (HIL). Cela permet de transférer les données entre l'IPL et l'HIL. Ensuite, ces données sont produites par les poids des liens reliant l'IPL et l'HIL. Dans ce contexte, les calculs sont effectués avec un biais déterminé appliqué à la variable HIL, et les résultats obtenus sont accumulés dans la variable OPL. Ici, on choisit un ANN de type LMBP [17]. Les poids des connexions sont ajustés au cours de l'apprentissage en mesurant l'erreur afin d'atteindre la sortie cible. Le LMBP est utilisé pour entraîner les ANN avec l'erreur quadratique moyenne (MSE) comme fonction de performance. L'algorithme LMBP utilise les dérivées obtenues pour mettre à jour les poids, ce qui confère des propriétés d'apprentissage efficace et d'une convergence accélérée.
Chaque neurone dans un réseau de perceptron multicouche possède une fonction de sommation et une fonction d'activation. Cependant, il existe certains poids numériques (wpk) qui relient ces neurones entre les niveaux. Lorsque les entrées sont multipliées par les poids conformément à l'équation (5). La fonction sigmoïde non linéaire est considérée avec la minimisation de l'erreur quadratique moyenne donnée dans l'équation (6).
(5)
(6)
SBOA a optimisé les paramètres du filtre et de la commande du système proposé
De nos jours, les algorithmes métaheuristiques jouent un rôle clé dans la résolution de problèmes d'ingénierie. Figure 5 présente une classification des algorithmes. L'algorithme SBOA23 imite le comportement du SB pour survivre dans son habitat naturel. Les phases d'exploration et d'exploitation du comportement de chasse du serpentaire de Beaudouin sont modélisées dans l'algorithme SBOA. La phase d'exploration de l'algorithme imite le comportement du SB lorsqu'il attrape des serpents, tandis que sa phase d'exploitation reproduit son comportement d'esquive face à des prédateurs tels que les aigles. Le processus de résolution commence par la génération de valeurs aléatoires comprises dans leurs limites respectives pour chaque SB de la population, puis par le calcul de la valeur de la fonction objectif en exécutant le modèle Simulink, chaque valeur de SB étant traitée comme un paramètre de conception. En fonction de la valeur objectif, chaque SB modifie sa position afin de représenter une solution améliorée. La modélisation mathématique de l'algorithme SBOA est décrite ci-dessous :
Phase de préparation initiale
Tout d'abord, il est nécessaire d'identifier les solutions initiales qui ont été utilisées pour commencer la recherche d'une minimisation typique de la fonction objectif F(Y). Dans ce cas, la population aléatoire initiale Y = [Y1, Y2, ......YN] de SB est formée à partir de N solutions initiales. L'équation (7) initialise la population et Yi représente une solution.
(7)
Les bornes inférieure et supérieure des variables de choix sont notées lb et ub. Dans [0, 1], r est un nombre aléatoire. N est la dimension du problème. En outre, la valeur d'ajustement de la solution Fi = F(Yi) est utilisée pour quantifier la qualité de Yi.
Stratégie de chasse des SB
Trouver, manger et attaquer la proie sont les trois phases principales de la chasse du SB. L'ensemble du processus de chasse a été divisé en trois intervalles de temps égaux : t < 1/3T, 1/3T < t < 2/3T et 2/3T < t < T. Ici, t représente l'itération courante et T le nombre maximal d'itérations. Ces intervalles correspondent aux trois phases de la prédation du SB : recherche, consommation et attaque de la proie. Ces divisions reposent sur les données biologiques relatives aux phases de chasse de l'oiseau et sur la durée de chaque phase. Par conséquent, chaque phase du SBOA est modélisée comme suit :
À la recherche de proies
Le SB doit chasser ses proies à distance sécurisée durant cette phase. Afin de recueillir des données suffisantes sur toute la zone de recherche, la première étape des algorithmes d'optimisation exige une exploration plus poussée. En utilisant les positions des deux autres SB comme guide, le SB peut examiner d'autres régions possibles. Afin d'augmenter davantage la diversité de l'algorithme, des processus de mutation différentielle sont donc introduits. Lorsque t est inférieur au tiers de T, l'équation (8) est utilisée pour mettre à jour la position de chaque individu. Il cherchera, consommera et attaquera les proies. Par conséquent, chaque phase de SBOA est modélisée comme suit :
(8)
Où Yi(t) est la position actuelle du ième SB, Ynewi(t) donne la solution mise à jour, deux individus Yi(t) et Yr2(t) ont été choisis aléatoirement parmi la population existante. Un vecteur aléatoire, pour la stratégie de sélection des paramètres, R1 de dimension 1 × N est choisi aléatoirement dans l'intervalle [0, 1] afin d'assurer un équilibre efficace entre les phases d'exploitation et d'exploration. La solution mise à jour n'est prise en compte que si elle conduit à une amélioration de la valeur de l'objectif par rapport à la solution précédemment obtenue.
Consommation de proies
Le SB utilise une agilité vive et se déplace en planant autour du serpent après avoir repéré une proie potentielle. La patience de la proie sera exploitée pour abaisser sa garde, en l'observant et en l'attirant tandis qu'ils tournent l'un autour de l'autre. D'autres SB ajustent leur posture près de la proie. Ce comportement augmente considérablement les chances de succès de la chasse. L'équation (9) illustre comment le mouvement brownien (B) est utilisé pour simuler le déplacement aléatoire du SB lorsque 1/3 T < t < 2/3 T.
(9)
Le B est implémenté dans MATLAB à l'aide de la distribution normale standard par le vecteur généré aléatoirement r(1,M). Les SB utilisent ensuite l'équation (10) pour mettre à jour leurs positions.
(10)
Où Ybest(t) est la meilleure solution et le facteur d'échelle exponentiel exp((t/T)4) est choisi pour augmenter la capacité d'exploitation à chaque itération.
Attaque des proies
Après une consommation continue, la victime épuisera son énergie. Les oiseaux-serpents devraient lancer l'assaut maintenant. Ici, différentes méthodes d'attaque, telles que des pas incessants ou des bonds longs et sporadiques sur une courte période, sont simulées à l'aide de la technique du vol de Lévy. Les caractéristiques sont décrites par les équations (11) et (12). Les options candidates se situent actuellement près de la réponse optimale, car les oiseaux-serpents s'approchent rapidement de leur proie. Ce plan sera mis en œuvre lorsque t > 2/3T.
(11)
(12)
Où LF représente la stratégie de marche de Lévy et
, avec s = 0,01 et φ = 1,5, des valeurs fixes. μ et υ sont des valeurs arbitraires comprises dans [0, 1]. Ici, τ désigne la fonction gamma et η = 0,5 ; le facteur de pondération 0,5 a été adopté, comme proposé dans la formulation initiale de SBOA, afin de réguler la taille du pas et de maintenir un équilibre entre de grands sauts exploratoires et un affinement local. Tableau 4 présente les valeurs utilisées pour le RNA et l'algorithme SB dans le système développé.
Stratégie d'échappement pour SB
Lorsqu'ils chassent d'autres proies dans la nature, les SB risquent d'être chassés eux-mêmes. Les aigles, les faucons, les renards et les chacals sont leurs principaux adversaires auxquels ils doivent faire face. Ils doivent utiliser une variété de techniques d'évasion pour se protéger, eux ou leurs proies, lorsqu'ils perçoivent un danger. Cet algorithme simule les méthodes de fuite en modélisant les modes de déplacement (D2) et le camouflage (D1). Face à des adversaires, les oiseaux-serpents s'intègrent d'abord à leur environnement pour rester en sécurité. Ici, les oiseaux-serpents s'efforcent d'éviter les optima locaux dans les algorithmes en mettant à jour leurs positions autour de la proie (l'individu le meilleur). Ils utiliseront des tactiques de vol ou de course rapide pour rester en sécurité s'ils ne parviennent pas à échapper à l'ennemi. À titre de référence, un individu aléatoire Yrand est choisi comme chef afin d'éviter d'être limité à un optimum local. Les équations (13) et (14) sont utilisées par les oiseaux-serpents pour mettre à jour leurs positions. Le diagramme de flux de l'algorithme SBOA est présenté dans Figure 6.
(13)
(14)
Représentation des variables de conception
Les variables du problème dans cette étude comprennent les paramètres de conception, tels que les gains du régulateur PI des convertisseurs DC-DC, ainsi que les résistances et inductances des filtres et des convertisseurs. L'équation (15) donne la représentation des variables de commande. Les bornes des paramètres de commande sont indiquées dans le Tableau supplémentaire 2.
(15)
Fonction de fitness (FF)
La minimisation de THD est choisie comme fonction objective (Obj) dans cette étude. La valeur maximale de FF correspond à la valeur minimale de Obj, donnée par les équations (16) et (17).
(16)
(17)