UNIVERSITÉ DU QUÉBEC EN ABITIBI-TÉMISCAMINGUE

MÉMOIRE PRÉSENTÉ À L'ÉCOLE DE GÉNIE

DE L'UNIVERSITÉ DU QUEBEC EN ABITIBI TÉMISCAMINGUE

COMME EXIGENCE PARTIELLE

DE LA MAÎTRISE EN INGÉNIERIE

PAR

KHALED OURAMDANE

MODÉLISATION ET OPTIMISATION DU SYSTÈME D'EXCITATI ON

DU GROUPE TURBINE-ALTERNATEUR DU SIMULATEUR ANALOGIQUE
D'HYDRO-QUEBEC

MAI2015
 Mise en garde

La bibliothèque du Cégep de l’Abitibi-Témiscamingue et de l’Université du Québec en Abitibi-

Témiscamingue a obtenu l’autorisation de l’auteur de ce document afin de diffuser, dans un but
non lucratif, une copie de son œuvre dans Depositum, site d’archives numériques, gratuit et
accessible à tous.

L’auteur conserve néanmoins ses droits de propriété intellectuelle, dont son droit d’auteur, sur
cette œuvre. Il est donc interdit de reproduire ou de publier en totalité ou en partie ce
document sans l’autorisation de l’auteur.

Warning
The library of the Cégep de l’Abitibi-Témiscamingue and the Université du Québec en Abitibi-
Témiscamingue obtained the permission of the author to use a copy of this document for non-
profit purposes in order to put it in the open archives Depositum, which is free and accessible to
all.

The author retains ownership of the copyright on this document. Neither the whole document,
nor substantial extracts from it, may be printed or otherwise reproduced without the author's
permission.
 Résumé

Le présent mémoire porte sur la modélisation et l'optimisation du système

d'excitation du groupe Turbine-Alternateur du Simulateur analogique d'Hydro-
Québec. L'important rôle que joue la partie excitatrice dans l'ensemble des
composantes du simulateur implique non seulement que cet élément assure la
protection des autres éléments du groupe (la surtension par exemple aux bornes de
l'alternateur), mais aussi sa propre protection (surintensité du circuit d'excitation par
exemple). Il intervient aussi dans l'amélioration des performances dynamiques du
simulateur tout en assurant sa fiabilité et sa durabilité.

Ce document présente en premier lieu : (1) une introduction générale, qui illustre
l'intérêt de ce projet d'un point de vue socio-économique et éducatif, les objectifs
visés ainsi que la méthodologie adoptée; (2) un état de l'art sur les systèmes
d'excitations des génératrices synchrones; (3) les modèles mathématiques de chaque
partie du groupe-turbine-alternateur sont développés, plus particulièrement le modèle
du système d'excitation, avant de faire une introduction aux régulateurs
conventionnel ( 4) souvent introduits au niveau de l'excitatrice; (5) une méthode
d' optimisation basée sur les algorithmes génétiques est présentée dans le but
d'optimiser les paramètres du système d'excitation.

Les modèles élaborés de chaque partie du groupe turbine-alternateur sont assemblés.

En prenant les valeurs des paramètres de la littérature courante, (6) une validation
sous matlab/simulink est effectuée, présentant ainsi les résultats de la régulation et
d'optimisation sur le système global.
 11

Abstract

This thesis focuses on the modeling and optimization of the excitation system of the
group turbine-altemator Hydro-Québec's digital Simulator. The important role plays
by the exciter part in the who le components of the simulator implies that this element
not only insures the protection of other elements of the group (overvoltage for
example across the terminais of altemator), but also its own protection (overcurrent of
the excitation circuit for example ). It also envolves in the improvement of the
dynamic performances of the simulator while assuring its reliability and durability.

This document first presents a general introduction (1), which illustrates the interest
of this project of a socioeconomic and educational point view, the objectives and
methodology used; (2) the state of art in excitation systems of synchronous generators
is presented; (3) the mathematical models of each part of the group slog altemator are
developed, particularly the excitation systems models, before making an introduction
for conventional regulators (4) often introduced at the exciter level. (5) The
optimization method based on the genetic algorithms is presented in order to optimize
the excitation system' s parameters.

Developed models of each part of the turbine-altemator group are assembled. The
validation of present work is made using matlab/simulink software and parameters
values obtained from the literature. Results of the control and optimization of the
global system's model are presented and discussed (6).
 111

Remerciements

Ce mémoire a été réalisé à l'Université du Québec en Abitibi-Témiscamingue

(UQAT) avec l'aide et le soutien de nombreux intervenants clés qui ont grandement
contribué à sa réalisation. Aussi, je tiens à les remercier et à leur exprimer toute ma
reconnaissance ;

En premier lieu, mon directeur de recherche, le Professeur René W amkeue (ing.

Ph.D), professeur titulaire à l'UQAT, pour m'avoir proposé ce sujet, pour son soutien
permanant et les nombreux conseils et informations apportés. Ce dernier a
grandement contribué à rehausser la qualité technique de ce travail.

M. Jean-Jacques Beaudoin, Professionnel, chargé des Laboratoires à l'UQAT,

pour sa disponibilité et son savoir-faire technique, mts en contribution pour la
compréhension de nombreux concepts pratiques.

Mes coéqmpters du groupe de recherche GREENG, pour leurs nombreuses

suggestions lors de nos présentations.

Je dédie particulièrement ce mémoire à ma sœur ainée Ouiza que Dieu veuille sur
son âme, à toute ma famille , à mon épouse et à Koceyla, et également à tous ceux qui
n'ont pas été cités. Je remercie le Seigneur de m'avoir permis d'accomplir ce projet.
 lV

Table des matières

Résumé .............. ........................ ........................ ....................... ........................ .............. i

Abstract ......................................................................................................................... ii

Remerciements ............................................................................................................. iii

Table des matières ........................................................................................................ iv

Liste des Figures ....................... ........................ ....................... ........................ ............. x

Liste des tableaux ...................................................................................................... xvii

Liste des symboles ................. ........................ ............................................... ........... xviii

C HAPITRE 1 INTRODUCTION GÉNÉRALE .................................................. 1

1. 1 Motivations et obj ectifs du travail... ................... ........................ ............. 2

1.1.1 Motivations socio-économiques .......................................................... 3

1.1.2 Objectif général .................... ....................... ........................ ................ 7
1.1.3 Obj ectifs spécifiques ........................................................................... 7

1.2 Méthodologie et organisation du proj et .................................................. 8

1. 3 Conclusion .......... ....................... ........................ ........................ ........... 10

CHAPITRE 2 ETAT DE L'A RT SUR LES SYSTÉMES D'EXCITATION .. 11

2.1 Introduction .............................. ....................... ........................ .............. 11

 v

2.2 Exigences d'un système d'excitation ................. ........................ ........... 11

2.3 Éléments d'un système d'excitation ..................................................... 13

2.4 Différents types de système d'excitation .............. ........................ ........ 15

2.4.1 Système d'excitation à courant continu (DC) ................................... 15

2.4.2 Systèmes d'excitation à courant alternatif (AC) ............................... 16
2.4.3 Systèmes d'excitation statique .......................................................... 18

2.5 Commande et circuits de protection des systèmes d'excitation ............ 20

2.5.1 Régulateurs AC et DC ....................................................................... 21

2.5.2 Circuits de st abilisation d'un système d'excitation ................... ........ 21
2.5.3 Stabilisat eur de puissance (PS S) ....................................................... 22
2. 5.4 Compensat eur de charge ...... ........................ ....................... .............. 23

2. 6 Conclusion ............................................................................................ 24

CHAPITRE 3 MODÉLISATION DES SOUS-ENSEMBLES DU GROUPE

TURBINE-ALTERNATEUR .................................................................................. 25

3.1 Introduction ......... ....................... ........................ ........................ ........... 25

3.2 Modélisation de la génér atrice synchrone ............................................. 26

3.2.1 Modèle de la génératrice synchrone dans le repère abc .................... 26

3.2. 2 Modèle de la génératrice synchrone dans le repère de Park (dqO) .... 28
3 .2.3 Développem ent du modèle hybride de la génér atrice synchrone ...... 34

3. 3 Modèle de la t urbine hydraulique et servom oteur ................................. 38

 Vl

3.4 Modélisation du système d'excitation statique .......................... ........... 40

3.4.1 Modèle par la méthode DQ ............................................................... 40

3.4.2 Modèle IEEE type ST. ...................................................................... 55

3. 5 Conclusion ............................................................................................ 62

CHAPITRE 4 ANALYSE DES RÉGULATEURS CONVENTIONNELS

D'UN GROUPE TURBINE-ALTERNATEUR ..................................................... 63

4.1 Introduction ........................................................................................... 63

4.2 Régulateur automatique de tension ....................................................... 63

4.2.1 Régulation et stabilisation de la tension d 'excitation ........................ 64

4.2.2 Régulation de la tension de sortie de la génér atrice synchrone : ....... 66

4. 3 Régulation de la vitesse de rotation ...................................................... 74

4.4 Stabilisateur de puissance ..................................................................... 76

4. 5 Conclusion .......... ....................... ........................ ........................ ........... 78

CHAPITRE 5 OPTIMISATION PAR ALGORITHME GÉNÉTIQUE ......... 79

5.1 Introduction .............................. ....................... ........................ .............. 79

5.2 Technique d' optimisation par algorithmes génétiques ......................... 80

5 .2.1 Principe de base d'un AG standard ................ ........................ ........... 81

 Vll

5.2.2 Codage des variables dans un AG .................. ........................ ........... 83

5.2.3 Opérateurs d'un algorithme génétique ........... ........................ ........... 83
5.2.3.1 Opérateur de sélection.......... ....................... ........................ ........... 84
5.2.3.2 Opérateur de croisement ................................................................. 84
5.2.3.3 Opérateur de Mutation .................................................................... 87

5.3 Optimisation mono objectif ........... ............................................... ........ 88

5.4 Conclusion ............................................................................................ 90

CHAPITRE 6 TESTS ET VALIDATION DE LA RÉGULATION (AVR,

PID) ET DE L'OPTIMISATION PAR ALGORITHME GÉNÉTIQUE DES
PARAMÉTRES D'EXCITATION DU GROUPE TURBINE-ALTERNATEUR
..................................................................................................................................... 91

6.1 Introduction ......... ....................... ........................ ........................ ........... 9 1

6.2 Paramètre des différents blocs .... ........................ ....................... ........... 92

6.3 Régulation de la tension d'excitation .............. ........................ .............. 94

6.4 Régulation de la tension de sortie de la génératrice .............................. 95

6.5 Régulation des t ensions dans le syst ème global.. ....................... ........... 97

6.6 Tests et simulations sur le systèm e global ............................................ 99

6.6.1 Test d'am orçage des tensions de sortie ............................................. 99

6.6.2 Tests -réponse à une variation de tension ........................... ............ 100

6.7 Régulation de la vitesse de rotation et de la puissance mécanique ..... 105

 V111

6. 8 Stabilisation de puissance ........... ....................... ........................ ......... 108

6. 9 Optimisation des paramètres du système d'excitation ........................ 112

6.9 .1 Optimisation des paramètres de 1'AVR ........................................... 114

6.9.2 Optimisation des paramètres du régulateur PID .............................. 115
6.9.3 Tests et Simulation sur le système global après optimisation ......... 116
6.9.3.1 Réponse à un échelon unitaire ..................................................... 117
6.9.3.2 Réponse à une variation de tension ............................................. 119

6.10 Conclusion .......... ....................... ........................ ........................ ......... 121

CONCLUSION GÉNÉRALE ................................................................................ 122

Références 125

Annexe 1 129

Annexe2 130

Annexe3 135

Annexe 4 135

Annexe 5 138

Annexe 6 140

Annexe 7 142

Annexe 8 142
 lX

Annexe 9 143

Annexe 10 144

Annexe 11 149

Annexe 12 150
 x

Liste des Figures

Figure 1-1 Schéma représentatif du simulateur analogique d'Hydro-Québec. 3

Figure 1-2 Production d'énergie aux États-Unis en 2013 - au Canada et au

Québec en 2012 ................................................................................ 4

Figure 1-3 Comparaison des prix moyens de l'électricité pour les clients
résidentiels des grandes villes en Amériques du Nord [4] ................ 5

Figure 1-4 Émission de gaz à effet de serre de différentes filières de

production [4] ................................................................................... 6

Figure 1-5 Schéma descriptif du projet ........................ ........................ ........... 10

Figure 2-1 Schéma fonctionnel du système d'excitation d'une génératrice

synchrone [ 1] ............................................................... 13

Figure 2-2 Système d'excitation de type DC ................ ........................ ........... 16

Figure 2-3 Excitation de type AC avec redresseur non commandé ................ 17

Figure 2-4 Excitation type AC avec redresseur commandé ........................... 17

Figure 2-5 Système d'excitation de type Brushless ........................................ 18

Figure 2-6 Excitatrice statique par redresseur commandé à source de tension

[1] .................................................................................................... 19
 Xl

Figure 2-7 Système d'excitation muni de circuits de protection .......... ........... 20

Figure 2-8 Schéma de stabilisation et régulation de 1' excitation [ 1] [2] ........ 22

Figure 2-9 Schéma du compensateur de charge [1] ....... ..... ........ . .... . ...... 23

Figure 3-1 Enroulements d'une machine synchrone triphasée (repère abc) . . 28

Figure 3-2 Schéma de la machine synchrone en diphasé (référentiel dqO) .... 29

Figure 3-3 Circuits électriques équivalents de la génératrice synchrone

triphasée dans le repère de Park (en pu) ......................................... 33

Figure 3-4 Illustration des entrées-sorties du modèle hybride de la génératrice

....................................................................................................... 34

Figure 3-5 Schéma fonctionnel des états (flux) .............................................. 36

Figure 3-6 Schéma fonctionnel des sorties ( l.f/s, ~ ) ....................................... 37

Figure 3-7 Schéma du convertisseur commandé lié au réseau ............ ........... 40

Figure 3-8 Circuit équivalent du redresseur commandé dans le repère dq .... 4 1

Figure 3-9 Circuit électrique équivalent du redresseur commandé lié au réseau

dans le repère dq [1 9]. ..................................................................... 43

Figure 3-10 Convertisseur avec régulateurs ................ ........................ .............. 46

Figure 3-11 Régulation par PI ........................................................................... 46

 Xll

Figure 3-12 Circuit électrique équivalent du convertisseur régulé [20] ............ 47

Figure 3-13 Tension de sortie Vj(a = 10°) ........................................................ 50

Figure 3-14 Tension de sortie Vj (a= 15°) ........................................................ 51

Figure 3-15 Calcul des paramètres Kpi et Kii ...................................................... 52

Figure 3-16 Calcul des paramètres Kpv et Kiv·················· ························ ··········· 52

Figure 3-17 Tension de sortie Vj régulée .......................................................... 53

Figure 3-18 Modèle de 1'amplificateur ............................................................. 56

Figure 3-19 Système d'excitation et transformateur de stabilisation [1]. .......... 56

Figure 3-20 Fonctionnement des limitations .................. ........................ ........... 58

Figure 3-21 Fonctions de déclenchement (Gate) ............ ........................ ........... 60

Figure 3-22 Structure du modèle détaillé du système d'excitation .................. 60

Figure 3-23 Modèle IEEE du système d'excitation statique type STlA ......... 6 1

Figure 4-1 Principe de la régulation de tension dans un groupe turbine-

alternateur [1] [25] [26] .................................................................. 64

Figure 4-2 Système d'excitation type STl de IEEE ....................................... 65

Figure 4-3 Régulation et stabilisation de la tension d'excitation ........ ........... 65

 X111

Figure 4- 4 Régulation de la tension de sortie de la génératrice ............ ........... 66

Figure 4-5 Schéma fonctionnel d'un correcteur PID classique ...................... 66

Figure 4-6 Schéma de régulation de tension simplifié ................................... 67

Figure 4-7 Principe de la régulation de la vitesse ........................................... 74

Figure 4-8 Régulateur de vitesse type PID .................. ........................ ........... 75

Figure 4-9 Schéma bloc du Servomoteur ....................................................... 75

Figure 4-10 Système d'excitation type ST1 avec PSS ...................................... 76

Figure 4-11 Stabilisateur de puissance de type IEEE PSS1A .......................... 77

Figure 4-12 Schéma de la régulation globale ........ ............. ..... .............. 77

Figure 5-1 Organigramme d'un AG standard.............. ........................ ........... 82

Figure 5-2 Croisement en 1 point (codage binaire) ........................................ 86

Figure 5-3 Croisement en 2 points (codage binaire) ...................................... 86

Figure 5-4 Croisement uniforme (codage binaire) ......................................... 86

Figure 5-5 Représentation d'une mutation de bit dans une chaîne (codage
binaire) ............ ..... ..... ............. ..... ............. ..... ............ 88

Figure 6-1 Régulation et stabilisation de l'excitation .... ................................. 94

 XlV

Figure 6-2 Tension d'excitation VFn .... ........................ ....................... ........... 94

Figure 6-3 Régulation de la tension génératrice V8 ........................................ 95

Figure 6-4 Tension génératrice Vs .................................................................. 95

Figure 6-5 Tension d'excitation correspondante VFD····································· 96

Figure 6-6 Représentation schématique du système global (Génératrice +

Excitation) ....................................................................................... 97

Figure 6-7 Tension de sortie de la génératrice Vsabc- ...................................... 98

Figure 6-8 Tension d'excitation Vpn correspondante .... ....................... ........... 98

Figure 6-9 Tension de sortie de la génératrice Vsahc- ...................................... 99

Figure 6-10 Tension d'excitation Vpn correspondante ...................................... 99

Figure 6-11 Puissances au rotor. ............. ....................... ........................ ......... 100

Figure 6-12 Tension d' excitation Vpn et tension de sortie Vsabc ············· ········· 101

Figure 6-13 Puissances au niveau du rotor. ................. ........................ ............ 101

Figure 6-14 Puissances au niveau stator. ......................................................... 102

Figure 6-15 Tension d'excitation VFn et tension de sortie VSabc ············· ········· 102

Figure 6-16 Puissances au niveau du rotor. ..................................................... 103

 xv

Figure 6-17 Puissances au niveau du stator.. .................. ........................ ......... 103

Figure 6-18 Régulateur de vitesse type PID .................................................... 105

Figure 6-19 Vitesse de rotation de la turbine mret puissance mécanique Pm·· 106

Figure 6-20 Ouverture de vannage g et débit d'eau q associés .............. ......... 107

Figure 6-21 L'écart de vitesse co- cor•f ........................................................... 107

Figure 6-22 Schéma de la régulation globale .................................................. 108

Figure 6-23 Variation de la tension d'excitation VFD avec et sans PSS .......... 109

Figure 6-24 Variation des puissances au niveau du rotor avec et sans PSS .... 109

Figure 6-25 Tension VFD······················· ······················· ························ ············ 110

Figure 6-26 Puissances au rotor ....................................................................... 110

Figure 6-27 Tension de stabilisation Vstah ...................... ........................ ......... 111

Figure 6-28 Tension d' excitation VFD- avant et après optimisation ............... 114

Figure 6-29 Tension Vs - avant et après optimisation.... ................................. 115

Figure 6-30 Tension génératrice Vsahc (après optimisation) ............................ 117

Figure 6-31 Tension d'excitation VFD associée (après optimisation) ..... ......... 117

Figure 6-32 Variation des puissances au rotor ............... ........................ ......... 118
 XVl

Figure 6-33 Puissances au stator ............. ....................... ........................ ......... 118

Figure 6-34 Tension d'excitation VFD et Tension génératrice Vsabc ................ 119

Figure 6-3 5 Variation des puissances au rotor ................................................ 120

Figure 6-36 Variation des puissances au stator ..................................... 120

 XVll

Liste des tableaux

Tableau 1: Paramètres de la génératrice synchrone Labvolt ...................... ........... 92

Tableau 2: Paramètres du système d'excitation statique STl ............................... 93

Tableau 3 : Paramètres de la turbine hydraulique ................................................. 93

Tableau 4: paramètres du régulateur automatique de tension AVR ................... 114

Tableau 5 : paramètres du régulateur PID .... ....................... ....................... ......... 115

Tableau 6 :paramétrés optimaux du système d'excitation STl.. ........................ 116

 XV111

Liste des symboles

ST Statique

AC Courant alternatif

DC Courant continu

AVR Régulateur Automatique de tension

PSS Stabilisateur de puissance

KJ Gain de compensation de 1' excitatrice ST

Tf Constante de temps du compensateur

KA Gain de l'amplificateur de l'excitatrice ST

TA Constante de temps de l'amplificateur

Va, Vb, Vc V/pu Tensions au stator repère abc

Vd,Vq V/pu Tensions des amortisseurs d et q

Vj V/pu Tension de champ

ia, ib, ic A/pu Courants au stator repère abc

id, iq A/pu Courants axe d et q

iJ A/pu Courant de champ

lf'a' lf'b' lf'c Wb/pu Flux au stator repère abc

lf/d ' lf/q

Wb/pu Flux axe d et q

If/! Wb/pu Flux de champ

 Xl X

P(B) Matrice de Park

1
F (B) Matrice de Park inverse

.O./pu Matrice résistance au stator

.O./pu Matrice résistance au rotor

.o./pu Résistances d'armature axe d et q

.o./pu Résistances de champ, d'amortisseurs D et Q

abc
Vs 'Vs
dqO
v Vecteurs tension au stator

abc dqO
Vr 'Vr v Vecteurs de tension au rotor

.abc .dq
ls , l s A Vecteurs courant au stator

. abc . dqO
lr ' lr A Vecteurs courant au rotor

Il Fabc Il Fdq O wb Vecteurs flux au stator

't" s ''t" s

, ,abc , • .dq O
lYr , tpr wb Vect eurs flux au r otor

rad/s Vitesse de synchronism e

e rad Angle

H Matrices des inductances mutuelles

pu Matrices des réact ances mutuelles

pu Réactances propres d' axe d et q

X md, Xmq pu Réactances de magnétisation d 'axe d et q

Xkjl pu Réactance de kany

 xx

Xj pu Réactance de fuite

