04.02.2013 Views

Eoliennes à vitesse variable basées sur des machines

Eoliennes à vitesse variable basées sur des machines

Eoliennes à vitesse variable basées sur des machines

SHOW MORE
SHOW LESS

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)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!