48

Système inertiel de stockage d’énergie associé à des générateurs éoliens

I.Hamzaoui1, F.Bouchafaa1, A.Talha1
1
Laboratory of Instrumentation, Faculty of Electronics and Computer, University of Sciences and Technology
Houari Boumediene, BP 32 El-Alia 16111 Bab-Ezzouar Algiers, Algeria.
Email: [email protected], [email protected], [email protected].

Résumé- Le type d’unité de production décentralisée qui fait l’objet de cet article concerne les sources à énergie renouvelable et plus
particulièrement l’énergie éolienne. L’objectif étant d’une part de proposer une solution facilitant l’intégration des éoliennes
existantes et d’autre part de proposer un système de génération d’énergie à base d’éolienne à vitesse variable pouvant participer à la
gestion du réseau. En raison du caractère très fluctuant et aléatoire du vent, une éolienne seule ne peut pas faire le réglage de la
production. Afin d’assurer l’équilibre production –consommation, on propose une solution constituée un système de stockage
de type inertiel qui jouant le rôle d’un tampon entre la source et le consommateur. Cette solution possède une bonne dynamique, un
bon rendement et une durée de vie similaire à l’éolienne. Les résultats de simulation de l’association d’un SISE avec un générateur
éolienne, pour augmenter son taux de pénétration dans le réseau et améliore son contribution au service système sous
Matlab/Simulink son présentés.

Mots clés : Générateurs Eoliens, Système Inertiel de Stockage d’Energie (SISE) – Volant inertiel - MADA - Commande Vectorielle

1. Introduction

1
types de stockage [1][2].
Les générateurs éoliens sont des générateurs dont la Le système étudier est connecté au réseau
source primaire d’énergie est le vent. Il est bien connu '1 (Fig.1). Comporte:
que le vent a des caractéristiques très fluctuantes et - Un aérogénérateur à vitesse variable utilisant une
imprévisibles et il est impossible de prédire sa valeur pour MADA et pilotée par le rotor au moyen des
un moment donné. convertisseurs MLI à deux niveaux, un multiplicateur de
Comme la puissance mécanique débitée par la turbine est vitesse et une turbine à trois pales.
très fluctuante, la puissance générer par le générateur
EO
- Un SISE Comportant un volant d’inertie, une machine
éolien sera aussi fluctuante. Cela pose beaucoup de
asynchrone à cage (MAS) et des convertisseurs
problèmes pour les gestionnaires du système énergétique
(AC/DC/AC) MLI à deux niveaux.
pour deux raisons. Il faut assurer la balance entre la
puissance générer et la puissance consommée. De plus, la En cet article, nous commençons par une brève
puissance consommée est difficilement prévisible et présentation du modèle d’éolien, la machine asynchrone
variable. A cause de ces restrictions, les générateurs à double alimentation et sa commande indépendante des
éoliens actuels ne peuvent pas fonctionner sans être puissances actives et réactive. Notre étude est concentrée
associés à une source classique d’énergie [2][4][5]. sur l’association d’un système inertiel de stockage
IC

En cet article l’idée de stockage est exploitée, afin de d’énergie au générateur éolienne afin de constituer un
régler la puissance débitée sur le réseau. Le stockage ensemble capable de fournir des services système ou de
inertiel est une solution appropriée aux aérogénérateurs, fonctionner en ilote.
où il offre de meilleurs avantages par rapport aux autres

Préseau , Qréseau
Ps , Qs
C arbre

Ω t Ω méc

Pr , Qr

PStockage

Fig.1. Schéma global du système étudié

 49

2. Modèle de la turbine
0.5

La puissance aérodynamique (extraite du vent) apparait

Coefficient de puissance (Cp)

0.4
au niveau du rotor de la turbine dépend du coefficient de
puissance Cp. Elle est donner par: 0.3

1
= C p ( λ , β ) ρπ R 2 v 3
0.2

Paero (1)
2 0.1

Ou ρ la densité de l’air, R le rayon de la pale et v la 0

vitesse du vent. 0 2 4 6 8 10
Rapport de vitesse (Landa)
12 14

Cp est le coefficient de puissance, dont la limite théorique