Vsref v Tension de base

fsref A Courant de base

Srej VA Puissance de base

Cùrej radis Vitesse de référence

If/rej wb Flux de référence

Zef Q Impédance de base

H pu Constante d'inertie

TM pu Couple mécanique

D Coefficient d'amortissement

Te pu Couple électrique

Vs, Vr pu Vecteur tension au stator et rotor

fs, fr pu Vecteur courant au stator et rotor

ll's' ll'r pu Vecteur flux au stator et rotor

Cùn pu Pulsation propre

g Vannage

q Débit d' eau

Pm Puissance mécanique

A, Gain de la turbine
 XXl

Tw s Temps d'établissement du débit nominal

Ta Constante de temps du servomoteur

Ka Gain du servomoteur

u Commande du servomoteur

CPL Charge idéale

PcPL Kw Puissance à travers CPL

Vsa, Vsb, Vsc v Tensions de source triphasée en repère abc

Req• Leq• Ceq Q,H,F hnpédances équivalentes à 1' entrée

rF, LF, CF Q,H,F hnpédances liées au filtre à la sortie

IcPL A Courant à travers la charge CPL

Vj v Tension aux bornes de la charge CPL

Vbus,a, Vbus,b, Vbus,c v Tensions de bus d'entrée en repère abc

Vbus,d, Vbus,q v Tensions de bus d'entrée en repère dq

Edc, EdcJ v Tensions de sortie du convertisseur

la,h)c A Courants d'entrée en repère abc

IJ, Iq A Courants d'entrée en repère dq

Ide A Courant à la sortie du convertisseur

À rad Déphasage

Sabc Matrice de passage dans le repère abc

 XXll

T
Sabc Matrice de passage inverse, repère abc

a rad Angle d'amorçage des thyristors

~ rad Angle de phase des tensions de bus

~1 rad Angle de phase dans le repère dq

T[B(t)] Matrice de transformation

Sdq Matrice de passage dans le repère dq

Sd/ Matrice de passage inverse dans le repère dq

Vsd 'Vsq v Tension de source dans le repère dq

fdsJqs A Courants de source dans le repère dq

Vj* v Consigne de la tension de sortie

* A Consigne du courant de sortie

Ide

Kc Facteur de rectification de tension

VFD pu Tension d 'excitation

KE pu Gain de 1' excitatrice

TE Constante de temps de 1' excitatrice

KLR Facteur de rectification du courant

VAmax' VAmin pu Limitation de tension

Tc,Tcl ,TB,TBI Constante de temps du compensateur

Vstab pu Tension de stabilisation

 XX111

VREF pu Tension génératrice de référence

Vc pu Tension corrigée

/pD pu Courant d'excitation

PID Régulateur classique

Kp Gain proportionnel du PID

KI Gain intégral du PID

KD Gain dérivée du PID

Rp Statisme

w pu Vitesse de rotation

Wref pu Vitesse de rotation de référence

AG Algorithme Génétique

f Fonction objective

Pc Probabilité de croisement

Pm Probabilité de mutation

F Fonction d'adaptation
 1

CHAPITRE 1

INTRODUCTION GÉNÉRALE

La plus part des centrales hydroélectriques au Québec innovent et œuvrent afin

d'assurer l'intérêt de tous en matière de production d'énergie électrique. A cet effet,
le projet de construction d'un simulateur de centrale hydroélectrique groupe turbine-
alternateur a été lancé chez Hydro-Québec en collaboration avec l'École de Génie de
l'Université du Québec en Abitibi-Témiscamingue, dans le but de former du
personnel technique destiné à œuvrer dans ce secteur. Le projet intitulé :
« Modélisation et Optimisation du système d'excitation du groupe turbine-alternateur
du simulateur d'Hydro-Québec »vise l'étude d'un élément important et fondamental
dans le fonctionnement du groupe turbine-alternateur qui est essentiellement le
système d 'excitation ou excitatrice.

Cependant, le système d'excitation est élaboré en accord avec les exigences des
cahiers de charges, et pour se faire, on utilise des correcteurs en vue de contrôler les
tensions de sortie des machines. En effet, l 'amplitude de ces tensions dépend
intégralement du champ tournant principal; l'idée est donc, en fonction de 1"etat de la
machine synchrone utilisée, d'adapter le champ inducteur pour garder une amplitude
constante en sortie. Pour cela, la principale contrainte est généralement son immunité
vis-à-vis des perturbations extérieures (charge, délestage) et intérieures (variations
des paramètres).
 2

En outre, l'utilisation de l'excitatrice au niveau de l'inducteur présente quelques

exigences par rapport au fonctionnement de la machine synchrone ainsi que sur le
réseau électrique, à savoir :

La fiabilité de l'excitatrice et sa capacité à amortir les oscillations de puissance

au niveau de la génératrice [ 1].
L'amélioration de la qualité des s1gnaux, qm présente souvent des taux de
distorsion assez importants.
L'assurance des limitations nécessaires en cas de surcharges; le courant
d'excitation de 1' inducteur principal peut être limité par la tension de sortie de la
machine via les circuits limiteurs de l'excitatrice [1] [2].

1.1 Motivations et objectifs du travail

En plus des apports pédagogiques que le simulateur d'Hydro-Québec (Figure 1-1)

peut offrir en t ermes d 'éducation et d 'apprentissage dans le domaine de génie et plus
précisément dans la production d 'énergie électrique en centrales hydroélectriques,
une formation beaucoup plus facile et prat ique aux étudiants en génie. La mise en
place du simulateur est également motivée par des enjeux socio-économiques assez
importants.
 3

Réservoir d 'eau supérieur

Conduit e forcée

Réservoir d'eau

Figure 1-1 Schéma représentatif du simulateur analogique d 'Hydro-Québec

1.1.1 Motivations socio-économiques

Le Québec est la terre de l'hydroélectricité. Grâce à sa géomorphologie et à ses

ressources hydroélectriques, le Québec a accès à une forme d'énergie propre,
renouvelable, abondante et constante.

Avec plus de 130 000 cours d'eau et 500 000 lacs, qui représentent au total 12%
de son territoire, le Québec est l'une des régions les plus riches en eau du monde, et
contient plus de 40% des ressources hydrauliques du Canada [3]. Cette spécificité
 4

constitue un réel atout dont le gouvernement québécois a su tirer profit, en faisant du

développement de la filière hydroélectrique une priorité de sa stratégie énergétique.
L'aménagement de la Baie-James en est l'exemple le plus illustre, puisque celle-ci
abrite aujourd'hui le plus grand complexe hydroélectrique au monde.

96 % de l'électricité produite au Québec provient de l'énerg ie hydroélectrique

(Figure 1-1) [3] [4]. L'entreprise d'État, Hydro-Québec, g ère d'immenses ensembles
de bmTages, qu'elle peut ouvrù· ou fe1mer instantanément pour ,Yuster l'offi·e à la
demande. L'élect1icité est donc abondante, propre et surtout constante, contrairement
aux autres énergies renouvelables généralement in te1m itten tes.

États-Unis 2013l•l Canada 20121bl Québec 20121betc)

Production d'électricité: 4 060 TWh Production d'électricité : 615 TWh Production d'électricité : 200 TWh

Hydroufique Nucléaire
6,5% 2%
\ Autres
Autres , Éolien 1%
Nucléaire 1% 1 4% Èolien
Mazout
19 % 1 1% 1%
Biomasse
Mazout 1,5 % Gaz
1% naturel
9%

Charbon_
10%
Hydraulique
Biomasse 61 %
Gaz 2% 1
naturel Charbon 1 •
28% l Eolien Hydraulique
39%
2% 96%
Énergies renouvelables • 13 % Énergies renouvelables • 65% Énergies renouvelables > 97 %

Figtue 1-2 Production d'énergie aux États-Unis en 2013- au Canada et au

Québec en 2012

L'électricité est tellement bon marché au Québec que 70% des foyers se chauffent
par ce moyen en hiver. Le prix de l'électricité est en outre le même partout, par
 5

péréquation (Figure 1-2), que l'on se situe dans la communauté inuit du nord ou en
plein centre de Montréal, le tarif pratiqué est plat; ni dégressif, ni progressif.

C~>nts par kWh

30

25

20

15

10

0
~
0
>-
~
8
.!!
~
f
c:
~
~
iii
:t:
~
Cl
tl
sc:
e0
...
"'~
B
0 u
~
~
"'
u~

~
•
9
Uj
~
~
0
~
i
~
ï:
c:
li
'Ec:
0
z IL
c: :: ~ :1

"'
Cl)

Not" : Pour une consommation mensuelle de 1 000 kWh, tarifs en vigueur le 1"' avril 20 14

Figure 1-3 Comparaison des prix moyens de l' électricité pour les clients
résidentiels des grandes villes en Amériques du Nord [4]

Grâce à ses ressources hydrauliques, le Québec bénéficie donc d'un actif non
négligeable, qui prend une valeur additionnelle dans un contexte de lutte contre les
gaz à effet de serre (Figure 1-3). Avec cette prépondérance de l'hydroélectricité dans
son mix énergétique, le Québec produit deux foi s moins de gaz à effet de serre par
habitant que l'ensemble du Canada.
 6

, ... ... ... --

---... ... ..
/
\
Hydraulique
au fil de l'eau ,
'14
........... .. , '
-------
Nucléaire 1 6

Éolien 19
',
{ Hyd~auliqu~
', ... _
avec reservotr,' 10
, '
........ _ .; ...
Solaire
photovoltaïque 38

Thermique:
gaz naturel 422
~------------------~
Thermique:
charbon 957

0 200 400 600 800 1000

(g éq. C021kWhl

Hydro-Ouébec

Figure 1-4 Émission de gaz à effet de serre de différentes filières de production [4]

Par conséquent, l ' immense importance économique qu'occupe l'hydroélectricité

au Québec, encourage d'avantage le principal fournisseur en électricité, Hydra-
Québec à œuvrer pour l'établissement du simulateur analogique, qui certainement,
contribuera non seulement dans le développement du secteur hydroélectrique
provincial, mais aussi assurera le relais entre universitaire et industriel en termes de
progrès et d'innovation (sachant que plusieurs essais et tests pourront être effectués
sur le simulateur analogique indépendamment du réseaux électrique).
 7

1.1.2 Objectif général

Le but principal visé, étant d'assurer une meilleure qualité et fiabilité des
alimentations électriques tout en garantissant une solide protection des équipements
électriques; soit de la génératrice synchrone, la turbine hydraulique, et éventuellement
du réseau électrique au niveau du simulateur d'Hydro-Québec. Cela en agissant
seulement sur la partie excitation du groupe.

1.1.3 Objectifs spécifiques

Les objectifs qui sont spécifiquement visés, dans cette étude peuvent être illustrés
comme suit:

Élimination des perturbations extérieures (charges, délestages, ... etc) via des
contrôleurs ou régulateurs bien spécifiques introduits au niveau de l 'excitatrice.

Optimisation des paramètres du système d'excitation du groupe turbine-

alternateur du simulateur d'Hydro-Québec.

Assurance d'une meilleure stabilité, rapidité et précision au mveau de la

tension génératrice en optimisant les différents paramètres de 1' excitatrice par la
méthode d'optimisation par algorithme génétique.

Amélioration de l' amortissement des oscillations de pmssance de la

génératrice.
 8

1.2 Méthodologie et organisation du projet

Notre travail est structuré suivant quatre parties principales élaborées selon le
cheminement ci -dessous :

• Revue de littérature :

Un État de l'art sur les systèmes d'excitation des génératrices synchrones sera
présenté. Les différents types de systèmes d'excitations seront alors abordés, leurs
utilisations, leurs fonctionnements et le rôle important qu'ils jouent dans un groupe
turbine-alternateur (chapitre 2).

• Modélisation des sous ensembles du simulateur :

Avant de procéder aux essais et différents tests sur 1' ensemble du simulateur
analogique, groupe turbine-alternateur, la modélisation de chaque bloc est
indispensable. Le modèle de la génératrice est alors présenté, puis adapté à la
problématique du thème, à son tour la partie hydraulique du simulateur définie
principalement par la turbine hydraulique est donnée sous le modèle très fréquent
dans la littérature [5]. Et enfin, la modélisation du système d 'excitation, qui présente
la partie la plus importante dans notre travail, avant de faire l'assemblage et la mise
en commun des différents modèles, dans un seul modèle global, cela permettra alors
d'effectuer des essais et des simulations sur l'ensemble du système, ce qui conclut le
chapitre 3.
 9

• Régulation et contrôle des grandeurs de sortie :

Cette étape consiste, en la régulation via des contrôleurs bien choisis (PI, PD,
PID, A VR, ... etc) qui seront introduits dans la partie excitatrice du système global
définissant le groupe turbine alternateur, afin d'améliorer les performances
dynamiques de l'ensemble du système. Le fonctionnement et le rôle de chaque
régulateur est alors illustré et présenté dans le chapitre 4.

• Optimisation des paramètres du système d'excitation:

Une méthode d'optimisation appropriée sera choisie (algorithme génétique) et

appliquée, afin de générer les résultats les plus optimaux de chaque paramètre du
système d 'excitation du simulateur analogique du groupe turbine-alternateur. La
présentation de cette méthode d 'optimisation est définie dans le chapitre 5.

Notre travail alors peut être représenté par le schéma descriptif de la Figure 1-5
montrée ci-dessous où sont résumées, les parties essentielles du projet.

La validation et les résultats obtenus sont alors présentés dans le dernier chapitre
(6) avant de finir par une conclusion générale.
 10

Puissance
électrique

-Adaptation d ' un algoritbme -Protection (Circuits limiteurs).

d'optimisation approprié (A.G). - Compensateur de cbar&e
-Assurance des conditions - Régulation de tension.
optimale - Stabilisateur de puissance

Puissance
mécanique

Figure 1-5 Schéma descriptif du projet

1.3 Conclusion

Ce chapitre nous a donné un aperçu introductif de notre travail de recherche, le

situant dans le contexte socio-économique et éducatif, les motivations de ce projet et
les objectifs visés ont été énumérés ainsi que la méthode adoptée dans la résolution de
chaque pattie du thème.
 11

CHAPITRE2

ETAT DE L'ART SUR LES SYSTÉMES D'EXCITATION

2.1 Introduction

Ce chapitre est consacré à la partie la plus importante dans notre travail, qui est
définie par le système d'excitation du groupe turbine-alternateur. Avant d'aborder les
éléments essentiels qui constituent un système d'excitation, nous allons discuter les
exigences et nécessités pour le fonctionnement des systèmes d'excitations. Les
différents types de systèmes d'excitations seront par la suite présentés suivant leurs
utilisations. Le rôle et l'importance du système d'excitation dans un groupe turbine-
alternateur est en principe la possibilité d'effectuer des opérations de commande, de
régulation et de stabilisation dans le but d 'assurer la protection des éléments dans
chaque bloc du groupe. Les différents circuits de commande et de protection seront
alors détaillés dans cette partie montrant alors 1'utilité et la nécessité et le
fonctionnement de chaque élément d'un système d'excitation; excitateur, régulateurs,
stabilisateur, compensateur, ... etc.

2.2 Exigences d'un système d'excitation

Les extgences de performances d'un système d'excitation sont principalement

déterminées par la génératrice synchrone et aussi le circuit d' alimentation. La
 12

condition de base est que le système d'excitation doit ajuster le courant de champ afin
d'exciter la génératrice synchrone et maintenir la tension de sortie selon la consigne.

La fonction de base de l'excitation est alors de fournir le courant continu aux

bobines créant un champ au rotor de la machine synchrone; de plus, l'excitation
assure le contrôle de la tension et de la puissance réactive [1] [6] [7], ainsi que
certaines fonctions protectrices du réseau (respect des limites des machines), le tout
par le biais de la tension d'excitation.

D'un point de vue d'installation électrique, le système d'excitation devrait

contribuer au contrôle efficace de la tension et à l'amélioration de la stabilité du
système [1] [8]. Cela implique aussi le maintien de la stabilité transitoire en cas de
perturbations et de faibles signaux.

Historiquement, le rôle du système d'excitation dans l 'amélioration des

performances du système de production n'a pas cessé de se développer. Auparavant,
les systèmes d'excitation intervenaient manuellement dans la régulation et le maintien
constant de la tension et la puissance réactive en sortie des génératrices. Quant le
contrôle de la tension génératrice a été automatisé la première fois, il était lent. Au
début des années 1920, l'intégration des régulateurs continus dans les systèmes
d' excitations a contribué dans le développement de ceux-ci, permis d' obtenir de
meilleurs résultats, et suscité l'intérêt d'un grand nombre d'industriels. Depuis ce
moment, le système d' excitation a connu une évolution continue. À partir des années
1960, l'importance et le rôle du système d'excitation ont été élargis en utilisant des
signaux de stabilisation auxiliaires. En plus de la régulation de tension de sortie,
l' excitation permettait d'amortir les oscillations de puissance du système à l' aide de
stabilisateur de puissance (power system stabiliz er) [9] [10]. Les systèmes
d'excitation modernes sont actuellement capables de fournir une réponse
pratiquement instantanée même avec des tensions élevées (HT) [ 1] [ 11].
 13

Pour remplir son rôle de manière satisfaisante, le système d'excitation doit

répondre aux exigences suivantes [ 1] [ 12] :
Respecter les critères souhaités du cahier de charge.
Fournir les limitations nécessaires afin d' assurer la protection de l'excitatrice
elle même, de la génératrice, et des autres équipements (turbine, réseau, ... etc).
Respecter les exigences tout en assurant une flexibilité d'opération
(régulation, stabilisation, ... etc)
Assurer la fiabilité et la validité des équipements.

2.3 Éléments d'un système d'excitation

Le schém a ci-dessous (Figure 2-1 ) représente les différentes parties et les

éléments d'un système d'excitation appliqué pour les génératrices synchrones.

®
Limite urs et c ircuits ..,_
prot ecteurs

®
Senseur de t ension et
compensat eur de c harge ~

0 CD
Réf
~
~
~
Régulat eur ____, Excitat rice Génératrice r ___.. Syst ème
d'ali mentati on

Stabilisateur de puissance ~

Figure 2-1 Schéma fonctionnel du système d'excitation d'une génér atrice

synchrone [ 1]
 14

- (1) L'excitatrice: fournit la puissance aux bobines qui créent le champ au rotor.
Typiquement, un générateur d'une puissance de 600MW nécessite une excitatrice
d'environ 1MW.

- (2) Le régulateur : permet de contrôler et stabiliser le signal d'entrée de

1' excitatrice.

- (3) Le senseur de tension : mesure la tension aux bornes de la machine, la filtre

et la redresse en une valeur continue. Ensuite, cette mesure est comparée à la
référence (tension voulue aux bornes de la machine).

- Le compensateur de charge : est un élément additionnel et facultatif qui permet

de conserver une tension constante en un point particulier distant des bornes de la
machine.

- (4) Le stabilisateur de puissance (PSS) : il fournit un signal d'entrée additionnel

au régulateur du système d'excitation afin d'amortir les oscillations de puissance à
basse fréquence du réseau électrique.

- (5) Les limiteurs et les circuits protecteurs : sont toute une gamme de moyens
qui s'assurent que les limites du circuit d'excitation et du générateur ne sont pas
outrepassées. Par exemple, limitation du courant d'excitation, limiteur de
surexcitation, limiteur de tension et limiteur de sous-excitation.
 15

2.4 Différents types de système d'excitation

Dépendamment de la source d'alimentation, les systèmes d'excitations pour les

génératrices synchrones se distinguent selon trois types essentiels [1] [13]:
./ Systèmes d'excitation à courant continu (DC)
./ Systèmes d'excitation à courant alternatif (AC)
./ Systèmes d'excitation statiques

Dans cette section, ces différents types d'excitation seront décrits, suivant la
forme et la structure générale de chaque.

2.4.1 Système d'excitation à courant continu (DC)

Cette catégorie d'excitation utilise une machine à courant continu comme source
d'excitation, fournissant le courant au niveau du rotor de la génératrice synchrone à
travers des bagues collectrices. La machine à courant continue est entrainée par un
moteur de l'arbre de la génératrice, et peut être soit auto-excitée, soit excitée
séparément.

L'excitation à courant continu est la première utilisée de l'année 1920 jusqu 'en
1960 où elle a été remplacé par l'excitation de type AC [1].

Actuellement, ce type d'excitation est rarement utilisé pour fin d'étude de la

stabilité car il exige une modélisation complexe. Il est souvent remplacé par
l'excitatrice de type AC ou l'excitatrice statique ST [1].
 16

Génératrice
Excitation OC synchrone
Bague collectrice
Ampli

Rhéostat

L _____ -------

Régulateur
de tension

Figure 2-2 Système d'excitation de type DC

2.4.2 Systèmes d'excitation à courant alternatif (AC)

L' excitatrice à courant alternatif utilise un alternateur associé a un redresseur

stationnaire ou tournant pour produire le courant de champ excitant la machine
synchrone, le redresseur utilisé peut aussi être soit commandé ou non. Dans le cas
d'un redresseur non commandé, le réglage s'effectue sur les enroulements de champs
au niveau de l'excitatrice AC.

Les figures 2-3 et 2-4 montrent l'excitatrice de type AC muni d'un redresseur
stationnaire non commandé et commandé respectivement.
 17

Génératrice
Excitation AC Redresseur synchrone
non commandé Bague collectrice

commandé
Régulateur
oc
Référence AC
Régulateur
AC Entrées

Figure 2-3 Excitation de type AC avec redresseur non commandé

Génératrice
Excitation AC Redresseur
synchrone
commandé Bague collectrice

]G) cTT
Régulateur
oc
Référence AC
Régulateur
AC Entrées

Figure 2-4 Excitation type AC avec redresseur commandé

 18

Dans le cas de 1'excitation à courant alternatif avec redresseur tournant, les bagues
collectrices sont éliminées et la sortie du redresseur à diodes est directement reliée à
la machine synchrone. La Figure 2-5 illustre l'excitation AC de type Brushless où
1' on observe la structure tournante composée de 1' excitatrice principale AC et du
redresseur à diodes, tel que la sortie du redresseur de 1' excitatrice pilote active le
champ à la structure tournante [1] [13].

