See discussions, stats, and author profiles for this publication at: [Link]

net/publication/368637070

Etude et simulation des composants de système hybride éolien-

photovoltaïques avec stockage (supercondensateur _ batteries )

Article in IEE Proceedings A Science Measurement and Technology · February 2023

CITATIONS READS

0 1,929

4 authors, including:

Amina Benhamou El Madjid Berkouk

University M'Hamed Bougara of Boumerdes National Polytechnic School of Algiers
36 PUBLICATIONS 396 CITATIONS 316 PUBLICATIONS 3,724 CITATIONS

SEE PROFILE SEE PROFILE

All content following this page was uploaded by Amina Benhamou on 28 February 2023.

The user has requested enhancement of the downloaded file.

ème Conférence Internationale des Energies Renouvelables (CIER-2016)
Proceedings of Engineering and Technology – PET
Vol.14, pp.167 - 173

Etude et simulation des composants de système

hybride éolien-photovoltaïques avec stockage
(supercondensateur _ batteries )
ELHAMRANI Ramadan 1, BENHAMOU Amina*2,BERKOUK Madjid 3
1
Laboratoire életronique de puissance-Ecole militaire polytechnique
Bordj-el-Bahri-Alger-Algérie
1
ramadan_elhamrani@[Link]
2*
Boumerdes University-Boumerdes-Algeria
2
amina_benhamou@[Link]
3
école nationale polytechnique-elHarrach –Alger-Algeria

Résumé: l’objectif de ce travail est Les systèmes de production hybrides, travaillant en parallèle avec le réseau électrique,
d’énergie hybrides PV-éolien, se présentent comme l’une des options appelés aussi connectés réseau. Ces systèmes contribuent à
technologiques pouvant permettre d’opérer de façon progressive la satisfaire la charge du système électrique du pays. Les systèmes
transition entre les sources renouvelables (solaire et éolien). En hybrides du deuxième groupe fonctionnent en régime isolé ou
particulier pour les zones isolées d’Afrique Subsaharienne qui
en mode autonome. Ils doivent répondre aux besoins des
disposent d’un potentiel important en ressource solaire et pour
lesquelles la production d’énergie est généralement assurée par les consommateurs situés dans des sites éloignés du réseau
groupes Diesels . électrique.
Pour permettre d’exploiter au mieux ces sources d’énergie, le système Le système d’énergie hybride peut fonctionner avec la
doit être correctement dimensionné pour qu’il puisse fournir présence ou non d’une source d’énergie conventionnelle qui
l’énergie d’une manière continue. Le dimensionnement d'un tel peut être un générateur diesel, une micro turbine à gaz…etc,
système implique la détermination des caractéristiques nominales des aussi à la présence ou non d’un dispositif de stockage qui
sources et du système de stockage, pour des ressources énergétiques permet d’assurer une meilleure satisfaction des charges
et des contraintes d’utilisation données. électriques pendant les périodes d'absence d'une ressource
Pour ce fait, une modélisation et simulation des différents éléments
primaire à convertir en électricité. La gestion de la capacité de
de système hybride solaire et éolien , sont illustrés dans ce travail,
après une recherche parmi les différents concepts et techniques de stockage limitée, ainsi que leur coût élevé représente un
ces composants. inconvénient, et afin de minimiser le recours au stockage, on a
Mots clés— énergie solaire, énergie éolienne, système hybride, eu recours aux systèmes d’énergie hybride. [3]
stockage supercondensateur, hacheur buck-boost-simulation. La structure du système hybride peut contenir un système
photovoltaïque, une éolienne, un générateur hydraulique, ou
I. INTRODUCTION bien une combinaison de plusieurs sources renouvelables. Le
La croissance économique est devenue équivalente à une choix des sources se base sur une analyse technico- économique.
consommation énergétique de plus en plus importante. Les Dans notre travail, on s’est intéressé à faire le bon choix
systèmes de génération éolienne et photovoltaïque sont les plus technique des composants de notre système hybride éolien
prometteurs parmi les technologies de production d’énergie photovoltaïque, avec stockage et isolé du réseau. (fig 1).
renouvelables, par leurs nombreux avantages. Cependant,
chacune des technologies mentionnées ci-dessus a ses propres
inconvénients [1]. Néanmoins, le système d’énergie hybride est
utilisé pour surmonter l’intermittence, l’incertitude, et la faible
disponibilité de chaque source d'énergie renouvelable, ce qui
rend le système plus fiable.
Les systèmes d'énergie hybride (éolien –photovoltaïque) bien
adaptés à une production connectée au réseau ou bien
décentralisée peuvent contribuer à résoudre plusieurs problèmes
d’approvisionnement énergétique [2].
Le système d’énergie hybride a connu des nouveaux
développements ces dernières années pour ce qui concerne les
méthodes de gestions de la production, ainsi que l'optimisation
du dimensionnement par des logiciels spécialisés (HOMER, Fig. 1 schéma de système hybride éolien –photovoltaïque
SOMES, RAPSIM...).