Fig.2. Coefficient de puissance Cp en fonction de rapport de vitesse λ
(limite de Betz) est de 0.59, dépend de plusieurs
paramètres [4],[5],[8]: L’équation (6) donne l’expression de la puissance
• Le rapport de vitesse λ défini par: maximale obtenue en utilisant la stratégie MPPT
RΩ t (Maximum Power Point Tracking) qui permet de régler
λ = (2)
automatiquement la vitesse spécifique à sa valeur
v
• L’angle de calage β. optimale (λopt), de manière à obtenir le coefficient de
puissance maximale (Cpmax):
Pour la turbine utilisée, le coefficient de puissance est
1 ρπ R t C p max
5

approché par la formule suivante [8]: PMPPT = Ω 3t (6)

2 λ 3opt
C2 C5
C P ( λ , β ) = C1 ( − C 3 − C 4 ). exp( ) + C6λ (3)
λi λi La maximisation de la puissance et représenter sous forme
de schéma bloc simplifier donner dans la figure 3.

1
Avec 1 = 1 0 . 035 Plusieurs correcteurs peuvent être utilise, nous avons
−
λi λ + 0 . 08 . β β 3
+ 1 choisi un régulateur de type PI pour effectuer notre
modèle.
C 1 = 0 .5176 , C 2 = 116 , C 3 = 0 .4 , C 4 = 5, C 5 = 21, C 6 = 0 .0068 .
'1 3. Modèle de la MADA
Le couple aérodynamique (éolien) est le rapport de la
puissance à la vitesse mécanique de la turbine Ωt : Les équations électriques classiques du MADA dans le
modèle de PARK sont écrites comme suit:
Paero 1
EO
C aer = = C p ( λ , β ). ρπ R 2 v 3 (4)  d
Ωt 2Ω t v ds = R s .i ds + dt Φ ds − ωs Φ qs
La vitesse mécaniques de la turbine et celle de la  d
v qs = R s .i qs + Φ qs + ωs Φ ds (7)
génératrice sont liées par le gain du multiplicateur, ainsi  dt
que pour le couple aérodynamique et celle issu du  d
v
 dr = R .i
r dr + Φ dr − ω r Φ qr
multiplicateur sont exprimés par la relation suivante:  dt
 v qr = R r .i qr + d Φ qr + ω r Φ dr
C g  1  C aero  (5)  dt
Ω  = G Ω   M
Ce = −P L (Φ qs .i dr − Q ds .i qr )
IC

 t  mec 
 s
La figure 2 représente la variation de Cp en fonction du
rapport de vitesse λ pour un angle de calage nul. Les puissances actives et réactives statorique sont définies
Il est clair de cette figure qu'il y a une valeur de λ Pour par:
lequel le coefficient de puissance est maximisé,
 Ps = v ds ids + v qs iqs (8)
maximisant ainsi la puissance pour une vitesse de vent 
donnée. Q s = v qs ids − v ds iqs

λ R .Ω
Ωturbine Ωmec
β v
t 1
G

Cp
1 1 1 Ω mec
v Cp ρ Sv 3 1
2 Ωt G Cg Jp + f
C aero

Cem − ref
Ωturbine−ref Ω ref
λCp max .v 1
R G

Fig.3. Schéma bloc de la MPPT extraite avec asservissement de la vitesse

 50

Les puissances actives et réactives statorique sont définies type inertiel est envisage, un volant d’inertie couplé
par: mécaniquement à une machine asynchrone et pilotée par
un convertisseur de puissance comme représenté à la
 Pr = vdr idr + vqr iqr (9) figure 1.
 L’énergie Ev stockée dans le volant d’inertie Jv à pour