Génératrice
Excitatrice Structure tournante synchrone
pilote

Triphasé
AC

Régulateur Control manuel

Entrées

Figure 2-5 Système d'excitation de type Brushless

2.4.3 Systèmes d'excitation statique

Dans ce type de système d'excitation, tous les composants sont statiques ou

stationnaires. Le redresseur statique fournit le courant d'excitation de la génératrice
synchrone à travers des bagues collectrices, et celui-ci est alimenté par la génératrice
synchrone à travers un transformateur d'excitation. Cependant il existe trois types
d'excitation statique :
 19

(a) avec redresseur commandé à source de tension

(b) avec redresseur non commandé à source combinée
( c) avec redresseur commandé à source combinée

Génératrice
Transformateur Redresseur
synchrone
d'excitation commandé Bague collectrice

Régulateur
oc
Référence AC
Régulateur
AC Entrées

Figure 2-6 Excitatrice statique par redresseur commandé à source de

tension [1]

L' excitatrice statique par redresseur commandé à source de tension (Figure 2-6)
est la plus utilisée des trois, voici énumérés quelques unes de ses caractéristiques :
./ Pas de partie tournante à part la machine .
./ Utilisation d 'un générateur de tension aux bornes, d'un autotransformateur
d' excitation et d'un pont de thyristor pour le redressement de la tension AC d'entrée .
./ Possibilité d'utilisation de plusieurs types de correcteurs et régulateurs (PI,
PID, Flous, ... etc) [14].
./ Moins coûteuse et facile à entretenir
./ Plus adaptée pour les petites centrales hydroélectriques
 20

2.5 Commande et circuits de protection des systèmes d'excitation

La commande d'un système d 'excitation ne dépend pas simplement du régulateur

de tension, mais aussi des fonctions de limitation et de protection qui permettent de
bien répondre aux exigences de performance présentées en section 2.2 [1] [12].

Mesure de tension
Compensateur de charge

Stabilisateur
de puissance ~

Mesure de
tension **

4 Rég ulateur
Ajus~~
1----
DC
tensi on
DC .1
Î
Excitatrice 1

1
:1 Court circuit 1 * :~ Génératrice ~ ---...
,... Régu late ur
1----
AC
Aju st.___....
tens ion AC

Exc. Syst
Circuits de stab.

Limiteur de ***
Surexcitation

Limiteur de
Sousexcitation

Limiteur V/Hz
et protection

* Court-circuit : Applicable pour l'excitatrice type ac et statique seulement

** Q uelques systèmes ont un régulateur de en boucle ouverte
*** Limite ur surexcitation : peut être utilisé avec un régulateur de

Figure 2-7 Système d'excitation muni de circuits de protection

 21

La Figure 2-7 présente toutes les différentes fonctions et circuits de protection

d'un système d'excitation. Cependant l'utilisation de ces fonctions n'est pas requise
par les différents types de système d'excitation, mais elles sont distinctes d'un
système à 1' autre [8].

2.5.1 Régulateurs AC et DC

Le régulateur à courant alternatif (AC) contrôle la tension de sortie au niveau du

stator de la génératrice synchrone, de plus il permet l'application d'autres contrôle et
de fonctions protectrices, comme illustré dans la Figure 2-7.

Le régulateur à courant continu (DC) contrôle à son tour, la tension d'excitation.

Souvent, il est définit par un régulateur automatique de tension (AVR) [15].

Par ailleurs, plusieurs types de régulateurs peuvent être utilisés dans le cas du
courant alternatif, citant les correcteurs conventionnels (PI, PID, ... ), qui seront
présentés dans le chapitre 4.

2.5.2 Circuits de stabilisation d'un système d'excitation

Les systèmes d'excitation comportent des circuits de stabilisation, qm

