Compensateurs Actifs en Basse Tension
Compensateurs Actifs en Basse Tension
THSE
prsente pour obtenir le grade de
Mohamad Alaa Eddin Alali Contribution lEtude des Compensateurs Actifs des Rseaux Electriques Basse Tension
(Automatisation des systmes de puissance lectriques)
Membres du jury
Directeur de thse : Directeur de thse : Rapporteur interne : M. Shahrokh Saadate, professeur, Universit Henri Poincar, Nancy M. Francis Braun, professeur, Universit Louis Pasteur, Strasbourg M. Eric Ostertag, professeur, Universit Louis Pasteur, Strasbourg
Rapporteur externe : M. Jean-Marie Kauffmann, professeur, Universit Franche-Comt, Belfort Rapporteur externe : M. Brayima Dakyo, professeur, Universit du Havre, Le Havre Examinateur : Membre invit : Membre invit : M. Yves-Andr Chapuis, M.C, Universit Louis Pasteur, Strasbourg M. Pascal Boos, Ingnieur, Socit Socomec, Benfeld M. Mohamed Machmoum, M.C, IRESTE, Saint-Nazaire
THSE
prsente pour obtenir le grade de
Mohamad Alaa Eddin Alali Contribution lEtude des Compensateurs Actifs des Rseaux Electriques Basse Tension
Soutenue publiquement le 12 Septembre 2002
Membres du jury
Directeur de thse : Directeur de thse : Rapporteur interne : M. Shahrokh Saadate, professeur, Univ. Henri Poincar, (Electrotechnique) M. Francis Braun, professeur, Univ. Louis Pasteur, (Electrotechnique) M. Eric Ostertag, professeur, Univ. Louis Pasteur, (Automatique)
Rapporteur externe : M. Jean-Marie Kauffmann, professeur, Univ. Franche-Comt, (Electrotechnique) Rapporteur externe : M. Brayima Dakyo, professeur, Univ. du Havre, (Electrotechnique) Examinateur : Membre invit : Membre invit : M. Yves-Andr Chapuis, M.C, Univ. Louis Pasteur, (Microlectronique) M. Pascal Boos, Ingnieur, Socit Socomec, (Electrotechnique) M. Mohamed Machmoum, M.C, IRESTE, Saint-Nazaire, (Electrotechnique)
AVANT PROPOS
Le travail prsent dans ce mmoire a t effectu en codirection entre les deux laboratoires : GREEN (Groupe de Recherche en Electrotechnique et en Electronique de Nancy) de lUniversit Henri Poincar, UHP-Nancy1, et LEPSI (Laboratoire dElectronique et de Physique des Systmes Instrumentaux) de lUniversit Louis Pasteur, ULP-Strasbourg1. Je remercie monsieur Jean-Marie KAUFFMANN, professeur lUniversit de Franche-Comt Belfort, qui a accept dtre rapporteur de mon travail et qui ma fait lhonneur de prsider mon jury. Je tiens remercier mes encadreurs sans qui cette thse naurait jamais eu lieu : messieurs Shahrokh SAADATE professeur lUHP, Francis BRAUN professeur lULP et Yves-Andr CHAPUIS matre de confrences lULP. Leur encadrement ma offert la possibilit de raliser une thse multidisciplinaire ce qui ma permis dlargir mon spectre de connaissances scientifiques et douvrir la voie des collaborations futures en forme de co-tutelles. Je tiens remercier galement aussi la direction des deux laboratoires daccueil, messieurs Abderrezak REZZOUG, professeur lUHP, directeur du GREEN et Jean-Louis RIESTER professeur lULP, directeur du LEPSI, de mavoir offert la possibilit de raliser ma thse au sein de leurs laboratoires. Je remercie monsieur Eric OSTERTAG, professeur lUniversit Louis Pasteur Strasbourg, davoir accept dtre rapporteur de ma thse et avec qui jaurai le plaisir de travailler sur des sujets portant sur des thmatiques scientifiques qui nous intressent. Je remercie monsieur Brayima DAKYO, professeur lUniversit du Havre, davoir accept dtre rapporteur de ma thse et qui a enrichi le dbat avec ses remarques trs pertinentes. Je tiens remercier aussi messieurs Pascal BOOS, ingnieur la socit SOCOMEC Benfeld et Mohamad MACHMOUM, matre de confrences lIRESTE Saint-Nazaire, davoir accept leur participation au jury de thse. Je suis trs reconnaissant la socit de SOCOMEC de Benfeld du partenariat industriel quils nous ont confis. Je remercie spcialement toutes les personnes qui ont particip ce projet messieurs : Michel BOMMER, Pascal BOOS, Stphane MERCIER, Emmanuel GRUBER et Thierry DANOBER. Finalement, je tiens remercier galement monsieur Jean LEVEQUE, matre de confrences au GREEN de lUHP, de mavoir soutenu tout au long de la thse.
SOMMAIRE
INTRODUCTION GNRALE ............................................................................................. 1 CHAPITRE I PERTURBATIONS DES RSEAUX LECTRIQUES ET PRINCIPE DE COMPENSATION ACTIF..................................................................................................... 5
I. INTRODUCTION ............................................................................................................................7 II. CARATERISTIQUES DES PERTURBATIONS ELECTRIQUES .........................................................7
II.1 Perturbations harmoniques en courant et en tension ................................................................. 8
II.1.1 Origine des harmoniques............................................................................................................................. 8 II.1.2 Consquences des harmoniques .................................................................................................................. 8 II.1.3 Rglementation ......................................................................................................................................... 10
IV.3 La combinaison parallle-srie actifs (UPQC) ......................................................................... 20 IV.4 Combinaison hybride active et passive ..................................................................................... 20
IV.4.1 Le filtre actif srie avec des filtres passifs parallles............................................................................... 21 IV.4.2 Le filtre actif srie connect en srie avec des filtres passifs parallles .................................................. 21 IV.4.3 Le filtre actif parallle avec un filtre passif parallle .............................................................................. 21
V. CONCLUSION.......................................................................................................................23
II.3.3 Le systme base de PLL......................................................................................................................... 39 II.3.3.1 Structure gnrale............................................................................................................................. 39 II.3.3.2 Principe de fonctionnement .............................................................................................................. 39 II.3.3.3 Calcul de la rfrence poursuivre Vd-ref.......................................................................................... 40 II.3.4 Rgulation de la tension continue ............................................................................................................. 40 II.3.5 La commande de londuleur...................................................................................................................... 41 II.3.5.1 Commande par hystrsis................................................................................................................. 41 II.3.5.2 Commande par modulation de largeur dimpulsion ......................................................................... 42 II.3.6 La Rgulation du courant du filtre actif parallle...................................................................................... 43 II.3.6.1 Structure de la boucle de rgulation ................................................................................................. 43 II.3.6.2 Mthode de rgulation par PI ........................................................................................................... 44 II.3.6.3 Mthode de rgulation par RST........................................................................................................ 45
II.4 Conclusion..................................................................................................................................... 45
III.3.1.2 Placement de ples du rgulateur RST............................................................................................ 52 III.3.1.3 Analyse des Performances .............................................................................................................. 53 III.3.2 Rgulateur RST appliqu la rgulation du courant du F.A.P ................................................................ 55 III.3.2.1 Description polynomiale du rgulateur RST................................................................................... 55 III.3.2.2 Placement de ples du rgulateur RST............................................................................................ 55 III.3.2.3 Analyse des Performances .............................................................................................................. 56
V.2 Cas de la Compensation dans un milieu perturb..................................................................... 71 V.3 Cas de la compensation dun filtre actif parallle structure classique ................................. 73
V.3.1 Compensation des courants harmoniques sans rejet des composantes dues aux dcoupages ................... 73 V.3.2 Compensation des courants harmoniques avec rejet des composantes dues aux dcoupages................... 75
II.3.2.2 Identification des creux de tension biphass..................................................................................... 92 II.3.2.3 Identification des creux de tension triphass .................................................................................... 93
II.6.2 Compensation des harmoniques et des dsquilibres de tension ............................................................ 113 II.6.2.1 Rgulateur PID ............................................................................................................................... 113 II.6.2.2 Rgulateur RSTam ........................................................................................................................... 115 II.6.3 Evaluation des mthodes didentification et de rgulation...................................................................... 116
III. CONCLUSION ET PERSPECTIVES SUR LA DPOLLUTION UNIVERSELLE DES RSEAUX ELECTRIQUES .......................................................................................................................118
III.1 Conclusion sur le filtre actif srie............................................................................................ 118 III.2 Dpollution universelle : Combinaison Parallle-Srie Actifs (UPQC) ............................... 119
III.2.1 Structure gnrale et emplacement de lUPQC ..................................................................................... 119 III.2.2 Rsultats de simulation .......................................................................................................................... 120
CHAPITRE IV FAISABILIT EXPRIMENTALE DUN FILTRE ACTIF PARALLLE UNIVERSEL ................................................................................................................................... 125
I. INTRODUCTION ........................................................................................................................127
I.1 Problmatiques de recherche...................................................................................................... 127 I.2 Contextes et objectifs de recherche ............................................................................................ 127
IV.2 Calcul direct dans le plan z des rgulateurs numriques ...................................................... 149
IV.2.1 Calcul du rgulateur RSTam(z) du courant............................................................................................. 149 IV.2.2 Etudes des performances du rgulateur RSTam(z) en boucle ferme ..................................................... 149 IV.2.2.1 Etude des performances en terme de poursuite ............................................................................. 149 IV.2.2.2 Etude des performances en compensation des perturbations ........................................................ 150
IV.[Link] Effets des perturbations de la tension du rseau ........................................................................................150 IV.[Link] Effets des perturbations de londuleur.......................................................................................................151
ANNEXE II ..................................................................................................................................180
1. Compensation des creux de tension biphass ............................................................................ 180
1.1 Rgulation de tension par PID ................................................................................................................... 180 1.2 Rgulation de tension par RST .................................................................................................................. 181 1.3 Rgulation de tension par RSTam ............................................................................................................... 182
2.
2.1 Rgulation de tension par PID ................................................................................................................... 184 2.2 Rgulation de tension par RST .................................................................................................................. 185 2.3 Rgulation de tension par RSTam ............................................................................................................... 186
Introduction Gnrale ___________________________________________________________________________ Lutilisation croissante dans lindustrie de systmes commands base dlectronique de puissance entrane de plus en plus de problmes de perturbation au niveau des rseaux lectriques. Ainsi, on assiste une augmentation rgulire, de la part des utilisateurs, des taux dharmonique et de dsquilibre des courants, ainsi qu une importante consommation de la puissance ractive. La circulation de ces mmes courants perturbs va galement provoquer des harmoniques et des dsquilibres de tension, lesquels vont se superposer la tension nominale du rseau lectrique. De plus, des incidents du type coups de foudre ou un brusque dmarrage dune machine tournante forte puissance peuvent causer une chute soudaine et importante de tension. On nommera ce type dincident : creux de tension. Ces perturbations ont bien entendu des consquences nfastes sur les quipements lectriques, lesquelles peuvent aller dun fort chauffement ou dun arrt soudain des machines tournantes jusqu la destruction totale de ces quipements. Plusieurs solutions de dpollution des rseaux lectriques ont t dj proposes dans la littrature. Celles qui rpondent le mieux aux contraintes industrielles daujourdhui sont les compensateurs actifs parallle, srie et combinaison parallle-srie actifs (aussi appels UPQC). Le compensateur actif parallle peut tre utilis pour compenser les courants harmoniques et dsquilibrs de mme que la puissance ractive. Le compensateur srie peut compenser la fois les tensions harmoniques et dsquilibres, et les creux de tension. LUPQC correspond une solution universelle de compensation pour un, plusieurs ou tous les types de perturbation de courant et de tension que lon observe sur le rseau. Progressivement, et de manire endiguer laugmentation des problmes de perturbation sur les rseaux lectriques, des normes de qualit lectrique de plus en plus contraignantes seront imposes aux fournisseurs et aux consommateurs industriels. Face ces nouvelles rglementations mais aussi face la demande de plus en plus spcifique des utilisateurs dnergie lectrique, les compensateurs actifs doivent sadapter en rendant plus complexes leurs structures et leurs mthodes de rgulation. Dautant plus que ces solutions doivent se soumettre aux limites physiques et technologiques inhrentes aux dispositifs de filtrage actif. Dans le cadre de cette problmatique de recherche, une thse de doctorat a t initie au sein des laboratoires du GREEN et du LEPSI en vue doptimiser les performances de contrle des compensateurs actifs tout en tenant compte de lensemble des contraintes et limitations matrielles inhrentes leurs structures. Ainsi, des tudes avances ont t ralises sur les diffrents types de rgulateur existants et sur leur faisabilit en terme dimplantation numrique. Un partenariat industriel est venu complter ce travail de recherche en vue de valider la faisabilit exprimentale de nos investigations. Il sagissait de la socit Socomec (Benfeld, Alsace) qui souhaitait dvelopper un prototype industriel de filtre actif parallle dit universel, capable de sadapter et de compenser le plus largement possible lensemble des perturbations en courant, pouvant apparatre sur le rseau lectrique. De manire atteindre ces objectifs de recherche, ce mmoire sera divis en quatre chapitres : Le premier chapitre sera consacr la description des perturbations en courant et en tension pouvant intervenir dans un rseau lectrique. Les origines, les consquences et les normes inhrentes ces perturbations seront analyses. Dans ce chapitre, nous prsenterons galement les solutions traditionnelles et modernes de dpollution. Finalement, nous dtaillerons les avantages, inconvnients et les perspectives de dveloppement des solutions de type filtre actif, aujourdhui trs utilises dans lindustrie. Dans le second chapitre de ce mmoire, nous traiterons du problme des perturbations en courant et de leur compensation partir dune structure de filtre actif parallle. La structure du
Introduction Gnrale ___________________________________________________________________________ compensateur sera rendue plus complexe par l'introduction d'un filtre du troisime ordre en sortie de londuleur MLI. Nous tudierons diffrentes mthodes didentification et de rgulation (PI, RST) en vue doptimiser les performances de filtrage du compensateur sur la plus large bande de frquence possibles. Nous prsenterons les limites de ces modles de rgulation suite la dgradation des performances de compensation inhrente au dphasage entre les courants identifis et ceux effectivement injects sur le rseau lectrique. Nous proposerons alors une nouvelle mthode de rgulation, base sur un rgulateur RST, permettant de limiter largement ces effets de dphasage et datteindre dexcellentes performances de compensation avec une trs bonne robustesse. La validation des mthodes proposes sera effectue par simulation partir du logiciel SABER. Dans le troisime chapitre, nous traiterons la fois du compensateur actif srie, appliqu ici pour rduire les perturbations en tension, et de la solution combine srie-parallle UPQC, propose pour compenser lensemble des perturbations du rseau lectrique basse tension. Nous tudierons principalement les diffrentes mthodes didentification et de rgulation en vue doptimiser les performances de compensation du filtre actif srie. Ainsi, la mthode de rgulation amliore, base sur le rgulateur RST que nous aurons introduite dans le cas du filtre actif parallle, sera reprise et applique pour rguler ici la tension. Nous observerons de nouveau lefficacit de cette mthode travers les trs bons rsultats de compensation observs, cette fois-ci, dans le cas des perturbations en tension. La validation des mthodes proposes sera effectue par simulation partir du logiciel SABER. Faisant suite notre synthse dtude sur le filtre actif srie, nous prsenterons, dans ce chapitre, la pertinence de notre approche dans la perspective de conception dun compensateur universel (UPQC). Dans le quatrime et dernier chapitre de cette thse, nous aborderons la faisabilit exprimentale dun filtre actif parallle universel partir du cahier des charges de notre partenaire industriel. Dans un premier temps, un modle de la structure du filtre actif parallle universel sera ralis en mode continu. Ce dernier devra tenir compte des nouvelles contraintes de conception et dapplication du filtre actif dans son futur environnement dutilisation. Luniversalit en mode continu du filtre actif propos sera teste travers des simulations effectues partir les logiciels SABER et Matlab- Simulink- Power System Blockset. Ensuite, nous aborderons la discrtisation des rgulateurs impliqus dans la rgulation des courants perturbs, en vue de complter ltude de faisabilit exprimentale du filtre actif parallle universel. Deux mthodes de conception seront proposes et analyses dans ce travail : la conception indirecte en mode discret partir du mode continu et la conception directe en mode discret. Finalement, nous conclurons ce travail en tablissant les limites thoriques et technologiques de notre approche de faisabilit exprimentale.
I. INTRODUCTION
Comme tout gnrateur dnergie lectrique, un rseau de puissance fournit de lnergie aux appareils utilisateurs par lintermdiaire des tensions quil maintient leurs bornes. Il apparat vident que la qualit de cette nergie dpend de celle de la tension au point de livraison. Cette tension subit gnralement beaucoup de perturbations de deux origines distinctes : les perturbations de tension causes par le passage, dans les rseaux lectriques, des courants perturbateurs comme les courants harmoniques, dsquilibrs et ractifs, les perturbations de tension causes par des tensions perturbatrices comme les tensions harmoniques et dsquilibres et les creux de tension.
Dans la premire partie de ce chapitre, nous tudierons les caractristiques gnrales des perturbations lectriques. Ainsi, nous dtaillerons les origines, les consquences matrielles et les limites tolres et imposes par les normes internationales de ces perturbations. Dans une deuxime partie, nous prsenterons et comparerons des solutions de compensation de ces perturbations gnralement proposes dans la littrature. Nous pourrons ainsi adopter la compensation active la plus adquate. La dernire partie de ce chapitre sera consacre la prsentation des diffrents montages des filtres actifs. Nous voquerons ainsi les principes et la mise en uvre de ces structures modernes de compensation.
La qualit de la tension peut tre affecte, soit du fait de certains incidents inhrents la nature physique et aux sujtions lies lexploitation du rseau, soit du fait de certains rcepteurs. Ces dfauts se manifestent sous forme de diffrentes perturbations affectant un ou plusieurs des quatre paramtres prcdemment dfinis. On a donc quatre possibilits distinctes de perturbations : les fluctuations de la frquence 50 Hz : elles sont rares et ne sont observes que lors de circonstances exceptionnelles, par exemple certains dfauts graves du rseau, au niveau de la production ou du transport;
Chapitre I : Perturbations des rseaux lectriques et principe de compensation actif _________________________________________________________________________ les variations de lamplitude : il ne sagit pas des variations lentes de tension qui sont corriges par les transformateurs de rglage en charge, mais de variations rapides de tension ou de creux de tension se prsentant souvent sous forme d-coups brusques. Les creux de tension peuvent tre soit isols, soit plus ou moins rptitifs, de forme rgulire ou non; la modification de la forme donde de la tension : cette onde nest alors plus sinusodale, et peut tre considre comme reprsentable par une onde fondamentale 50 Hz, associe soit des harmoniques de frquence multiple entier de 50 Hz, soit mme parfois des ondes de frquence quelconque; la dissymtrie du systme triphas, que lon appelle dsquilibre.
On peut, en plus, mentionner un type particulier de perturbations difficile classer puisquil concerne tout la fois lamplitude et la forme donde : ce sont les variations transitoires damplitudes dont la dure est infrieure 10 ms. Afin de bien analyser les pollutions des rseaux lectriques en basse tension et, par consquent, de trouver les meilleures mthodes de dpollution, on va distinguer deux types de perturbations, savoir les perturbations de courant et celles de tension. Les courants perturbateurs comme les courants harmoniques, les courants dsquilibrs et la puissance ractive sont majoritairement mis par des charges non linaires, base dlectronique de puissance, et/ou dsquilibres. La puissance ractive peut tre aussi consomme par des charges linaires inductives comme les moteurs asynchrones qui sont largement prsents dans les sites industriels. Les perturbations de tension comme les creux, les dsquilibres et les harmoniques de tension trouvent gnralement leurs origines dans le rseau lectrique lui-mme mais parfois galement dans les charges. Ces types de perturbation ont des effets trs nfastes sur les quipements lectriques. Ces effets peuvent aller des chauffements ou de larrt des machines tournantes jusqu la destruction totale de ces quipements.
Chapitre I : Perturbations des rseaux lectriques et principe de compensation actif _________________________________________________________________________ de tlcommunication, les dfauts de fonctionnement de certains quipements lectriques et le risque dexcitation de rsonance : Lchauffement : Les pertes totales par effet Joule sont la somme de celles du fondamental et des harmoniques :
2 2 I .R = I h.R h =1
(I.1)
avec I le courant total, Ih le courant harmonique de rang h qui reprsente le fondamental pour h=1, et R la rsistance traverse par le courant I. Les harmoniques augmentent aussi les pertes fer (pertes par courants de Foucault). Ils prennent de limportance dans les matriels utilisant les circuits magntiques (moteurs, transformateurs...). Le vieillissement des isolants est souvent d une contrainte en tension conscutive la prsence de la tension harmonique, et donc une augmentation locale du courant de fuite, ou encore lchauffement exagr dans les conducteurs. Linterfrence avec les rseaux de tlcommunication : Le couplage lectromagntique entre les rseaux lectriques et de tlcommunication peut induire dans ces derniers des bruits importants. Dans le cas de rsonances, une partie des rseaux de tlcommunication peut tre rendue inutilisable [DE98, SZ84]. Les dfauts de fonctionnements de certains quipements lectriques : En prsence des harmoniques, la tension (ou le courant ) peut changer plusieurs fois de signe dans une demipriode; par consquent, tout appareil dont le fonctionnement est bas sur le passage par zro des grandeurs lectriques (appareils utilisant la tension comme rfrence) peut tre perturb. Le risque dexcitation de rsonance : Les frquences de rsonance des circuits forms par des inductances du transformateur et des cbles sont normalement leves. Ce nest pas le cas lorsque des batteries de capacit sont raccordes au rseau pour relever le facteur de puissance; les frquences de rsonance peuvent devenir assez faibles, et concider ainsi avec celles des harmoniques engendrs par les convertisseurs statiques. Dans ce cas, il y aura des phnomnes damplification dharmoniques.
Diffrentes grandeurs sont dfinies pour chiffrer ces perturbations. Parmi celles-ci les plus utilises sont : le taux harmonique de rang h : Sh = Ch / C1 (I.2)
ou Ch reprsente la composante harmonique de rang h, C1 reprsente la composante fondamentale, le taux global de distorsion harmonique : THD =
2
C 2h C 21
(I.3)
Chapitre I : Perturbations des rseaux lectriques et principe de compensation actif _________________________________________________________________________ En gnral, les harmoniques pris en compte dans un rseau lectrique sont infrieurs 2500 Hz, ce qui correspond au domaine des perturbations basses frquences au sens de la normalisation. Les harmoniques de frquence plus leve sont fortement attnus par l'effet de peau et par la prsence des inductances de lignes. De plus, les appareils gnrant des harmoniques ont, en grande majorit, un spectre d'mission infrieur 2500 Hz, c'est la raison pour laquelle le domaine d'tude des harmoniques s'tend gnralement de 100 2500 Hz, cest--dire des rangs 2 50 [DE98]. le facteur de puissance :
Normalement, pour un signal sinusodal le facteur de puissance est donn par le rapport entre la puissance active P et la puissance apparente S. Les gnrateurs, les transformateurs, les lignes de transport et les appareils de contrle et de mesure sont dimensionns pour la tension et le courant nominaux. Une faible valeur du facteur de puissance se traduit par une mauvaise utilisation de ces quipements. Dans le cas o il y a des harmoniques, une puissance supplmentaire appele la puissance dformante (D), donne par la relation (I.4), apparat comme le montre le diagramme de Fresnel de la Fig. I.1.
2 D = 3.V1 I h h=2 50
(I.4)
(I.5)
Q g j1 P j
Fig. I.1 : Diagramme de Fresnel des puissances II.1.3 Rglementation Daprs le contrat EMERAUDE dEDF, les deux parties (fournisseur et rcepteur) doivent sengager respecter les normes limitant les perturbations harmoniques. De son ct, EDF sengage ce que les taux individuels de tension harmonique, exprims en pourcentage de la tension fondamentale V1 pour les rseaux HTA (1 50 kV), ne dpassent pas les seuils donns dans le tableau ci-dessous. 10
Chapitre I : Perturbations des rseaux lectriques et principe de compensation actif _________________________________________________________________________ Harmoniques impairs non multiple de 3 Rang Vh/V1 (%) 5 6 7 5 11 3,5 13 3 17 2 19, 23 et 25 1,5 multiple de 3 Rang Vh/V1 (%) 3 5 9 1,5 15 et 21 0,5
Tab. I.1 : Engagement EMERAUDE sur les harmoniques de tension (rseaux HTA) Concernant les rseaux HTB (plus de 50 kV), EDF sengage ne pas dpasser les seuils donns dans le suivant : Harmoniques impairs non multiple de 3 Rang Vh/V1 (%) 5 et 7 2 11 et 13 1,5 17 et 19 1 23 et 25 0,7 Harmoniques pairs multiple de 3 Rang Vh/V1 (%) Rang Vh/V1 (%) 3 2 2 1,5 9 1 4 1 15 et 21 0,5 6 24 0,5 le THD global de tension ne dpassant pas 3 %.
Tab. I.2 : Engagement EMERAUDE sur les harmoniques de tension (rseaux HTB) Les rgles de limitation des courants harmoniques recommandes aux clients par EDF travers le contrat EMERAUDE sont donnes dans le tableau ci-dessous: Rangs impaires 3 5 et 7 9 11 et 13 > 13 Ih/I1 (%) 4 5 2 3 2 Rangs paires 2 4 >4 Ih/I1 (%) 2 1 0,5
Tab. I.3 Limitation EMERAUDE des courants harmoniques Les limitations en tension harmonique que les clients de EDF doivent respecter sont : Pour un harmonique pair : V h 0 ,6%, V1 Pour un harmonique impair : V h 1%, V1 Pour le taux de distorsion global de tension : THD < 1,6 %.
Il est dusage de dire que, dans les installations industrielles, les tensions harmoniques dont le THD est infrieur 5% ne produisent pas deffet notable. Entre 5% et 7% on commence observer des effets, et pour plus de 10% les effets sont quasi certains [DE98]. 11
Chapitre I : Perturbations des rseaux lectriques et principe de compensation actif _________________________________________________________________________ Concernant la puissance ractive, EDF autorise ses clients en consommer, sans tre factur, jusqu 40% de la puissance active absorbe. Cela se traduit, pour des charges linaires, par un facteur de puissance cos j 0,928 ou par un angle de phase j 21,8.
Concernant les dispositifs triphass dlectronique de puissance, principalement les ponts redresseurs, le fonctionnement en prsence de dsquilibre entrane l'apparition de composantes harmoniques non caractristiques, notamment des harmoniques de rang multiple de 3. Lapparition de ces courants harmoniques peut poser des problmes, comme la gnration dune anti-rsonance lors du filtrage de lharmonique de rang 5. Outre les effets classiques des harmoniques, ces frquences non caractristiques peuvent conduire, dans certains cas, au blocage de la commande. La consquence des composantes inverses sur les machines tournantes est la cration dun champ tournant en sens inverse du sens de rotation normal, d'o un couple de freinage parasite et des pertes supplmentaires qui provoquent lchauffement de la machine. Concernant l'effet du dsquilibre homopolaire, il faut signaler le risque d'chauffement du conducteur neutre dans un rseau BT qui, lorsque le conducteur est d'un diamtre trop faible, peut provoquer une rupture du conducteur ou un incendie. II.2.3 Rglementation Le dsquilibre en tension est caractris par le taux de dsquilibre de tension tv donn par le rapport des amplitudes des tensions inverse et directe : tv = Vi Vd (I.6)
Daprs le contrat EMERAUDE, EDF sengage fournir, aux clients raccords aux rseaux HTA (1 50 kV) et HTB (50 130 kV), une tension dont le taux de dsquilibre moyen tm ne dpasse pas 2 %.
12
Chapitre I : Perturbations des rseaux lectriques et principe de compensation actif _________________________________________________________________________ Le fournisseur EDF, travers le contrat EMERAUDE, permet aux clients de salimenter au rseau lectrique sans rserve pour des charges perturbatrices infrieures ou gales : 500 kVA pour une puissance de court-circuit de 40 MVA de 1 50 kV, 4 MVA pour une puissance de court-circuit de 400 MVA de 63 90 kV 15 MVA pour 1500 MVA de puissance de court-circuit 225 kV.
Cependant EDF demandera aux clients qui ne sont pas couverts par ces champs de ne pas provoquer un taux de dsquilibre suprieur 1 %.
Les creux de tension sont caractriss par leur amplitude et par leur dure. Ils sont monophass, biphass ou triphass suivant le nombre de phases concernes. II.3.2 Consquences des creu x de tension Les creux de tension sont susceptibles de perturber le fonctionnement de certaines installations industrielles et tertiaires. En effet, ce type de perturbation peut causer des dgradations de fonctionnement des quipements lectriques qui peuvent aller jusqu la destruction totale de ces quipements. Le Tab. I.4 rsume les consquences nfastes causes par les creux de tension sur quelques matriels industriels et tertiaires sensibles [ME90]. II.3.3 Rglementation Lengagement dEDF, travers le contrat EMERAUDE, se prsente sous forme de seuils, les creux de tension tant caractriss par leur profondeur et leur dure, avec des limites de 30% et de 600 ms, comme le montre la Fig. I.2.
13
Uf
Un 0.7 Un
600 ms
Dure du phnomne
Fig. I.2 Engagement EMERAUDE sur les creux de tension Types dappareils Consquences nfastes
Eclairage Moins de luminosit, extinction et rallumage (lampes arc) Systmes base dlectronique Arrt du dispositif de puissance Dispositifs de protection Ouverture des contacteurs Ralentissements, dcrochage, surintensit au retour de la Moteurs asynchrones tension Moteurs synchrones Perte de synchronisme, dcrochage et arrt du moteur Variateurs de vitesse pour un En mode onduleur : destruction des protections moteur courant continu En mode redresseur : ralentissement de la machine Variateurs de vitesse pour un Ralentissement, dcrochage, surintensit au retour de la tension, destruction ventuelle de matriel au niveau du moteur asynchrone convertisseur Tab. I.4 Consquences des creux de tension sur quelques quipements lectriques sensibles
14
X = L .w = 3 .R
Charge rsistive
Rseau triphas
X= 1 = 3 .R c.w
Fig. I.3 : Montage de Steinmetz pour le rquilibrage III.1.1.2 Compensation de la puiss ance ractive La puissance ractive est majoritairement consomme par les moteurs asynchrones et plus rcemment par des dispositifs base dlectronique de puissance. Diffrentes mthodes de compensation sont utilises pour relever le facteur de puissance. La plus simple consiste placer des batteries de condensateurs en parallle avec le rseau. Linconvnient de cette mthode rside dans le fait que la puissance ractive fournie par les condensateurs est constante et quelle ne sadapte pas lvolution du besoin. III.1.1.3 Compensation des couran ts harmoniques Plusieurs solutions existent pour limiter la propagation et leffet des harmoniques dans les rseaux lectriques : 15
Chapitre I : Perturbations des rseaux lectriques et principe de compensation actif _________________________________________________________________________ laugmentation de la puissance de court-circuit du rseau et lutilisation de convertisseurs peu polluants qui ont pour effet de diminuer la distorsion harmonique, lutilisation de dispositifs de filtrage pour rduire la propagation des harmoniques produits par des charges non linaires.
Le filtrage consiste placer en parallle sur le rseau dalimentation une impdance de valeur trs faible autour de la frquence filtrer et suffisamment importante la frquence fondamentale du rseau. Parmi les dispositifs de filtrage les plus rpandus, on distingue le filtre passif rsonnant (Fig. I.4) et le filtre passif amorti ou passe-haut (Fig. I.5).
L
R L
R
C
Le filtre rsonnant est un filtre trs slectif. Il peut se connecter en parallle avec dautres filtres rsonnants. Le filtre passe-haut compense les harmoniques suprieurs ou gaux sa frquence propre. Il peut se connecter en parallle avec dautres filtres rsonnants. Ces dispositifs sont utiliss pour empcher les courants harmoniques de se propager dans les rseaux lectriques. Ils peuvent aussi tre utiliss pour compenser la puissance ractive. Malgr leur large utilisation dans lindustrie, ces dispositifs peuvent prsenter beaucoup dinconvnients : manque de souplesse sadapter aux variations du rseau et de la charge, quipements volumineux, problmes de rsonance avec limpdance du rseau. III.1.2 Dpollution des tension s perturbatrices Les tensions perturbatrices dans un rseau lectrique basse tension sont principalement les creux de tension, les tensions harmoniques et les tensions dsquilibrs. Ces deux dernires sont gnralement causes par la circulation des courants harmoniques et/ou dsquilibrs. Pour dpolluer les rseaux lectriques de ces deux types de perturbation, on peut limiter la circulation des courants perturbateurs en utilisant les solutions traditionnelles prsentes prcdemment dans le cas des perturbations de courant. Quant aux creux de tension, la solution la plus frquente dans les milieux sensibles (hpitaux, sites industrielsetc.) est dutiliser des groupes lectrognes qui se substituent au rseau lectrique. Mais la limitation de la puissance de ces groupes ainsi que la qualit mdiocre de lnergie lectrique fournie restent un problme.
16
17
~
Iinj
Fig. I.6 : Filtre actif parallle IV.1.2 Mise en uvre Les premiers principes de fonctionnement des filtres actifs parallles ont t prsents dans la littrature ds le dbut des annes 1970 [BI69, SA71]. En 1976, une premire famille de filtre actif parallle a t conue partir donduleurs transistors de puissance commands en MLI [GY76]. Ces premiers dispositifs de puissance taient alors proposs pour la seule compensation des courants harmoniques. Cependant, cette poque, il tait encore difficile de concevoir ce type de systme pour des applications industrielles. En effet, dans ces annes l, il tait presque impossible de trouver, sur le march, des interrupteurs de puissance capables de fonctionner aux frquences de commutation et aux niveaux de puissance exigs par la ralit industrielle. Cette
18
Chapitre I : Perturbations des rseaux lectriques et principe de compensation actif _________________________________________________________________________ barrire technologique sera franchie, ds 1977, lors de la conception dun premier prototype de filtre actif parallle base de thyristors commutations naturelles pour la compensation de courant harmonique [MO77]. Cependant, lapplication des onduleurs base de thyristor a tout de suite pos le problme de la gnration non dsire de composantes injectes sur le rseau la frquence de commutation. La mme raison a galement empch lutilisation de compensateurs statiques parallles thyristors, lesquels avaient t conus pour la compensation conjointe de la puissance ractive et des courants dsquilibrs [GR77]. Au cours des annes 1980, des progrs importants dans le domaine des semiconducteurs ont permis de dvelopper de nouveaux composants de puissance associant hautes frquences de commutation et fortes puissances. Profitant de ces avances, et de lavnement des interrupteurs de puissance du type GTO et IGBT, de nombreux onduleurs de puissance, commands en MLI, ont pu tre conus en vue de rpondre aux contraintes industrielles de conception des filtres actifs parallles [KA83, AK86]. Ainsi, ces derniers ont commenc tre commercialiss et installs travers le monde, et plus spcialement dans les pays les plus industrialiss comme le Japon [AK96, AK97]. Ces premiers prototypes ne compensaient alors que les perturbations harmoniques de courant. Suite ces premiers dveloppements, dautres types de filtre actif parallle ont pu tre raliss pour compenser la fois la puissance ractive, et/ou les harmoniques et les dsquilibres de courant [CH2000]. Actuellement, les filtres actifs parallles sont essentiellement installs par les consommateurs industriels. Lvolution future de ces dispositifs de puissance pourrait autoriser le fournisseur dnergie prendre un rle plus important, en lui permettant de les installer lui-mme. Cette approche permettrait damortir la propagation des harmoniques causes par la rsonance, laquelle peut tre observe entre les inductances du rseau et les batteries de condensateur installes pour amliorer le facteur de puissance [AK96, AK97]. De mme, des filtres actifs parallles installs par le fournisseur auraient aussi pour objectif de rduire la distorsion harmonique de tension en amont, ct rseau lectrique.
Chapitre I : Perturbations des rseaux lectriques et principe de compensation actif _________________________________________________________________________ IV.2.2 Mise en uvre Le filtre actif Srie (F.A.S) est une solution adapte la compensation des tensions perturbatrices, harmoniques, dsquilibres et des creux de tension. Ces perturbations trouvent gnralement leurs origines dans le rseau lui-mme mais peuvent parfois tre provoques par les charges elles-mmes. En 1976, une famille de filtres actifs sries avec des onduleurs transistors contrls en MLI a t prsente [GY76]. Un compensateur statique pour le rquilibrage de la tension du rseau lectrique a t propos en 1985 [MU85]. La gnration, par londuleur thyristors, des composantes non dsires, ayant une faible frquence de commutation, a dcourag lavancement pratique de cette solution. En 1990, un compensateur de tension base donduleur de tension triphas a t propos pour compenser le dsquilibre de tension de rseau lectrique [MOR90]. Depuis, beaucoup de travaux ont t consacrs lamlioration de la qualit de compensation du dsquilibre de la tension du rseau [CA92, CA94]. Dautres articles ont introduit le problme des creux de tension et la solution de compensation par le compensateur actif srie [BH94, HAD97, TO99].
Iinj
F.A.S
F.A.P
Chapitre I : Perturbations des rseaux lectriques et principe de compensation actif _________________________________________________________________________ IV.4.1 Le filtre actif srie ave c des filtres passifs parallles Le rle du filtre actif srie dans ce cas est dempcher les courants harmoniques de circuler vers le rseau et de les obliger passer par les filtres passifs raccords leurs frquences comme le montre la Fig. I.9 [WA92].
Fig. I.9 : Filtre actif srie et filtre passif parallle IV.4.2 Le filtre actif srie conn ect en srie avec des filtres passifs parallles Le principe de fonctionnement de cette configuration, prsente en Fig. I.10, est la mme que la prcdente avec lavantage de rduire encore le dimensionnement du filtre actif srie car le courant qui le traverse est plus faible [WA92]. De plus, le filtre actif srie est labri dun ventuel courtcircuit de la charge.
~
Filtre passif parallle
Fig. I.10 : Filtre actif srie connect en srie avec un filtre passif parallle IV.4.3 Le filtre actif parallle avec un filtre passif parallle Le rle du filtre actif parallle dans cette configuration, montre en Fig. I.11, est la compensation des courants harmoniques basses frquences mis par la charge polluante. Le filtre passif accord sur une frquence leve, limine les harmoniques hautes frquences y compris ceux crs par le filtre actif parallle. Ce type de filtrage a dj t appliqu la compensation des courants harmoniques mis par un cycloconvertisseur de forte puissance [WA92].
21
~
Filtre passif parallle
22
V. CONCLUSION
Dans ce chapitre, nous avons prsent diffrents types de perturbations affectant londe de tension du rseau lectrique. Comme nous avons pu le constater, les harmoniques et les dsquilibres de courant et de tension, la puissance ractive et les creux de tension ont des effets nfastes sur les quipements lectriques. Ces effets peuvent aller des chauffements et de la dgradation du fonctionnement jusqu la destruction totale de ces quipements. Plusieurs solutions traditionnelles et modernes de dpollution ont t prsentes. La solution classique base de filtres passifs est souvent pnalise en termes dencombrement et de rsonance. De plus, les filtres passifs ne peuvent pas sadapter lvolution du rseau et aux charges polluantes. Rcemment, en plus du filtrage des harmoniques, les filtres actifs parallles et sries, et leur combinaison, sont tudis pour la compensation de tous les types de perturbation susceptible dapparatre dans un rseau lectrique basse tension. En effet, profitant des progrs raliss dans le domaine de llectronique de puissance, ces solutions peu encombrantes noccasionnent aucune rsonance avec les lments passifs du rseau et font preuve dune grande flexibilit face lvolution du rseau lectrique et de la charge polluante. Le filtre actif parallle peut tre install pour compenser toutes les perturbations de courant comme les harmoniques, les dsquilibres et la puissance ractive. Le filtre actif srie peut tre install pour compenser toutes les perturbations de tension comme les harmoniques, les dsquilibres et les creux de tension. La combinaison parallle-srie actifs est la solution universelle pour compenser toutes les perturbations en courant et en tension. Dautres solutions de dpollution consistent associer des filtres actifs et passifs pour trouver un compromis entre lefficacit de compensation et le cot. Ces solutions peuvent tre des solutions futures trs prometteuses. Les filtres hybrides tant conus essentiellement pour dpolluer les rseaux lectriques des courants harmoniques, notre objectif de recherche portant sur la dpollution de toutes sortes de perturbations, seules les solutions actives de filtrage seront analyses dans ce mmoire. Cest ainsi que dans lobjectif damliorer la qualit de lnergie lectrique qui doit tre conforme aux nouvelles contraintes normatives, nous tudierons, dans les chapitres suivants, des mthodes de rgulation avances que nous appliquerons dans le cas des structures de filtres actifs parallle, srie et combin.
23
24
25
26
I. INTRODUCTION
Les objectifs de recherche que nous nous sommes fixs dans ce chapitre devront nous permettre damliorer les performances de compensation du filtre actif parallle dans le cas de conditions svres, ou mme extrmes, de fonctionnement des rseaux lectriques. Les perturbations de courant que nous tudierons seront de type aussi bien harmonique, ractif, que dsquilibr. Dans le but datteindre une partie de ces objectifs, nous aborderons une structure de filtre actif parallle associant un filtre du troisime ordre, galement appel filtre en T, en sortie de londuleur. Cette solution permet dempcher, sans entraner deffets secondaires, les composantes dues aux commutations de se propager sur le rseau lectrique. Cependant, la structure avec un filtre en T entrane une plus grande complexit de mise en uvre du filtre actif parallle. Il nous a paru important daborder ces problmes qui restent encore peu abords dans la littrature scientifique. Finalement, ce chapitre sera dcompos en quatre parties. Dans une premire partie, nous tudierons la structure gnrale du filtre actif parallle, laquelle a t divise en deux : la partie puissance et la partie contrle-commande. Dans la partie puissance, nous introduirons les trois principaux blocs de cette structure, savoir londuleur de tension, llment de stockage dnergie et le filtre de sortie. Dans la partie contrle-commande, nous aborderons les cinq principaux blocs, soit lidentification des courants perturbs, la rgulation de la tension continue et du courant du filtre actif parallle, la commande de londuleur et le systme base de PLL. Le systme base de PLL sera intgr dans la mthode didentification des courants perturbs afin de gnraliser lapplication de cette dernire dans les cas les plus dfavorables de la tension du rseau. Dans la seconde partie de ce chapitre, nous tudierons lensemble des contraintes et objectifs de rgulation du filtre actif parallle. Nous aborderons, plus particulirement, les mthodes de rgulation du courant dans le cas dun rgulateur de forme RST. Ensuite, nous analyserons les amliorations et les dgradations, en termes de compensation, observes partir de ce mode de rgulation. Nous identifierons, plus spcifiquement, un problme de dphasage entre les courants de rfrence et ceux injects, entranant une diminution des performances du filtre actif. Finalement, en vue de rsorber ces difficults, nous proposerons une toute nouvelle mthode de rgulation, dveloppe partir dune structure RST classique et que lon nommera : RSTam. Les troisime et quatrime parties de ce chapitre viendront corroborer les mthodes de rgulation proposes dans cette thse en vue dtudier la faisabilit de la commande du filtre actif parallle. Dans un premier temps, nous approcherons des contraintes conomiques de conception du filtre actif parallle travers lestimation de la puissance apparente de ce dernier dans les trois cas habituels de compensation. Puis, dans un second temps, nous raliserons, par le logiciel SABER, plusieurs cas de simulation comportementale du filtre actif parallle suivant la commande et les perturbations choisies.
27
Partie Puissance
Filtre de sortie Onduleur Elment de stockage
Courant de la charge
Tension du rseau
Systme PLL
Partie ContrleCommande
Chapitre II : Filtre actif parallle, stratgie de commande, rgulation et dimensionnement ___________________________________________________________________________ se fait par lintermdiaire dun condensateur Cdc de tension Vdc. Le filtre de sortie est un filtre passif habituellement du premier ordre (Lf,Rf) employ pour connecter londuleur de tension au rseau lectrique.
T2
T3
Cdc
Vdc
T5
T6
Fig. II.2 : Onduleur de tension triphas Cette structure du filtre actif parallle ne permet pas la fermeture simultane des semiconducteurs dun mme bras sous peine de court-circuiter le condensateur de stockage. Par contre, ils peuvent tre tous les deux ouverts (pendant un temps mort). La continuit des courants est alors assure par la mise en conduction dune des diodes dun mme bras. En pratique, nous commandons les deux semi-conducteurs dun mme bras de faon complmentaire : la conduction de lun entrane le blocage de lautre. En ralit, le mode, o les semi-conducteurs d'un mme bras sont tous les deux ferms, n'existe que durant les commutations. Afin d'viter un court-circuit cause du dlai de blocage des interrupteurs, il faut insrer sur un mme bras, un temps d'attente, galement appel temps mort, entre la commande de blocage d'un interrupteur et la commande d'amorage de l'autre [WA92]. Avec l'hypothse des commutations instantanes, ce mode de fonctionnement ne sera pas pris en compte et par consquent, aucun risque de court-circuiter le condensateur nest craindre. II.2.1.2 Tension fournie par lond uleur Louverture et la fermeture des interrupteurs de londuleur de la Fig. II.2 dpendent de ltat des signaux de commande (S1,S2,,S3,), comme dfini ci-dessous : 1 T ferm et T ouvert S1 = 0 T1 ouvert et T4 ferm 1 4 1 T ferm et T ouvert S2 = 0 T2 ouvert et T5 ferm 2 5 S3 = 1 T3 ferm et T6 ouvert 0 T ouvert et T ferm 3 6 Ainsi, on pourra exprimer huit cas possibles de tension de sortie du filtre actif Vf (rfres au neutre n de la source), comme le montre le Tab. II.1 [XU94]. II.2.1.3 Reprsentation vectoriell e Dans le plan biphas (a,b), en considrant V f le vecteur correspondant aux tensions de londuleur, les huit cas possibles du vecteur V f sont donns par la Fig. II.3. (II.1)
29
Chapitre II : Filtre actif parallle, stratgie de commande, rgulation et dimensionnement ___________________________________________________________________________ N du cas 0 1 2 3 4 5 6 7 S3 0 0 0 0 1 1 1 1 S2 0 0 1 1 0 0 1 1 S1 0 1 0 1 0 1 0 1 Vf3 0 -Vdc/3 -Vdc/3 -2Vdc/3 2Vdc/3 Vdc/3 Vdc/3 0 Vf2 0 -Vdc/3 2Vdc/3 Vdc/3 -Vc/3 -2Vdc/3 Vdc/3 0 Vf1 0 2Vdc/3 -Vdc/3 Vdc/3 -Vdc/3 Vdc/3 -2Vdc/3 0
Vf ( 2 )
Vf ( 3 )
Vf ( 6 )
Vf ( 0 )
Vf (7 )
V f*
Vf ( 1 )
Vf ( 4 )
Vf ( 5 )
Fig. II.3 : Reprsentation vectorielle des tensions gnres par londuleur O V f reprsente la tension de rfrence que doit produire londuleur pour pouvoir crer les courants perturbateurs identifis. Ceci signifie que londuleur nest capable de fournir des tensions gales aux tensions de rfrence que si le vecteur form par ces derniers reste lintrieur de lhexagone montr dans la Fig. II.3 [AK90, XU94]. II.2.2 Systme de stockage d nergie Le stockage de lnergie du ct continu se fait souvent par un systme de stockage capacitif reprsent par un condensateur Cdc qui joue le rle dune source de tension continue Vdc, comme le montre la Fig. II.2. Le choix des paramtres du systme de stockage (Vdc et Cdc) se rpercute sur la dynamique et sur la qualit de compensation du filtre actif parallle. En effet, une tension Vdc leve amliore la dynamique du filtre actif. De plus, les ondulations de la tension continue Vdc, causes par les courants engendrs par le filtre actif et limites par le choix de Cdc, peuvent dgrader la qualit de compensation du filtre actif parallle. Ces fluctuations sont dautant plus importantes que lamplitude du courant du filtre est grande et que sa frquence est faible. Pour cette raison, nous 30
*
(II.2)
avec Vs la tension simple du rseau, Ih le courant harmonique du rang h et a langle dallumage des thyristors du pont Graetz. Puisque laugmentation de la valeur de la tension continue amliore la commandabilit du filtre actif (Cf. Tab. II.1 et Fig. II.3) et sachant que le choix de cette tension se rpercute en grande partie sur le choix des interrupteurs, la tension continue Vdc doit tre choisie comme la plus grande tension respectant les contraintes des interrupteurs. La deuxime mthode, plus simple, se base sur la mesure du courant harmonique Ih du rang le plus faible. La capacit Cdc se calcule de la faon suivante [ALA2000-1]:
Cdc = Ih e .Vdc .w h
(II.3)
avec wh, la pulsation la plus faible des harmoniques compenser. II.2.3 Filtre de sortie Le filtre de sortie est un filtre passif utilis pour connecter londuleur de tension au rseau lectrique. Le filtre de sortie est dimensionn pour satisfaire les deux critres suivants : assurer la dynamique du courant :
d d I h _ ch = I inj dt dt
(II.4)
empcher les composantes dues aux commutations de se propager sur le rseau lectrique.
Deux types du filtre de sortie peuvent tre employs : un filtre de sortie du premier ordre et un filtre de sortie du troisime ordre.
II.2.3.1 Filtre du premier ordre Ce type de filtre est le plus utilis dans la littrature [AK90, WA92, XU94]. Il est compos dune inductance Lf de rsistance interne Rf , comme le montre la Fig. II.2. Un filtre de ce type ne permet pas de satisfaire simultanment les deux critres de dimensionnement du filtre de sortie. En effet, seule une valeur relativement faible de Lf peut raliser une bonne dynamique du filtre actif en satisfaisant lgalit (II.4). Malheureusement, une valeur faible de Lf permet la majorit des 31
Cf es Vs Rf Vf
Onduleur de tension
Fig. II.4 : Schma quivalent par phase du filtre de sortie du troisime ordre en T Les quations qui modlisent le filtre de sortie sont :
I (s) =
inj
(II.5)
avec Vf la tension de sortie de londuleur, B1(s)/A(s) la fonction de transfert du filtre de sortie correspondant au systme dorigine et B2(s)/A(s) la fonction de transfert correspondant au modle de perturbation. Ces perturbations sont causes par la tension du rseau lectrique e(s) que lon considre dsormais gale la tension de raccordement Vs pour les rseaux lectriques de forte puissance de court-circuit.
A(s) = a1 s 3 + a2 s 2 + a3 s + a4 B1(s) = b11 s + b12 B2 (s) = -( b21 s 2 + b22 s + b23 )
(II.6)
avec
32
(II.7)
(II.8)
La frquence de rsonance fcp du systme dorigine est donne par la relation II.(9) :
f cp = 2.p . 1 ( L f 1 . L f 2 )C f Lf 1 + Lf 2
(II.9)
33
identification partir de la dtection du courant de la charge polluante, identification partir de la dtection du courant de la source, identification partir de la dtection de la tension de la source.
La premire mthode est la plus approprie au filtre actif parallle install par le consommateur pour compenser les courants perturbateurs causs par ses charges polluantes [AK96]. Cest pourquoi cette mthode de dtection sera employe dans ce travail. II.3.2 La mthode didentific ation des courants perturbs II.3.2.1 Gnralits sur les mtho des didentification Les diffrentes mthodes didentification de courant perturbateur peuvent tre regroupes en deux familles dapproche. La premire utilise la transforme de Fourier rapide dans le domaine frquentiel, pour extraire les harmoniques du courant. Cette mthode est bien adapte aux charges o le contenu harmonique varie lentement. Elle donne aussi lavantage de slectionner les harmoniques individuellement et de ne choisir de compenser que les plus prpondrants. Il est noter que cette mthode ncessite une grande puissance de calcul afin de raliser, en temps rel, toutes les transformations ncessaires pour extraire les harmoniques [NAK89]. La deuxime famille est base sur le calcul des puissances instantanes dans le domaine temporel. Certaines de ces mthodes se basent sur le calcul des puissances harmoniques de la charge non linaire [BENCH98]. Dautres peuvent tre utilises pour compenser la fois les courants harmoniques et la puissance ractive, en se basant sur la soustraction de la partie fondamentale active du courant total [BENCH98]. Rcemment, des nouvelles mthodes didentification ont t prsentes pour donner le choix de compenser un, plusieurs ou voire mme tous les types de courants perturbateurs. En effet, en se basant sur la rgulation de la tension continue et sur celles du rseau lectrique aux points de raccordement, nous pouvons compenser la fois tous les courants perturbateurs, tout en offrant la possibilit de rguler la tension de la charge [CH2000]. Cette mthode, qui ne peut tre implante que numriquement, ne garantit pas une compensation parfaite de la puissance ractive, de mme que la rgulation de tension nassure pas toujours une bonne qualit la tension de la charge. Une autre mthode, appele la mthode de dtection synchrone et reposant sur la transforme de Park, a t propose [GU2000]. Cette mthode se base essentiellement sur le calcul de la pulsation
34
II.[Link] Gnralits sur les puissances instantanes En prsence des harmoniques, la puissance apparente est compose de trois parties : active (P), ractive (Q) et dformante (D) comme le montre la relation II.10 (Cf. Fig. I.1). S = P2 + Q2 + D2
(II.10)
Cette mthode exploite la transformation a-b pour obtenir les puissances relles et imaginaires. Notons par (Va,Vb) et (Ia, Ib) les composantes orthogonales du repre a-b associes respectivement aux tensions de raccordement du filtre actif parallle (Vs) et aux courants absorbs par les charges polluantes (Is). La transformation a-b-0 triphase permet d'crire, la relation des tensions suivante :
2 . 3 1 2 1 0 1 2 1 2 3 2 1 2 Vs 1 1 - Vs 2 2 3 Vs 3 2
V0 Va = Vb
(II.11)
I0 Ia = I b
(II.12)
Les composantes avec lindice (0) reprsentent les squences homopolaires du systme triphas de courant et de tension. La puissance active instantane, note P(t), est dfinie par la relation suivante :
P( t ) = Vs 1 I s 1 + Vs 2 I s 2 + Vs 3 I s 3 = p( t ) + p0 ( t ) p( t ) = Va Ia + Vb I b p0 ( t ) = V0 I 0
(II.13)
35
p V sa V sb I a . q = V sb V sa I b
(II.15)
Dans le cas gnral, chacune des puissances p et q comporte une partie continue et une partie alternative, ce qui nous permet d'crire lexpression ci-dessous :
p p = p + ~ ~ q = q + q
avec :
(II.16)
p une puissance continue lie la composante fondamentale active du courant et de la tension, q une puissance continue lie la composante fondamentale ractive du courant et de la tension, ~ et q des puissances alternatives lies la somme des composantes perturbatrices du ~ p courant et de la tension.
Afin disoler les puissances actives et ractives conventionnelles, il faut connatre avec prcision les frquences des pulsations des puissances instantanes formes partir de la relation (II.15). Pour cela, nous baserons notre tude sur lanalyse frquentielle des puissances instantanes donnes dans diffrents cas de courant et tension ( Annexe-I). Il est noter quon considre dsormais que le rseau tudi est compos de trois fils ce qui empche les composantes homopolaires dy circuler. II.[Link] Sparation des puissances perturbatrices Aprs avoir identifi les pulsations des puissances instantanes, le filtre de puissance charg disoler les puissances active et ractive conventionnelles peut tre dimensionn. Un circuit constitu d'un filtre passe bas avec un soustracteur peut tre employ, comme le prsente la Fig. II.6.
36
X -
~ X
Fig. II.6 : Schma reprsentant le principe de sparation des puissances Lordre de ce filtre passe-bas dfinit la dynamique et lefficacit de la mthode didentification. Des filtres de puissance du quatrime ou cinquime ordre ont t proposs [AK86]. Dans notre tude, nous avons choisi un filtre passe-bas du deuxime ordre en vue de simplifier lapproche dimplantation numrique de ce dernier. En effet, un ordre plus lev entranerait des temps de calcul plus longs ce qui peut tre prjudiciable dans notre tude. La relation suivante donne lexpression gnrale dun filtre passe-bas du deuxime ordre :
w0 2 2 s + 2xw 0 s + w 0
2
(II.17)
La frquence de coupure, f0=w0/2p, est choisie pour que le filtre de puissance puisse bloquer ~ toute composante perturbatrice des puissances instantanes ( X ). Elle doit aussi permettre aussi le passage des composantes continues reprsentant les puissances active et ractive la frquence fondamentale ( X ). Cette frquence est donc choisie selon le type de la charge, soit :
60 Hz pour un courant de charge quilibr avec un temps de rponse du filtre de 20 ms, 20 Hz pour un courant de charge dsquilibr avec un temps de rponse du filtre de 60 ms.
II.[Link] Calcul des courants perturbateurs En inversant la relation (II.15), nous pouvons recalculer les courants dans le repre a-b comme le montre lquation (II.18) :
Vsa Ia 1 I = 2 2 b V sa + V sb Vsb
- Vsb p Vsa q
(II.18)
En considrant les quations (II.16) et (II.18), nous pouvons sparer le courant dans le repre (a-b) en trois composantes, active et ractive la frquence fondamentale et les harmoniques. Ceci conduit :
Ia 1 Vsa - Vsb p 1 Vsa - Vsb 0 1 Vsa - Vsb ~ p I = V + V + V ~ Vs Vs q D b D b D b b 1 s4 24a 0 1 s4 24a 3 1 s4 2Vsa4 q 4 4 44 4 4 4 4 3 4 4 44 3
courant actif courant ractif courants harmoniques
(II.19)
avec D = Va2 + Vb2 suppos constant dans l'hypothse dune tension sinusodale quilibre du rseau lectrique. Cette expression montre donc que l'identification des diffrentes composantes du courant dans le repre (a-b) revient sparer des termes continus des puissances relle et imaginaire
37
(II.20)
avec Ipa, Ipb courants perturbateurs calculs dans le repre a-b partir des courants ractifs et harmoniques de la relation (II.19). II.[Link] Algorithme d'identification Finalement, lalgorithme didentification que nous avons dfini peut tre reprsent par le synopsis de la Fig. II.7. Par souci de simplification, nous dcrirons les matrices des relations (II.11 et II.12), aprs avoir enlev les composantes homopolaires, et (II.20) respectivement par les blocs nots C32 et C23.
Vs1,2,3
Systme PLL
Is1,2,3
Vdc-ref
-
Vdc
+ Rgulateur de la tension continue Rgulateur du courant du F.A.P
qd
Vd C32
Isa Isb
* *
+ +
Filtre Passe-Haut
% p
2
pc
+ 2
V da + V db
Vdb
Vda
* *
Vdb
+ Iah -
C23
Vd1,2,3
Vda
Isb Isa
*
Vsb
+ +
Filtre Passe-Haut
% q
2
Iref-1,2,3
* *
C23
V da + V db
2
Vda
+ +
Ibh
Vda
Vda
*
Vdb
+ +
db
Fig. II.7 : Schma de principe de la mthode didentification des puissances instantanes En plus de lalgorithme gnral de la mthode didentification, nous reprsentons sur la Fig. II.7, trois des blocs appartenant la partie de contrle-commande et que nous aborderons dans ce qui suit : le bloc rgulation de tension continue, le bloc rgulation du courant du filtre actif et le bloc systme PLL. 38
+-
RST
wd
1/s
qd
Vsd Vsq Filtre passe-bas Vsa (qd ) Vs1 (qs) Vs2 (qs) Vs3 (qs) Vsb (qd )
P( -q d )
Vd
Sin Cos
Vd 1 = Vd .Sin ( wt + d d )
Ct32
-1
C23 Vd
Vd 2 = Vd .Sin ( wt + d d - 2p / 3 )
Vd 3 = Vd .Sin ( wt + d d + 2p / 3 )
Fig. II.8 : Schma global du systme base de PLL II.3.3.2 Principe de fonctionneme nt Le principe de la PLL est bas sur lemploi dun rgulateur RST spcifique, que nous dvelopperons dans la troisime partie de ce chapitre, garantissant un bon niveau de fiabilit et une rponse rapide [BR99]. Son fonctionnement se base sur la transformation de Park dans le repre d,q. A partir de la Fig. II .8, les tensions Vs123(qs), mesures au point de raccordement du filtre actif parallle, subissent dans un premier temps la transformation a-b. Les tensions ainsi obtenues sont exprimes dans le repre de Park par une rotation P( -q d ). Langle de cette rotation, issu de lintgration de lestimation de la pulsation w d , dtermine par le rgulateur RST, reprsente langle de la composante directe de tension dd. Nous souhaitons que langle de rotation estim q soit gal
d
langle rel de la tension directe q d . Cela peut tre ralis en choisissant la valeur de Vd-ref qui satisfait cette condition. La PLL sera verrouille lorsque langle estim q d sera gal q d .
39
(II.21)
V=Vd qs=qd, avec qd=wdt+dd et wd=ws (pulsation de la tension du rseau), dd (langle de la composante directe de la tension du rseau).
On calcule les composantes de tension directe et en quadrature dans le repre de Park, soit respectivement Vsd et Vsq, en appliquant dabord la transformation (a-b), puis celle de Park (d,q) avec un angle de rotation q d . Cette approche de calcul est regroupe dans la relation suivante :
sin( q d - q d ) t Vsd sin( Dq d ) V = [C32 ] .P( -q d ).Vs 123 ( q d ) = 3Vd ) 3Vd - cos( Dq d ) sq - cos( q d - q d
(II.22)
Lasservissement de phase sera ralis lorsque Dqd = 0. Dans ce cas, langle de la tension du rseau qd et langle estim q d sont gaux et nous pouvons ainsi crire :
(II.23)
De lquation (II.23), on remarque que pour obtenir, dans le mme temps, langle et lamplitude de la composante directe de tension (dd,Vd), la tension directe Vsd doit tre nulle. Cela implique que la tension directe de rfrence Vd-ref soit nulle. Aprs avoir verrouill, nous utilisons ces paramtres (dd,Vd) pour reconstituer la composante directe triphase de tension (Vd123) qui serait gale la tension du rseau Vs123, si celle-ci tait saine. Si la tension du rseau est dsquilibre et pollue, les composantes de tension issues de la transformation de Park le seront galement, et de ce fait le fonctionnement de la PLL sera perturb. Pour cela, lors du choix des paramtres du rgulateur, il faut tenir compte de ces perturbations afin que le systme base de PLL fonctionne avec une bonne prcision dans tous les cas de figure. II.3.4 Rgulation de la tensio n continue La tension moyenne Vdc aux bornes du condensateur doit tre maintenue une valeur fixe. La principale cause susceptible de la modifier sont les pertes dans le filtre actif (interrupteurs et filtre de sortie). La rgulation de la tension moyenne aux bornes du condensateur de stockage dnergie doit se faire par ladjonction des courants fondamentaux actifs dans les courants de rfrence. La ~ sortie du rgulateur Pc sajoute, un signe prs, la puissance active perturbatrice P et donne lieu un courant fondamental actif corrigeant ainsi Vdc (Cf . Fig. II.7). La puissance Pc reprsente la 40
(II.24)
Notons que la relation (II.24) est non linaire. Pour des faibles variations de la tension Vdc autour de sa rfrence Vdc-ref, elle peut tre linarise travers les relations suivantes : d Pc = Cdc .Vdc - ref ( Vdc ) dt (II.25) V ( s ) = Pc ( s ) dc Vdc - ref Cdc s A partir de la relation (II.25), et en prenant en compte le rgulateur proportionnel (Kc), la boucle de rgulation de la tension continue peut tre reprsente par le schma de la Fig. II.11. Le choix du paramtre Kc aura pour objectif dobtenir un temps de rponse minimal afin de ne pas nuire la dynamique du filtre actif.
Vdc-ref
+-
Kc
Pc
Vdc
Fig. II.9 : Boucle de rgulation de la tension continue II.3.5 La commande de lond uleur Le but de la commande de londuleur est de permettre la meilleure reproduction des courants perturbs de rfrence, travers les ordres de commande appliqus aux drivers interrupteurs de puissance. Les deux principales familles de commande des convertisseurs statiques sont : la commande par hystrsis, la commande par modulation de largeur dimpulsion (MLI). II.3.5.1 Commande par hystrsis La commande par hystrsis, appele aussi commande en tout ou rien, est une commande non linaire qui utilise lerreur existant entre le courant de rfrence et le courant produit par londuleur [BENCH98]. Lerreur est compare un gabarit appel bande dhystrsis. Ds que lerreur atteint
41
Iref
Iinj
Fig. II.10 : Principe de commande des courants par hystrsis II.3.5.2 Commande par modulati on de largeur dimpulsion Afin de contourner les problmes prcdents, nous introduirons une deuxime famille de commande de londuleur : la commande par modulation de largeur dimpulsion (MLI). La technique de commande par MLI rsout le problme de la matrise de la frquence de commutation en fonctionnant avec une frquence fixe facile filtrer en aval de londuleur. La plus simple et la plus connue des modulations de largeur dimpulsion est sans doute la MLI chantillonnage naturel, dite MLI intersective. Cette technique de commande met en uvre dabord un rgulateur qui dtermine la tension de rfrence de londuleur (modulatrice) partir de lcart entre le courant mesur et sa rfrence. Cette dernire est ensuite compare avec un signal triangulaire (porteuse frquence leve fixant la frquence de commutation). La sortie du comparateur fournit lordre de commande des interrupteurs. Le schma de principe est donn par la Fig. II.11.
Iref
Iinj
Rgulateur
Fig. II.11 : Principe de commande des courants par MLI Dautres techniques de MLI existent galement dans la littrature comme la MLI chantillonnage rgulier o on peut distinguer deux mthodes : la MLI chantillonnage rgulier symtrique o la rfrence est chantillonne chaque priode de la porteuse,
42
(II.26)
avec Lf = Lf1 + Lf2. r Notons par D I f la diffrence entre le courant de rfrence et le courant mesur partir de la relation suivante : r r r D I f = I ref - I inj Des quations (II.26) et (II.27), nous obtenons lexpression ci-dessous :
Lf
r r d r d r D I f = ( Vs + L f I ref ) - V f dt dt
(II.27)
(II.28)
Le premier terme de la partie droite de la relation (II.28) peut tre dfini comme tension de rfrence (Vf-ref), ce qui nous donne lexpression suivante :
r r d r V f - ref = Vs + L f I ref dt
(II.29)
r r Lcart entre V f - ref et V f produit alors une erreur sur le courant. Selon la relation (II.29), la tension de rfrence est compose de deux termes frquences diffrentes. Le premier reprsente la 43
Iref
Iinj
Rgulateur
Vf-ref
G(s) Onduleur
Vf
1/Lf s
Iinj
Fig. II.12 : Schma de la rgulation des courants du F.A.P Dans ce schma, G(s) reprsente londuleur qui peut tre modlis par la relation suivante [AK97] :
1 G( s ) = K 1 + t s V K = dc 2V p
(II.30)
avec Vdc la tension ct continu de londuleur, Vp lamplitude de la porteuse triangulaire et (t) reprsentant le retard caus par le calcul des courants perturbateurs. Dans un premier temps, nous pourrons considrer le gain K comme tant gal 1 dans le cas o Vdc = 2Vp. De plus, nous introduirons, dans la boucle de rgulation du courant, un retard (t) correspondant une priode dchantillonnage (Te) : t = Te. Ce retard correspond au temps de calcul, d une implantation numrique de la partie contrle-commande du filtre actif parallle. En effet, lors dune implantation numrique, les ordres de commande envoys sur londuleur ne peuvent tre pris en compte quaprs une premire priode dchantillonnage ncessaire llaboration de ces derniers. Nous verrons plus spcifiquement lensemble des contraintes dimplantation du filtre actif parallle dans le Chapitre IV de ce prsent mmoire. II.3.6.2 Mthode de rgulation pa r PI Un rgulateur doit satisfaire aux objectifs gnraux de la rgulation ainsi quaux contraintes lies au rejet des perturbations. Dans le cas dun rgulateur aussi simple que le rgulateur PI, et cause du manque de degr de libert que nous impose ce dernier, seule la poursuite des rfrences pour un point de fonctionnement fixe pourra tre envisage. A cela sajoute le problme du compromis trouver dans le cas dune structure de filtre actif parallle avec un filtre de sortie du premier ordre (choix de linductance) [ALA2000-1, ALA2001]. En effet, un rgulateur PI ne peut pas assurer, avec les seuls paramtres de correction KP et KI , les objectifs de rgulation en boucle ferme pour un systme dorigine dordre suprieur un. Par consquent, pour rsoudre des problmes plus complexes de contrle-commande de filtre actif parallle, il est souhaitable daborder des rgulateurs plus avancs, dont le nombre de paramtres se dtermine en fonction de lordre du systme dorigine. Le choix dune mthode de rgulation base sur un rgulateur RST peut rpondre ces critres de rgulation plus svre .
44
II.4 Conclusion
Finalement, nous avons pu prsenter et dfinir, dans la premire partie de ce chapitre, la plupart des lments constituant la structure gnrale du filtre actif parallle. Cette structure a dabord t partage en deux parties, lune dite partie puissance, et lautre dite partie contrle-commande. Nous avons exprim les critres qui nous ont guid dans le choix de chacun des lments des parties concernes. Ainsi, nous avons pu fixer le choix des lments de la partie contrle-commande, tels que la mthode didentification des courants perturbs, la commande de londuleur et la rgulation de la tension continue, de mme que ceux de la partie puissance, comme londuleur, llment de stockage et le filtre de sortie. Par contre, les lments de la partie contrle-commande, tels que le systme base de PLL et la rgulation du courant du filtre actif, nont pu tre compltement dfinis. En effet, ces derniers lments ncessitaient la connaissance de la mthode de rgulation qui leur sera applique en vue de rpondre aux exigences de compensation du filtre actif parallle. Cette tude fera lobjet de la troisime partie de ce chapitre, intitule tude des boucles de rgulation.
45
T(s)
+R(s)
1/S(s)
B(s)
++
1/A(s)
Sortie Y
+ +
dex
46
Y(s) =
(T.B)(s) S(s) (B. R)(s) Y (s) + d (s) + d (s) (S. A + B. R)(s) c (S. A + B. R)(s) in (S. A + B. R)(s) ex
(II.31)
Le principe du placement de ples consiste spcifier un polynme de stabilit arbitraire D(s) et calculer les polynmes S(s) et R(s) de sorte que lon ait :
D(s) = A(s). S(s) + B(s). R(s)
(II.32)
Cette quation est dite quation diophantine, ou quation de Bezout. Supposons que la rfrence Yc et les perturbations internes din soient de type constantes. Daprs lquation (II.31), la fonction de transfert (Y/din)(s) sera nulle en rgime permanent si nous imposons S(0)=0, ce qui implique que S(s) doit se fractionner en :
S(s) = s S(s)
(II.33)
Pour garantir un gain statique unitaire de la fonction de transfert en poursuite (Y/Yc)(s), lquation (II.31) montre quil suffit maintenant de choisir un polynme T(s) qui vrifie la contrainte T(0)=R(0). A cette condition prs, le polynme T(s) est ici totalement libre (sous rserve que (T/S)(s) soit propre). Pour simplifier le transfert en poursuite, nous choisirons, dans un premier temps, T(s) multiple quelconque de D(s). III.2.1 Rsolution de lquatio n de Bezout Il convient dabord danalyser les degrs des polynmes du rgulateur RST. Le systme dorigine (B/A)(s) et le rgulateur (R/S)(s) tant propres, le degr de D(s) est dtermin par celui du produit (A(s). S(s)), ce qui implique donc :
deg(D)=deg(A)+deg(S)
(II.34)
Lgalit des coefficients des puissances de (s) dans chacun des deux membres de lquation (II.32) se traduira par un systme dquations linaires, en nombre gal (deg(D) +1). Pour que ce systme ait une solution, quels que soient les coefficients du polynme D(s), il faut que le rgulateur prsente au moins autant de degrs de libert que dquations. Or, nous disposons de :
deg(R)+1 coefficients pour R(s), deg(S) coefficients pour S(s), car S(0)=0.
deg(R) deg(A)
47
(II.37)
Alors, lidentit selon les puissances dcroissantes de (s), se traduit par le systme linaire suivant, dit systme de Sylvester :
1 a 1 0 1 M O O O O M an O 1 a1 0 O O O M O O M O an 0
0 s0 d0 s d 0 O 0 1 1 M M b1 O O M O O O M M M O O 0 sn d n = bn O 0 r0 d n + 1 O b1 r1 M O M M M 0 O M M M bn rn d 2 n + 1
(II.38)
48
S(z) = 0 pour z = 1
Le systme de Sylvester est alors videmment diffrent, mais les contraintes sur les degrs de A(z), B(z), R(z), S(z) et T(z) restent inchanges. Le nombre de ples placer pour un systme discret de degr (n) reste toujours au moins gal (2n). III.2.2 Stratgie de placement de ples robuste Pour raliser un compromis entre les performances et la robustesse, le polynme D(s) va tre factoris en D(s) = P(s).C(s) o P(s) est de degr n. Cela permet de rpartir les modes en boucle ferme en ples de filtrage (zros de C(s)) et ples de commande (poursuite) (zros de P(s)). Dans lobjectif de raliser cette factorisation, deux paramtres de synthse sont employs : un horizon de commande : Tc un horizon de filtrage : Tf Nous dduisons alors les ples de commande partir des ples en boucle ouverte, en effectuant dans lordre, comme le montre la Fig. II.14, les oprations suivantes :
Im (2)
d
Fig. II.14 : Construction des zros de P(s) partir de ceux de A(s)
(1) Les zros instables de A(s) sont remplacs par leurs symtriques par rapport laxe imaginaire, (2) Les zros complexes peu amortis sont ramens lamortissement spcifi,
49
Chapitre II : Filtre actif parallle, stratgie de commande, rgulation et dimensionnement ___________________________________________________________________________ (3) Les zros trop lents ( droite de la verticale dabscisse 1/Tc) sont ramens sur cette verticale, (4) Les zros restants demeurent inchangs.
Nous reconstituons alors le polynme P(s) partir des nouveaux zros obtenus. Les zros du polynme C(s) se dduisent partir du paramtre de synthse spcifi Tf, selon la mme procdure que pour P(s). Pour un rgulateur strictement propre, un ple nul (s=0) sajoute aux zros de A(s) et les n+1 zros de C(s) se dduisent dans ce cas de ceux de s.A(s) suivant toujours la mme procdure. Langle du secteur, dfini par les droites d et d, dans la Fig. II.14 fixe lamortissement des ples placs. Dans ce travail, un secteur de (90) (245) est choisi pour fixer lamortissement des ples placs 0,707. Nous pouvons associer Tc et Tf un certain nombre de rgles simples que nous allons numrer : Si le temps de rponse en asservissement est trop grand, diminuer Tc, Pour augmenter la robustesse et la marge de stabilit, augmenter Tc ou Tf ou les deux, Nous avons souvent intrt dstabiliser les horizons Tc et Tf pour accrotre la robustesse, Si la commande est trop sensible aux bruits des capteurs, augmenter Tf, Si la commande est trs sollicite en mode asservissement, augmenter Tc, ou prvoir un pr-filtrage adoucisseur de la rfrence. Remarque - Dans le cas discret, les ples de A(z) seront dabord convertis en ples en (s), puis modifis selon les horizons Tc et Tf, et reconvertis enfin en ples en z, pour dterminer P(z) et C(z).
50
un premier polynme correspondant l'action de rgulation : Rr(s) et Sr(s), un second polynme correspondant l'action de rejet des perturbations : Rf(s) et Sf(s) (ces deux polynmes agissant comme des filtres aux perturbations).
Cette approche de partitionnement va entraner la modification de l'quation de Bezout (II.32), qui va s'exprimer prsent par la relation suivante :
D(s) = A(s) S r (s). S f (s) + B(s) Rr (s). R f (s) 14243 4 4 14243 4 4 S(s) R(s)
(II.39)
Afin de filtrer les perturbations constantes, il suffit dliminer le gain statique du deuxime terme de la relation (II.31) : Vd(s)/din(s), en ajoutant un zro au polynme Sf(s) tel que S(0) = 0, ce qui est repris dans la relation suivante :
Sf (s) = s
(II.40)
En considrant comme prpondrantes, les composantes perturbatrices externes de frquence 100 Hz et 300 Hz, le rgulateur va devoir filtrer ces dernires pour ne pas les retrouver dans la commande u(s) (Cf. Fig. II.13). A partir de la relation (II.39) et afin de satisfaire lensemble de ces nouvelles contraintes, des zros cals sur les frquences (f1=w1/2p) de 100 Hz et ( f2=w2/2p) de 300 Hz sont rajouts au polynme Rf(s). La largeur et la profondeur des puits dattnuation peuvent tre imposes au concepteur en introduisant des paramtres de rglage r1 et r2. Ces deux paramtres 51
(II.41)
Les perturbations de tension de frquences suprieures f2 ntant pas prpondrantes, elles pourront tre filtres par une pente de 20 dB/dcade du rgulateur adapt, comme nous le verrons plus loin. III.3.1.2 Placement de ples du r gulateur RST Nous reprsentons sur la Fig. II.15 le schma de ralisation du rgulateur RST que nous avons appliqu dans la suite de notre tude.
Commande u
Vd-ref
T(s)/R(s)
+-
R(s)/S(s)
- 3Vs
Vd
Fig. II.15 : Boucle de rgulation structure RST du systme base de PLL Le systme dorigine B(s)/A(s) du systme base de PLL est modlis par la relation suivante :
B 3Vs (s)=s A
(II.42)
En nous basant sur la relation (II.39) , nous en dduisons le degr des polynmes R(s) et S(s) donn par les relations suivantes :
deg(R)=deg(Rr)+deg(Rf) deg(S)=deg(Sr)+deg(Sf)
Daprs les donnes prcdentes, nous savons que le polynme Rr(s) est dordre 1, ce qui correspond lordre du systme dorigine du systme base de PLL prsent en (II.42), et que le polynme Rf(s) est dordre 4 (relation II.41). Par consquent, le degr global du polynme R(s) est gal 5, tel que :
deg(R)=1 + 4 = 5
De la mme faon, en choisissant un rgulateur spcifique strictement propre, nous montrerons, en nous basant sur la relation (II.36), que le degr global du polynme S(s) est gal 6 :
deg(S)=1 + 5 = 6
Finalement, nous en dduisons le degr global du polynme D(s) qui est alors gal 7 :
deg(D) = 7
52
Fig. II.16 : Diagramme de Bode du rgulateur spcifique RST De lanalyse de ces rsultats, nous pouvons galement dduire que le filtrage des frquences suprieures 300 Hz est ralis par une pente de (20) dB/dcade, ce qui a t permis grce au choix dun rgulateur strictement propre. Cas de simulation - Le cas dun rseau dsquilibr a t choisi afin dexaminer la robustesse en termes de rapidit et de prcision du systme base de PLL lors dun changement de paramtre du systme dorigine et de variations des perturbations. La Fig. II.17 montre les rsultats de simulation de la tension dsquilibre du rseau lectrique (Vs123), de lamplitude de la composante directe de tension (Vd), de langle de la composante directe de tension (dd) reprsent sous la forme des sinus triphass, et de la composante directe triphase de tension (Vd123). Le dsquilibre de la tension est caractris par les amplitudes de tension suivantes : Vs1=0,65Vs, Vs2=0,8Vs, Vs3=0,95Vs avec (Vs=230 V). Le dsquilibre nest appliqu que sur quatre priodes et ne dbute qu linstant t = 0,12 s.
53
(V)
Vs123
(V)
Vd
Sin wd (V)
Vd123 (V)
t (s)
Fig. II.17 : Fonctionnement du systme base de PLL sous une tension dsquilibre Des rsultats de simulation de la Fig. II.17, nous observons que le systme base de PLL, que nous avons dvelopp dans cette partie, est robuste face aux variations de la tension Vs (systme dorigine) et aux perturbations. La rapidit et la prcision du systme base de PLL peuvent tre observes par un verrouillage instantane de langle dd, une frquence inchange de 50 Hz et par une amplitude exacte de Vd. On observe un petit retard de lextraction de Vd, ce dernier tant caus par le filtre passe-bas du quatrime ordre (Cf. Fig. II.8). En effet, langle dd est obtenu travers une boucle de rgulation alors que lamplitude Vd est lissu dun simple filtrage du quatrime ordre. Le retard caus par ce filtre peut tre minimis encore en trouvant un compromis plus judicieux entre la rapidit et lefficacit du filtre passe-bas.
54
B1(s) /A(s)
Fig. II.18 : Schma de la rgulation en poursuite par RST du courant du F.A.P De la Fig. II.18, la relation suivante, caractrisant la poursuite de la consigne, peut tre crite :
I inj (s) = (T. B1 )(s) I ref (s) (g.S. A + B1 .R)(s)
(II.43)
avec 1/g(s)=1/(t.s+1) le modle reprsentant le retard de calcul. La frquence de coupure de la boucle de rgulation doit tre choisie de manire pouvoir satisfaire les deux critres suivants :
dune part, la frquence de coupure doit tre plus grande que lharmonique le plus lev du signal de rfrence, de manire ne pas influencer la compensation des harmoniques,
dautre part, la frquence de coupure doit tre suffisamment loigne de la frquence de commutation pour que lon puisse modliser londuleur par un gain simple, comme lexprime la relation (II.30). III.3.2.2 Placement de ples du r gulateur RST A partir de la Fig. II.18 prcdente, le systme dorigine dcrit par lexpression (II.5), est donn par la nouvelle relation suivante :
Btot ( s ) 1 B (s) = 1 Atot ( s ) t .s + 1 A( s )
(II.44)
55
deg(D)=2n=8
Le degr des polynmes P(s) et C(s) est de 4 :
Les ples du systme dorigine vont subir les changements dcrits sur la Fig. II.14. Puis, ils sont projets sur labscisse (-1/Tc) pour former le polynme P(s), et sur labscisse (-1/Tf) pour le polynme C(s). Aprs avoir form le polynme de stabilit arbitraire (D(s)=P(s).C(s)), les polynmes R(s), S(s) et par consquent T(s) se dduisent en appliquant la relation (II.32). Le degr des polynmes R(s) et S(s) est le mme et est gal 4 pour un rgulateur propre, comme le donne la relation (II.35) :
deg(R)=deg(S)=n=4
T(s) est choisi comme un simple gain tel que (Iinj(s)/Iref(s)) = 1 quand la variable symbolique (s) tend vers 0.
III.3.2.3 Analyse des Performance s Le signal de rfrence poursuivre, pouvant contenir des courants harmoniques allant jusquau rang 40, la mthode de rgulation par RST telle que prsente ne peut pas satisfaire les performances requises par les normes internationales. Ces mauvais rsultats sexpliquent en grande partie par linjection, aprs rgulation, dun dphasage entre les courants de rfrence (Iref) et ceux injects en sortie de la boucle de rgulation (Iinj). De plus, nous avons mesur que ce dphasage augmentait :
dune part avec lordre de la fonction de transfert en boucle ferme, dautre part avec la frquence des signaux de rfrence poursuivre.
Dans le but de rsoudre ce problme de dphasage, nous allons dabord identifier leffet de ces dphasages, puis nous prsenterons une mthode de rgulation amliore base sur un rgulateur RST, que lon nommera RSTam.
56
Xref (b)
Xinj
(c)
Xrel
Xidal
t (s)
Fig. II.19 : Leffet du dphasage sur la qualit de compensation On voit bien, sur les dernires courbes de la Fig. II.19, quel point le problme de dphasage peut dgrader la qualit de compensation du filtre actif parallle. Ce problme n'est pas seulement inhrent au rgulateur RST. On pourrait observer, en effet, les mmes effets ngatifs en utilisant n'importe quelle autre mthode de rgulation linaire connue, comme les rgulateurs PI, PID, H2, H,etc. L'application d'un rgulateur de ce type n'est donc pas une solution allant de soi, il nous a donc paru opportun de continuer pousser notre tude sur les rgulateurs RST. III.4.2 Approche de conceptio n du rgulateur amlior RTSam Les principes du rgulateur RST sont gnralement utiliss quand la rfrence suivre est constitue de signaux constants ou une frquence seule et relativement faible. Dans ce dernier cas, le dphasage entre les rfrences identifies et la sortie de la boucle ferme est acceptable. Par contre, si la rfrence suivre est compose de signaux plusieurs frquences, le dphasage nest plus ngligeable. En effet, le dphasage augmente avec la frquence, comme le montre le diagramme de phase donn par la Fig. II.20 et le Tab. II.2 (voir rsultats de synthse plus loin dans le paragraphe III.5.1).
57
(II.45)
(II.46)
avec D(s)=(Atot.S+Btot.R)(s) le polynme de stabilit arbitraire. Afin de minimiser lerreur e(s) et de la faire tendre vers zro, le numrateur de la relation (II.46) doit tre nul. Puisque les polynmes R(s) et S(s) ont dj t choisis et fixs, cest le polynme T(s) qui devra remplir cette condition. De plus, lordre du polynme T(s) devra tre choisi de faon ce que la fonction de transfert T(s)/R(s), prsente sur la Fig. II.18, soit ralisable ou causale. La stabilit du systme, dfinie par le rgulateur RST classique, nest pas modifie si lon utilise un polynme T(s) stable [DELA96]. Dans le cas tudi, le quatrime ordre ayant t choisi pour les polynmes R(s) et S(s), le polynme T(s) est du quatrime ordre, ce qui est le plus grand ordre utilisable pour T(s) tout en satisfaisant le critre de causalit de T(s)/R(s). Ainsi, le polynme T(s) scrira partir de lexpression suivante :
T ( s ) = t1 s 4 + t 2 s 3 + t3 s 2 + t4 s + t5
(II.47)
En reportant la relation (II.47) dans (II.46) et en remplaant s=jw, le numrateur de la relation (II.46) va pouvoir scrire sous la forme suivante:
(T .Btot - D)(jw ) = Re ( w ) + jI m ( w )
(II.48)
Afin de faire tendre la relation (II.48) vers zro, les deux parties relle Re(w) et imaginaire Im(w) doivent tre gales zro. Cette contrainte permet davoir deux quations pour chaque frquence de chaque zro plac. En prenant en compte le nombre des paramtres inconnus (t1,t2,t3,t4,t5) de la relation (II.47), uniquement trois frquences peuvent tre choisies. Les zros sont placs proches des ples lents de la boucle ferme. Le gain et la phase de la fonction de transfert en boucle ferme Iinj(s)/Iref(s) sont donns par la Fig. II.20 et le Tab. II.2.
58
Fig. II.20 : Fonction de transfert (Iinj/Iref) avec les rgulateur RST et RSTam
50 1 -1
250 350 550 650 850 950 1150 1250 1450 1550 1 -8 1 -11 1 -1 1 -18 1,1 -4 0.99 0.98 0.98 0.97 -21 1 -6 -27 -30 -36 0,9 -10 1 0,98 -10 -11 0,96 -40 0,9 -10 0,95 -46 0,9 -7 0,94 -49 0,9 -5
1 1 -0,01 -0,5
Tab. II.2 : Gain et phase du systme en boucle ferme utilisant RST et RSTam
59
60
(II.49)
avec Vs la tension du rseau au point de raccordement et In-ch le courant de la charge non linaire (Cf. Annexe-I, A-I.1). La puissance apparente du filtre actif (Sf) compensant le courant harmonique Ih, inject par un pont redresseur triphas thyristors (pont de Graetz), est donne par lquation suivante :
2 S f = Dch = [Link] .I h
(II.50)
Ce courant harmonique Ih qui doit tre cr par le filtre actif peut scrire de la faon suivante :
2 I h = I n - ch - I 2 fon
(II.51)
avec Ifon le courant fondamental consomm par la charge non linaire. Le courant fondamental et le courant de la charge peuvent tre crits en fonction du courant direct de la charge non linaire, Id, de la faon suivante :
I n - ch =
2 6 I d , I fon = Id 3 p
(II.52)
En reportant les relations (II.51) et (II.52) dans celles de (II.49) et (II.50), on obtient le rapport des puissances (th) donn par lexpression suivante :
th = Sf Sch = 0.24 I d -a ( 2 )I d 3
(II.53)
61
(II.54)
avec Ud la tension de la charge non linaire ct continu. Des expressions (II.53) et (II.54), on obtient lexpression finale de th donne par la relation suivante : p2 -9 .cos a 0.3 cos a th = p (II.55)
La Fig. II.21 suivante montre la variation du rapport des puissances, (th), du filtre actif parallle par rapport celle de la charge non linaire, en fonction de langle dallumage des thyristors (a).
th
a (deg)
Fig. II.21 : Rapport des puissances du filtre actif parallle et de la charge non linaire pour la compensation des courants harmoniques A partir de cette figure, on remarque que pour a=0 la puissance maximale du filtre actif est de Sf30% Sch. Cette puissance diminue avec laugmentation de langle dallumage (a) grce la diminution du courant harmonique.
62
(II.56)
t hr =
(II.57)
En reportant les relations (II.51) et (II.52) dans celles de (II.54) et (II.57), on obtient lexpression du rapport des puissances (thr) suivante :
t hr = cos a 1 -
9
p
2
cos a
(II.58)
La Fig. II.22 donne la reprsentation graphique du rapport des puissances (thr) en fonction de langle dallumage des thyristors a de la charge non linaire.
thr
a (deg)
Fig. II.22 : Rapport des puissances du filtre actif parallle et de la charge non linaire pour la compensation du courant harmonique et de la puissance ractive A partir de la Fig. II.22, on montre que pour un angle damorage a=0, on trouve le mme rapport de puissance que celui obtenu dans le cas de compensation prcdent (Sf30% Sch). Par contre, dans ce cas de compensation, le filtre actif atteint ici un rapport maximum de puissance de Sf52% Sch pour un angle de a42. Les rsultats graphiques de la Fig. II.22 montrent galement que le filtre actif parallle, et du point de vue conomique, peut tre appliqu pour compenser la fois le courant harmonique et
63
ractive
Dans cette troisime tude, nous nous intresserons au calcul du rapport des puissances apparentes dans le cas dune compensation des courants harmoniques et dsquilibrs et de la puissance ractive. Les courants harmoniques et la puissance ractive sont consomms par un pont de Graetz tandis que le courant dsquilibr est caus par une charge linaire dsquilibre connecte en parallle avec celle non linaire. Le rseau tudi tant de trois fils, le dsquilibre de courant est reprsent uniquement par la composante inverse du courant Ii de la charge linaire. Le nouveau rapport des puissances (thri) calcul dans ce cas peut scrire de la faon suivante :
t hri = Sf Sch = ( [Link] .I h )2 + ( [Link] .I fon . sin a )2 + ( [Link] .I i )2 [Link] .I n - ch
(II.59)
(II.60)
La Fig. II.23 donne la reprsentation graphique du rapport des puissances (thri) en fonction de langle dallumage des thyristors a de la charge non linaire. Ce rapport de puissance est donn pour plusieurs valeurs du taux inverse de courant (X=Ii/Ifon), Ifon tant toujours le courant fondamental de la charge non linaire. A partir de la Fig. II.23, on remarque que pour X=0, on retrouve la mme courbe que dans le cas de compensation prcdent. Les rsultats graphiques de la Fig. II.23 montrent galement que la puissance du filtre actif parallle augmente de faon quasi linaire avec laugmentation du taux inverse du courant. Du point de vue conomique, le filtre actif parallle peut ainsi tre install pour compenser un faible taux du courant inverse (X<0,1). Il peut aussi tre appliqu pour compenser uniquement le courant harmonique tout en ayant la possibilit de compenser, dans la limite de la puissance restante, une partie du dsquilibre dans un cas de sous-charge.
64
thri
a (deg)
Fig. II.23 : Rapport de puissances pour la compensation des courants harmonique et dsquilibrs et de la puissance ractive
65
Rseau Electrique
Is1 Is2
Rc Lc
es
Rs
Ls
Filtre de Sortie
Lf2,Rf2 Lf1,Rf1 Cf Rf
Cdc
Vs1,2,3
PLL
Is1,2,3 Iinj1,2,3 Vdc
Fig. II.24 : Structure gnrale du filtre actif parallle au sein du systme tudier La modlisation globale de la structure, prsente sur la Fig. II.24, a permis dtudier la validit fonctionnelle du filtre actif parallle travers plusieurs simulations. Les valeurs des lments caractrisant cette structure ont t regroupes dans le Tab. II.3.
66
Tab. II.3 : Valeurs des lments construisant la structure tudie Finalement, ces simulations nous permettront dvaluer la prcision et la rapidit du systme base de PLL, et de la mthode de rgulation (RSTam) proposes dans ce chapitre. Ainsi, nous avons labor plusieurs cas de simulation que nous regrouperons dans trois cas dtude :
le cas de la compensation gnrale des harmoniques et dsquilibre de courant, et de la puissance ractive, le cas de la compensation dans un milieu perturb, le cas de la compensation pour une structure classique de filtre actif parallle (filtre de sortie du premier ordre avec un rgulateur PI).
Ces simulations ont t ralises partir du code numrique SABER et nous ont permis danalyser les performances du filtre aussi bien dans le domaine temporel que spectral.
67
In-ch1 (A)
IL1 (A)
t (s)
(A)
IL1
(A)
f (Hz)
Fig. II.26 : Analyse spectrale de la compensation du courant harmonique V.1.2 Compensation des cou rants harmoniques et de la puissance ractive Dans ce cas, on tudie la compensation du courant harmonique et de la puissance ractive consomms par la mme charge non linaire avec une charge linaire toujours dconnecte. La rapidit de compensation du filtre actif est tudie lors dun changement du courant de la charge non linaire caus par une variation de langle dallumage entre a=5 et a=42. La Fig. II.27 prsente lanalyse temporelle de la phase 1 du courant de la charge (In-ch1), des courants identifi (Iref1) et inject (Iinj1) superposs, du courant du rseau aprs filtrage (IL1) et de la tension continue (Vdc). On remarque que la boucle de rgulation de la tension continue (Cf. Fig. II.11) arrive poursuivre sa rfrence de (Vdc-ref = 840 V) dans un temps de rponse de 20 ms avec un petit dpassement de 3% et un taux dondulation de 0,4%. Un courant sinusodal du 68
In-ch1 (A)
(A)
Iref1 Iinj1
IL1 (A)
t (s)
(A)
f (Hz)
69
ractive
Le filtre actif permet ici de compenser le courant harmonique ainsi que la puissance ractive, causs par la charge non linaire (a=42), de mme que les squences ngatives du courant de la charge linaire dsquilibre Rch1,2,3. Les rsultats de la Fig. II.29 montrent lanalyse temporelle du courant de la charge non linaire (In-ch123) et du courant de la charge dsquilibre (IL-ch123). Ils donnent galement les courants de rfrence des trois phases identifis (Iref-123) ainsi que ceux injects au rseau (Iinj-123) superposs et le courant circulant en amont du rseau (IL123). Un facteur de puissance unitaire tmoigne de la bonne compensation de la puissance ractive. Le courant de la ligne aprs compensation est un courant quilibr sinusodal quasiment dpourvu des harmoniques y compris ceux dus aux dcoupages.
(A) IN-ch1,2,3 IL-ch1,2,3 Iref1, Iinj1 Iref2, Iinj2 Iref3, Iinj3 IL1,2,3
t (s)
Fig. II.29 : Analyse temporelle de la compensation des courants harmonique et dsquilibr et de la puissance ractive
In-ch1 IL1
(A)
f (Hz)
70
(A)
(A)
(A)
(A)
(A)
(A)
f (Hz)
Fig. II.31 : Analyse spectrale montrant la compensation de la composante inverse du courant La Fig. II.30 montre le spectre harmonique de la phase 1. Le THD global du courant de la charge non linaire est de 27,08%, et de 3,88% ct rseau lectrique, ces deux THD tant reports au mme courant fondamental de la charge non linaire. La Fig. II.31 montre lanalyse spectrale des courants des charges perturbatrices (Is1,2,3) et du rseau lectrique aprs filtrage (IL1,2,3). Le taux inverse du courant de la charge dsquilibre est de 10,7% et de 0,85% ct rseau lectrique, ces deux taux inverses tant rapports au mme courant fondamental de la charge non linaire.
71
t (s) Fig. II.32 : Analyse temporelle de la PLL sous une tension dsquilibre du rseau lectrique
Vs1,2,3 (V)
In-ch1,2,3 (A)
IL1,2,3 (A)
t (s)
Fig. II.33 : Analyse temporelle de la compensation du courant harmonique sous une tension dsquilibre du rseau lectrique
72
In-ch1 IL1
(A)
f (Hz)
Fig. II.34 : Analyse spectrale de la compensation du courant harmonique sous une tension dsquilibre du rseau lectrique
aux dcoupages
Le premier cas simul est un filtre actif parallle sens compenser, en dconnectant la charge linaire, le courant harmonique engendr par le pont de Graetz thyristors de 100 kVA avec un angle dallumage de a=42. La Fig. II.35 prsente lanalyse temporelle du courant de la phase 1 de la charge (In-ch1), des courants identifi (Iref1) et inject (Iinj1) superposs et du courant du rseau aprs filtrage (IL1).
73
In-ch1 (A)
IL1 (A)
t (s)
Fig. II.35 : Analyse temporelle de la compensation du courant harmonique On remarque daprs cette figure que lallure du courant du rseau aprs filtrage prend la forme dune sinusode. Ce courant est masqu par les composantes dues la frquence de commutation que le filtre de sortie, du premier ordre, na pas pu empcher de se propager sur le rseau lectrique. La Fig. II.36 prsente lanalyse spectrale, phase 1, du courant de la charge non linaire et de celui du rseau aprs filtrage. Le THD du courant ct charge est de 27,1%, il est de 4,8% ct rseau sans compter les harmoniques dus la frquence de commutation. Le THD du courant du rseau, prenant en compte les harmoniques dus aux dcoupages, est de 13,9%. Un THD acceptable de 4,8% est suprieur 3,66% obtenu pour un cas similaire simul en paragraphe ( V.1.2) en employant le rgulateur RSTam. Cette diffrence est due au dphasage caus, en boucle ferme, par le rgulateur employ de type PI comme le montre le Tab. 4 :
F (Hz) G PI PH 50 1 -0,2 250 1 -0,6 350 1 -5 550 1 -8 650 0,98 -9 850 0,97 -12 950 0,97 -13 1150 0,96 -16 1250 0,95 -17 1450 0,94 -21 1550 0,93 -25
Tab. II.4 : Gain et phase en boucle ferme utilisant PI Avant danalyser ces rsultats, il est intressant de comparer lefficacit, en termes du filtrage des composantes haute frquence, des deux filtres de sorties tudis. La Fig. II.37 montre lanalyse spectrale comparant le contenu harmonique hautes frquences des courants IL1 obtenus dans ce paragraphe et dans celui tudi au ( V.1.2). La Fig. II.37 explique le THD lev (13,9%) obtenu pour le cas dun filtre de sortie du premier ordre. En effet, ce filtre nayant que la mme valeur des inductances employes dans le cas dun filtre de sortie en T, ne peut empcher correctement la propagation des composantes dues aux dcoupages sur le rseau lectrique.
74
In-ch1
(A)
IL1
(A)
f (Hz)
In-ch1
(1er ordre)
(A)
f (Hz)
Fig. II.37 : Analyse spectrale des effets des filtres de sortie du troisime du premire ordre V.3.2 Compensation des cou rants harmoniques avec rejet des composantes dues
aux dcoupages
Afin dobtenir le mme taux de rejet des composantes hautes frquences obtenu par le filtre en T, une deuxime simulation a t effectue. Dans ce cas, un filtre de sortie du premire ordre dont linductance est de Lf=1000 mH doit tre employ. Le THD du courant de rseau aprs filtrage est de 12,8% comme le montre la Fig. II.38. Le taux lev de THD obtenu dans ce cas peut tre expliqu par une dynamique trs lente du filtre actif parallle cause par une valeur leve de linductance Lf. En effet, une valeur leve de Lf ne satisfait pas lquation II.4 de commandabilit de londuleur de tension. Cette analyse prouve quil est difficile au filtre de sortie du premire ordre de faire le bon compromis entre un bon rejet des composantes hautes frquences et une bonne dynamique du filtre actif parallle.
75
In-ch1
(A)
IL1
(A)
f (Hz)
Fig. II.38 : Analyse spectrale de la compensation du courant harmonique (cas dune valeur leve dinductance)
76
VI. CONCLUSION
Dans ce chapitre, nous avons prsent le filtre actif parallle comme une solution de dpollution des rseaux lectriques des courants perturbateurs comme le courant harmonique, dsquilibr et ractif. La mthode didentification des puissances instantanes a t choisie pour offrir la possibilit didentifier un, plusieurs ou tous les courants perturbateurs du rseau lectrique. Cette mthode ntant pas adapte au cas dune tension perturbe du rseau, un systme base de PLL a t propos pour dlivrer la composante directe de la tension du rseau lectrique. Le systme base de PLL, bas sur lemploi dun rgulateur RST spcifique, dlivre la composante directe de tension avec rapidit, prcision et robustesse. Un filtre de sortie du troisime ordre, appel filtre en T, a t propos afin dempcher les composantes dues aux dcoupages de se propager sur le rseau lectrique. Ce type de filtre de sortie rsout les problmes secondaires causs par lemploi dun filtre passif auxiliaire associ un filtre de sortie du premier ordre. Afin de raliser les objectifs de rgulation (robustesse en stabilit et en performance), un rgulateur RST du courant du filtre actif parallle a t tudi. Le problme de dphasage, entre les courants identifis et ceux injects par le filtre actif, occasionn par le RST et par tous les rgulateurs linaires, a t introduit. Ce problme de dphasage dgrade la qualit de compensation du filtre actif parallle limitant son exploitation, alors des progrs raliss dans le domaine de lautomatique permettent de le surmonter. Afin de rsoudre le problme du dphasage tout en gardant les avantages de robustesse raliss par le rgulateur RST, un rgulateur amlior, appel RSTam, a t propos. Ce rgulateur se base sur le mme principe de placement de ples que le RST classique, plaant les zros afin de minimiser lerreur entre la consigne et la sortie de la boucle de rgulation et de minimiser par consquent le dphasage. Une comparaison, en termes de suivi de consigne, entre les rgulateurs RST et RSTam a t effectue afin de valider thoriquement, dans un premier temps, la nouvelle mthode de rgulation propose. Afin de valider la fonctionnalit de lensemble des mthodes proposes (systme base de PLL, filtre de sortie en T et le rgulateur RSTam), plusieurs simulations ralises sur SABER ont t effectues. La rapidit et la qualit de compensation des courants harmoniques, des courants dsquilibrs et de la puissance ractive ont t observes. Enfin, la difficult dassurer une bonne qualit de compensation, partir dun filtre de sortie du premier ordre avec un rgulateur PI du courant, a t montre travers deux simulations. Finalement, ltude du filtre actif parallle, que nous avons aborde dans ce chapitre, nous a permis datteindre de trs bonnes performances de compensation des perturbations en courant, et ceci dans les conditions extrmes de fonctionnement que nous nous tions fixes. Une partie de nos objectifs de recherche a donc t atteint. Cependant, nous devons prsent complter ce travail par une tude spcifique sur la compensation des perturbations en tension dans les mmes conditions dapplication que prcdemment.
77
78
79
80
I. INTRODUCTION
Lobjectif de recherche de ce chapitre est de contribuer lamlioration des performances de compensation des perturbations en tension partir dune structure de filtre actif srie. En vue de dpolluer les rseaux lectriques de toutes sortes de perturbation en courant et en tension, une structure combine parallle-srie (UPQC) sera aussi aborde et dveloppe dans ce chapitre. Finalement, ce chapitre sera dcompos en quatre parties : Dans une premire partie, nous tudierons plusieurs structures de filtre actif srie. La description complte de la structure la plus utilise sera prsente. Puis, nous effectuerons un classement des mthodes didentification des tensions perturbatrices dans lobjectif de montrer les avantages et les limites dapplication des unes et des autres. Cela nous conduira proposer une nouvelle mthode didentification base sur un systme base de PLL dj intgr dans la mthode didentification des courants perturbateurs. Cette mthode aura pour objectif de rsoudre les problmes en terme dadaptabilit et de calcul dj rencontrs avec dautres mthodes. La prcision et la robustesse de cette nouvelle mthode didentification sera analyse travers plusieurs simulations effectues sous SABER. Dans une deuxime partie, la boucle de rgulation de tension du filtre actif srie sera tudie dans le but de compenser les tensions perturbatrices de tout rang. Les avantages et les limites des rgulateurs PID et RST, appliqus dans ce cas, seront analyss. De plus, les principes du rgulateur amlior RSTam, dj utilis pour la rgulation du courant du filtre actif parallle, seront galement proposs. Une synthse comparative entre ces trois rgulateurs sera ensuite effectue en terme de poursuite en amplitude et en phase dans le but de dsigner le rgulateur permettant la compensation la plus large possible. Dans une troisime partie, nous effectuerons une tude dvaluation de la puissance apparente du filtre actif srie lors de la compensation de plusieurs cas de creux de tension. Loptimisation de la puissance du filtre actif srie sera ensuite tudie en comparant les mthodes didentification avec et sans le systme base de PLL. Dans une quatrime partie, des simulations comparatives du filtre actif srie, employant les trois mthodes de rgulation (PID, RST, RSTam), seront effectues afin danalyser les performances de ce dernier en termes de robustesse et de qualit de compensation. Finalement, la dernire partie de ce chapitre sera consacre au dveloppement des performances de lUPQC install pour compenser tout type de perturbation pouvant apparatre dans un rseau lectrique. Afin datteindre ces objectifs, lUPQC va profiter de tous les dveloppements dj raliss lors de ltude des filtres actifs parallle et srie. Des rsultats de simulation, effectus sous SABER, valueront la qualit de compensation du filtre actif combin.
Partie Puissance
Transformateurs dinjection avec protection By-Pass Filtre de sortie Onduleur Elments de stockage Alimentation DC
Tension du rseau
Partie ContrleCommande
Fig. III.1 : Structure gnrale du filtre actif srie Afin de mieux comprendre cette structure, les parties puissance et contrle-commande vont tre dtailles.
la source de tension aux bornes de londuleur ne doit jamais tre court-circuite, les deux interrupteurs dun mme bras doivent donc avoir deux commandes complmentaires, la source de courant ct alternatif ne doit jamais tre en circuit ouvert, le courant circulant doit toujours trouver un chemin libre do la mise en antiparallle de diodes avec les interrupteurs.
Londuleur de tension employ se compose de six interrupteurs repartis, dune faon symtrique, sur trois bras interrupteurs commands la fermeture et l'ouverture (IGBT, GTO) avec des diodes en antiparallle, comme le montre la Fig. III.2.
82
Cdc1 Rsf1 Rsf2 Rsf3 Lsf1 Lsf2 Lsf3 Vinj1 Csf1 Csf2 Csf3 Vinj3 - Vdc/2 Vinj2
Cdc2
Fig. III.2 : Partie puissance du filtre actif srie II.1.2 Systme de stockage d nergie Le choix dun systme de stockage deux condensateurs avec un point milieu (Cdc1,Cdc2), comme le montre la Fig. III.2, donne la possibilit de mesurer, aux bornes des condensateurs, une tension continue Vdc sans fluctuation. Cela permet aussi davoir trois phases indpendantes du filtre actif srie fonctionnant comme un demi-pont. Par contre, ce choix dans le cas du filtre actif parallle aurait dgrad la dynamique de ce dernier, lequel dpendait, entre autre, de la tension continue Vdc. En effet, un systme de stockage sans point milieu permet londuleur davoir une tension de sortie dont la valeur maximale est de Vf=2/[Link] (Cf. Tab. II.1), au lieu dune valeur maximale de Vdc/2 dans le cas dune structure avec un point milieu. Afin dassurer que londuleur dispose dune tension maximale sa sortie, il faut que les deux tensions continues ( Vdc/2) soient maintenues dans une plage bien spcifie. En effet, la valeur minimale des tensions continues dtermine directement lamplitude de la tension de sortie que le filtre actif srie est capable de gnrer. Pour cela plusieurs solutions sont envisageables : une alimentation DC base dun pont redresseur diodes, triphas [CA92] ou monophas [DI97], inclus dans la structure, comme la montr la Fig. III.1 : le pont redresseur tant aliment par le rseau lectrique, chaque condensateur aura la valeur crte de la tension simple du rseau. Ce montage prvoit des valeurs leves des capacits de stockage pour assurer, dans le cas dun rseau perturb, un faible taux dondulation aux bornes des condensateurs [TO99]. Il faut galement sassurer, dans ce cas, du dimensionnement en courant des diodes du pont redresseur, un systme DC indpendant base des batteries alimentant les condensateurs de stockage dnergie [VI2000]. Cette solution assure au filtre actif srie une dure dautonomie assez large avec une tension de sortie relativement stable pendant la dcharge, pour les cas des filtres actifs de puissance leve et de tension continue relativement basse, lutilisation des super-condensateurs ou des supraconducteurs peut tre envisage [VI2000].
Dans le premier cas propos, la chute de tension du condensateur pendant la priode de dcharge ne doit pas influencer la qualit de compensation du filtre actif srie. Cela est possible si la diffrence entre la tension nominale (Vdc-n) et celle la fin de la priode de dcharge (Vdc-f) est assez petite.
83
i .P x Dt 3 100
(III.1)
avec P, la puissance active de la charge protger et i le nombre des phases perturbes. La variation dnergie dans un condensateur pendant la compensation du creux de tension est gale :
1 Cdc (V 2 - n - V 2 - f ) dc dc 2
(III.2)
La chute de tension aux bornes de chaque condensateur pendant la priode de dcharge, nous permet de dimensionner la valeur de la capacit. Donc, pour pouvoir compenser un creux de tension de profondeur x%, la chute de tension aux bornes de chaque condensateur doit rester infrieure (1-x)%. Cela revient dire quil faut assurer que la tension continue aux bornes du condensateur soit suprieure la tension de rfrence (tensions perturbatrices identifies). II.1.3 Filtre de sortie Le filtre de sortie dun filtre actif srie est gnralement un filtre passif du deuxime ordre (Lsf, Rsf, Csf), comme le montre la Fig. III.3. Lensemble onduleur et filtre de sortie se comportant comme une source de tension, le filtre de sortie sert dune part attnuer les composantes dues aux commutations de londuleur, dautre part connecter le filtre actif au rseau lectrique.
Vsrie Arrive Rseau Lt, Rt IL Vinj Csf Lsf, Rsf Vsf Dpart Charge
Onduleur
Fig. III.3 : Schma de modlisation du filtre de sortie A partir de la Fig. III.3, avec Vsf la tension de sortie de londuleur, Vinj la tension aux bornes du condensateur Csf, et IL le courant de la charge protger, les quations qui modlisent le filtre de sortie sont :
Vinj (s) = B1 (s) B (s) Vsf (s) + 2 I L (s) A(s) A(s)
(III.3)
84
(III.4)
et :
(III.5)
(B1/A)(s) reprsente la fonction de transfert du filtre de sortie correspondant au systme dorigine, et (B2/A)(s) la fonction de transfert correspondant au modle de perturbation. Ces perturbations sont causes par le passage du courant de la charge protger IL.
La frquence de rsonance fcs du systme dorigine est donne par la relation (III.6) :
f cs = 1 2.p . Lsf .Csf
(III.6)
Le filtre de sortie du deuxime ordre est dimensionn pour rejeter les composantes dues la frquence de commutation de londuleur qui vaut 12 kHz dans ce travail. Pour cela un rejet de 55 dB, permettant de ngliger les composantes hautes frquences passant du ct rseau lectrique, est obtenu pour une frquence de 500 Hz, comme le montre la Fig. III.4.
Fig. III.4 : Diagramme reprsentant le gain du filtre de sortie La rsistance damortissement Rsf est choisie pour que le facteur de qualit Q, donn par lquation (III.7), soit Q=50. 85
(III.7)
II.1.4 Transformateurs dinj ection et protection du F.A.S Les transformateurs dinjection de tension (Lt,Rt) permettent dinjecter en srie avec le rseau lectrique et avec le taux de transformation souhait la tension produite par londuleur, comme le montre la Fig. III.3. Par consquent, ces transformateurs, considrs comme des transformateurs de courant, permettent de rduire, par leurs taux de transformation, le courant de la charge traversant londuleur de tension. A cela sajoute que ces transformateurs, grce leurs caractristiques magntiques de saturation, contribuent la protection du filtre actif srie contre les dfauts ct charge [MOR99]. Les dsavantages de lemploi de ces transformateurs sont leurs cots levs et la chute de tension occasionne par leurs inductances de fuite (Lt,Rt). Le choix dun filtre actif srie sans transformateur dinjection ncessite lemploi de trois onduleurs de tension monophass avec trois sources continues indpendantes. Cette solution augmentant lnergie ncessaire pour compenser les mmes creux de tension, un filtre actif srie avec un onduleur multicellulaire et un seul bus continu a t propos [VI2000]. Dans ce travail, trois transformateurs dinjection de tension seront employs et dimensionns pour supporter le courant de la charge protger et la tension maximale injecte au rseau lectrique. Le taux de transformation sera unitaire. Le filtre srie est protg contre un court-circuit en aval ct charge par un systme by-pass, constitu par une impdance variante, par deux thyristors en antiparallle en srie avec une petite rsistance, et par des transformateurs circuit magntique saturable [MOR99].
Vref
Vinj
Rgulateur
Modulatrice
Fig. III.5 : Principe de commande de londuleur par MLI La tension de sortie de londuleur Vsf est gale :
Vsf = K .Ventre Vdc K = 2. Vp
(III.8)
86
Chapitre III : Filtre actif srie et combinaison parallle-srie actifs ___________________________________________________________________________ avec V p la valeur crte de la tension de la porteuse, et Vdc la tension ct continu de londuleur.
II.2.2 Identification des tensi ons perturbatrices La mthode didentification sert calculer les tensions perturbatrices qui sont injectes par londuleur, en opposition de phase, pour dpolluer la tension aux bornes de la charge protger. Ces tensions perturbatrices reprsentent les tensions dsquilibres et harmoniques, et les creux de tension. Plusieurs mthodes didentification ont dj t proposes dans la littrature. La plupart des mthodes sont bases sur le calcul des composantes symtriques dans le repre de Park (d,q), ncessitant une bonne connaissance du rseau lectrique en prsence dharmoniques de tension. Dautre part, la complexit de ces calculs rend plus difficile leur implantation. Cest pour ces raisons que nous proposons une nouvelle mthode didentification des tensions perturbatrices base sur lutilisation du systme base de PLL dj utilis dans lidentification des courants perturbateurs. II.2.3 Rgulation de la tensio n du filtre actif srie Aprs avoir identifi les tensions perturbatrices et afin dassurer une injection rapide, prcise et robuste de ces tensions, plusieurs mthodes de rgulation de la tension de sortie de londuleur ont t proposes dans la littrature. Dans ce chapitre, et afin de raliser les objectifs de rgulation dtaills prochainement, deux rgulateurs PID et RST ont t tudis. Les problmes ainsi que les limites de ces deux rgulateurs ont t analyss. De plus, nous avons galement appliqu notre mthode de rgulation RSTam, dj employe dans la rgulation du courant du filtre actif parallle.
87
[Vabc ] = [Vab
Vbc Vca ]
(III.9)
Ces tensions peuvent tre dcomposes en trois squences positives (Vp) (appeles aussi directes), ngatives (Vn) (appeles aussi inverses) et homopolaires (Vo), comme le dcrit lquation (III.10) :
1 Vab 1 1 V0 V = 1 1 a a 2 V bc p 3 Vn 1 a 2 a Vca
(III.10)
avec a=ej120.
Considrant que le rseau tudi est compos de trois fils, les composantes homopolaires Vo deviennent nulles. Afin de calculer la composante ngative, le premier paramtre de compensation doit tre dfini. Il sagit du facteur ngatif Un, appel aussi taux inverse, qui est le rapport entre les composantes ngative et positive de tension :
Un = Vn Vp
(III.11)
Afin de calculer la composante positive ncessaire pour rguler la tension de la charge protger, le deuxime paramtre de compensation doit tre dfini. Il sagit du facteur damplitude Mp qui est le rapport entre la composante positive et la tension de rfrence Vd-des :
Mp = Vp Vd -des
(III.12)
Les tensions composes ncessaires pour compenser les tensions ngatives sont :
[U c ] = -Vn [Tn ] [Tn ] = 1 a a 2
(III.13)
En reportant les valeurs de Vn des quations (III.11) et (III.12), dans lexpression (III.13), ce qui nous donne : 88
[U c ] = -U n M p [Tn ]
les tensions ngatives simples (Vn1,2,3) peuvent tre obtenues par lquation suivante :
(III.14)
[V
n1 ,2 ,3
]=
[T ] Uc = -U n M p n 1 -a 1 -a
(III.15)
Les tensions composes ncessaires pour rguler la tension de la charge afin dobtenir une amplitude Vd-des sobtiennent par :
[M c ] = -V p Tp Tp = 1 a 2 a
[ ]
[ ] [
(III.16)
[M c ] = (1 - M p ) [Tp ]
Les tensions positives simples Vp1,2,3 peuvent tre obtenues par lquation suivante :
(III.17)
[V
p 1 ,2 ,3
] = 1 Ma p
(1 - M ) [T ] 1-a
p 2 p
(III.18)
A partir des quations (III.15) et (III.18), on obtient les tensions ncessaires pour rguler et compenser le dsquilibre de tension de la charge protger : Vref 1,2 ,3 = [Vn1,2 ,3 ] + V p1,2 ,3 Tp T U [T ] V - M p p 2 + n n ref 1 ,2 ,3 = 2 1 -a 1-a 1-a
] [ ]
[ ]
(III.19)
Ces mthodes proposes pour calculer les composantes ngative et positive ncessaires pour rquilibrer et rguler la tension de la charge protger, ne prvoient pas la prsence dans la tension du rseau des composantes homopolaire et harmoniques. De plus, ces mthodes didentification ne conservant pas langle de la composante positive de tension, lors de la rgulation de tension de la charge, peuvent occasionner un surdimensionnement du filtre actif srie, comme cela est expliqu par la suite. II.3.1.2 Mthodes bases sur le ca lcul des perturbations dans le repre (d,q) Ces mthodes utilisent la dcomposition des tensions perturbatrices en squences positive (Vp), ngative (Vn), et homopolaire (Vo). En passant par la transformation de Park nous obtenons pour chaque composante symtrique des valeurs constantes (Vpd,q-mes,Vnd,q-mes,Vod,q-mes) [TO99, SAN 2000]. En fixant nos composantes dsires dans le repre direct et en quadrature aux grandeurs (Vpdp n n des,V q-des) pour la composante positive, aux grandeurs (V d-des,V q-des) pour la composante ngative, et aux grandeurs (Vod-des,Voq-des) pour la composante homopolaire, les tensions de rfrence dans le repre de Park sont dtermines de la faon suivante :
89
(III.20)
Les tensions perturbatrices (tensions de rfrence) triphases sont ensuite calcules grce la transformation suivante :
Vref 1 p n C23 P( w ,t ) V d - ref + C23 P( -w ,t ) V d - ref + V V n V p q - ref q - ref ref 2 = V + [1 1 1] (cos w t V o - ref + sin w t V o - ref ) ref 3 d q wt avec : P(w ,t) = cosw t - sin w t sin cos w t
(III.21)
avec C23 la transformation de Concordia inverse (Cf. (II.20), chapitre II). Ces mthodes se compliquent en prsence des composantes harmoniques, car il faut appliquer la transformation de Park directe et inverse pour chaque frquence et procder au filtrage des signaux frquence autre que celle de la transformation de Park. A cela nous ajoutons le fait que ces mthodes ont besoin dune trs bonne connaissance de la tension du rseau et quelles pourraient occasionner un surdimensionnement du filtre actif srie, sil ny avait pas de conservation de langle de la composante directe de tension lors de la rgulation de la tension de la charge. II.3.1.3 Analyse des performance s des mthodes didentification tudies Les mthodes didentifications tudies sont gnralement mises en uvre pour compenser les creux et les dsquilibres de tension. La prsence dans la tension du rseau lectrique des composantes harmoniques complique la ralisation de ces mthodes qui deviennent dans ce cas dpendantes du spectre harmonique de la tension du rseau. Afin de surmonter les problmes introduits par ces mthodes didentification, ce mmoire propose une nouvelle mthode didentification base sur lemploi dun systme base de PLL.
90
-+
Transformation 2 3
Vr3
+ + -
Vdif1
+ +
Vdif2
Vref1
+ +
Vdif3
Vref2
+ +
Vref3
Fig. III.6 : Mthode didentification base sur lemploi dun systme base de PLL Afin de mesurer lefficacit de cette mthode didentification, trois cas de creux de tension, monophas, biphas et triphas, ont t simuls et analyss. II.3.2.1 Identification des creux d e tension monophass Dans ce cas, un creux de tension monophas de 80% de la tension nominale (Vn) sur la phase 1 a t simul. Les rsultats de simulation de la tension du rseau lectrique (Vs1,2,3), de la composante directe de tension (Vd1,2,3), de la diffrence de ces deux tensions (Vdif1,2,3), de la tension rgulatrice (Vr1,2,3), et de la tension de rfrence (Vref1,2,3) sont donnes sur la Fig. III.7. A partir de la Fig. III.7 et dans lobjectif de retrouver la tension nominale, les tensions de rfrence (Vref1=0,8Vn) et (Vref2=Vref3=0) sont dtermines. La bonne prcision de cette mthode didentification peut tre constate.
91
(V)
Vd1,2,3
(V)
Vdif1,2,3
(V)
Vr1,2,3
(V)
Vref1,2,3
t (s)
Fig. III.7 : Identification dun creux de tension monophas de 80% (Lgende - phase 1 : ., phase 2 : - - - - -, phase 3 : -..-..-..-) II.3.2.2 Identification des creux d e tension biphass Dans ce cas, un creux de tension biphas de 50% de la tension nominale (Vn) sur les phases (phase 1 et phase 2) a t simul. Les mmes grandeurs que prcdemment sont prsentes sur la Fig. III.8.
Vs1,2,3
(V)
(V)
Vd1,2,3
(V)
Vdif1,2,3
(V)
Vr1,2,3
(V)
Vref1,2,3
t (s)
Fig. III.8 : Identification dun creux de tension biphas de 50% (Lgende - phase 1 : ., phase 2 : - - - - -, phase 3 : -..-..-..-) A partir de la Fig. III.8 et dans lobjectif de rtablir la tension nominale, les tensions de rfrence, (Vref1= Vref2=0.5Vn) et (Vref3=0), sont calcules, ce qui prouve la bonne prcision de cette mthode didentification.
92
(V)
(V)
Vd1,2,3
(V)
Vdif1,2,3
(V)
Vr1,2,3
(V)
Vref1,2,3
t (s)
Fig. III.9 : Identification dun creux de tension triphas de 30% (Lgende - phase 1 : ., phase 2 : - - - - -, phase 3 : -..-..-..-) Dans lobjectif de reconstruire la tension nominale, la Fig. III.9 montre travers les tensions de rfrence calcules (Vref1=Vref2=Vref3=0,3Vn), la bonne prcision de cette mthode didentification. Le faible retard de rtablissement observ sur les courbes de Vdif1,2,3 est un problme numrique propre aux conditions initiales de la simulation. La rapidit de cette mthode est traite par la suite dans ce chapitre (Cf. II.6).
93
B2(s)
Vref
+-
Rgulateur
G(s) Onduleur
B1(s)
++
1/A(s)
Vinj
II.4.1 Objectifs gnraux Les objectifs de rgulation sont la robustesse du rgulateur en stabilit et en performance. La robustesse en stabilit permet de conserver la stabilit du systme en boucle ferme dans le cas o les paramtres du systme dorigine changent. La robustesse en performance se mesure par une poursuite rapide et prcise des rfrences. La robustesse des lois de commande (rgulateurs) a t examine par la suite dans ce chapitre par le biais de simulations. II.4.2 Application des rgula teurs classiques Dans lobjectif dassurer une bonne qualit de compensation des perturbations de tension dans des cas de fonctionnements extrmes du rseau lectrique, les rgulateurs classiques de type PID et RST ont t intgrs dans la boucle de rgulation de tension du filtre actif srie. II.4.2.1 Rgulateur RST pour la sservissement de tension du F.A.S Le rgulateur RST, dj introduit (Cf. chapitre II, III.2), a pour objectif dassurer la robustesse du filtre actif srie en stabilit et en performance. En effet, le changement des valeurs des lments passifs du filtre de sortie, caus par le vieillissement, ainsi que le modle des transformateurs dinjection qui na pas t pris en compte dans lquation (III.3), entrane une modification des paramtres du systme dorigine (B1/A). De plus, le changement de la rfrence poursuivre, en frquence et en amplitude, ne doit pas perturber la prcision et la rapidit de la poursuite. Cela peut tre assur par la robustesse en performance du rgulateur. En rfrence la Fig. II.15, du chapitre II, et la Fig. III.10, la boucle de rgulation utilisant une structure RST peut tre donne par la Fig. III.11.
94
B2(s)
Vref
T(s)/R(s)
+-
R(s)/S(s)
B1(s)
++
1/A(s)
Vinj
Fig. III.11 : Schma de la rgulation de la tension du F.A.S par structure RST A partir de la Fig. III.11, on obtient lquation caractrisant la poursuite :
Vinj (s) = (T. B1 )(s) (S. B2 )(s) Vref (s) + I L (s) (S. A + B1 .R)(s) (S. A + B1 .R)(s)
(III.22)
Le premier terme de la relation (III.22) reprsente la fonction de transfert en boucle ferme en termes de poursuite, tandis que le deuxime terme reprsente la fonction de transfert entre les perturbations (causes par le courant de la charge protger IL ) et la sortie Vinj. La frquence de coupure de la boucle de rgulation doit tre choisie pour satisfaire deux critres :
dune part, elle doit tre plus grande que lharmonique le plus lev prsent dans le signal de rfrence, pour ne pas influencer la compensation des tensions harmoniques, dautre part, elle doit tre suffisamment loin de la frquence de commutation pour que lon puisse modliser londuleur par un gain simple (unitaire dans ce cas).
En ajoutant aux objectifs de rgulation une compensation en boucle ferme des perturbations engendres par IL et en prenant en compte les critres du choix de la frquence de coupure de la boucle de rgulation, les polynmes R(s), S(s) et T(s) sont choisis suivant les procdures dcrites au chapitre II, III.2.2. Etant donn notre systme dorigine, reprsent par les quations (III.3) et (III.4), du deuxime ordre (n=2), et en choisissant un rgulateur propre pour diminuer leffet de dphasage entre les grandeurs injectes et rfrences (dj introduit dans le chapitre II), le polynme de stabilit arbitraire sera dordre 4 :
deg(D)=2.n=4.
Le degr des polynmes P(s) et C(s) est de 2 :
Les ples du systme dorigine subissent les changements dcrits dans le chapitre II sur la Fig. II.14. Ils sont ensuite projets sur labscisse (-1/Tc) pour former le polynme P(s), et sur labscisse (-1/Tf) pour former le polynme C(s).
95
deg(R)=deg(S)=n=2
Le polynme T(s), quant lui, est choisi comme un gain simple, tel que (Vinj/Vref)(s) = 1, quand la variable symbolique (s) tend vers 0. Le problme de cette mthode, comme toute mthode linaire, est le dphasage de tension entre la rfrence poursuivre (Vref) et la sortie de la boucle (Vinj). Ce dphasage augmente, dune part avec lordre du systme en boucle ferme (cest pourquoi on a choisi un rgulateur propre et non pas strictement propre), dautre part avec la frquence des signaux de rfrence poursuivre (Cf. chapitre II, III.4.1). II.4.2.2 Rgulateur PID pour las servissement de la tension du F.A.S Les rangs des harmoniques de tension poursuivre pouvant aller au del de 20 et comme lordre du systme en boucle ferme (Vinj/Vref) est de 4, la mthode de la structure RST, cause du problme de dphasage, ne peut satisfaire les performances requises par les normes internationales. Afin de rsoudre ce problme, dune faon partielle, un rgulateur de type PID a t utilis. Le schma de rgulation de la tension du filtre actif srie est donn par la figure suivante :
IL
B2(s)
Vref
+-
PID(s)
B1(s)
++
1/A(s)
Vinj
Fig. III.12 : Schma de la rgulation de la tension du F.A.S par PID En compensant le dnominateur du systme dorigine A(s) par le numrateur de la fonction de transfert du rgulateur NPID et en se basant sur la Fig. III.12, on dtermine lquation caractrisant la poursuite :
K .bo .N B2 ( s ).DPID ( s ) Ti .Td V + I V = inj N K bo N ref A( s ).DPID + B1 ( s ).N PID L 2 s + s+ Td Ti .Td avec H PID = N PID DPID Td s 1 = K 1 + + Ti s T 1+ d N s
(III.23)
96
97
(III.24)
(III.25)
avec D(s)=(A.S+B1.R)(s), le polynme de stabilit arbitraire. Afin de minimiser lerreur e(s) et de la faire tendre vers zro, le numrateur de lquation (III.25) doit tre nul. Puisque les polynmes R(s) et S(s) ont dj t choisis, le polynme T(s) doit remplir cette condition. De plus, lordre du polynme T(s) doit tre choisi de faon ce que la fonction de transfert T(s)/R(s), prsente sur la Fig. III.11, soit ralisable (causale). La stabilit du systme, dfinie par une structure RST classique, nest pas modifie si lon utilise un polynme T(s) stable [DELA96]. Dans le cas tudi, les polynmes R(s) et S(s) tant du deuxime ordre, le polynme T(s) est galement du deuxime ordre, soit le plus grand ordre applicable T(s) pour satisfaire le critre de causalit de T(s)/R(s). Le polynme T(s) devient ainsi :
T ( s ) = t1s 2 + t2 s + t3
(III.26)
En reportant lquation (III.26) dans celle de (III.25) et pour s=jw, le numrateur de la relation (III.25) devient :
(T .B1 - D)(jw ) = Re ( w ) + jI m ( w )
(III.27)
Afin de faire tendre lquation (III.27) vers zro, les deux parties relle Re(w) et imaginaire Im(w) doivent tre nulles. Cette contrainte permet davoir deux quations pour chaque frquence correspondant aux zros placs. Ces frquences, au nombre de deux, doivent tre choisies en prenant en compte le nombre des paramtres inconnus dans la relation (III.27), soit t1,t2 et t3. On choisit de placer les zros proches des ples lents de la boucle ferme. II.4.4 Synthse des rgulateu rs RST, PID et RSTam II.4.4.1 Comparaison des rgulat eurs en termes de poursuite Afin de comparer les trois mthodes de rgulation, PID, RST et RSTam, en termes de poursuite, en gain et en phase, nous avons tabli des rsultats que nous avons regroups dans la Fig. III.13 et le Tab. III.1. Dans la Fig. III.13, le gain et la phase des fonctions de transfert en boucle ferme (Vinj/Vref) sont tracs. Dans le Tab. III.1, ces mmes rsultats ont t regroups sous forme de tableau pour des multiples de la frquence fondamentale Ff allant de Ff 25 Ff. En gnral, dans un rseau lectrique, les tensions harmoniques de rang plus lev sont ngligeables. En terme de poursuite en gain, on remarque, partir du Tab. III.1 et de la Fig. III.13, que le gain est satisfaisant dans les trois mthodes sur toute la bande de frquence considre.
98
Fig. III.13 : Fonction de transfert (Vinj/Vref) avec les rgulateur PID, RST et RSTam F (Hz) PID RST RSTam G PH G PH G PH 50 1 -0.7 1 -1,3 1 250 1 -3 1 -6 1 350 1 -4 1 -9 1 -0,5 550 1 -7 1 -14 1 0,2 650 1 -8 1 -17 1 0,7 850 1 -11 1,05 -22 1,06 1,5 950 1 -12 1,05 -24 1,07 1,5 1150 1 -15 1,1 -29 1,2 0,5 1250 1 -16 1,11 -32 1,2 -0,8
-0,2 -0,5
Tab. III.1: Gain et phase du systme en boucle ferme utilisant PID, RST et RSTam
99
Fig. III.14 : Gain de la fonction de transfert perturbations-sortie (Vinj/IL) A partir de la Fig. III.14, on remarque que les perturbations la frquence fondamentale de IL sont attnues de 18 dB dans le cas du PID et par un gain de 55 dB dans les cas des rgulateurs RST et RSTam. Cela veut dire que ces perturbations nont aucune influence sur la boucle de rgulation dans le cas des rgulateurs RST et RSTam. Par contre, une forte perturbation de ce type peut influencer la qualit de poursuite dans le cas du rgulateur PID. La compensation des perturbations est assure par les rgulateurs prsentant un gain relativement lev la frquence de perturbation (50 Hz), comme le montre la Fig. III.15. En effet, les perturbations sont gnralement attnues proportionnellement au gain du rgulateur.
Fig. III.15: Gain des fonctions de transfert des rgulateurs employs 100
101
(III.28)
Cette mme tension du rseau a t dcompose en trois squences directe (Vd), inverse (Vi) et homopolaire (Vo). La puissance apparente de la charge protger aprs compensation (Sch), pour un courant (IL) et une tension (Vch) dont langle est (dch), est donne par lquation (III.29).
Sch = 3 Vch I L
(III.29)
La puissance apparente du filtre actif srie (Sf) est compose de trois puissances apparentes : directe (Sd-inj), inverse (Si) et homopolaire (So). La puissance directe (injecte) (Sd-inj) est due la tension directe injecte par le filtre actif (Vd-inj). Cette tension est la diffrence entre la tension directe dsire (Vd-des, dd-des) et la tension directe du rseau lectrique donne par le systme base de PLL (Vd, dd). La puissance inverse (Si) est due la tension inverse injecte par le filtre actif (Vi) tandis que la puissance homopolaire (So) est due la tension homopolaire injecte par le filtre actif (Vo), comme le dcrivent les quations suivantes :
r r S d - inj = 3 Vd - des - Vd I L Si = 3 Vi I L S o = 3 Vo I L
(III.30)
En se basant sur lquation (III.30), la puissance totale apparente du filtre actif srie peut tre donne par lquation (III.31).
Sf =
(S
d - inj
) + (S ) + (S )
2 2 i o
(III.31)
102
(III.32)
Des quations (III.29) et (III.31), le rapport, tmo-ph, entre la puissance du filtre actif srie et la puissance de la charge protger est donn par lquation suivante :
Sf S ch
t mo - ph =
r r Vd - des - Vd r = Vch
Vi Vo + + V V ch ch
(III.33)
A partir de cette quation, deux cas avec et sans lutilisation de la PLL, peuvent tre abords. En effet, le systme base de PLL donne la possibilit de synchroniser sur la composante directe de tension en conservant langle dd. Alors que sans la PLL (ou dautres techniques fournissant cet angle), un autre angle de synchronisation dd-des doit tre choisi lors de la rgulation de la tension de charge. II.5.1.1 Cas nutilisant pas un sys tme base de PLL Dans ce cas le rapport des puissances tmo-ph peut scrire :
t mo - ph V V V = 1 - d ( d d - d d - des ) + i + o V V Vch ch ch
2 2 2
(III.34)
II.5.1.2 Cas utilisant un systme base de PLL Dans ce cas, langle de la composante directe donn par la PLL et langle dsir tant gaux, dd-des=dd, le rapport des puissances tmo-ph devient :
t mo - ph
V V V = 1 - d + i + o V V V ch ch ch
(III.35)
A partir des quations (III.34) et (III.35), les courbes du rapport des puissances, du filtre actif srie et de la charge protger, avec et sans synchronisation sont donnes sur la Fig. III.16. 103
......... b = 60 + + + b = 30 ****
1
b = 0 b = 0
tmo-ph 0.6
0.4 0.2 0 0
2
1 1,2 2
0.2
0.4
0.6
0.8
Fig. III.16 : Rapport de puissances de F.A.S et de charge protger pour le cas des creux monophass avec synchronisation (2) et sans synchronisation (1) Dans cette figure, (1-k) est la profondeur du creux de tension. A partir de cette figure, on observe que la puissance maximale injecte par le filtre actif srie, avec et sans synchronisation, est (Sf=1,18 Sch) pour le cas extrme (1-k= 0) et (b=180). Par contre, l'avantage de lutilisation de la PLL, du point de vue nergtique, est plus significatif pour lensemble des angles de dphasage (0<b<180). II.5.2 Compensation des creu x de tension biphass La deuxime tude de la puissance apparente concerne la compensation des creux de tension biphass. Dans ce cas, le filtre actif srie injecte au rseau lectrique deux composantes de tension directe et inverse. Sachant que les angles de dphasage de tension sont dans ce cas b1=b2=b et b3=0 et que les coefficients damplitude ont le valeurs k1=k2=k et k3=1, les tensions injectes au rseau sont :
2 k V Vd = 3 (cos b + j sin b + 1) V = k V (1 - cos b - j sin b ) i 3
(III.36)
Des quations (III.29) et (III.31), le rapport, tbi-ph, entre la puissance du filtre actif srie et la puissance de la charge protger est donn par lquation suivante :
r r Vd - des - Vd r = = Sch Vch
Sf
t bi - ph
Vi + V ch
(III.37)
En se basant sur cette quation, deux cas, avec et sans lutilisation de la PLL, peuvent tre discuts.
104
(III.38)
II.5.2.2 Cas utilisant un systme base de PLL Dans ce cas, langle de la composante directe donn par la PLL et langle dsir tant gaux, dd-des=dd, le rapport des puissances tbi-ph devient :
t sbi - ph
V V = 1 - d + i V V ch ch
(III.39)
A partir des quations (III.38) et (III.39), les courbes du rapport des puissances, du filtre actif srie et de la charge protger, avec et sans synchronisation sont donnes sur la Fig. III.17.
1.5
1
2 1
****
tbi-ph
0.5
1,2 2 2
0 0
0.2
0.4
0.6
0.8
Fig. III.17 : Rapport de puissances de F.A.S et de charge protger pour le cas des creux biphass avec synchronisation (2) et sans synchronisation (1) Dans ce cas et partir de la figure prcdente, la puissance maximale du filtre actif avec synchronisation est (Sf=1,2 Sch) pour le cas extrme (1-k= 0,5) et (b=180). Sans synchronisation, la puissance maximale du filtre actif srie est de (Sf=1,5 Sch) pour le cas extrme (1-k= 0) et (b=180). II.5.3 Compensation des creu x de tensions triphass Dans la troisime tude, la puissance apparente du filtre actif srie, compensant les creux de tension triphass, a t calcule. Dans ce cas, le filtre actif srie injecte uniquement au rseau lectrique la composante directe de tension. Sachant que les angles de dphasage de tension sont
105
(III.40)
Des quation (III.29) et (III.31), le rapport, ttri-ph, entre la puissance du filtre actif srie et la puissance de la charge protger est donn par lquation suivante :
t tri - ph
(III.41)
En se basant sur cette quation, deux cas, avec et sans lutilisation de la PLL, peuvent tre abords. II.5.3.1 Cas nutilisant pas un sys tme base de PLL Dans ce cas le rapport des puissances tbi-ph peut scrire :
t tri - ph = 1 -
Vd ( d d - d d - des ) Vch
(III.42)
II.5.3.2 Cas utilisant un systme base de PLL Dans ce cas, langle de la composante directe donn par la PLL et langle dsir tant gaux, dd-des=dd, le rapport des puissances ttri-ph devient :
V t stri - ph = 1 - d V ch
(III.43)
A partir des quations (III.42) et (III.43), les courbes du rapport des puissances, du filtre actif srie et de la charge protger, avec et sans synchronisation sont donnes sur la Fig. III.18. Dans ce cas et partir de la figure prcdente, la puissance maximale injecte par le filtre actif srie avec synchronisation est Sf =Sch, pour le cas extrme 1-k=1. Sans synchronisation, la puissance du filtre actif srie est de Sf=2 Sch, pour le cas extrme 1-k= 0 et b=180. II.5.4 Analyse conomique d e la compensation A partir des trois cas de creux de tension analyss, on peut remarquer que : il est indispensable de synchroniser la tension de la charge avec la composante directe de tension lors de la rgulation de la tension de charge. En effet, ne pas synchroniser peut augmenter la puissance du filtre actif srie jusqu' deux fois celle de la charge protger, pour des cas extrmes, la puissance apparente injecte par le filtre actif srie peut tre trs leve pour les cas extrmes de creux de tension. Cela, du point de vue conomique, peut conduire limiter la compensation des creux de tension jusqu 30%, comme prvu par le contrat EMERAUDE.
106
ttri-ph 1
0.8 0.6 0.4 0.2 0
1,2
0.2
0.4
0.6
0.8
Fig. III.18 : Rapport de puissances de F.A.S et de charge protger pour le cas des creux triphass avec synchronisation (2) et sans synchronisation (1)
107
es
Rs
Ls
Vs2
Lt,Rt
Vch2 Vch3
Vdc
Cdc1 Cdc2
Filtre de Sortie
Lsf,Rsf
Vinj1
Vinj2
Vinj3
Vs1,2,3
PLL
Vinj1,2,3 Vd-des Vdc
Fig. III.19 : Structure gnrale du filtre actif srie au sein du systme tudier Les valeurs des lments caractrisant la structure tudie sont regroupes dans le Tab. III.2. Rseau lectrique 1 MVA, 230 V, ucc=4% (Rs, Ls) Charges protger 100 kVA, cos j = 0,8, 230 V Filtre actif srie Vdc (source de tension continue) (840 V) (Lsf, Rsf) (500 mH, 31 mW) (Csf) (200 mF) Vd-des 325 V Frquence de commutation 12 kHz Transformateurs dinjection de tension 200 kVA, 230 V, ucc=2% (Rt, Lt) (1,2 mW, 17 mH) Tab. III.2: Valeurs des lments constituant la structure tudie Plusieurs simulations travers deux analyses, temporelle et spectrale, ont t effectues sur SABER afin de valider le fonctionnement du filtre actif srie sous des conditions relles, et parfois exagres, des rseaux lectriques. Lefficacit et la rapidit de la mthode didentification, base de PLL, et celle de rgulation de tension, RSTam, seront analyses travers plusieurs simulations du filtre actif srie compensant des creux de tension, des dsquilibres de tension et des harmoniques de tension. Les avantages ainsi que les limites des autres mthodes de rgulation de tension (RST, (0,25 mW, 19,4 mH)
108
Charge Protger
IL1
Vs1
Vch1
Chapitre III : Filtre actif srie et combinaison parallle-srie actifs ___________________________________________________________________________ PID) seront compars avec la mthode RSTam dans lobjectif dun choix final de la mthode de rgulation de la tension du filtre actif srie.
II.6.1 Compensation des creu x de tension Dans cette partie de simulation, on va analyser la robustesse en termes de rapidit et de prcision du filtre actif srie compensant les creux de tension monophass, biphass et triphass. II.6.1.1 Compensation des creux de tension monophass Dans ce cas, un filtre actif srie, compensant un creux de tension monophas dune profondeur de 80% et dune dure de 4 priodes sur la Phase 1, a t tudi. Les analyses temporelles et spectrales dun filtre actif srie, employant la mthode didentification base dun systme PLL ainsi que les trois mthodes de rgulation de tension (PID, RST et RSTam), seront illustres. Nous comparerons la robustesse en rapidit et en prcision dans le cas de ces trois mthodes de rgulation. En effet, avec le changement du point de fonctionnement du rseau lectrique (avant, pendant et aprs le creux de tension), la mthode didentification ( travers Vref1) et les mthodes de rgulation ( travers Vinj1) doivent poursuivre avec prcision et rapidit des valeurs variant entre 0% et 80% de la tension nominale Vn. II.[Link] Rgulation de tension par PID La Fig. III.20 montre lanalyse temporelle du filtre actif srie employant le rgulateur PID et compensant ce creux de tension monophas. Dans cette figure, la tension du rseau (Vs1,2,3), le courant de la charge protger (IL1,2,3), les tensions superposes de la Phase 1 identifie (Vref1) et injecte (Vinj1) et la tension de la charge aprs compensation (Vch1,2,3) sont prsentes.
Vs1,2,3 (V)
IL1,2,3 (A)
(V)
Vref1 Vinj1
Vch1,2,3 (V)
t (s)
Fig. III.20 : Analyse temporelle de la compensation du creux de tension monophas (PID) (Lgende - phase 1 : ., phase 2 : - - - - -, phase 3 : -..-..-..-) Cette figure montre la rapidit des mthodes didentification, dun systme base de PLL, et de rgulation PID face un changement brutal de la consigne. En effet, suivant le changement du point de fonctionnement du rseau lectrique, la tension identifie (Vref1) et celle injecte (Vinj1) rpondent rapidement ces variations brutales. Cette rapidit de compensation du creux de tension permet la circulation dun courant quilibr malgr ces changements. 109
Vref1 (V)
Vinj1 (V)
Vch1 (V)
f (Hz)
Fig. III.21 : Analyse spectrale de la compensation du creux de tension monophas (PID) Cette figure, o la tension identifie (Vref1) permet la tension du rseau (Vs1) datteindre la valeur crte de la tension de charge (Vch1) dfinie par (Vd-des=325 V), tmoigne de la bonne prcision de la mthode didentification propose. De plus, une tension de charge protger aprs compensation trs proche de la tension impose par (Vd-des), confirme la possibilit dutiliser le rgulateur PID dans ce type de compensation. II.[Link] Rgulation de tension par RST Les analyses temporelle et spectrale du filtre actif srie employant le rgulateur RST et compensant ce creux de tension monophas ont t effectues sur les Fig. III.22 et Fig. III.23. La Fig. III.22 montre lanalyse temporelle. A partir de cette figure, la rapidit de rponse des tensions identifie (Vref1) et injecte (Vinj1), ces variations brutales du point de fonctionnement du rseau lectrique, est une preuve de la robustesse en rapidit des mthodes employes. Le courant du rseau lectrique est toujours quilibr grce cette rapidit de compensation du creux de tension. La Fig. III.23 montre lanalyse spectrale de la Phase 1 afin dvaluer la prcision des mthodes didentification et de rgulation utilises. Cette figure dmontre galement la bonne prcision des mthodes didentification et de rgulation proposes. La petite diffrence entre les valeurs crtes souhaites et obtenues de la tension de charge peut tre explique par un dphasage de 1,287 (Cf. Tab. III.1).
110
IL1,2,3 (A)
(V)
Vref1 Vinj1
Vch1,2,3 (V)
t (s)
Fig. III.22 : Analyse temporelle de la compensation du creux de tension monophas (RST) (Lgende - phase 1 : ., phase 2 : - - - - -, phase 3 : -..-..-..-)
Vs1 (V)
Vref1 (V)
Vinj1 (V)
Vch1 (V)
f (Hz)
II.[Link] Rgulation de tension par RSTam Le comportement du filtre actif srie compensant toujours le mme creux de tension monophas, et employant cette fois-ci le rgulateur RSTam sera simul sur les Fig. III.24 et Fig. III.25. La Fig. III.24 montre lanalyse temporelle. Face aux mmes changements brutaux de la consigne, la rapidit de la rponse des mthodes didentification et de rgulation de tension RSTam sont montres sur la Fig. III.24. Le courant de la charge protger est donc quilibr grce cette rapidit de compensation.
111
IL1,2,3 (A)
(V)
Vref1 Vinj1
Vch1,2,3 (V)
t (s)
Fig. III.24 : Analyse temporelle de la compensation du creux de tension monophas (RSTam) (Lgende - phase 1 : ., phase 2 : - - - - -, phase 3 : -..-..-..-) La Fig. III.25 montre lanalyse spectrale de la Phase 1 afin dvaluer la prcision des mthodes didentification et de rgulation proposes.
Vs1 (V)
Vref1 (V)
Vinj1 (V)
Vch1 (V)
f (Hz)
Fig. III.25 : Analyse spectrale de la compensation du creux de tension monophas (RSTam) Dautres rsultats de simulation du filtre actif srie compensant des creux de tension biphas et triphas sont reports dans lannexe-II. Dans cette annexe, ltude des performances de la mthode didentification base du systme PLL et des mthodes de rgulation de tension (PID, RST et RSTam) a t effectue.
112
(V)
Vref1 Vinj1
(V)
Vref2 Vinj2
(V)
Vref3 Vinj3
Vch1,2,3 (V)
t (s)
Fig. III.26 : Analyse temporelle de la compensation du dsquilibre et des harmoniques de tension (PID) A partir de cette figure, on peut observer que le dphasage entre les tensions de rfrences (Vref1, Vref2, Vref3) et celles injectes (Vinj1, Vinj2, Vinj3) peuvent tre ngliges conformment au Tab. III.1. Lallure sinusodale quilibre de la tension de charge aprs compensation prouve le bon fonctionnement de la mthode didentification ainsi que celui de la mthode de rgulation employes. La Fig. III.27 montre lanalyse spectrale de la Phase 1 de la tension du rseau (Vs1) et celle de la charge aprs compensation (Vch1). Le taux individuel de distorsion harmonique tant respectivement pour les trois phases de 19,4%, 15,89% et 14,16%, contre 0,94% aprs compensation. 113
(V)
Vch1
(V)
f (Hz)
Fig. III.27 : Analyse spectrale de la compensation des harmoniques de tension (PID) Afin dvaluer avec prcision la qualit du rquilibrage du filtre actif srie intgrant le rgulateur PID, une autre analyse spectrale des tensions du rseau (Vs1, Vs2, Vs3) et de la charge (Vch1, Vch2, Vch3) est donne sur la Fig. III.28.
(V) Vs1
(V)
Vs2
(V)
Vs3
(V)
Vch1
(V)
Vch2 Vch3
(V)
f (Hz)
Fig. III.28 : Analyse spectrale de la compensation du dsquilibre de tension (PID) A partir de cette figure on peut dterminer avec prcision le taux de dsquilibre de tension qui est de 9%, ct rseau, et de 0,047%, ct charge, aprs compensation. Cela confirme la bonne qualit de compensation du dsquilibre de tension dj observe sur lanalyse temporelle.
114
(V)
(V)
(V)
(V)
t (s)
Fig. III.29 : Analyse temporelle de la compensation du dsquilibre et des harmoniques de tension (RSTam) A partir de cette figure, on peut observer que le dphasage entre les tensions de rfrences (Vref1, Vref2, Vref3) et celles injectes (Vinj1, Vinj2, Vinj3) est vraiment ngligeable, conformment au Tab. III.1. La tension de la charge protger a une allure sinusodale et quilibre. Langle de la tension de la charge, aprs compensation et rgulation, est celui de la composante directe de tension dd=3,2 pour la Phase 1. Cela tmoigne dune trs bonne prcision de la mthode didentification base de PLL et dun dphasage quasi nul ralis grce la mthode de rgulation RSTam. La Fig. III.30 montre lanalyse spectrale de la Phase 1 de la tension du rseau (Vs1) et de la tension de la charge aprs compensation (Vch1). Le taux de distorsion harmonique individuel tait respectivement pour les trois phases de 19,4%, 15,89% et 14,16% contre 0,042% aprs compensation. Afin de dterminer avec prcision la qualit du rquilibrage du filtre actif srie intgrant le rgulateur RSTam, une deuxime analyse spectrale des tensions du rseau (Vs1, Vs2, Vs3) et de la charge (Vch1, Vch2, Vch3) est montre sur la Fig. III.31. La Fig. III.31 montre que le taux de dsquilibre de tension qui est de 13,6%, ct rseau, devient 0,015%, ct charge, aprs compensation. Cela confirme la bonne qualit de compensation du dsquilibre observe sur lanalyse temporelle.
115
(V)
Vch1
(V)
f (Hz)
(V)
Vs2
(V)
(V)
(V)
(V)
Vch3
f (Hz)
Fig. III.31 : Analyse spectrale de la compensation du dsquilibre de tension (RSTam) II.6.3 Evaluation des mthod es didentification et de rgulation Tous les rsultats de simulation de la compensation des perturbations de tension (creux, dsquilibre et harmonique) montrent que la mthode didentification base de PLL, propose dans ce mmoire, est une mthode robuste sadaptant rapidement et avec prcision aux variations brutales du point de fonctionnement du rseau lectrique. Quant aux mthodes de rgulation de tension du filtre actif srie, les rsultats de simulation de la compensation des creux de tension, monophas, biphas et triphas et du dsquilibre de tension, montrent que les trois mthodes de rgulation (PID, RST, RSTam) sont valides. En effet et conformment au Tab. III.1, pour un gain quasi unitaire avec un dphasage ngligeable 50 Hz, les trois mthodes de rgulation peuvent tre utilises pour compenser ces types de perturbations. 116
117
DEPOLLUTION
Dans ce chapitre, nous avons prsent le filtre actif srie. Ce dernier est une solution de dpollution des rseaux lectriques des perturbations en tension, comme les tensions harmoniques, les tensions dsquilibres et les creux de tension. Nous avons choisi un filtre de sortie du deuxime ordre dans la structure du filtre actif srie, de manire empcher les composantes dues la frquence de commutation de se propager sur le rseau lectrique. Une nouvelle mthode didentification a t propose de manire identifier toutes sortes de perturbations de tension du rseau lectrique. Cette mthode didentification est base sur lemploi dun systme base de PLL, dj intgr dans la mthode didentification des courants perturbateurs. Elle rduit considrablement le nombre de calculs par rapport aux autres mthodes proposes et ne ncessite pas une connaissance prcise du rseau lectrique. Une tude comparative, avec et sans PLL, montre que la puissance apparente ncessaire au filtre actif srie est minimale sil y a une synchronisation en phase avec la composante directe de la tension lors de la rgulation de tension de la charge protger. Afin de raliser les objectifs de rgulation que nous nous tions fixs, nous avons tudi trois rgulateurs. Un rgulateur RST a t employ dans la structure du filtre actif srie de manire compenser les creux de tension et les dsquilibres de tension. Cependant, le problme de dphasage, dj observ dans la rgulation du courant, a t de nouveau identifi. Ce problme de dphasage dgrade la qualit de compensation du filtre actif srie dans le cas de perturbations dues aux tensions harmoniques. Afin de rsoudre ce problme, les principes du rgulateur RSTam, dj introduit dans le chapitre II, et du rgulateur PID, ont t appliqus. Une tude comparative en termes de suivi de consigne entre les rgulateurs PID, RST et RSTam a t effectue afin de mesurer le champ dapplication de chaque mthode. Cette tude a montr que les trois mthodes de rgulation taient valides pour la compensation des creux et dsquilibres de tension. De plus, le rgulateur PID est utilisable uniquement pour compenser les tensions harmoniques de rangs et damplitudes relativement faibles. Par contre, la mthode propose RSTam prvoit la compensation de toutes les tensions harmoniques sans contraintes, ni sur leur amplitude, ni sur leur frquence. Enfin, plusieurs simulations ont t effectues sur SABER afin de valider la fonctionnalit des mthodes du filtre actif srie compensant toutes les perturbations de tension. Suite au travail men prcdemment dans ce mmoire et portant sur la contribution aux amliorations des performances des filtres actifs parallles et sries, il nous a paru intressant dtudier les amliorations des performances de lUPQC regroupant ces deux types des filtres actifs.
118
Filtre de sortie
Onduleur
Elment de stockage
Onduleur
Filtre de sortie
Fig. III.32 : Structure gnrale de lUPQC Suivant lobjectif de dpollution du filtre actif srie, deux emplacements de lUPQC sont possibles : Le filtre actif srie est install en amont du filtre actif parallle, (Cf : Fig. I.8), afin de dpolluer la source des tensions perturbatrices, Le filtre actif srie est install en aval du filtre actif parallle, Fig. III.33 ci-dessous, afin disoler la charge de la source perturbe.
Dans ce travail, le filtre actif srie sera install ct charge afin de ne pas le surdimensionner. En effet, un filtre actif srie install ct poste de transformation devra tre dimensionn pour supporter la totalit des courants de tous les utilisateurs connects au mme point de raccordement. Afin de rduire le nombre de calculs effectuer, la combinaison parallle-srie actifs va bnficier du fonctionnement commun de certains blocs. En effet, le systme base de PLL, utilis pour lidentification des courants et des tensions perturbateurs, le systme de stockage dnergie ainsi que la rgulation de la tension continue, vont fonctionner dune faon commune pour les deux parties de lUPQC.
119
Is1 Is2
Rc Lc
es
Rs
Ls
Vdc
Cdc1
Vinj2
Vinj3
Cf Rf
Cdc2
Csf
Vs1,2,3
PLL
Is1,2,3 Iinj1,2,3 Vdc
Fig. III.33 : Structure tudie de lUPQC Les valeurs des lments de la structure tudie sont regroupes dans le Tab. III.3. Rseau lectrique 1 MVA, 230 V, ucc=4% (Rs, Ls) Charges protger pont Graetz de 100 kVA, a=5 (Rc, Lc) (100 kVA, ucc=4%) Filtre actif parallle (Lf1, Lf2), (Rf1, Rf2, Rf), (Cf) Frquence de commutation (0,25 mW, 19,4 mH) 100 kVA (14,6 mW, 155,6 mH)
(90 , 70) mH, (5, 5, 16) mW, (200) mF 12 kHz Filtre actif srie
(Lsf, Rsf, Csf) (500 mH, 31 mW, 200 mF) (Rt, Lt) (200 kVA, ucc=2%) (1,2 mW, 17 mH) Vd-des 325 V Frquence de commutation 12 kHz Source continue commune Vdc (Cdc1, Cdc1) 840 V (8,8, 8,8) mF
Tab. III.3 : Valeurs des lments constituant la structure adopte (UPQC) Afin de valider le comportement de lUPQC, un point de fonctionnement du rseau lectrique rassemblant la plupart des perturbations, a t simul sur SABER. 120
Is1,2,3 (A)
Vch1,2,3 (V)
IL1,2,3 (A)
t (s)
Fig. III.34 : Analyse temporelle du fonctionnement de lUPQC A partir de cette figure, on remarque que durant les premires 80 ms, o le F.A.S ne fonctionne pas, le creux de tension influence le courant de la charge et par consquent la forme de la tension continue. Aprs le dbut de la compensation du creux de tension, la forme du courant de la charge et de la tension continue redevient correcte. La Fig. III.34 montre aussi que la tension continue poursuit sa rfrence de 840 V avec rapidit et prcision, grce la boucle de rgulation continue. Les condensateurs tant rechargs en ligne par la filtre actif parallle, leurs valeurs peuvent tre considrablement rduites par rapport aux calculs dcrits par les quations III.1 et III.2 [TO99]. En effet, en se basant sur ces deux quations et pour un creux de tension monophas dune profondeur de 50% et dune dure de 60 ms, chaque condensateur doit avoir une valeur de Cdc1= Cdc2=186 mF. Cette rduction (1/21 dans notre cas) est limite par le courant fondamental, fourni par le rseau lectrique, ncessaire pour recharger ces condensateurs. La tension de la charge redevient quilibre aprs 80 ms ainsi que le courant du rseau est sinusodal durant toute la dure de compensation. Afin de montrer le fonctionnement de chaque partie de lUPQC, deux autres figures, Fig. III.35 et Fig. III.36, sont donnes. La Fig. III.35 montre lanalyse temporelle de la compensation du creux de tension par la partie srie de lUPQC. Dans cette figure, la tension du rseau lectrique (Vs1,2,3), la tension injecte de la Phase 1 (Vinj1), les tensions superposes, injectes aux deux autres phases (Vinj2, Vinj3) et la tension aux bornes de la charge protger (Vch1,2,3) sont donnes.
121
Vinj1, (V)
(V)
Vinj2 Vinj3
Vch1,2,3 (V)
t (s)
Fig. III.35 : Analyse temporelle de la compensation du creux de tension Monophas par la partie srie de lUPQC A partir de cette figure, on remarque que le filtre actif srie arrive compenser ce creux de tension et rguler la tension de la charge avec prcision et ceci presque instantanment. La tension ngligeable injecte par le filtre srie sur les Phase 2 et Phase 3, pour toute la dure de simulation, et sur la Phase 1, avant 80 ms, est due aux dcoupages malgr le filtrage de sortie du deuxime ordre. La Fig. III.36 montre lanalyse temporelle de la compensation par la partie parallle de lUPQC, du courant harmonique et de la puissance ractive du pont Graetz. Dans cette figure, le courant de la charge polluante (Is1,2,3), les courants superposs, Phase 1, de la rfrence (Iref1) et du filtre actif parallle (Iinj1) et le courant du rseau (IL1,2,3) sont donns.
Is1,2,3 (A)
IL1, (A)
t (s)
Fig. III.36 : Analyse temporelle de la compensation du courant harmonique et de la puissance ractive par la partie parallle de lUPQC
122
123
124
125
126
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
I. INTRODUCTION
I.1 Problmatiques de re cherche
Si de nos jours, les solutions de type filtrage actif rentrent de plus en plus dans les plans dinvestissement des industriels franais, cest aussi grce lengagement des socits de conception qui ont russi rduire la fois les contraintes de dveloppement, et les cots lachat de ces dispositifs. A ltranger, comme rcemment en France, on compte aujourdhui plusieurs socits qui fabriquent et commercialisent des filtres actifs. Ces appareils sont essentiellement des filtres actifs parallles appliqus la rduction des cots occasionns par la fourniture en nergie lectrique. Ces cots tant principalement dus aux perturbations harmoniques qui transitent entre le site industriel client et le rseau lectrique du fournisseur. Pour aller dans le sens de lamlioration de ces produits, les concepteurs sintressent de plus en plus la ralisation de filtres actifs universels qui pourraient rsoudre la plupart des problmes provoqus par les perturbations en courant et en tension. De plus, ces dispositifs doivent aussi rpondre des contraintes trs industrielles, telles que la souplesse dutilisation, lapplication sur des gammes de puissance, de tension et de frquence extensibles et la capacit se connecter sur nimporte quel rseau ou systme prexistant. Toutes ces exigences de fonctionnement entranent le dveloppement de nouvelles techniques de compensation par filtrage actif, comme celles que nous avons dj introduites lors des premiers chapitres de cette thse. Pourtant, la conception dun filtre actif de ce type doit prendre en compte des critres de faisabilit exprimentale trs complexe mettre en uvre.
127
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
simulations effectues laide des logiciels SABER et Matlab (toolbox Simulink & Power System Blockset) . Dans la troisime et dernire partie, nous tudierons la faisabilit exprimentale du filtre actif parallle en mode discret. Dans cette tape, deux mthodes de conception de rgulateurs partir du plan z seront abordes. La premire mthode se basera sur la discrtisation des modles conus en continu en utilisant la commande de Matlab C2d (conversion continuous-time models to discrete time). La deuxime mthode permettra de concevoir directement les rgulateurs dans le plan z. Dans cette deuxime mthode, le schma global de la boucle de rgulation du courant, telle quelle sera implante, sera construite et les performances du rgulateur RST en z seront testes en termes de poursuite et de compensation des perturbations. Finalement, nous effectuerons une tude comparative des rgulateurs conus en s et en z. Cette tude comparative consistera chercher un compromis entre les diffrentes contraintes numriques dans lobjectif dtudier la faisabilit exprimentale du prototype de filtre actif parallle universel propos par la socit partenaire.
128
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
Is1 Is2
Rc Lc
es
Rs
Ls
Vs1,2,3
PLL
Is1,2,3 Iinj-1,2,3 Vdc
Fig. IV.1 : Structure gnrale du systme tudi Le filtre de sortie du filtre actif est un filtre passif du troisime ordre en T, reprsent par deux inductances (Rf1,Lf1, Rf2,Lf2) et une capacit Cf en srie avec une rsistance damortissement Rf. Ce filtre sert dune part, connecter londuleur de tension au rseau lectrique, dautre part empcher les composantes dues aux commutations de se propager sur le rseau lectrique.
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
La tension d'alimentation est de 240/415 V dans le cas du rseau lectrique et de 265/375 V dans le cas du groupe lectrogne. Ces tensions de rseau ont t choisies de faon se placer dans le cas le plus dfavorable du cahier des charges concernant la zone d'incommandabilit du filtre actif parallle, la tension continue tant fixe (Cf. Fig. II.3). Dans le cas du rseau lectrique, deux frquences centrales (50 Hz et 60 Hz) seront tudies, alors que dans le cas d'un groupe lectrogne, trois frquences extrmes (45 Hz, 55 Hz et 65 Hz) seront tudies.
130
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
Iref
T(s)/R(s)
+-
R(s)/S(s)
++
1/g(s)
++
B1(s)
++
1/A(s)
Iinj
Fig. IV.2 : Schma de la rgulation par RSTam du courant du F.A.P De la Fig. IV.2, la relation suivante, caractrisant la poursuite de la consigne, peut tre dduite :
I inj (s) =
(T. B1 )(s) (g.S. B2 )(s) I ref (s) + V (s) + (g.S. A + B1 .R)(s) (g.S. A + B1 .R)(s) s (S. B1 )(s) (g.S. B1 )(s) + V (s) + P (s) s (g.S. A + B1 .R)(s) (g.S. A + B1 .R)(s) ond
(IV.1)
Avec 1/g(s)=1/(t.s+1) reprsentant le modle du temps de dcalage prendre en compte lors de limplantation numrique (t=80ms). En effet, le calcul du processus seffectue lintrieur de la priode dchantillonnage et nest appliqu quau dbut de la priode dchantillonnage suivante. Le premier terme de la relation (IV.1) reprsente la fonction de transfert en boucle ferme en termes de poursuite, tandis que le deuxime terme reprsente la fonction de transfert entre les 131
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
perturbations et la sortie Iinj. Ces perturbations, causes par la tension du rseau, et que lon appellera tension naturelle, vont tre compenses par lajout de la tension du rseau Vs la sortie du rgulateur (troisime terme de lquation). A partir de la relation (II.29), on montre que lajout de cette tension, quon va appeler tension compensatrice, est ncessaire pour empcher lappel, par le filtre de sortie, dun courant fondamental au rseau lectrique. Le quatrime terme de la relation (IV.1) montre la fonction de transfert en boucle ferme entre les perturbations causes par londuleur Pond(s) et le courant du filtre actif Iinj(s). En effet, on sait que des perturbations peuvent tre engendres par londuleur de tension du filtre actif. Ces perturbations sont provoques par le temps mort des interrupteurs de londuleur. Ces perturbations sont aussi lies dautres contraintes dimplantation comme nous le montrerons dans la suite de ce mmoire.
I inj (s) =
(IV.2)
avec (B2/A)(s) reprsentant le modle de perturbation. Le gain de cette quation est donn sur la figure suivante :
Fig. IV. 3 : Gain de la fonction de transfert (Iinj/Vs)(s) A partir de cette figure, nous constatons que la moindre tension, fondamentale ou harmonique, du rseau ajoute une composante perturbatrice non ngligeable au courant inject. 132
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
III.2.2 Perturbations de lond uleur- Rsultats exprimentaux Habituellement, les perturbations engendres par un filtre actif parallle sont essentiellement dues la tension perturbe du rseau lectrique. Cependant, des essais exprimentaux pratiqus sur le prototype dvelopp par notre partenaire industriel, pour une charge non-linaire dconnecte, montre que ce dispositif produit des courants harmoniques de tout rang. On observe sur la Fig. IV.4, le courant de la phase 1 du rseau (IL1) et la tension entre les phases 1 et 2 de londuleur (Vf1-2). A noter que le trs faible courant du rseau (30 mA) est consomm par des lments passifs de protection du dispositif.
IL
Vf1-2
Fig. IV.4 : Mesure du courant du rseau et de la tension de londuleur (rsultats exprimentaux) Le spectre de la tension compose de londuleur Vf1-2 est prsent sur la Fig. IV.5 : A partir de la Fig. IV.5, nous constatons que londuleur produit une tension dont le spectre contient des tensions harmoniques de tout rang. Notons que ces rsultats exprimentaux nous ont t fournis par notre partenaire industriel. Afin dvaluer leffet de ces harmoniques sur la qualit de compensation des courants harmoniques, le Tab. IV.4 calcule les courants harmoniques correspondant aux tensions harmoniques simples prpondrantes partir du signal Vf1-2 montres sur la Fig. IV.5. En effet, partir de la Fig. IV.2 nous trouvons que le gain de la fonction de transfert en boucle ouverte (Iinj/Pond) est le mme que celui reprsent en Fig. II.5. A partir du Tab. IV.4 , on observe des courants harmoniques impairs mais galement, et ceci de manire moins habituelle, des courants harmoniques pairs. Les courants impairs sont provoqus par un temps mort des IGBTs qui est limit 2 ms. Les courants pairs sont probablement dus un problme de centrage de la MLI vectorielle employe dans cette manipulation. En effet, cause du temps mort, le centrage de la MLI sur une priode ne peut pas seffectuer avec prcision, ce qui occasionne la gnration, partir de londuleur, des harmoniques pairs non-conventionnels. Il est noter que la prsence des harmoniques multiples de trois dans le spectre de la tension Vf1-2 est due un dsquilibre de la commande des interrupteurs de londuleur.
133
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
(V)
50
f (Hz)
950
Fig. IV.5 : Spectre des mesures de la tension de londuleur (rsultats exprimentaux) Frquence (Hz) 100 250 350 550 650 Tension simple de londuleur (V) 1,15 0,98 0,9 2,88 0,9 Courant correspondant (Iinj) (A) 8,66 3 1,99 4,3 1,15
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
termes de poursuite, la fonction de transfert en boucle ferme peut tre donne par la relation suivante :
(T. Btot )(s) I inj (s) = (S. A + B .R)(s) I ref (s) tot tot 1 1 = avec g(s) s + 1 Btot ( s ) = 1 B1( s ) Atot ( s ) ts + 1 A( s )
(IV.3)
Pour un systme dorigine global, Btot(s)/Atot(s), dordre n=4 et en adoptant les mmes dmarches que dans les paragraphes III.3.2.2 et III.4.2 du chapitre II, on trouve les ordres des polynmes R(s), S(s) et T(s) :
deg(D)=2n=8, deg(R)=deg(S)=deg(T)=n=4.
Afin dvaluer lefficacit de la mthode de rgulation RSTam en termes de poursuite, le diagramme de Bode de la fonction de transfert (Iinj/Iref)(s) est donn sur la figure suivante :
Fig. IV.6 : Fonction de transfert (Iinj/Iref)(s) avec le rgulateur RSTam(s) De cette figure, on remarque que le dphasage maximal dans la bande passante est de 10 1250 Hz, ce qui tmoigne de la bonne poursuite de cette mthode de rgulation.
135
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
III.3.2 Etude des performanc es en compensation des perturbations III.3.2.1 Compensation des pertur bations de la tension du rseau A partir de la relation (IV.1), on peut exprimer linfluence, en boucle ferme, de la tension du rseau (tension naturelle) sur le courant inject Iinj par lquation suivante :
I inj (s) = (g.S. B2 )(s) Vs (s) (g.S. A + B1 .R)(s)
(IV.4)
Lquation (IV.5) montre leffet de lajout de la tension du rseau (tension compensatrice) sur le mme courant inject :
I inj (s) = (S. B1 )(s) Vs (s) (g.S. A + B1 .R)(s)
(IV.5)
A partir des deux expressions (IV.4) et (IV.5), nous avons pu tablir le diagramme de Bode qui est donn sur la Fig. IV.7. Ainsi, on peut observer linfluence de la tension du rseau (tension naturelle) et la tension ajoute (tension compensatrice) sur la boucle de rgulation.
Fig. IV.7 : Fonctions de transfert (Iinj/Vs)(s) dans les cas des tensions naturelle et compensatrice A partir de la Fig. IV.7, on remarque que la tension du rseau (tension naturelle), ne peut pas influencer la boucle de rgulation car elle est compense par la (tension compensatrice) (ayant, jusquau rang 17, les mmes gains et en opposition de phase). De plus, grce au rgulateur employ, ce type de perturbation est compens en boucle ferme, sur toute la bande de frquence, par un gain suffisamment ngatif en dB. III.3.2.2 Compensation des pertur bations de londuleur Quant leffet des perturbations de londuleur (Pond), elles sont modlises par la relation (IV.6) : 136
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
(g.S. B1 )(s) P (s) (g.S. A + B1 .R)(s) ond
I inj (s) =
(IV.6)
Fig. IV.8 : Fonction de transfert des perturbations de londuleur (Iinj/Pond)(s) A partir de la Fig. IV.8, on trouve que les perturbations de l'onduleur (Pond) nont aucune influence sur la boucle de rgulation. En effet, grce aux rgulateurs utiliss, les perturbations de rangs prpondrants, se situant dans la plage de frquence entre 100 Hz et 1200 Hz, sont rduites par des gains compris entre 17 dB et 33 dB. Cette compensation des perturbations est obtenue grce au rgulateur RSTam dont le gain est relativement lev sur la mme plage de frquence, comme le montre la Fig. IV.9. En effet, dans une boucle de rgulation, la rduction des perturbations est proportionnelle au gain du rgulateur.
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
Le gain dun rgulateur de type RSTam peut tre augment en diminuant les horizons Tc et Tf. Cependant, cette augmentation est limite par la robustesse du rgulateur et la saturation de la commande lors de limplantion numrique. Un compromis doit tre trouv entre ces diffrentes contraintes de manire rendre compte de la faisabilit de conception du filtre actif.
138
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
(THD)Is (THD)IL-23 (THD)IL (THD)Vs-av (%) (THD)Vs-ap (%) (%) (%) (%)
23 31 32,8 33.12 33,.12 31 30,5 0,9118 1,1 1,6 1,7 1,7 2,26 1,45 0,95 1,15 1,9 2,2 2,26 2,5 1,6 4,7 2 8 2,3 1,14 0,65 0,58 0,34 0,58 0,34 0,62 0,35 7,8 8
Groupe lectrogne :
Tab. IV.5 : Rsultats de simulation tablis partir du cahier des charges de la socit Socomec 139
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
Les notations de (THD)Is, (THD)IL et (THD)IL-23 reprsentent respectivement les taux de distorsion harmonique des courants de la charge, du rseau aprs filtrage et des 23 premiers harmoniques du courant du rseau aprs compensation. De mme, les notations de (THD)Vs-av et (THD)Vs-ap reprsentent respectivement les taux de distorsion harmonique des tensions de raccordement avant et aprs compensation. III.4.2 Synthse des rsultats En guise de synthse, nous avons repris point par point les exigences du cahier des charges en les comparant nos rsultats de recherche : Les rsultats de simulation montrent que le filtre actif parallle dpollue efficacement aussi bien les diffrents rseaux lectriques proposs que les groupes lectrognes dont la puissance varie de 100 kVA 2 MVA. La tension rseau peut varier de 6% avec deux frquences centrales (50 Hz et 60 Hz) sans consquences sur le fonctionnement du filtre actif. De mme que les frquences des groupes lectrognes varient jusqu 10% pour une variation de tension de 15%, comme lexige le cahier des charges. Le THD du courant aprs filtrage, pour les 23 premiers harmoniques, montre que les rsultats obtenus sont meilleurs que ceux requis par le cahier des charges avec un rapport qui varie entre 1/6 et 1/2 suivant les contraintes imposes. Ces rsultats confirment lefficacit de la mthode RSTam propose dans le cadre de la rgulation du courant du filtre actif parallle. Le THD de la tension au point de raccordement aprs la compensation est plus petit que celui obtenu avant compensation. Ce rsultat satisfait la contrainte exige de ne pas dgrader la tension dalimentation en compensant les courants harmoniques. Lefficacit de la PLL en termes de prcision et de rapidit a t prouve en appliquant aussi bien un changement de la tension du rseau, quun changement de la frquence du rseau.
Le temps de rponse de la poursuite de la tension continue est largement plus petit que le temps 100 ms demand. Cependant, celle-ci est accompagne dun petit dpassement transitoire en tension qui est de 7%. Aucun problme de rsonance nest observ entre le filtre de sortie et le rseau lectrique et/ou un filtre passif install. Le filtre actif peut fonctionner seul ou avec plusieurs modules en parallle et compenser des puissances plus leves, tout en restant loin des limites de THD requises. La mthode didentification des puissances instantanes, avec le systme base de PLL, sadapte aux changements de lamplitude et de la frquence du courant et de la tension du rseau.
Toutes ces performances de filtrage, obtenues par simulation, nous permettent daffirmer que le cahier des charges est largement satisfait, de dmontrer le caractre universel du filtre actif parallle dfini en prambule et de valider sa faisabilit exprimentale en mode continu.
140
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
Aucun risque datteindre les limites du cahier des charges ne semblait tre envisageable ce niveau du travail, mme en tenant compte des dgradations que pourrait engendrer une implantation numrique. Cependant, dans le but de prendre en compte toutes les contraintes dimplantation numrique, nous nous devons de dvelopper un modle discret des boucles de rgulation du filtre actif parallle.
141
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
lapproche, dite indirecte, qui consiste en une transposition du mode continu au mode discret des rgulateurs. Cette voie a t dveloppe partir de la commande C2d de discrtisation du logiciel Matlab, la synthse dite directe en z qui permettra dobtenir directement le modle discret des rgulateurs, dans le plan z.
(IV.7)
avec Te=80 ms la priode dchantillonnage correspondant une frquence dchantillonnage de Fe=12,5 kHz qui est la frquence de commutation de londuleur. Cette mthode de transformation est applique sur les deux rgulateurs :
Le rgulateur RST spcifique du systme base de PLL Le rgulateur RSTam de la rgulation du courant du filtre actif parallle.
IV.1.1 Discrtisation de la bo ucle de rgulation du systme base de PLL IV.1.1.1 Schma numrique
En se basant sur le schma de rgulation en continu du systme base de PLL (Cf. chapitre II, Fig. II.15), le schma quivalant en discret peut tre donn comme suit :
Vd-ref T(z)/R(z)
+-
R(z)/S(z)
Commande u
B(z)/A(z)
Vd
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
avec B(z)/A(z) reprsentant le systme dorigine qui consiste en un intgrateur, qui aura sous forme discrtise lexpression suivante :
Te z -1
(IV.8)
IV.1.1.2 Calcul du rgulateur RST (z) La commande C2d et la mthode Tustin sont appliques sur les fonctions de transfert T(s)/R(s) et R(s)/S(s) dj calcules en s (Cf. chapitre II, III.3.1). A partir de la Fig. IV.10, les diagrammes de Bode du rgulateur (R/S) et du systme en boucle ferme Vd-ref/Vd, la fois en mode continu et en mode discret, sont superposs sur les Fig. IV.11 et Fig. IV.12. De la Fig. IV.11, on observe que jusqu la frquence de 2 kHz les rgulateurs discrtiss et continus sont identiques, en gain et en phase. Les perturbations dues la tension perturbe du rseau lectrique peuvent donc tre rejetes. Un cart en gain et en phase commence apparatre en sapprochant de la frquence de Nyquist FN : FN=Fe/2. Cette frquence FN ajoute une contrainte supplmentaire par rapport aux contraintes de rgulation dj abordes dans le domaine continu. De la Fig. IV.12, on remarque que les gains en boucle ferme des systmes discrtis et continu sont bien superposs pour toutes les frquences y compris celle de Nyquist. Un cart de phase peut tre en revanche observ. Cet cart ninfluence pas le fonctionnement du systme base de PLL car les rfrences poursuivre sont constantes. Dans le cas du systme base de PLL, on constate que le passage du mode continu au mode discret seffectue de faon satisfaisante, ce qui tendrait montrer que cette approche de conception pourrait tre gnralise.
Discret
Discret
143
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
Discret
Discret
Fig. IV.12: Diagramme de Bode en boucle ferme Vd-ref/Vd en mode continu et discret IV.1.2 Discrtisation de la bo ucle de rgulation du courant IV.1.2.1 Schma numrique Le schma global du systme en boucle ferme de la rgulation du courant inject par le filtre actif, tel quil sera implant, est donn sur la Fig. IV.13. Le terme (1/z) de la boucle de rgulation reprsente le retard dune priode dchantillonnage d au temps de calcul par le processeur numrique. En effet, les calculs sont effectus lors dune premire priode dchantillonnage et ne peuvent tre utiliss qu la priode dchantillonnage suivante.
Vs
Gff
Pond
B2(z) Iinj
Iref
GVI
T(z)/R(z)
+-
R(z)/S(z)
+ +
1/z
GMLI
+ +
B1(z)
1/A(z)
GIinj
Fig. IV.13 : Schma dimplantation de la boucle de rgulation du courant Les gains, dans ce schma, reprsentent le taux de rduction des valeurs implantes par rapport celles relles du rseau lectrique, avec : 144
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________ - GIinj le gain de mesure du courant du filtre actif, - GVI le gain issu des mesures des courants et des tensions du rseau lectrique, - Gff le gain de la tension du rseau, ajoute pour empcher lappel dun courant fondamental au rseau lectrique par le filtre de sortie,
- GMLI = Vdc/3 le gain modlisant londuleur et reprsentant le coefficient K prsent dans la relation (II.30). A partir de la Fig. IV.13, la relation caractrisant la poursuite de la consigne est donne par lexpression suivante :
I inj (z) =
(T. B1 .GMLI .GVI )(z) (B2 .z.S)(z) I ref (z) + Vs (z) + (z.A.S + B1 .[Link] .GIinj )(z) (z.A.S + B1 .[Link] .GIinj )(z) + (B1 .S.G
MLI
.G )(z)
ff
Vs (z) +
(B .z.S)(z)
1
(IV.9)
Pond (z)
Les termes constituant la relation (IV.9) reprsentent la poursuite et la compensation des perturbations de la tension du rseau et de londuleur comme cela a t expliqu dans la relation (IV.1). A partir de la Fig. IV.13, on trouve que la fonction de transfert du modle dorigine (B1/A)(s), dfinie en s par les quations (II.6) et (IV.3), devient :
(IV.10)
A(z) = a z 3 + a z 2 + a z + a 1 2 3 4 2 B1 (s) = b1 z + b2 z + b3
(IV.11)
Le zro supplmentaire apparu dans la relation (IV.11), par rapport lquation (II.6), est d au retard de discrtisation du systme dorigine introduit par la commande C2d de Matlab et la mthode dchantillonnage avec bloqueur dordre zro (Zoh). Conformment lquation (IV.3), lquation (IV.11) montre que lordre du systme dorigine global est aussi n=4 en discret. Le retard dune priode dchantillonnage modlis par 1/z est d au fait que les courants de rfrence calculs, Iref, ne sont injects au rseau quaprs une priode dchantillonnage. IV.1.2.2 Calcul du rgulateur RST (z) La procdure utilise dans le cas du systme base de PLL peut aussi tre applique ici pour la conception du rgulateur de courant. En se basant sur le calcul des polynmes R, S et T effectu prcdemment dans ce chapitre dans le plan s (Cf. III.3.1), la transformation en z des fonctions de transfert (T/R et R/S) se fait laide de la commande C2d et de la mthode Tustin. Malheureusement, si cette mthode de discrtisation (C2d) peut garantir la stabilit dans le plan z dun rgulateur dj stable en mode continu, elle ne peut pas assurer la stabilit de la boucle ferme en termes de poursuite (Iinj/Iref)(z) [DELA96]. En effet, cette boucle ferme constitue par
145
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
ce rgulateur discrtis avec le modle dorigine global dfini par lquation IV.10 nest pas stable comme le montre la Fig. IV.14.
Fig. IV.14 : Ples et zros du systme en boucle ferme dans le plan z A partir de la Fig. IV.14, on remarque que trois ples en boucle ferme sont lextrieur du cercle unitaire, ce qui rend le systme en boucle ferme instable dans le plan (z). Les zros extrieurs au cercle unit, montrs dans cette figure, ne jouent aucun rle dans ltude de la stabilit. IV.1.2.3 Calcul dun rgulateur R ST(z) pour un systme boucl stable Le problme de cette technique de passage en numrique est la limitation impose par la frquence de Nyquist FN qui est gale la moiti de la frquence dchantillonnage Fe (FN= Fe/2). Cette limitation est trs restrictive par rapport ltude en mode continu o lavantage dune bande passante large permettait une meilleure application de la mthode de minimisation derreur. En effet, une bande passante limite par la frquence FN ne permet pas de relever suffisamment la courbe de phase sans influencer la poursuite en gain des rfrences. Il est rappeler quen mode continu, la frquence de coupure de la fonction de transfert, en terme de poursuite, est limite par la frquence de commutation qui est gale la frquence dchantillonnage dans notre cas. Ce passage serait possible en boucle ferme, en termes de stabilit, si la frquence dchantillonnage Fe tait au moins cinq fois suprieure la frquence de coupure en boucle ferme Fc (Fe= 5Fc). Afin de conserver les mmes avantages offerts par le rgulateur continu, il faudrait donc satisfaire la relation (Fe= 10Fc). Cela nous conduirait deux solutions :
diminuer la frquence de coupure moins de 2000 Hz, ce qui dgrade normment la compensation cause du dphasage, augmenter la frquence dchantillonnage et par consquent la frquence de commutation dun mme rapport.
La dernire solution prsente deux inconvnients, lun provenant de londuleur et lautre du DSP. En effet, laugmentation de la frquence de commutation dans un rapport de 5 10, pour une puissance de londuleur denviron 30 kVA, va augmenter significativement les pertes de londuleur 146
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
et par consquent diminuer son rendement. A lheure actuelle, cette augmentation de frquence nest pas ralisable techniquement. De plus, la diminution de la priode dchantillonnage est limite par la performance de calcul des DSP gnralement appliqus pour ce type dapplication. En effet, le modle contrle-commande du filtre actif requiert un minimum de temps de calcul qui ne peut tre rduit quavec une option plus performante de processeur. La frquence minimale dchantillonnage permettant, dans notre cas, le passage du mode continu au mode discret est de Fe =41 kHz. La Fig. IV.15 nous permet de comparer les diagrammes de Bode des systmes en boucle ferme dans le cas discrtis, avec Fe =41 kHz, et dans le cas continu, dj montr dans la Fig. IV.6.
Fig. IV.15 : Fonction de transfert (Iinj/Iref) avec les rgulateurs RSTam continu et discrtis A partir de la Fig. IV.15, on relve que le modle discret de la fonction de transfert du systme suit, jusqu la frquence de coupure, le mme gain et la mme phase que le modle continu. A la frquence de coupure, un dpassement important est observ dans le cas du modle discrtis. Ce dpassement est la consquence de la limitation de frquence de Nyquist (gale dans ce cas Fe/2=20.5 kHz). Il faut donc aller jusqu Fe=80 kHz pour obtenir un comportement similaire au rgulateur continu. Il est noter galement que le passage du mode continu au mode discret ne posait aucun problme dans le cas du systme base de PLL. En effet, dans ce cas et en mode continu, la frquence de coupure de la fonction de transfert en boucle ferme, montre dans la Fig. IV.12, tait 85 fois infrieure la frquence dchantillonnage. IV.1.3 Synthse de la discrti sation par la commande C2d La premire mthode de passage du mode continu au mode discret effectue par la commande C2d de Matlab montre que cette mthode est applicable si la frquence de coupure de la fonction de transfert en continu est trs loigne de la frquence dchantillonnage. Cela est ralisable dans le cas de la poursuite des signaux continus (cas de la PLL) ou des signaux de frquence trs faible (cas de la compensation des creux de tension 50 Hz par le filtre actif srie).
147
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
Si la rfrence poursuivre contient des harmoniques, cette mthode de discrtisation nest pas utilisable dans les limites technologiques actuelles. Ceci nous conduit utiliser une autre mthode de passage en mode discret, il sagit de concevoir le rgulateur directement dans le plan z.
148
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
deg(D)=2n=8,
Cela implique que les degrs des polynmes R(z) et S(z) soient gaux 4 :
deg(R)=deg(S)= n=4
Le principe de la mthode de minimisation derreur, dj dcrite dans ce mmoire (Cf. chapitre II, III.4.2), reste le mme. Par contre, le calcul de T(z) devient plus compliqu effectuer. Ce calcul, conformment ce que nous avons obtenu dans le plan s, nous donne un polynme T(z) dordre 4 :
deg(T)=n=4
Finalement, les polynmes R(z), S(z) et T(z) ont directement pu tre calculs dans le plan z, nous permettant ainsi dliminer les problmes dinstabilit observs en boucle ferme lors de la discrtisation par la commande C2d de Matlab. IV.2.2 Etudes des performanc es du rgulateur RSTam(z) en boucle ferme Dans cette partie, nous allons tudier les performances du rgulateur RST(z), en boucle ferme, en termes de poursuite et en termes de compensation des perturbations. IV.2.2.1 Etude des performances e n terme de poursuite A partir de la Fig. IV.13 et de la relation (IV.9), la fonction de transfert en boucle ferme entre le courant de rfrence et le courant inject peut tre donne par la relation suivante :
I inj (z) =
(T. B1 .GMLI .GVI )(z) I ref (z) (z.A.S + B1 .[Link] .GIinj )(z)
(IV.12)
149
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
La Fig. IV.16 reprsente le diagramme de Bode du systme en boucle ferme, tabli partir de la relation (IV.12). Cette dernire est constitue par le rgulateur RSTam, labor dans le plan z partir du modle dorigine global dcrit dans lquation (IV.10).
Fig. IV.16 : Fonction de transfert (Iinj/ Iref)(z) avec le rgulateur RSTam calcul en z La Fig. IV.16 montre que nous pouvons compenser, de faon quasi parfaite, les harmoniques 5, 7, 11 et 13. Les harmoniques de rang plus lev, appartenant la bande passante, sont bien poursuivis en phase, grce la minimisation derreur, mais ils sont par contre amplifis en gain. Cette amplification est due la limitation de frquence de Nyquist qui est de 6,25 kHz dans notre cas. IV.2.2.2 Etude des performances e n compensation des perturbations Dans cette partie, nous allons tudier les perturbations dues la tension du rseau et londuleur ainsi que leurs effets sur la qualit de rgulation en boucle ferme. IV.[Link] Effets des perturbations de la tension du rseau A partir de la relation (IV.9), nous pouvons exprimer linfluence de la tension du rseau Vs (tension naturelle) sur le courant inject Iinj par lquation suivante :
I inj (z) =
(IV.13)
Lquation IV.14, ci-dessous, montre leffet de lajout de la tension du rseau (tension compensatrice) sur le mme courant inject :
I inj (z) =
(IV.14)
A partir des deux expressions IV.13 et IV.14, nous avons pu tablir le diagramme de Bode donn sur la Fig. IV.17. Ainsi, on peut observer linfluence de la tension du rseau (tension naturelle) et celle de la tension ajoute (tension compensatrice) sur la boucle de rgulation. 150
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
Fig. IV.17 : Fonctions de transfert (Iinj/Vs)(z) naturelle et compensatrice De la Fig. IV.17, on remarque que la tension du rseau (tension naturelle), ne peut pas influencer srieusement la boucle de rgulation car elle est compense par la tension compensatrice pour les tensions harmoniques conventionnelles. En effet, cette dernire a les mmes gains et se trouve en opposition de phase avec la tension naturelle du rseau. Sils existent, les harmoniques plus levs de cette tension seront compenss en boucle ferme par un gain ngatif en dB. IV.[Link] Effets des perturbations de londuleur Lapparition des perturbations engendres par londuleur nous a amen modliser leurs effets en boucle ferme par la relation (IV.15) en nous basant sur la Fig. IV.13 et la relation (IV.9).
I inj (z) =
(IV.15)
Afin dtudier leffet de ce type de perturbation, en boucle ferme et sur la sortie (Iinj), le gain du diagramme de Bode, reprsentant la relation (IV.15), est donn par la Fig. IV.18. Le gain du diagramme de Bode de la Fig. IV.18 montre linfluence des perturbations engendres par londuleur Pond sur la boucle de rgulation. A partir de cette figure, on trouve que les perturbations Pond sont lgrement amplifies en boucle ferme ce qui augmente la dgradation de rgulation dj releve dans la boucle rgulation en termes de poursuite de la Fig. IV.16.
151
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
152
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
Cette diffrence de performance peut sexpliquer par la diffrence de gain des rgulateurs lors dune conception en mode continu ou en mode discret. En effet, les frquences des ples du systme boucl en mode continu sont limites par la frquence de commutation (Fc=12,5 kHz), alors que les frquences des ples en mode discret se limitent la frquence de Nyquist FN qui est la moiti de la frquence de commutation (FN=Fc/2=Fe/2=6,25 kHz). Cette limitation en frquence limitera le gain du rgulateur et par consquent son influence sur les perturbations comme le montre la figure suivante.
Fig. IV.19 : Fonction de transfert du rgulateur conu en discret A partir de la Fig. IV.19, nous remarquons que le rgulateur prsente un gain trs faible pour les frquences correspondant aux perturbations prpondrantes cites sur le Tab. IV.4. Une limitation en frquence plus leve (12,5 kHz) a permis en mode continu, en utilisant la mme
153
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
mthode de rgulation, davoir un gain beaucoup plus important, compensant toutes les perturbations qui risquent de dgrader la qualit de filtrage (Cf. Fig. IV.9). Suite ces analyses, nous constatons que les performances du rgulateur RSTam, proposes en mode continu pour rpondre des contraintes nouvelles de compensation, subissent dans le domaine discret des dgradations dues principalement la limitation en frquence. En effet, le dpassement remarqu en Fig. IV.16 sur la courbe du gain en poursuite ainsi que la lgre amplification des perturbations observe sur les Fig. IV.18, sont dus la limitation de frquence de Nyquist. Afin de rsoudre ce problme, nous proposons deux solutions :
ajouter des contraintes supplmentaires sur le polynme S(z) pour faire tendre les quations IV.13, IV.14 et IV.15 vers zro pour des harmoniques spcifiques, augmenter raisonnablement la frquence de commutation de londuleur pour rsoudre les problmes de dpassement et damplification des perturbations.
La premire solution, dont le principe a dj t utilis dans le cas du systme base de PLL (Cf. chapitre II, III.3.1), implique lajout de filtres frquences prcises. Cela peut rsoudre les problmes engendrs par les perturbations (Vs et Pond) tout en augmentant en mme temps lordre global de la fonction de transfert en termes de poursuite. Cette augmentation de lordre du systme occasionne une amplification du problme de dphasage qui peut se rsoudre par laugmentation de lordre du polynme T(z) (Cf. chapitre II, III.4). Malheureusement, cette solution va contribuer laugmentation du dpassement sur la mme plage de frquence o les harmoniques doivent tre compenss (Cf. Fig. IV.16). La deuxime solution, qui consisterait repousser les limites de commutation de londuleur, malgr les contraintes technologiques, peut constituer une solution srieuse aux problmes numriques rencontrs dans notre tude. Il nous a donc apparu intressant dtudier la frquence minimale de commutation de londuleur permettant la discrtisation de notre rgulateur.
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
frquence 1500 Hz, le gain commence augmenter entranant le type de dpassement dj observ. Ce dpassement peut influencer les harmoniques suprieurs 1500 Hz qui sont supposs tre ngligeables.
Fig. IV.20 : Fonction de transfert (Iinj/ Iref)(z) avec le rgulateur RSTam calcul en z pour Te=45 ms V.2.2 Etude des performanc es du rgulateur en compensation des perturbations Linfluence en boucle ferme des mmes types de perturbation, causs par la tension du rseau Vs et par londuleur Pond, est tudie dans cette partie du chapitre. V.2.2.1 compensation des perturb ations causes par la tension du rseau A partir des deux expressions (IV.13) et (IV.14), nous avons pu tablir le diagramme de Bode qui est donn sur la Fig. IV.21.
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
Ainsi, on peut observer linfluence de la tension du rseau (tension naturelle) et la tension ajoute (tension compensatrice) sur la boucle de rgulation. A partir de la Fig. IV.21, on remarque que la tension naturelle du rseau, ne peut pas influencer la boucle de rgulation. En effet, cette dernire est compense dune part par la tension compensatrice pour les tensions harmoniques conventionnelles (ayant les mmes gains et tant en opposition de phase), dautre part par un gain ngatif en dB. V.2.2.2 compensation des perturb ations causes par londuleur Afin dtudier leffet de ce type de perturbation en boucle ferme et sur la sortie (Iinj), le gain du diagramme de Bode, reprsentant la relation (IV.15), est donn par la Fig. IV.22.
Fig. IV.22 : Fonction de transfert des perturbations de londuleur(Iinj/Pond)(z) A partir de la Fig. IV.22, on trouve que les perturbations prpondrantes causes par l'onduleur (Pond) sont compenses en boucle ferme par un gain de 9 dB.
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
noccasionne pas damplification des perturbations dues londuleur mais ne les rduit pas non plus.
157
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
VI. CONCLUSION
Dans ce chapitre, nous avons dfini et abord la faisabilit exprimentale dun filtre actif parallle universel partir des contraintes dun cahier des charges prcis fourni par une socit partenaire (la socit franaise Socomec) . Ce filtre actif devait sadapter aux changements du rseau lectrique (ou groupe lectrogne) en termes de puissance, de frquence, et de tension. Afin de satisfaire le cahier des charges, le systme base de PLL a t inclus dans la mthode didentification des harmoniques. De plus, un filtre de sortie du troisime ordre (filtre en T) a t employ et le principe de la nouvelle mthode de rgulation du courant du filtre actif RSTam a t appliqu. Des tudes de faisabilit exprimentale de ce filtre actif en mode continu et discret ont t effectues. Les perturbations qui peuvent dgrader la qualit de compensation du filtre actif ont t identifies travers des mesures effectues sur un prototype du filtre actif dvelopp par la socit Socomec. Dans un premier temps, lensemble du filtre actif parallle a t tudi en mode continu, en prenant en compte, dans les modles de dfinition, des perturbations et des contraintes matrielles inhrentes la faisabilit exprimentale du filtre actif dvelopp par notre partenaire de recherche. La mthode de rgulation RSTam propose a dmontr sa robustesse en termes de stabilit et de performance. La prcision du systme base de PLL, face aux changements des rfrences et aux perturbations en amplitude et en frquence, a t galement examine travers des simulations. La rapidit de rponse du systme retrouver la tension continue de rfrence, aprs une dconnexion brutale de la charge, a t dmontre. Aucune rsonance entre le filtre de sortie, le rseau lectrique et un ventuel filtre passif install sur le rseau, nest apparu lors des simulations effectues. Lensemble de ces rsultats ont permis de valider en mode continu lapproche de faisabilit exprimentale du filtre actif parallle universel. Cependant, cette approche de faisabilit exprimentale ne pouvait tre complte sans une tude en mode discret du rgulateur de courant appliqu au filtre actif parallle universel. En effet, en boucle ferme, nous avons t contraint par la frquence de Nyquist, qui est la moiti de la frquence de commutation de londuleur, ce qui nous a oblig repenser notre modlisation de rgulation. Nous avons alors propos deux approches de conception en mode discret. La premire approche sapplique lorsque les contraintes du rgulateur se limitent poursuivre une rfrence unique et constante, tout en rejetant, lintrieur de la boucle de rgulation, quelques harmoniques. Ainsi, nous avons conu le rgulateur en mode continu, puis nous avons effectu la transformation du mode continu au mode discret partir de la seule commande de discrtisation C2d du logiciel de calcul Matlab. Ce passage de mode a pu tre ralis, pour la mme frquence de commutation prvue (12,5 kHz), et sans perte dinformation, pour le rgulateur RST spcifique au systme base de PLL. Cependant, dans le cas de la rgulation du courant du filtre actif, les rfrences poursuivre pouvant aller jusqu 2000 Hz, la mthode du passage du mode continu au mode discret ne garantit plus la stabilit en boucle ferme pour la frquence de commutation prvue. Afin de rsoudre ce problme, deux rgulateurs, bass sur les mmes principes de conception du rgulateur RSTam en mode continu ont t conus directement dans le plan z. Le premier a t calcul pour la frquence dchantillonnage requise dans le cahier des charges (12,5 kHz). Cependant, ce rgulateur na permis une poursuite quasi parfaite, en gain et en phase, des harmoniques de rfrence que jusquau rang 13, tout en amplifiant lgrement les perturbations causes par londuleur. 158
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
Nous avons remarqu travers toutes les tudes effectues en mode continu et en mode discret que les perturbations dues la tension, dite naturelle, du rseau sont compenses dans tous les cas par lajout de la tension, dite compensatrice, du rseau. Le deuxime rgulateur a t conu en discret pour une frquence dchantillonnage de 22 kHz. Ce rgulateur permet de poursuivre parfaitement en gain et en phase les harmoniques de rfrence jusquau rang 25 tout en rduisant raisonnablement les perturbations de londuleur. Laugmentation de la frquence dchantillonnage et par consquent la frquence de commutation jusqu 22 kHz limite, actuellement et cause des problmes technologiques, linstallation de ce filtre actif universel pour dpolluer des charges de puissances trs leves.
159
Chapitre IV : Faisabilit exprimentale dun filtre actif active parallle universel ___________________________________________________________________________
160
161
162
163
On peut aussi envisager de travailler spcifiquement sur la rgulation du systme PLL, comme par exemple : dintroduire lamplitude de tension directe (Vd) du systme PLL dans la boucle de rgulation pour viter le retard occasionn par le filtre passe-bas employ. dessayer dliminer toutes les perturbations, dans la boucle de rgulation de la PLL, au lieu den liminer deux aux frquences 100 Hz et 300 Hz. Cela va amliorer la prcision de la mthode didentification de tensions base sur un systme PLL. Cette amlioration peut tre envisage par lapplication dun rgulateur H.
En lectronique et microlectronique, nous pensons videmment la puissance de traitement de certains processeurs qui pourrait nous permettre datteindre nos objectifs de vitesse dchantillonnage, de manire repousser les limites thoriques et pratiques de certains rgulateurs. Dans ce sens, des tudes doivent tre envisages en microlctronique, en particulier dans lutilisation de circuits intgrs spcifiques, tels que les ASIC (Application Specific Integrated Circuit) et PLD (Programmable Logic Device), pour rpondre aux contraintes de calcul actuelles. Des recherches peuvent tre ainsi dveloppes dans le domaine de ladquation algorithme architecture en vue de limplantation de structures de commande plus complexes.
165
166
REFERENCES
167
168
Rfrences ___________________________________________________________________________
[AK83] H. Akagi, Y. Kanazawa and A. Nabae, Generalized theory of the instantaneous reactive power in three-phase circuits , Proceeding 1983 International power electronics conference. Tokyo, Japan, PP. 1375-1386, 1983. [AK86] H. Akagi, A. Nabae and S. Atoh, Control strategy of active power filters using multiple voltage-source PWM converters , IEEE Trans. on Industry applications, vol. IA-22, pp. 460-465, 1986. [AK90] H. Akagi, Y. Tsukamuto, A. Nabae, Analysis and design of an active power filter using quad-series voltage source PWM converters , IEEE Trans. on Industry applications, vol. 26, No. 1, pp. 93-98, 1990. [AK96] H. Akagi, New trends in active filters for power conditioning , IEEE Trans. on Industry applications, vol. 32, No. 6, pp. 1312-1322, November/December 1996. [AK97] H. Akagi, Control strategy and site selection of a shunt active filter for damping of harmonic propagation in power distribution systems , IEEE Trans. on power delivery, vol. 12, No. 1, pp. 354-363, January 1997. [AL2000] G. Alarcon, C. Nunez, V. Cardenas, M. Olivier, Design and implementation of a 3-phase series active filter to compensate voltage disturbance CIEP-2000, October 2000, Acapulco, Mexico, P. 93-98. [ALA98] M.A.E. Alali, A.S. Snyder and al, A robust damping controller for power systems using linear matrix inequalities , IEEE power engineering society, winter meeting, New York 1998. vol . 1, PP. 519-524, 1998. [ALA2000-1] M.A.E. Alali, S. Saadate, Y.A. Chapuis, F. Braun, Energetic study of a shunt active conditioner compensating current harmonics, power factor and unbalanced , EPE-PEMC 2000, Kosic, Slovak Republic, vol. 5, pp. 211-216, September 2000. [ALA2000-2] M.A.E. Alali, S. Saadate, Y.A. Chapuis, F. Braun, Energetic study of a series active conditioner compensating voltage dips, unbalanced voltage and voltage harmonics , IEEE-CIEP-2000 Acapulco, Mexico, P. 80-86, Octobre 2000. [ALA2000-3] M.A.E. Alali, S. Saadate, Y.A. Chapuis, F. Braun, Control and analysis of series and shunt active filters with SABER , IPEC 2000, Tokyo, Japan, PP. 1467-1472, April 2000. [ALA2001] M.A.E. Alali, S. Saadate, [Link], [Link], Stratgies de dpollution active des rseaux lectriques, parti I : Commande gnralise et analyse de performances dun compensateur actif parallle , revue internationale de gnie lectrique, vol : 4, -n3-4/2001, PP. 293-316, 2001. [ALA2001-1] M.A.E. Alali, S. Saadate, [Link], [Link], Stratgies de dpollution active des rseaux lectriques, parti II : Compensateur actif srie dharmoniques, de dsquilibre et de creux de tension des rseaux lectriques , revue internationale de gnie lectrique, vol : 4, -n3-4/2001, PP. 317-332, 2001.
169
Rfrences ___________________________________________________________________________
[ALA2001-2] M.A.E. Alali, S. Saadate, Y.A. Chapuis, F. Braun, Advanced corrector with FPGA-based PLL to improve performance of series active filter compensating all voltage disturbances , EPE-2001 Graz, Austria, August 2001. [ALA2001-3] Mohamad .A.E. Alali, Optimisation des performances de filtres actifs parallles-sries partir des correcteurs avancs , JCGE01, Belfort/Nancy, pp. 215-222, 2001. [AR95] M. Aredes, J. Hfner, K. Heumann, A combained series and shunt active power filter , IEEE/KTH Stockholm power tech conference, Stockholm, Sweden, June 18-22, 1995. [BENCH98] L. Benchaita, Etude, par simulation numrique et exprimentation, dun filtre actif parallle structure courant avec une nouvelle mthode de contrle-commande Thse de lInstitut National Polytechnique de Lorraine, Nancy 30 Octobre 1998. [BEN98] H. Benquassimi, Conception de convertisseurs alternatifs continu mono-tage absorption sinusodale Thse de lInstitut National Polytechnique de Grenoble, 18 Juin 1998. [BH94] V. B. Bhavaraju, P. Enjeti, A fast active power filter to correct line voltage sags , IEEE Trans. on Industrial Electronics, vol. 41, NO. 3, P. 333-338, June 1994. [BI69] B. Bird, J.F Marsh, P.L McLellan, Harmonic reduction in multiple converter by triple frequency current injection , Inst. Elect. Eng. Proc, vol. 116, No. 10, pp. 1730-1734, April 1969. [BO99] M.H.J. Bollen, Understanding power quality problems: voltage sags and interruptions , IEEE Trans. on Industrial Electronics, vol. 41, NO. 3, P. 333-338, June 1994. [BR99] N. Bruyant, Etude et commande gnralises de filtres actifs parallles : compensation globale ou slective des harmoniques rgime quilibr ou dsquilibr Thse de doctorat de luniversit de Nantes, Saint-Nazaire, 25 Novembre 1998. [CA92] A. Campos, G. Joos, P. Ziogas, J. Lindsay, Analysis and design of a series voltage compensator for three-phase unbalanced sources , IEEE Trans. on Industrial Electronics, vol.39, N2, P. 159-167, April 1992. [CA94] A. Campos, G. Joos, P. Ziogas, J. Lindsay, Analysis and design of a series voltage unbalance compensator based on a three-phase VSI operating with unbalanced switching function , IEEE Trans. on power electronics, vol.9, N3, P. 269-274, May 1994. [CH2000] A. Chandra, B. Zingh, B.N. Zingh, K. Al-Haddad, An improved control algorithm of shunt active filter for voltage regulation, harmonic elimination, power factor correction and balancing of nonlinear loads , IEEE Trans. on power electronics, vol.15, No. 3, pp. 495-507, May 2000.
170
Rfrences ___________________________________________________________________________
[DES95] G. Desquilbet, C. Foucher, P. Fauquembergue, Statistical analysis of voltage dips Notes EDF, 96NR00102, 1995. [DE98] T. Deflandre, P. Mauras, Les harmoniques sue les rseaux lectriques Edition Eyrolles, 1998. [DeCH97] F. De Chateauvieux, Creux de tension et coupures brves Notes EDF, Groupe de travail UIE power quality: Guide sur la qualit de llectricit pour lalimentation des installations industrielle, partie 2, HR22/96/045.A, 1996. [DELA96] P. De larmina, Automatique, Commande des Systmes Linaires , 2nd Edition, HERMES, 1996. [DI97] J. W. Dixon, G. Venegas, L. Moran, A series active power filter based on sinusoidal current-controlled voltage-source inverter , IEEE Transaction on Industrial Electronics, vol. 44, NO. 5, PP. 612-620, October 1997. [FU96] H. Fujita, H. Akagi, The unified power quality conditioner: the integration of series- and shunt-active filters , IEEE Trans. on Power Electronics, vol. 13, No. 2, pp. 1312-1322, March 1998. [GY76] L. Gyugyi and E.C. Strycula, Active AC power filters , IEEE/IAS Annual Meeting, pp. 529-535. 1976. [GR77] M. Grandpierre, B. Trannoy, A static power device to rebalance and compensate reactive power in three phase network: design and control , IEEE Ind. Appl. Soc. Annual meeting, pp. 127-135, 1977. [GR90] W. Grady, M. Samotyj, and A. Noyola, Survey of active power line conditioning methodologies , IEEE Trans. on Power Delivery, vol. 5, pp. 1536-1542, 1990. [GU2000] S. Guffon, Modlisation et commandes a structure variable de filtres actifs de puissance Thse de lInstitut National Polytechnique de Grenoble, 24 Juillet 2000. [HAD97] K. Haddad, G. Joos, Distribution system voltage regulation under faults conditions using static series regulators , IEEE Industry Applications Society, Annual meeting, New Orleans, PP. 1383-1389, October 5-9, 1997. [KA83] H. Kawahira, T. Nakamura, S. Nakazawa, Active power filtres , JIEE IPEC, Tokyo, Japon, pp. 981, 1983. [KAU97] V. Kaura, V. Blasko, Operation of a phase locked loop system under distorted utility conditions , IEEE Trans. on Industry Applications, vol. 33, No. 1, pp. 58-63, January/February 1997. [LA99] O.L. Lapena, L.G. de Vicuna, J. Matas, M. Lopez and M. Castilla, A sliding mode control design of an active power filter with a low-pass filter as input stage , EPE-99 Lausanne, Suisse, 1999.
171
Rfrences ___________________________________________________________________________
[LO95] C. Lott, Filtrage actif parallle des harmoniques de courant sur des rseaux industriels : tude thorique et ralisation dune maquette GTO Thse de doctorat de lInstitut National Polytechnique de Lorraine, Nancy 17 Octobre 1995. [ME90] P. Meynaud, Qualit de la tension dans les rseaux lectriques creux de tension, flicker et harmoniques Techniques de lingnieur, trait Gnie lectrique, D 4 260, pp. 1-25, 1990. [ME96] P. Meynaud, P. Pene, G. desquilbet, Dsquilibre de tension Notes EDF, 96NR00003, 1996. [MO77] [Link], and al, Active filters for ac harmonic suppression , IEEE/PES winter meeting, A77 026-8, 1977. [MOR90] L. Morn, P. D. Ziagas, G. Joos, Three-phase solid-state voltage compensator system , Can. J. Elect. Comp. Eng., vol. 15, no. 1, PP. 27-35, 1990. [MOR99] L.A. Morn, I. Pastorini, J. Dixon,R. Wallace, A fault protection scheme for series active power filters , IEEE Trans. on power electronics, vol.14, N5, P. 928-938, September 1999. [MU85] M. Muljadi, R. Shiferl, T. A. Lipo, Induction machine phase balancing by unsymmetrical thyristor voltage control , IEEE Trans. on Industry applications, vol. IA-21, no. 4, PP. 669-678, May/June 1985. [NAK89] T. Nakajima, E. Masada, An active power filter with monitoring of harmonic spectrum , EPE-89, 3rd European conference on power electronics and applications, Aachen, Germany, 1989. [PE90] F.Z. Peng, H. Akagi and A. Nabae, A study of active power filters using quad-series voltage-source PWM converters for harmonic compensation , IEEE Trans. on power electronics, vol. 5, N1, pp. 9-15, January 1990. [SA71] [Link], [Link], A new method to eliminate ac harmonic currents by magnetic compensation-consideration on basic design , IEEE Trans. Power Appl. Syst, , vol. 90, No. 5, pp. 2009-2019, 1971. [SAL96] A. Salem Nia, Contribution ltude thorique et exprimental dun filtre actif parallle commande analogique et numrique temps rel Thse de lInstitut National Polytechnique de Lorraine, Nancy 24 Septembre 1996. [SAN2000] A. Sannino, J. Svensson, A series connected voltage source converter for voltage sag mitigation using vector control and a filter compensation algorithm , 35th IAS annual meeting of the 2000 IEEE Industry Applications Conference, 8-12October 2000 Rome, Italy, P. 2476-2481. [SZ84] B. Szabados, Rpercution sur le rseau lectrique de lalimentation des convertisseurs de grande puissance REG, N. 12, pp. 791-796, 1984.
172
Rfrences ___________________________________________________________________________
[TO99] R. Tounsi, Dveloppement dun contrle commande pour un compensateur srie de creux de tension. Validation par simulation du fonctionnement avec des charges industrilles Thse de lInstitut National Polytechnique de Toulouse, 30 Octobre 1999. [VI2000] A.J. Visser, T. Mouton, J.H.R. Enslin, Directed-coupled cascaded multilevel sag compensator IEEE 31st annual power electronics specialists conference, 18-23 June 2000 Galway, Ireland. [WA92] Miao-Xin Wang, Filtrage actif de puissance : Etudes et ralisation dun filtre actif commande numrique temps rel Thse de lInstitut National Polytechnique de Toulouse, 18 Dcembre 1992. [XU94] J. Xu, Filtrage actif parallle des harmoniques des rseaux de distribution dlectricit Thse de lInstitut National Polytechnique de Lorraine, Nancy 20 Janvier 1994.
173
Rfrences ___________________________________________________________________________
174
ANNEXES
175
176
Annexes ___________________________________________________________________________
(A-I.1)
Nous remarquons que les courants harmoniques sont de rang (6K 1) (K =1, 2, 3, ....etc.) et que la valeur efficace de chaque courant harmonique Ih est inversement proportionnelle au rang correspondant (Ih = I1/h). Dans cette annexe, les cas de la tension du rseau et du courant de la charge polluante les plus probables parvenir, vont tre analyss et les pulsations des puissances instantanes pour chaque cas vont tre calcules.
(A-I.2)
Le courant de la charge polluante est celui reprsent par lquation A-I.1. En faisant la transformation a-b, nous obtenons les courants et les tensions suivants :
V a sin w t = 3Vs - cos w t V b sin w t 1 sin 5( w t - a ) 1 sin 7( w t - a ) Ia - - + .... = 3 I 1 I b cos w t 5 cos 5( w t - a ) 7 - cos 7( w t - a )
(A-I.3)
(A-I.4)
177
Annexes ___________________________________________________________________________
Dans ce calcul, nous allons analyser la puissance relle p(t) puis gnraliser les rsultats pour la puissance imaginaire q(t). La puissance relle peut tre calcule partir des quations (A-I.3) et (A-I.4) comme le montre la relation (A-I.5):
P( t ) = 3Vs I 1 cos a + 14 4 2 3
P
(A-I.5)
Le premier terme (continu) de cette relation reprsente la puissance active conventionnelle, tandis que les autres termes reprsentent la puissance alternative cause par les courants perturbateurs (harmonique dans ce cas). A partir de cette relation, nous pouvons tracer la transformation du spectre de la puissance active instantane :
Composantes harmoniques directes
1 5
7 11
12
13 17
18
19
Rang harmonique
Fig. A-I.1 : Transformation du spectre de la puissance relle instantane (courant harmonique et tension sinusodale quilibre) De la Fig. A-I.1 nous constatons que la puissance conventionnelle active devient une composante continue, alors que les courants harmoniques (5,7), (11,13), (17,19), causent des puissances perturbatrices aux pulsations respectives (6w, 12w, 18w). Le mme spectre peut tre obtenu pour la puissance imaginaire instantane.
(A-I.6)
avec Vd, Vinv reprsentent respectivement les composantes directes et inverses de tension. Le courant est celui reprsent par lquation (A-I.1). En effectuant la transformation a-b, nous obtenons les courants, donns par la relation (A-I.4), et les tensions suivants :
(A-I.7)
178
Annexes ___________________________________________________________________________
Dans ce cas, nous allons aussi analyser la puissance relle p(t) et puis gnraliser les rsultats pour la puissance imaginaire q(t). La puissance relle peut tre calcule partir des quations (II.15), (A-I.4) et (A-I.7) comme le montre la relation A-I.8:
3V I 3V I 3Vinv I 1 cos( 2wt - a ) + d 1 cos( 6wt - 5a ) - inv 1 cos( 4wt - 5a ) 5 5 P( t ) = 3Vd I 1 cos a - (A-I.8) 3Vinv I 1 3Vd I 1 14 4 2 3 cos( 8wt - 7a ) + ..... cos( 6wt - 7a ) + P 7 4 7 1444444444444 2444444444444 3 4 ~ P
Le premier terme (continu) de cette relation reprsente la puissance active conventionnelle, tandis que les autres termes reprsentent la puissance alternative cause par les courants perturbateurs harmoniques et les tensions directes et inverses du rseau lectrique. A partir de cette relation, nous pouvons tracer la transformation du spectre de la puissance active instantane reprsente par la Fig. II.6 (courant harmonique et tension directe) et celle suivante :
Composantes harmoniques inverses 0 1 2 5 7 8 11 13 14 17 16 19 20 h Rang harmonique
10
Fig. A-I.2 : Transformation du spectre de la puissance relle instantane (courant harmonique et tension sinusodale inverse) De la relation A-I.8, nous constatons que la partie directe de la tension du rseau lectrique avec les courants harmoniques (fondamental inclus) ont le mme effet prsent dans le premier cas analys. Par contre, nous pouvons constater de la Fig. A-I.2 que la partie inverse de la tension avec le courant fondamental crent une puissance perturbatrice la pulsation 2w. Alors que cette mme tension inverse avec les courants harmoniques inverses de rangs (5,11,17), et ceux directs de rangs (7,13,19.), causent respectivement des puissances perturbatrices aux pulsations (4w,10w,16w) et (8w,14w,20w). Le mme spectre peut tre obtenu pour la puissance imaginaire instantane. Pour le cas inverse, dune tension quilibre avec un courant harmonique dsquilibr, les mmes spectres prsents dans les figures A-I.1 et A-I.2 peuvent tre obtenus. De la mme manire, dans le cas o la tension du rseau est dsquilibre contenant des harmoniques, avec un courant harmonique dsquilibr, nous pouvons observer que les composantes perturbatrices de courant et de tension du mme rang occasionnent des composantes continues qui peuvent sajouter celles reprsentant les puissances conventionnelles active et ractive. De ces analyses, il parat comme vident que la tension du rseau lectrique ne doit contenir que la composante directe la frquence fondamentale. Sinon, les composantes perturbatrices de tension (harmoniques et dsquilibres) introduisent des composantes non dsires des puissances instantanes. Ce qui peut avoir pour consquence de donner une valeur errone des courants perturbateurs identifis.
179
Annexes ___________________________________________________________________________
ANNEXE II
Dans cette annexe, les analyses temporelle et spectrale du filtre actif srie compensant les creux de tension biphas et triphas seront montres. Dans chaque cas simul, le filtre actif srie emploie la mthode didentification des perturbations de tension, base sur le systme base de PLL, et une des mthodes de rgulation de tension (PID, RST et RSTam).
(V)
Vref2 Vinj2
(V)
Vref1 Vinj1
Vch1,2,3 (V)
t (s)
Fig. A-II.1: Analyse temporelle de la compensation du creux de tension biphas (PID) (Lgende - phase 1 : ., phase 2 : - - - - -, phase 3 : -..-..-..-) Cette figure montre la rapidit des mthodes didentification, base dun systme base de PLL, et de la rgulation PID face un changement brutal de la consigne. En effet, suivant le changement du point de fonctionnement du rseau lectrique, les tensions identifies (Vref1, Vref2) et celles injectes (Vinj1, Vinj2) rpondent rapidement ces variations brutales. Cette rapidit de compensation du creux de tension permet la circulation dun courant quilibr malgr ces changements du point de fonctionnement du rseau lectrique.
180
Annexes ___________________________________________________________________________
Afin dvaluer la prcision des mthodes didentification et de rgulation employes, une analyse spectrale des Ph1 et ph2 a t effectue sur la Fig. A-II.2.
Vs1 (V)
Vref1 (V)
Vinj1 (V)
Vch1 (V)
f (Hz)
Fig. A-II.2 : Analyse spectrale de la compensation dun creux de tension biphas (PID) Cette figure o les tensions identifie (Vref1) et injecte (Vinj1) sont presque gales, permettant la tension de la charge (Vch1) atteindre la valeur crte (Vd-des=325 V), tmoigne de la bonne prcision des mthodes didentification et de rgulation tudies.
Remarque :
Les petites oscillations remarques sur les tensions injectes sur les Phase 2 et Phase 3 sont dues, entre autres, des problmes de calcul d au simulateur, car le dbut et la fin du creux de tension des Phase 2 et Phase 3 n'arrivent pas au moment du passage par zro, comme c'est le cas de la Phase 1.
181
Annexes ___________________________________________________________________________
Vs1,2,3 (V)
(V)
Vref2 Vinj2
(V)
Vref1 Vinj1
Vch1,2,3 (V)
t (s)
Fig. A-II.3: Analyse temporelle de la compensation du creux de tension biphas (RST) (Lgende - phase 1 : ., phase 2 : - - - - -, phase 3 : -..-..-..-)
Vs1 (V)
Vref1 (V)
Vinj1 (V)
Vch1 (V)
f (Hz)
182
Annexes ___________________________________________________________________________
Vs1,2,3 (V)
(V)
Vref1 Vinj1
(V)
Vref2 Vinj2
Vch1,2,3 (V)
t (s)
Fig. A-II.5 : Analyse temporelle de la compensation du creux de tension biphas (RSTam) (Lgende - phase 1 : ., phase 2 : - - - - -, phase 3 : -..-..-..-) Face aux mmes changements brutaux de la consigne, la rapidit de la rponse des mthodes didentification et de rgulation de tension RSTam proposes sont dmontres sur la Fig. A-II.5. Afin dvaluer la prcision des mthodes didentification et de rgulation proposes, une analyse spectrale de la Phase 1 a t effectue sur la Fig. A-II.6.
V s1 (V)
16
V ref1
(V)
V inj1 (V)
V ch1 (V)
f (Hz)
183
Annexes ___________________________________________________________________________
(V)
Vref1 Vinj1
(V)
Vref2 Vinj2
(V)
Vref3 Vinj3
Vch1,2,3 (V)
t (s)
Vref1 (V)
Vinj1 (V)
Vch1 (V)
f (Hz)
Fig. A-II.8 : Analyse spectrale de la compensation du creux de tension triphas (PID) 184
Annexes ___________________________________________________________________________
La Fig. A-II.7 montre la rapidit, des mthodes didentification et de rgulation PID face un changement brutal de la consigne. Afin dvaluer la prcision des mthodes didentification et de rgulation employes, une analyse spectrale a t effectue sur la Fig. A-II.8. Cette figure tmoigne de la bonne prcision des mthodes didentification et de rgulation proposes.
(V)
Vref1 Vinj1
(V)
Vref2 Vinj2
(V)
Vref3 Vinj3
Vch1,2,3 (V)
t (s)
Vref1 (V)
Vinj1 (V)
Vch1 (V)
f (Hz)
185
Annexes ___________________________________________________________________________
Face aux changements brutaux de la consigne, la rapidit de la rponse des mthodes didentification et de rgulation de tension RST, proposes, peuvent tre observes sur la Fig. A-II.9. La Fig. A-II.10 montre lanalyse spectrale effectue sur la Phase 1 afin dvaluer la prcision des mthodes didentification et de rgulation employes.
(V)
Vref1 Vinj1
(V)
Vref2 Vinj2
(V)
Vref3 Vinj3
(V)
Vch1,2,3
t (s)
Fig. Fig. A-II.11: Analyse temporelle de la compensation du creux de tension triphas (RSTam)
Vs1 (V)
Vref1 (V)
Vinj1 (V)
Vch1 (V)
f (Hz)
Fig. A-II.12 : Analyse spectrale de la compensation du creux de tension triphas (RSTam) 186
Annexes ___________________________________________________________________________
Face aux mmes changements brutaux de la consigne, la rapidit de la rponse des mthodes didentification et de rgulation de tension RSTam proposes peuvent tre observes sur la Fig. A-II.11. La Fig. A-II.12 montre lanalyse spectrale de la Phase 1 afin dvaluer la prcision des mthodes didentification et de rgulation proposes. Cette figure tmoigne de la bonne prcision des mthodes didentification et de rgulation proposes.
187
Annexes ___________________________________________________________________________
La Fig. [Link].1 prsente lanalyse temporelle du courant de la phase 1 de la charge (Is1), des courants superposs identifi (Iref1) et inject (Iinj1) et du courant du rseau aprs filtrage (IL1). A partir de la Fig. [Link].1, on voit bien lefficacit de la mthode de rgulation propose, base du rgulateur RST amlior. En effet, comme le confirme galement la Fig. IV.6, il ny a pratiquement pas de dphasage entre le courant identifi et celui inject. Lallure du courant du rseau aprs filtrage prouve la trs bonne qualit de filtrage du filtre actif parallle. La Fig. [Link].2 prsente lanalyse spectrale du courant de la phase 1 de la charge et celui du rseau aprs filtrage. Le THD du courant ct charge qui tait est de 23% passe 0,91% ct rseau ce qui tmoigne de la bonne qualit de filtrage, dj constate par lanalyse temporelle. Remarque Il est noter que la puissance de court-circuit du rseau lectrique doit tre plus leve que la puissance apparente de la charge. Dans ce chapitre, le choix de la puissance du rseau a pour objectif dtudier la robustesse du rgulateur RSTam travers le changement de limpdance de court-circuit du rseau. En effet, le changement de limpdance de court-circuit du rseau lectrique 188
Annexes ___________________________________________________________________________
fait changer le systme dorigine reprsent par le filtre en T en srie avec limpdance du rseau (Cf. [Link].4). Pour cette raison, la puissance du rseau lectrique est quelquefois proche de la puissance de la charge polluante.
Is1 (A)
(A)
Iref1 Iinj1
IL1 (A)
t (s)
(A)
IL1
(A)
f (Hz)
189
Annexes ___________________________________________________________________________
2.1.2 Simulation par Matlab- Simulink- Power System Blockset Le mme cas tudi prcdemment est simul ici par Matlab- Simulink- Power System Blockset. Les Fig. [Link].3 et Fig. [Link].4 montrent les analyses temporelle et spectrale du courant de la charge Is1 et celui du rseau IL1. Les THD, du courant de la charge 24% et du rseau aprs compensation 1,2%, sont donns sur la Fig. [Link].5.
% du fondamental (A)
% % du fondamental
190
Annexes ___________________________________________________________________________
Fig. [Link].5 : Le THD des courants de la charge et du rseau aprs filtrage (a=0) En comparant ces rsultats de simulation avec ceux obtenus travers SABER, nous trouvons que la modlisation par Matlab- Simulink- Power System Blockset des parties construisant la structure tudie est correcte et les rsultats de simulation donns par ce logiciel sont donc crdibles. Remarque Tout le reste des simulations qui vont tre prsentes dans ce chapitre a t effectu sous SABER.
Daprs lanalyse temporelle de la Fig. [Link].6, nous remarquons que la bonne qualit de filtrage se maintient malgr laugmentation de langle dallumage (a) qui dgrade gnralement la qualit de compensation. Daprs lanalyse spectrale de la Fig. [Link].7, Le THD du courant est de 31% ct charge, et se trouve rduit 1,08% ct rseau.
191
Annexes ___________________________________________________________________________
Is1 (A)
IL1 (A)
t (s)
(A)
IL1
(A)
f (Hz)
192
Annexes ___________________________________________________________________________
La Fig. [Link].8 prsente lanalyse temporelle de la tension continue (Vdc), du courant de la charge (Is1), du courant inject (Iinj1) et du courant du rseau (IL1) aprs filtrage. A partir de la Fig. [Link].8, on remarque que la boucle de rgulation de la tension continue arrive poursuivre une rfrence de Vdc-ref =840 V avec un temps de rponse de 20 ms, un dpassement de 7% et un taux dondulation de 1,5%. Ce temps de rponse est largement infrieur aux 100 ms requis par le cahier des charges. Quant au courant inject par le filtre actif, il suit toujours le courant de rfrence en gardant quasiment la mme rapidit et robustesse.
Vdc (V)
Is1 (A)
(A)
Iref1 Iinj1
IL1 (A)
t (s)
193
Annexes ___________________________________________________________________________
La Fig. [Link].9 prsente lanalyse temporelle de la tension de la phase 1 au point de raccordement (Vs1), du courant de la charge (Is1), des courants superposs identifi (Iref1) et inject (Iinj1) et du courant du rseau aprs filtrage (IL1). Nous remarquons, en observant le courant du rseau aprs filtrage, que la qualit du filtrage est toujours maintenue, malgr la dformation de la tension dalimentation du rseau. Cela sexplique par lajout de la tension du rseau (tension compensatrice) dans le schma de rgulation et par la compensation en boucle ferme de ce type de perturbation grce au rgulateur RSTam (Cf. Fig. IV.7). Ce rsultat dmontre galement la bonne performance du systme base de PLL dans un nouveau cas de fonctionnement extrme du rseau lectrique. La Fig. [Link].10 prsente lanalyse spectrale, phase 1, de la tension de raccordement, du courant de la charge et de celui du rseau aprs filtrage. La tension au point de raccordement contient principalement des harmoniques de tension de rangs 5, 7, 11, 13, 17, 19, 23 occasionnant un THDv global de 8%. Le THD du courant ct charge est de 30,5%, de 1,45% ct rseau. Cela montre que linfluence dune telle dformation de tension du rseau est largement matrise.
Vs1 (V)
Is1 (A)
(A)
Iref1 Iinj1
IL1 (A)
t (s)
194
Annexes ___________________________________________________________________________
Vs1 (V)
Is1 (A)
IL1 (A)
f (Hz)
(A)
Vs1-PLL (V)
(A)
Iref1 Iinj1
IL1 (A)
t (s)
Fig. [Link].11 : Analyse temporelle de la compensation (frquence 60 Hz) La Fig. [Link].11 prsente lanalyse temporelle, de la phase 1, du courant de la charge (Is1), de la tension au point de raccordement en sortie de la PLL (Vs1-PLL), des courants superposs identifi (Iref1) et inject (Iinj1) et du courant du rseau aprs filtrage (IL1). Cette figure tmoigne de lefficacit du systme base de PLL sous un tel changement en frquence du rseau lectrique. 195
Annexes ___________________________________________________________________________
Is1
(A)
IL1
(A)
f (Hz)
Fig. [Link].12 : Analyse spectrale de la compensation (frquence 60 Hz) La Fig. [Link].12 prsente lanalyse spectrale du courant de la phase 1 de la charge et du celui du rseau aprs filtrage. Le THD du courant est de 31% ct charge, et de 2,26% ct rseau, ce qui constitue une bonne qualit de compensation malgr laugmentation relative du THD du courant aprs filtrage.
196
Annexes ___________________________________________________________________________
(A)
IL1 (A)
(1000 kVA)
(A)
t (s)
Fig. [Link].13 : Analyse temporelle de la compensation harmonique en prsence des lments passifs en rseau lectrique Daprs lallure des courants aprs compensation, aucun risque de rsonance ne serait craindre. Cela sexplique par un dpassement ngligeable la frquence de coupure dans le diagramme de gain de la fonction de transfert entre le courant de rfrence et le courant inject (Cf. Fig. IV.6).
Annexes ___________________________________________________________________________ filtre actif parallle. Le cas tudi est modlis par un rseau lectrique et une charge non linaire dont les caractristiques sont donnes dans le tableau suivant :
Rseau lectrique Charge non-linaire
2 MVA, 240 V, 50 Hz
Quatre modules de filtrage sont employs pour compenser chacun 1/4 du courant total identifi. La Fig. [Link].14 prsente lanalyse temporelle, de la phase 1, du courant de la charge (Is1), de la tension au point de raccordement donne par la PLL (Vs1-PLL), de la superposition du courant total identifi (Iref1-tot) et celui inject par les quatre modules (Iinj1-tot). Cette figure prsente aussi le courant du rseau aprs filtrage (IL1). Daprs la Fig. [Link].14, lallure du courant du rseau aprs filtrage tmoigne de la bonne qualit de compensation du filtre actif parallle. Un faible dphasage apparat entre le courant total identifi et celui inject. Ce phnomne est d au montage par module des filtres que nous avons adopt ici. En effet, en choisissant autant de filtres de sortie et de boucles de rgulation de courant que le nombre de modules de filtre, les dphasages occasionns dans chaque boucle de rgulation sajoutent, affectant la forme du courant total inject. Cette stratgie de montage a t adopte, dune part pour profiter de la marge de compensation, dautre part, pour viter lemploi dun filtre de sortie commun qui entranerait un dimensionnement adapt la puissance de la charge polluante. La Fig. [Link].15 prsente lanalyse spectrale, de la phase 1, du courant de la charge et de celui du rseau aprs filtrage. Le THD du courant ct charge est de 32%, de 3,5% ct rseau, ce qui prouve la bonne qualit de filtrage remarque de lanalyse temporelle.
Is1 (A)
Vs1-PLL (V)
(A)
Iref1_tot Iinj1_tot
IL1 (A)
t (s)
Fig. [Link].14 : Analyse temporelle de la compensation harmonique (cas de quatre filtres actifs en parallle)
198
Annexes ___________________________________________________________________________
Is1
(A)
IL1
(A)
f (Hz)
Fig. [Link].15 : Analyse spectrale de la compensation harmonique (cas de quatre filtres actifs en parallle)
199
Annexes ___________________________________________________________________________
La Fig. [Link].16 prsente lanalyse temporelle de la phase 1. A partir de la Fig. [Link].16, on constate que le filtre actif arrive difficilement suivre le courant harmonique de la charge polluante. Ce phnomne est d en partie laugmentation de 15% de la tension dalimentation, mais galement la diminution de la frquence fondamentale (50 Hz) de 10%. Nous remarquons, par contre, que le systme base de PLL reste robuste face aux changements des paramtres du systme dorigine (Vs) et la variation des perturbations rejeter (en frquence et en amplitude).
Is1 (A)
Vs1-PLL (V)
(A)
Iref1 Iinj1
IL1 (A)
t (s)
Fig. [Link].16 : Analyse temporelle de la compensation harmonique (cas dune tension de 265 V et dune frquence de 45 Hz) La Fig. [Link].17 prsente lanalyse spectrale du courant de la charge et de celui du rseau aprs filtrage. Le THD du courant ct charge est de 30%, de 5,9% ct rseau, ce qui est loin des (10%) requis par le cahier des charges.
200
Annexes ___________________________________________________________________________
Is1
(A)
IL1
(A)
f (Hz)
Fig. [Link].17 : Analyse spectrale de la compensation harmonique (cas dune tension de 265 V et dune frquence de 45 Hz)
La Fig. [Link].18 prsente lanalyse temporelle de la phase 1. Daprs cette figure, en observant lallure du courant rseau aprs compensation, on estime que le filtre actif retrouve une bonne qualit de filtrage, par rapport au cas prcdent, malgr les conditions extrmes de fonctionnement. En effet, malgr laugmentation de la tension dalimentation de 15%, et de la frquence fondamentale (60 Hz) denviron 10%, on montre bien que le filtre actif parallle compense correctement ces perturbations. Ces bons rsultats peuvent sexpliquer par une monte moins rapide des harmoniques par rapport au cas prcdent. On remarque aussi la bonne prcision du systme base de PLL dans ce cas extrme. La Fig. [Link].19 prsente lanalyse spectrale du courant de la charge et celui du rseau aprs filtrage. Le THD du courant ct charge est de 30,5%, de 2,69% ct rseau, ce qui est largement loin des 10% requis par le cahier des charges.
201
Annexes ___________________________________________________________________________
Is1 (A)
Vs1-PLL (V)
(A)
Iref1 Iinj1
IL1 (A)
t (s)
Fig. [Link].118 : Analyse temporelle de la compensation harmonique (cas dune tension de 265 V et dune frquence de 65 Hz)
Is1
(A)
IL1
(A)
f (Hz)
Fig. [Link].19 : Analyse spectrale de la compensation harmonique (cas dune tension de 265 V et dune frquence de 65 Hz
202
Rsum
Dans un rseau lectrique basse tension, les harmoniques, la puissance ractive, les dsquilibres de courant et de tension, et les creux de tension ont des effets nfastes sur les quipements lectriques. Dans cette thse, nous avons abord les structures de filtres actifs parallles, sries et combines (UPQC) en vue dtudier la compensation de tous les types de perturbations susceptibles dapparatre dans le rseau lectrique. De plus, dans lobjectif damliorer la qualit de lnergie lectrique, en la rendant conforme aux nouvelles contraintes normatives, nous avons galement t amen dvelopper et appliquer des mthodes dautomatique avances. Dans ce cadre, cette thse rapporte de plusieurs structures de rgulation, notamment dune version amlior du RST, cette dernire repoussant les limites dapplication du RST classique. De plus, dans lobjectif de raliser un filtre actif parallle, une tude approfondie de faisabilit exprimentale a t effectue, prenant en compte les contraintes industrielles du produit aussi bien au niveau de sa conception que de son application.
Mots cls
filtre actif parallle, filtre actif srie, UPQC, harmoniques, dsquilibre, puissance ractive, creux de tension, mthodes automatiques, RST, faisabilit exprimentale, problmatiques industrielles.
Abstract
In an electrical network, unbalanced and harmonic currents generated by nonlinear and unbalanced loads can cause harmonics and unbalanced voltage. These voltage perturbations along with voltage sag can strongly degrade customer power quality. Nowadays, the active filtres such as shunt, series and combined shunt-series are studied as a flexible solution to compensate all current and voltage perturbations. Therefore, in order to improve power system quality, we have studied advanced control methods applied on the structures of the above mentioned active filtres. For such requirements, we have developed a new RST control method in order to extend the application spectrum of RST controller. Thus, an universal shunt active filter has been developed using a continuous and discret time experimental study.
Key words
Shunt actif filter, series actif filtre, UPQC, harmonics, unbalance, reactive power, voltage sag, control, poles placement, experimental feasibility, industrial constraints.