Qr = vqr idr − vdr iqr expression:
1
=
2
Le couple électromagnétique est exprimé par: E V J vΩ v (14)
2
Cem = p (φds .iqs − φqs ids ) (10)
Pour calculer l’inertie du Volant, on se base sur une
puissance à fournir durant un temps ∆t: on souhaite que le
Où p et le nombre de paire de pôle stockage inertiel fournisse la puissance nominale PmasN
durant un temps ∆t l’énergie nécessaire vaut alors:
4. Commande de puissance active et réactive
∆Ev= PmasN ∆t. Sachant que : ∆Ev=1 /2 Jv.∆Ωv2
Un contrôle vectoriel de cette machine (MADA) a été et que ∆Ω v2 = ∆Ω vMAX
2
− ∆Ω vMIN
2
il vient:
conçu en orientant le repère de PARK pour que le flux
statorique suivant l’axe q soit constamment nul 2P ∆t
Jv = masN (15)
φ ds = φ s et φ qs = 0 [1][5][6][7]. (Ω 2
vMAX − Ω 2vMIN )
Alors le couple et simplifier: Nous constatons d’après cette équation que lors de la
M (11) fabrication d’un volant d’inertie pour un SISE, il ya deux
C em = −p φ s .i qr condition initiales dont il faut tenir compte : La vitesse
Ls
maximale de rotation du volant et la capacité du volant
Le couple électromagnétique et la puissance actives d’inertie. La machine asynchrone est choisie en fonction

1
dépendra seulement du courant rotororique d'axe q. de ces avantages en termes de simplicité et de robustesse
En négligeant la résistance du stator Rs, on aura: des parties tournantes, son modèle dans le repère de
vds = 0, vqs = vs . Ainsi, les courants rotoriques sont liés '1 PARK, avec orientation du flux rotorique (Φqr=0 et
Φdr=Φr), peut être décrit par les équations suivantes:
aux puissances active et réactive statoriques comme suit:
 dΦ dr Rr MR r
 M  dt = − L Φ dr + (ω s − pωω) qr + L i ds
Ps = −vs L iqr (12)  r r
 s  dΦ qr = − R r Φ − (ω − pωω) + MR r i

EO
v Φ vM  dt Lr
qr s dr
Lr
qs
Q = s s
− s i dr  di (16)
 s Ls  ds MR r Mp Ω R sr 1
Ls  = 2
Φ dr + Φ qr − i ds + ω s i qs + v ds
 dt σL s L r σL s L r σL s σL s

Une régulation sur les courants rotorique permet de  di qs MR r Mp Ω R sr 1

 dt = σL L2 Φ qr − σL L Φ dr − σL i qs − ω s i ds + σL v qs
contrôler l'échange de puissance entre le stator et le  s r s r s s

 C pM
em − ref =
réseau. Φ dr i qs
 Lr
L'expression des tensions rotorique en fonction des
courants rotorique permet de compléter la description du En distingue deux zones de fonctionnement pour la
machine électrique présenté à la figure 5. [1], [5]
IC

système à réguler dont le schéma fonctionnel est présenté

sur la figure 4.
 M2 didr M2
vdr = R r i dr + (Lr − ) − gωs (Lr − )iqr ;
(13)
 Ls dt Ls
PmasN (w)
 2 2
v = R i + (L − M ) diqr + gω (L − M )i + gω Mvs . puissance
 qr r qr r
Ls dt
s r
Ls
dr s
ωs Ls CemN (Nm)

Couple
Où g est le glissement du générateur.