II. CLASSIFICATION DES SYSTEMES D’ENERGIE HYBRIDE A. Le générateur photovoltaique:

Les systèmes d’énergie hybride sont classifiés suivant La caractéristique (I/V) d’un GPV est non linéaire (fg.2). Le
plusieurs critères. module PV produit une puissance à un point appelé point de
Il existe deux principaux groupes des systèmes d’énergie fonctionnement qui appartient à la courbe caractéristique. Les
hybride. Dans le premier groupe, on trouve les systèmes coordonnées de ce point sont la tension et le courant de
Copyright IPCO-2017
ISSN 2356-5608
ème Conférence Internationale des Energies Renouvelables (CIER-2016)
Proceedings of Engineering and Technology – PET
Vol.14, pp.167 - 173
fonctionnement. De même que pour la cellule, le MPP
correspond au point de fonctionnement pour lequel le GPV.
Figure 5 : Modèle d’une cellule réelle.

( Sont donnés par le fabriquant.

a) Courant de saturation ( )
Pour avoir l’expression du courant de saturation de
référence (Io,n), on remplace par les conditions du point du
circuit ouvert (Voc,n, 0) dans l’équation suivante :

[Link]éristique I/V de module PV

: Courant de saturation aux conditions standards ;

: Tension du circuit ouvert aux conditions standards ;

: Potentiel thermique aux conditions standards.

Dans les conditions de fonctionnement, le courant de

saturation est donné par :
[Link]éristique de PV au MPPT.

La figure 3 présente un exemple de la caractéristique (P/V) et

(I/V) d’un module PV. Elle illustre l’importance de faire : Énergie de gap du semi-conducteur ( =1.12
fonctionner le système au MPP pour tirer le maximum de la eV pour le silicium)
puissance disponible [4].
Une cellule photovoltaïque est modélisée de manière
simple par une source idéale de courant qui produit un courant
Iph proportionnel à la puissance lumineuse (figure 4), en
parallèle avec une diode qui correspond à la jonction p-n.

Figure [Link] de l’intensité en fonction de la tension et la puissance en

Vcel fonction de l’intensité
V
D
Pour étudier l’influence de l’éclairement (G) et la température
(Tc) sur les caractéristiques du panneau on fait varier l’un et on
fixe l’autre à sa valeur de référence, les résultats de simulation
sont représentés dans les figures 7 et 8

Figure 4 : Modèle d’une cellule idéale

Le circuit représenté ci-dessus est régi par la relation suivante:

Avec : Figure [Link]éristique de la

Ce modèle reste théorique et ne prend pas compte du puissance de panneau PV et le courant en fonction de l’intensité.
comportement d’une cellule photovoltaïque dans des conditions D’après les résultats de simulation, nous remarquons que
réelles, où on observe une perte de tension en sortie ainsi que le courant de court-circuit est fortement lié à l’irradiation solaire,
des courants de fuites. On modélise donc cette perte de tension il est directement proportionnel avec ce facteur. Contrairement à
par une résistance en série Rs, et les courants de fuite par une Icc, la tension du circuit ouvert Voc est légèrement affectée par la
résistance en parallèle Rp (figure 5). variation de l’irradiation, elle augmente un peu avec
l’augmentation de cette dernière, [5]on peut conclure donc que
la puissance maximale du panneau est fortement liée à
Iph Id Ip I Rs l’irradiation solaire.
VD Rp V
Copyright IPCO-2017
ISSN 2356-5608
 ème Conférence Internationale des Energies Renouvelables (CIER-2016)
