Eoliennes à vitesse variable basées sur des machines
Eoliennes à vitesse variable basées sur des machines
Eoliennes à vitesse variable basées sur des machines
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
Etude d’éoliennes <strong>à</strong> <strong>vitesse</strong> <strong>variable</strong><br />
<strong>basées</strong> <strong>sur</strong> <strong>des</strong> <strong>machines</strong> asynchrones<br />
(MAS-MADA)<br />
Equipe : Conversion et qualité de l’énergie<br />
M. Machmoum, Machmoum,<br />
F. Poitiers, L. Moreau et M.E. Zaïm, Zaïm,<br />
R. Le Doeuff<br />
LARGE-GE44, LARGE GE44, Bd de l’Université, BP 406, 44602 Saint Nazaire Cedex<br />
email : mohamed.machmoum<br />
mohamed. machmoum@polytech polytech.univ univ-nantes. nantes.fr fr
Plan<br />
• Rappels <strong>sur</strong> la conversion de l’énergie éolienne : objectifs<br />
et génératrices utilisées<br />
• Machine asynchrone <strong>à</strong> cage : fonctionnement en<br />
autonome<br />
• Machine asynchrone <strong>à</strong> cage <strong>à</strong> <strong>vitesse</strong> <strong>variable</strong><br />
• Machine asynchrone <strong>à</strong> double alimentation connectée <strong>sur</strong><br />
le réseau.<br />
• Conclusions et Perspectives
C p<br />
P m<br />
45<br />
40<br />
35<br />
30<br />
25<br />
20<br />
15<br />
10<br />
5<br />
Rappels <strong>sur</strong> la conversion d’énergie<br />
0<br />
2 3 4 5 6 7 8 9 10 11<br />
16000<br />
14000<br />
12000<br />
10000<br />
8000<br />
6000<br />
4000<br />
2000<br />
λ<br />
8 m.s -1<br />
4 m.s -1<br />
7 m.s -1<br />
6 m.s -1<br />
5 m.s -1<br />
λ<br />
0<br />
1000 1500 2000 2500 3000 3500 4000<br />
Ω<br />
par éolienne<br />
Courbe C p(λ) (spécifique <strong>à</strong> chaque éolienne).<br />
C p : coefficient de puissance<br />
λ : Rapport d ’avance (<strong>vitesse</strong> de l’extrémité<br />
<strong>des</strong> pales <strong>sur</strong> la <strong>vitesse</strong> du vent).<br />
Dimension<br />
Densité de<br />
d ’une pale<br />
l ’air<br />
Vitesse<br />
’air<br />
Vitesse du<br />
vent<br />
1 2<br />
2 3<br />
Pm= ( ) 1<br />
2 1<br />
Puissance mécanique disponible en<br />
fonction de la <strong>vitesse</strong> du générateur<br />
pour différentes <strong>vitesse</strong>s de vent.<br />
p<br />
R<br />
C R V<br />
KV ρπ<br />
Ω
P m<br />
16000<br />
14000<br />
12000<br />
10000<br />
8000<br />
6000<br />
4000<br />
2000<br />
Objectifs de cette conversion<br />
8 m.s -1<br />
4 m.s -1<br />
7 m.s -1<br />
6 m.s -1<br />
5 m.s -1<br />
0<br />
1000 1500 2000 2500 3000 3500 4000<br />
Ω<br />
En <strong>vitesse</strong> fixe : le maximum<br />
théorique n’est pas atteint<br />
En <strong>vitesse</strong> <strong>variable</strong> : La puissance<br />
maximale est exploitée pour toutes<br />
les <strong>vitesse</strong>s du vent (régulation de<br />
λ optimal pour avoir Cp<br />
maximum)
MAS<br />
MASDA<br />
MAS<br />
R-L<br />
Génératrices utilisées<br />
RESEAU DE<br />
DISTRIBUTION<br />
C1<br />
U<br />
RESEAU DE<br />
DISTRIBUTION<br />
Machine asynchrone en autonome.<br />
* Nécessité de capacités d’autoexcitation<br />
pour magnétiser la machine.<br />
Machine asynchrone <strong>à</strong> cage débitant<br />
<strong>sur</strong> un réseau.<br />
* Obligation de fonctionner au<br />
voisinage du synchronisme.<br />
Autorise le fonctionnement <strong>à</strong><br />
<strong>vitesse</strong> <strong>variable</strong> ( MAS ou<br />
MADA).
3<br />
GS<br />
Génératrices utilisées<br />
Machine synchrone <strong>à</strong> aimants permanents<br />
(MSAP)<br />
* bon rendement, faibles puissances,<br />
adaptée aux faibles <strong>vitesse</strong> mais prix<br />
élevé.<br />
Machine <strong>à</strong> réluctance <strong>variable</strong> (MRV)<br />
permettant de supprimer totalement ou<br />
partiellement le multiplicateur de<br />
<strong>vitesse</strong>.<br />
* Etu<strong>des</strong> actuellement en cours.
Machine asynchrone <strong>à</strong> cage en<br />
fonctionnement autonome<br />
TURBINE<br />
EOLIENNE<br />
Ω m Γ<br />
R<br />
V abc<br />
MAS<br />
(modèle de<br />
PARK)<br />
MAS<br />
i abc<br />
R L CHARGE<br />
Charge +<br />
excitation<br />
V abc<br />
Modèle diphasé associé <strong>à</strong> une transformation étoile-triangle côté charge<br />
méthode simple d’étude de régimes déséquilibrés, gain en temps de calcul et<br />
modèles indépendants machine - charge/excitation.<br />
Paramètres : Machine LS, 3 kW, p=1, 220/380 V, J=0.117 kg.m 2 , f=5.10 -3 N.m.s -1
E<br />
350<br />
300<br />
250<br />
200<br />
150<br />
100<br />
50<br />
E<br />
Machine asynchrone <strong>à</strong> cage en<br />
fonctionnement autonome<br />
R r /g<br />
C=30µF<br />
C=40µF<br />
l g<br />
0<br />
0 0.5 1 1.5 2 2.5 3<br />
C=35µF<br />
I m<br />
I M<br />
M<br />
Points de fonctionnement pour différentes<br />
valeurs de capacités<br />
M<br />
Valeur de la mutuelle M (H)<br />
E<br />
0.5<br />
0.45<br />
0.4<br />
0.35<br />
0.3<br />
0.25<br />
0.2<br />
I C<br />
C<br />
R<br />
x : points expérimentaux<br />
-: Approximation polynomiale<br />
0.15<br />
0 0.5 1 1.5 2 2.5 3 3.5 4<br />
Courant magnétisant Im (A)<br />
Prise en compte de la saturation<br />
I m
Vsa (V)<br />
Vsc (V)<br />
400<br />
200<br />
0<br />
-200<br />
-400<br />
0<br />
400<br />
200<br />
0<br />
-200<br />
-400<br />
0<br />
Machine asynchrone <strong>à</strong> cage en<br />
fonctionnement autonome<br />
Voltage build-up<br />
2 s<br />
1 2 3 4 5 6 7 8<br />
time (s)<br />
Stator voltage of phase c<br />
Tensions stator <strong>des</strong> phases a et c<br />
Stator voltage of phase a<br />
2 s 2 s<br />
1 2 3 4<br />
Time (s)<br />
5 6 7 8<br />
Simulation Relevé expérimental<br />
Processus d ’amorçage en régime équilibré
Vsc (V)<br />
400<br />
300<br />
200<br />
100<br />
-100<br />
-200<br />
-300<br />
0<br />
Machine asynchrone <strong>à</strong> cage en<br />
fonctionnement autonome<br />
Voltage fall-down<br />
Stator voltage of phase c<br />
-400<br />
20 22 24 26 28 30<br />
Time (s)<br />
32 34 36 38 40<br />
Ca=Cb Ca= Cb=Cc=28.5 =Cc=28.5 µF, , Ω=3000 =3000 tr/min, essai <strong>à</strong> vide<br />
Déconnexion brusque d’un condensateur ou de la charge<br />
Ca=Cb Ca= Cb=Cc=33 =Cc=33 µF, , Ω=3000 =3000 tr/min, Ra=Rb Ra= Rb=Rc Rc=1k =1kΩ<br />
100 V/div<br />
1s/div<br />
100 V/div<br />
1s/div<br />
Simulation Relevé expérimental
350<br />
300<br />
250<br />
200<br />
150<br />
100<br />
50<br />
Machine asynchrone <strong>à</strong> cage en<br />
fonctionnement autonome<br />
0<br />
4 4.5 5 5.5 6 6.5 7<br />
t(s)<br />
Tension<br />
stator (V)<br />
260<br />
240<br />
220<br />
200<br />
180<br />
160<br />
140<br />
28µF<br />
∆V<br />
Influence d’une charge purement résistive<br />
Vitesse (rad/s)<br />
Tension (V)<br />
Fréquence (Hz)<br />
Influence d’un impact de <strong>vitesse</strong><br />
BILAN<br />
32µF<br />
∆V<br />
36µF<br />
∆V<br />
120<br />
200 400 600 800 1000 1200 1400 1600<br />
P (W)<br />
* Impossibilités d’utiliser <strong>des</strong> charges<br />
exigeantes en tension et fréquence.<br />
* Risque de démagnétisation.<br />
* Fonctionnement <strong>sur</strong> une plage de<br />
<strong>vitesse</strong> restreinte.<br />
autres solutions !
• Avantages de la <strong>vitesse</strong> <strong>variable</strong><br />
Vitesse <strong>variable</strong><br />
• Optimisation de l’énergie captée<br />
• Contrôle du transfert de puissance et énergie propre envoyée <strong>sur</strong><br />
le réseau.<br />
• Inconvénients de la <strong>vitesse</strong> <strong>variable</strong><br />
• Utilisation de <strong>machines</strong> spéciales<br />
• Convertisseur de puissance « complexe »
Génératrice asynchrone<br />
<strong>à</strong> <strong>vitesse</strong> <strong>variable</strong><br />
• nécessité d’utiliser un convertisseur de fréquence (ac ( ac/dc dc dc/ac dc ac)<br />
• ac/dc ac dc<br />
redresseur <strong>à</strong> dio<strong>des</strong> associé <strong>à</strong> un convertisseur dc/dc dc dc<br />
(<strong>machines</strong> synchrones)<br />
redresseur <strong>à</strong> MLI (<strong>machines</strong> asynchrones)<br />
• gestion du transfert de puissance entre le redresseur <strong>à</strong> MLI et<br />
l’onduleur (régulation du bus continu) et le placement au point<br />
de puissance optimum de l’éolienne.<br />
2 approches possibles : approche source de courant ou<br />
approche source de tension.
MAS associée <strong>à</strong> un redresseur MLI<br />
Régulation de Puissance (1)<br />
• Vdc régulée par l ’onduleur, le redresseur fixe le couple<br />
• Commande vectorielle de la MAS <strong>sur</strong> modèle diphasé :<br />
Cref = Pref/Ωrotor et Ψref constant<br />
Vitesse<br />
du vent<br />
Génératrice<br />
asynchrone<br />
Vitesse de rotation<br />
rotor<br />
Algorithme de<br />
placement au point<br />
optimal de<br />
fonctionnement<br />
Puissance<br />
de référence<br />
÷<br />
Flux de référence<br />
Couple<br />
de référence<br />
Redresseur MLI<br />
Commande<br />
vectorielle<br />
Vdc
MAS associée <strong>à</strong> un redresseur MLI<br />
Régulation de Puissance (2)
MAS associée <strong>à</strong> un redresseur MLI<br />
Régulation de Vdc (1)
MAS associée <strong>à</strong> un redresseur MLI<br />
Régulation de Vdc (2)<br />
Vitesse déclenchement<br />
• t=1.15s Pref=Pnom/4 et <strong>à</strong> t= 2s Pref=Pnom/2 Equilibre puissance as<strong>sur</strong>é - bon choix du<br />
rendement de l’ensemble convertisseur machine.<br />
• Bonne régulation de Vdc <strong>à</strong> l’aide d ’un PI + compensation.
Avantages :<br />
Machine asynchrone <strong>à</strong> double alimentation<br />
(MADA) <strong>sur</strong> le réseau (fortes puissances)<br />
MASDA<br />
P (f)<br />
• transfert bidirectionnel de la puissance rotorique.<br />
fonctionnement hypo. ou hypersynchronisme, hypersynchronisme,<br />
donc large plage de <strong>vitesse</strong>.<br />
• Convertisseur rotorique de faible puissance.<br />
• meilleure utilisation de la machine : puissance élevée au del<strong>à</strong> du synchronisme<br />
possibilité de fonctionner <strong>à</strong> <strong>vitesse</strong> élevée du vent.<br />
• Mode « MPTM » au <strong>des</strong>sous de Ωs s et couple constant au del<strong>à</strong> de Ωs.<br />
gP<br />
(gf)<br />
C1<br />
U<br />
(f)<br />
RESEAU<br />
P(1+g)<br />
(f)
Machine asynchrone <strong>à</strong> double alimentation<br />
(MADA)<br />
Tensions statoriques et<br />
rotoriques<br />
⎧ d<br />
⎫<br />
V ds = R sI ds + ψ ds − θ� ⎪ sψ qs<br />
⎪ dt<br />
⎪<br />
⎨ ⎬<br />
⎪ d<br />
V qs = R sI qs + ψ qs + θ� sψ ⎪<br />
ds<br />
⎪⎩ dt<br />
⎪⎭<br />
⎧ d<br />
⎫<br />
V dr = R r I dr + ψ dr − θ� ⎪ rψ qr<br />
⎪ dt<br />
⎪<br />
⎨ ⎬<br />
⎪ d<br />
V qr = R sI qr + ψ qr + θ� rψ ⎪<br />
dr<br />
⎪⎩ dt<br />
⎪⎭<br />
• Contrôle de la puissance active et réactive<br />
Modèle diphasé :<br />
Flux statoriques et rotoriques :<br />
⎧ψ= L I + M I ⎫<br />
ds s ds dr<br />
⎪<br />
ψ qs L sI qs M I<br />
⎪<br />
⎪ = + qr ⎪<br />
⎨ ⎬<br />
⎪ψdr = L r I dr + M I ds ⎪<br />
⎪ψqr L r I qr M I ⎪<br />
⎩ = + qs ⎭<br />
Equation m écanique :<br />
d Ω<br />
Γ m = Γ e + J + f Ω<br />
dt<br />
ψ<br />
Γ<br />
V<br />
e<br />
ds<br />
ds<br />
=<br />
= ψ<br />
=<br />
−<br />
0<br />
p<br />
s<br />
M<br />
Ls<br />
et<br />
et<br />
I<br />
qr<br />
V<br />
ψ<br />
ψ<br />
qs<br />
qs<br />
ds<br />
=<br />
=<br />
V<br />
s<br />
0<br />
≅<br />
ω<br />
s<br />
ψ<br />
ds
V q<br />
V d<br />
Machine asynchrone <strong>à</strong> double alimentation<br />
(MADA)<br />
M<br />
gω<br />
ψ<br />
1<br />
1<br />
s ds<br />
L s<br />
K inv<br />
+ pT<br />
K inv<br />
+ pT<br />
inv<br />
inv<br />
+<br />
V qr<br />
+<br />
V dr<br />
-<br />
ν 2<br />
ν 3<br />
Modèle interne de la MADA<br />
ν 1<br />
-<br />
+<br />
Rr +<br />
1<br />
p( Lr<br />
−<br />
2<br />
M<br />
)<br />
L<br />
2<br />
M<br />
gωs( Lr<br />
− )<br />
L<br />
gωs( Lr<br />
−<br />
2<br />
M<br />
)<br />
L s<br />
1<br />
Rr + p( Lr<br />
−<br />
2<br />
M<br />
)<br />
L<br />
s<br />
s<br />
s<br />
I qr<br />
I dr<br />
2<br />
V s<br />
L ω<br />
s s<br />
−<br />
−<br />
MV<br />
L<br />
s<br />
MV<br />
L<br />
s<br />
s<br />
s<br />
+<br />
ν 4<br />
+<br />
P<br />
Q
Machine asynchrone <strong>à</strong> double alimentation<br />
(MADA)<br />
( ) ( 4 )<br />
t ( )<br />
t t<br />
⎧ ⎪xd=<br />
i d,v dr ,xdw = Q ref υ3υ ,u= v d,<br />
⎨<br />
⎪⎩ e= Q− Q ref ,y = id Qref<br />
υ3<br />
⎧ ⎛ x�ˆ ⎞ d ⎛ xˆd ⎞ ⎛B1,d ⎞<br />
⎪( θ d) : ⎜ ⎟= ( Ad− LdCy,d) ⎜ u Ldy ˆx<br />
ˆx<br />
⎟+ ⎜ ⎟ +<br />
⎪ ⎜�⎟ dw<br />
dw 0<br />
⎪ ⎝ ⎠<br />
⎝ ⎠ ⎝ ⎠<br />
⎨<br />
−Kd<br />
⎪ �������<br />
⎛ ˆx d ⎞<br />
⎪u= vd =−⎡⎣K1,d<br />
K2,d⎤⎦⎜<br />
ˆx<br />
⎟<br />
⎪⎩ ⎝ dw ⎠<br />
objectif : synthétiser un observateur θ donnant u<br />
tq e=0 quelle que soit w.<br />
Les poles de (A d -L d C y,d ) sont placés en respectant<br />
un horizon de filtrage -ω f =-1/T f et de commande<br />
donnés.<br />
Même démarche pour l’axe q .<br />
ω<br />
u<br />
Q ref<br />
ν4 ν3 V d<br />
Y : me<strong>sur</strong>es, e :signaux d ’erreurs <strong>à</strong> réguler<br />
u : commande, w : perturbations et références<br />
Kinv<br />
1+<br />
pT<br />
P<br />
inv<br />
V dr<br />
-K d<br />
+<br />
-<br />
Xˆ , ˆ<br />
d Xdω 1<br />
2<br />
M<br />
Rr+ p( Lr−<br />
)<br />
L<br />
θd observer<br />
s<br />
I dr<br />
MVs<br />
−<br />
L<br />
s<br />
+<br />
+<br />
Contrôle de la puissance réactive - Schéma standard<br />
Q<br />
+<br />
-<br />
e<br />
y
Test 1 :<br />
Machine asynchrone <strong>à</strong> double alimentation<br />
(MADA)<br />
Régulateurs mono<strong>variable</strong>s utilisés : IP, RST, Observateur<br />
• Poursuite : échelon de puissance active ∆P= P=-5kW 5kW <strong>à</strong> t=1s<br />
suivi d’un échelon de puissance réactive ∆P= P= 2kVar <strong>à</strong> t=2s<br />
• Sensibilité aux perturbations : échelon de <strong>vitesse</strong> ∆Ω=(318 ∆Ω=(318-260)<br />
260) rd/s rd/s<br />
<strong>à</strong> t=3s<br />
Test 2 :<br />
• mêmes échelons en puissances avec variation aléatoire de la <strong>vitesse</strong><br />
• Robustesse vis-<strong>à</strong>-vis vis vis de la résistance rotorique
vsa<br />
isa*20 isa*20<br />
P<br />
Q<br />
Ω<br />
Machine asynchrone <strong>à</strong> double alimentation<br />
(MADA)<br />
Réponses <strong>à</strong> <strong>des</strong> impacts de puissance active, de puissance réactive et de <strong>vitesse</strong><br />
• meilleur découplage avec le RST qu’avec le PI<br />
• temps de réponse : RST : tr=45ms , PI : tr=130ms<br />
• échelon Ω : RST<br />
∆P=150W P=150W (3%) et ∆Q=100Var Q=100Var (5%)<br />
PI<br />
∆P=1300W P=1300W (26%) et ∆Q= Q= 1000Var (50%) pour le PI.
vsa<br />
isa*20 isa*20<br />
P<br />
Q<br />
Ω<br />
Machine asynchrone <strong>à</strong> double alimentation<br />
(MADA)<br />
Performances du système suite <strong>à</strong> <strong>des</strong> variations aléatoires de la <strong>vitesse</strong><br />
• Performances altérées avec le PI.<br />
• Bonnes performances avec le RST : ∆P P et ∆Q Q < 5%.<br />
• Courant stator quasi constant.
Rr<br />
P<br />
Q<br />
Ω<br />
Machine asynchrone <strong>à</strong> double alimentation<br />
(MADA)<br />
Variations aléatoires de la <strong>vitesse</strong> + test robustesser<br />
• Performances médiocress avec le PI.<br />
• Bonnes performances avec le RST : ∆P P et ∆Q Q < 5%.<br />
• Courant stator quasi constant et robustesse via <strong>à</strong> vis de Rr.<br />
Rr
Vvent<br />
Synoptique du simulateur (macroscopique) de<br />
la turbine éolienne<br />
Coefficient de puissance<br />
Cp de l'éolienne<br />
1R<br />
λ<br />
V<br />
Ω<br />
=<br />
K<br />
Ω1<br />
vent<br />
Dispositif mécanique (éolienne)<br />
Calcul du couple<br />
généré par l'éolienne<br />
1 C p<br />
Γ = ⋅π ⋅ ρ ⋅ ⋅R ⋅V<br />
m<br />
2 λ<br />
Γ r<br />
3 2<br />
Ω =Ω<br />
me<strong>sur</strong>é<br />
Γ m<br />
Calcul<br />
du<br />
couple<br />
Pme<strong>sur</strong>é<br />
2<br />
MADA<br />
f<br />
-<br />
Γ<br />
r<br />
Equation<br />
mécanique<br />
1<br />
Js f<br />
+<br />
Ω2<br />
Multiplicateur<br />
de <strong>vitesse</strong><br />
Ω1ref<br />
M<br />
1<br />
K<br />
AC/AC<br />
Vitesse de référence (rad/s).<br />
Ω2ref<br />
Vitesse me<strong>sur</strong>ée (avec l’observateur) (rad/s).
Evolution <strong>des</strong> travaux<br />
� Analyse du fonctionnement de l’ensemble : simulateur - MCC - MADA.<br />
� Validation expérimentale <strong>des</strong> lois de commande mono-<strong>variable</strong>s établies en<br />
simulation pour la commande de la MADA.<br />
� Etude d'une commande robuste de la MADA par une approche multi-<strong>variable</strong><br />
�� Prise en compte de la saturation<br />
�� Comparaison de la MADA avec la MAS <strong>à</strong> cage en prenant en compte les les<br />
considérations de rendement.<br />
�� Thèse en cours : Machines lentes (MRV)
MADA (7.5 kW) MCC (10kW)<br />
Dispositif de<br />
production<br />
d’énergie <strong>à</strong><br />
étudier<br />
Banc d’essai MADA<br />
Moteur<br />
d’entraînement<br />
<strong>des</strong>tiné <strong>à</strong> simuler<br />
le comportement<br />
d’une éolienne<br />
Dispositif de commande<br />
(Onduleur +PC<br />
+ carte d’interface DSP)