0 Ω vN 2Ω vN Ω v ( rad / s

Fig.5. Allure de la puissance et du couple en fonction de la vitesse

- Pour 0 ≤ Ω v ≤ Ω vN , le couple nominale de la machine

est disponible, mais la puissance maximale est variable,
en fonction de la vitesse ( Pmas = k Ω v ), et plus petite que
la puissance nominale. Cette zone ne présente pas
beaucoup d’intérêt pour les SISE.
Fig.4. Schéma bloc de la Commande en puissance de la MADA - Pour Ω v > Ω vN , la puissance est maximum et
correspond à la puissance nominale de la machine, le
5. Système inertiel de stockage d’énergie couple électromagnétique est alors inversement
Dans l’objectif de faire participer une éolienne à vitesse proportionnel à la vitesse de rotation ( C em = k / Ω v ).
variable aux services systèmes, le stockage d’énergie de C’est la zone de fonctionnement utilisée dans les SISE
 51

parce qu’ici la puissance de la machine est disponible 2 X 1500 X 2.5

J= ≈ 0.1Kg.m 2
pour n’importe quelle vitesse. (314) 2 − (157) 2
Le stockage inertiel avec machine asynchrone sera donc
utilisé dans la plage de vitesse suivant Les figures 7 et 8 illustrent le fonctionnement du système
Ω vN ≤ Ω v ≤ 2 Ω vN . Permettant ainsi un fonctionnement de stockage inertiel. La puissance de référence de la
figure 7 est inversée à chaque fois que la vitesse du volant
à puissance nominale constante. A partir d’une puissance
atteint une limite haute ou basse en 2,5s. Ainsi, on
de référence Pv-ref, on peut déduire le couple
demande à la MAS de fournir où de stocker l’énergie
électromagnétique de référence de la machine, Cem-ref,
nominale. Une puissance positive correspond à une
entrainant le volant d’inertie par une mesure de la vitesse
puissance consommée par la machine et une puissance
de rotation, Ωv-mes.
négative correspond à une puissance fournie par la
Pv − ref machine.
C em − ref = (17)
Ω v − mes 2000

1500
Le couple électromagnétique de référence, devra être
1000
limité au couple nominale pour la plage de vitesse
500
compris entre 0 et la vitesse nominale, au delà de la

Pv (Kw)
0
vitesse nominale, la valeur du couple diminuera de
-500
manière à garder le produit Cem-ref.Ωv constant. La
-1000
réduction du couple est effectuée par le défluxage de la
-1500
machine au delà de la vitesse de synchronisme. La loi de
-2000
défluxage est comme suit: 0 1 2 3 4 5
Temps (s)
6 7 8 9 10

Fig.7. Puissance absorbée et débitée par la MAS

Pv−ref L r 1

1
320

Φ dr−ref = (18) Vitesse de rotation du volant (rad/s) 300

pMiqs Ω v −mes
280

L’estimation du flux est donnée par l’équation suivante:

Φ dr − estimé =
1+
M
L r
S
i ds − mes
'1
(19)
260

240

220

200
R r
EO
180

La figure 6 présente le schéma fonctionnel de la 160

commande du SISE. Les courants ids-ref et iqs-ref sont 140

0 1 2 3 4 5 6 7 8 9 10
déterminés par le régulateur de flux pour le courant d’axe Temps (s)
d, et par la référence du couple électromagnétique pour le
courant d’axe q. Le couple électromagnétique étant Fig.8. Vitesse du volant (rad/s)
calculé à partir de l’équation (17), le courant en La fonction principale du SISE est de lisser la puissance
quadrature est déterminé par inversion l’équation du fournie par le générateur éolien qui peut causer plusieurs
couple (16). problèmes dans le réseau. Pour réduire au minimum les
IC

fluctuations de cette puissance, le SISE doit assurer la

6. Résultat de simulation compensation des variations de la puissance de vent.
La valeur du coefficient d’inertie a été calculée pour une La puissance de référence du SISE est déterminée par la
plage de vitesse comprise entre 1500 et 3000 tr/min, et différence entre la puissance générée par le générateur
une puissance nominale de 1.5kW durant un temps éolienne et la puissance qu’il faut délivrer au réseau ou
correspondant à 2.5s: sur les charges isolées.

Ω v− mes

Φ dr − estimé i ds − mes

U C − mes
Φ i ds − ref
dr − ref
vds vqs

Pv − ref iqs − ref i qs − ref

i qs − mes

Fig .6. Schéma bloc de la commande du système de Stockage inertiel

 52

Pv − ref = Preseau − Péol . Selon le principe illustré par la correspondant au mode de fonctionnement hypo
synchrone et augmente sur les périodes où la MADA
figure 9.
fonctionne en régime hyper synchrone, et sa valeur est
presque constante dans des plages où la vitesse mécanique
de la MADA est proche de 157rd/s.
La figure 15 nous montre le transit de puissance entre le
SISE et le réseau. Le sens d’écoulement de cette
puissance varie en fonction de la puissance générée par
l’aérogénérateur. Si cette dernier est inférieur à -7.5Kw, le
transfert d’énergie du réseau vers le SISE (stockage), du
SISE vers le réseau (déstockage) afin de compenser le
déficit de puissance de l’éolien.
Les composantes du flux rotorique direct et quadrature de
la MAS sont illustrées sur la figure 16. La composante en
Fig.9. La puissance lisser par SISE quadrature est toujours nulle, ce qui justifie la commande
à flux rotorique orienté.
Dans notre application, la puissance de référence du
7. Conclusion
réseau est fixée à la valeur de -7.5kW, et pour maintenir
le facteur de puissance du réseau unitaire la puissance Dans cet article, nous avons étudié un générateur éolien à
réactive de référence du réseau et fixée à une valeur nulle. vitesse variable associé à un système de stockage inertiel.
La figure 10 présente les différentes grandeurs turbine- Le SISE est contrôlé par une puissance de référence
MADA. Comme la fréquence des grandeurs statorique de obtenue en fonction de la puissance générer et de la
la MADA et imposée par le réseau 50Hz, le glissement puissance à envoyer au réseau, dans le but était de lisser

1
varie sur une plage de 12% en passant par des valeurs la puissance injectée dans le réseau afin de participer aux
nulles. Il est à noter que le profil du vent est choisi de services auxiliaires (réglage de la tension, la fréquence,
façons à obtenir justement ces différentes valeur du '1 compensation …). Les résultats de simulation montrent
glissement, c'est-à-dire différents modes de clairement le stockage et la restitution d’énergie.
fonctionnement (hypersynchrone, hyposynchrone,
8. Références
synchronisée) de la MADA. La figure 10-c montre
l’allure de la puissance mécanique sur l’arbre (générer [1] G.O.Cimuca, "Système inertiel de stockage d’énergie associé à des
éoliennes", thèse de doctorat, Ecole Nationale Supérieure d’Arts et
par la turbine) qui varie entre -7kW et -8kW, sur les
EO
Métiers, Centre de Lille, 2005.
phases de fonctionnement en mode hyposynchrone entre [2] M.khaterchi, J.elhadj, M.elleuch, "DPC for Three-Phase Inverter to
-7 kW et -7.5kW, et sur les périodes de fonctionnement en Improve the Integration of Wind Turbine Associated to Flywheel
mode hyper synchrone entre -7.5 kW et -8kW. La figure Energy Storage System into the Grid" the 6th international Muli-
Conference on systems, Signals and Devices, (SSD) 2009.
10-d présente les puissances active et réactive échangées [3] A.Boyette, "Contrôle-Commande d’un générateur asynchrone à
entre le réseau et le rotor. Cette relation s’établit double alimentation avec système de stockage pour la production
directement du fait que le convertisseur (AC/DC/AC) ne éolienne ", Thèse de doctorat de l’université Henri Poincaré, Nancy
contient aucun élément de stockage intermédiaire, leur 2006.
[4] A.Davigny, "Participation aux services système de fermes
sens d’écoulement dépend du signe du glissement. Pour
IC

d’éoliennes à vitesse variable intégrant du stockage inertiel d’énergie",

glissement négative, les puissances sont transmises du thèse de doctorat, Ecole Doctorale Sciences Pour l’ingénieur, Lille,
rotor au réseau (mode génératrice hypersynchrone). Pour 2007.
glissement positive, les puissances sont transmises du [5] L.Leclerq, "Apport du stockage inertiel associé a des éoliennes dans
un réseau électrique en vue d'assurer des services systèmes", Ecole
réseau au rotor (mode génératrice hyposynchrone). Pour
Doctorale Sciences Pour l’ingénieur, Lille, 2004.
un glissement nul, le rotor absorbe une puissance active [6] K.Ghedamsi "Contribution à la modélisation et la commande d’un
qui représente les pertes joule rotorique. convertisseur direct de fréquence Application à la conduit de la machine
La figure 11 présente les puissances actives et réactives asynchrone ", Thèse de doctorat, Ecole nationale polytechnique 2008.
[7] S.EL-Aimani. "Modélisation de différentes technologies d’éoliennes
du réseau. La figure 12 montre que la tension est imposée intégrées dans un réseau de Moyenne Tension", Thèse de Doctorat
par le réseau et l’amplitude de courant injecté au réseau d’Etat en Electronique et Génie Electrique. Ecole Centrale de Lille
est imposée par la commande. La tension et le courant du (ECL), 2005.
réseau sont en opposition de phase pour tous les modes [8] N. Laverdure, "Sur l’intégration des générateurs éoliens dans les
réseaux faibles ou insulaires" , Thèse de Doctorat, Institut National
de fonctionnement, ceci justifier que l’écoulement de
Polytechnique de Grenoble, décembre 2005.
puissance est toujours de l’aérogénérateur vers le réseau. [9] M.Aktarujjaman, M.A.Kashem, M.Negnevitsky, G.Ledwich
La figure 13 présente la tension et le courant de la MAS, "Smoothing output power of a doubly fed wind turbine with an Energy
et sur une période correspondant au mode de Storage System", Conférence 2006
[10] A.ARBAOUI, "Aide à la décision pour la définition d’un système
fonctionnement hyposynchrone de la MADA, la MAS
éolien adéquation au site et a un réseau faible", Thèse de Doctorat, École
fonctionne en génératrice (phase de déstockage), le Nationale Supérieure d’Arts et Métiers Centre de Bordeaux,2006
courant et en avance de phase par rapport à la tension [11] G.Cimuca, M. Radulescu,C.Saudemont,B.Robyns,"Losses and
fig.13-a, et sur une période correspondant au mode de Efficiency of a Flywheel Energy Storage System with Permanent-
Magnet Synchronous Machine Associated to a Variable-Speed Wind
fonctionnement hyper synchrone de la MADA, la MAS Generator"
fonctionne en moteur (phase de stockage) le courant et en
retard par rapport a la tension fig.13-b.
La figure 14 montre la variation de vitesse du volant
d’inertie. Cette dernière diminue sur des plages
53

200
0.8
180
0.7

e(rad/s) 160 0.6

0.5
u

140

n
t
iq

e
0.4

m
n

(b)
ca

e
120
0.3

s
ssemé

lis
(a)

G
100 0.2

0.1
ite

80
V

0
60
-0.1

40 -0.2
0 2 4 6 8 10 12 14 16 18 20 0 2 4 6 8 10 12 14 16 18 20
Temps (s) Temps (s)

-2000
2000

-3000
1500

-4000 1000
Fig .16. Composante direct et quadrature du
w )

(C)
e(K

500

W
-5000

)
in

r(K
rb

0
Ptu

-6000
(d)

P
-500
-7000
-1000

-8000
-1500

-9000 -2000
0 2 4 6 8 10 12 14 16 18 20 2 4 6 8 10 12 14 16 18 20
Temps (s) Temps (s)

FigFig.9. Lesdifférentes
.10. Les différentesgrandeurs
grandeursturbine-MADA
turbine-MADA

2000
200
400

0 150 300
Pres (Kw), Qres(KVAR)

Vs(V) et Is X 2.5 (A)

200
100
-2000 Qres
100
50
Is (A)

-4000 0

0
Pres -100
-6000
-50 -200

1
-8000 -300
-100

-400
-10000 -150
2 4 6 8 10 12 14 16 18 20 0 2 4 6 8 10 12 14 16 18 20 4.12 4.125 4.13 4.135 4.14 4.145 4.15 4.155 4.16 4.165 4.17
Temps (s) Temps (s) Temps (s)

'1
EO

Fig .11. Les puissances actives et réactives du réseau Fig .12. Tension et courant du réseau
IC

300 250 250

350
200
(a) 2000 Vmas
200 (b) Vmas
200 1000
150 150
Vmas (V) ,Imas x 20 (A)

Vmas (V) ,Imas x 20 (A)

Vmas (V) ,Imas x 20 (A)

300
100 0 100
100
50 -1000 50
250
m(rad/s)

0 0 -2000 0
Pv-ref

-50 -50
200 -3000
-100
-100 -100
W

-4000
-150 -150
-200150
-5000
-200
Imas Vmas -200

-6000
Imas Imas
-300100 -250 -250
0 1 2 3 4 5 6 7 8 9 10 8.1 8.15 8.2 8.25 8.3 2.2 2.25 2.3 2.35 2.4 2.45
Temps (s) -7000
Temps (s) Temps (s)
50
0 2 4 6 8 10 12 14
Fig.13.
16
Les grandeurs
18 20
électriques
-8000
0 2
statorique
4 6
de la8 MAS
10 12 14 16 18 20
Temps (s) Temps( s)

Fig .14. Vitesse du volant d’inertie Fig .15. La puissance de stockage et sa référence

1.8 0.3

1.6
0.2

1.4
0.1
1.2
Qdr (wb)

Qqr (wb)

1 0

0.8 -0.1

0.6
-0.2
0.4

-0.3
0.2

0 -0.4
0 2 4 6 8 10 12 14 16 18 20 0 2 4 6 8 10 12 14 16 18 20
Temps(s) Temps (s)