Proceedings of Engineering and Technology – PET
Vol.14, pp.167 - 173
Pour la température nous constatons que le courant de Lors du fonctionnement de la MADA en moteur, la puissance
court-circuit est peu sensible aux variations de cette dernière. Ps est fournie par le réseau [8].
Nous remarquons aussi que la tension à vide est inversement Dans le modèle mécanique proposé, l’inertie totale J est
proportionnelle avec la température ce qui engendre une petite celle de la génératrice Jg et l’inertie de la turbine Jt ramenée au
diminution de la puissance. rotor de la génératrice.

J J g JtG2
Simulation du générateur photovoltaïque
L’évolution de la vitesse mécanique dépend du couple
Dans notre travail, nous avons utilisé un générateur mécanique appliqué au rotor de la génératrice Tmec qui est la
photovoltaïque constitué de 12 panneaux en série et 5 panneaux résultante d’un couple électromagnétique produit par la
en parallèle afin de produire la puissance dimensionnée du génératrice Tem, d’un couple de frottement visqueux Tvis et d’un
générateur qui est de 3.5 KWatt. Les résultats de simulation sont couple du multiplicateur de vitesse Tg.
présentés dans la figure 9 :

Tmec J d: dt
4000 20
3500
3000 15
Tmec Tg Tem Tvis
Pgpv (W)

2500
Igpv (A)

2000 10
1500
1000
Le couple issu du frottement est modélisé par le
5
500 coefficient de
0
0 50 100 150 200 250 300 0 t frottement visqueux f
0 50 100 150 200 250 300
Vgpv (V)
Vgpv (V) Tvis f :
nous aurons :
Figure 9. Caractéristiques P-V et I-V du générateur photovoltaïque étudié. Tg Tem J d: f : dt
B. Eolienne à vitesse variable à base de la Machine A partir des équations précédentes la REM et le schéma
Asynchrone à Double Alimentation (MADA) : bloc de l’arbre mécanique sont représentés comme suit :
La Machine Asynchrone Doublement Alimentée (MADA) a
suscité un intérêt particulier surtout en tant que génératrice
Te
m
dans le domaine de l’énergie éolienne. Pour les éoliennes
Tg Tme 1
c
utilisant la MADA, le stator de celle-ci est directement T
vis
J
couplé au réseau alors que son rotor est connecté au réseau sf
à travers une interface composée de deux convertisseurs
(a) REM (b) Schéma bloc
statiques (convertisseur coté MADA et convertisseur coté
Fig. [Link] et schéma bloc du modèle de l’arbre mécanique
réseau) [6]
La REM et le schéma bloc représentant le modèle de
l’ensemble de la chaîne cinématique de la turbine (Fig. 12), sont
obtenus en regroupant les composantes respectives développées
précédemment (Fig. 9, 10 et 11)

Figure 10. schéma de l’éolienne MADA

Elle possède un stator identique à celui d’une machine
asynchrone classique ou d’une machine synchrone. La
différence réside dans le rotor composé d’enroulements Fig 12. la chaine cinématique de la turbine.
triphasés disposés de la même manière que les enroulements
Cette figure montre que la turbine produit un couple
statoriques et connectés en étoile [7], ce qui permet d’accéder
aérodynamique Taero à partir de la puissance cinétique Pv du vent
aux grandeurs rotoriques.
et par action de l’angle d’orientation des pales.
La vitesse mécanique de l’arbre de la génératrice : résulte
de l’application du couple moteur disponible à la sortie du
multiplicateur Tg auquel s’oppose le couple électromagnétique
Tem. Le contrôle de cette dernière peut donc être effectué soit par
action sur l’angle d’orientation des pales soit à travers le couple
électromagnétique de la génératrice. [9]