interviennent souvent en cas d'instabilité du système causé par un important temps de
réponse (le cas présent dans le système d'excitation de type DC ou AC) [6], ou
lorsque le gain du régulateur de tension au stator est très faible. Le circuit de
 22

stabilisation est caractérisé par une compensation additionnelle à la boucle de

régulation pour améliorer la performance dynamique du système. La Figure 2-8
illustre la compensation la plus fréquente, qui est définie par une fonction de transfert

sous la forme sK F . Les paramètres de la fonction de stabilisation peuvent alors

1+sTF
s'ajuster pour assurer le bon fonctionnement du système.

vref + VFD
Excitation
,."" L ..."" "" Vers la
,_ et AVR
Génératrice
v;
Compensation

sKF
l +sTF
.. "

Figure 2-8 Schéma de stabilisation et régulation de l'excitation [1] [2]

2.5.3 Stabilisateur de puissance (PSS)

Le stabilisateur de pmssance (PSS ) utilise des s1gnaux de stabilisation

additionnels qui agissent sur le systèm e d 'excitation dans le but d 'améliorer les
performances dynamiques du syst èm e électrique global. Généralement les grandeurs
utilisées à 1'entrée d'un st abilisateur de puissance sont génér ées par : la vit esse de
rotation de l'arbre, la fréquence et la puissance [2] [1 6].

L'utilisation d'un st abilisat eur de puissance dans les syst èmes d 'excitation, est
souvent requise, car cet élément permet essentiellement 1' am ortissement des
oscillations de puissance surtout dans le cas de faibles signaux, ce qui assure ainsi
une m eilleure st abilité.
 23

2.5.4 Compensateur de charge

Normalement, c'est le régulateur automatique de tension (AVR) qui assure le

contrôle de la tension de sortie de la génératrice. Mais parfois, un compensateur de
charge est aussi ajouté pour réguler la tension [1] [17]. Comme montré dans la Figure
2-9, le compensateur de charge est composé d' une résistance Re et une réactance X c

en série, qui sont ajustables. En utilisant cette impédance, le courant des armatures 11

mesuré et la tension E1 , le calcul de la tension de compensation Vc est donné selon la

formule suivante :

(2.1 )

Génératrice
Bus HT

[>-------+-.
synchron e Transformateur

3 0+~
3 ~
-~---------<]] [ Au système
d'alimentation

Rég ulateur
Exc itatrice
de tensio n

Figure 2-9 Schém a du compensateur de charge [ 1]

 24

2.6 Conclusion

Ce chapitre nous a perm1s de sav01r plus sur l'importance, le rôle et le

fonctionnement d'un système d'excitation dans un groupe turbine-alternateur. Cela
nous introduit à une étape importante dans notre travail, qui consiste en la
modélisation du système d'excitation du simulateur analogique du groupe turbine-
alternateur chez Hydro-Québec, ce qui est sujet dans notre prochain chapitre où le
choix du type d'excitation est fait et le modèle mathématique est ensuite élaboré.
 25

CHAPITRE 3

MODÉLISATION DES SOUS-ENSEMBLES DU GROUPE TURBINE-

ALTERNATEUR

3.1 Introduction

Dans ce chapitre, la modélisation de chaque bloc (génératrice, turbine et système

d'excitation) du groupe turbine-alternateur est présentée. L'objectif de cette partie est
d'illustrer chaque modèle indépendamment, et plus particulièrement présenter le
modèle mathématique du système d'excitation du groupe turbine alternateur.

Premièrement, le modèle mathématique de la génératrice synchrone est présenté,

dans le repère triphasé abc puis dans le repère de Park [1]. Le modèle de la machine
dit hybride le plus approprié dans notre travail de recherche, est développé.

Après la présentation du modèle de la turbine hydraulique, le système d 'excitation

de type statique ST est modélisé suivant deux méthodes qui seront bien détaillées,
selon tout élément qui constitue 1'excitation. Cela permettra par la suite, d'effectuer
différents essais et tests au niveau du système global après assemblage des trois
modèles, pour montrer l'importance de cette partie excitatrice.
 26

3.2 Modélisation de la génératiice synchrone

La machine synchrone considérée dans notre travail représente un alternateur

synchrone à pôles saillants, elle est principalement composée de trois enroulements
au stator, un enroulement au rotor et deux enroulements amortisseurs.

Après l'application de la transformée de Park aux équations de la génératrice, le

modèle mathématique est exprimé en unité réduite (pu) par les équations électriques
de tension et de flux ainsi que les équations mécaniques [ 1].

Dans notre travail, il s 'agit en partie de réaliser un contrôle des grandeurs de

sortie de la génératrice au niveau du système d'excitation. Le modèle hybride de la
génératrice qui est le plus approprié, sera alors développé.

3.2.1 Modèle de la génératrice synchrone dans le repère abc

Les relations entre les flux et courants s 'écrivent sous la forme matricielle
suivante:

(3.1)

Avec: aoc
lf's = [ lf'a lf'b lf'c r; lf':Wc = [If'f lf'D lf'Q r; (3.2)
z-aoc [.
i ab= = [i
s = za ib ir .
c ' r f iD iQr
 27

(3.3)

Et les équations de tensions sont écrites comme suit :

(3.4)

Telles que :

(3.5)

J [-Rs 0 ](-abcls J+-d [[ -Lalx (B)

Vsa1x - 3,3 ss sr (B)l (.a/x
Lalx ls J] (3.6)
(
va: - 0 3•3 Rr i':c dt ( -L~c ( B) Y L rr ia:
 28

Figure 3-1 Enroulements d'une machine synchrone triphasée (repère abc) .

3.2.2 Modèle de la génératrice synchrone dans le repère de Park (dqO)

La transformation de Park définie une transformation triphasée-diphasée

suivie d'une rotation. Elle permet de passer du repère fixe (abc) vers le repère mobile
(dqO), afin de simplifier le modèle du système vu la complexité de l'étude du
système.

Les enroulements au niveau du stator (repère abc) sont alors remplacés par deux
enroulements (dqO) en quadrature.
 29

d
2.
a
ji
/
/
/
/
/
o/
(("-/

Fi gure 3-2 Schéma de la m achine synchrone en diphasé (référentiel dqO)

La matrice de passage P est notée comme suit :

2rr 4rr
cos(B) cos(B - -) cos(B - -)
3 3
P(B) = . e-2rr
3_ - sin(B) - sm( . e -4rr
-) - sm( -) (3.7)
3 3 3
1 1 1
2 2 2

1
La matrice inverse F sera alors :
 30

cos(B) -sin(B) 1
. e - -2rr)
2rr
p-1(8) = cos(B--) -sm( 1 (3.8)
3 3
2rr
cos(B+-) . e +-)
-sm( 2rr
1
3 3

L'application de la transformation de Park se fait comme suit :

[Idqo ]=P[Iabc ]=[id iq i0 r (3.9)

J= p ['f'abel =
[ 'f' dqO [!pd !f'q !f'o r
[~qO J= p [V:,bc ]= [Vd Vq VO r
Et les équations obtenues sont données par :

./ Équations de flux :

(3 .10)

Telles que :

1
00 ·L
' sr
=PLaJx(B)=
sr
lL~ ifO
0
(3.11)

L oo 0 0

./ Équations des tensions :

 31

A partir de la formule des tensions dans le repère abc dans l'équation (3.6). Les
tensions au niveau du stator peuvent s'écrire :

v:qo =
[ vrdqO
l [-R l l l
0
5

3,3
03,3
R
r
s
[idqO
-dqO
lr
+ p [lf/dqO
s +m
dqO
lf/,
m
W [lf/dqO
s l
dqO
lf/,
(3.12)

(3.13)

La conversion des équations en valeur relative (pu) permet de faciliter le calcul et

l'évaluation du modèle mathématique plus aisément. Après avoir introduits les

valeurs de références (V:.réf ' I sréf ' S sréf' m rëf' •• •) de chaque grandeur, les résultats

alors obtenus sont donnés en valeur réduite (pu) dans les équations suivantes :

./ Expression des flux en pu :

(3.14)

Avec:
xd 0
x ss = 0 o0 J·x
' sr
[xdf
= o 0 (3.15)
[0 0 X QQ 0 0
 32

Telles que les réactances mutuelles sont calculées comme suit :

(3.16)
x ff =xmd +xkfi +x1 ; xQQ = x mq +xQ;

x df =xfd =xmd ; xQq = x qQ = x mq ;

./ Expression des tensions en pu :

(3.17)

(3.1 8)
 33

id ra xa xkfï xl rf if
~

Vd(t)

1 eq
1
VI
+
ûlml//q

i q ra xa
~
lmq = -lq +IQ
rQ
xmq
Vq(t)

1
~+
e,
XQ

Figure 3-3 Circuits électriques équivalents de la génératrice synchrone

triphasée dans le repère de Park (en pu)

Pour obtenir une modélisation complète de la m achine, nous adj oignons aux
équations électriques, l ' ' equation mécanique de la machine régie par la relation
suivante :

dw (t)
2H m = T - T - Dw (pu) (3 .19)
dt m e m

Le couple électromagnétique en pu est calculé par :

T = (!ffdi - !f/ id )
e q q (3 .20)
 34

3.2.3 Développement du modèle hybride de la génératrice synchrone

Le but principalement visé étant de contrôler, stabiliser et optimiser les valeurs

des grandeurs (tension, courant) de sortie au niveau du réseau électrique. Le modèle
hybride de la machine synchrone qui est le plus adapté dans ce cas, est alors formulé
à partir des équations de tensions et de flux obtenues précédemment.

Tensions au rotor V,. ..

~
Modèle
..
~
Tensions au stator V:
hybride
Courants au stator 1s .. de la GS ...
~
Courants au rotor 1r

Figure 3-4 Illustration des entrées-sorties du modèle hybride de la

génératrice

Les expressions de tension et flux sont résumés dans :

1) Équation de t ension :

~ 3,3] [lf/s] (3.21)

3,3 If-Ir

A partir de 1'équation (3 .21 ), on déduit :

 35

(3.22)
(3.21) ::::>

(3.23)

2) Équation de flux :

[!f/s] =[~:s~
lf'r sr
~sr][~s]
rr r
(3.24)

(3.25)
( . ) ::::> {!f/s = -XsJs + XsJr
3 24
!ffr =-xsr T ls + XrJ r (3.26)

À partir de l'équation (3.23) on aura:

(3.27)

De l'équation (3.26), on obtient:

-X (3.28)
fr - rr
-l (
l.f/, +X sr T] s )
1r =Xrr -l~~~
r r
+X r r - l xsr ri s (3.29)
 36

En remplaçant (3.29) dans (3.27), on aboutit à:

(3.30)

(3.31)

En prenant le flux comme variable d'état.

dtpr
dt l.f/r

Figure 3-5 Schéma foncti onnel des états (flux)

~ Calcul des sorties :

À partir de l'équation (3.22)

 37

(3.32)

En remplaçant (3.29) dans (3.32), on obtient :

lf.ls
=-Xss1s +Xsr x rr- l l.f/, +Xsrx rr- l xsr ri s (3.33)

~ -xsrx rr- l lf', +(xsr x rr- l xsrr - x ss )1s

lf/s -
(3.34)
'------v------'
C D

Figure 3-6 Schém a fonctionnel des sorties ( l.f/s , ~ )

 38

3.3 Modèle de la turbine hydraulique et servomoteur

Dans une turbine hydraulique, les vannes sont actionnées par un servomoteur.
L'eau pénètre dans une conduite forcée pour augmenter sa vitesse. À la sortie de la
conduite, elle fait tourner une turbine qui entraîne un alternateur. L'eau est ensuite
évacuée par une galerie de fuite.

La modélisation de tous les composants de la turbine hydraulique aboutit à un

modèle non linéaire décrit par les équations de la dynamique du débit et la puissance
mécanique [5] :

(3.35 ~

Le modèle linéaire est obtenu par la linéarisation autour d'un point de

fonctionnement (g0 ,q0 ,P0 ) [18] (annexe 2), on aura alors les équations suivantes.

(3.36)
 39

En utilisant les équations suivantes au point de fonctionnement :

(3.37)

On obtient facilement le modèle linéaire de la turbine hydraulique avec le

servomoteur comme suite :

(3.38)

Où q, g,Pm ,4,Tw,~ ,Ka , et u représentent respectivement le débit d'eau dans la

turbine, le vannage, la puissance mécanique, gain de la turbine, le temps

d'établissement du débit nominal, la constante de temps du servomoteur, le gain et la
commande du servomoteur [5].
 40

3.4 Modélisation du système d'excitation statique

Dans notre travail, nous avons opté pour un système d'excitation de type statique
à source de tension avec redresseur commandé présenté dans le chapitre 2 (Figure 2-
6). Celui-ci correspond parfaitement aux attentes du simulateur analogique d' Hydra-
Québec.

Dans cette partie, nous développons le modèle mathématique utilisé de

1' excitatrice statique suivant deux méthodes : le modèle par méthode DQ et le modèle
type ST de IEEE.

3.4.1 Modèle par la méthode DQ

Cette méthode consiste en la représentation du redresseur commandé de

l' excitatrice statique (Figure 2-6), par le circuit électrique dans le cas d'une
alimentation triphasée de la Figure 3-7 [ 19] [20].

-------1
1 rF

~------- A --------~ 1 +
Source Bus ACBus 1

CPL

Pont à thyristors

Figure 3-7 Schéma du convertisseur commandé lié au réseau

 41

Le redresseur est représentée par un pont à thyristor, la ligne de transmission est

définie par les impédances équivalentes (Req, Leq et Ceq), et le filtre par (rF, LF et CF).
La sortie du filtre est souvent reliée a une charge idéale CPL afin d'évaluer les
grandeurs de sortie du redresseurs commandé.

Tel que le courant de sortie est donné par la formule suivante :

1CPL -- PCPL

vi (3.39)

Se basant sur le circuit équivalent dans la Figure 3-9, l'effet des inductances
équivalentes (Leq) en alternatif entraine une chute de tension après redressement, qui
3wL
définie par est résistance variable r = __e_q ; où w est la fréquence de source .
.u 7r

La modélisation du redresseur commandé dans le repère dq est simplement

développée sous forme d'un transformateur (Figure 3-8) en utilisant les équations de
transformation qui vont suivre.

r
3CùL
= --
eq r-----------,
~ 7r
1 r~ 1
---......
--:---;1,------, +
+
+ --+ + Id 1
Ide vbus,d
1
la Repère dq
Vbus,a
vbus,b lb Ede
c==J +
vbu,,c Je vbus ,q

1 1· s
L --~ ~ -------1
Pont à thyristors Transformateur

Figure 3-8 Circuit équivalent du redresseur commandé dans le repère dq

 42

./ Modèle mathématique dans le repère abc :

(3.40)

Telle que:

(3.41)

Où ~ :Angle de phase des tensions de bus et a: Angle d'amorçage des thyristors.

S abc : Matrice de passage dans le repère abc .

./ Modèle mathématique dans le repère dq :

Le passage du modèle mathématique du redresseur commandé dans le repère abc

s'effectue par la matrice de transformation T[ B(t)] définie par la formule suivante :

2 2
cos( B(t)) cos( B(t)- ; ) cos( B(t)+ ; ) (3.42)
T[B(t)]=
2 2
-sin( B(t)) -sin ( B(t)- ; ) -sin ( B(t) + ; )

7r (3.43)
Où B(t) = wt--+<A,
2

Tel que w représente la fréquence du signal de source.

Et ~1 : Angle de phase des tensions dans le repère dq.

 43

En combinant ainsi, les équations (3.40)- (3.42), on obtient :

(3.44)

La matrice de passage S dq dans le repère dq, est donnée par :

s dq =
2../3[ cosC<A -<P+a)
.--
. 1 -<P+a) ]T
-sm(</J (3.45)
7[

Finalement, en utilisant le modèle par la méthode DQ, se basant sur la matrice de

transformation T[ B(t)] définie en (3.42) et le circuit équivalent du redresseur de la

Figure 3-8, nous pouvons représenter le redresseur commandé lié au ré seau de la

Figure 3-7 par le circuit équivalent suivant dans la Figure 3-9.

rf1 rF

+ --+ +
Id
v sd +
vbus,d c.q t IVC<f11Vrus,q vf

vsq +
-

Figure 3-9 Circuit électrique équivalent du redresseur commandé lié au

r éseau dans le repère dq [1 9]
 44

./ Modèle d'état rmal du convertisseur commandé:

A travers le circuit électrique équivalent du redresseur commandé Figure 3-9, et

après l'application de la loi de kirchhoff pour les grandeurs électriques de tension et
de courant, nous aboutissons aux équations suivantes [19] [20] :

(3.46)

Tel que:

Variables d'états: X= [Ids ] qs V:bus ,d V:bus ,q ] de V1 Jr

La variable de commande d 'états : U = [Vm PcPL r

La variable de sortie : Y= V1

Le système d'équation (3.45) est non linéaire, la linéarisation se fait alors par le

terme du premier ordre des séries de Taylor autour d'un point d'équilibre ( X 0 , U 0 )

[18] (annexe 2).

 45

Le système d'état après la linéarisation sera alors défini par:

c5X:A(X0 ,U 0 )c5X +B(X0 ,U 0 )c5U

{ c5Y -C(X ,U )c5X +D(X ,U )c5U
0 0 0 0

Où 1'état, la commande et la sortie sont respectivement :

oX = [of ds ofqs oVb"",d o~us,q ofde oV1 r

ou= [ovm oPcPLr
oY =oV1

Les matrices A, B, C et D sont données comme suit :

1
(û -- 0 0 0
L.q L.q
- R.q 1
- (û 0 -- 0 0
L.q L,q (3.47)
0 0 (û ~ 2~ 0
c.q 7!C,q
A(Xo,Uo) =
0 - (û 0 0 0
c .q

0 0 ~2~ 0 - ( !L+!E_ J
J!Leq LF LF LF
PCPL
0 0 0 0
CF CF.V2j0 6x6

~ cos(Âu+a0 ) 0
Leq C(X0 ,U0 ) = (0 0 0 0 0 1 L6
~ sin(Âu + a 0 ) 0
Leq D(X0 ,U0 ) = [0 0 L2
B(X0 ,U0 ) = 0 0
0 0
0 0

0
 46

./ Modèle d'état du convertisseur commandé avec régulation de l'excitation

[20] [21]

Afin de réguler la tension d'excitation Vj, des régulateurs classiques PI peuvent

être utilisés. Le circuit électrique du convertisseur statique de la Figure 3-7 avec
régulation devient comme montré à la Figure 3-10 [20].

~------- À --------~
Somce Bus AC Bus

Contrôle + v·j
des sig naux

Figure 3-10 Convertisseur avec régulateurs

Le schéma fonctionnel du régulateur PI est donné par :

Figure 3-11 Régulation par PI

 47

A partir du schéma de régulation PI (Figure 3-11 ), on déduit les équations suivantes

(3.48)

Le circuit électrique équivalent du convertisseur statique avec régulation de la

tension d'excitation, dans le repère dq est alors le suivant :

,.1

-
1q =0
ae.,v~u.d

Figure 3-12 Circuit électrique équivalent du convertisseur régulé [20]

A partir du circuit équivalent, on peu écrire :

(3.49)

En remplaçant v* (3.49) dans l' équation (3.48), on obtient

 48

(3.50)

En intégrant les équations additionnelles (3.48) et (3.50) obtenues après

régulations, Les équations alors en (3.45), deviennent [20] :

1 1
L 1ds + w1qs - -L - Vbu.s,d + -L - Vsd
P(1ds ) = - R,q
q q q

R
1 1
P(1qs ) = -wlds -~1
L qs -L- Vbus ,q +L - Vsq
eq eq eq (3.51)
1
p(Vbu.s,d ) =C- 1ds -wV:bus,q -
~2../3
-- - 1de
:Tr.C
q
2 ~

1
p(v;,us,q ) = -wv;,u.s ,d + c lqs
eq

V ) = -1-1 __1_ PCPL

p ( .( C de C V
F F .(

p(XJ = - V.r +V/

p(X) = - l d - K pv V.r+ K IV. X +K pv V.r*
: C V

Le système d 'équation dans (3.5 1) est non linéaire. Après linéarisation autour d'un
point d 'équilibre (annexe 2). Le systèm e d' état sera définit comme suit :

5 X : A,eg (X 0 ,U0 )5 X + Breg (X0 ,U0 )5U

{ 5Y - Creg(X ,U )5X + Dreg(X ,U )5U
0 0 0 0

Où 1' état, la commande et la sortie sont écrits comme suit :

 49

oX =[olds sv.bus,q oldc ox,r

oU =[ovm sv·f
oY = sv1

Et les matrices 4 eg , Breg , Creg et D reg sont données par :

 50

./ Exemple de simulation du modèle par méthode DQ :

Afin d'illustrer la dynamique du convertisseur commandé se basant sur le modèle

obtenu par la méthode DQ, nous allons présenter quelques tests de simulations
obtenus dans matlab (l'exemple d'application numérique et le scripte associé sont
donnés en annexe 3 et 5 respectivement).

Pour une tension d'entrée V.. = 230 f/,_msl phase, et une fréquence f =60hz.

1) Redresseur lié au réseau, sans régulation de la tension d'excitation Vf (Figure

3-9)
• Pour une variation de charge à la sortie Pcpl de 2 kW à 4 kW, opérée à t =

0.5 s, avec un angle d'amorçage des thyristors a = 10 degrés:

Tens ion Ex citatrice

190
4 kYI
185 L KW

en 1\
±o! 180

\ 1~h--
0
6
>c 175
.Q
en
v
c
Q)
170
1-- ~

165

160
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9
Temps(s)

Figure 3-13 Tension de sortie Vj(a = 10°)

 51

• Pour une variation de charge à la sortie Pcpl de 2 kW à 4 kW, opérée à t =

0.5 s, avec un angle d'amorçage des thyristors a = 15 degrés:

Tension Excitatrice
150
4 kw
145
2 kw
140
:§' ~
0 135
- ~

\
6
>c 130
.Q
(J)
(\ 1\ /'\..
\1- \ ~v ~ -
'-"'""
cQ)
f-
125

120
v
115
v
110
0.3 0.4 0.5 0.6 0. 7 0.8 0.9
Temps(s)

Figure 3-14 Tension de sortie Vj(a = 15°)

Remarque:

Selon les Figures 3-13 et 3-14, plus la charge Pep! augmente, plus la tension
de sortie Vjdiminue.

Plus l'angle d'amorçage est grand plus la tension de sortie Vj diminue, et le

système devient de plus en plus lent et instable (Figure 3-14).

Les résultats obtenus exigent forcément un contrôle de la sortie pour remédier

à l'instabilité du système.
 52

2) Redresseur lié au réseau, avec régulation de la tension d'excitation Vf

Après la synthèse des correcteurs PI (Figure 3-11) comme suit, à partir des
schémas en Figure 3-15 et Figure 3-16 [20] :

Figure 3-15 Calcul des paramètres Kpi et Kii

Le système étant du 2éme ordre sous la forme s 2 + 2.;-wn.s + wn2

Tel que .; et (J)n r eprésentent le coeffi cient d' amortissem ent réduit et la pulsation

propre du syst ème respectivem ent. On déduit alor s les paramètres du PI :

(3 .53)

Figure 3-16 Calcul des paramètres Kpv et Kiv

 53

(3.54)

• Pour une variation de charge à la sortie Pcpl de 3 kW à 4 kW, opérée à t =

0.3s, avec un angle d'amorçage des thyristors a = 10 degrés, avec une

consigne de tension V1 • = 300 volts.

Tension Excitatrice
302

4 kw
301
3 kw

2 300

~r
Ci
:2
>c 299
v
0
Vl
c
~ 298

297

296
0.1 0.1 5 0.2 0. 25 0.3 0.35 0.4 0.45 0.5
Te mps(s)

Figure 3-17 Tension de sortie Vj régulée

D'apr ès la Figure 3-17, le Régulateur PI m aintient la tension d'excitation Vj à la

valeur désirée Vj * = 300 Volts, m ême avec une variation de la charge extérieure Pcpl.
 54

./ Remarques

La modélisation du système d'excitation par la méthode DQ présente des résultats

satisfaisants. Cette méthode permet la régulation de la tension de sortie du
convertisseur classique [20] [21], cela dit, la complexité des calculs et des équations
du système d'état qui représentent le convertisseur statique ainsi que l'intervention de
plusieurs paramètres, implique souvent, un temps de simulation assez important. Ce
type de modèle incluant le convertisseur statique est souvent utilisé dans
l'amélioration de la stabilité due à une charge CPL [21]. Il introduit une étude
intéressante dans les systèmes plus complexe de production, tels que les génératrices
synchrones.
 55

3.4.2 Modèle IEEE type ST

Les modèles mathématiques des systèmes d'excitation sont essentiels pour

l'évaluation des performances exigées, pour l'optimisation et la conception ainsi que
la coordination de commande supplémentaire et de circuits de protection [1] [12], et
aussi important pour les études de stabilité de systèmes liées à la planification et à
l'exploitation des systèmes d'énergie électriques [6] [9] [11]. Le détail du modèle
requis dépend souvent du but de 1'étude à faire [22]. Dans notre travail, le contrôle et
les dispositifs de sécurité qui ont un impacte sur 1' étude de la stabilité transitoire et
les petits signaux sont principalement définis par le régulateur de tension, le
stabilisateur de puissance et la commande et stabilisation de l'excitation [1] [9].

Dans cette partie, la modélisation du système d 'excitation statique de type ST1

standardisé par IEEE est présentée [1] [23]; le modèle de chaque élément est illustré
indépendamment puis assemblé pour présenter le modèle globale de l'excitatrice.

~ Modèle de chaque élément du système d'excitation type ST [1]

./ Modèle de l'amplificateur:

Le pont à thyristor de l'excitatrice statique est souvent modélisé par un

amplificateur tout en prenant compte des circuits de stabilisations et de limitations qui
en découlent. L'amplificateur peut être de type magnétique ou électronique, il est
caractérisé par un système de premier ordre incluant un gain KA et une constante de
temps TA (Figure 3-18).
 56

VIM4X

vlM'/
Figure 3-18 Modèle de l'amplificateur

La tension de sortie de l ' amplificateur est limitée par une saturation, spécifiée par
les bornes VRMAX et VRMIN dites les limites de saturation.

La limitation de tension de sortie de l'amplificateur varie souvent en fonction de

la tension externe de la génératrice (saturation dynamique) .

./ Modèle du circuit de stabilisation :

Excitation

Armature
Champs ~----~----------~+

r~
Transfo de
stabilsation

L...--~+
+
V:rr
L...-----~~-------------------------10-

Figure 3-19 Système d 'excitation et transformateur de stabilisation [1]

 57

Dans ce type de système d'excitation, les transformateurs série sont les plus
utilisés tel que montré dans la Figure 3-19 ci-dessus. Les équations liées au
transformateur de stabilisation sont données dans le domaine de Laplace comme suit :

{V~ =R1i~
2
= + s~i1 .+ sMi
JS z + s4.z + sMz
2 2
2.

1
(3.55)

Où les indices 1 et 2 désignent le primaire et le secondaire du transformateur, et R,

L et M représente la résistance, l'induction de fuite et l'induction mutuelle,
respectivement. Étant donné que le secondaire du transformateur est connecté à un
circuit à haute impédance, alors i2 est négligé, on aura donc :

~ : (1\ + s4 )i1
{ Tl;_ - sMz 1
(3.5 6)

Alors:
V2 _ sM1 _ sKF
(3.57)
~ R1 +s4 1+ sTF

Où K =M/ et T
F Il\ F
=4/
~~
 58

./ Les limitations :

Dans la modélisation des systèmes d'excitation de type ST, les circuits limiteurs
jouent un rôle important, cependant il est nécessaire de distinguer entre deux sortes de
limitation, la limitation après bloc (windup), et la limitation en bloc (non-windup),
pour ce faire nous illustrons un exemple d'une fonction d'intégration de l'entrée dans
chacun des deux cas, présentée dans la Figure 3-20 ci-dessous :

Lx
Lx

u
·1
1
s
lv/.
J
y u 1
s
y

LN LN
(a) (b)

Figure 3-20 Fonctionnement des limitations

A partir de la Figure 3-20, on a :

dv
(a) => - = U
(3.58)
dt

La sortie y est calculée selon les cas suivants :

(3.59)

(3.60)

Y -L
- N (3.61)
 59

dy
(b) => -=u (3.62)
dt

La sortie y est calculée comme suit :

dy
Si LN --< y --< L x => -=u (3.63)
dt

Si y'2L et dy >-0 => dy =0 et y=Lx

x dt dt (3.64)

Si y<S,LN et dy--< 0 => dy =0 et y=LN

dt dt (3.65)

Pour les limitations après bloc (externe), la variable v n'est pas limitée, par
conséquence la variable de sortie y suit la valeur de v, jusqu'à ce que v atteigne les
bornes Lx et LN. Alors que pour les limitations en bloc (interne), la variable de sortie y
est limitée et prend la valeur des bornes dés que u change de signe .

./ Les fonctions de déclenchement «Gate » :

Des circuits de déclenchement sont souvent utilisés pour contrôler un ou deux

signaux d'entrée. La Figure 3-21 illustre les fonctions de basse valeur (Low Value,
LV gate) et de grande valeur (High Value, HV gate) :
 60

:1 .~i.) :1 .~.)
u u
~y ~y
v v

Si u ~v, y=u Si u :::; v, y=u

Si u --<v, y=v Si u >-v, y=v

Figure 3-21 Fonctions de déclenchement (Gate)

•!• Modèle globale du système d'excitation statique:

La Figure 3-22 ci-dessous représente la structure générale d'un modèle détaillé de

systèm e d'excitation ayant une correspondance pour chaque équipem ent physique.
Pour simplifier les ét udes sur ce modèle, plusieurs techniques sont adaptées
dépendamment de 1' étude souhaitée.

Étages d'amplification Excitation

Autres
Boucles mineures
signaux
de stabilisation

Boucle majeure de
stabili sation

Figure 3-22 Structure du m odèle détaillé du système d' excitation

 61

./ Modèle du système d'excitation statique type ST1A

Le modèle du système d'excitation statique de type ST1A est définit par un

redresseur commandé à source de tension tel que présenté en chapitre 2 (Figure 2-26).
Le courant d'excitation du redresseur commandé est fourni par la génératrice
synchrone à travers le transformateur de puissance, les valeurs de la tension
d'excitation alors sont proportionnelles à celles de la tension de sortie de la
génératrice, la constante Kc est définie au niveau des limitations de VFD (saturation
dynamique). La limitation du courant d'excitation est aussi prise en compte par
1'introduction de la constante hR et le gain KLR· La Figure 3-23 ci-dessous présente ce
type de système d 'excitation en toutes ses parties principales [ 1] [23] [24] [25].

* Entrées altern ées

Figure 3-23 Modèle IEEE du système d' excitation statique type ST1 A
 62

3.5 Conclusion

Cette section du document, nous a présenté toute la partie modélisation

mathématique, et mise en équation de chaque élément constitutif d'un groupe turbine-
alternateur; génératrice synchrone, turbine hydraulique et système d'excitation.
Cependant, 1'élément essentiel dans notre travail, défini par 1' excitatrice, comporte
des parties de commande, régulation et stabilisation, le chapitre suivant va présenter
ces différentes parties.
 63

CHAPITRE 4

ANALYSE DES RÉGULATEURS CONVENTIONNELS D'UN GROUPE

TURBINE-ALTERNATEUR

4.1 Introduction

Ce chapitre traite les différents équipements qui interviennent principalement au

niveau du système d'excitation d'un groupe turbine-alternateur. Ces éléments dits
régulateurs conventionnels sont utilisés afin d'assurer le bon fonctionnement du
groupe. On distingue le régulateur automatique de tension AVR qui permet de
contrôler la tension d'excitation, le régulateur classique PID de la tension terminale
génératrice, le régulateur de vitesse (speed govemor) qui lui contrôle la vitesse de
rotation et enfin le stabilisateur de puissance (PSS) qui intervient en principe dans
l'amortissement des oscillations de puissance au niveau du réseau.

4.2 Régulateur automatique de tension

La sortie de l'AVR dépend uniquement de la tension terminale. Il agit sur

l'excitatrice du générateur afin de maintenir la tension terminale constante [1] [25].
La Figure 4-1 représente le schéma de principe de la régulation de tension dans une
génératrice.
 64

vrej ~Ir
v. v
+
R
_____..
FD
~.
, , Régulation Excitation .."" Générateur ..
~t. -

vc

Capteur

Figure 4-1 Principe de la régulation de tension dans un groupe turbine-

alternateur [ 1] [25] [26]

4.2.1 Régulation et stabilisation de la tension d'excitation

Le régulateur de tension automatique A VR et 1'excitatrice constituent

principalement le système d'excitation du générat eur [ 1] [26]. Comme vu dans le
chapitre 2, il existe plusieurs t ypes de systèmes d'excitation. Le régulateur de tension
a pour entrée la t ension de réf érence qui constitue la consigne, la tension terminale
provenant de la machine synchrone, une t ension de stabilisation provenant du
stabilisateur de puissance et enfin un retour après amortissement (compensation). La
différence entre la t ension terminale et la tension de référence passe par un
compensateur de phase puis par le régulateur principal. C i-dessous est représenté le
système d'excitation statique de t ype ST1 standardisé par IEEE.
 65

AVR
~---------------------~

1 VAMAX
1

1
Compensation Excitation V1 VJIMIN
VAMIN 1
retard-avance
Amplification :
~----------------------

sKF
~----------------~ ~----------------------~
l+sTF
Compensation

Figure 4-2 Système d' excitation type STl de IEEE

Afin de contrôler la tension d'excitation VFn, le système simplifié (Tc, Tc1 , TB, et
TB 1 sont négligés) peut être représenté comme suit [1] [2] :

r---------
1
V FD Vers
la GS

Amplification Excitati on
_ _ _ _ _ _ _ _ _ __J

Excitati on et AVR

Compensati on

Figure 4-3 Régulation et stabilisation de la tension d' excitation

 66

4.2.2 Régulation de la tension de sortie de la génératrice synchrone :

Afin de corriger la tension de sortie de la génératrice Vs, un régulateur de type

PID classique est le plus souvent intégré au système vu son efficacité de donner des
réponses plus performantes tellement il associe les trois actions principales
(proportionnelle, intégrale et dérivée) [25] [26].

Système d'excitation

~ef .:. Î: V:rr Régulation Vc l +~ KA VA KE 1 VFD KG vs...

,.. ':,/ ,..
PlO i~ 2: ~ ,.. ,..
~l'- 1 - l' l+TAs l+TEs 1 l+Ta s
1 1
Amplification Excitation Génératrice
1 v; 1
1
sKF ...,
1
1 1
1 l+sTF "' 1
1 1
Compensation
1
--------------------------
1

VSr _K_r_ ...,

~

1+ Tr S
Capteur

Figure 4- 4 Régulation de la tension de sortie de la génératrice

K de(t )
d dt

Figure 4-5 Schéma fonctionnel d'un correct eur PID classique

 67

La tension de sortie corrigée est donnée par :

Vc(t)=KP[e(t)+_!_e(t)dt+TD de(t)l (4.1)

0 dt

Tel que e(t) représente l'écart, Kp le gain proportionnel, T1 la constante intégrale

de temps et T D la constante dérivée du temps.

La fonction de transfert du correcteur PID s'écrit alors comme suit :

1
G Pw (s) =KP(1 +-+Tv s) (4.2)
~s

./ Synthèse des paramètres des correcteurs :

vrej~+
~ fl.rr ...,. Régul ati on vc
,....
KI VFD ...
,
KG vs ...
,.
PID 1+T1 s 1+Ta s
J~-

Génératrice

VSr ~ ""
1+ T, s
Capteur

Figure 4-6 Schéma de régulation de tension simplifié

 68

Méthode par compensation de pôles :

Après simplification du système global précédent (Figure 4-6), le schéma

fonctionnel obtenu montré ci-dessus nous permet de calculer aisément les paramètres
du correcteur PID, en utilisant la méthode par compensation de pôles.
La fonction transfert du correcteur PID peut s'écrire sous la forme :

(4.3)

La fonction de transfert (Excitation + Génératrice) est exprimée comme suit :

(4.4)

Après compensation de pôles par les zéros du correcteur PID, on aura :

(4.5)

En boucle ouverte, la fonction de transfert est exprimée comme :

K
GBO ( s) = ___§Q_
T (4. 6)
Is
 69

En boucle fermée :

- GBO ( s) H ( s)
FBF
(S ) - (4.7)
l+GBO ( s ) H ( s )

Tel que H ( s) définit la fonction de transfert du capteur qui s'écrit comme suit:

K
H(s)=-' (4.8)
l+T,s

En remplaçant les équations (4.6) et (4.8) dans (4.7), on déduit la fonction de transfert
du système globale en boucle fermée :

(4.9)

La forme canomque d'une fonction de transfert du deuxième ordre est la

suivante :
 70

(4.10)

Tel que E et mn représentent respectivement le coefficient d'amortissement et la

pulsation propre du système.

Par superposition entre les deux fonctions de transferts FBF(s) équation (4.9) et
F(s) équation ( 4.10) on aura :

0 2&
-- (4.11)
K BOK r aJn

TI~ 1
= ( 4.12)
aJ2
K BOKr n

(4.13)

2
TT
( 4.12) => K BO = I raJn (4.14)
K r

En remplaçant (4.13) dans (4.14), on obtient:

 71

On déduit alors les gains :

( 4.15)
ou

v' Calcul des différents paramètres de chaque fonction de transfert :

Fonction de transfert simplifiée de la Génératrice synchrone [ 1] :

( 4.1 6)

Le calcul de K 0 et T0 se fait comme suit :

G(s ) = K 1+s.J;_d /3.1f/ds = G(s)*efd

0
(1 +sT;0 )(1 +sT;)
0

(4.17)

avec

Tm et r: sontnégligées [1].
0
 72

Fonction de transfert de l'amplificateur:

( 4.18)

La relation entre la tension continue VsDc obtenue après redressement par le

convertisseur commandé, et la tension alternative entre phase Vsah reçue par le réseau
à travers le transformateur d'excitation est donnée comme suit:

VSoc = 1.35Vsab_eff cos( a)

V1 : 1.35V2 cos ~a)= 1.35KeVc cos (a) (4.19)

{ V:, - V:b cos( a)- KtbVc cos( a)

Où: v1 est la tension d'excitation (en courant continu), v; la valeur efficace de la

tension secondaire du transformateur utilisé, ~ tension corrigée à la sortie du

régulateur PID, et V:b définit la tension génératrice maximale (Vbase). Et Ke = Yvc et

Ktb = Vive représente respectivement le ratio entre excitation/transformateur et le

ratio entre génératrice/transformateur [2].

Le gain du K 1 de l'amplificateur est alors calculer comme suit :

(4.20)
 73

On déduit alors,

(4.21)

Et la constante de temps

(4.22)

Tels que:

a : Angle d 'amorçage des thyristors

T : Période du signal filtré

Te : Période du signal d ' entée

 74

4.3 Régulation de la vitesse de rotation

Le schéma ci-dessous (Figure 4-7) montre le principe de la régulation de la vitesse

de rotation et la puissance mécanique d'un groupe turbine hydraulique-alternateur. La
boucle de régulation est composée généralement d'un régulateur (governor), d'un
servomoteur et d'une turbine.

w rej +
... werr ...
Governor ___.. Servomoteur r---. Turbine
w ..
,.. Générateur _____.
,_
w

Figure 4-7 Principe de la régulation de la vitesse

Le régulateur (governor) est souvent représenté par un régulateur classique PID

[5] (Proportionnel, intégral et dérivée) dans la Figure 4-8. Le PID agit sur la position
du vannage de la turbine en fonction de la variation de la vitesse par rapport à sa
valeur nominale. Le statisme a pour rôle le partage de la charge totale du réseau entre
différents générateurs s' il y a lieu. Il permet aussi d'assurer le bon fonctionnement de
plusieurs turbines dans un réseau électrique. Un exemple de calcul des paramètres
Kp, Ki et Kd du régulateur PID classique de vitesse est donné en (Annexe 4).
 75

Kp

Servomoteur l-------1r--~ Vannage Turbine

Statisme

Figure 4-8 Régulateur de vitesse type PID

Le modèle du servomoteur est décrit par une fonction de transfert simple du

K
premier ordre __a _ semblable à celle d'un moteur. Il génère la vitesse de rotation
Tas+l
de la turbine en commandant 1' ouverture du vannage. Il est représenté dans le schéma
suivant [27] :

gmax
Ka
Tas+I
Servomoteur

Figure 4-9 Schéma bloc du Servomoteur

 76

4.4 Stabilisateur de puissance

Le stabilisateur de puissance (PSS) permet en principe d'amortir les oscillations de

puissances. Il fonctionne avec le régulateur automatique de tension (AVR) en
ajoutant une entrée supplémentaire au système d'excitation. La Figure 4-10 illustre la
configuration générale d'un PSS avec un système d'excitation de type ST1 de IEEE.
Le stabilisateur de puissance utilise comme entrée l'écart entre la vitesse de la
machine et la vitesse nominale w- Wref et a pour sortie la tension de stabilisation Vstah

[10] [24].

.----------l Stabilisateur de puissance PSS 1 4 - - - w- wref

( V, vRMAX
VAMAX
- KJ FD)
(1 + sTe )(1 + sTe, )
(l+ sT8 )(l + sT )
Compensation v,vRMJN
VAMIN
retard-avance
Amplification

sKF
~---------~ ~r-----~
l +sTF
Compensation

Figure 4-10 Syst ème d'excitation type ST1 avec PSS

Il existe plusieurs types de stabilisateurs utilisés dans la commande des réseaux

électriques. Le plus important est le m odèle conventionnel à une entrée type IEEE
PSSlA [24] [28] présenté en Figure 4-11 , il est constitué d'un gain Kstah qui détermine
le niveau d'amortissement du PSS, d'un filtre passe-haut qui va permettre au PSS
 77

d'agir seulement en cas de variation rapide de vitesse et éliminer l'action de ce dernier

lorsque la vitesse de la machine est constante et enfin un compensateur de phase.

w-wref
•.J K
1 stab
1 ...,... Tws -----. ~ns+1 ------. Tzas
T ns + 1
2 ,..... ~tab
1
Tws+1 ~as+1 +1

Figure 4-11 Stabilisateur de puissance de type IEEE PSSlA

Le schéma suivant (Figure 4-12) montre la régulation globale dans les différents
blocs du système (groupe Turbine hydraulique - Alternateur synchrone avec leurs
système d'excitation STl).

Excitation

Vannage
Machine
Synchrone et
Turbine

Figure 4-12 Schéma de la régulation globale

 78

4.5 Conclusion

L'action de contrôle des grandeurs de sortie d'un groupe turbine-alternateur est

donc obtenue par deux niveaux de régulations principales, soit les régulateurs de
tension et de vitesse (voir Figure 4-12), ajoutant à cela un stabilisateur de puissance
(PSS) qui agit pour fin de l'amortissement des oscillations de puissances
électromécaniques [9]. Cependant l'interaction entre le régulateur de tension (AVR)
et le stabilisateur de puissance (PSS) qui opèrent simultanément au niveau du système
d'excitation, peut compromettre la performance dynamique entre leurs deux fonctions
(contrôle de la tension et amortissement de la puissance) [24] [28]. Pour remédier à
tout inconvénient pouvant survenir à ce niveau, un travail d'optimisation des
paramètres de l'excitatrice sera effectué.

Dans le chapitre 5, nous allons alors présenter la méthode d'optimisation la plus

adaptée dans ce cas-ci, dite optimisation par algorithme génétique afin d 'assurer le
contrôle optimale des grandeurs de sortie, et améliorer les performances dynamiques
du groupe-turbine alternateur du simulateur analogique.
 79

CHAPITRES

OPTIMISATION PAR ALGORITHME GÉNÉTIQUE

5.1 Introduction

Ce chapitre présente la partie optimisation, où est définit la méthode et la

technique d'optimisation utilisée, appelée méthode d'optimisation par algorithmes
génétiques.

Ce choix est justifié par la nécessité de trouver l'optimum global d'un problème
d'optimisation qui comprend des variables continues ou discrètes, et qui est
représenté sous forme d'un système (système d'excitation statique ST) dont le
modèle est non linéaire, non dérivable ou ne possédant pas d'expressions analytiques.
L'optimisation par algorithmes génétiques est la méthode la mieux adaptée dans ce
cas-ci, car elle prend en compte toutes les différentes contraintes d'optimisations, et
de plus de sa robustesse, elle offre souvent des résultats très fiables.
 80

5.2 Technique d'optimisation par algorithmes génétiques

Les algorithmes génétiques (AG) sont des méthodes d'optimisation. Ils tirent leur
nom de 1'évolution biologique des êtres vivants dans le monde réel. Ces algorithmes
cherchent à simuler le processus de la sélection naturelle dans un environnement
défavorable en s'inspirant de la théorie de l'évolution proposée parC. Darwin. Tels
que les « individus » les mieux adaptés tendent à vivre assez longtemps pour se
reproduire alors que les plus faibles ont tendance à disparaître (the survival of the
fittest); Michalewicz et al., 1999; Tomassini, 1999) [29].

Par analogie avec l'évolution naturelle, les AG font évoluer un ensemble de

solutions candidates, appelé une « population d'individus ». Un « individu » n'est
autre qu'une solution possible du problème à résoudre. Chaque individu de cette
population se voit attribuer une fonction appelée fonction d'adaptation (fztness) qui
permet de mesurer sa qualité ou son poids; cette fonction d'adaptation peut
représenter la fonction objectif à optimiser. Ensuite, les meilleurs individus de cette
population sont sélectionnés, subissent des croisements et des mutations et une
nouvelle population de solutions est produite pour la génération suivante. Ce
processus se poursuit, génération après génération, jusqu'à ce que le critère d'arrêt
soit atteint, comme par exemple le nombre maximal de générations.

Les techniques de recherche et d'optimisation sont en général classées en trois

catégories (Coello Coello et al, 2002) : énumératives, déterministes et stochastiques.
Les AG font partie de la troisième catégorie et quatre caractéristiques les distinguent
des autres techniques d 'optimisation (Goldberg, 1989; 1994) [30] [31] :

ils utilisent un codage des paramètres et non les paramètres eux-mêmes.

ils travaillent sur une population d'individus (ou de solutions).
 81

ils n'utilisent que les valeurs de la fonction à optimiser, pas sa dérivée, ou une
autre connaissance auxiliaire.
ils utilisent des règles de transition probabilistes et non déterministes.

5.2.1 Principe de base d'un AG standard

Un AG standard nécessite en premier le codage de l 'ensemble des paramètres du

problème d'optimisation en une chaîne de longueur finie. Le principe d'un AG est
simple, il s'agit de simuler l'évolution d'une population d'individus jusqu'à un critère
d'arrêt. On commence par générer une population initiale d'individus (solutions).
Puis, à chaque génération, des individus sont sélectionnés pour la reproduction selon
le principe de la survie du plus adapté. Ensuite, des individus « enfants » (ou les
descendants) sont générés en appliquant les deux opérateurs génétiques suivants : le
croisement et la mutation. Ces enfants sont placés dans une nouvelle population P(t)
et vont se substituer, en tout ou en partie, à la population de la génération précédente.
De nouvelles populations d'individus vont ensuite se succéder, d'une génération (t) à
la génération (t+ 1), chaque génération représentant une itération jusqu'à l 'atteinte du
critère d'arrêt. La Figure 5-1 présente l'organigramme d'un AG standard [31].
 82

Début

Création de la population
initiale de la génération
t = 0

Évaluation de la fonction
d'adaptation de chaque individu
(définition de la fonction objective f)

Sélection des individus

(paramètres à optimiser)

Application des opérateurs de

croiseme nt et de muta tion
Non

G énération de l a nouvelle Critère d 'arrêt

population P (t) a tteint?
fnùn ?

Oui

Meilleur résulta t

Fin d e
l'optimisa tion

Figure 5-l Organigramme d'un AG standard

 83

5.2.2 Codage des variables dans un AG

Le codage est une partie très importante des algorithmes génétiques. Il permet de
représenter l'individu sous la forme d'un chromosome. Ce chromosome est constitué
de gènes qui prennent des valeurs dans un alphabet binaire ou non. Certains auteurs
n'hésitent pas à faire le parallèle avec la biologie et parlent de génotype en ce qui
concerne la représentation binaire d'un individu, et de phénotype pour ce qui est de sa
valeur réelle correspondante dans 1'espace de recherche [29].

Le choix du codage est délicat. Il doit permettre de coder toutes les solutions et
permettre la mise en œuvre des opérateurs de reproduction. C'est ainsi que le bon
déroulement des algorithmes génétiques sera assuré. Plusieurs type de codages sont
utilisés, on citera à titre d'exemple: codage réel, codage binaire, Gray [31].

Le codage binaire a pour principe de coder la solution selon une chaîne de bits qui
peuvent prendre les valeurs 0 ou 1. Ce type de codage est le plus utilisé. Il présente
l'avantage d' être facile et permet de coder toutes sortes d'objets : des réels, des
entiers, des chaînes de caractères, etc.

5.2.3 Opérateurs d'un algorithme génétique

Trois type d'opérateurs sont utilisés dans les AG; la sélection, le croisement et la
mutation. Et il existe plusieurs méthodes pour chacun des opérateurs [31] [32].
 84

5.2.3.1 Opérateur de sélection

La sélection a pour objectif d'identifier les individus qui doivent se reproduire sur
la base de la fonction d'adaptation, de telle sorte que les individus les mieux adaptés
sont sélectionnés alors que les moins adaptés sont écartés. Cet opérateur favorise
alors les meilleurs éléments selon le critère à optimiser (minimiser ou maximiser), qui
dans notre cas, est simplement définit par la valeur de la fonction objectif.

Cependant, il existe plusieurs méthodes de sélection, les plus connues sont : « la

roue de la fortune » ou par roue biaisée, et « la sélection par tournoi » (tournament
selection). Nous citerons aussi la sélection par la méthode « élitiste », et la sélection
stochastique.

5.2.3.2 Opérateur de croisement

Les opérateurs de croisement dans un AG sont conçus par analogie à la propriété

naturelle d'hybridation de 1'ADN. Cet opérateur combine les gènes de deux individus
parents pour donner deux nouveaux chromosomes d 'individus enfants (descendants)
possédants des caractéristiques issues des deux parents. La zone de croisement est
généralement choisie aléatoirement dans les chromosomes, et les méthodes de
croisement sont liées au codage mais leur principe est identique. Le croisement
permet d'enrichir la diversité de la population en manipulant la structure des
chromosomes (Chipperfield et al. 1994), comme il favorise aussi l'exploration de
l' ensemble des solutions possibles, il s ' agit alors de la phase d'intensification de la
recherche.
 85

De nombreux types de croisement existent dans la littérature. Ils préservent plus

ou moins l'identité génétique des parents et permettent un déplacement dans tout
l'espace des solutions. Les types de croisement les plus usuels sont : croisement en 1
point, en 2 points, et le croisement uniforme.

Le croisement utilisé avec les chaînes de bits est le croisement à découpages de

chromosomes qui se fait en trois étapes [31]. Au départ, deux solutions P 1 et P2
appelées parents sont sélectionnées en fonction de leurs performances. Puis, on tire
aléatoirement une position inter-gènes dans chacun des parents (point de croisement
Pt). Enfin, on échange les deux sous-chaînes de chacun des chromosomes, ce qm
produit deux autres solutions El et E2 appelées enfants (Figure 5-2).

Pour chaque couple, un nombre aléatoire P est tiré dans 1' intervalle [0, 1] et
comparé ensuite avec la probabilité de croisement Pc :

si P > Pc, le couple ne subit pas de croissement et un clonage de chromosome

aura lieu. Les deux enfants produits sont ainsi une copie de leurs parents.
si P < Pc, le croisement a lieu et un échange des parties des chromosomes des
parents va produire deux enfants par couple de parents.

Ce mécanisme présente l'inconvénient de privilégier les extrémités des individus

ce qui peut générer des fils plus ou moins proches de leurs parents. Pour éviter ce
problème, on peut étendre ce principe en découpant le chromosome en 3 sous-chaînes
comme illustré sur la Figure 5.3 ou encore en plusieurs sous chaînes.

Le croisement peut, aussi, se faire d'une manière uniforme. Chaque gène des
individus enfants est créé en copiant le gène correspondant d'un parent, choisi grâce à
un masque [32]. Pour chaque couple de parents, un masque binaire est généré
 86

aléatoirement. Comme illustré sur la Figure 5.4, pour chaque gène, lorsque le masque
a pour valeur 1, le gène du premier enfant est la copie du gène du premier parent.
Lorsqu'il a pour valeur 0, le gène du premier enfant est la copie du gène du second
parent. Cet opérateur de croisement favorise les recombinaisons complexes de gènes
entre les individus appariés.

Pt

lülülllllülllll•
1
1
1
1 ..
Croisement
El 1 1 1 ü 1 ü 1 1 1 ü 1 ü 1 1 1 ü 1
1
1
1
1
1

E2 1 1 1 1 1 ü 1 1 1 1 1 ü 1 1 1 1 1

Figure 5-2 Croisement en 1 point (codage binaire)

Ptl Pt2

Pl Il 1 ü i1
1
ü 1 1 11 f1
1
1
ü Il 1 1
Croisement
El Il lü lü 1 1 lü lü 1 1 1 1
1
1
1
1

'
'
1
1
1
1
1 1

P2ll Il lü II lü lü II lü 1 E2ll 1 1 lü 1 1 1 1 lü 1 1 lü 1

Figure 5-3 Croisement en 2 points (codage binaire)

Pl II 1 ü 1 ü 1 1 II 1 ü II 1 1 1 El 1 1 1 ü 1 ü Il 1 1 1 ü 1 1 1 ü 1
Croisement

E2 1 1 1 1 1 ü 1 1 1 ü 1 ü 1 1 1 1 1

~sque 1 1 1ü 1ü 1 1 1ü 1ü 1ü 1

Figure 5-4 Croisement uniforme (codage binaire)

 87

Le croisement est géré par le seul coefficient de croisement. Un tirage aléatoire est
réalisé et s'il est inférieur au coefficient de croisement alors un croisement a lieu
entre deux individus de la population parente. Classiquement, le coefficient de
croisement est compris entre 0, 7 et 0,9 pour un algorithme codé en réel. Dans le cas
contraire les enfants sont des copies des parents.

5.2.3.3 Opérateur de Mutation

Les mutations sont primordiales dans le processus de convergence vers 1'optimum

global. Elles permettent d'aller explorer l'espace des solutions de manière aléatoire.
Ainsi, on évite de rester piégé par un optimum local ; il s'agit de la phase de
diversification de la recherche. La mutation consiste à créer un individu fils a partir
d'un individu parent en faisant évoluer aléatoirement tout ou une partie des gênes de
ce dernier. Théoriquement, les mutations permettent à elles seules de faire converger
1' algorithme.

Elles sont gérées par le seul coefficient de mutation. Son principe d' application
est le même que celui du coefficient de croisement, ça consiste à modifier, avec une
probabilité Pm faible, certains bits des chromosomes. Nous tirons tout d 'abord pour
chaque bit un nombre aléatoire P dans 1'intervalle [0, 1]. Puis, nous le comparons avec
une probabilité de mutation Pm donnée :

Si P > Pm, le bit ne subira pas de modification.

Si P < Pm , la mutation est appliquée au bit correspondant.
 88

Comme pour les croisements, de nombreuses méthodes de mutation ont été

développées dans la littérature mais 1'une des plus efficaces est celle qui consiste à
muter chaque paramètre de la fonction à optimiser avec une probabilité dépendant des
informations contenues dans les gènes des individus. Dans le cas du codage binaire,
chaque bit est remplacé selon la probabilité Pm par son inverse. C'est ce qu'illustre la
Figure 5-5.

P 10 10 11 1 10 11 10 10 1 Mutation Elolol1 11 11o1o1

Figure 5-5 Représentation d'une mutation de bit dans une chaîne (codage
binaire)

5.3 Optimisation mono objectif

L' optimisation est dite mono objectif, à cause qu' on a une seule fonction
objective à optimiser. Où pour chaque individu de chaque génération, la fonction
d'évaluation est calculée puis utilisée plus tard pour choisir les individus qui vont
générer la population suivante.

Pour calculer le coût d'un point de l'espace de recherche, on utilise une fonction
d'évaluation ou d'adaptation F. L'évaluation d'un individu ne dépendant pas de celle
des autres individus, le résultat fournit par la fonction d'évaluation va permettre de
sélectionner ou de refuser un individu pour ne garder que les individus ayant le
meilleur coût en fonction de la population courante : c'est le rôle de la fonction (F).
Cette procédure permet de s'assurer que les individus performants seront conservés,
alors que les individus peu adaptés seront progressivement éliminés [33].
 89

La fonction d'adaptation, associe une valeur pour chaque individu. Cette valeur à
pour but d'évaluer le degré d'adaptation d'un individu à son environnement. Les
individus peuvent être aussi comparés entre eux. Cette fonction, propre au problème,
est souvent simple à formuler lorsqu'il y a peu de paramètres. Au contraire, lorsqu'il
y a beaucoup de paramètres ou lorsqu'ils sont corrélés, elle est plus difficile à définir.
Dans ce cas, la fonction devient une somme pondérée de plusieurs fonctions .

La fonction d'adaptation doit expnmer le plus fidèlement possible, la

problématique posée sous forme mathématique. Sa définition peut être simplement
analytique, ou elle peut éventuellement faire appel au jugement de l'utilisateur. En
raison de l'analogie avec la théorie de l'évolution (survie des individus les mieux
adaptés à leur environnement), les algorithmes génétiques sont naturellement
formulés en t erme de maximisation. Ils servent donc à déterminer le m aximum d'une
fonction F réelle à une ou plusieurs variables. Le problèm e d'opt imisation sur
l'espace de rechercheE est formulé comme suit:

maxF(x) (5.1 )
XEE

Si le problèm e à résoudre est un problème de minimisation d'une fonction f.

minf(x) (5.2)
x EE

Ceci équivaut au problèm e de m aximisation de F(x) , que l'on définit comme suit:

1 (5.3)
F(x) - - --
1+ [(x)
 90

5.4 Conclusion

Les algorithmes génétiques sont des méthodes méta heuristiques d'optimisation

globale basées sur des concepts de génétique et de sélection naturelle. Le composant
principal des AG est le gène qui se compose d'une chaîne de caractères (souvent
binaire). Les gènes s'enchaînent et forment les chromosomes. Ces derniers forment
les individus dans l'espace de recherche. Ainsi, les AG travaillent sur une population
d'individus, où chacun de ces derniers représente une solution possible pour le
problème donné. Dans chaque itération de l'AG, la performance de chaque individu
de la population courante est calculée. Les opérateurs de génétiques, sélection,
croisement et mutation, sont appliqués successivement pour créer une nouvelle
population jusqu'à l'approche rigoureuse de la solution optimale [30] [31].

Les algorithmes génétiques offrent plusieurs avantages :

- Ils ne demandent pas d' informations à priori ou des propriétés particulières de la

fonction objectif du problème.

- Leurs performances par rapport aux algorithmes classiques sont bien remarquées
lorsque par exemple les espaces de recherches sont importants [34].

- Outre leur facilité de programmation et de manipulation, ils sont facilement

adaptables à tout type de problème d'optimisation.

- Ils peuvent être utilisés avec profit pour traiter des problèmes n 'étant pas
optimisables efficacement par des approches purement mathématiquement.
 91

CHAPITRE6

TESTS ET V ALIDATI ON DE LA RÉGULATION (AVR, PID) ET DE

L'OPTIMISATION PAR ALGORITHME GÉNÉTIQUE DES PARAMÉTRES
D'EXCITATION DU GROUPE TURBINE-ALTERNATEUR

6.1 Introduction

Ce chapitre présente les différents tests et résultats d'application qui valide notre
travail. Les résultats des différents niveaux de régulations sont donnés (régulation de
tension par A VR, régulation de la vitesse de rotation, stabilisation de puissance par
PSS ... ), ainsi que les simulations de plusieurs tests sur le modèle du système global
assemblé et implémenté sous formes de blocs simulink et script matlab.m (voir
annexe); test d'amorçage des tensions, test de réponse à un saut de tension, ... etc. Les
résultats de l'application de la méthode d'optimisation par Algorithme Génétique
(sous matlab toolbox pour AG [35]) sont aussi présentés, et comparés aux cas avant
optimisation.

Les paramètres de chaque élément du groupe turbine-alternateur sont donnés en

tableaux (paramètres de la machine Lab-volts, paramètres de la turbine hydraulique
idéale, et paramètres du système d' excitation type STl).
 92

6.2 Paramètre des différents blocs

./ Paramètres de la génératrice (Lab-volt) :

Paramètres S!mboles Valeurs

Puissance apparente s 120 VA

Tension nominale au stator v 208 v

Courant maximal au stator I 0.33 A

Fréquence f 60 Hz

Vitesse de rotation N 120 Rpm

- Résistances (nu) :

Résistance au stator r, 0.0198

Résistance de fuite rr 0 .1983

- Réactances (J!u) :

X<f 0.350
Synchrones axe d et q
Xq 0.34 17

x'd 0.050
Transitoires axe d et q
x' q 0 .34 17

x"d 0 .0295
Sub-transitoires axe d et q
x" q 0.0295

Fuite au stator xr 0.0295

Kany xto -0.008 1

Tableau 1: Paramètres de la génératrice synchrone Labvolt

 93

./ Paramètres initiaux de l'excitatrice statique type ST1A:

Limitation des tensions VAmax= 300 V Amin= -300 VRmax = 6,4 VRmin = -6,4

Limitation des courants Yrmax= 4,4 Yrmin= -4,4 - -

Rectification Kc = 0.04 KLR = 1 - -

KA =1.62 Kf =0,0503 KE = 1 Tb = 1
Régulateur automatique de tension
(AYR)
TA =0,0087 Tf =0.3774 TE = 0 Tc = 1

Kp =1,13.10.4 Kr = 6,3496 Kct = 0,0082

Régulateur classique PID

Rt = 0,0159 L1= 0,006 M = 0,8.10.3

Paramètres du Transformateur de
stabilisation (en pu)

Tableau 2: Paramètres du syst ème d' excitation statique ST1

./ Paramètres de la turbine hydraulique (Idéale) :

Puissance nominale Pm 127.4MW

Vannage -gmax -Ouverture nominale : 5 pu
-gO -Ouverture à vide : 0.02 pu

Débit d'eau nominal q 85 m 3/s

Tableau 3 : Paramètres de la turbine hy draulique

 94

6.3 Régulation de la tension d'excitation

Un régulateur automatique de tension A VR est utilisé afin de contrôler et

stabiliser la tension d'excitation VFD, comme illustré par le schéma suivant :

~-------------------

1 Vers
laGS
1

Amplification Excitation 1

~------------------J
Excitation et AVR

Compensation

Figure 6-1 Régulation et stabilisation de l' excitation

En utilisant les param ètres initiaux du régulateur automatique de tension

A VR, on visualise la tension d'excitation comme suit :

1.5

ss
(/)
Cl.
E
ro
.r:
u
il)
-o
c 0. 5
0
(ii
c
2

0
0 2 3 4 5
Temps(s)

Figure 6-2 Tension d'excitation V FD

 95

6.4 Régulation de la tension de sortie de la génératrice

L'introduction d'un régulateur classique PID est faite dans le but d 'assurer une
meilleure régulation de la tension de sortie de la génératrice synchrone Vs, où cette
dernière est représentée par une fonction de transfert.

Système d'excitation

vref ... L verr ..... vs·...

~
+ KA VA KE Ka
.....
1 V FD
Régulation
;/ PID l+TA s
~ l+T0 s
J,- 1 - J l+TEs 1
1 1
Amplification Excitation Génératrice
1 v; 1
1
sKF
---
1
1 1
1 l+sTF 1
1 1
Compensation 1
1
-------------------------~

Vs, _!!__,___
l +J;.s
Capteur
Figure 6-3 Régulation de la tension génératrice Vs

Les variations des tensions de sortie VFD et Vs sont données comme suit :
T ension Génér atrice
1.6

1.4

1.2
s
s r"\.....
v
U)

>.,
-.;
<::
.E 0 .8
2<::
0 0.6
ën
<::
2
0 .4

0.2

0
0 2 3 4 5 6 7 8 9 10
Tem ps(s)

Figure 6-4 Tension génératrice Vs

 96

Tension d excitation

1.2

S'
§
>
§ 0.8
~
(::;
1il
0.6
1
"C
v
J
c:
0
·;;;
c:
~ 0.4
v
0.2

0
0 2 3 4 5 6 7 8 9 10
Temps(s)

Figure 6-5 Tension d'excitation correspondante VFD

./ Remarque :

Réponse concordante entre la tension génératrice Vs et la tension d'excitation VFD

(Figure 6-4 et 6-5).
La régulation par PID classique utilisée révèle quelques dépassements, une bonne
précision et une rapidité de réponse satisfaisante (Figure 6-4).
 97

6.5 Régulation des tensions dans le système global

Le modèle hybride de la machine synchrone est assemblé au système d'excitation

statique de type STl. Une représentation schématique est illustrée ci-dessous pour
montrer les différents blocs. Et les résultats de simulations sont donnés après
implémentation sous matlab simulink.

v{f) Génératrice
Synchrone
V:
r
vq
vh

r
Is I q
Ih
Modèle hybride

+:!
Système
d'excitation
statique type ST vfrej
(Modèle IEEE)
1- - - - - - - - - - - - -1
1 Régulation de la 1
1 1
L . . - - - - - --1 tension i-----{ V.ref
: génératrice :
L------------1

Figure 6-6 Représentation schématique du système global (Génératrice +

Excitation)
 98

Ten sion Génératrice

Tension Génératrice
1 .5 ,----,-----.----.-----.--.------r--I

::::J
30
.0

>"'
"'
!il
.!:
ê
~
<::
.2
(J)
<::

"'
1-

-1

-1.5 c___ _ r __ ____[_ ___[__ __[__ _L __ _ L __ ___,c___ _ r __ __L_----J:

0 2 3 4 5 6 7 8 9 10
Temps(s)

Figure 6-7 Tension de sortie de la génératrice Vsabc

Tensi on d excitation

0 .6

0.5 /
S'
c.
"="
>
c: 0 .4 1
1
0

~
·u
x 0 .3
Q.)
-.::::>
c:
0
·;;; 0 .2
c:
Q.)
1-
0 .1

-0 .1
0 2 3 4 5 6 7 8 9 10
Temps(s)

Figure 6-8 Tension d' excitation VFD correspondante

Les Figures 6-7 et 6-8 montrent un bon amorçage des tensions de sortie, et une
concordance entre la tension génératrice Vsabc et la tension d' excitation VFD·
 99

6.6 Tests et simulations sur le système global

6.6.1 Test d'amorçage des tensions de sortie

En déconnectant les courants d 'entrée Jsabc (courants à vide), les tensions VFD et
Vsabc sont visualisées et données comme suit :

Tension Génératrice

S'
su
..0

"'
>.,
""ffi
c:
1~
~
c:
.2
.,
"'
c:
1-

-1

2 3 4 5 6 7 8 9 10
Temps(s)

Figure 6-9 Tension de sortie de la génératrice Vsabc

T en s ion d excitat ion

Figure 6-10 Tension d'excitation VFD correspondante

 100

Les variation des puissances sont observées et affichées :

Puissance rotorique acti-.e

2r---c----r---,c---c----r---c----r---~---c--~

S'
sa: 1.5 1f

~ ~~~---r--~----+---~--~----~--~---+--~

Jo:v'L--~----L_--~__
0 2 3
_ L_ _ _ _L __ _

4 5
~ ~----~--~
6
_ __ L_ _

7 8 9 10
Temps(s)
Puissance rctorique réactive

S'
0.1
{\
-ê
a
<1.)
()
c
<Il
(/)
0.05
\
"
(/)
·:;
0.. 0

0 2 3 4 5 6 7 8 9 10
Temps(s)

Figure 6-11 Puissances au rotor

6.6.2 Tests - réponse à une variation de tension

Afin de démontrer la variation de la tension génératrice ams1 que sa tension

d'excitation correspondante, des tests de réponse à un échelon ont été effectués. Les
puissances relatives sont aussi visualisées .

./ Augmentation de tension (échelon croissant de 0.5 pu à 1 pu, à t=5s)

 101

Tension d excitation

S'
a. 0.6
;:::-
>
c
0
0.4 1
~
"" 0
à'l
-o 0.2
c
v
0
ïii
cQ)
1-
0
Î
0 2 3 4 5 6 7 8 9 10
Temps(s)

Tension Génératrice

2 3 4 5 6 7 8 9 10
Temps(s)

Figure 6-12 Tension d'excitation VFD et tension de sortie Vsahc

- Les puissances relatives sont présentées comme suit :

Puiss anc e rotorique acti-.e

S' 1.5
-ê
c.. 1
1
Q)
0
c
"'"'
"'
0.5
1
·::;
c..
0
....-
0 2 3 4 5 6 7 8 9 10
Temps(s)
Puiss ance rot orique réacti-.e
0.1
S'
5 0.05
{\
Q)
0
c
"'"' 1 \
"' ~

'
·::;
c.. 0

0 2 3 4 5 6 7 8 9 10
Temps(s)

Figure 6-13 Puissances au niveau du rotor

 102

Puissance statorique active

'3' -0.1
_e,
(/)
[]_
-0.2

\\.
Q)
u
c
!U
(/)
-0.3
(/)
·:;
[]_ -0.4

0 2 3 4 5 6 7 a 9 10
Temps(s)
Puissance statorique réactive

'3'
a.
ùl 0

"
a
Q)
u
c -0.5

'
!U
(/)
(/)
·:;
[]_ -1

0 2 3 4 5 6 7 a 9 10
Temps (s)

Figure 6-14 Puissances au niveau stator

./ Teste: Diminution de t ension : (échelon décroissant de 1 pu à 0.5 pu, à t =5s)

Tens ion d excit ation

'3'
a.
~ 0.6
>
c
0

~ 0.4
1 \
"()
~
0.2 1 \._
"tJ
1
§
(j)
c 0
Q)
1-
0 2 3 4 5 6 7 a 9 10
Temps(s)
Tens ion Génératrice
'3' 1.5
_e,
u
..c
!U
>
Q) 0.5
ro
-~ 0
E
2 -0.5
c
0
(j) -1
c
Q)
1-
2 3 4 5 6 7 a 9 10
Temps(s)

Figure 6-15 Tension d'excitation VFD et tension de sortie Vsabc

 103

Les puissances correspondantes sont affichées comme suit :

Puissance rctcrique active

~ 1.5
'!:::'
Il.
Q)
( 1\
u
ffiffJ 1 1\
ffJ
·s
Il.
0.5

0 v
0 2 3 4 5
' 6 7 8 9 10
Temps(s)
Puissance rctcrique réactive

:::l
o. 0.1 .....
'!:::'
a 0.05
1\
~
Q)
u
c
ro 0
l'..
ffJ
ffJ 1\ . /.,-
·s -0.05 l"
Il.
-0.1
0 2 3 4 5 6 7 8 9 10
Temps(s)

Figure 6-16 Puissances au niveau du rotor

Puissance statoriqu e active

~ -0.1 ~--~--~----~--~----~--~----+----+----+---~
:::J
s ffJ
Il. -0.2 l----~--~----~--~----~--~----+----+----+---~

1-0.3 f-'\
k-~~--~--~----~--~----~---+----+----+------+
·:;
[J_ -0.4 f-----~--~----~--~----+---~----+-----+----+------1-

0 2 3 4 5 6 7 8 9 10
Temps(s)
Puissance statorique réactive

S"
o.
v; 0

~
a Q)
u /
c
ro -0.5
~ 1
ffJ
.!!!
:::J
Il. -1

0 2 3 4 5 6 7 8 9 10
Temps(s)

Figure 6-17 Puissances au niveau du stator

 104

./ Remarques et discussions :

Les tests effectués montrent que la régulation de la tension génératrice

triphasée Vsahc présente une bonne performance en termes de précision et de stabilité
(Figure 6-9, 6-12 et 6-15).
Cependant on remarque un léger retard avant d'atteindre la valeur de référence
au niveau de la tension d'excitation VFD (Figure 6-10).
La parfaite concordance entre la tension génératrice et la tension d'excitation
est bien évidente.
Les puissances électriques au mveau du rotor et stator évoluent auss1 en
correspondance avec la régulation de la tension d'excitation ainsi que la tension
génératrice (Figure 6-11, 6-13, 6-14, 6-16 et 6-17).
 105

6. 7 Régulation de la vitesse de rotation et de la puissance mécanique

Afin de réguler la vitesse de rotation et la puissance mécanique, un PID classique

est choisi (voir chapitre 2), selon le schéma fonctionnel suivant :

Kp

Servomoteur r-----r--.t Vannage Turbine

Statisme

Figure 6-18 Régulateur de vitesse t ype PID

La m éthode de calcul des paramètres du régulateur PID (par compensation de

pôles) est donnée en (annexe 4 ).

Formulation de la vitesse :

w(t) = ~w(t ) + w0
1 t
{ ~w(t) = -
2H
j (T - T )dt- Kd ~w(t)
m e
0
 106

En utilisant les paramètres d'une turbine hydraulique idéale les plus usuelles
(Tableau 3), les résultats de simulation sont donnés sur les figures ci-dessous :

Vitesse de rotation

J\A
...
'§: 0.8
~
<1.)
(f)
0.6
(f)
<1.)
0.4
5
0.2

0
0 0.5 1.5 2 2.5 3 3.5 4 4.5 5
Temps(s)

Puissance Mécanique
2.5

2
S'
0..
Q) 1.5
()
c:
<Il
1\.
(f)
(f)
·s
0..
v-
0.5

0
0 0.5 1.5 2 2.5 3 3.5 4 4.5 5
Temps(s)

Figure 6-19 Vitesse de rotation de la turbine mr et puissance mécanique Pm

Selon la Figure 6-19, la vitesse de rotation ainsi que la puissance mécanique

suivent parfaitement la consigne grâce au régulateur PID, avec quelques
dépassements.
 107

L'ouverture du vannage g et le débit d'eau q correspondants sont comme suit:

Réglage du vannage

2 i - - Vannage g
- - Débitq
~

15 1.5
·<!.)
0
v<1.)
Ol
1\
<a
c
c ~-
<a
> If\
0.5 v

0
0 0.5 1.5 2 2.5 3 3.5 4 4.5 5
Temps(s)
Figure 6-20 Ouverture de vannage g et débit d 'eau q associés

L' erreur de vitesse est illustrée comme suit :

Erreur de vitesse
0.6

0 .4

0 .2

0 1\.
~ v
~
-0 .2

~ -0 .4
-0 .6

-0.8

-1
0 0.5 1.5 2 2.5 3 3 .5 4 4.5 5
Temps(s)
Figure 6-21 L'écart de vitesse co - co,.1

L'ouverture de vannage g et le débit d'eau q montre une bonne symétrie avec

la puissance et la vitesse de rotation (Figure 6-19 et 6-20). Et l'erreur de vitesse
s' annule assez vite (0.5s) Figure 6-21, ce qui justifie une répons e rapide.
 108

6.8 Stabilisation de puissance

Dans le but de remédier à l'amortissement des oscillations de puissances, et

d'avantage la stabilisation de la tension d'excitation, un PSS lA de type IEEE est
ajouté au niveau de la boucle de régulation des tensions (voir chapitre 2).

La régulation globale dans les différents blocs du système (groupe Turbine

hydraulique- Alternateur synchrone avec leurs système d'excitation STl), sont alors
présenté par le schéma suivant :

Excitation

Vannage

Machine
L-----~ Statisme Synchrone et
Turbine

Figure 6-22 Schéma de la régulation globale

 109

Le rôle du stabilisateur de puissance PSS est démontré dans les figures suivantes :

Tension d excitation

s ,, - - Pss Off

,r~\ v
8'+- 0.8 - - Pss On
>c
0 0.6
:;::::;
Cil
:t=
(.)
x 0.4
Q)

""0
c 0.2
.Q
(/)
c
Q) 0
f-

-0.2
0 0.5 1.5 2 2.5 3 3.5 4 4.5 5
Temps(s)
Figure 6-23 Variation de la tension d'excitation VFn avec et sans PSS

Puissance rotorique active

4
S' 3 1. ~ - - Pss Off

1~
0..

- - PssOn
't'
0...
2
1 '\.
Q)
(.)
c
(1l
(/)
.!!l
:::s li - """'
0
0...

-1
0 0.5 1. 5 2 2.5 3 3.5 4 4. 5 5
Temps(s)
Puissance rotorique réactive
0.6
- - Pss Off

"'
S'
8 0.4
a {' ~ - PssOn

-
Q)
(.) 0.2

'
c
(1l
Vl
.!!l 0 VI .....
:::s
0...

-0.2
'=~
0 0.5 1. 5 2 2.5 3 3.5 4 4. 5 5
Temps(s)

Figure 6-24 Variation des puissances au niveau du rotor avec et sans PS S

 110

ZOOM:

Tension d excitation
0.9
s ~ - PssOff

-
8 0.8
>
c
0
0.7
1"'\.
t "'\.\
- PssOn

~
(.)
x 0.6 1 ~
Q)
""0
c
1
\ / ---......
0.5
~
.Q
<J)
c
Q)
f- 0.4
/, '
0.5 1.5 2 2.5 3
Temps(s)

Figure 6-25 Tension VFD

Puissance rotorique active

s
,g, 3 1-----1--k-~---+-----+----f-------1 - - Pss Off
o._ - - Pss On
g
ro
2 ~-Jl-~~~-4-----+----~-~==~======~
Vl
-~
:::J
(l_

0.5 1. 5 2 2.5
Temps(s)
Puissance rotorique réactive
s Î'
,g, 0.4 1 - PssOff ~
~
"'
a [- Ps sOn
~ 0.2
c
ro
<J)
0 l# '\ __,_
' -
Vl
·s
(l_

0.5 1. 5 2 2. 5
Temps(s)

Figure 6-26 Puissances au rotor

 Ill

Sortie du Pss
0.4
::::J
o._
:ï5'
0.3
"'
û)
>
c
0
~
0.2
IJ)

:.0
"'
û) 0.1

"
(!)
"Cl
c
0
.in 0
c
(!)
1-
-0.1
0 0.5 1.5 2 2.5 3 3.5 4 4.5 5
Temps(s)

Figure 6-27 Tension de stabilisation Vstab

./ Remarques :

Le stabilisateur de puissance PSS permet de réduire les oscillations au niveau des

puissances (Figure 6-24, 6-26), tel qu'une baisse d'amplitude et gain en rapidité
au niveau de la tension d 'excitation (Figure 6-23, 6-25).
Le PSS améliore aussi la stabilité et la rapidité du système (6-27), grâce à la
tension de stabilisation Vstab qui intervient pour réduire les oscillations plus
rapidement.
 112

6.9 Optimisation des paramètres du système d'excitation

./ Optimisation par Algorithme Génétique [26]:

Afin d'optimiser les différents paramètres du système d'excitation statique type

ST utilisé dans notre travail, nous adoptons l'algorithme génétique connu pour sa
précision et sa robustesse dans ce type d'optimisation.

Après l'initialisation des paramètres de l'excitatrice (paramètres du régulateur

automatique de tension AVR, paramètres du PID), dans le but d'améliorer les
performances dynamiques du système, on applique l'indice d'optimisation ITSE
(Integral Time Square Error) définit par la fonction objectif f à minimiser, écrite
comme suit:

Où le carré de 1' erreur de tension e ( t) est intégré.

Application des opérateurs de sélection, croisement et mutation de l'algorithme

génétique dans matlab toolbox : [26] [35]
Sélection des paramètres : permet de spécifier le choix de sélection des
parents pour la nouvelle génération dans l'algorithme génétique
(SelectionFcn). Dans notre travailla méthode de sélection stochastique a été
utilisé (@selectionstochunif).
Croisement des paramètres : démontre la combinaison de deux individus
parents pour former les individus enfants pour la nouvelle génération dans
 113

l'algorithme génétique (CrossoverFcn). Le type de croisement utilisé est le cas

uniforme (@crossoverscattered.
Mutation des paramètres : spécifie le changement aléatoire du chromosome
d'un individu dans une population pour former l'individu enfant muté de la
nouvelle génération (MutationFcn). La mutation gaussienne est appliquée
dans notre validation (@mutationgaussian).

Les résultats obtenus après l'application de la méthode d'optimisation par

algorithme génétique sont donnés dans les tableaux et les figures ci-après.
 114

6.9.1 Optimisation des paramètres de l'AVR

Paramètres du Valeurs minimales Valeurs maximales Valeurs optimales

régulateur AVR (Algo Génétique)

KA 10 400 10

TA 0.01 0.1 0.1

KF 0.02 0.1 0.02

TF 0.3 1.5 1.5

Tableau 4: paramètres du régulateur automatique de tension AVR

Affichage des résultats :

Tension d excitation
1. 5
~ -- Avant optimisation 1

- - A prés optimisation
s
8
en
D..
E ~

'v/
ro
..r:::
()
Q)
-o
c
- ~ 0.5
c
2

0
0 2 4 6 8 10
Temps(s)

Figure 6-28 Tension d' excitation VFn - avant et après optimisation

 115

6.9.2 Optimisation des paramètres du régulateur PID

Paramètres du Valeurs minimales Valeurs maximales Valeurs optimales

régulateur PID (Algo Génétique)

Kp 0 100 1,0844
Kr 0 100 3,5137
Kn 0 100 0,0476

Tableau 5: paramètres du régulateur PID

Résultats obtenus :

Tens ion Génératri ce

1.6
1\ - Avant optimisation ~
1.4
- - A prés optimisation
~ 1.2
êiî
>
Q)
roc ' J'\-
E 0.8
2c 0.6
.Q
(/)

~ 0.4
.....,
0.2

0
0 2 3 4 5 6 7 8 9 10
Temps(s)

Figure 6-29 Tension Vs - avant et après optimisation

 116

./ Discussion :

- Les Figures 6-28 et 6-29 montrent clairement l'efficacité de l'optimisation par

Algorithme Génétique (AG), tel qu'on obtient de meilleurs résultats en
performance; rapidité de réponse (court temps de réponse), système stable
(élimination des dépassements), parfaite précision (suivi adéquat de consigne).

6.9.3 Tests et Simulation surie système global après optimisation

KA= 10 Kf= 0,02 KE = 1 Tb= 1

Régulateur automatique de tension
(AVR)
TA =0,1 Tf =1,5 TE = 0 Tc = 1

Kp =1.0844 K1 = 3,5137 Kct = 0,0476 -

Régulateur classique PID

Limitation des tensions VAmax= 300 VAmin= -300 VRmax = 6,4 VRmin = -6,4

Limitation des courants Vrmax= 4,4 Vrmin = -4,4 - -

Rectification Kc = 0.04 KLR = 1 - -

Tableau 6: paramétrés optimaux du système d'excitation STl

En utilisant les paramètres optimaux du système d'excitation statique type STl

IEEE obtenus par optimisation par AG (tableau 6), on a effectué quelques s imulations
sur le modèle global implémenté dans matlab simulink.
 117

6.9.3.1 Réponse à un échelon unitaire

Tensico Génératrice

Temps(s)

-1

2 3 4 5 6 7 8 9 10
Temps (s)

Figure 6-30 Tension génératrice Vsahc (après optimisation)

Tension d excitation
0.7

0.6
::;,
0..
f'\......
;;:::-' 0.5
>
c
0

~
·u
0.4

(i3
"0
0.3
c
0
·v; 0.2
cQ)
f-
0.1

o ~--_L _ _ _ _L __ __ L_ _ _ _L __ __ L_ _ _ _L __ __ L_ _ _ _L __ __ L_ _~

0 2 3 4 5 6 7 8 9 10
Temps(s)
Figure 6-31 Tension d'excitation VFD associée (après optimisation)
 118

./ Affichage des puissances :

Puissance rotorique acti~.e

2~---r---,----~--~----~---c----y----r----r---~

: :> ('....
~ 1.5 ~~~--~----~--~----+---~----+---~----~---+
Il.
Q)
u
c
l1l
fJ)

·s 0.5
fJ)

Il.
0
0 2 3 4 5 6 7 8 9 10
Temps(s)
Puissance rotorique réactive
0.4
S' A.
~ 0.3
aQ) 0.2
1\
u
c
l1l
fJ) 0.1
\
fJ)
·s
Il. 0
\
0 2 3 4 5 6 7 8 9 10
Temps(s)

Figure 6-32 Variation des puissances au rotor

Puissance statorique acti~.e

S" -0.1
o.
ù)
Il. -0. 2
Q)
u
c
l1l -0.3 1\
'
fJ)
fJ)
·s
Il. -0.4
0 2 3 4 5 6 7 8 9 10
Temps(s)
Puissance stat oriqu e réacti~.e

S'
sfJ) 0
a
Q)
u
c
l1l
-0.5
fJ)

-~
:::J
Il. -1
\
0 2 3 4 5 6 7 8 9 10
Temps(s)

Figure 6-33 Puissances au stator

 119

6.9.3.2 Réponse à une variation de tension

Afin de démontrer la variation de la tension génératrice ams1 que sa tension

d'excitation correspondante après 1' optimisation des paramètres de 1' excitatrice par
AG, un test réponse à un saut de tension est effectué. Les puissances relatives sont
aussi visualisées.
Augmentation de tension (échelon croissant de 0.5 pu à 1 pu, à t=0.5s)

Tension d excitation
:::J
o.
§; 0.6
5 If
~ 0.4
(3
>< ,......,
~ 0.2
c
0
-~ 0
Q)
1-
0 2 3 4 5 6 7 8 9 10
Temps(s)
Tension Génératrice
s 1. 5
s
_g
<0
>
Q) 0.5
(ij
-~ 0
§ -0.5
5
·v; -1
c
Q)
1-
2 3 4 5 6 7 8 9 10
Temps(s)

Figure 6-34 Tension d'excitation VFD et Tension génératrice Vsabc

 120

./ Affichage des puissances associées :

-
Puissance rotorique active

'§: 1.5
f
<1)
u
c
rn o.5
.!!2
::J If'
~ 0
0 2 3 4 5 6 7 8 9 10
Temps(s)
Puissance rotorique réactive

'§: 0.2
1:::"
0
~ 0.1
c
\
ro

~
VI
VI
0
r\ \
0 2 3 4 5 6 7 8 9 10
Temps(s)

Figure 6-35 Variation des puissances au rot or

Pui ssance st atorique act ive

s -0.1 1---- +-- -+- -+-- -+- --+- - +-- -+- -+-- -+-- ----4
BVI
~ -0 . 2 ~-+---+--+---+---+--+---+--+---+----+

i
·::;
-0.3 H
\t-----+-
. -+-- -+- -t--- +-- -+- -+-- -+- -----1- --+

~ -0.4 1---- +-- -+- -+-- -+- --+- -+-- -+- - +-- -+-- ----4
0 2 3 4 5 6 7 8 9 10
Temps(s)
Pui ssanc e s t at orique réac tive

s
B 0
VI
a \.
<1)
u
c -0.5
ro
\.
VI
VI
·::;
~ -1

0 2 3 4 5 6 7 8 9 10
Temps(s)

Figure 6-36 Variation des puissances au st ator

 121

./ Remarques et discussions des résultats :

La réponse optimale obtenue en utilisant les AG est très satisfaisante, en

comparant les résultats obtenus avant et après l'optimisation, dans tous les résultats
de régulation de la tension génératrice Vsabc et la tension d'excitation VFn associée sur
les Figures ( 6-7 et 6-8 1 6-30 et 6-31 ), on déduit :
Plus de rapidité (temps de réponse moindre) après l'optimisation.
Parfaite précision dans les réponses de Vsabc et VFD observées sur les Figures (6-
30 à 6-32).
Les résultats obtenus, après optimisation par AG montrent une stabilité meilleure
(élimination de dépassement).

6.10 Conclusion

C e chapitre illustre les résultats essentiels obtenus pour la validation de chaque

partie du travail visé; citant principalement la régulation et l' optimisation par
algorithme génétique des param ètres de l 'excitatrice après adoption du m odèle
st atique (ST). Nous pouvons alors, confirmer l 'efficacité et la fiabilité de la m éthode
d' opt imisation par AG dans l ' amélioration des performances dynamiques du système.
 122

CONCLUSION GÉNÉRALE

Ce projet d 'étude vtse particulièrement un élément primordial dans le bon

fonctionnement du simulateur analogique du groupe turbine-alternateur, que la
compagnie Hydro-Québec prévoit de construire, mais aussi un élément crucial dans
une centrale hydroélectrique en général. Le système d'excitation opère un rôle très
important dans un groupe-turbine-alternateur, entre assurance de la protection des
éléments du groupe et amélioration des performances dynamiques de l'ensemble.
Notre travail a visé cet important élément, traitant alors sa modélisation, et
1' optimisation de chaque paramètre qui intervient à ce niveau.

En premier lieu, nous avons abordé les motivations de ce travail d'un point de vue
socio-économique et éducatif tout en justifiant les objectifs et l'importance d'une
étude portée sur le système d'excitation dédié au simulateur analogique du groupe-
turbine-alternateur d'Hydro-Québec. Au second chapitre, l 'état de l'art sur les
systèmes d'excitations des génératrices synchrones est présenté, introduisant alors les
différents éléments constitutifs d'un système d'excitation, ainsi que les différents
types d'excitatrices pouvant intervenir dans un groupe turbine-alternateur. Le chapitre
3 de ce document, établie la partie modélisation de chaque bloc du simulateur
analogique, à savoir la génératrice synchrone, la turbine hydraulique et le système
d' excitation, le type et le modèle de ce dernier est élaboré par deux méthodes, avant
de faire une introduction aux régulateurs conventionnels qui interviennent au niveau
de 1'excitatrice statique de type ST, le fonctionnement de chaque partie de la
régulation est alors explicité dans le chapitre 4. En vue d'améliorer, ajuster et de
perfectionner les performances dynamique du système global défini après assemblage
des modèle des sous-ensembles du groupe turbine-alternateur, une méthode
 123

d'optimisation des paramètres du système d'excitation est choisie et présentée dans le

chapitre 5, cette méthode est basée sur les algorithmes génétiques (AG) appartenant à
la famille des algorithmes évolutionnistes qui sont de plus en plus utilisés de nos
jours dans divers domaine d'étude. L'adoption des AG à notre travail dans le but
d'optimiser les paramètres du système d'excitation de type statique ST du simulateur
analogique nous a permis de valider et d'atteindre les objectifs de notre travail.

Dans ce travail, plusieurs défis et objectifs ont été relevés, qui sont énumérés
comme suit:

Établissement d'un modèle d'état du système d'excitation statique par la

méthode DQ.

- Développement du modèle d'état de 1'excitatrice statique par méthode DQ avec

régulation de la tension d'excitation (chapitre 4).

- Développement du modèle du système d 'excitation statique de type STl avec

tous ses éléments, se basant au modèle de référence établi sur IEEE STlA.

Élaboration d'un modèle mathématique hybride des machines synchrones à

pôles saillants.

Application d'un travail de régulation de la tension par un régulateur

automatique AVR, au niveau du système d'excitation statique ST.

- Application d'un stabilisateur de puissance PSS au niveau de l ' excitatrice

statique STen vue d' élimination des oscillations de puissances.

- Mise en point d'un modèle global qui assemble les trois principaux modèles des
blocs essentiels du groupe turbine-alternateur sous forme de blocs simulink.

- Optimisation des différents paramètres du système d' excitation de type statique

en utilisant la méthode d'optimisation par algorithme génétiques.
 124

Le travail sur les systèmes d'excitations en général, implique encore plus de

recherche dans le domaine de la production d'énergie hydroélectrique, vue son
Immense importance dans 1' ensemble des composantes des centrales
hydroélectriques. Notre travail en particulier, effectué sur l'excitatrice du simulateur
analogique prévu chez Hydro-Québec présente une sérieuse perspective et un grand
acquis en vue de la conception expérimentale de cet élément.
 125

Références

[1] P.Kundur. "Power System Stability and Control" , New York: McGraw-Hill
Inc, 1994.

[2] Wei, and Zheng Xu, Member, IEEE, "Excitation System Parameters Setting
for Power System Planning". Power Engineering Society Summer Meeting
IEEE 2002.

[3] Mm.Gérard Cornu, Jean-Luc Fichet, Hervé Maurey et Rémy Pointereau.

"Rapport d'information n ° 463 (2013-2014)", fait au nom de la commission
du développement durable, Sénat le 16 avril 2014.

[4] Hydro-Quebec en bref, octobre 2014, (http: //www.hydroquebec.com/a-

propos-hydro-quebec/qui-sommes-nous/hydro-quebec-bref.html)

[5] Working Group on Prime and Energy Supply "Hydraulic turbine and turbine
control models for system dynamic studies ". 1992 IEEE tram. Power System.

[6] Haj agos, L.M. ; Basler, M.J. "Recommended practice for excitation system
models for power system stability studies". Changes to IEEE 421.5 Power
Engineering Societ y General Meeting 2005 IEEE.

[7] Dr.Dhiya Ali Al-Nimma, Dr. Majid Salim Matti, "Reactive Power Control of
an Alternator with Static Excitation System Connected to a Network", Mosul
University, Vol.l8, N°3, June 2010.

[8] M.S. Ghazizadeh, F. M. Hughs, "A Generator Transfer Function R egulator

for Jmproved Excitation Control", IEEE Trans on. Power Systems, Vo1.13,
N°2, May 1998, pp. 437- 441.

[9] P.M.Anderson, A.A Fouad. "Power System Control and Stability", IEEE
Press, 2003.

[10] Prepared by the IEEE Working Group on Computer Modelling of Excitation

Systems."Excitation system models for power system stability studies", IEEE
committee report, February 1981.

[11] Vladimir Chuvychin, Roman Petrichenko, Nikolay Gurov." Optimization of

Excitation System Parameters for Kegums Hydra Power Plant of Latvia".
 126

Environment and Electrical Engineering (EEEIC), 2011 , lOth International

Conference on.

[12] Q.Zhao.J.Jiang,"Robust Controller Design for Generator Excitation Systems",

IEEE Trans on Energy Conversion,Vol.lO, N°2, June 1995, pp.201-207.

[13] Technical Report IEC International Electrotechnical Commission, ref.34-16-2,

"Excitation systems for synchronous machines". First edition 1991-02.

[14] N.Srikanth (M.Tech), Atejasri-M.Tech Dept of IEEE, GIET, Rajahmundry,

A.P-India, "Enhancing power system stability by using thyristor controlled
series compensator", Vo1.2, Issue 5, September-October 2012, pp 1817-18-24.

[15] Bhumaiah Jula, M.Tech, Md.ASIF Assoc. Prof, Krishna Enugala, M.Tech.
"Simulation for protection of Huge Hydra Generator from Short Circuit
Faults", Vol.4, November 2012.

[16] J.Machowski and S.Robak, J.W.Bialek, J.R.Bumby, "Excitation control for

use with synchronous generators ". lEE Proc-Gener. Trans.distrub, Vo1.145,
N° 5 september 1998.

[17] C.H.Cheng, Y.Y.Hsu, "Damping of Generator Oscillation Using an

Adaptative Static Var Compensateur", IEEE Trans on. PowerSystems, Vol. 7,
no.2, May 1992, pp.718-724.

[18] F. Rotella, "Commande des systèmes linéaires à plusieurs entrées". Notes de

cours, Ecole Nationale d'Ingénieurs de Tarbes, 2002, rotella@enit.fr.

[19] K. Chaij arumudomrung, K-N. Areerak, and K-L. Areerak School of Electrical
Engineering, Suranaree University of Technology Nakhon Ratchasima,
Thailand. ''Modeling of Three-phase Controlled Rectifier using a DQ
Method'', 2010 International Conference on Advances in Energy Engineering.

[20] Koson.Chaijarumudomrung, Kongpan.Areerak, and Kongpo.Areerak. School

of Elect rical Engineering, Suranaree University of Technology Nakhon
Ratchasima,Thailand."The Controller Design of Three-Phase Controlled
R ectifier using an A rtificial Intelligence Technique" European Journal of
Scientific Research ISSN 2011.

[21] Koson.Chaijarumudomrung, Kongpan.Areerak and Kongpo.Areerak School of

Electrical Engineering, Suranaree University of Technology Nakhon
Ratchasima, Thailand. "The Stability Study of AC-DC Power System with
 127

Controlled Rectifier Including Effect of Voltage Control" European Journal of

Scientific Research ISSN 1450-216x vo1.62 no.4 (2011).

[22] Vadim Slenduhhov, Jako Kilter Tallinn University of Technology (Estonia)

v.slenduhhov@gmail.com, jako.kilter@ttu.ee "Modeling and Analysis of the
Synchronous Generators Excitation Systems", pamu 2013, p 252-257.

[23] The Institute of Electrical and Electronics Engineers, Inc, New York, USA.
Copyright © 1992. "IEEE Recommended Practice for Excitation System
Modelsfor Power System Stability Studies" .

[24] T.Zabaiou, Member IEEE, Louis-A Dessaint, senior member IEEE, Patrice
Brunelle. "Development of new library of IEEE Excitation system and its
validation with PSS/E'. 2012 IEEE.

[25] Anant Oonsivilai and Padej Pao-La-Or. "Application ofAdaptive Tabu Search
for Optimum PID Controller tuning AVR System". Institute of Engineering,
Suranaree University of Technology, vol. 3, june 2008.

[26] Cao Yuan, Ma Jin, "Research on PID Parameters Optimization of

Synchronous Generator Excitation Control System", Critical infrastructure
(CRIS), 2010 5th International Conference. 20-22 Sept. 2010.

[27] Working Group on Prime Mover and Energy Supply Models f or System
Dynamic Performance Studies "Hydraulic Turbine and turbine control models
for system dynamic studies", Transactions on Power Systems, vo1.7, no.1,
February 1992.

[28] Meng Shen, Student Member, IEEE, Vaithiananthan Venkatasubramanian,

Member, IEEE, Nicholas Abi-Samra, and Dejan Sobajic, "A New Framework
for Estimation of Generator Dynamic Parameters". IEEE transactions on
power systems, vol.l5, no.2, may 2000.

[29] E.Lutton, "Darwinisme artificiel", INRIA Rocquencourt Equipe Complex

Proj et Fractales, 2004.

[30] J.Philippe Rennard, "Genetic Algorithm Viewer: Démonstration d'un

algorithme génétique", www.rennad. org/alif, Avril 2000.

[3 1] D. E. Goldberg, "Genetic algorithms in search, optimisation, and machine

learning". London: Addison-Wesley, 1989.
 128

[32] L.Davis, Handbook of Genetic Algorithms.Van Nostrand Reinhold New

York, 1991.

[33] Z.Michalewicz, GeneticAlgorithms. Springer-Verlag: New York, 1992.

[34] Sidhartha Panda, N. P. Padhy, "Comparison afP article Swarm Optimization

and Genetic Algorithmfor TCSC-based Controller Design". World Academy
of Science, Engineering and Technology 3 2007.

[35] A.J.Chipperfield and P.J.Fleming. "The matlab genetic algorithm toolbox".

Department of Automatic Control and Systems Engineering, University of
Sheffield. From lEE colloquium on applied control techniques using
MATLAB, Digest no.1995/014.
 129

Annexe 1

Conversion des grandeurs en (pu)

valeur de la grandeur actuelle

va1eur d e 1a gran deur en pu = -------=--------
valeur de la base de la grandeur

Stator Circuit de champ Unité

Tension
usrtf = .[ivN u = ssréf
M
fréf 1
fréf

Courant J fréf [A]

J sréf = ; i JN

Puissance 3 [VA]
s sf"(jf =2Uro/ rof =3VNJN

Impédance s [0]
z = usréf u fréf
sréf 1 z fréf = --=----z
sréf I fréf I

Inductance [H]
L = z sréf L = z fréf
sréf fréf
(J)réf (J)réf

Flux [Wb)
= 1 = u sréf
l.f/sréf L sréf sréf
(J)réf
 130

Annexe2

./ Linéarisation autour d'un point de fonctionnement ou d'équilibre:

Soit une fonction y= f(x, z) non linéaire, définie autour du point d'équilibre par

On a la série de Taylor :

Tronquée après les termes du premier ordre;

5y=y-ye
5x=x-xe
{
5z= z - ze

La loi linéarisée est alors définie par :

1) Le modèle linéaire de la turbine hydraulique dans le chapitre 3 est obtenu selon

cette méthode, voici dessous, le développement des équations :
 131

A partir du modèle non linéaire suivant (3.35) :

Selon le développement de Taylor;

On aura les deux fonctions:

dq
dt = J;.(q,g)
{
PM = fz(q,g)

(En régime stationnaire ( q 0 , g 0 ) , la fonction dq est nulle), on obtient enfin :

dt
 132

;o
3

=> Avec Pa= 2

4 go

2) Obtention du système linéarisé du modèle DQ du convertisseur statique

(chapitre 3) :

Soit, en général, un système non linéaire définit par les fonctions d'état et de sortie
suivantes :

Xn_=fn (X (t), ...,X (t),~ (t), ... ,Um(t)) => {X(t): f(X(t), U(t))
1 1

~ - g 1 (X1 (t), ... ,X1 (t) , U1 (t), ... , Um (t)) Y(t) - g(X(t ), U(t ))

La linéarisation du modèle d'état autour du point d' équilibre (X0 ,U 0 ) est

développée selon Taylor au 1er ordre des fonctions f et g :

 133

Équations d'état

xl R!J;(Xo,Uo) + !;!lxu (Xl-Xo)+ .. +:; 1 (Xn-Xo)+ :D 1 (Ul-Uo)+ .. +a~ 1 (Um-Uo)

O· o n Xo.Uo 1 Xa.Uo m Xo.Uo

Équations de sortie

En utilisant les relations suivantes :

5X =X - X 0
5U = U - U 0
{
5Y = Y - Y0

Les équations peuvent alors s'écrire :

- Équations d'états

oX1 ~ J;(X ,U )+ oX
0 0 ê!f; 1 oX1 + ... + oJ;
oXn
1 oXn + ê!f;
oU1
1 oU1 + ... + ê!f;
oUm
1 oUm
1 Xo .Uo Xo .Uo Xo .Uo Xo .Uo
 134

-Équations de sortie :

La forme matricielle est :

c5X(t):A(X0 ,U0 )c5X +B(X0 ,U0 )c5U

{ c5Y(t) -C(X ,U )c5X +D(X ,U )c5U
0 0 0 0

Avec :

% of;. % of;.
oX1 0Xn oU1 oUm
A(Xo,Uo) = -of-
r 1 ;B(Xo,Uo) = lfT 1

oX Xo ,Uo ou Xo,Uo
of. lfn of. lfn
oX1 oXn Xo,Uo
oU1 oUm Xo Po

og1 og1 og1 og1

oX1 0Xn oU1 oUm
C(X o, Uo) = -og-
T 1 .,D(Xo,Uo) -- -og-
T 1 =
oX XoPo oU XoPo
ogp ogp ogp ogp
oX1 oXn XoPo
oU1 oUm Xo P o
 135

Annexe3

Paramètres utilisés dans les simulations du modèle du convertisseur statique par la

méthode DQ:

Sans régulation de la tension Vf [19] :

Req = 0.150; Leq =30pH; Ceq = 2nF;

[CF= lOOOpF; rF= 0.030; LF = 8mH;

Avec régulation de la sortie Vf [20] :

Re = 0.150; Leq = 30 pH; Ceq = 2nF;

[C: = 8000 pF; rF = 0.030; LF = 8mH;

Annexe 4

Régulation de la vitesse de rotation : Exemple de calcul des paramètres du

régulateur PID ;

./ Méthode par compensation de pôles :

 136

1
=> GBF(s) = 1
1+-s
Ka

On déduit alors les paramètres du régulateur :

1
T1 = -K =>K1 = K a
a

Td =Ta => Kd = Ta
K = 1 tel que Tc : constante de temps
p Tc .Ka .RP

Régulation de la puissance mécanique :

A partir du modèle d'état de la turbine hydraulique équation (3.38)

dq 2 2~
-=---q+--g
dt TWGO TWGO

dg =-K g+K u (a)

dt a a

Pm =3q-2~g

./ méthode par retour d' état:

p(x)=Ax + Bu
A partir de la formulation d'un système d'état : (b)
{ y=Cx+Du
 137

Et par superpositions des deux systèmes (a) et (b), les matrices d'états sont:

C = 3 et D = 2~ où la commande représente le

vannage g, le débit q et la puissance mécanique P rn les états.

Le calcul de la matrice d'état K c se fait comme suite :

La loi de commande : uc = -K_ cx+v(t)

Prenant un pôle désiré p =1, on aura ~ = S -1

2
Sachant que P(
A+BK_c
) =S+A+BK c =S---+
-
2
J: G
~
J: G
K c
-
w 0 w 0

Pour Pd = P(A + BK_c ) =>

=> K
_c
= _!_(
~
1- TwGo)
2
 138

Annexe 5

Code matlab du modèle par DQ du convertisseur statique non commandé

% Dé claration d e s équa ti o ns :

fun c t i on px fcn(x , params , U)

%#eml
ids = x(1) ; iqs = x(2) ; vbusd = x(3) ; vbusq x ( 4) ; i de x ( 5) ;
vf = x ( 6) ;

pcpl = U(1) ; vsd = U(2) ; vsq = U(3) ;

Req = params(1) ; Leq = params(2) ; Ce q = params(3) ; r mu =

params(4) ;
rf params(5) ; Lf = params(6) ; Cf= params(7) ; w = params(S) ;

px= zeros ( 6 ,1);

px(1 ) -( Req/ Leq ) *i ds + w* i qs - (1 /Leq)*vbusd + (1/Leq ) *vsd ;

px(2 ) - w* i ds - (Req/ Leq ) * i qs - (1/Leq )*vb u sq + (1/Leq ) *vsq ;
px(3) (1/Ceq )*ids + w*vbusq - (3*sqrt ( 2 )/(pi *Ceq ) )*ide ;
px(4) - w*vbusd + (1/Ceq ) *iqs ;
px(5) ( 3*sqrt (2)/(pi * Lf) ) *vbu sd - ( (rf+rmu)/Lf)*idc -
(1/Lf) *vf ;
px(6) = (1/ Cf ) *idc - (1/ Cf ) * (pcp l/vf);

% initia li sa ti o n de l ' espace de travail

cl ear a ll ; c l ose a ll ;
cl c
gl oba l pa r ams U

%paramètres de l ' excitatrice

% Donn ées
f 60 ;
Req 0 . 15 ;
Leq = 30e - 6
Ceq = 2e- 9;
w = 2*pi *f;
rmu = (3*w* Leq )/pi
rf = 0 . 03 ;
 139

Lf = 8e- 3 ;
Cf = 1 0000 e - 6 ;
landha O = p i/ 3 ;
a l phaO 15*pi /180 ;
params = [Req Leq Ceq rmu rf Lf Cf w];

Vs 23 0 ;

Vm sqrt( 2 )*Vs ;
vsd sqrt(3/2)*cos( l andha0 + a lphaO)*Vm;
vsq = sqrt (3 /2)*sin(landha0 + a lphaO)*Vm;

% conditions initiales :
XO = [14.25 ; 0 ; 130 . 5 ; 3 7 4 . 2 ; 10 . 55 ; 180 . 55] ; %système non
commandé

% Affichage des résultats

p l o t(t , Vf , ' lineWidth ' , 2) ; h o l d on

g rid on
xlab el( ' Temps(s) ' ) ;
y l abe l( ' Te n sion Vf ' ) ;
titl e ( ' Tens i o n Excitatrice ' )
axis ([0 .4 0 . 9 11 0 1 50 ])

% ----------------------------- FIN ---------------------------- %

 140

Annexe 6

Code matlab du modèle par DQ du convertisseur statique avec régulation de

l'excitation:

% Déclaration des équat i ons :

fun c t i o n px = fcn(x , params , U, cont ro l)

%#eml

ids = x( 1 ) ; iqs = x(2) ; vbusd x( 3 ) ; vbusq x ( 4) ; ide x ( 5) ;

vf x ( 6) ;
Xv = x( 7 ) ; Xi x(8) ;

pcpl = U(1) ; vsd = U(2) ; vsq U(3) ; vf r ef u (4) ;

Req = params(1) ; Leq = params(2) ; Ceq = params(3) ; r mu =

params (4) ;
rf= params(5 ); Lf = params ( 6 ); Cf= params(7 ) ; w = params (B) ;

kpv = cont ro l(1); kpi = contro l( 2 ); kiv = cont ro l( 3 ); kii =

contro l(4);

px= zeros (8 ,1) ;

px (1) - (Req/ Leq ) *i ds + w* i qs - (1 / Leq ) *vbusd + (1 /Leq ) *vsd;

px ( 2 ) - w* ids - (Req/ Leq ) * i qs - (1 /Leq ) *vbu sq + (1 / Leq ) *vsq ;
px(3 ) (1/Ceq ) *ids + w*vbusq - (3*sqrt ( 2 ) /(pi*Ceq ) ) *ide;
px(4 ) - w*vbusd + (1/Ceq ) * i qs ;
px(5 ) - ( (rf+kp i) / Lf ) * i dc + (( kiv *kp i) / Lf) *Xv + (kii/ Lf)*Xi -
(( kpv*kpi)/Lf)*vf + ( ( kpv*kpi ) /Lf)*vf ref ;
px(6 ) (1 /Cf)*idc - (1 /Cf)* (pcpl /f) ;
px(7 ) - vf + vf r ef ;
px(B ) - i de - kpv * vf + kiv*Xv + kpv* vf ref
 141

% paramètres de convertisseur commandé :

% Données
f 60 ; % fréquence
Req 0 .1 5 ;
Leq = 30e- 6
Ceq = 2 e - 9 ;
w = 2*p i *f ; % pulsati on p ropre
rmu = (3*w*Leq)/pi % r ési s t a nce v a riabl e AC- DC
rf 0 . 03;
Lf = 8e- 3 ;
Cf = 8 00 0e- 6 ;
landhaO = pi/3 ; % Déphasage
a l phaO p i /18 ; % angle d ' a morçage des thyr istors
params = [Re q Leq Ceq r mu rf Lf Cf w] ;

% c o ntrôleur PI
vf ref = 300 ;
e = 0 .8;
wni 2 *pi*50 ;
wnv = 2 *pi*30 ;

% contr ôle de tensi o n

kpv 2* e *wnv*Cf ;
kiv = (wnvA2 ) *Cf ;

% cont rô l e de cou ra nt
kp i = 2*e*wni*Lf - rf ;
k ii = (wn i A2 ) *Lf ;
contr o l = [ kpv kp i k i v ki i] ;

Vs = 230 ;
Vm = sqrt (2 ) *Vs ;
vsd sqrt ( 3 / 2 ) *cos (la n dha0 + a l p h aO ) *Vm;
v sq = sqrt( 3/2 ) *s in(la n dh a0 + a l p h aO ) *Vm;

% conditions i nitia l es :
XO = [14. 25 ; 0 ; 1 30 . 5 ; 37 4. 2 ; 1 0 . 55 ; 300 ; 0 . 009 ; 2 . 7837e- 4 ]; %
sys t ème comma ndé .

% Affi c h age des rés ult a t s

p l ot ( t , Vf , ' li n e Width ' , 2 ) ; h o l d on
grid on
x l abe l( ' Te mps(s ) ' );
y l abe l( ' Te n s i o n Vf( p u) ' );
tit l e ( ' Te n sion Excitatr ice ' );
ax i s ( [ 0 .1 0 . 5 296 3 0 2 ] );

% ----------------------------- FIN ----------------------------%

 142

Annexe 7

Dessin simulink du système d'excitation type STl- Régulateur de tension

Vt.I'Rmaxllc.lf

Vf

Saturnt10n
O)llamic1

Amort1,.ment

Annexe 8

Dessin simulink du Stabilisateur de puissance PSSlA

 143

Annexe 9

Dessin simulink du régulateur de vitesse de la turbine hydraulique

Scope 1

1---1.-..... ŒJ
q

l--~--.t·werr Wc g ate l----'-------.----1~

Scop e 2
PID control Bloc Servomoteur

Goto

Ga1n
M odel_turb _hyd CD
Pm
L-~---------~ R p ~--------~

Ga tn3

L---------------<é -K-~--------------,

B
From 1

Add Gain 2
w
 144

Annexe 10

Dessin simulink du modèle hybride de la machine synchrone

Ph i _S

Got o

______. Phi_r

r
Phi_s

Phi_s
Vs 1
Vs

IS IS

Fl ux Statori que Ph1 _ s TenSon S: atori que_ Vs

,.....,.....,. l s
L:.J IS Phi_r
ls

Phi_r Ir
-CD
Ir

,...,....... Vr

~
Vr
IS
Vr

Fl ux Roto n ue P h1 r
Co urant rot o nqu e_l r
 145

Présentation des blocs séparément:

Phi_r

ls Produ ct1

Calcul du flux statorique Phi_s

Ph i_r
Ph i_r
Product2
Phi_r

IS

ls Prod uct

Sc ope 10

Yr

Vr Prod uct1

Calcul du flux rotorique Phi_r

 146

1/VvTl

ls
2
ls

Calcul de la tension statorique Vs

Phi_r

Ph1 r Product

ls

ls Product1

Calcul du courant rotorique Ir

 147

hela
Id
Constan t lntegrator Id

la
la
Park lq

lb lq

lb

lh
le
lh
le
abc ••• dq

functi o n [Id , Iq , Ih] = Park(theta , I a , I b , Ic)

%#eml
P = (2/3) *[ cos(theta) , cos(theta - (2 *p i) /3) ,
cos (theta+(2*p i)/3) ; ...
- sin(theta) , - sin(theta- (2*p i)/ 3 ) , -
sin (t h eta+(2*pi)/3 ); ...
0 . 5 , 0 . 5 , 0 . 5 ];

Id P( l,:) * [I a I b I c ] ' ;
Iq P (2 ,:) * [ Ia Ib Ic ]';
Ih P (3 ,:) * [ Ia Ib Ic ]';
e nd

Bloc transformée de park

 148

E}--+[I}-+ heta
Va
ŒJ
Va
Conslant1 lntegrator1

ŒJ • Vd
Vd lnvPark Vb ... CD
Vb
CD • Vq
Vq
Vc Œ)
CD
Vh
• Vh
Vc
dq --abc

function [Va , Vb , Vc ] = i nv Park (theta , Vd , Vq , Vh)

i nv P = [cos(the ta) , - si n( theta) , 1; ...

cos(the t a - (2*pi) /3 ) , - s in(theta - (2 *p i )/ 3 ) , 1 ; ...
cos (t h eta+( 2*p i) /3 ) , - s i n( t h e ta+ ( 2*p i ) /3 ) , 1 ] ;
Va i nv_ P (1 , :) * [Vd Vq Vh] ' ;
Vb inv_ P (2 , :) * [Vd Vq Vh ] ' ;
Vc inv_ P (3 , :) * [Vd Vq Vh] ' ;

e nd

Blo c transformée de Park inverse

 149

Annexe 11

Dessin simulink du modèle global du groupe turbine-alternateur

Ç9-------+L]
C tock T o Wotk~ ~~,ce1

T ra n s~in11p a rk
T ra n s_p a rk1

êj
G ot o14

~.....--------.
Attri b uti ons

T o Wo rk~ a ce4
~ Fro m4
1~'
Scope5
~ V•I.ab
G ro u nd
C a lcul d u cou p l e Te

{VT]

G o to3
M od ul e
V ref
Exc1l e t

M odul e2

P art ie mé ca n i q u e _ T u rb i ne e l Gov e rn eur j.-- - - + -( w_W ref

G o lo 1

To IN:J r~ a ce1 5
 150

Annexe 12

Scripte matlab.m du programme global

% initiali sation de l ' espace de tra v ail

%=======================================

cl e a r all ; c l o s e all
cl c

% Do nnées

% Caracté ristiques d e b a se de la machine

s base 120 ; % p u i ssance d e ba s e (VA)

V base 208 *sqr t (2 ) ; % te ns i o n de b as e (V)
I b a se 0 . 3 3*s q rt( 2 ); % c o urant de b a s e (A )
Z b a se V b a se/I bas e ; % impédance de b ase (Ohms )
f 60 ; % fr é qu e nce d e b ase (he rtz)
wn 2* pi * f ; % p ul sation é l ectr i que de base
p 2; % Nombre de pa i res d e pô l e
N 60*f/p ; % v itesse de r o t at i o n
wm 1; % vitesse de synchronisme

% Paramètres de la machine (en pu )

ra 0 . 0 1 98 ; % rés i s t a n ce strat or i q u e
rD 0 . 012 ; % résistance de l ' amortisseur D
rQ 0 . 0 07 3 ; % résistance de l ' amortisseur Q
rh 1; %r a ; % r és i s t a n ce h omopo l a ire
rf 0 .1 983 ; % r és i s t a n ce de f uit e
xa 0 . 0 1 0 1; % réactance stratori q ue
xD 0 . 669 ; % réactance de l ' a mortisseur D
xQ 0 .1 352 ; % réacta n ce de l ' a mo rti sseur Q
xh 1; %x a ; % réacta n ce homopo l a ire
xd 0 . 35 ; % réactance mutue l le à l ' axe d
xq 0 . 3 4 17 ; % réacta n ce mu tue ll e à l ' axe q
xk f 1 - 0 . 00 8 1; % réacta n ce de ka n y
xf 0 . 02 95 ; % réactance de fuite

% réactance strator i que axe d

xmd = x d - xa ;
% réactance stratori que axe q
xmq = x q - xa ;
% réac t ances mu t u e ll es
x fd xmd ;
x Dd xfd ;
xqQ xmq ;
 151

Suite du scripte matlab (page 2)

xfD xmd + x kfl ;

xDf xfD ;
xff x md + xkfl + xf ;
xDD x md + xkfl + xD ;
xQQ xmq + xQ ;

% Matrice des résistances du stator

r s = di ag ([ra ra rh]) ;
% Matrice des résistances du rotor
r r = diag ([rf rD rQ ]) ;
% Matrice globale des résistances du stator et rotor
R = blkdiag( - r s , r_r) ;

% Matrice des réactances du stator

Xss = diag ( [xd xq xh]) ;
% Matrice des réactances mutuelles entre le stator et le rotor
Xsr = [xfd xDd 0 ; 0 0 xqQ ; O 0 0] ;
% Matrice des réactances du rotor
Xr r = [xff xfD O; xDf xDD 0 ; 0 0 xQQ ];
% Matrice g lobal e des réactances
X= [ - Xss Xsr ; - Xsr ' Xrr ];

% calcul des ma trices du modé l e :

% Flux r o t o ri q u e
A_ phi = - wn*r_ r*inv(Xrr );
B_ phil - wn*r_ r*inv(Xrr ) *Xsr ' ;
B_ phi 2 = wn *eye (3 );

%Flux stato rique

C Xsr*inv (Xrr );
D = Xsr*inv (Xrr ) *Xsr '- Xss ;

% t e n s i on s t atoriqu e :
J = [ 0 - 1 0 ; 1 0 0 ; 0 0 0 ];
% Courant a u rotor
Al i nv (Xrr) ;
A2 = i nv(Xrr)*Xsr ' ;

% ---------------------- Partie Mécanique -----------------------

% Paramètres r e liés à l a turbine et governeu r:

At= 1. 8 ; % gain de la turbi ne

Tw = 0 .1; % temps de démarrage de l ' eau
beta = 0 . 001; % es t e de la turbine
 152

Suite du scripte (page 3)

Ka = 30; % este gain statique du servomoteur

Ta= 0 . 1 ; % este de temps du servomoteur
gmin = 0 . 0 1 ; % ouve r ture minimale du vannage
gmax = 1 0 ; % ouverture maximale du vannage
Vgmin = 0 . 0 1 ; % vitesse d ' ouverture minimale
Vgmax = 1 0 ; % vitesse d ' ouverture maximale
Rp 0.8; % 0 . 02 statisme
gO = 0 . 02; % vannage initial
GO = At* gO ;
H 0 . 05 ; % constant d ' inertie
d = 1;

% Régulation de la vitesse

Wr e f = 1 ; % vitesse de référence
Pr e f = 0 . 6 ; % puissance de référence
Pe = - 0 . 35 ; % puissance électrique
d r ef = 1 ; % mode de fonctionnement du régulateur
wO = 0 . 1;
% par a mètres des rég u l ateur :
% cal cu l des par ams PID
tau c = 0 . 01;
Kp_w 1 /(tau_ c*Rp*Ka) ;
Ki w Ka ;
Kd w Ta ;

% Gain d u r e t o ur d ' é t a t
K_ c = (1 /At)*( 1 - (Tw*G0 /2)) ;

% %- -------------------- Partie Exc i tat i o n-----------------------

%par a mè tres de l ' exc ita tri ce :

% Réf érence artic l e " Exc i tat i o n s y s p a r am setting fo r p o we r

sys t em p l a ning " [2 ].

% Do nnées par amèt res de l ' exc itatri ce ST1

fe = 60 ; % fréq ue nce
Te = 1 /fe ; % péri od e
Te _pu = Te*2*p i;
 153

Suite du scripte (page 4)

% params transformateur de stabilisation (pu)

Rt = 0 . 0159 ; % resistance
L1 = 0 . 6e- 2 ; % leakage inductance
M = O. Be - 3 ; % mutual inductance

% angle d ' amorçage des thyristors

a l pha = pi/ 6 ;

% tensio ns de commande :
Vs rms = 208 ;
V_ base = Vs_ rms*sqrt( 2);
Vab = sqrt(3)* Vs_ rms ;
Vab_pu = Vab/V_base ;
VAmax 300 ;
VAmin - VAmax ;
Vl max 4.4;
Vl mi n - Vlmax ;
VRmax 6. 4;
VRmin - VRmax ;

% facteur de rectification .
Kc = 0 . 04;
KLR = 1;

Vref = 1;
Vstab = 0 ;

% Gains et este de temps;

Tc = 1;
Tb = 1;
Kexc 1;
Texc = 0;

% Facteur de transformation
Ke = 12 ; % Ratio excitation/transfo
Ktb = 1 0 ; % Ratio transfo/GS

% Amplification
KA= (1. 35*Ke )/Ktb
TA = Te pu/ 1 2 ;
K f M/Rt;
T f = L1/ Rt ;
 154

Suite du scripte (page 5)

% FT Génératrice
xmd1 = 0 . 33 99 ;
Rf 0 .1 983 ;
xf 0 . 0 2 95 ;

Kg xmd1 /Rf ;
Tg xf/Rf ;

% Capteur
Kr 1;
Tr = 6e - 3 ;

% Correcteur PID par compensation de pôles:

e = 0 . 9; %coefficient d amortissement
w n = 1; % pulsation propre du système
Kp ( (TA+Tg ) *Tr *w n A2)/(Kr *KA*Kg) ;
Ki 1/(TA+Tg) ; ~
Kd (TA*Tg )/(TA+Tg ) ;

% par a mètres d u Stabil isateur de pu issance -------------------- %

Ks t ab = 10;
T w
~
0. 2;
T1 n 1;
T1d 1;
T2 n 0;
T2 d 0;
%- -------------------------- FIN ------------------------------ %
 155

Partie Optimisation par AG :

% Optimisation par algor ithme g énétique

% Optimisation des paramètres de l' AVR

par am [KA TA K f T_f ];

f_obj optim_ for Ef d(pa r am ); % fonction Objective

% Optimisation AG

opt i ons =
gaoptimset( 1 0 utputFcn 1 , @opt i m_ f or Ef d , 1 SelectionFcn 1 , @se l ect i o nst o
chu n if , .. .
1
CrossoverFcn 1 , @crossoverscatt ered , . . .
1
MutationFcn 1 , @mutati ongaus s i a n) ;
l b = [ 1 0 0 . 0 1 0 . 02 0 . 3 ] ;
ub = [ 400 0 . 1 0 . 1 1. 5 ] ;
x= ga(@optim_ f or Efd , 4 , [ ], [], [ ], [ ],lb , u b );

% fonction ob j ective

function f obj = opt im_ Efd_ Vs (parametre )

1
ass i g n i n( b a se 1 ,
1
KA 1 , paramet r e (1) );
1
ass i g n i n( b a se 1 ,
1
TA 1 , paramet r e (2 ) );
ass i gn in( 1
base 1 ,
1
K_ f 1 , parametre (3 ) ) ;
ass i gn in( 1
base 1 ,
1
T_ f 1 , parametre (4) ) ;

[to ut , xout , yout ] s im ( 1 Con t r o l Efd 1 ) ;

f obj=O ;
for i = 1: 50 0 1
f ob j = f obj + t ou t (i , 1) * (yout(i , 1) - 1) A2 ;
end
 156

% Optimisa t ion des paramètres PID

para= [Kp,Ki,Kd];
f Ob = optim_pid(pa r a ); % f o n c ti on Objective

% Optimisation par algor i thme g énétique

options =
gaop timset( ' OutputFcn ' ,@ optim_pid , ' SelectionFcn ' ,@ s electio n stochun
if , ...
' CrossoverFcn ' , @crossoverscatte red, ...
' MutationFcn ' , @mutationgaussian );
lb = z eros ( 3, 1) ';
ub = [10 0 100 1 00 ];
x= ga(@op tim_pid,3, [], [], [], [],lb ,ub );

% f o nction objective

funct i on f Ob = opti m_pi d(para)

ass ignin( ' base ' , ' Kp ' , p a r a (1) );

ass ignin( ' base ' , ' Ki ' , p a r a (2 ) );
a ss i g n i n( ' base ' , ' Kd ' , para( 3 ) );

[to ut , xout , yout] s i m( ' Cont ro l _Vs _ op t ' );

fOb =O ;
f o r i =1: 5 00 1
fOb = f Ob + t o ut( i, 1) * (yout(i , 1 )-1) A2 ;
end

%---------------------------- FI N ------------------------------- %
 157

Blocs simulink associés au scrypte matlab (partie optimisation par AG) :

~
Clock T o Wo~ace

Efd f-------i~---+\Tî
Efd

Vq

1-----j~·vr_ref scope 1

vr_refl ST1

Ef d t-----t~---+!
Efd1

Vq 1

t------t~Vf_ref Scope2

Vf_ref2 S T2

Bloc simulink - Optimisation des paramètres de 1'A VR

 158

con sant2

Add

Control de tension Efd avant optimisation (bloc ST)

Constant2

Add

Amorti ~me nt

Control de la tension Efd après optimisation (bloc STl)

 159

,--------______.~ ~
Clock ToWork::;:lace1
ToWo rk:;:lace2

~--------------~ CJ
f-----~Vstab Sco pe

Con9:a nt1 " (') f-----<.--•Œ.J

VI ( ~

0. 01 f-----~VfO ,-----------~CJ
Con3ant2 Sco p e1

Eld f----+~{

f-----I~Vs_rel Efd

Step

Con tro l_ VS

ToWork::;:lace4

r---------------.18
f--- -- li>IV""" Scope2

Con3ant3 VI(•Jf----l.--~3)
Vl{ ~1

00 1 t-----ll>lvto ,--------------li>ICJ
Con 3 a nl4 Scope3

Efd f--- -+-+{

f-- -- li>lvs_rel Efd 1

St ep1

Controi _V81

ToWor~ace 5

Bloc simulink - Optimisation des param ètres du PID

 160

Vf

Mlorti~ ert

_ Gain1
1
Capte !.X
Kr
Trs-1

Régulation de la tension Vs avant optimisation

_ G ain1
1
Capteur
VI

Régulation de la tension Vs après optimisation