Qualite de L'energie Electrique
Qualite de L'energie Electrique
Qualite de L'energie Electrique
i
Rsum :
Nous prsentons dans ce travail la problmatique de la Qualit de lEnergie Electrique, les
phnomnes principaux qui la dgradent, leurs origines ainsi que leurs consquences sur les
quipements et lappareillage lectriques. Pour rpondre aux besoins de notre tude, nous
avons eu recours un ensemble doutils et de mthodes dont la Mthode du Vecteur
dEspace qui sinscrit dans lanalyse des perturbations lectriques. Enfin, nous avons procd
la ralisation dune application EPQA (Electrical Power Quality Analysis) ddie la
surveillance et lanalyse de la qualit de lnergie lectrique.
Mots Cls: Qualit de lEnergie Electrique, Mthode du Vecteur dEspace, Application,
surveillance.
Abstract:
In this work, we have presented the electrical power quality issues, the main phenomena
which degrade it, their origins and their consequences on electrical equipments. For our study
we have resort to the whole tools and methods, which the space vector method is a new
method for analysis of the electrical power disturbances. Finely we have proceeded to the
development of an application dedicated to the monitoring and analysis of the electrical
power issues.
Key words: Power Quality, Vector Space Method, software, monitoring
ii
Table ues matieies
Table des matires ............................................................................................................................. ii
Introduction Gnrale .................................................................................................................... viii
Chapitre I ............................................................................................. Perturbation de la tension
1. Introduction .................................................................................................................................. 3
1.1 Classification des perturbations lectriques............................................................4
1.2 Catgories des phnomnes lectromagntiques....................................................6
2. Creux de tension et coupure brves ...................................................................................... 7
2.1 Dfinition......................................................................................................................7
2.2 Origines........................................................................................................................8
2.3 Consquences...............................................................................................................9
2.4 Analyse des creux de tension.....................................................................................10
2.4.1 Signatuie uun cieux ue tension..............................................................................................10
2.4.2 Piopagation sui le iseau........................................................................................................11
2.4.S Piopagation tiaveis les tiansfoimateuis.............................................................................11
2.4.4 Paiameties uteiminant le type uu cieux ue tension............................................................12
2.4.S Classification ues cieux ue tension..........................................................................................13
3. Surtension ................................................................................................................................... 15
3.1 Dfinition.....................................................................................................................15
3.2 Effets des perturbations.............................................................................................16
4. Variation et fluctuation de la tension................................................................................ 17
4.1 Dfinition.....................................................................................................................17
4.2 Effets des fluctuations................................................................................................18
5. Dsquilibre ............................................................................................................................... 18
5.1 Dfinition.....................................................................................................................18
5.2 Effets de dsquilibre.................................................................................................19
6. Harmoniques et interharmoniques .................................................................................. 19
6.1 Dfinition.....................................................................................................................19
6.2 Effets des harmoniques..............................................................................................20
7. Conclusion ................................................................................................................................... 21
Chapitre II ................................................ Mthode d'analyse des perturbations lectriques
1. Introduction ............................................................................................................................... 22
2. Analyse des creux de tension, coupures brves et surtensions ................................ 22
2.1 Estimation de lamplitude..........................................................................................22
2.1.1 valeui efficace...........................................................................................................................22
2.1.2 Signal analytique.......................................................................................................................23
2.1.S Composante fonuamentale......................................................................................................23
2.2 Segmentation..............................................................................................................23
2.3 Mthodes de classification.........................................................................................24
2.4 Caractrisation............................................................................................................26
3. Harmoniques .............................................................................................................................. 27
3.1 Outils danalyse...........................................................................................................27
3.1.1 Srie de Fourier......................................................................................................................27
3.1.2 Transformation de Fourier...................................................................................................28
3.2 Caractrisation............................................................................................................28
4. Dsquilibre ............................................................................................................................... 29
iii
5. Conclusion ................................................................................................................................... 29
Chapitre III .................................................................................... Mthode du Vecteur d'Espace
Mthode du Vecteur dEspace ....................................................................................................... 28
1. Introduction ............................................................................................................................... 30
2. Transformation du Vecteur dEspace ................................................................................ 30
3. Caractristiques et reprsentation du Vecteur dEspace ........................................... 32
3.1 Perturbations au niveau du fondamental................................................................32
3.2 Perturbations harmoniques......................................................................................33
4. La transformation du Vecteur dEspace pour lanalyse des perturbations
lectriques ........................................................................................................................................... 34
4.1 Vecteur dEspace et composante homopolaire en cas de creux de tension..........34
4.1.1 Cieux ue tension monophas...................................................................................................34
4.1.2 Cieux ue tension biphas.........................................................................................................39
4.1.S Cieux ue tension tiiphas........................................................................................................42
4.1.4 Cieux ue tension avec suitension............................................................................................43
4.2 Analyse des creux de tension par la transformation du Vecteur dEspace...........45
4.2.1 Extiaction ues tiaits caiactiistiques ues tensions mesuies..............................................46
4.2.2 Segmentation............................................................................................................................47
4.2.S Classification.............................................................................................................................50
4.2.4 Caiactiisation..........................................................................................................................53
5. Conclusion ................................................................................................................................... 53
Application ddie la surveillance en temps rel .............................................................. 55
Chapitre IV ............................................ Application ddie la surveillance en temps rel
1. Introduction ............................................................................................................................... 55
2. Organigramme gnral de lapplication ......................................................................... 57
2.1 Paramtres de la tension...........................................................................................58
2.1.1 Algoiithme uu piogiamme......................................................................................................58
2.1.2 Nouele utilis............................................................................................................................59
2.2 Charge :........................................................................................................................59
2.2.1 Algoiithme uu piogiamme......................................................................................................59
2.2.2 Nouele ue isolution ues quations uiffientielles...............................................................60
2.3 Redresseur :................................................................................................................60
2.4 Acquisition :.................................................................................................................61
2.5 Analyse :.......................................................................................................................62
3. Prsentation des menus et des fentres de lapplication (EPQA): ........................... 63
3.1 Menu :...........................................................................................................................63
3.2 Fentre principal.........................................................................................................63
3.3 Paramtres de la tension :...........................................................................................65
3.4 Charge:.........................................................................................................................66
3.5 Redresseur:.................................................................................................................67
3.6 Acquisition:.................................................................................................................68
4. Validation de lapplication .................................................................................................... 70
4.1 Redresseur diodes :.................................................................................................70
4.2 Acquisition..................................................................................................................76
5. Exemples dapplications : ...................................................................................................... 77
6. Conclusion ................................................................................................................................... 86
iv
Conclusion Gnrale ........................................................................................................................ 87
Bibliographie : ................................................................................................................................... 88
Annexe .................................................................................................................................................. 89
v
Liste des Figures :
Figure II. 1 Algorithme de classification creux de tension bas sur la mthode des six tensions...25
Figure II. 2 Algorithme de classification des creux de tension bas sur les composantes
symtriques...................................................................................................................................26
vi
Figure IV. 17 Valeurs obtenues par Simplorer et EPQA....................................................................71
Figuie Iv. 18 Couiant ue souice uonn pai Simploiei....................................................................72
Figure IV. 19 Courant de source donn par EPQA..............................................................................72
Figure IV. 20 Redresseur diodes D1, D5, D6 dconnectes..............................................................73
Figure IV. 21 Tension et courant redresss donns par Simplorer....................................................73
Figure IV. 22 Tension et courant redresss donns par EPQA...........................................................73
Figure IV. 23 Valeurs obtenues par Simplorer et EPQA.....................................................................74
Figure IV. 24 Courant de source de la phase 1 donn par Simplorer................................................74
Figure IV. 25 Courant de source de la phase 1 donn EPQA.............................................................74
Figure IV. 26 Courant de source de la phase 2 donn par Simplorer..............................................75
Figure IV. 27 Courant de source de la phase 2 donn par EPQA.......................................................75
Figure IV. 28 Courant de source de la phase 3 donn par Simplorer................................................75
Figure IV. 29 Courant de source de la phase 3 donn par EPQA.......................................................76
Figure IV. 30 Courant dune source mesur par EPQA ..................................................................76
Figure IV. 31 Courant de sources mesur par DAQVIEW...............................................................77
Figure IV. 32 Cration de la tension Exemple1...................................................................................77
Figure IV. 33 Rsultats de la simulation Exemple : 01.......................................................................78
Figure IV. 34 Charge linaire Exemple 1............................................................................................78
Figure IV. 35 La reprsentation de Fresnel Exemple 1......................................................................79
Figure IV. 36 Vecteur dEspace Exemple 1..........................................................................................79
Figure IV. 37 Analyse du courant de source Exemple 1.....................................................................80
Figure IV. 38 les courants de source dune charge RLC quilibre exemple 1................................80
Figure IV. 39 Charge non linaire exemple 2......................................................................................81
Figure IV. 40 Rsultats de simulation Exemple 2...............................................................................81
Figure IV. 41 Analyse des courants de source Exemple 2...................................................................82
Figure IV. 42 Vecteur despace Exemple 2..........................................................................................82
Figure IV. 43 Reprsentation de Fresnel Exemple 2...........................................................................83
Figure IV. 44 Analyse du courant de source Exemple 3.....................................................................83
Figure IV. 45 Vecteur dEspace (systme de courant) Exemple 3......................................................84
Figure IV. 46 Paramtres du dsquilibre (courant) Exemple 3.......................................................84
Figure IV. 47 Reprsentation de Fresnel (courant) Exemple 3..........................................................85
Figure IV. 48 Analyse des tensions de source Exemple 3....................................................................85
Figure IV. 49 Vecteur (tension) despace exemple 3...........................................................................86
Figure IV. 50 Reprsentation de Fresnel (tension) Exemple 3...........................................................86
vii
Liste des Tableaux :
Introduction
Gnrale
Intiouuction uniale
ENP 2u12 1
ENP 2u12 2
Chapitre I
Perturbation de la tension
Chapitie I Peituibation ue la tension
ENP 2u12 S
1. Introduction
La Qualit de lnergie lectrique se rapporte une varit de phnomnes
lectromagntiques qui caractrisent la tension et le courant un moment donn et un
emplacement donn sur le rseau, ltude de cette dernire revient ltude de la qualit
de la tension et du courant. Dans ce premier chapitre nous allons prsenter la qualit de la
tension, ses caractristiques, les phnomnes qui la dgradent, leurs origines, et leurs
consquences sur les quipements et appareillages lectriques.
Dans la distribution de lnergie lectrique, la tension constitue le support de la
puissance et le courant reprsente lnergie relle qui transite dun point un autre. Un
systme de tension triphas est caractris par quatre grandeurs. Toute variation enregistre
de ces grandeurs constitue une perturbation lectrique. Ces grandeurs peuvent tre rsumes
de la manire suivante :
Amplitude : Cest gnralement la premire exigence impose au distributeur
dnergie. La variation tolre pour cette grandeur est de 10% de sa valeur nominale
pendant un temps de 10 ms quelques minutes.
Frquence : Dans le cas idal, les trois tensions sont alternatives et sinusodales dune
frquence constante de 50 ou 60 Hz selon le pays. La norme EN 50160 fixe 1% la
variation de la frquence pendant 95% dune semaine et [-6%, +4%] en cas de
perturbations importantes. Les marges de variations de la frquence sont plus importantes
pour les rseaux autonomes.
Forme donde : la forme idale des trois tensions est une sinusode, mais on considre
en gnral que le signal est compos dun fondamental de frquence nominale auquel se
superpose des signaux diffrentes frquences appels harmoniques, la tension peut
contenir des signaux non priodiques appels bruits.
Equilibre : la symtrie impose aux trois tensions davoir la mme amplitude et le
mme dphasage relatif qui est de
2
3
pour le systme triphas.
Chapitie I Peituibation ue la tension
ENP 2u12 4
ENP 2u12 S
g
i
m
e
t
a
b
l
i
Coupures brves
<0.1 pu Courts-circuits Arrts dquipements,
pertes de production
Surtensions
1.1 1.8
pu
Courts-circuits,
dbranchement des
charges
importantes
Dclenchements,
dangers
pour les personnes et
pour
les matriels
Dsquilibre
-
Charges asymtriques
ou
monophases
Echauffements des
machines tournantes,
vibrations
Variations rapides de
tension (Flicker)
0.1 7 %
Charges fluctuantes
(fours
arc, moteur dmarrage
frquent, soudeuses,
oliennes)
Papillotements
de lclairage
Harmoniques
0 20% Charges non linaires
(structures
dlectronique de
puissance, arcs
lectriques)
Echauffements,
vieillissements, pertes
supplmentaires
Inter-harmoniques
0 2%
Charges non linaires et
fluctuantes (fours arc,
soudeuses, oliennes)
Papillotements
de lclairage
Bruit
0 1% Fours arc, charges non
linaires
Echauffements, pertes,
vieillissements
Variations de la
frquence
-
Dsquilibre entre la
production et la
consommation
Dysfonctionnements des
quipements lectriques
Tableau I. 1 Vue densemble des principales perturbations lectriques [2].
Chapitie I Peituibation ue la tension
ENP 2u12 6
ENP 2u12 7
5.1 DC offset
5.2 harmoniques
5.3 Inter-harmoniques
5.4 Bruit
0-
100
me
harmonique
bande large
tat d'quilibre
tat d'quilibre
tat d'quilibre
tat d'quilibre
0 0,1 %
0 20 %
0 2 %
0,1 %
6. Fluctuation de la tension < 25 kHz Intermittent 0,1- 7 %
7. Variation de frquence <10 sec
Tableau I. 2 Classification des phnomnes lectromagntiques [3].
2. Creux de tension et coupure brves
2.1 Dfinition
Un creux de tension figure(I.1) est une baisse brutale de la tension en un point
dun rseau, une valeur comprise entre 1% et 90% (CEI 61000-2- 1, CENELEC EN
50160), ou entre 90% et 10% (IEEE 1159) de la tension nominale suivie dun rtablissement
de la tension aprs un court laps de temps compris entre la demi-priode fondamentale du
rseau et une minute, un creux de tension est caractris par :
Sa profondeur (U),
Sa dure (t) dfinie comme la dure pendant laquelle la tension est infrieure
90 %. On parle de creux de tension x% si la valeur efficace passe en dessous de
x% de la valeur sa valeur nominale.
ENP 2u12 8
de dure suprieure. Les coupures brves et les coupures longues sont diffrentes tant du
point de vue de lorigine que des solutions mettre en uvre pour sen prserver ou pour en
rduire le nombre.
Dans le cas dun systme triphas, les paramtres U et t sont en gnral diffrents
sur les trois phases. Cest la raison pour laquelle un creux de tension doit tre dtect et
caractris sparment sur chacune des phases. Un systme triphas est considr comme
subissant un creux de tension si au moins une phase est affecte par cette perturbation [4].
2.2 Origines
Les creux de tension et les coupures brves sont principalement causs par des
phnomnes conduisant des courants levs qui provoquent travers les impdances des
lments du rseau une chute de tension damplitude dautant plus faible que le point
dobservation est lectriquement loign de la source de la perturbation. Les creux de tension
et les coupures brves ont diffrentes causes :
Des dfauts sur le rseau de transport (HT), de distribution (BT et MT) ou sur
linstallation elle-mme.
Lapparition des dfauts provoque des creux de tension pour tous les utilisateurs. La
dure dun creux est en gnral conditionne par les temporisations de fonctionnement des
organes de protection. Lisolement des dfauts par les dispositifs de protections (disjoncteurs,
fusibles) provoquent des coupures (brves ou longues) pour les utilisateurs aliments par la
section o le dfaut est prsent. Bien que la source dalimentation ait disparue, la tension du
rseau peut tre entretenue par la tension rsiduelle restitue par les moteurs asynchrones ou
synchrones en cours de ralentissement (pendant 0,3 1 s) ou la tension due la dcharge des
condensateurs branchs sur le rseau. Les coupures brves sont souvent le rsultat du
fonctionnement des automatismes de rseau tels que les renclencheurs rapides ou lents, les
permutations de transformateurs ou de lignes. Les utilisateurs subissent une succession de
creux de tension ou de coupures brves lors de dfauts arc intermittents, de cycles de
dclenchement-renclenchement automatiques permettant llimination des dfauts fugitifs
ou encore en cas de renvois de tension permettant la localisation du dfaut.
la commutation de charges de puissance importante (moteurs asynchrones, fours arc,
machines souder, ) par rapport la puissance de court-circuit.
Les coupures longues sont le rsultat de lisolement dfinitif dun dfaut permanent
par les dispositifs de protection ou de louverture volontaire ou intempestive dun appareil.
Les creux de tension ou coupures se propagent vers les niveaux de tension infrieurs travers
Chapitie I Peituibation ue la tension
ENP 2u12 9
les transformateurs. Le nombre de phases affectes ainsi que la svrit de ces creux de
tension dpendent du type de dfaut et du couplage du transformateur. Le nombre de creux
de tension et de coupures est plus lev dans les rseaux ariens soumis aux intempries que
dans les rseaux souterrains. Mais un dpart souterrain issu du mme jeu de barres que des
dparts ariens ou mixtes subira aussi des creux de tension dus aux dfauts affectant les
lignes ariennes.
Les transitoires de dure infrieure la demi-priode fondamentale sont causs, par
exemple, par la mise sous tension de condensateurs ou lisolement dun dfaut par un fusible
ou par un disjoncteur rapide BT.
ENP 2u12 1u
[
p
.
u
]
Va
Vb
Vc
0.5
1
30
210
60
240
90
270
120
300
150
330
180 0
Biagiamme ue Fiesnel
Chapitie I Peituibation ue la tension
ENP 2u12 11
)
(Z
S
+ Z
)
E (I. 1)
V
A
tension mesure en point A.
En avale le creux de tension nest pas attenue. En effet si le court-circuit se produit
avant le point A, on aura VA = VD = 0.
ENP 2u12 12
B
a
A
C
b
c
Dy
E
a
E
b
E
c
I
a
I
b
I
c
V
a
V
b
V
c
I
n
Chapitie I Peituibation ue la tension
ENP 2u12 1S
3) Endroit du dfaut
Il influe sur la propagation du creux de tension (aval et amont).
4) Type de mesure
Deux types de mesures sont utiliss :
Mesures entre phase et neutre : elles donnent gnralement plus dinformations
Mesures entre phases : elles permettent de mesurer les perturbations telles quelles
sont subies par les charges.
2.4.5 Classification des creux de tension
Selon la classification universelle reconnue de [6], on distingue sept principaux
types de creux de tension, on les rsume dans le tableau ci-dessous :
E
a
E
b
E
c
I
a
I
b
I
c
V
a
V
b
V
c
Chapitie I Peituibation ue la tension
ENP 2u12 14
ENP 2u12 1S
3. Surtension
3.1 Dfinition
Toute tension applique un appareil lectrique qui dpasse la crte de la tension
normale de fonctionnement de cet appareil est appele surtension. Les surtensions sont de
trois natures :
temporaires ( frquence industrielle),
de manuvre,
dorigine atmosphrique.
Les surtensions temporaires :
Elles sont la mme frquence que la frquence du rseau. Elles ont plusieurs
origines :
Un dfaut disolement :
Lors dun dfaut disolement entre une phase et la terre dans un rseau neutre
impdant ou isol, la tension des phases saines par rapport la terre peut atteindre la tension
compose. Des surtensions sur les installations BT peuvent provenir des installations HT par
lintermdiaire de la prise de terre du poste HT/BT.
La ferro-rsonance :
Il sagit dun phnomne oscillatoire non linaire rare, souvent dangereux pour les
quipements, se produisant dans un circuit comportant un condensateur et une inductance
saturable.
La rupture du conducteur de neutre :
Les appareils aliments par la phase la moins charge voient leur tension augmenter
(parfois jusqu la tension compose).
La surcompensation de lnergie ractive :
Les condensateurs shunt produisent une augmentation de la tension depuis la source
jusquau point o ils se trouvent. Cette tension est particulirement leve en priode de
faibles charges.
Les surtensions de manuvre :
Elles sont provoques par des modifications rapides de la structure du rseau
(ouverture dappareils de protection).
Chapitie I Peituibation ue la tension
ENP 2u12 16
Figure I. 8 Surtension.
3.2 Effets des perturbations
Leurs consquences sont trs diverses selon le temps dapplication, la rptitivit,
lamplitude, le mode (commun ou diffrentiel), la raideur du front de monte et la frquence:
Claquage dilectrique : cause de destruction de matriel sensible (composants
lectroniques),
dgradation de matriel par vieillissement,
coupure longue entrane par la destruction de matriel (perte de facturation
pour les distributeurs, de production pour les industriels),
perturbations des circuits de contrle commande et de communication courant
faible,
contraintes lectrodynamiques et thermiques (incendie) causes par :
La foudre essentiellement,
Les surtensions de manuvre.
0 20 40 60 80 100 120 140
-300
-200
-100
0
100
200
300
400
Temps en [ms]
T
e
n
s
i
o
n
e
n
[
V
]
Chapitie I Peituibation ue la tension
ENP 2u12 17
ENP 2u12 18
=
|u
1
|
|u
1 d
|
ct u
0
=
|u
1 0
|
|u
1 d
|
(I. S)
Dans les conditions normales dexploitation, la norme EN 50-160 impose que sur
chaque priode dune semaine, 95 % des valeurs efficaces moyennes sur 10 minutes de la
composante inverse de la tension dalimentation ne dpassent pas 2 % de la composante
directe.
ENP 2u12 19
ENP 2u12 2u
y(t) =
0
+
k
n
k=1
2sin(2nk
0
t +
k
) (I. 6)
c]]
= _(
k
)
2
k=0
(I. 7)
Le taux de distorsion harmonique (THD) donne une mesure de la dformation du signal
IE = _(_
1
]
2
k=2
(I. 8)
La norme EN 50-160 impose que dans des conditions normales dexploitation, le
taux global de distorsion ne dpasse pas les 8%, et pendant chaque priode dune semaine,
95 % des valeurs moyennes sur 10 mn des taux individuels ne dpassent pas les valeurs
indiques dans le tableau suivant :
Harmoniques non multiple de 3 Harmoniques multiple de 3 Harmoniques paire
Rang Seuil (%) Rang Seuil (%) Rang Seuil (%)
5 6 3 5 2 2
7 5 9 15 4 1
11 3.5 1521 0.5 624 0.5
13 3
17 2
19,23,25 0.5
6.2 Effets des harmoniques
Leurs consquences sont lies laugmentation des valeurs crtes (claquage
dilectrique) et efficaces (chauffement supplmentaire) et au spectre en frquence
(vibration et fatigue mcanique) des tensions et des courants. Leurs effets ont toujours un
impact conomique du fait du surcot li :
une dgradation du rendement nergtique de linstallation,
un surdimensionnement des quipements,
Chapitie I Peituibation ue la tension
ENP 2u12 21
Chapitre II
1. Introduction
Dans ce chapitre, nous allons prsenter les diffrentes mthodes danalyse des
principales perturbations lectriques. Pour les perturbations importantes au niveau de
lamplitude (creux de tension, coupures et surtensions), les mthodes destimation
damplitude, de dtection, de classification et de caractrisation sont envisages. Pour ce qui
concerne les perturbations au niveau de la forme donde (harmoniques), la transforme de
Fourier (FFT, DFT) est envisage. Enfin, les techniques destimation et de reprsentation
statistique pour les variations de tension et de frquence sont galement prsentes.
2. Analyse des creux de tension, coupures brves et surtensions
Les creux de tension, les coupures brves et les surtensions se caractrisent par des
variations importantes de tension. Les coupures et les surtensions sont analyses par les
mmes mthodes danalyse que les creux de tension.
2.1 Estimation de lamplitude
Elle se fait par lintermdiaire de la valeur efficace du signal estim dans la plupart
des cas. Il y a notamment dautres mthodes utilises essentiellement bases sur lestimation
de lenveloppe du signal ou de son fondamental.
2.1.1 Valeur efficace
La dtermination de lamplitude en temps rel (de 4 6 priodes fondamentales)
dun signal, se fait par lestimation de sa valeur efficace sur une fentre glissante (relation
II.1), de langueur multiple dune demi-priode. A noter que, si la longueur de la fentre nest
pas un multiple entier de la demi-priode, cela introduira des oscillations dans lamplitude
estime. Pour une fentre de longueur kI
0
avec k H
X
c]]
(t) =
_
1
kI
0
_ x()
2
J
t+
k1
0
2
t-
k1
0
2
(II. 1)
Pour un signal sinusodal, son amplitude est estime par la formule suivante :
A
x
(t) = 2 X
c]]
(t) (II. 2)
Cette mthode est performante et fiable dans le cas sinusodal. Cependant elle ne
lest pas en cas de prsence dharmoniques qui introduisent des erreurs dans lamplitude
estime. En effet, dans ce cas elle ne reflte pas les variations brusques de lamplitude.
Chapitie II Nthoue uanalyse ues peituibations lectiiques
ENP 2u12 2S
2
.
Le signal analytique dun signal de forme sinusodale x(t) = A(t) cos(t + ) est
donc une exponentielle complexe de mme frquence, mme amplitude et mme phase :
x
u
(t) = A(t)c
](ot+ q)
La valeur absolue du signal analytique constitue finalement une estime de lamplitude
instantane du signal x(t) A(t) = |x
u
(t)|
Le signal analytique permet de dtecter instantanment les perturbations au niveau
de lamplitude et destimer avec prcision leur dure. Cette technique est intressante dans le
cas de perturbations de trs courte dure (moins de deux priodes fondamentales).
2.1.3 Composante fondamentale
La composante fondamentale dun signal est dfinie comme suit :
X
]ond
(t) =
2
I
0
__ x() c
-]
2n
1
0
:
J
t+
k1
0
2
t-
k1
0
2
_ (II. 4)
Cette technique prsente lavantage destimer lamplitude du signal correctement et sans
erreurs, mme en prsence dharmoniques. Elle permet donc disoler le signal du reste des
bruits et dautres perturbations harmoniques.
2.2 Segmentation
Afin danalyser une perturbation, elle doit tre dtecte et isole. Le processus de
dcomposition des signaux mesurs en segments constants ou faiblement variables est appel
segmentation . Les mthodes de segmentation peuvent agir dans le domaine temporel,
dans le domaine frquentiel, ou bien tre bases sur un modle du signal [1].
Les mthodes de segmentation dans le domaine temporel utilisent les outils de la
dcision statistique tels que les tests dhypothses squentiels pour dtecter les changements
dans les paramtres du signal.
Chapitie II Nthoue uanalyse ues peituibations lectiiques
ENP 2u12 24
Figure II. 1 Algorithme de classification creux de tension bas sur la mthode des six tensions [1].
Mthode des composantes symtriques :
La mthode des composantes symtriques dtermine le type de creux de tension par
une tude dans le plan complexe des tensions directe et inverse. Les impdances directes et
inverses sont supposes gales. La tension homopolaire est suppose nulle.
Cette mthode analyse les creux de tension tels quils sont subits par les charges, elle
est dveloppe et valable seulement pour les creux de tension de type A, C et D. Lalgorithme
de classification est prsent la figure II.2.
Perturbation mesure
Amplitude des six
Tensions
Equilibre des six tensions Creux de tension non quilibr
Creux de tension biphas
Creux de tension monophas
Creux de tension de type E
Pas de creux de tension
Creux de tension triphas
vmax vmin < Seuil Equ
vmin = vab vbc ou vca
vmin = va, vb ou vc
Nombie vmin > 1
vmin > Seuil Tii
vmin < Seuil Tii
vmax vmin >Seuil Equ v
max
v
min
< Seuil v
max
v
min
> Seuil
v
min
> Seuil Tiiphas
v
min
< Seuil Tiiphas
Chapitie II Nthoue uanalyse ues peituibations lectiiques
ENP 2u12 26
Figure II. 2 Algorithme de classification des creux de tension bas sur les composantes symtriques.
2.4 Caractrisation
La dernire tape dans lanalyse dun creux de tension consiste dterminer la
svrit du creux et son impact sur les quipements lectriques. Deux paramtres sont
observer, lamplitude et la dure. On peut donc reprsenter un creux de tension sur un repre
deux dimensions, amplitude (en pourcentage) en ordonne et dure en abscisse. Ces
courbes obtenues dites courbes de sensibilit sont semblable aux limites fixes par les
courbes de la compatibilit lectromagntique, en effet diffrentes bandes existent pour situer
les zones risque pour diffrents appareils. Des indices sont utiliss pour traduire les rsultats
sous forme algbrique. Parmi lesquels on cite :
La perte de tension : est dfinie comme lintgrale de la tension dans lintervalle du
creux de tension :
I
)
3.14
I
3
=1
(II. 9)
I
ct I
k=1
(II. 1u)
O o
0
reprsente la valeur moyenne du signal (composante continue DC offset), et
Chapitie II Nthoue uanalyse ues peituibations lectiiques
ENP 2u12 28
`
1
1
1
1
o
k
=
2
I
0
_ x(t) cos(k
0
t) Jt
1
0
2
-
1
0
2
b
k
=
2
I
0
_ x(t) sin(k
0
t) Jt
1
0
2
-
1
0
2
0
0
=
2n
I
0
(II. 11)
Ce signal peut, de la mme manire, tre reprsent par une srie de signaux exponentiels
complexes :
< x >
k
=
1
T
0
_ x(t)e
j kmt
ut
T
0
0
(II. 12)
3.1.2 Transformation de Fourier
Les signaux non priodiques peuvent eux aussi sous certaines conditions tre
dcomposs en une somme dexponentielles complexes :
x(t) =
1
2n
_ X()c
]ot
J o X() = _ x(t)c
-]ot
Jt
+
-
+
-
Avec x(t) tant le signal temporelle et X() sa transformer de Fourier dans le domaine
frquentiel.
3.2 Caractrisation
Il existe plusieurs indices de quantification dharmoniques, le plus utilis est dit
taux global de distorsion (Total Harmonic Distortion) |2] :
Taux global de distorsion (THD)
IE = _(_
u
k
u
1
]
2
k=2
(II. 1S)
Taux de distorsion (Facteur de distorsion)
= _
u
k
2
k=2
u
k
2
k=1
(II. 14)
Taux individuel dharmoniques
E
k
=
u
k
_ u
k
2
k=1
selon la CEI (II. 1S)
Chapitie II Nthoue uanalyse ues peituibations lectiiques
ENP 2u12 29
E
k
=
u
k
u
1
selon la CIuRE (II. 16)
Facteur de puissance F
p
, de dphasage cos et de dformation :
`
1
1
1
1
1
1F
p
=
P
S
cos =
P
1
S
1
=
F
p
cos
(II. 17)
Remarque :
Dans le cas dun signal purement sinusodal, le facteur de dformation est gal 1,
dans ce on confond le facteur de puissance et le cos.
4. Dsquilibre
Le dsquilibre dun systme triphas est quantifi par le facteur de dsquilibre,
dont la dfinition pour [10] :
Des mesures entre phase et neutre ou phase et terre
I
% = 1uu_
I
u
+ o`I
b
+ oI
c
I
u
+ oI
b
+ o`I
c
_ (II. 18)
0 I
u
, I
b
, I
u
, I
c
, sont les fondamentaux des trois phases.
Des mesures entre phases
I
% = 1uu_
I
ub
+ oI
bc
I
ub
+ o`I
bc
_ (II. 19)
0 I
ub
, I
bc
, sont les fondamentaux des deux tensions Vab et Vbc.
5. Conclusion
Les diffrentes mthodes danalyse des principales perturbations lectriques sont
prsentes dans ce chapitre. Ces mthodes diffrent selon la dgradation des paramtres qui
caractrisent la tension telle que lamplitude, la frquence, et la forme donde. Dans le
chapitre suivant, nous allons prsenter une nouvelle mthode trs intressante et plus
performante pour lanalyse de ces perturbations qui est base sur la transformation du
Vecteur dEspace.
Chapitre III
ENP 2u12 Su
1. Introduction
Une nouvelle mthode pour lanalyse des perturbations de la tension est prsente
dans ce chapitre. Elle est base sur la transformation du Vecteur dEspace (VE), qui dcrit un
systme triphas par une grandeur complexe, le Vecteur dEspace, et une grandeur relle la
composante homopolaire. Dans le cas dun systme quilibr la composante homopolaire est
nulle, et la trajectoire de lextrmit du VE prend une forme circulaire dans le plan
complexe. Dans le cas o le systme est perturb par une quelconque perturbation, on
observe une modification de la forme du VE, et/ou la composante homopolaire. Ces
modifications sont utilises afin de dtecter, identifier et classer les perturbations mesures.
2. Transformation du Vecteur dEspace
La transformation du Vecteur dEspace est drive de celle des composantes
symtriques. Ce paragraphe prsente un bref historique sur les transformations des
composantes symtriques, et leur lien avec le Vecteur dEspace.
Les composantes symtriques ont t introduites en 1918 par Fortescue. La
transformation de Fortescue, applicable aux grandeurs sinusodales, permet lanalyse des
systmes triphass asymtriques. Plus tard, Lyon reprend la transformation de Fortescue et
lapplique des grandeurs dpendant du temps quelles que soient les conditions
(sinusodales ou non). Cette transforme gnrale et instantane est donne par :
_
x
d
(t)
x
(t)
x
o
(t)
_ =
1
S
_
1 o o
2
1 o
2
o
1 1 1
__
x
u
(t)
x
b
(t)
x
c
(t)
_ (III. 1)
O o = c
]
2n
3
, et ] = -1
Si x
u
,x
b
et x
c
sont valeurs relles, les composantes x
d
et x
I
sont valeurs
complexes conjugues lune de lautre. Ces grandeurs sont redondantes, elles portent alors la
mme information. Dans ce cas seules les deux composantes x
d
et x
0
permettent de dcrire
compltement le systme original. Afin dviter la redondance des composantes directe et
inverse, la transformation de Lyon a tait modifie par Clarke [1] comme suit :
_
x
u
(t)
x
[
(t)
x
o
(t)
_ =
2
S
`
1 -
1
2
-
1
2
u
S
2
-
S
2
1
2
1
2
1
2
/
_
x
u
(t)
x
b
(t)
x
c
(t)
_ (III. 2)
Chapitie III Nthoue uu vecteui uEspace
ENP 2u12 S1
La matrice de transformation ainsi obtenue est cette fois-ci valeurs relles. Les
composantes x
u
et x
u
forment le Vecteur dEspace qui scrit de la manire suivante :
x(t) = x
u
(t) + ]x
[
(t) (III. S)
Et la troisime composante reprsente la composante homopolaire, en combinant (III.2) et
(III.3) la transformation du Vecteur dEspace peut galement tre reprsente sous la forme
matricielle suivante :
_
x(t)
x
o
(t)
_ =
2
S
_
1 o o
2
1
2
1
2
1
2
_ _
x
u
(t)
x
b
(t)
x
c
(t)
_ (III. 4)
En comparant (III.4) avec (III.1), on remarque que la transformation du Vecteur
dEspace calcule les mmes composantes que la transforme de Lyon, mis part la
composante inverse x
1
1
2
o
2
1
2
o
1
2
/
_
x(t)
x
o
(t)
_ (III. S)
Daprs (III.5) on remarque que pour pouvoir reconstruire les grandeurs triphases dorigine,
le Vecteur dEspace et la composante homopolaire sont tous deux ncessaires. Le Vecteur
dEspace peut dcrire seul compltement le systme triphas dorigine seulement si la
composante homopolaire est nulle.
Cette dernire ne peut donc pas tre nglige pour lanalyse des problmes de la Qualit de
lnergie lectrique puisque les perturbations lectriques se traduisent souvent par des
modifications de la composante homopolaire.
Chapitie III Nthoue uu vecteui uEspace
ENP 2u12 S2
= x
p
c
]o
c
t
+x
n
c
-]o
c
t
(III. 7)
O x
p
= |x
p
|c
]q
p
et x
n
= |x
n
|c
]q
n
sont respectivement les phaseurs positif et ngatif.
Exemple : dun cas o les phases et les amplitudes des deux phaseurs sont diffrentes :
Figure III. 1 Reprsentation du Vecteur dEspace.
On montre partir de lquation (III.7), que le grand rayon et le petit rayon de lellipse
dpendent directement des amplitudes des deux phaseurs, alors que langle dinclinaison de
lellipse est fonction de leurs phases [11]:
`
1
1
1
1
r
mu]
= |x
p
| + |x
n
|
r
mn
= |x
p
| - |x
n
|
mn
=
1
2
(
n
+
p
)
(III. 8)
r
maj
r
m|n
m
x
p
I
m
I
m
R
e
I
m
R
e
q
|nc
R
e
x
n
Chapitie III Nthoue uu vecteui uEspace
ENP 2u12 SS
= ( x
pk
(t)c
]ko
c
t
+x
nk
(t)c
-]ko
c
t
) (III. 9)
k=
k=0
O k est le rang des harmoniques prsents dans les grandeurs triphases.
La figure III.2 montre les dformations caractrisant la forme du Vecteur dEspace dans le
plan complexe.
Figure III. 2 Reprsentation du Vecteur dEspace par la somme de deux vecteurs tournants.
=
I
m
I
m
I
m
Re Re
o
-S
o
Chapitie III Nthoue uu vecteui uEspace
ENP 2u12 S4
ENP 2u12 SS
Figure III. 3 Creux de tension de type B.
Les modifications dans le Vecteur dEspace et la composante homopolaire peuvent
tre exprimes thoriquement partir des quations des creux de tension. Par exemple, un
creux de tension de type B sur la phase a est dcrit par les quations suivantes :
`
1
1
1
1
:
u
(t) = (1 -J)I cos(
0
t + )
:
b
(t) = I cos(
0
t + -
2n
S
)
:
c
(t) = I cos(
0
t + +
2n
S
)
J : tant la profondeur du creux de tension.
En appliquant la formule dEuler aux trois tensions de phase et en utilisant la transformation
du Vecteur dEspace (III.4), on peut exprimer le Vecteur dEspace et la composante
homopolaire pour ce type de creux de tension de la manire suivante :
_
x(t)
= [1 -
J
S
c
]
c
]
o
t
+
J
S
c
] (
Sn
S
- n
2n
S
- )
c
-]
o
t
x
o
(t) = -
J
S
Icos [
u
t + +(n +1)
2n
S
O n = 1,2,3 correspond la phase a, b et c respectivement.
Figure III. 4 Vecteur dEspace pour un creux de tension monophas de type B.
d
a
c
b
c
c
b
R
Im
r
m|n
r
maj
q
|nc
Im
R
a
c
b
r
m|n
r
maj
q
|nc
Im
R
a
b
a q
|nc
r
maj
r
m|n
Phase : a
Phase: b
Phase: c
Chapitie III Nthoue uu vecteui uEspace
ENP 2u12 S6
ENP 2u12 S7
nc
=
n
2
,
n
6
ct
5n
6
Pour un creux de tension sur les phases a, b et c respectivement.
La composante homopolaire pour ce type de creux reste gale 0. En effet, ce type
de creux de tension provient de mesures effectues un niveau de tension diffrent du niveau
de tension o le dfaut a t produit, et la composante homopolaire nest pas modifie.
ENP 2u12 S8
Les inclinaisons possibles de lellipse restent les mmes que dans le cas des creux de
tensions de type B et D. La composante homopolaire reste galement nulle dans ce cas.
Figure III. 6 Creux de tension de type F.
Exemple : chute de tension de profondeur d = 40% sur la troisime phase et un dphasage de
+9 et -9 respectivement sur la phase 1 et 2, plus une chute de tension de 17% sur leur
amplitudes.
Les rsultats de simulation sont lannexe (figure 3).
La comparaison entre les rsultats prdtermins par les formules thoriques et ceux donns
par lapplication est donne dans le tableau III.3.
Rsultats prdtermins par les formules thoriques Rsultats de simulation
r
mu]
(I) r
mn
(I)
nc
() |x
o
| r
mu]
(I) r
mn
(I)
nc
() |x
o
|
269.643 186.6762 -30 0 310.917 217.718 -30.040 0.530
Tableau III. 3 Rsultats de lapplication de la mthode du Vecteur dEspace.
Vecteur dEspace pour les creux de tension monophass :
Le tableau III.4 rsume les caractristiques du Vecteur dEspace et de la composante
homopolaire dans le cas des creux de monophass, o n = 1, 2,3 correspond la phase a, b et
c respectivement.
b
d
a
c
Chapitie III Nthoue uu vecteui uEspace
ENP 2u12 S9
nc
=
n
3
pour une chute de tension sur les phases a et b et
nc
=
2n
3
en cas de creux de
tension sur les phases a et c. Enfin, ce type de creux de tension se caractrise par une
composante homopolaire nulle.
Figure III.1 Creux de tension de type C.
a
b
c
d
Chapitie III Nthoue uu vecteui uEspace
ENP 2u12 4u
ENP 2u12 41
Exemple : creux de tension de profondeur d = 50% sur les phases 1 et 2. Les rsultats de
simulation sont lannexe (figure 4).
La comparaison entre les rsultats prdtermins par les formules thoriques et ceux donns
par lapplication est donne dans le tableau III.6.
Rsultats prdtermins par les formules thoriques Rsultats de simulation
r
mu]
(I) r
mn
(I)
nc
() |x
o
| r
mu]
(I) r
mn
(I)
nc
() |x
o
|
259.250 155.550 -60 51.830 259.280 155.561 -59.963 51.860
Tableau III. 6 Rsultats de lapplication de la mthode du Vecteur dEspace.
Creux de tension de type G :
Le creux de tension type G est caractris par des modifications des tensions des
trois phases le long daxes diffrents. Le Vecteur dEspace correspondant parcourt galement
une ellipse dans le plan complexe avec des rayons dpendants de la profondeur de creux et
des inclinaisons identiques celles en cas de creux de tension de type C et E. Et la
composante homopolaire reste nulle.
Figure III.3 Creux de tension de type G.
Exemple : creux de tension de profondeur d = 12 .5% sur les phases 1 et 2, avec un
dphasage supplmentaire de 10.
Les rsultats de simulation sont lannexe (figure 6).
La comparaison entre les rsultats prdtermins par les formules thoriques et ceux donns
par lapplication est donne dans le tableau III.7.
b
a
c
d
Chapitie III Nthoue uu vecteui uEspace
ENP 2u12 42
ENP 2u12 4S
ENP 2u12 44
ENP 2u12 4S
ENP 2u12 46
enregistres sont dcoupes en segments durant lesquels le systme triphas est sain ou subit
un creux de tension. Les segments pour lesquels un creux de tension est dtect sont ensuite
analyss et le type exact des creux est dtermin. La caractrisation est la dernire tape,
durant laquelle la gravit des creux de tension est estime.
4.2.1 Extraction des traits caractristiques des tensions mesures
Le Vecteur dEspace et la composante homopolaire sont dtermins par la
transformation du Vecteur dEspace partir des tensions mesures (quation III.4). Les
paramtres du Vecteur dEspace utiliss pour lanalyse des creux de tension sont les rayons et
linclinaison de la forme (cercle ou ellipse) quil parcourt dans le plan complexe. Ces
paramtres ne peuvent pas tre correctement estims gomtriquement en utilisant
directement cette forme. En effet, les harmoniques et le bruit prsents le plus souvent au sein
des tensions mesures dforment le cercle ou lellipse que devrait parcourir le Vecteur
dEspace. Ces dformations peuvent introduire des erreurs importantes dans lestimation des
paramtres du Vecteur dEspace et par la suite, dans les tapes suivantes danalyse des creux
de tension.
Afin dviter leffet nfaste des harmoniques et du bruit, on cherche estimer les
phaseurs positif et ngatif du Vecteur dEspace la frquence fondamentale
(x
p
et x
n
). A partir de ces deux grandeurs, les caractristiques de la forme du Vecteur
dEspace dans le plan complexe sont correctement estimes. Afin dobtenir x
p
et x
n
, la
transforme de Fourier est applique au Vecteur dEspace. Pour la frquence fondamentale,
on obtient :
X() = _ (x
p
c
]o
c
t
+ x
n
c
-]o
c
t
-
)c
]ot
Jt = _ x
p
c
-2n(]+]
c
)t
Jt
-
+ _ x
n
c
-2n(]-]
c
)t
Jt (III. 1u)
-
X() = x
p
o( +
o
) +x
p
o( -
o
) (III. 11)
O
o
et
o
sont respectivement la frquence et la pulsation fondamentale et o est
limpulsion de Dirac.
Lquation III.11 montre que la transforme de Fourier permet disoler facilement x
p
et x
n
grce aux impulsions de Dirac
o
. Les modules de x
p
et x
n
sont utiliss pour calculer les
rayons de lellipse et leurs phases pour estimer son inclinaison (relation III.8).
Chapitie III Nthoue uu vecteui uEspace
ENP 2u12 47
Figure III.6 Spectre du Vecteur dEspace pour la frquence fondamentale.
Lapplication de la transforme de Fourier lanalyse du fondamental de la tension
homopolaire et du Vecteur dEspace permettent une estimation prcise de leurs paramtres et
garantissent le fonctionnement et la fiabilit de la mthode, mme en cas de distorsions
harmoniques importantes ou de mesures bruites.
Les diffrentes transformes de Fourier sont estimes sur des signaux chantillonnes et sont
donc calcules numriquement laide de transformes de Fourier discrtes (TFD). On peut
pour cela utiliser un algorithme rapide tel que celui de la FFT (Fast Fourier Transform). Une
illustration de lanalyse dun cas rel par EPQA V1.0.0 (tensions du rseau SONELGAZ
mesures dans le laboratoire LRE de lENP) est prsente lannexe (figure 8).
4.2.2 Segmentation
Le Vecteur dEspace et la composante homopolaire permettent de simplifier ltape
de dtection des creux de tension par rapport aux mthodes classiques prsentes au chapitre
II en utilisant seulement deux variables pour la segmentation et en vitant les algorithmes de
comparaison des segments. Les deux variables utilises pour la segmentation sont le petit
rayon de la forme parcourue par le Vecteur dEspace dans le plan complexe, et lamplitude de
la composante homopolaire.
|x()|
x
p
o
-
o
x
n
Chapitie III Nthoue uu vecteui uEspace
ENP 2u12 48
Figure III.7 Algorithme de segmentation et destimation de la zone de creux de tension [1].
a) Estimation de Rmin :
Le Vecteur dEspace est estim chaque chantillon du signal en utilisant la
transforme instantane (III.4). Ensuite, la transforme de Fourier discrte est applique ce
Vecteur dEspace chaque chantillon et sur une fentre glissante dune priode
fondamentale. Les valeurs des phaseurs correspondants la frquence fondamentale positive
(xp) et ngative (xn) sont alors releves, et le petit rayon est dtermin partir de lquation
(III.8).
Le calcul de la TFD sur une fentre glissante dune priode implique que les
variations brusques de lamplitude des tensions triphases seront reprsentes par des
variations dau moins une priode dans r
mn
.
b) Segmentation de r
mn
Les variations dans la valeur de r
mn
sont dtectes par une double application de
lalgorithme Cusum. Cet algorithme a t choisi grce sa simplicit et ses bonnes
performances. Lalgorithme Cusum part de lhypothse que le signal dintrt est r
mn
de
type alatoire, chantillons indpendants et identiquement distribus. Sa loi est suppose
tre de type gaussien, et sa moyenne est constante par morceaux. Son principe gnral est
deffectuer un test dhypothse squentiel ralis chaque chantillon, afin de dcider si la
valeur moyenne des tous derniers chantillons a dvi de sa valeur courante. Ce test est bas
Donne relles
Estimation de Rmin
Segmentation de Rmin
Estimation de Rmin Estimation de Xo
Sgmentation de Xo
Pas de creux dtect
Zone de creux de tension
Nombre de changement >=1 Nombre de changement = 0
Nombre de changement >=1
Nombre de changement = 0
Chapitie III Nthoue uu vecteui uEspace
ENP 2u12 49
ENP 2u12 Su
4.2.3 Classification
Lintroduction de la mthode du Vecteur dEspace permet dutiliser de nouvelles
techniques de classification diffrentes de celles dj vu au chapitre II. Nous allons prsenter
lune de ces mthodes. Elle permet une classification exhaustive des creux de tension en
diffrentiant galement les creux de tension avec signature semblables.
Lalgorithme de classification des creux de tension est prsent en figure III.15. Il
est constitu des tapes suivantes :
La transformation du Vecteur dEspace est applique aux grandeurs triphases pour le
segment analys (quation III.4).
Les phaseurs du Vecteur dEspace sont estims la frquence fondamentale positive
et ngative en utilisant la transforme de Fourier discrte. Ils sont ensuite utiliss pour
calculer les caractristiques du Vecteur dEspace dans le plan complexe (facteur de forme SI,
angle dinclinaison
nc
, et rayons r
mn
et r
mn
de lellipse).
nc
r
mu]
, x
o
x
o
Chapitie III Nthoue uu vecteui uEspace
ENP 2u12 S1
les creux non quilibrs des creux triphass ou avec des surtensions doit donc tre
dtermine. La dfinition des creux de tension vue au chapitre I permet daffirmer que la
profondeur minimale dun creux est de d = 0,1p.u, ce qui correspond un facteur de forme
maximal allant de 0,86 0,93 en fonction du type de creux. De ce fait, un creux de tension
avec un facteur de forme suprieur 0,93 est considr comme un creux de tension triphas
ou avec surtension, et un creux de tension avec un facteur de forme infrieur cette valeur est
considr comme un creux de tension monophas ou biphas.
Linclinaison de lellipse diffrencie les creux monophass des creux biphass et
dtermine la ou les phase(s) en chute. En effet, langle entre le grand rayon de lellipse et
laxe rel dfini la figure III.16 permet de dterminer les phases en chute. Ce rsultat est
illustr la figure III.15, o les creux monophass sont dsigns par Mono et les creux
biphass par Bi, suivi des phases affectes par le creux.
Figure III.16 Type de creux de tension en fonction de langle dinclinaison de lellipse.
Langle dinclinaison de lellipse prend rarement des valeurs exactement multiples de 30.
Pour cette raison, lindex k est introduit afin de dterminer le type de creux et de faciliter
lintgration de lalgorithme dans un logiciel ddi. La relation entre lindex k, langle
dinclinaison de lellipse et les types de creux de tension sont donns au tableau III.12.
nc
u _ 1S Su _ 1S 6u _ 1S 9u _ 1S 12u _1S 1Su _1S
Mono a
Mono a
Bi ab
Mono b
Bi bc
Mono c
Bi ac Bi ab
Mono b
Bi bc
Mono c
Bi ac
Chapitie III Nthoue uu vecteui uEspace
ENP 2u12 S2
k
u 1 2 S 4 S
Type de creux
Bi ob Hono b Bi ob Hono o Bi oc Hono c
Tableau III. 12 Type de creux de tension en fonction de linclinaison de lellipse [1].
Dterminer le type du creux de tension en utilisant la tension homopolaire et le grand
rayon de lellipse. Comme dans le cas de facteur de forme, des seuils doivent tre introduits
pour la composante homopolaire et le grand rayon de lellipse. Ces seuils peuvent tre
calculs partir de la profondeur du creux estime partir de la valeur de r
mn
, ou
dtermins par des observations du rseau analys.
Diffrencier les creux de tension triphass et les creux de tension avec surtension, Si
la perturbation mesure est classifie comme un creux de tension triphas ou comme un creux
de tension avec surtension, la composante homopolaire est utilise pour la dtermination
exacte du type de creux. Si lamplitude de la tension homopolaire est en dessous dun certain
seuil positif, le creux de tension est classifi comme triphas. Dans le cas contraire, le creux
de tension mesur se caractrise par des surtensions. Dans ce cas, la phase de la composante
homopolaire indique le type de creux : monophas (type H) ou biphas (type I), ainsi que les
phases affectes par le dfaut comme indiqu la figure III.16.
Figure III.9 Type de creux de tension en fonction de la phase de la composante homopolaire [1].
Mono b
Bi bc
Mono c
Bi bc
Mono a
Bi ab
Chapitie III Nthoue uu vecteui uEspace
ENP 2u12 SS
4.2.4 Caractrisation
Deux approches pour estimer la gravit des creux de tension par la transformation du
Vecteur dEspace sont possibles. La premire approche consiste valuer la gravit des creux
de tension par leur profondeur. Pour des creux de tension de type de A G, la gravit des
creux est estime par le petit rayon de lellipse r
mn
, qui dpend directement de la
profondeur des creux de tension d (voir Tableau III.1-3). Pour des creux de tension de type H
et I, le petit rayon de lellipse reste gal la tension nominale. Dans ce cas, lamplitude de la
composante homopolaire est utilise pour estimer la gravit des creux de tension, car elle est
proportionnelle leur profondeur (voir Tableau III.11).
La deuxime approche consiste estimer la gravit des creux de tension par la
surface de la forme parcourue par le Vecteur dEspace dans le plan complexe. En effet, elle
est directement proportionnelle limpact des creux de tension sur les quipements. Par
exemple, pour une dure et une profondeur de creux de tension donnes, les creux de tension
triphass sont les plus svres, suivis par les creux de tension biphass puis monophass. Les
creux de tension avec surtensions se caractrisent par la mme surface pour la forme
parcourue par le Vecteur dEspace quen conditions normales. En effet, ces types de creux de
tension ne se propagent pas et ne sont pas vus par les charges gnralement connectes en
triangle au niveau MT [1], et ne sont donc pas considrs comme svres. Cette approche est
donc trs approprie pour lanalyse des creux de tension triphass.
5. Conclusion
Lanalyse des problmes de la Qualit de lnergie lectrique par la transformation
du Vecteur dEspace est prsente dans ce chapitre. Des approches spcifiques pour les
diffrents types de perturbations (creux de tension, surtensions, dsquilibre et harmoniques)
ont t dveloppes. La transformation du Vecteur dEspace extrait les traits caractristiques
des creux de tension et facilite toutes les tapes ncessaires leur analyse. Ainsi, la
segmentation des creux de tension est effectue de manire plus fiable et ncessite lanalyse
de moins de grandeurs. Le type de creux de tension est dtermin grce aux caractristiques
du Vecteur dEspace et de la composante homopolaire. Enfin, pour estimer la gravit des
creux de tension, deux techniques spcialement destines aux creux de tension dans les
systmes triphass sont proposes. La transformation du Vecteur dEspace est galement
performante pour lanalyse dautres types de perturbations. Lanalyse des surtensions est
effectue dune manire analogue celle des creux de tension. Le degr de dsquilibre est
estim par le rapport des composantes spectrales du Vecteur dEspace aux frquences
Chapitie III Nthoue uu vecteui uEspace
ENP 2u12 S4
Chapitre IV
Application ddie la
surveillance en temps rel
Chapitie Iv Application uuie la suiveillance en temps iel
ENP 2u12 SS
1. Introduction
Dans ce chapitre, nous prsentons lapplication ralise travers une implantation
dun organigramme fonctionnel. Chaque menu, sous menu et fentre est dcrit, visualis et
discut. Lapplication est dveloppe sous MATLAB R2009b, qui est un logiciel qui permet
la programmation oriente objet.
Il comporte plusieurs librairies (toolboxes) ddies aux diffrents domaines de
lingnierie moderne. Comme il permet aussi de dvelopper des applications sous un
environnement Windows en deux modes programming et GUIDE (guide user interface
Development Environment).
Environnement de dveloppement des interfaces graphiques (GUIDE)
Pour accder cet outil, lutilisateur de MATLAB peut insrer la commande
guide qui affiche lcran suivant :
La fentre menu et outils comporte :
Les boutons New , Open et Save qui permettent la cration, louverture et
lenregistrement des projets.
Les boutons Align Object , Property Inspector et Object Browser qui
permettent respectivement lalignement des objets, linspection des proprits de
lobjet slectionn et lexploration des diffrents objets contenus dans le projet.
Chapitie Iv Application uuie la suiveillance en temps iel
ENP 2u12 S6
ENP 2u12 S7
ENP 2u12 S8
ENP 2u12 S9
(t) = (1 -J
)I sin(2n
0
t +
+ (1 -i)
2n
S
)
O i est lindice de la phase i, i = 1: S, (i est choisi selon le type du systme
dalimentation),
J
(t), dintensit d %.
Lutilisateur dispose de deux choix du type de bruit : distribution normale ou uniforme,
lintensit d en pourcent est relative lamplitude de chaque phase.
2.2 Charge :
2.2.1 Algorithme du programme
Dbut
Chargement de
donnes
Reconnaissance de type du systeme
Tracer (appel des fonctions
de calcul)
Enrigstrement
Fin
Choix du
couplage
Paramtres de
la charge
Chapitie Iv Application uuie la suiveillance en temps iel
ENP 2u12 6u
+b(t
, u
) i 1: N
Le pas b est dynamique dans le programme, il est en fonction de la rapidit du
systme.
2.3 Redresseur :
ENP 2u12 61
2.4 Acquisition :
Dbut
Inspction du materiels
Collection des proprits des cartes
installes
Choix de la carte
Fs,T,Voltage range
Acquisition priodique
chixs des canaux
Dbut
d'acquisition
Affichage de
donnes
Priodique
Fin
OUI
NON
Chapitie Iv Application uuie la suiveillance en temps iel
ENP 2u12 62
2.5 Analyse :
_ =
2
S
_
1 o o
2
1
2
1
2
1
2
_ _
x
u
(t)
x
b
(t)
x
c
(t)
_
Les spectres des tensions sont calculs par la transforme de Fourier rapide, et les
phaseurs positifs ngatifs sont reprs par leurs frquences + f0 et f0 respectivement.
Lalgorithme de dtection des creux de tension est bas sur la comparaison des
valeurs crtes des trois tensions sur une fentre glissante de langueur multiple dune demi-
priode fondamentale.
La mthode utilise pour lestimation de la valeur efficace est prsent dans le
chapitre II 2.1.1
Dbut
Chargement de
donnes
Calcul du spectre
Controle d'affichage choix
des cannaux
Taille >3
Activation de la
slection du
systeme
Appels des fonctions : analyse du systeme,
(Application de la mthode du vecteur espace)
Acquisition
Chapitie Iv Application uuie la suiveillance en temps iel
ENP 2u12 6S
ENP 2u12 64
O Di : reprsente la profondeur du creux de la phase i en %.
Phi : reprsente le dphasage supplmentaire de la phase i en .
Fr : reprsente la frquence fondamentale du systme,
Xeff : reprsente lamplitude maximale du fondamental des trois phases.
Xref : reprsente la valeur max des amplitudes estimes du systme,
Le panneau de contrle daffichage offre la possibilit de choisir les signaux visualiser :
Le canal I permet la visualisation :
Des signaux dans le domaine temporel,
Le systme triphas slectionn,
Le Vecteur dEspace (partie relle et imaginaire),
La composante homopolaire (en cas dquilibre, elle est nulle),
Le spectre du Vecteur dEspace,
Le canal II permet la visualisation de :
Des signaux dans le domaine frquentiel,
Du Vecteur dEspace dans le plan complexe,
Diagramme de Fresnel du systme choisi.
Le THD (Taux global de distorsion) du canal slectionn est galement donn.
Dans le dernier panneau Rsultats , on trouve :
Rmin, Rmaj : le rayon min et max de lellipse respectivement,
SI : le coefficient de forme de lellipse,
Xp, Xn : les amplitudes estimes des phaseurs positif et ngatif respectivement,
Inc. : reprsente linclinaison de lellipse.
Si le panneau (choix du systme) est dsactiv, tous les champs relatifs aux rsultats
danalyse du systme triphas sont initialiss.
Chapitie Iv Application uuie la suiveillance en temps iel
ENP 2u12 6S
ENP 2u12 66
3.4 Charge:
ENP 2u12 67
ENP 2u12 68
Enfin lutilisateur peut enregistrer ses rsultats de simulation en format (*.mat, *.xls).
3.6 Acquisition:
La fentre prsente la figure IV.12 reprsente la partie pratique de lapplication, elle
assure la communication avec plusieurs types de cartes dacquisition laide dun adaptateur
appropri, ce dernier permet lapplication de communiquer avec le hardware. Toutes les
cartes compatibles avec MATLAB sont prises en charge par cette application si elles ont
moins de 16 canaux sauf les cartes port srie.
ENP 2u12 69
Acquisition priodique : cette option nous permet deffectuer des acquisitions priodes
selon la priode choisie, soit par minute ou bien par heure,
Canaux de slection et facteurs de calibrage : on choisit les canaux et les facteurs de
calibrage qui correspondent lacquisition voulue,
En cliquant sur le Bouton Appliquer , les paramtres utilisateurs sont enregistrs
et la fentre se ferme automatiquement. Pour commencer lacquisition il suffit dappuyer sur
dmarrer pour lancer lacquisition.
Lors dune acquisition, lapplication affiche le premier canal, lutilisateur peut visualiser
galement les autres entres en appuyant sur superposer puis slectionner une autre entre.
Et lenregistrement en format (*.mat ou *.xls) peut se faire en cliquant sur enregistrer .
En cas dchec de lacquisition, le programme informe lutilisateur sur le type
dincident, avec les messages derreur.
ENP 2u12 7u
4. Validation de lapplication
Cette partie consiste comparer les rsultats obtenus par EPQA et ceux de
Simplorer. En effet, afin de valider notre application, une comparaison avec un logiciel de
rfrence est primordiale.
La premire comparaison consiste valider les programmes de simulations en
comparant les rsultats obtenus par EPQA et ceux obtenus par Simplorer, cela nous
permettra dvaluer la prcision dEPQA du point de vue intgration numrique.
La deuxime comparaison consiste valider le programme dacquisition en
comparant les donnes acquises par EPQA et DAQVIEW.
4.1 Redresseur diodes :
Premier cas :
ENP 2u12 71
ENP 2u12 72
ENP 2u12 7S
On reprend le mme montage avec la mme charge, mais les diodes D1, D5, D6 dconnectes
ENP 2u12 74
ENP 2u12 7S
ENP 2u12 76
ENP 2u12 77
ENP 2u12 78
ENP 2u12 79
ENP 2u12 8u
Figure IV. 38 les courants de source dune charge RLC quilibre exemple 1.
ENP 2u12 81
Exemple 2 : charge non linaire (redresseur alimentant une charge R-L) alimente
par un systme de tension quilibr. Nous procderons par les mmes tapes que
prcdemment pour la gnration des tensions quilibres, puis on charge le fichier des
tensions par la fentre Redresseur . Nous choisissons la configuration suivante :
Toutes les diodes sont connectes,
Mode en charge, R=20, L=100mH.
Les rsultats de simulation sont prsents dans la Figure.IV.39.
ENP 2u12 82
ENP 2u12 8S
ENP 2u12 84
Figure IV. 45 Vecteur dEspace (systme de courant) Exemple 3.
Figure IV. 46 paramtres du dsquilibre (courant) Exemple 3.
Lestimation du dsquilibre en amplitude est perdue car le modle utilise suppose
que le signal est sinusodal voir la relation (II.2). A cet effet, le champ Xref de la figure.IV.46
propose une autre amplitude estime partir de la moyenne des maximums sur une fentre
glissante dune longueur dune priode fondamentale.
Le dsquilibre en phase est estim par la mthode du passage par zro, la phase 1 est
toujours prise comme origine des phases.
Chapitie Iv Application uuie la suiveillance en temps iel
ENP 2u12 8S
ENP 2u12 86
ENP 2u12 87
Conclusion
Gnrale
Conclusion uniale
ENP 2u12 88
Dans ce travail, nous avons prsent les diffrentes mthodes danalyse des
perturbations lectriques. Ceci nous a permis de conclure que la mthode du Vecteur
dEspace est la plus adquate. En vue du dveloppement dune application pour le suivi en
temps rel de la Qualit de l'Energie Electrique, nous avons utilis la mthode du Vecteur
dEspace en se basant sur le Vecteur dEspace et la composante homopolaire.
Lapplication propose permet de simuler le dsquilibre dun systme
dalimentation et ses consquences sur les charges linaire et non linaire (redresseur triphas
double diodes).Cette application permet aussi lanalyse de donnes obtenues soit par la
simulation soit par lacquisition.
Nous avons compar les rsultats obtenus en utilisant lapplication EPQA dune part
avec ceux obtenus par le logiciel de simulation Simplorer, et dautre part avec ceux acquis
partir de DAQVIEW. Cette comparaison nous a permis de conclure que notre application
donne de bons rsultats avec une erreur relative maximale de 0,1%, ce qui nous a donc permis
de valider EPQA.
Le travail effectu sur la surveillance de la Qualit de lEnergie Electrique nest pas
encore achev. A cet effet, nous proposons denvisager des travaux ultrieurs, notamment
dans le dveloppement de lapplication EPQA. En effet, lapplication dveloppe a la
possibilit dtre enrichie par dautres modules et options comme la Tlsurveillance. Pour
cela, nous avons ouvert la communication avec le module Analyse .
Dautres projets pourront tre faits en association avec les industriels pour essayer de
leur apporter des solutions aux problmes lis la QEE auxquels ils sont confronts surtout si
on tient compte des efforts de normalisation qui imposent aux professionnels du domaine
nergtique des rgles de plus en plus strictes.
ENP 2u12
Bibliographie :
ENP 2u12
Annexe
Figure1Creuxdetension(surlaphase2)detypeB.
Figure2Creuxdetension(surlaphase1)detypeD.
ENP 2u12
Figure3CreuxdetensionmonophasdetypeF.
Figure4CreuxdetensiondetypeC.
ENP 2u12
Figure5CreuxdetensiondetypeE.
Figure6CreuxdetensiondetypeG.
ENP 2u12
Figure7Creuxdetensiontriphas(typeA).
Figure8ApplicationduVecteurdEspaceuncasrel.