Figure [Link] du rotor de la MADA Simulation du modèle de la turbine: Dans la simulation

numérique de la partie mécanique du système éolien on a utilisé
Copyright IPCO-2017
ISSN 2356-5608
ème Conférence Internationale des Energies Renouvelables (CIER-2016)
Proceedings of Engineering and Technology – PET
Vol.14, pp.167 - 173
les modèles mathématiques. L’ensemble des modèles de la
turbine, multiplicateur de vitesse et arbre mécanique ont été
simulés en conservant l’architecture réelle. La vitesse du vent v,
l’angle d’orientation des pales et le couple électromagnétique
Tem fournis par la MADA constituent les entrées du système
tandis que le couple de sortie du multiplicateur Tg et la vitesse
mécanique ! sont les sorties.

L’angle d’orientation des pales est fixé à une valeur

constante (égale à 2). Le couple électromagnétique de la MADA Figure 16 (c) coefficient de puissance (d) rapport de vitesse de vent.
correspond à un point de fonctionnement optimal quant à la
puissance générée (Maximum Power Point Tracking MPPT).

Figure 17 (e) couple de multiplicateur (f) couple électromagnétique de MADA.

Fig.13. Profil du vent appliqué à la turbine pendant 250 secondes. C. Système de stockage :
La figure 13, représente le profil du vent, figure 14 (a) et Pour un bon fonctionnement du système d’énergie hybride, le
(b) montrent les vitesses de la turbine et la vitesse mécanique de système de stockage joue un rôle crucial, il permet d’avoir une
l’arbre, on remarque que pendant 250s, le système éolien passe continuité de service et une meilleure qualité d’énergie fournie.
par deux de fonctionnement hypo et hyper synchrone. Nous rappelons quelques paramètres électriques employés pour
caractériser une batterie, ces derniers sont:
 La capacité nominale (Q) : C’est le nombre maximal
d'ampères-heures (Ah) qui peut être extrait de la batterie,
pour des conditions de décharge données.
 L'état de charge « SOC » (State of charge) : C’est le
rapport entre la capacité présente « q »et la capacité
nominale « Q » : (0 ≤ SOC ≤1). Si SOC = 1 la
batterie est totalement chargée et si SOC = 0 la batterie est
Figure 14.(a) vitesse de la turbine , (b) vitesse de l’arbre totalement déchargée.
Les figures 15 (a) et (b) présentent respectivement la variation  Le régime de charge (ou décharge) : C’est le paramètre
du rapport de vitesse et le coefficient de puissance Cp en qui reflète le rapport entre la capacité nominale d'une
fonction du temps. On vérifie, que les valeurs de batterie et le courant auquel celle-ci est chargée (ou
coefficient de puissance n’atteignent pas la valeur théorique déchargée). Il est exprimé en heures.
maximale donnée par Betz (0,59). La puissance aérodynamique  La durée de vie : C’est le nombre de cycles
disponible au niveau de la turbine Paero et le couple charge/décharge que la batterie peut soutenir avant de
aérodynamique Taero sont illustrés respectivement dans les perdre 20% de sa capacité nominale. [9]
figures 16.c et 16.d. Les figures 16.e et 16.f montrent
respectivement le couple issu du multiplicateur et le couple
appliqué par la MADA qui agissent directement sur l’arbre.

Figure 18. Évolution de la tension Vbatt, durant la charge et la décharge.

La courbe de décharge (Q-V) est donnée par la figure 19.

Figure 15 (a) coefficient de puissance (b) rapport de vitesse de vent.

Copyright IPCO-2017
ISSN 2356-5608
ème Conférence Internationale des Energies Renouvelables (CIER-2016)
Proceedings of Engineering and Technology – PET
Vol.14, pp.167 - 173

Figure 19. Courbe de décharge (Q-V) de la batterie [11].

