Calcul D'écoulement de Puissance Dans Les Réseaux Électriques Lors de L'intégration Des Fermes Éoliennes Et Des Dispositifs FACTS
Calcul D'écoulement de Puissance Dans Les Réseaux Électriques Lors de L'intégration Des Fermes Éoliennes Et Des Dispositifs FACTS
Calcul D'écoulement de Puissance Dans Les Réseaux Électriques Lors de L'intégration Des Fermes Éoliennes Et Des Dispositifs FACTS
UNIVERSITE DE M’SILA
FACULTE DE TECHNOLOGIE
DEPARTEMENT DE GENIE ELECTRIQUE
MEMOIRE DE FIN D’ETUDES EN VUE DE L’OBTENTION DU
DIPLÔME
DE MASTER EN GENIE ELECTRIQUE
INGENIERIE DES SYSTEMES ELECTROMECANIQUE
THEME
N°d’ordre : 055
Remerciements :
Nous tenons à remercier tout premièrement Dieu le tout
puissant pour la volonté, la santé et la patience, qu'il nous a
donné durant toutes ces longues années.
A toute ma famille
D’électrotechnique
Promotion 2012/2013
Sommaire
Sommaire
Sommaire
Table de figure
Table de tableau
Liste des symboles et abréviations
Introduction Générale............................................................................................................................ 1
Chapitre I :Fonctionnement du Réseaux Electrique
I.1.Introduction .................................................................................................................................... 3
I.2. Définition d’un réseau électrique .................................................................................................. 3
I.3. Classification des réseaux électrique............................................................................................. 4
I.3.1.Description des réseaux électrique .............................................................................................. 4
I.3.1.1.Le réseau de transport THT :.................................................................................................... 4
I.3.1.2.Le réseau de répartition HT :.................................................................................................... 4
I.3.1.3.Le réseau de distribution : ........................................................................................................ 5
I.3.2.Les postes de transformation ....................................................................................................... 5
I.3.2.1.Poste maçonné .......................................................................................................................... 5
I.3.2.2.Poste aérien accroché (ACC) ................................................................................................... 6
I.3.2.3. Postes sous enveloppe préfabriqué .......................................................................................... 7
I.4. Problématique du fonctionnement des réseaux ............................................................................. 7
I.5. L’acheminement de la puissance................................................................................................... 8
I.6. Modélisation des éléments du réseau ............................................................................................ 9
I.6.1. Modélisation d’une ligne de transport ........................................................................................ 9
I.6.2. Modélisation du transformateur ................................................................................................. 9
I.6.3. Modélisation d’une charge ........................................................................................................ 9
I.7. Puissance transmissible dans une ligne ....................................................................................... 10
I.7.1. Chute de tension dans une ligne ............................................................................................... 11
I.8. Conclusion ................................................................................................................................... 11
Chapitre II :Calcul l’écoulement de puissance dans le réseau électrique
II.1.Introduction ................................................................................................................................. 13
II.2. L’écoulement des puissances ..................................................................................................... 13
II.2.1. Définition : .............................................................................................................................. 13
II.2.2. Classification des jeux de barre : ............................................................................................ 13
II.3. Méthodes itératives del’écoulement de puissance: .................................................................... 14
II.3.1.Equation du réseau électrique : ................................................................................................ 14
II.3.2. Méthode de Gauss-Seidel : ..................................................................................................... 15
II.4. le Calcule d’écoulement de puissance: ...................................................................................... 15
II.4.1. Dans les nœuds PQ ................................................................................................................. 16
II.4.2. Dans les nœuds PV ................................................................................................................. 17
II.4.3. Dans le nœud balancier (slack) ............................................................................................... 18
Sommaire
II.5. Puissance dans une ligne ............................................................................................................ 18
II.5.1. La puissance transit dans une ligne :....................................................................................... 18
II.5.2. Les pertes de puissance dans un lignes ................................................................................... 19
II.6. Calcul l’écoulement de puissance par la méthode de Gauss-Seidel : ........................................ 22
II.6.1. Simulations de -6- jeux de barres: .......................................................................................... 22
II.6.2. Simulations de -14- jeux de barres : ....................................................................................... 24
II.6.3. Simulations de -30- jeux de barres : ....................................................................................... 26
II.7. Conclusion: ................................................................................................................................ 29
Chapitre III Modélisation sur les systèmes éoliens
III.1.Introduction ............................................................................................................................... 30
III.2. Généralités sur les éoliennes .................................................................................................... 30
III.2.1. Définition de l’énergie éolienne ............................................................................................ 31
III.2.2. Développement de l’énergie éolienne ................................................................................... 31
III.2.3. Les différents types d’éoliennes ............................................................................................ 31
III.2.3.1. Eoliennes à Axe Horizontal : .............................................................................................. 31
III.2.3.2. Eoliennes à Axe Vertical : .................................................................................................. 32
III.2.4. Fonctionnement selon type de vitesse : ......................................................................... 33
III.2.4.1. Fonctionnement de vitesse fixe : ........................................................................................ 33
III.2.4.2. Fonctionnement de vitesse variable :.................................................................................. 33
III.2.5.Principaux composants d’une éolienne : ................................................................................ 33
III.3. Généralités sur la MADA : ....................................................................................................... 35
III.3.1. Description de la MADA ....................................................................................................... 35
III.3.2.Les formes application des machines asynchrones à double alimentation :........................... 35
III.3.2.1.Application moteur : ............................................................................................................ 35
III.3.2.2.Application génératrice : ..................................................................................................... 35
III.3.3.Domaine d’applications de la GADA :................................................................................... 35
III.3.4.Hypothèses simplificatrices pour la modélisation de la MADA : .......................................... 36
III.3.5.Le principe de machine asynchrone à double alimentation :.................................................. 36
III.4.System de modélisation de la MADA : ..................................................................................... 37
III.4.1. Modèle mathématique de la GADA dans le référentiel généralise : ................................... 38
III.4.2.Equations magnétiques de GADA :........................................................................................ 39
III.4.3.La transformée de Park : ...................................................................................................... 39
III.4.4.Equations électriques et mécanique de la GADA dans leréférentiel (d,q) : .......................... 40
III.4.5.Mise sous forme d’équation d’état ......................................................................................... 42
III.5. simulation de GADA ................................................................................................................ 47
III.5.1. Interprétations des résultats ................................................................................................... 47
III.6.Conversion d’énergie cinétique du vent en énergie mécanique: ............................................... 48
III.6.1 Loi de BETZ ........................................................................................................................... 48
III.6.2.modèle de la turbine d’éolienne.............................................................................................. 49
III.6.2.1Modèle du multiplicateur d’éolienne .................................................................................... 52
III.6.2.1Modèlisation de l’arbre de la turbine .................................................................................... 52
III.7.Caractéristique de puissance ...................................................................................................... 54
III.7.1.Commande en dessous de la puissance nominale (optimisation de la puissance).................. 54
III.7.2.Principe de la MPPT ............................................................................................................... 54
III.7.3.Commande indirecte de vitesse .............................................................................................. 54
Sommaire
III.7.4. Le schéma globale de l’association GADA turbine : ............................................................ 57
III.7.5. Résultats de simulation (en dessous de la puissance nominale) ............................................ 57
III.7.6. Interprétations des résultats ................................................................................................... 59
III.8. Control de puissance active et réactive par le GADA .............................................................. 59
III.8.1. Commande direct ................................................................................................................... 59
III.8.2. Resultats de simulation .......................................................................................................... 59
III.9. Amélioration de la tension de réseau électrique par intégration de puissance éolienne........... 59
III.9.1. Description du réseau étudié.................................................................................................. 63
III.9.2.Resultats de simulation ........................................................................................................... 59
III.10. CONCLUSION ...................................................................................................................... 59
Chapiter IV Les systèmes de transmission à courant alternatif flexible FACTS
IV.1. Introduction .............................................................................................................................. 66
IV.2. Classification des dispositifs FACTS ....................................................................................... 66
IV.3. Les différents types des systèmes FACTS ............................................................................... 67
IV.3.1. Compensateurs parallèles ...................................................................................................... 67
IV.3.1.1. Compensateurs parallèles à base de thyristors ................................................................... 69
IV.3.1.2. Compensateurs parallèles à base de thyristors GTO .......................................................... 71
IV.3.2. Compensateurs série .............................................................................................................. 72
IV.3.2.1. Compensateurs série à base de thyristors ........................................................................... 72
IV.3.2.2. Compensateurs série à base de thyristors GTO .................................................................. 73
IV.3.3. Compensateurs hybrides (série – parallèle)........................................................................... 74
IV.3.3.1. Compensateurs hybrides à base de thyristors ..................................................................... 74
IV.3.3.2. Compensateurs hybrides à base de thyristors GTO............................................................ 74
IV.4. Les avantages de la technologie des dispositifs FACTS .......................................................... 75
IV.5. Amélioration des performances des réseaux électriques par les dispositifs FACTS ............... 76
IV.5.1. Introduction ................................................................................................................................ 76
IV.5.2. La variation de la tension et la puissance sur les réseaux de transport d'énergie
électriques .......................................................................................................................................... 76
IV.5.2.1. Une ligne résistive .............................................................................................................. 77
IV.5.2.2.Une ligne inductive ............................................................................................................. 78
IV.5.2.3.Une ligne inductive relie à deux grands réseaux : ............................................................... 79
IV.6. Exploitation d'un réseau électrique: ......................................................................................... 80
IV.7. Caractéristique d’un compensateur statique TCSC:................................................................. 81
IV.7.1. Etude d’un TCSC ................................................................................................................. 81
IV.7.2. Définition de TCSC:.............................................................................................................. 81
IV.7.3. La structure d’un TCSC ........................................................................................................ 82
IV.7.4. Schéma block globale de TCSC ......................................................................................... 84
IV.7.5. Les avantage de TCSC ......................................................................................................... 84
IV.7.6.Résultat de simulation du TCSC ............................................................................................ 85
IV.7.7.Interprétation .......................................................................................................................... 85
IV.8. Compensateur statique shunt de puissance réactive (SVC) ..................................................... 86
IV.8.1. Définition du SVC ................................................................................................................. 86
IV.8.2.Différents schémas de SVC ................................................................................................... 86
IV.8.3.Modèle de compensateur statique de puissance réactive SVC .............................................. 86
IV.8.4. Les avantages de SVC ........................................................................................................... 88
Sommaire
IV.9. Conclusion ................................................................................................................................ 90
Conclusion générale .......................................................................................................................... 91
Annexe
Références bibliographiques
Résumé
Table de figure
Tableau II.1 Spécification des différentes grandeurs des jeux de barres. [13].......................... 14
Tableau II.2 présentée les données du réseau de -6- jeux de barres .......................................... 23
Tableau II.3 présentée les différents données du réseau de -6- jeux de barres ......................... 23
Tableau II.4 Les tension obtenues par la méthode de GAUSS-SEIDEL .................................. 24
Tableau II.5 Les données des branches du réseau de -14- jeux de barres ................................. 25
Tableau II.6 Les différents données du réseau de -14- jeux de barres ...................................... 25
Tableau II.7 Résultats de simulation du Gauss-Seidel par réseau de -14- jeux de barres ......... 26
Tableau II.8 Les données du réseau de -30- jeux de barres ...................................................... 25
Tableau II.9 Les différents données du réseau de -30- jeux de barres ...................................... 28
Tableau II.10 Résultats de simulation du Gauss-Seidel avec réseau de -30- jeux de barres
................................................................................................................................................... 28
Table de figure
Introduction Générale
Les problèmes d'approvisionnement en énergie rencontrés à l'échelle mondiale ne sont pas
uniquement dus à la destruction de l'environnement et aux changements climatiques qui en
découlent, mais bien aussi à une croissance de la consommation et à la tarification des
ressources fossiles. Une consommation mondiale en hausse, en particulier en ce qui concerne
l'énergie électrique et une tendance à la libéralisation sont à l’origine de changements
considérables, tant dans le domaine des réseaux d'alimentation que pour le remplacement du
parc des centrales électriques. Les conséquences, pour le moins préoccupantes, se répercutent sur
la disponibilité et la qualité de l’énergie. C'est pourquoi les solutions recherchées doivent être
réalisables techniquement parlant tout en répondant aux exigences en matière de durabilité de
productions décentralisées à base d’énergies renouvelables (éolien, photovoltaïque, entre
autres).[1]
Ces productions seront intégrées au réseau, assurer la disponibilité de divers services pour
les gestionnaires de réseau. En d’autres termes, le réseau doit évoluer vers un réseau flexible et
intelligent qui intègre au mieux les énergies locales et/ou renouvelables pour minimiser les
impacts engendrés par l’insertion de productions décentralisées et/ou par la recherche de
nouvelles architectures. Ces deux solutions devraient permettre l’augmentation du taux de
production décentralisée dans le réseau dans les meilleures conditions économiques et de
sécurité.
L’énergie éolienne fut la première des sources d’énergies exploitées par l’homme,
oubliée depuis long temps, cette source d’énergie inépuisable a connu un développement
important depuis le début des années 1990 partout dans le monde. C'est le cas de
l’Allemagne, leader mondial avec une puissance éolienne installée de 22 GW, suivi des
états unis d’Amérique (USA) avec 16 GW et l’Espagne avec 15 GW [2].
Cependant, la technologie des turbines éoliennes semble avoir atteindre la limite en
matière de taille des hélices. En effet, les plus grandes éoliennes actuellement sur le
marcher, ont une hélice de plus 120 m de diamètre et une tour qui dépasse les 100 m de haut,
pour produire une puissance électrique de 6 MW environ [2] .
Les éléments proposés qui permettent ce contrôle amélioré des systèmes sont les dispositifs
FACTS « Flexible AlternatingCurrent Transmission System ». Les dispositifs FACTS font en
général appel à de l'électronique de puissance, des microprocesseurs, de l'automatique, des
télécommunications et des logiciels pour parvenir à contrôler les systèmes de puissance. Ce sont
des éléments de réponse rapide. Ils donnent en principe un contrôle plus souple de l'écoulement
de puissance. Ils donnent aussi la possibilité de charger les lignes de transit à des valeurs près de
leur limite thermique, et augmentent la capacité de transférer de la puissance d'une région à une
autre. Ils Limitent aussi les effets des défauts et des défaillances de l'équipement, et stabilisent le
comportement du réseau.[3]
Le sujet de ce mémoire concerne, en particulier, calcul d’écoulement de puissance dansles
réseaux électriques lors de l’intégration des fermes éoliennes et des dispositifs FACTS. Pour
atteindre ces objectifs de recherche, ce mémoire est organisé en quatre chapitres:
Dans le premier chapitre, nous allons présenter le fonctionnement des réseaux électriqueset
modélisé les différentes composantes des réseaux électrique.
Le second chapitre a été consacré à l’étude de l'écoulement de puissance, où différentes
méthodes itératives classiques à savoir les méthodes de Gauss, Gauss Seidel pour déterminer les
tensions au niveau des jeux de barres d’un réseau électrique
Le troisième chapitre est consacré à des rappels sur les systèmes éoliens, aussi la
modélisation et la simulation de GADA, et avec la simulation de captage du maximum de
puissance possible MPPT, et le contrôle de puissance active et réactive par la commende
vectorielle directe. Dans ce même chapitre nous proposons, une comparaison entre des réseaux
de test sans et avec l’intégration des fermes éoliennes.
Le quatrième chapitre est consacré à l'étude théorique des différents systèmesFACTS en
générale, à savoir leurs classifications et leur principe de fonctionnement, et aussi étudie la
variation de la tension et la puissance sur les réseaux de transport d'énergie électriques, et on
présente les différentes configurations du SVC et TCSC sont présentées ainsi que sa simulation,
sur le logiciel MATLAB.
Nous terminerons ce présent travail par une conclusion générale et quelques perspectives de
recherche envisagées.
Les réseaux de répartition sont à haute tension (HT) de l'ordre de 50 à 150 kV (on
rencontre le plus souvent le 63 et le 90 kV), ont pour but d'assurer à l'échelle régional
la fourniture d'électricité.
I.3.1.3.Le réseau de distribution :
Les réseaux de distribution sont généralement basés sur une structure arborescente de
réseau à partir d'un poste source, l'énergie parcourt l'artère ainsi que ses dérivations
avant d'arriver aux postes de transformation MT/BT. Les réseaux de distribution ont pour
but d'alimenter l'ensemble des consommateurs. Il existe deux sous niveaux de tension [5] :
- Les réseaux à moyenne tension MT (de 3 à 33 kV).
- Les réseaux à basse tension BT (de 220 à 380 V), sur lesquels sont raccordés les
utilisateurs domestiques.
La distribution ce fait souvent soit par des lignes aériennes ou des câbles sous terrain. Les
réseaux de distribution ont principalement une structure radiale.
Contrairement aux réseaux de transport , les réseaux de distribution présentent une
grande diversité de solutions techniques à la fois selon les pays concernés, ainsi que selon
la densité de population..
a)-les Réseaux de distribution à moyenne tension(MT)
On appelle réseaux MT les réseaux couvrants la plage de 1 à 36 kV.
b)-les Réseaux de distribution à basse tension(BT)
Permet d’acheminer l’énergie électrique du réseau de distribution MT aux points de faible
consommation. Ce réseau alimente un nombre très élevé de consommateurs du domaine
domestique. La structure de ce réseau est de type aérien ou souterrain.
I.3.2.1.Poste maçonné
C’est une construction maçonnée qui contient un transformateur, un jeu de barres pour
donner un départ vers une autre zone de consommation et un répartiteur ou une armoire pour
donner plusieurs voies de distribution BT avec un dispositif de protection. (Voir la figure I.3).
Une ligne de transmission, reliant deux nœuds i et k, est habituellement modélisée par un
circuit équivalent, comme le montre la figure (I.6) dont les quatre paramètres de la ligne de
transmission sont réparties uniformément le long de la ligne.
I.6.2.Modélisation du transformateur
Généralement le modèle d’un transformateur est donné par une réactance inductive.
Le modèle équivalent d’une charge est représenté par une admittance Y définie par :
(I.1)
: Admittance de la charge.
(I.2)
Cette valeur maximale est atteinte pour θ = 45°.Si l’on maintient constante, ce qui implique
que n’est plus nulle, mais ne remet pas en cause l’expression (I.2), on a, pour = :
(I.3)
Pour améliorer la capacité de transfert des réseaux, on voit donc qu’il est essentiel de disposer du
plus grand nombre possible de points à tension fixée.
Concéderons la figure (I.4) qui représente une ligne de transport d’impédance complexe
, et que la tension n’est tenue qu’à l’extrémité 1 (coté générateur), l’extrémité 2(Coté charge)
absorbant une puissance .
Si le réseau n’est pas trop chargé, le diagramme de tension donné par la figure (I.9) conduit à
assimiler la chute de tension ∆V à :
∆V≈ - (I.4)
(I.5)
(I.6)
θ (I.8)
Dans ces conditions, les relations (I.10) et (I.11) illustrent le fait que :
- La chute de tension dépend principalement de la puissance réactive consommée par
l’extrémité réceptrice ;
- L’angle de transport θ dépend principalement de la puissance active transmise.
I.8. Conclusion
Dans ce chapitre nous avons présenté un aperçu sur le fonctionnement des réseaux d’énergie
électrique on commençant par définition des réseaux, leurs classifications des réseaux électrique
et leurs architectures et l’acheminement de puissance à travers ces réseaux. On a vu brièvement
les facteurs qui limitent la puissance transmissible dans les lignes de transport et la chute
de tension dépend de puissance réactive consommée.
Chapitre II :
Tableau II.1 Spécification des différentes grandeurs des jeux de barres. [13]
(II.1)
Où :
(II.2)
[ ]
(II.3)
[ ]
∑
∑
[ ] (II.4)
∑
3. Matrice complexe.
5. Chaque élément de la diagonale (k, k) est la somme des admittances de toutes les
branches liées au jeu de barres k.
: L’impédance des jeux de barres
(II.5)
(II.6)
(II.7)
(II.8)
(II.9)
[ ] (II.10)
Pour un réseau de N jeux de barres, la tension calculée à n’importe quel jeu de barre k, où
GAUSS
() ( )
[ ( ) ∑ ] (II.11)
GAUSS-Seidel
() () ( )
[ ( ) ∑ ∑ ] (II.12)
( ) (II.14)
Dans chaque itération la tension est déterminée pour tous les jeux de barres sauf le <<SWING
BUS>>. Le processus est répété plusieurs fois jusqu’à ce que les tensions dans deux itérations
successives ne se diffèrent pas l’une de l’autre de plus que la tolérance. Ce processus de résoudre
les équations linéaires est connu sous le nom de la méthode itérative de GAUSS Seidel. Si le
même ensemble des valeurs des tensions est utilisé dans une itération complète (au lieu de
substituer les nouvelles valeurs chaque fois qu’elles sont déterminées pour calculer la tension du
jeu de barre suivant), le processus serait appelé la méthode itérative de Gauss. [8][35]
En se basant sur la solution de l’itération (k), le calcul des nouvelles tensions à l’itération
(k+1) dépend du type de nœud concerné.
Il faut noter que le nœud 1 est choisi comme nœud balancier dont l'amplitude et la phase de
tension sont connues. Cela explique pourquoi on commence les calculs pour i = 2.
(II.16)
D'où, on tire:
( ) ( )
[∑ ] (II.17)
Il faut remarquer que les meilleures valeurs de tension sont utilisées dans le calcul de la
puissance réactive .Une fois est déterminée, on utilise l'équation (II.15) pour calculer la
nouvelle valeur de du nœud de génération i et en tirer la phase θ i. Habituellement, une limite
maximum et/ou minimum :
(II.19)
Q doit être spécifié. Si excède la limite spécifiée , cette dernière est alors choisie à la
place de la valeur calculée pour être utilisée dans (II.15). Si l'amplitude de tension | | de la
nouvelle valeur calculée de Vi est plus grande que la valeur spécifiée d'amplitude de tension du
nœud de génération i, la nouvelle tension est corrigée en la multipliant par le facteur
correspondant au rapport de l'amplitude de tension spécifiée par l'amplitude de tension calculée,
tout en gardant la phase de tension de la valeur calculée. En d'autres termes, seule l'amplitude
doit être corrigée.
En résumé, le processus itératif de Gauss-Seidel commence en prenant des valeurs initiales
pour les amplitudes et phases de tension inconnues (sauf pour le cas du nœud balancier) et
calcule les nouvelles valeurs, qui sont ainsi des valeurs corrigées. Au niveau de chaque nœud, et
à chaque itération, la tension corrigée est remplacée dans l'équation (II.15) comme une valeur
estimée afin de calculer une nouvelle valeur. Ce processus est répété jusqu'à convergence.[8]
Une fois que toutes les tensions aux nœuds , ,..., Vn sont calculées, la puissance du nœud
balancier est alors calculée:
(II.20)
(II.21)
B/2 B/2
( ) (II.22)
Où
admittance de la ligne ij,
admittance shunt total de la ligne ij.
La puissance complexe transité le long de la ligne ij s'écrit:
(II.23)
En substituant l'expression de (1.17) dans (1.18), on obtient:
( ) ( ) (II.24)
( ) ( ) (II.25)
Début
I=1,……, N Branches, I S
Itération = 1,Maxiter
Non Oui
Qi Qi Qi Qi
(max) (min)
Oui Non
Calculer l’erreur △V
2 1
Oui
Iter Maxiter
et V˃
Non
Accélération la convergence de la tension
Fin
Charge P,Q
Charge P,Q
5
4 P, Q
P, Q
6
Charge P,Q
Les données du réseau à -6- jeux de barres sont présentées par le tableau (II.2) suivant:
La tension de base est égale 220KV, la puissance de base est égale 100MVA
Y=
14
10 8 11 12
13
9 7 6
4 5
3 2 1
Les données du réseau à -14- jeux de barres sont présentées par le tableau (II.5) suivant:
Nœuds de Nœuds de R (pu) X (pu) B/2(pu)
départ d’arrivée
1 2 0.01938 0.05917 0.0264
1 5 0.05403 0.22304 0.0246
2 3 0.04699 0.19797 0.0219
2 4 0.05811 0.17632 0.0170
2 5 0.05695 0.17388 0.0173
3 4 0.06701 0.17103 0.0064
4 5 0.01335 0.04211 0.0
4 7 0.0 0.20912 0.0
4 9 0.0 0.55618 0.0
5 6 0.0 0.25202 0.0
6 11 0.09498 0.19890 0.0
6 12 0.12291 0.25581 0.0
6 13 0.06615 0.13027 0.0
7 8 0.0 0.17615 0.0
7 9 0.0 0.11001 0.0
9 10 0.03181 0.08450 0.0
9 14 0.12711 0.27038 0.0
10 11 0.08205 0.19207 0.0
12 13 0.22092 0.19988 0.0
13 14 0.17093 0.34802 0.0
Tableau II.5Les données des branches du réseau de -14- jeux de barres
Les résultats de simulationobtenus par la méthode de Gauss-Seidel sont donnés dans lestableaux
suivants :
Le nombre d’itération = 69
Nœud de Jeu Les tensions L’angle
de barres
1 1.06 0
2 1.045 -4.9575
3 1.0028 -12.614
4 1.0147 -10.230
5 1.0214 -8.7728
6 1.029 -14.81
7 1.0212 -13.476
8 1.0610 -13.476
9 1.0010 -15.216
10 0.9982 -15.451
11 1.0100 -15.266
12 1.0128 -15.723
13 1.0069 -15.772
14 1.0145 -16.561
Tableau II.7 Résultats de simulation du Gauss-Seidel pourle réseau de -14- jeux de barre
23 3 1 0 0 0 0.032 0.016 0 0
24 3 1 0 0 0.043 0.087 0.067 0 0
25 3 1 0 0 0 0 0 0 0
26 3 1 0 0 0 0.035 0.023 0 0
27 3 1 0 0 0 0 0 0 0
28 3 1 0 0 0 0 0 0 0
29 3 1 0 0 0 0.024 0.009 0 0
30 3 1 0 0 0 0.106 0.019 0 0
Tableau II.9lesdifférents données du réseau de -30- jeux de barres
Interprétation
Dans la simulation de -30-jeux des barres on remarque les tensions aux niveaux des bus (14 et
15) sont (0.9765 pu =214.83Kv) et (0.9647 pu =212.23Kv) respectivement, sont très inferieur à
la tension nominale, dont l’intérêt de l’amélioration par le système de l’éolien.
II.7. Conclusion:
L’objectif principal visé dans ce chapitre c’est d’introduire l’étudiant au problème de calcul
d’écoulement de charge dans un réseau électrique. Nous avons vu en modélisant tous les
éléments du réseau que le problème d’écoulement de puissance et déterminer les valeurs de
tension et la phase dans chaque nœud par utiliser la méthode de GAUSS-SEIDEL, pour ce la des
techniques d’amélioration sont nécessaire, tel que les dispositifs FACTS et les éoliennes dont
nous utiliserons dans la suite de ce travail, afin d’améliorer le transit de puissance dans les
réseaux électriques.
Les éoliennes se divisent en deux grandes familles : les éoliennes à axe vertical qui tendent à
disparaître, les éoliennes lentes à axe horizontal utilisées surtout pour le pompage et les
éoliennes rapides à axe horizontal plus particulièrement utilisées dans la production d’énergie
électrique.
Les éoliennes à axe horizontal beaucoup plus largement employées, elles nécessitent
souvent un mécanisme d’orientation des pales, présentant un rendement aérodynamique plus
élevé, démarrent de façon autonome et présentent un faible encombrement au niveau du sol.
Dans ces types d’éolienne l’arbre est parallèle au sol. Le nombre de pales utilisé pour la
production d’électricité varie entre 1 et 3, le rotor tripale étant le plus utilisé car il constitue un
compromis entre le coefficient de puissance, le coût et la vitesse de rotation du capteur éolien, ce
type d’éolienne a pris le dessus sur celles à axe vertical car elles représentent un coût moins
important.[23]
agricoles. Leurs nuisances sonores sont de plus relativement faibles, le niveau sonore d’une
éolienne est de 50Db à 150 mètres et devient imperceptible au-delà de 400 mètres. Dans la
plupart des cas, le bruit du vent est supérieur à celui engendré par l’éolienne.
Le fonctionnement à vitesse fixe ou à vitesse variable constitue l’une des plus importantes
spécifications des aérogénérateurs, Nous allons respectivement citer en ce qui suit les avantages
de l’un par rapport à l’autre :
Le mât :
Généralement un tube d’acier ou éventuellement un treillis métallique, doit être le plus haut
possible pour éviter les perturbations près du sol. Toutefois, la quantité de matière mise en œuvre
représente un coût non négligeable et le poids doit être limité. Un compromis consiste
généralement à prendre un mât de taille très légèrement supérieure au diamètre du rotor de
l’aérogénérateur.[21]
La nacelle :
Regroupe tous les éléments mécaniques permettant de coupler le rotor éolien au générateur
électrique arbre lent et rapide, roulements, multiplicateur. Le frein à disque, différent du frein
aérodynamique, qui permet d’arrêter le système en cas de surcharge. Le générateur qui est
généralement une machine synchrone ou asynchrone et les systèmes hydrauliques ou électriques
d’orientation des pales (frein aérodynamique) et de la nacelle (nécessaire pour garder la surface
balayée par l’aérogénérateur perpendiculaire à la direction du vent). A cela viennent s’ajouter le
système de refroidissement par air ou par eau, un anémomètre et le système électronique de
gestion de l’éolienne.[20]
Le rotor :
Formé par les pales assemblées dans leur moyeu. Pour les éoliennes destinées à la
production d’électricité, le nombre de pales varie classiquement de 1 à 3, le rotor tripale (concept
danois) étant de loin le plus répandu car il représente un bon compromis entre le coût, le
comportement vibratoire, la pollution visuelle et le bruit.
III.3.2.1.Application moteur :
III.3.2.2.Application génératrice :
C’est surtout le mode de fonctionnement en génératrice qui nous intéresse car il correspond
parfaitement dans notre cas, à la MADA utilisée dans une éolienne, la puissance est fournie
par le dispositif d’entraînement, dans notre cas la turbine éolienne ; une partie de la
puissance transitant par le stator est réabsorbée par le rotor.
Pour cette machine, le rotor est bobiné et comporte des enroulements triphasés couplés en
étoile associés à des bagues collectrices qui permettent d’alimenter ses enroulements. Le stator
est identique à celui d’une machine asynchrone à cage. [28]
[ ] [ ][ ] [ ]
{ (III.1)
[ ] [ ][ ] [ ]
[ ] [ ] [ ] [ ]
[ ] [ ] [ ]= [ ]
Phase.
[ ]= [ ]; [ ]=[ ]
Où
, sont respectivement les inductances principales de fuite statoriques et rotoriques ;
Ms, Mr : sont respectivement les inductances mutuelles statoriques et rotoriques.
La matrice des inductances mutuelles (ou matrice de couplage rotor- stator) est donnée par :
( ) ( ) ( )
[ ] [ ] . ( ) ( ) ( ) (III.3)
[ ( ) ( ) ( ) ]
Appelée souvent transformation des deux axes , cette transformation a été présente par
PARK en 1929 lors de l’étude de la machine asynchrone comme moyen de simplifier
les équations de celle-ci . fait correspondre aux variables réelles leurs composantes :
* D’axes directs (indice d) ;
* D’axes en quadrature (indice q). [20]
L'idée de Park repose sur le fait qu'un champ tournant crée par un systéme triphasé peut
l'étre aussi à l'idéntique par deux courants Id ,et Iq ,et deux tension . la transformation
qui traduit ce passage est donné par :
[ ] [ ( )] [ ] [ ] [ ( )] [ ]
[ ( )] est la matrice de Park :
( ) ( ) ( )
[ ( )] √ * ( ) ( ) ( ) (III.4)
[ √ √ √ ]
Pour la transformation inverse , c'est-à-dire pour le passage du système biphasé au système
triphasé, on utilise la matrice de Park inverse suivante :
Université de m’sila 2012/2013 Page 39
Chapitre IIIModélisation sur les systèmes éoliens
( ) ( )
√
[ ( )] √ * ( )) ( ) (III.5)
√
[ ( ) ( ) ]
√
( ) ( ) ( )
[ ] √ * ( ) ( ) ( ) *[ ] (III.6)
[ √ √ √ ]
Courant :
( ) ( ) ( )
[ ] √ * ( ) ( ) ( ) *[ ] (III.7)
[ √ √ √ ]
(III.8)
{
Les composantes des flux statoriques et rotoriques sont données par :
(III.9)
{
Avec
, =
Avec respectivement:
, , et les tensions statoriques et rotoriques directs en système diphasé.
, , et lescourantsstatoriques et rotoriques directs en système diphasé.
, , et les flux statoriques etrotoriques directs en système diphasé.
, sont les pulsations rotoriques ; des grandeurs électriques statoriques et rotoriques.
L’expression du couple électromagnétique est donnée par les relations suivantes :
( )
{ (III.10)
( )
(III.11)
Avec respectivement :
: Couple électromagnétique.
: Couple résistant.
: Le coefficient de frottement visqueux de la GADA ;
: La vitesse de rotation de l’axe de la GADA ;
: L’inertie des parties tournantes.
Les puissances active et réactive statoriques sont représentées par :
( )
{ (III.12)
( )
: La puissance active statoriques.
: La puissance réactive statoriques.
Les puissances active et réactive rotoriques sont représentées par :
Université de m’sila 2012/2013 Page 41
Chapitre IIIModélisation sur les systèmes éoliens
( )
{ (III.13)
( )
: La puissance active rotoriques.
: La puissance réactive rotoriques.
Pour la génératrice asynchrone à double alimentation les variables de contrôle sont les
tensions statoriques et rotoriques. En considérant les flux statoriques et les courant rotoriques
comme des vecteurs d’état, alors le modèle de la MADA est décrit par l’équation d’état suivante:
[ ̇] [ ][ ] [ ][ ] (III.14)
Pour la machine asynchrone à double alimentation les variables de contrôle sont les tensions
statoriques et rotoriques. En considérant les
Courant statoriques et les flux rotoriques comme des vecteurs d’état, alors le modèle de la
MADA est décrit par l’équation d’état suivante :
[ ] [ ] (III.15)
[ ] [ ] (III.16)
A partir d’expressions (II.9), on déduit les expressions des flux rotoriques :
{ (III.17)
{ (III.18)
( ) ( ) (III.19)
( ) ( )
{ ( )
( )
( )
(III.20)
{
La matrice d’état s’écrit :
( )
( )
[ ] (III.21)
[ ]
La matrice des entrées s’écrivent:
[ ] (III.22)
[ ]
On rappelle que :
( )
{
Pour facilité la simulation, nous avons écris la matrice [A]comme suit :
[ ] [ ] [ ] [ ] (III.23)
[ ] ;[ ] [ ] [ ] [ ]
[ ]
A partir des équations (II, 18) et (II, 19) (II, 20) le modèle peut être exprimé dans l'espace
d'état sous la forme suivante :
[ ̇] [ ][ ] [ ][ ]
( )
( )
[ ]= [ ]+
[ ]
(III.24)
[ ]
[ ]
60 80
courant rotorique Iqr [A]
20 40
0 20
-20 0
-40 -20
-60 -40
0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1
Temps (s) Temps (s)
60 100
40 80
20 60
0 40
-20 20
-40 0
-60 -20
0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1
Temps (s) Temps (s)
0 0.2
flux rotorique fdr (wb)
0
-0.5
-0.2
-0.4
-1
-0.6
-1.5 -0.8
0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1
temps (s) temps (s)
0
0.5
flux statorique fds (wb)
-0.5 0
-1 -0.5
-1.5 -1
-2 -1.5
0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1
temps (s) temps (s)
20
0
couple électromagnétique [N*m]
-20
-40
-60
-80
-100
-120
-140
-160
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
temps [s]
1500 500
1000
0
500
-500
-1000
-500
-1000 -1500
-1500
0 0.2 0.4 0.6 0.8 1 -2000
0 0.2 0.4 0.6 0.8 1
temps [s]
temps [s]
4 4
x 10 x 10
2.5 3.5
puissance active statorique [wat]
2 3
puissance réactive statorique [var]
1.5 2.5
1 2
0.5 1.5
0 1
-0.5 0.5
-1 0
-1.5 -0.5
-2 -1
0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1
temps [S] temps [S]
III.6. Conversiond’énergiecinétiqueduventenénergiemécanique:
III.6.1 LoideBETZ
La théorie du moteur éolien à axe horizontal a été établie par BETZ. Il suppose que le
moteuréolien est placé dans un air animé à l’infini amont d’une vitesse et à l’infini aval d’une
vitesse : [20][25]
( )
En déduit :
(III.29)
ρ. ( ) ( ) (III.31)
Un vent théoriquement non perturbé traverserait cette même surface Ssans diminution de
Vitesse soit à la vitesse, la puissance correspondante(mécanique théorique) serait alors:
(III.32)
On définit le rapport entre les deux puissances Puissance extraite du vent et celle qui est
théoriquement disponible qui est appelé le coefficient de puissance :[28]
( )( ( ) )
(III.33)
Ce coefficient présent un maximum de 16/27 soit 0,59. C'est cette limite théorique appelée
limite de Betz qui fixe la puissance maximale extractible pour une vitesse de vent donnée.
Cette limite n'est en réalité jamais atteinte et chaque éolienne est définie par son propre
coefficient de puissance exprimé en fonction de la vitesse relative.[24][28]
L’éolienne capte l’énergie cinétique du vent et la convertit en un couple qui fait tourner les
pales du rotor. Trois facteurs déterminent le rapport entre l’énergie du vent et l’énergie
mécanique récupérée par le rotor: la densité de l’air, la surface balayée par le rotor et la vitesse
du vent. La densité de l’air et la vitesse du vent sont des paramètres climatologiques qui
dépendent du site.[22]
Le ratio de la vitesse est défini comme le rapport entre la vitesse linéaire des pales et le
vent : (III.36)
Avec
: Le coefficient de performance ou coefficient de puissance.
: Le ratio de la vitesse.
: L’angle de l’orientation de la pale.
: La longueur de la pale.
: La vitesse de rotation de la turbine.
Ainsicommeilestillustrédanslafigure (III-19),ilestpossibledeconstatéquele
coefficientdepuissanceévolueenformedecloche,lemaximumdecelle-ci( )est
atteintpourunevaleurdurapport d’avance( )etdel’angled’orientationdespales
( )Aveccettevaleurlaturbinefonctionneaveclerendementmaximumthéorique.Dans
lasuitedenotretravailunecommandeadaptéeseraélaboréepouratteindrecepointdefonctionnement.[2
1]
0.6
beta = 0°
0.55 beta = 5°
beta = 10°
0.5 beta = 15°
beta = 20°
0.45 beta = 25°
Coefficient de puissance (cp)
beta = 30°
0.4
0.35
0.3
0.25
0.2
0.15
0.1
0.05
0
0 2 4 6 8 10 12 14 16 18
Ratio de vitesse spécifique (lambda)
La puissance mécanique récupérée par une turbine éolienne peut s’écrire sous la forme:
( ) (III.37)
= * (III.38)
et
S= * (III.39)
avec
:Le couple aérodynamique
:La puissance aérodynamique
ρ : la densité de l’air.
v : la vitesse du vent.
S : la surface circulaire balayée par la turbine, le rayon du cercle est déterminé par la
longueur de la pale.
(III.40)
Avec :
J (III.42)
1 /G
* 1 /G
Zone I: Correspond aux faibles vitesses du vent, la puissance disponible dans cette zone est
inferieure à la puissance nominale de la turbine. L’objectif dans cette zone est d’extraire le
maximum de puissance du vent en appliquant des techniques appelées techniques d’extraction de
maximum de puissance (Maximum Power Point Tracking).
Zone II: Zone de transition entre la zone de charge partiel (Zone I) et la zone de charge
nominale (Zone II).Dans cette zone, la vitesse de rotation est maintenue constante pour des
raisons mécaniques (contraintes des pales) et acoustiques (bruits généré).
Zone III: Correspond aux vents forts, l’objectif dans cette zone est de limiter la puissance
produite à une valeur égale à la puissance nominale de l’éolienne pour éviter les surcharges .Cela
se fait par action sur l’angle de calage des pales.
III.7.2.PrincipedelaMPPT
La MPPT (Maximum PowerPoint Tracking) est une méthode de contrôle très fiable,
Robuste et simple à implanter. Elle consiste à déterminer la vitesse de la turbine qui permet
d’obtenir le maximum de puissance générée. Pour agir sur la vitesse de rotation de la turbine ,
on s’intéresse au sens de variation de la grandeur .Ce gradient vaut zéro lorsque l’on a atteint
le maximum de la puissance, point que l’on recherche .On le déduit à partir l’équation :
( ) (III.44)
La recherche du maximum de puissance peut se faire par action sur le couple (commande
indirect de la vitesse)
Le résultat de simulation montant le principe de MPPT est présenté sur la figure (III.23) :
4000
vent=5m/s
vent=6m/s
vent=7m/s
vent=7.7m/s
puissance optimale
3000
Puissance mécanique [W]
2000
1000
0
0 200 400 600 800 1000 1200 1400 1600 1800 2000
Vitesse de rotation du générateur [tr/min]
Figure III.23 Puissance éolienne en fonction de la vitesse de rotation pour différents vitesses du
vent.
( )
= (III.45)
: Le couple optimale.
Une estimation de la vitesse du vent est souvent utilisée pour contourner ce problème.
Connaissant la vitesse de rotation de la turbine, la vitesse du vent sera déduite de l’équation de la
vitesse spécifique comme suit.
(III.47)
( )
= (III.48)
- - =0 (III.52)
Avec :
(III.53)
D’où
= (III.54)
Cette expression donne le couple qui doit être imposé à la génératrice pour assurer le
fonctionnement optimale de l’éolienne. Le schéma block de cette structure de commande est
donné par la figure ci-dessous
8
Vitesse du vent (m/s)
7,5
7
6,5
6
5,5
5
4.5
4
0 5 10 15 20 25 30 35 40 45 50
Temps (s)
4000 1200
3500
1000
la puissance produit (w)
3000
2000 600
1500
400
1000
200
500
0 0
0 10 20 30 40 50 0 10 20 30 40 50
Temps (s) Temps (s)
0.48
La coefficient de puissance
0.475
0.47
0.465 Cp
Cp max
0.46
0.455
5 10 15 20 25 30 35 40 45 50
Temps
FigureIII.28 La courbe représenté lecoefficientdepuissance .
10
La vitesse spécifique
6
la vitesse spécifique
la vitesse spécifique optimale
4
0
0 5 10 15 20 25 30 35 40 45 50
Temps (s)
{ (III.55)
En introduisant ces expressions dans les équations des tensions rotoriques du système sont :
( ) ( )
{ (III.56)
( ) ( )
des courants rotoriques diphasés dans le système (III.56) disparaissent en régime permanent.
Nous pouvons donc écrire :
( )
{ (III.57)
( )
En revanche, le terme représente une force électromotrice dont l’influence n'est pas
négligeable, le système de contrôle devra donc compenser cette perturbation. Ainsi, on obtient
un modèle plus simple permettant le contrôle direct et indépendant des puissances active et
réactive en utilisant un seul régulateur sur chaque axe (Figure IV.2).[28][20][24]
-
+ +
PI
[GADA]
+
Partie de
+ +
PI Puissanc
- - e
0 1500
-1000
1000
-2000
500
-3000
-4000 0
0 1 2 3 4 5 0 1 2 3 4 5
temps (s) temps (s)
4 4
x 10 x 10
3 4
Figure III.31La puissance de référence (P active ; Q réactive)
2
4
3
x 10
2.5
2x 104
( messtatorique
1 Q mes P mes
/ ref ) (var)
Q réf
puissance active ( mes / ref ) (w) 2 P réf
2
0 1.5 0
1
1 -1
puissacce
0.5-2
-2
puissance réactive
0 0
-0.5-4
-3
-1 0 1 2 3 4 5 -1
0 1 2 3 4 5
temps (s) temps (s)
-1.5
-2
-2
-3 -2.5
0 1 2 3 4 5 0 1 2 3 4 5
Temps (s) Temps (s)
2000
Q mes P mes
puissance réactive ( mes / ref ) (var)
500
puissance active ( mes / ref ) (w)
Q réf P réf
1500 0
-500
1000 -1000
-1500
500 -2000
-2500
0
-3000
-3500
-500
-4000
-4500
1 1.5 2 2.5 3 1 1.5 2 2.5 3
Temps (s) Temps (s)
100 100
50 50
0 0
-50 -50
-100 -100
0 1 2 3 4 5 0 1 2 3 4 5
Temps (s) Temps (s)
100 100
Le courant Idr (A)
0 0
-50 -50
-100 -100
0 1 2 3 4 5 0 1 2 3 4 5
Temps (s) Temps (s)
Figure III.34 Lescomposantesducourantstatoriquesetrotorique
5
Le flux statorique Ids (w)
0.5
-0.5
0 -5
-1
-0.5 -10
-1.5 -15
-1
-20
-2 -1.5 -25
0 1 2 3 4 5 0 2 4 5
Temps (s) Temps (s)
100
Le couple électromagnétique (N.m)
50
-50
-100
-150
-200
-250
0 2 4 5
Temps (s)
Interprétation du résultat
Dans cette application, on fait l’intégration de puissance éolienne dans le réseau électrique
de -30 – jeux de barres représenté dans la figure (III.37), afin d’améliorer 1a tension au niveau
des jeux de barres (28 et 29).
Figure III.37 Schéma unifilaire du réseau électrique 30-bus avec intégration de la puissance
éolienne.[29]
III.9.2.Simulation:
Les résultats de simulation effectué sous programme sont présentés dans le tableau suivant :
Le nombre d’itération = 191
Nœud de Les tensions Les tensions Nœud de Les tensions Les tensions
Jeu de sans avec Jeu de sans avec
barres l’intégration l’intégration barres l’intégration l’intégration
1 1.0600 1.0600 16 1.0180 1.0279
2 1.0430 1.0430 17 1.0146 1.0254
3 1.0268 1.0316 18 1.0021 1.0128
4 1.0192 1.0192 19 0.9999 1.0109
5 1.0100 1.0100 20 1.0044 1.0154
6 1.0143 1.0218 21 1.0049 1.0170
7 1.0048 1.0092 22 1.0113 1.0255
8 1.0085 1.0167 23 1.0044 1.0168
9 1.0337 1.0436 24 0.9982 1.0165
10 1.0209 1.0321 25 0.9920 1.0245
11 1.0800 1.0894 26 0.9738 1.0069
12 1.0286 1.0375 27 0.9969 1.0379
13 1.0603 1.0690 28 1.0126 1.0241
14 1.0145 1.0241 29 0.9765 1.0183
15 1.0111 1.0215 30 0.9647 1.0070
On remarque après l’intégration l’énergie éolienne au niveau des jeux de barres de réseau
électrique,La tension au niveau de jeux de barres (28 et 29) augmente à (1.0183pu et 1.0070
pu).ce qui montre l’efficacité de l’éolienne pour l’amélioration de la tension.
III.10. Conclusion:
Dans cette chapitre, on a présenté un état de l’art de éolien, et on a rappelé la structure de
la GADA, ainsi que son modèle mathématique dans le repère (d,q), pour contourner cette
difficulté, on a fait appel à la transformation de Park Le modèle dynamique de la machine écrit
dans le repère de Park puis traduit sous forme d’équations d’état, a été simulé. Les résultats de
simulation sont satisfaisants.
Aussi on a présentéle contrôle de puissance active et réactive dans la GADA par la
commande vectorielle directe et l’intégration de puissance éolienne pour amélioration la tension
de réseau électrique de 30 jeux de barres.
Chapitre IV
Off). Ces éléments sont assemblés pour former les convertisseurs de tension ou de courant
afin d’injecter des tensions contrôlables dans le réseau.
3.Une troisième génération de FACTS utilisant des compensateurs hybrides et qui est adaptée à
chaque cas. Contrairement aux deux premières générations, celle-ci n’utilise pas de dispositifs
auxiliaires encombrants tels que des transformateurs pour le couplage avec le réseau [12].
STATCOM
SVC
TCVR
TCPAR
TCSC
SSSC
Figure IV.1 Diagramme schématique des applications des FACTS dans les réseaux.
Selon la demande des consommateurs, les lignes électriques doivent transporter des
puissances actives et aussi des puissances réactives.
Vers la fin des années 60 plusieurs équipements utilisant l’électronique de puissance ont fait
leurs apparitions. Ces derniers avaient l’avantage d’éliminer les parties mécaniques et d'avoir un
temps de réponse très court. Ces équipements étaient constitués essentiellement d’une inductance
en série avec un gradateur. Le retard à l’amorçage des thyristors permettait dérégler l’énergie
réactive absorbée par le dispositif.[10]
Afin d’éviter des pertes supplémentaires à cause de la transmission du courant réactif et pour
augmenter la stabilité des réseaux interconnectés il est nécessaire de compenser la puissance
réactive au niveau des sous stations d’interconnexion. L’avantage de ces dispositifs permet
d’éliminer les parties mécaniques de rendement mauvais et de réaction assez lente, et d’avoir un
temps de réponse très court .Ces équipements étaient constitués essentiellement d’une inductance
en série avec un gradateur, le retard à l’amorçage des thyristors permettait de régler l’énergie
électrique réactive absorbée par l’inductance du dispositif
•Système TSC
Le schéma du TSC consiste à introduire les condensateurs comme des sources de
puissance réactive. Chaque source monophasée est formée d’un condensateur (C) en série avec
deux interrupteurs statiques bidirectionnels (figure IV.3).
Ces compensateurs sont connectés en série avec le réseau et peuvent être utilisés comme une
impédance variable (inductive, capacitive) ou une source de tension variable. En général, ces
compensateurs modifient l’impédance des lignes de transport en insérant des éléments en série
avec celles-ci.
IV.5.1. Introduction
Comme déjà mentionné précédemment, notre travail consiste à optimiser la marge de la
stabilité statique de tension ainsi que la réduction des pertes actives.
De ce fait, il faut choisir, parmi les équipements FACTS, ceux qui ont des applications
spécifiques tel le maintien de la tension et le contrôle de l’écoulement de puissance : le SVC
étant le plus efficace dans la compensation de l’énergie réactive ainsi que le maintien de la
tension et le TCSC est le meilleur choix pour un contrôle approprié du transit de puissance et par
suite la réduction des pertes actives dans le réseau.
de ce fait, nous avons opté, dans la suite de notre travail, pour L’utilisation combinée du SVC et
du TCSC en vue d’assurer de meilleures performances dynamiques du système avec plus de
rentabilité et d’efficacité.
La régulation de la tension et la puissance maximale qu’une ligne peut transporter sont deux
de ses plus importantes caractéristiques. En effet, la tension d’une ligne doit demeurer assez
constante à mesure que la puissance active consommée par la charge varie. Ordinairement, la
variation de la tension de zéro à pleine charge ne doit pas dépasser 5 % de la tension nominale,
bien qu’on puisse tolérer parfois une régulation allant jusqu’à 10 % .Afin de connaitre la
variation de la tension et établir la puissance maximale transportable par une ligne, nous
étudierons successivement le comportement de trois type de lignes : [30]
1. une ligne résistive
2. une ligne inductive
3. une ligne inductive relie à deux grands réseaux
Une ligne possédant une résistance R, par phase, alimente une charge résistive consommant
une puissance variable , la tension de la source est constante, on suppose une charge
résistive, on à :
R
(IV.2)
( )
(IV.3)
La puissance maximale :
(IV.4)
1 Es
0,95 Es
La tension ER (pu)
0,5 Es
0
0 19% Pc 100%
la puissance Pc (%)
La puissance active que la ligne peut transporter .on atteint ce maximum lorsque la résistance de
la charge est égale à celle de la ligne .Il s’ensuit que la puissance égale 100%
(puissance maximale).
Si l’on permet une régulation maximale de 5% ( ), la ligne peut transporter une
puissance qui représenté seulement 19% de la puissance maximale. La ligne pourrait
transporter une puissance plus grande que mais la tension correspondant serait alors trop
basse.
Considérons maintenant une ligne dont la résistance est négligeable, mais qui possède une
réactance inductive donc on à la figure suivant :
√
(IV.8)
√
(IV.9)
La puissance maximale :
(IV.10)
1 Es
0,95 Es
La tension ER (pu)
0,707 Es
0
0 60 % Pc 100 %
La puissance Pc (%)
Ligne inductive donne donc une meilleure régulation qu’une ligne résistive.
La source doit fournir non seulement la puissance active consommée par la charge, mais
aussi la puissance réactive absorbée par la ligne. [33][34]
La figure (IV.16) donnée le circuit équivalent d’une ligne inductive relient deux régions S et
R dont les tensions et aux deux extrémités sont constantes. En ce qui concerne l’échange
de puissance active, donc on aura la figure suivant:[1]
X X
I I
P P
FigureIV.16Uneligneinductiverelie à deuxgrandsréseaux
Ɵ (IV.11)
100 %
La puissance P (%)
50 %
0
0 30 90 180
L'angle (degrés)
{ (IV.12)
( )
Cette équation montre qu'il est possible d'augmenter la puissance active transitée entre deux
réseaux soit en maintenant la tension des systèmes, soit en augmentant l'angle de transport entre
les deux systèmes, soit en réduisant artificiellement l'impédance de la liaison.
En jouant sur un ou plusieurs de ces paramètres, les FACTS permettent un contrôle précis des
transits de puissance réactive, une optimisation des transits de puissance active sur les
installations existantes et une amélioration de la stabilité dynamique du réseau. Ils permettent
aussi aux consommateurs industriels de réduire les déséquilibres de charges et de contrôler les
fluctuations de tensions crées par des variations rapides de la demande de puissance réactive et
ainsi d'augmenter les productions, de réduire les coûts et d'allonger la durée de vie des
équipements.[7]
Le TCSC est le premier contrôleur série développé dans les systèmes FACTS. Ce type de
dispositifs a été installé pour la première fois dans la sous-station de Kayenta par (Western Area
Power Administration (WAPA)) au début de l’année 1992. [11]
L’organe essentiel de commande du TCSC est le TCR qui est un contrôleur statique
utilisant une interface d’électronique de puissance pour fournir un contrôle rapide de l’énergie
réactive. Le principal composant contrôlable est une réactance inductive connectée en série avec
un thyristor bipolaire de commande. [31]
La variation de l’angle d’amorçage des thyristors, permet de varier la valeur de la réactance
inductive, qui cause un échange rapide de l’énergie réactive entre le TCR et le système.
(IV.13)
(IV.14)
( )
i j
Th1
L
RCT
Th2
La structure d’un TCSC est ainsi équivalente à un SVC de type FC-TCR, la différence étant
que le TCSC est connecté en série avec la ligne de transmission, alors que le SVC est connecté
en shunt à un nœud préalablement sélectionné. Une autre différence, est que le TCSC est
connecté directement à la ligne alors que le SVC exige un transformateur de liaison.
Le principe de base du TCSC est l’ajustement optimal et rapide de l’impédance de la ligne à
laquelle il est connecté. En contrôlant la réactance du TCSC, l’impédance de la ligne est
compensée, et donc l’écoulement de puissanceà travers la ligne est augmenté.[8]
XC
i R+jX J
jB/2 XL jB/2
Th
TCR
La figure (IV.20) présente l’implantation d’un TCSC dans une ligne électrique. La réactance
totale Est la somme algébrique de la réactance du TCSC et la réactance réelle de la ligne. Voir
l’équation [34] :
(IV.15)
Pourévitertoutesurcompensation,lerangdefonctionnementduTCSCdoitêtrelimitéentre : [34]
(IV.16)
La figure suivante représentée le modèle de la réactance de l’implantation d’un TCSC :
( )
-
+
( )
+
Régulateur
PID
-
1,04
1,03
1,02
La tension (pu)
1,01
0,99
0,98
0,97
0,96
0 0.5 1 1.5
Temps (s)
1.0025
Pmés
1.002 Préf
La puissavce transmise (pu)
1.0015
1.001
1.0005
0.9995
0.999
0 0.5 1 1.5
temps (s)
D'après les résultats obtenus on peut conclure que les TCSC permettent de contrôler des
puissances actives dans un réseau. Le choix optimal de ces dispositifs est un élément très
important pour le contrôle de puissance.
Le compensateur statique de puissance réactive SVC est le premier dispositif FACTS qui
apparaît dans les années 1970 pour répondre à des besoins de stabilisation de tension rendue
fortement variable du fait de charges industrielles très fluctuantes telles les laminoirs et les fours
à arc. Les SVC sont des FACTS de première génération. Ils utilisent des thyristors classiques,
commandables uniquement à l’amorçage. [33]
Le compensateur statique de puissance réactive SVC est un dispositif qui sert à maintenir la
tension en régime permanent et en régime transitoire à l'intérieur de limites désirées. Le SVC
injecte de la puissance réactive dans la barre où il est branché de manière à satisfaire la demande
de puissance réactive de la charge. [8]
a)symbol b) modèle
e Figure IV.26 Modélisation du SVC
Les valeurs des SVC sont exprimées sous forme de puissance réactive absorbée à la
tension nominale Un. La correspondance avec la susceptance est donnée par la relation On
définit la susceptance de SVC :
(IV.18)
( ) ( )
(IV.19)
( ) (IV.20)
( ) ( )
Ces essais au simulateur consistaient à produire des perturbations à la barre où il est installé
le compensateur statique SVC, afin de vérifier s'il était en mesure de contrôler rapidement les
variations de tension à la suite de ces perturbations.
0.8
Bsvc(réference)
0.6 Bsvc(mésuré)
0.4
susceptance (Bsvc)
0.2
-0.2
-0.4
-0.6
-0.8
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
temps(s)
0.6
0.4
0.2
-0.2
-0.4
-0.6
-0.8
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Temps (s)
1.08
Vm(réference)
1.06
Vm(mésuré)
La tension Vm (pu)
1.04
1.02
0.98
0.96
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
temps (s)
Figure IV.29 Tension primaire réel et tension primaire mesuré contrôlé par le SVC
IV.9. Conclusion
Les dispositifs FACTS présentés dans ce chapitre possèdent tous leurs propres
caractéristiques tant en régime permanent qu’en régime transitoire. Chaque type de dispositif
sera donc utilisé pour répondre à des objectifs bien définies.les gestionnaires de réseau sont
entrain d’exploiter le système de transport aux prés de ses limites thermiques et dynamiques,
alors que les consommateurs sont de plus en plus exigeants quant à la qualité de l’énergie et à la
continuité de service. Les équipements à base d’électronique de puissance, y compris leurs
commandes appropriées, offrent des solutions efficaces à ce problème.
Conclusion générale
L’industrie de l’énergie électrique a subi des changements profonds, suite aux différents
bouleversements politiques et économiques de plusieurs pays dans le monde. Ce nouvel
environnement du marché de l’électricité, a mis en évidence une clientèle plus importante et plus
exigeante vis-à-vis du fonctionnement des réseaux à travers la flexibilité de leur contrôle, la
fiabilité et la qualité de l’énergie offerte sur le marché.
Ce mémoire présente et explique calcul d’écoulement de puissance dans les réseaux
électriques lors de l’intégration des fermes éoliennes et des dispositifs FACTS, la technique de
Calcul d’écoulement de puissance dans les réseaux électriques étudiée dans ce travail nécessite
en premier lieu l'étude des méthodes itératives classiques (Gauss, Gauss-Seidel).
Aussi représenté la modélisation et la commande de puissance d’un système éolien à
base de génératrice asynchrone à double alimentation alimentée par la commande vectorielle
directe, pour contrôler la puissance transitée au réseau.
Nous avons traité le problème de calcul de l’écoulement de puissances et le contrôle
des tensions des réseaux d’énergie électrique en incorporant les dispositifs SVC. Une étude sur
les dispositifs FACTS à été réalisé et en plus une étude détaillée sur les dispositifs SVC
et TCSC (définition, modélisation, simulation).
Afin d’atteindre cet objectif, nous somme passés par les étapes suivantes :
Dansle premier chapitrenous avons présenté en bref un aperçu sur le fonctionnement des
réseaux électriques, on a montré les facteurs qui limitent la puissance transmissible dans une
ligne de transport et la chute tension provoquée par le transit de puissance.
La deuxième chapitre représenté le problème de calcul d’écoulement de puissance dans un
réseau électrique. Pour résoudre ce problème dont nous avons choisi la méthode de GAUSS et
GAUSS-SEIDEL.
Le troisième chapitre nous avons présenté généralité sur l’énergie éolienne. Et les
différents types d’éoliennes, leurs caractéristiques technologiques et les techniques de stockage
associées, en suite la modélisation de la machine asynchrone à double alimentation MADA. Et la
simulation du fonctionnement de la machine en mode génératrice GADA. Ensuite la
modélisation la turbine d’éolienne et l’optimisation du fonctionnement de la partie mécanique de
l’éolienne et l’élaboration d’une commande adéquate pour le captage du maximum de puissance
possible, et contrôle de puissance active et réactive puis les résultats de simulation.
Le quatrième chapitre a été consacré à la description des différents systèmes permettant
d’assurer le transport d’énergie électrique dans les meilleures conditions par utilisé les systèmes
FACTS. Les dispositifs FACTS peuvent être classés en trois catégories : les compensateurs
parallèles, les compensateurs séries et les compensateurs hybrides (parallèle-série), et aussi
modélisation de compensateur statique SVC et la capacité série commande par
thyristor(TCSC). Dans ce chapitre, on étudie le fonctionnement du SVC et TCSC et on présente
les éléments qui constituent ce dispositif. Différentes configurations du SVC et TCSC sont
présentées ainsi que sa modélisation.et représenté les résultats de simulation de ce dispositif.
L’ensemble de ces travaux peut être poursuivi et complété par des perspectives pouvant
contribuer à l’amélioration de l’ensemble des chaînes de conversion d’éolienne, dans le but est
intégrer une énergie éolienne au réseau électrique, parmi les perspectives envisageables :
- Insertion des systèmes FACT pour améliorer la transit de puissance des systèmes
électrique.
- L’intégration éventuelle d’un système de stockage d’énergie éolienne pour le lissage de la
puissance fournie au réseau.
- Introduire des techniques intelligentes comme les réseaux de neurone, les algorithmes
génétiques.
- Utilisation des autres équipements FACTS dans l’amélioration de la stabilité statique de
la tension et la transit de puissance des réseaux électriques à savoir : UPFC, SSSC,
STATCOM,…etc.
Annexe A
Les paramètres de la génératrice asynchrone double alimentée:
La valeur nominale
Puissance nominale
Tension nominale
La vitessenominale
Nombre de pairs de pôles
Paramètre de la machine
Resistance statorique
Resistance rotorique
Inductancestatorique
Inductance rotorique
Inductance mutuelle
Vitessesynchronisme
Constantes mécanique :
Moment Inertie
Coefficient de frottement
Annexe B
30
100
Bibliographie
Références bibliographiques
[1]Abdelaali Alibi, «Contrôle des Réseaux Electriques par les Systèmes FACTS: (Flexible AC
Transmission Systems)»Mémoire Magister,Université de Batna, 2009.
[2] DelendiLouardi, «Contrôle de l’écoulement de puissance active par système FACTS », Thèse
de Magister, Université de Batna, Décembre 2009.
[3] A.Abdelbaki&D.Belkhiri& Y.Lahouassa,«Calcul de courant de court-circuit dans les réseaux
électriques », Mémoire d’ingéniorat de l’université de M’sila, Département d’Electronique,
2012.
[4] T.Siabdellah&O.Aouidji&Z.Doghmane ,«Protections des réseaux électriques MT et HT
»Thèse de d’ingéniorat ,Université de M’sila , 2012.
[5] Mr. Zellagui Mohamed «étude des protections des réseaux électrique MT (30 &10 kV)»
mémoire magistère, Université Mentouri Constantine ,2009.
[6]Slimane Ramdane, «Programmation en MATLAB de la Stabilité Transitoire d’une
Modélisation Interactive d’un Réseau Electrique par Intégration Numérique», Mémoire
Magister, Université Amar Telidji –Laghouat, Juin 2008.
[7] Haimour Rachida « Contrôle des Puissances Réactives et des Tensions par les Dispositifs
FACTS dans un Réseau Electrique », Mémoire Magister, université d’Oran ,2008/2009.
[8] Oussama Mammeri« Différentes méthodes de calcul de la puissance réactive Dans un nœud a
charge non linéaire en présence d'un systeme de compensation de l'énergie», Mémoire Magister,
l’université de Batna,2012.
[9]K. Belacheheb, "Contribution à l'étude des systèmes de compensation FACTS en général,
UPFC en particulier, pour le contrôle du transit de puissance dans un réseau de transport ", Thèse
de doctorat de l'université de Henri Poincare Nancy I, 21 Juin 2001.
[10]EskandarGholipourShahraki«Apport de l'UPFC à l'amélioration de la stabilité transitoire des
réseaux électriques» Thèse de doctorat de l'université de Henri Poincare Nancy I, 13 Octobre
2003.
[11] BoussadiaFethi« Amélioration de la stabilité D’un réseau électrique par facts application au
réseau algérien» Mémoire magister l’université de Sétif, 2008.
[12]Boudjella Houari "Contrôle des puissances réactives et des tensions dans un réseau de
transport au moyen de dispositifs FACTS (SVC) " Mémoire de magistère, Université de
sidi bel abbés 23/01/2008.
Bibliographie
Présenté par :
FERAHTIA Saàd
Thème :
2012/2013
N°d’order : 055