En analysant la figure ci-dessus, on peut constater la Figure 21. Hacheur réversible en courant associé aux batteries
présence de trois points spécifiques sur la caractéristique (Q- Les deux interrupteurs « K11 » et « K10 » forment une
V) : ces trois points sont : la tension de la pleine charge (E0), la cellule de commutation. Leur fonctionnement sera donc
tension correspondante à la fin de la zone exponentielle (E exp) et complémentaire [12]. Si on associe deux fonctions de connexion
la tension correspondante à la fin de la zone nominale (E Nom).
et représentant respectivement l’état de conduction des
Les équations de charge et de décharge sont données interrupteurs « K11 » et « K10 », on obtient :
comme suit :

Selon l’état des deux interrupteurs « K11 » et « K10 »,

deux configurations sont possibles (figure 22):

La figure 20 , montre la caractéristique de décharge de la

batterie utilisée et l’évolution de la tension de cette dernière
pour différents courants de décharge.
Figure 22. Les deux configurations du convertisseur buck-boost
courant de décharge nominal (30A) réversible en courant.
14

13
Configuration 1 : K11 fermé et K10 ouvert
Voltage (V)

12

11

10

9 Configuration 2 : K11 ouvert et K10 fermé

0 1 2 3 4 5 6 7
t (heures)
E0 = 12.269, R = 0.0008, K = 0.0022475, A = 0.82078, B = 6

14
5A
13 10 A En analysant ces deux configurations, on peut conclure
Voltage (V)

15 A que les relations entre les grandeurs d’entrée (Vbat,Ibat) et celles

12
de sortie (Udc, Ibat,m) du convertisseur sont données par le
11
système d’équations ci-dessous
10

0 5 10 15 20 25 30 35 40 45
t (heures)

Figure 20. Caractéristique de décharge et tension de la batterie pour différents

courants de décharge.

Nous remarquons que plus le courant de décharge est Simulation du système de stockage
important, plus la batterie se décharge plus rapidement.
L’insertion du stockage dans le système hybride a été
[Link]élisation du hacheur buck-boost réversible en courant simulée dans les deux modes de fonctionnement (isolé et
Le hacheur buck-boost est utilisé pour le système de parallèle) et sous différentes contraintes telles que les variations
stockage par batteries. Il est considéré comme un étage de la charge et de l’éclairement. Dans cette section, le
d’adaptation, et il permet de réaliser les opérations de charge et générateur photovoltaïque est commandé de sorte à produire le
de décharge vu sa réversibilité en courant. Ainsi, il est élévateur maximum de puissance, la tension du bus continu est régulée à
une valeur de référence de 700V et la capacité du bus continu
de tension pour la décharge, et abaisseur de tension pour la
charge. [10] (Cdc) est supposée initialement chargée à 700V.

Le modèle du hacheur buck-boost associé aux batteries est Pour tester les performances du système de stockage,
représenté sur la figure 21. nous avons procédé à la variation de l’éclairement afin de varier
Ipv, et au changement du mode de fonctionnement (isolé à

Copyright IPCO-2017
ISSN 2356-5608
 ème Conférence Internationale des Energies Renouvelables (CIER-2016)
Proceedings of Engineering and Technology – PET
Vol.14, pp.167 - 173
parallèle) pour montrer le comportement du système de En se référant à la figure 24, nous constatons que le courant des
stockage pendant les deux phases [11] (charge et de décharge). batteries suit bien le courant de référence avec une très bonne
dynamique.
Les résultats de simulations sont présentés sur les figures ci-
dessous : À partir des figures (24, 25, 26), nous remarquons ce qui suit :
Paramètres de simulation : Pch=1500W. Cdc=4.7 mF. Du 0s à 0.2s : le mode de fonctionnement du système est le
PGPV=1000W. Pstock=2400W. mode isolé, le système de stockage se retrouve dans la phase de
706
décharge ( ) et la puissance demandée par la charge
Udc
Udc ref
est assurée par les batteries de stockage et le générateur PV, en
704
effet :
Udc (V)

702

700

698
Conclusion: le but de ce travail est l’intégration des énergies
696
0 0.2 0.4 0.6 0.8 1
renouvelables dans la production d’électricité mondiale, la
t (s) modélisation de la chaîne photovoltaïque avec tous ces
Figure 23. Régulation de la tension du bus continu composants, est effectuée en appliquant les principaux
algorithmes de maximisation de puissance MPPT, en particulier
la méthode de perturbation et d’observation (P&O).
10
Ibatt
Ibatt ref
5
Nous avons développé le modèle mathématique de la Machine
Ibatt (A)

0 Asynchrone à Double Alimentation MADA avec l’étude de la

-5
partie mécanique de la chaîne éolienne.
Nous avons intégré un système de stockage hybride batteries–
-10
0 0.1 0.2 0.3 0.4 0.5
t (s)
0.6 0.7 0.8 0.9 1 super condensateurs, dont la modélisation et l’asservissement
était détaillée avec son hacheur buck-boost Une description
Figure 24. Régulation du courant des batteries détaillée de la commande fera l’objet d’autres travaux,
1.5
Références :
1

[1]B. Multon, " Production d’énergie électrique par sources

Ipv (A)

0.5

renouvelables ", Techniques de l’Ingénieur, traité Génie

0
électrique D 4 005.
-0.5
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
[2] M. Amin, "Efficiency and Power Density Improvement of
t (s)
Grid-Connected Hybrid Renewable Energy Systems utilizing
Figure 25. Évolution du courant Ipv
High Frequency-Based Power Converte", Doctorat, Université
Internationale de Florida, 2012.
[3] I. Vichiu, ''Modélisation et analyse de l'intégration des
4
énergies renouvelable dans un réseau autonome'', Mémoire de
3 Doctorat, Université du Havre, Décembre, 2005.
2 [4] S. Astie,"Conversion photovoltaïque : du rayonnement
solaire à la cellule ", Techniques de l’Ingénieur, D 3 936.
Ich (A)

-1
[5]W. Omran, "Performance Analysis of Grid-Connected
-2 Photovoltaic Systems", Thèse de Doctorat, Université de
-3
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Waterloo, Ontario, Canada, 2010.
t (s)
[6]F. Poitiers, "Etude et commande de génératrices
Figure 26. Évolution du courant de charge coté bus continu
Asynchrones pour L’utilisation de l’énergie éolienne", Doctorat,
90.001 12.3 université de Nantes, 2003.
12.28
[7]B. Multon, X. Roboam, B. Dakyo, "Aérogénérateurs
90.0005
électriques", Techniques de l’Ingénieur, traité Génie électrique
Vbatt (V)
SOC (%)

12.26

90
D 3 960.
12.24 [8]J. M. Noël, "Énergie éolienne pour la fourniture d’électricité",
89.9995 12.22
Techniques de l’Ingénieur, traité Génie électrique BE 8 585.
[9] S. Elaimani, "Modélisation de différentes technologies
12.2
89.999
0 0.2 0.4 0.6 0.8 1
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 d’éoliennes intégrées dans un réseau de moyenne tension",
t (s)
t (s) Thèse de doctorat de l’Ecole Centrale de Lille, 2004. T.
Figure 27.Évolution de l’état de charge et la tension d’une batterie de [10] Scimone, " Energy Storage in Electric Power Generation
stockage. Plant from Renewable Sources ", Université Degli Studi di
La figure 27 montre l’évolution de la tension du bus continu, Catania, Italie, 2010.
nous remarquons que cette dernière est régulée à sa valeur de [11]P. F. Ribeiro, B. K. Johnson, et M. L. Crow, "Energy
référence. Storage Systems for Advanced Power Application",
Proceedings Of The IEEE, VOL. 89, NO. 12, December 2001.
[12]J. De Dinechin , "Le stockage énergétique", Cours A5-2
Energies Renouvelables. ENSTA, 2007.
Copyright IPCO-2017
ISSN 2356-5608
ème Conférence Internationale des Energies Renouvelables (CIER-2016)
Proceedings of Engineering and Technology – PET
Vol.14, pp.167 - 173

Copyright IPCO-2017
ISSN 2356-5608

View publication stats