Commande Intelligente d'un Système PV
Commande Intelligente d'un Système PV
Doctorat en Sciences
Spécialité Electronique
Présentée par
Nadjwa CHETTIBI
Directeur de Thèse
Pr. Adel MELLIT
Avant tout, je remercie Dieu le tout puissant qui m’a donné la force, la
patience et le courage pour accomplir ce modeste travail.
Keywords: Photovoltaic, Wind turbine, Fuel cell, PSO algorithm, Fuzzy logic, Artificial neural
network, ANFIS.
تصاعد في السنوات األخيرة االهتمام بدمج نظم إنتاج الطاقة الالمركزية القائمة على الموارد الشمسية و توربينات الرياح مع شبكة الكهرباء:ملخص
لقد برز مفهوم الشبكة الصغيرة. ولكن الطبيعة المتقطعة والمتقلبة إلنتاج الطاقة المتجددة تؤثر سلبا على استقرار وسيرورة عمل الشبكة الكهربائية.
تلعب واجهة الكترونيك. كحل واعد لهذه المشكلة، التي تشتمل علي مصادر الطاقة الموزعة باإلضافة الي أنظمة تخزين الطاقة وشحن االستهالك
. ألنها تمم خوارزميات المراقبة و التحكم التي من شأنها ضمان اإلدارة المثلى و الذكية للطاقة،الطاقة دورا رئيسيا في هذا النوع من نظم التوليد
الهجينة كما طبقنا استراتيجيات للتحكم واإلدارة مبنية على بعض تقنيات الذكاء االصطناعيPV درسنا مختلف أنظمة التوليد،في هذه األطروحة
وقد أظهرت. كالمنطق الغامض والشبكات العصبية االصطناعية وذلك بهدف تحسين المردود الطاقوي و نوعية الطاقة في ظل ظروف تشغيل مختلفة
. فعالية طرق التحكم التي اقترحناهاMatlab/Simulink المحاكاة بواسطة
الشبكة العصبية، المنطق الغامض، PSO خوارزمية، بطاريات الوقود، توربينات الرياح، ال كهروضوئيىة:كلمات مفتاحية
.ANFIS،االصطناعية
i
Table des Matières
Introduction Générale....................................................................................................................................................................... 1
Chapitre I
Chapitre II
II.1. Introduction................................................................................................................................................................................. 23
ii
II.2. Partie 1 : Modélisation mathématique des éléments des systèmes de production. ................................................... 23
II.2.1. Source photovoltaïque . .................................................................................................................................................. 23
II.2.1.1. Modélisation de la cellule PV .............................................................................................................................. 23
II.2.1.2. Simulation du module PV .................................................................................................................................... 25
II.2.2. Modélisation des composants d’une éolienne .......................................................................................................... 26
II.2.2.1. Turbine éolienne ................................................................................................................................................... 26
II.2.2.2. Multiplicateur de vitesse ...................................................................................................................................... 27
II.2.2.3. Arbre mécanique .................................................................................................................................................. 27
II.2.2.4. Génératrice synchrone à aimants permanents .............................................................................................. 28
II.2.3. Pile à combustible à oxyde solide................................................................................................................................. 30
II.2.3.1. Modélisation dynamique du stack de PAC-OS ............................................................................................. 30
II.2.3.2. Contraintes sur le fonctionnement de la PAC-OS ........................................................................................ 33
II.2.4. Modèle de Simulation de batterie ................................................................................................................................ 35
II.2.5. Modélisation des convertisseurs de puissance ......................................................................................................... 35
II.2.5.1. Convertisseur élévateur ............................................................................................................................................. 36
II.2.5.2. Convertisseur réversible en courant ....................................................................................................................... 37
II.2.5.3. Convertisseur triphasé à deux niveaux .......................................................................................................... 38
II.2.5.4. Convertisseur triphasé à trois niveaux ........................................................................................................... 40
II.3. Partie 2 : Aperçu sur les techniques d’intelligence artificielle ...................................................................................... 42
II.3.1. Logique floue .................................................................................................................................................................... 42
II.3.1.1. Ensemble flou .......................................................................................................................................................... 42
II.3.1.2. Fonction d’appartenance .................................................................................................................................... 42
II.3.1.3. Variable linguistique ............................................................................................................................................. 43
II.3.1.4. Contrôleur par logique floue ............................................................................................................................. 43
II.3.2. Réseaux de neurones artificiels .................................................................................................................................. 46
II.3.2.1. Neurone formel ..................................................................................................................................................... 46
II.3.2.2. Structures des RNAs ............................................................................................................................................. 47
II.[Link]. Réseaux de neurones non bouclés .................................................................................................. 47
II.[Link].1. Réseaux multicouches ................................................................................................................ 47
II.[Link].2. Réseaux à fonctions radiales de base...................................................................................... 48
II.[Link]. Réseaux de neurones bouclés ......................................................................................................... 49
II.3.2.3. Apprentissage des RNAs ....................................................................................................................................... 50
II.3.2.4. Règles d’apprentissage supervisé ........................................................................................................................ 51
II.[Link]. Algorithme de rétropropagation du gradient ................................................................................ 51
II.[Link]. Algorithme des moindres carrés ...................................................................................................... 52
II.[Link].1. Réseau ADALINE ......................................................................................................................... 52
II.[Link].2. Règle μ-LMS ................................................................................................................................. 53
II.[Link].3. Règle α-LMS ................................................................................................................................ 53
II.3.3. Système d’inférence neuro-flou adaptatif ANFIS .................................................................................................... 53
II.4. Conclusion .................................................................................................................................................................................. 54
Chapitre III
iii
III.2.1. Configuration du système PV étudié ......................................................................................................................... 55
III.2.2. Structure de commande du SPVCR ........................................................................................................................... 56
III.2.2.1. Commande du convertisseur élévateur ......................................................................................................... 56
III.2.2.2. Contrôle direct de puissance de l’onduleur .................................................................................................. 58
III.[Link]. Synthèse des régulateurs des puissances P et Q ........................................................................... 59
III.[Link]. Régulation de la tension du bus CC à base de PSO .................................................................... 61
III.2.3. Résultats de simulation ................................................................................................................................................ 64
III.3. Commande d’un système de génération éolien connecté au réseau électrique ...................................................... 66
III.3.1. Topologie du système de conversion éolien .......................................................................................................... 66
III.3.2. Contrôle du convertisseur côté génératrice............................................................................................................. 66
III.3.2.1. Contrôle directe du couple ............................................................................................................................... 66
III.3.2.2. Poursuite du PPM du GE .................................................................................................................................... 68
III.[Link]. Méthode de MPPT avec ratio de vitesse optimal ......................................................................... 69
III.[Link]. Méthode de commande MPPT adoptée ......................................................................................... 69
III.3.3. Commande à tension orientée de l’onduleur .......................................................................................................... 71
III.3.3.1. Contrôle de la tension du bus continu ............................................................................................................ 72
III.3.3.2. Optimisation des contrôleurs des courants d et q........................................................................................ 73
III.3.4. Synthèse des régulateurs PI à l’aide de PSO ............................................................................................................. 74
III.3.5. Résultats de simulation ................................................................................................................................................ 74
[Link] par LF d’un système de génération basé sur une PAC-OS........................................................................ 77
[Link] du système de production à PAC-OS ............................................................................................. 78
III.4.2. Système de commande du conditionneur de puissance ...................................................................................... 78
[Link]ôle du convertisseur élévateur ............................................................................................................... 78
III.4.2.2. Commande en puissance de l’onduleur ........................................................................................................ 80
III.4.3.Résultats de simulation .................................................................................................................................................. 81
III.5. Conclusion ............................................................................................................................................................................... 83
Chapitre IV
iv
IV.5.6. Gestion centrale d’énergie à base de la LF .............................................................................................................. 103
IV.6.Résultats de Simulation ........................................................................................................................................................... 105
[Link] .................................................................................................................................................................................. 113
Chapitre V
v
Liste des Figures
Figure I.1 : Coupe transversale d’une cellule PV typique. ........................................................................................................ 6
Figure I.2 : Exemples de trois types de cellules PV à base de Silicium : a) Monocristallin, b) Amorphe en films
minces, b) Multicristallin. ................................................................................................................................................................ 6
Figure I.3 : Cellule, module et champ PV . ................................................................................................................................... 6
Figure I.4 : Caractéristiques électriques du module BPMSX120 : a) I(V), b)P(V) pour T=25°C et G=1000W/m2..... 7
Figure I.5 : Couplage directe GPV-charge. .................................................................................................................................. 8
Figure I.6 : Configurations des systèmes PV avec : a) Un seule étage d’adaptation CC-CC, b) Un seule étage de
conversion CC-CA , c) Deux étages de conversion CC-CC et CC-CA . ............................................................................... 8
Figure I.7 : Schéma descriptif d’un aérogénérateur . ................................................................................................................ 9
Figure I.8: Types de turbines éoliennes a) Savonius, b) Darrieus, c) Evans, d) Tripale, e) Bipale, f) Mono-pale.. ..... 9
Figure I.9 : Diagramme de la puissance mécanique en fonction de la vitesse du vent. ................................................... 11
Figure I.10 : Chaines de conversion d’une éolienne à vitesse variable a) à deux convertisseurs, b) à trois
convertisseurs de puissance............................................................................................................................................................. 11
Figure I.11 : Composition d'un stack de PAC............................................................................................................................... 12
Figure I.12 : Principe de fonctionnement d’une PAC-OS. ....................................................................................................... 14
Figure I.13 : Structure d’un système de production à PAC : a) avec un onduleur seulement, b) avec un hacheur
et un onduleur..................................................................................................................................................................................... 14
Figure I.14 : Classification des systèmes de production hybrides. ......................................................................................... 16
Figure I.15 : Configuration d’un SEH (PV-éolien-stockage) : a) Connexion au bus CC, b) Connexion au bus CC
et CA, c) Connexion au bus CA. ...................................................................................................................................................... 19
Figure I.16 : Paradigme de contrôle centralisé d’un SEH. ........................................................................................................ 20
Figure I.17 : Creux de tension. ........................................................................................................................................................ 21
Figure II.1 : Circuit équivalant simplifié de la cellule PV. ........................................................................................................ 23
Figure II.2 : Bloc Simulink du GPV................................................................................................................................................. 25
Figure II.3 : Caractéristiques IPV (VPV) et PPV (VPV) du module BP MSX120 pour différents niveaux d’éclairement
à température constante T= 250C.................................................................................................................................................. 25
Figure II.4 : Caractéristiques IPV(VPV) et PPV(VPV) du module BP MSX120 pour différentes valeurs de température
et à G=1000W/m2. ............................................................................................................................................................................. 25
Figure II.5 : Modèle de la TE implémenté dans Matlab/Simulink ......................................................................... 27
Figure II.6 : Courbe de a) Puissance versus vitesse de rotation de la TE, b) Coefficient de puissance. ......................... 27
Figure II.7 : Représentation des repères fixe et tournant ......................................................................................................... 29
Figure II.8 : Exemple d’un canal de gaz à un seul orifice ....................................................................................................... 31
Figure II.9 : Modèle dynamique du stack de la PAC-OS ......................................................................................................... 32
Figure II.10 : Contrôle de la PAC-OS............................................................................................................................................. 34
Figure II.11 : Modèle de PAC-OS établi dans l’environnement Matlab/Simulink. ............................................................ 34
Figure II.12 : Evolution de : a) Pressions partielles (PH2,PO2,PH2O) dans la PAC-OS, b) Facteur d’utilisation de
carburant.. ............................................................................................................................................................................................ 34
Figure II.13 Variation de : a) Tension de sortie b) Courant de PAC-OS, c) Limites du courant..................................... 35
Figure II.14 : Modèle dynamique de batterie. ............................................................................................................................. 35
Figure II.15 : Circuit électrique d’un convertisseur élévateur ............................................................................................... 36
Figure II.16 : Circuit équivalent du convertisseur boost avec : a) Q1 passant, b) Q1 bloquant .................................... 36
Figure II.19 : Rapport cyclique d’un signal MLI ........................................................................................................................ 36
Figure II.17 : Circuit de puissance d’un convertisseur réversible abaisseur-élévateur .................................................. 38
vi
Figure II.18 : Circuit de puissance de l’onduleur à deux niveaux couplé au réseau ...................................................... 38
Figure II.19 : Circuit de puissance d’un onduleur 3L-NPC. .................................................................................................... 41
Figure II.21 : Différence entre la notion d’ensemble classique et d’ensemble flou. .......................................................... 42
Figure II.22 : Exemples de types de FA : (a) Fonction-L; (b) Fonction-Λ (triangulaire); (c) Fonction Gaussienne;
(d) Fonction-Π (trapézoïdale) . ...................................................................................................................................................... 43
Figure II.23 : FAs de la variable linguistique « Vitesse ». .......................................................................................................... 43
Figure II.24 : Structure d’un contrôleur par logique floue. .................................................................................................... 44
Figure II.25 : Neurone biologique et sa contrepartie artificielle ........................................................................................... 46
Figure II.26 : Perceptron à trois couches ..................................................................................................................................... 48
Figure II.27 : Structure d’un réseau RBFN .................................................................................................................................. 49
Figure II.28 : Structure d’un réseau d’Elman. ............................................................................................................................. 49
Figure II.29 : Types d’apprentissage : a) Supervisé, b) Non-supervisé ............................................................................... 50
Figure II.30 : Réseau ADALINE entrainé par la règle LMS. ..................................................................................................... 52
Figure II.31: Architecture d’un ANFIS .......................................................................................................................................... 54
Figure III.1 : Structure du SPVCR étudié. ...................................................................................................................................... 56
Figure III.2 : Position du point de fonctionnement sur la courbe PPV(VPV) .......................................................................... 57
Figure III.3 : Organigramme d’algorithme d’IncCond . ........................................................................................................... 57
Figure III.4 : Schéma de commande du convertisseur élévateur. .......................................................................................... 57
Figure III.5 : Schéma de DPC classique ......................................................................................................................................... 58
Figure III.6. Schéma de DPC adopté pour l’onduleur côté réseau. ........................................................................................ 59
Figure III.7 : Schéma de la boucle de régulation de puissance active (P) ........................................................................... 60
Figure III.8 : Déplacement d’une particule de l’essaim dans l’espace de recherche ........................................................ 61
Figure III.9 : Organigramme de l’algorithme de PSO .............................................................................................................. 63
Figure III.10 : Variation des puissances active et réactive sous un niveau d’éclairement variable. ............................. 65
Figure III.11 : Évolution de la tension du bus CC. ...................................................................................................................... 65
Figure III.12 : Variation du courant triphasé injecté au réseau. ............................................................................................ 65
Figure III.13 : Topologie du système de production éolien étudié. ....................................................................................... 66
Figure III.14 : Schéma classique de DTC d’une GSAP. ............................................................................................................. 67
Figure III.15 : Schéma de DTC adopté pour le pilotage du redresseur. ............................................................................... 68
Figure III.16 : Courbes de puissance Pt(ωt) d’une TE pour différentes vitesses de vent ................................................... 69
Figure III.17: Schéma de commande MPPT par TSR optimal ................................................................................................. 69
Figure III.18 : Structure de FFNN adopté pour la MPPT du GE............................................................................................... 70
Figure III.19 : Schéma de VOC de l’onduleur connecté au réseau ....................................................................................... 71
Figure III.20 : Orientation du repère dq selon le vecteur de tension du réseau ............................................................... 72
Figure III.21 : Bloc diagramme de la boucle de régulation de tension du bus CC. ........................................................... 73
Figure III.22 : Evolution de a) la puissance mécanique, b) la vitesse de rotation de la TE dans le premier cas de
simulation. ............................................................................................................................................................................................ 75
Figure III.23 : Evolution du : a) Couple de la TE b) Module du flux dans le premier cas de simulation. .................... 76
Figure III.24 : Evolution de : a) Composante d du courant réseau, b) Composante q du courant réseau, c)
Tension du bus CC, d) Puissances active et réactive dans le premier cas de simulation. ................................................ 76
Figure III.25 : Evolution de : a) Couple Te*, Te et Tt , b) Vitesse de rotation de la TE dans le deuxième cas de
simulation. ............................................................................................................................................................................................ 77
Figure III.26 : Evolution de : a) Composantes d et q du courant réseau , b) Tension du bus CC dans le deuxième
cas de simulation. ............................................................................................................................................................................... 77
Figure III.27 : Structure du système de GD basé sur une PAC-OS. ...................................................................................... 78
vii
Figure III.28 : Schéma de contrôle du convertisseur boost. .................................................................................................... 79
Figure III.29 : FAs des entrées et de la sortie du CLF de la tension du bus CC. ................................................................... 79
Figure III.30 : Commande en mode PQ de l’onduleur 3L-NPC. ............................................................................................. 80
Figure III.31. FAs des entrées et de la sortie du CLF de puissance P et Q. ............................................................................ 81
Figure III.32 : Evolution de la tension et du courant de sortie de la PAC-OS durant le premier cas de simulation 82
Figure III.33 : Evolution de la tension du bus CC durant le premier cas de simulation. ................................................. 82
Figure III.34 : Variation des puissances active et réactive de l’onduleur, du réseau et de la charge CA durant le
premier cas de simulation. ............................................................................................................................................................... 83
Figure III.35 : Evolution de a) Puissances active et réactive d’onduleur, b) Tension du bus-CC durant le
deuxième cas de simulation............................................................................................................................................................. 83
Figure IV.1. Bloc diagramme du MR (CA/CC) hybride. ............................................................................................................ 90
Figure IV.2: Contrôleur SN-RBFN du GPV implémenté dans l’environnement Matlab/Simulink. ............................... 92
Figure IV.3 : Schéma de commande prédictive du SGEE. ......................................................................................................... 93
Figure IV.4 : Prédiction de la consigne future du couple Te* ................................................................................................. 95
Figure IV.5 : Schéma de commande du convertisseur de la PAC-OS dans Matlab/Simulink......................................... 96
Figure IV.6 : Schéma de commande du convertisseur interfaçant le SSEB dans Matlab/Simulink. .............................. 97
Figure IV.7 : Orientation du FV du réseau. .................................................................................................................................. 98
Figure IV.8 : Bloc diagramme du schéma VFDPC-SVM proposé pour l’onduleur. .......................................................... 99
Figure IV.9 : Structure du réseau ENN adopté. ........................................................................................................................... 101
Figure IV.10 : Exemple d’ajustement en-ligne des poids du FFNN par AI dans Matlab/Simulink . .............................. 103
Figure IV.11 : FAs des variables a) Pnet , b) SOC, c) Pb* , d) Pfc*.................................................................................................. 105
Figure IV.12: Variation de : a) Irradiation solaire, b) Courant PV, c) Tension PV, d) Tension de référence en
sortie du SN-RBFN, e) Puissance PV de sortie. ............................................................................................................................ 106
Figure IV.13 : Evolution de : a) Vitesse de rotation, b) Couple électromagnétique et courant statorique direct, c)
module du flux statorique d) Tension du bus-CC. ................................................................................................................... 107
Figure IV.14 : Variation de : a) Demande des charges CA et CC, b) Puissances actives de l’onduleur et du réseau,
c) Puissances de sortie de PAC-OS et du SSEB, d) Puissances réactives de l’onduleur, de charge CA et du réseau
principal................................................................................................................................................................................................ 109
Figure IV.15 : Comparaison entre le SGE classique et le SGE flou. ........................................................................................ 109
Figure IV.16 : Evolution de la : a) Tension du bus CC, b) Puissances de charge total, du MR et du réseau
principal................................................................................................................................................................................................ 110
Figure IV.17 : Variation de : a) Puissances de charge CA, de l’onduleur et du réseau, b) Puissances de la PAC-OS
et du SSEB c) Tension du bus CC, d) Courant de sortie d’onduleur en présence de creux de tension, e)
Composantes αβ négatives du courant d’onduleur. ................................................................................................................ 111
Figure IV.18 : Evolution de : a) Tension PV, b) Profil du vent, c) Vitesse de rotation du GE, d) Erreur statique de
vitesse, e) Couple de la GSAP, f) Puissances active et réactive de l’onduleur...................................................................... 113
Figure V.1 : Structure adoptée du SEH connecté au réseau..................................................................................................... 115
Figure V.2 : Organigramme de l’algorithme de P&O ................................................................................................................ 116
Figure V.3 : Schéma bloc de commande à base d’ENN du convertisseur interfaçant le GPV ......................................... 117
Figure V.4 : Implémentation dans Matlab/Simulink de : a) Contrôleur ENN, b) la règle d’adaptation des poids de
la CS. ...................................................................................................................................................................................................... 119
Figure V.5: Schéma bloc de la commande du convertisseur boost de la PAC-OS. ............................................................ 120
Figure V.6: FAs assignées aux (a) Variables d’entrée (efc et Δefc), (b) Variables de sortie (KP,KI) du RGF du
contrôleur de PAC-OS....................................................................................................................................................................... 121
Figure V.7: Schéma bloc de commande du convertisseur réversible. ................................................................................... 122
viii
Figure V.8 : Schéma de VFOC adopté pour la commande de l’onduleur. ............................................................................ 124
Figure V.9 : Schéma bloc du CFMG implémenté dans Matlab/Simulink. ............................................................................ 125
Figure V.10 : FAs des entrées du CFMG : a) Surface de glissement S, b) Sa variation dS. ................................................ 125
Figure V.11 : FAs des entrées des régulateurs des gains KP et KI de contrôleur de courant direct. .............................. 127
Figure V.12: SGE adopté: a) Organigramme de sa stratégie de supervision, b) Son implémentation dans
Matlab/Simulink. ................................................................................................................................................................................ 130
Figure V.13 : Variation de : a) Tension PV, b) Courant PV, c) Tension de référence imposée par le ENN, d)
Puissance PV de sortie. ...................................................................................................................................................................... 131
Figure V.14 : Variation de : a) Puissance active d’onduleur, de charge et du réseau, b) Puissances de PAC-OS et
du SSEB, c) Puissance réactive d’onduleur, de charge et du réseau. ..................................................................................... 132
Figure V.15 : Variation de : a) Courant en quadrature d’onduleur, b) Gains du contrôleur PI du courant en
quadrature, c) Tension du bus CC. ................................................................................................................................................ 132
Figure V.16 : Variation de a) Irradiation solaire, b) Tension PV, c) Tension de référence imposée par le ENN, d)
Courant PV, e) Puissance PV de sortie. .......................................................................................................................................... 134
Figure V.17 : Evolution de : a) Courant de PAC-OS, b) Courant du SSEB, c) Puissance active d’onduleur, de
charge et du réseau. ........................................................................................................................................................................... 135
Figure V.18 : Variation de la : a) Température, b) Puissance de sortie du GPV. ................................................................. 136
Figure V.19 : Variation de : a) Puissance active d’onduleur, de charge et du réseau, b) Puissance réactive
d’onduleur, de charge et du réseau. .............................................................................................................................................. 137
Figure V.20 : Variation de : a) Courant direct d’onduleur, b) Courant en quadrature d’onduleur, c) Gains du
contrôleur PI de courant................................................................................................................................................................... 137
Figure V.21 : Evolution de a) Tension du bus CC, b) Puissance de PAC-OS et du SSEB, c) Courant d’onduleur et
tension du réseau (phase a). ............................................................................................................................................................ 138
Figure V.22 : Etapes de développement et d’implémentation d’un contrôleur sur une carte FPGA à l’aide de XSG .140
Figure V.23 : Structure du système PV-éolien étudié. ............................................................................................................... 140
Figure V.24: Bloc diagramme de la commande du convertisseur boost .............................................................................. 141
Figure V.25: Développement à base de XSG de a) Contrôleur P&O, b) L’ensemble contrôleur PI-générateur du
signal MLI. ............................................................................................................................................................................................ 142
Figure V.26: Bloc diagramme de la méthode de CCO adoptée............................................................................................... 142
Figure V.27 : Schéma de FOC appliqué au redresseur triphasé ............................................................................................. 144
Figure V.28: Blocs XSG du a) Schéma de FOC de la GSAP, b) Contrôleur PI du courant isd c) Générateur des
signaux MLI. ........................................................................................................................................................................................ 145
Figure V.29 : Schéma de VFOC à base de Backstepping. .......................................................................................................... 146
Figure V.30: Blocs XSG de a) Transformation des coordonnées (abc/αβ)et (αβ/dq), b) Transformation des
coordonnées (dq/αβ) et (αβ/abc), c) Estimateur du FV du réseau, d) Contrôleur ADALINE ,e) Schéma global de
VFOC par Backstepping. ................................................................................................................................................................... 149
Figure V.31: Performances de poursuite du PPM du GPV dans le premier cas de simulation....................................... 150
Figure V.32: Performances de poursuite du PPM du GE dans le premier cas de simulation. ........................................ 151
Figure V.33: Evolution de : a) Composantes q du courant statorique, b) Composantes d du courant statorique de
la GSAP, c) Tension du bus CC durant le premier test. ............................................................................................................. 151
Figure V.34: Evolution de : a) Puissances active b) Puissance réactive d’onduleur, c) Composantes αβ du FV du
réseau durant le premier test. ......................................................................................................................................................... 151
Figure V.35: Performances de poursuite du PPM du GPV dans le deuxième cas de simulation. .................................. 152
Figure V.36: Performances du schéma VFOC dans le deuxième cas de simulation.......................................................... 153
ix
Liste des Tableaux
x
Liste des acronymes
xi
MLI Modulation de Largeur d’Impulsion
MPPT Maximal Power Point Tracking
MR Micro-Réseau
MSAP Machine Synchrone à Aimants Permanents
P&O Perturbation et Observation
PAC Pile A Combustible
PAC-OS Pile A Combustible A Oxyde Solide
PAFC Phosphoric Acid Fuel Cell
PCC Point of Common Coupling
PEMFC Proton Exchange Membrane Fuel Cell
PI Proportionnel Intégral
PLL Phase Locked Loop
PPM Point de Puissance Maximale
PSO Particule Swarm Optimisation
PV PhotoVoltaïque
RBFN Radial Basis Function Network
RGF Régulateur de Gains Flou
RGNF Régulateurs de Gain Neuro-Flou
RNA Réseau de Neurones Artificiel
RP Rétro-Propagation
SEH Système d’Energie Hybride
SER Source d’Energie Renouvelable
SGE Système de Gestion d’Energie
SGEE Système de Génération d’Energie Eolien
SIF Système d’Inférence Floue
SN-RBFN Single Neuron Radial Basis Function Network
SOC State Of Charge
SOFC Solid Oxide Fuel Cell
SPVCR Système PV Connecté au Réseau
SPWM Sinusoidal Pulse Width Modulation
SSEB Système de Stockage d’Energie à base de Batteries
SVM Space Vector Modulation
TIA Technique de l’Intelligence Artificielle
TE Turbine Eolienne
THD Total Harmonic Distortion
TSK Takagi-Sugeno-Kang
TSR Tip Speed Ratio
VF-DPC Virtual Flux Direct Power Control
VFDPC-SVM Virtual Flux Direct Power Control based SVM
VFOC Virtual Flux Oriented Control
VOC Voltage Oriented Control
XSG Xilink System Generator
xii
Introduction générale
2
Introduction générale
3
Introduction générale
4
Chapitre I
I.1. Introduction
Vu l’accroissement de la demande énergétique d’une part et la limitation des
ressources à base d’énergie fossile d’autre part, les sources d’énergie renouvelable (PV,
éolienne, hydraulique, biomasse, etc.) ont bénéficié d’un regain d’intérêt ces dernières
décennies. Les sources PV et éolienne présentent des alternatives de production
d’électricité très prometteuses pour des applications autonomes ou connectées au
réseau.
A partir du rayonnement solaire incident, les panneaux PV génèrent de l’électricité
par le biais de l’effet photovoltaïque. Ils sont employés actuellement dans diverses
applications à savoir les relais de communication, les sites isolés, les centrales PV
connectés au réseau, etc. Les éoliennes sont des sources d’électricité qui récupèrent des
quantités de l’énergie cinétique du vent et les convertissent en énergie électrique à
l’aide de génératrices électriques. Les critères du coût, de rendement et de la gamme
de puissance influent sur le choix de type de la génératrice électrique et de la
configuration adéquate de la chaine de conversion éolienne.
L’hybridation de ces sources d’énergie avec ou sans système de stockage d’énergie
permet d’augmenter la fiabilité et d’obtenir un meilleur rendement. Dans ce chapitre,
nous dressons un état de l’art sur les systèmes PV hybrides. D’abord, des généralités
sur les cellules PV, les éoliennes et les PACs ainsi que leurs principes de
fonctionnement sont introduits. Ensuite, nous définissons le concept de « système
d’énergie hybride » et nous classifions les différents SEHs. Les exigences de contrôle
du flux énergétique et de la qualité d’énergie électrique dans les SEHs sont également
discutées.
I.2. Générateur photovoltaïque
I.2.1. Principe de fonctionnement
Les cellules photovoltaïques, souvent désignées par « photopiles », sont composées
de semi-conducteurs, notamment le silicium, qui possèdent la propriété d’absorption
de la lumière. C’est un dispositif de conversion de l’énergie de la radiation solaire à
une énergie électrique qui met en jeu l’effet photovoltaïque. Ce phénomène physique a
été découvert pour la première fois en 1839 par le physicien français Antoine Edmond
Becquerel [14, 15]. En 1954, les laboratoires Bell ont produit la première cellule PV à
base de silicium. Elle a trouvé sa première application dans les programmes spatiaux
de l’USA pour l’alimentation des satellites [15].
Typiquement, une cellule PV est constituée de deux couches semi-conductrices
dopées différemment (un dopage de type N et un dopage de type P) et qui forment
une jonction PN. Une différence de potentiel apparait donc aux bornes de la jonction
dans laquelle un champ électrique règne. Quand la photopile est illuminée, l’énergie
des photons du rayonnement solaire absorbé est transmise aux électrons des atomes
ce qui mène à la création des paires électron-trou. La collecte des charges libérées à
l’aide de deux électrodes (positive et négative) engendre alors un courant électrique
continu en sortie de la cellule PV [14-16]. La figure (I.1) représente une coupe
transversale d’une cellule PV expliquant son principe de fonctionnement.
5
Chapitre I : Etat de l’art sur les systèmes photovoltaïques hybrides
6
Chapitre I : Etat de l’art sur les systèmes photovoltaïques hybrides
Puissance PV (W)
2.5 80
2
60
1.5
40
1
0.5 20
Tension de circuit ouvert Vppm
0
Vco
0
0 5 10 15 20 25 30 35 40 45 0 5 10 15 20 25 30 35 40 45
Tension PV (V) Tension PV (V)
(a) (b)
Figure I.4 : Caractéristiques électriques du module BPMSX120 : a) IPV(VPV), b) PPV(VPV) pour
T=25°C et G=1000W/m2.
7
Chapitre I : Etat de l’art sur les systèmes photovoltaïques hybrides
Charge CA Charge CA
GPV Charge CC GPV GPV
/réseau /réseau
8
Chapitre I : Etat de l’art sur les systèmes photovoltaïques hybrides
Multiplicateur
Système de régulation
Frein de vitesse
électrique
Nacelle
Moyeu
Génératrice électrique
Pale Système d’orientation
Mat
Armoire du couplage
Fondation au réseau
9
Chapitre I : Etat de l’art sur les systèmes photovoltaïques hybrides
10
Chapitre I : Etat de l’art sur les systèmes photovoltaïques hybrides
A l’inverse des éoliennes à vitesse fixe, les éoliennes à vitesse variable requièrent
des convertisseurs de puissance en vue d’accomplir les tâches de contrôle du flux de
puissance, de régulation MPPT et de production d’une énergie de haute qualité [25].
Puissance
mécanique (W)
I II III IV
Pn
Contrôle MPPT
0 Vd Vn Vm
Figure I.9 : Diagramme de la puissance mécanique en fonction de la vitesse du vent [25,26].
Comme l’illustre la figure (I.10), deux configurations d’une chaine de conversion
éolienne, basée sur une GSAP (fonctionnant à vitesse variable), sont couramment
utilisées :
Configuration à deux convertisseurs tête-bêche (back to back).
Configuration avec trois convertisseurs de puissance (un redresseur à pont de
diode, un convertisseur boost et un onduleur triphasé).
CA-CC CC-CA
CDC
Charge CA
GSAP /réseau
(a)
CA-CC CC-CC CC-CA
Charge CA
GSAP
/réseau
(b)
Figure I.10 : Chaines de conversion d’une éolienne à vitesse variable : a) À deux convertisseurs, b)
À trois convertisseurs de puissance [27].
11
Chapitre I : Etat de l’art sur les systèmes photovoltaïques hybrides
12
Chapitre I : Etat de l’art sur les systèmes photovoltaïques hybrides
Efficacité
~50 ~40 40-50 45-55 50-60
électrique (%)
Céramique :
Acide Carbonate
Hydroxyde Polymère zirconium
Electrolyte phosphorique fondu
de potassium solide stabilisé à
liquide liquide
l’yttrium
Température de
fonctionnement 80-260 ~200 50-80 600-700 600-1000
(°C)
Gamme de
puissance des 1-100 50-1000 10-3-1000 100-105 5-105
applications (kW)
13
Chapitre I : Etat de l’art sur les systèmes photovoltaïques hybrides
Charge CA Charge CA
PAC PAC /réseau
/réseau
(a) (b)
Figure I.13 : Structure d’un système de production à PAC : a) Avec un onduleur seulement, b) Avec
un hacheur et un onduleur [29, 34].
14
Chapitre I : Etat de l’art sur les systèmes photovoltaïques hybrides
15
Chapitre I : Etat de l’art sur les systèmes photovoltaïques hybrides
Système hybride
Autonome
Connecté au réseau
16
Chapitre I : Etat de l’art sur les systèmes photovoltaïques hybrides
d’énergie générée par les SERs [46]. Le tableau (I.2) donne un aperçu sur les différentes
combinaisons de SEHs (comportant des sources PV) proposées dans la littérature.
La remarque qui peut être faite à partir du tableau (I.2) est que les SEHs peuvent
être groupés selon leur combinaison de ressources, comme suit [42]:
Des SEHs sans stockage d’énergie.
Des SEHs avec stockage d’énergie.
Des SEHs basés sur des sources d’énergie propre seulement.
Des SEHs comportant des sources d’énergie conventionnelles.
Le choix adéquat de la combinaison des ressources d’un SEH est fait en fonction
du potentiel énergétique (solaire, éolien, etc) de l’endroit de l’installation, des besoins
énergétiques des consommateurs (charges à CC et à CA) [42], du champ d’application
17
Chapitre I : Etat de l’art sur les systèmes photovoltaïques hybrides
(a)
18
Chapitre I : Etat de l’art sur les systèmes photovoltaïques hybrides
(b)
(c)
Figure I.15 : Configuration d’un SEH (PV-éolien-stockage) : a) Connexion au bus CC, b) Connexion
au bus CC et CA, c) Connexion au bus CA [40].
19
Chapitre I : Etat de l’art sur les systèmes photovoltaïques hybrides
Contrôleur centralisé
20
Chapitre I : Etat de l’art sur les systèmes photovoltaïques hybrides
Tableau I.3 : Limites de distorsion harmonique pour les systèmes de génération distribuée [9].
21
Chapitre I : Etat de l’art sur les systèmes photovoltaïques hybrides
I.9. Conclusion
Un système PV hybride est généralement constitué d’un champ PV, d’un
générateur éolien, d’une source d’appoint (PAC ou générateur diesel), d’un système
de stockage et d’un ensemble de convertisseurs de puissance. La quantité d’énergie
produite par un SEH à base de SERs (PV et éolienne) dépend des conditions
climatiques à savoir l’irradiation solaire, la température ambiante et la vitesse du vent.
Ce chapitre était consacré à l’introduction de quelques notions de base relatives aux
systèmes de production hybrides. D’abord, nous avons exposé le principe de
fonctionnement de la source PV, du générateur éolien et de la PAC-OS qui sont les
ressources énergétiques étudiées dans cette thèse. Ensuite, nous avons présenté les
différentes configurations des systèmes PV hybrides. Enfin, nous avons recensé les
problématiques liées à la gestion des flux énergétiques dans un SEH et à la qualité
d’énergie électrique fournie aux utilisateurs finaux.
Dans le chapitre suivant, nous nous intéresserons à la modélisation des sources de
production d’énergie et aux TIAs.
22
Chapitre II
II.1. Introduction
En effet, l’élaboration d’une loi de commande classique pour un système exige le
développement de son modèle mathématique. Ce dernier permet d’émuler le
comportement du système commandé d’une manière plus ou moins précise. La
complexité de la modélisation dépend des phénomènes et des dynamiques du système
à considérer et du degré de précision désirée [61]. Le modèle mathématique d’une
source énergétique permet alors d’examiner et d’analyser son comportement pour
différents modes d’opération. En outre, la modélisation des convertisseurs
d’électronique de puissance servit à leur simulation et à la synthèse de leurs
contrôleurs.
Ce chapitre est scindé en deux parties : la première partie est destinée à la
modélisation de différents types de sources d’énergie (PV, éolienne, PAC-OS) ainsi
que des éléments des chaines de conversion énergétique étudiées dans cette thèse.
Tandis que, la deuxième partie de ce chapitre se veut une introduction aux techniques
de l’intelligence artificielle à savoir la logique floue, les réseaux de neurones artificiels
et les réseaux adaptatifs à inférence floue. Cette synthèse a pour ambition de montrer
l’importance de ces TIAs pour la problématique de commande des SEHs dynamiques
et non-linéaires.
-
Figure II.1 : Circuit équivalant simplifié de la cellule PV.
23
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
Le circuit équivalent de ce modèle est illustré par la figure (II.1). Il comprend une
source du courant, une diode D et une résistance série RS. La caractéristique non
linéaire courant-tension d’une cellule PV peut être décrite par l’équation suivante :
𝑞
𝐼𝑃𝑉 = 𝐼𝑃𝐻 − 𝐼𝑠 [𝑒𝑥𝑝 (𝐾 𝑇 𝑎 (𝑉𝑃𝑉 + 𝐼𝑃𝑉 𝑅𝑠 )) − 1] (II.1)
Tel que
o IPV est le courant de sortie de la cellule PV.
o VPV est la tension aux bornes de la cellule PV.
o IPH est le photo-courant proportionnel au rayonnement solaire.
o Is est le courant de saturation de la diode.
o q est la charge de l’électron (q=1.6 × 10-19 C).
o K est la constante de Boltzmann (K=1.3806503 × 10-23 J/K).
o T est la température de la cellule.
o a est le facteur d’idéalité de la jonction PN.
o Rs est la résistance série.
Le photo-courant qui dépend de la température peut être exprimé comme suit :
𝐼𝑃𝐻 = 𝐼𝑃𝐻 (𝑇1 )(1 + 𝐾0 (𝑇 − 𝑇1 )) (II.2)
T1 est la température ambiante nominale qui est égale à 25°C. K0 est le coefficient de
température du courant de court-circuit. La relation linéaire entre le photo-courant et
l’irradiation solaire (G) peut être décrite par l’équation suivante :
𝐺
𝐼𝑃𝐻 (𝑇1 ) = 𝐼𝐶𝐶 (𝑇1,𝑛𝑜𝑚 ) 𝐺 (II.3)
𝑛𝑜𝑚
Où Gnom=1000W/m2 est l’irradiation nominale. ICC est le courant de court-circuit de
la cellule PV. Le coefficient K0 est déterminé par :
𝐼𝐶𝐶 (𝑇2 )−𝐼𝐶𝐶 (𝑇1 )
𝐾0 = (II.4)
𝑇2 −𝑇1
T2 =75°C. Le courant de saturation qui varie en fonction de la température de la
cellule peut être calculé comme suit :
3
𝑇 𝑎 −𝑞𝐸𝑔 1 1
𝐼𝑠 = 𝐼𝑠 (𝑇1 ) × (𝑇 ) exp ( 𝑎𝐾
( − ))
T T
(II.5)
1 1
24
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
𝑞 𝑉
𝐼𝑃𝑉 = 𝑁𝑃 𝐼𝑃𝐻 − 𝑁𝑃 𝐼𝑆 [𝑒𝑥𝑝 (𝐾𝑇𝐴 ( 𝑁𝑃𝑉 + 𝐼𝑃𝑉 𝑅𝑆 )) − 1] (II.9)
𝑆
Les équations de (II.1) à (II.9) sont donc utilisées pour construire le modèle du
panneau PV dans l’environnement Matlab/Simulink.
II.2.1.2. Simulation du module PV
Dans cette thèse, le module BP MSX120 constitué de 72 cellules en série est choisi
pour l’étude et la simulation. Leur caractéristiques électriques sont données dans
l’annexe A. Le schéma du modèle de GPV implémenté dans l’environnement
Matlab/Simulink est illustré par la figure (II.2). Ce modèle nous permet d’étudier
l’effet de la variation de l’éclairement et de la température sur les caractéristiques
IPV(VPV) et PPV (VPV) de la source PV.
1
T Ipv
Iph
i
2 W +
- 1
G Iph Current
+
Add R +
3 ID 1
1 VC
T
s
-
-
Vpv
T Id
2
IC
calcul Id
5 G=1000W/m2
G=800W/m2 120
G=1000W/m2
4.5 G=600W/m2
G=800W/m2
G=400W/m2
G=600W/m2
4 G=200W/m2 100 G=400W/m2
G=200W/m2
3.5
Puissance PV (W)
Courant PV (A)
80
3
2.5 60
2
40
1.5
1
20
0.5
0 0
0 5 10 15 20 25 30 35 40 0 5 10 15 20 25 30 35 40 45
Tension PV (V) Tension PV (V)
(a) (b)
Figure II.3 : Caractéristiques IPV(VPV) et PPV(VPV) du module BP MSX120 pour différents niveaux
d’éclairement et à une température constante (T= 250C).
140
4
120
T=70°C
3.5 T=50°C
T=30°C
100
3 T=10°C
Puissance PV (W)
T=0°C
Courant PV (A)
2.5 80
T=70°C
T=50°C
2 T=30°C
60
T=10°C
1.5 T=0°C
40
1
20
0.5
0 0
0 5 10 15 20 25 30 35 40 45 0 10 20 30 40 50
Tension PV (V) Tension PV (V)
(a) (b)
Figure II.4 : Caractéristiques IPV(VPV) et PPV(VPV) du module BP MSX120 pour différentes valeurs
de température et à G=1000W/m2.
25
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
La figure (II.5) illustre le modèle d’une TE à axe horizontal (de 10kW) implémenté
dans Matlab/Simulink. La caractéristique puissance-vitesse de rotation de la TE, sous
la vitesse nominale du vent (14 m/s), est illustrée par la figure (II.6.a). La courbe de son
coefficient de puissance versus le rapport de vitesse est représentée dans la figure
(II.6.b). Nous pouvons constater que le PPM de la TE est obtenu à la vitesse de rotation
optimale 𝜔𝑡−𝑜𝑝𝑡 =56.7rad/s (où, l’angle de calage β=0°). Cette valeur correspond ainsi
26
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
3 BETA
BET A Cp
Lamda
2 Di vi de1 Subsystem1 2
R 1.225 Pt
0.5
Product3
RO
Product2 1/2
pi (2)^2
1
pi R1 Avoi d di vi si on
2 u(1)^3 Tt
Di vi de
by zero
Wi nd speed Vwi nd
(m/s) beta^3
1
Generator speed (omega)
10000
Puissance Pt(W)
8000
6000
4000
2000
0
0 10 20 30 40 50 60 70 80 90 100
Vitesse de rotation (rad/s)
(a)
0.5
0.4
0.3
Cp
0.2
0.1
0
0 2 4 6 8 10 12 14
Rapport de vitesse
(b)
Figure II.6 : Courbe de a) Puissance versus vitesse de rotation de la TE, b) Coefficient de puissance.
27
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
28
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
Tel que 𝛟𝒓 est le vecteur du flux induit par les aimants permanents. LSS est la
matrice des inductances du stator. En remplaçant par l’équation (II.19) dans (II.18) , on
obtient le modèle de la MSAP dans le repère fixe (abc) :
𝑑𝒊𝑠 𝑑𝛟𝒓
𝒗𝑠 = 𝑅𝑠 𝒊𝑠 + 𝐿𝑠𝑠 𝑑𝑡
+ 𝑑𝑡
(II.20)
b 𝒅𝛉
d ωe= 𝒅𝒕𝒆
q
𝜃𝑒 a
29
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
Avec N0 est le nombre de cellules élémentaires mises en série dans la pile, E0 est le
potentiel aux conditions standards (température=25°C et pression =1atm) et R est la
constante universelle des gaz. T est la température absolue, F est la constante de
Faraday, Ifc est le courant de sortie de la PAC, r est la résistance interne représentant
30
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
les pertes ohmiques dans le stack. PH2, PO2, PH2O sont, respectivement, les pressions
partielles d’hydrogène, d’oxygène et d’eau.
a. Caractérisation des sorties des canaux
Quand on alimente un orifice considéré étouffé par une mixture de gaz de masse
molaire moyenne M, à une température constante, il va présenter la caractéristique
suivante :
𝑊
𝑃𝑢
= 𝐾√𝑀 (II.30)
Avec W est le débit de masse, K est la constante de la valve qui dépend de la surface
de l’orifice et Pu est la pression à l’intérieur du canal (voir figure (II.8)).
Canal
Valve
M, T, Pu W
Notons que 𝑞𝐻2 est le débit molaire d’hydrogène qui est égal à la dérivée de nH2.
31
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
Selon les relations basiques d’électrochimie, le débit 𝑞𝐻𝑟 2 peut être calculé par :
𝑁0 𝐼𝑓𝑐
𝑞𝐻𝑟 2 = = 2𝐾𝑟 𝐼𝑓𝑐 (II.36)
2𝐹
Kr est une constante définie comme le rapport (N0/4F). En substituant les équations
(II.31) et (II.36) dans l’équation (II.35), on obtient :
𝑅𝑇𝐾𝐻2 1
𝑃̇𝐻2 = 𝑉𝑎𝑛
(
𝐾𝐻2
(𝑞𝐻𝑖𝑛2 − 2𝐾𝑟 𝐼𝑓𝑐 ) − 𝑃𝐻2 ) (II.37)
On arrive finalement à l’expression suivante de pression partielle d’hydrogène :
1 1
𝑃̇𝐻2 = − 𝜏 (𝑃𝐻2 − 𝐾 (𝑞𝐻𝑖𝑛2 − 2𝐾𝑟 𝐼𝑓𝑐 )) (II.38)
𝐻2 𝐻2
Tel que (𝜏𝐻2 = 𝑉𝑎𝑛 ⁄𝑅𝑇𝐾𝐻2 ). On procède de la même manière pour calculer les
pressions partielles PO2 et PH2O. Les équations suivantes sont donc obtenues :
1 2
𝑃𝐻̇ 2 𝑂 = − 𝜏 (𝑃𝐻2𝑜 − 𝐾𝑟 𝐼𝑓𝑐 )
𝐻2 𝑂 𝐾 𝐻2 𝑂
(II.39)
1 1
𝑃𝑂̇ 2 = − 𝜏 (𝑃𝑂2 − (𝑞𝑂𝑖𝑛2 − 𝐾𝑟 𝐼𝑓𝑐 ))
{ 𝑂2 𝐾𝑂2
Tel que 𝑞𝑂𝑖𝑛2 est le débit molaire d’oxygène entrant au stack. Les constantes 𝜏𝑜2 et 𝜏𝐻2 𝑂
sont définies par 𝜏𝑜2 = 𝑉𝑐𝑎𝑡ℎ /𝑅𝑇𝐾𝑂2 et 𝜏𝐻2𝑂 = 𝑉𝑎𝑛 /𝑅𝑇𝐾𝐻2 𝑂 .
𝑰𝒇𝒄
𝟐𝑲𝒓
- 𝑲𝒓 𝒓
𝒒𝒊𝒏
𝑯𝟐 + -
𝒒𝒊𝒏
𝑶𝟐 +
𝟏⁄𝑲𝑯𝟐 𝟏⁄𝑲𝑯𝟐𝒐
𝟏⁄𝑲𝒐𝟐
𝟏 + 𝝉𝑯𝟐 𝒔 𝟏 + 𝝉𝑯𝟐𝒐 𝒔 𝟏 + 𝝉𝒐𝟐 𝒔
𝑷𝑯𝟐 𝑷𝑯𝟐 𝑶 𝑷𝑶𝟐
-
𝑹𝑻 𝑷𝑯 𝑷𝟎.𝟓
𝑶
𝑵𝟎 [𝑬𝟎 + 𝐥𝐧 ( 𝟐 𝟐 )]
𝟐𝑭 𝑷𝑯𝟐𝑶 + 𝑽𝒇𝒄
32
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
La réponse électrique rapide dans la PAC-OS est également modélisée par une
fonction de transfert de premier ordre ayant une constante de temps d’environ
𝑇𝑒 =0.8s [76].
Un indice de performance important a été défini dans [76] comme le facteur
d’utilisation de carburant. Il est égal au rapport entre le débit de carburant qui
entre en réaction et le débit d’hydrogène entrant dans la pile :
𝑖𝑛 𝑜𝑢𝑡 𝑟
𝑞𝐻 −𝑞𝐻 𝑞𝐻
𝑈𝑓 = 2
𝑖𝑛
𝑞𝐻
2
= 𝑖𝑛
𝑞𝐻
2
(II.42)
2 2
Donc, le débit d’hydrogène entrant au stack peut être contrôlé afin d’obtenir la
valeur optimale (85%) du facteur 𝑈𝑓 [76] :
𝑟𝑒𝑓 2𝐾𝑟 2𝐾
𝑞𝐻2 = 𝑈 𝐼𝑓𝑐 = 0.85𝑟 𝐼𝑓𝑐 (II.43)
𝑓,𝑜𝑝𝑡
𝑟𝑒𝑓
𝑞𝐻2 est la valeur de référence du débit d’hydrogène fourni à l’entrée du
processeur de combustible. En régime transitoire, le facteur 𝑈𝑓 doit avoir une
valeur dans la marge [0.8, 0.9]. Par la suite, le courant d’entrée de la PAC-OS doit
être maintenu dans l’intervalle défini comme suit [76]:
𝑖𝑛 𝑖𝑛
0.8 𝑞𝐻 𝑖𝑛 0.9 𝑞𝐻
𝐼𝑓𝑐_𝑚𝑖𝑛 = 2
≤ 𝐼𝑓𝑐 ≤ 𝐼𝑓𝑐_𝑚𝑎𝑥 = 2
(V.44)
2𝐾𝑟 2𝐾𝑟
Afin que la différence des pressions d’hydrogène et d’oxygène ne dépasse pas la
valeur de 4kPa , le rapport entre les débits d’hydrogène et d’oxygène (rH-O) doit
avoir la valeur optimale de 1.145 [76] :
𝑖𝑛
𝑞𝐻
𝑟𝐻_𝑂 = 𝑖𝑛
2
= 1.145 (II.45)
𝑞𝑂 2
33
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
𝒒𝒊𝒏
𝑯𝟐
𝟎. 𝟗
𝟐𝑲𝒓
𝑰𝒊𝒏 Ifc
𝒇𝒄 Limites 𝟏
𝟎. 𝟖
𝟏 + 𝑻𝒆 𝒔
𝒒𝒊𝒏
𝑯𝟐 𝟐𝑲𝒓
𝟏 𝒒𝒊𝒏
𝑯𝟐
Ifc 𝟐𝑲𝒓
𝑼𝒇,𝒐𝒑𝒕 𝒒𝒓𝒆𝒇 𝟏 + 𝑻𝒇 𝒔 𝟏
𝑯𝟐
𝒓𝑯_𝑶 𝒒𝒊𝒏
𝑶𝟐
qref H2
2*Kr 1/Uf
1
Gain1 1/Tf 4
Gain4 Tf
s
qin H2
Subtract5 Gain3 Integrator3
qin-H2 1 Vfc+
1/KH2 1
1/tH2 1
s
Subtract Gain6 PH2
+
2.5
1.5
Ph2
1
Po2
0.5 Ph2o
0
0 50 100 150 200 250 300 350 400 450 500
t(s)
(a)
0.86
0.855
Facteur Uf
0.85
0.845
0.84
0.835
0 50 100 150 200 250 300 350 400 450 500
t(s)
(b)
34
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
335.3 310
335.2
305
Vfc(V)
335.1
300
Ifcmin
286.6 280
Ifc(A)
286.5 275
286.4 270
286.3 265
0 500 1000 1500 0 50 100 150 200 250 300 350 400
t(s) t(s)
(b) (c)
En outre, l’état de charge (SOC : State Of Charge) de batterie est défini comme le
rapport entre la quantité de charge stockée dans la batterie, à l’instant t, et la capacité
nominale :
𝑡
∫0 𝐼𝑏𝑑𝑡
𝑆𝑂𝐶(𝑡) = 100(1 − 𝑄
) (II.47)
Tel que Vb est la tension de sortie de batterie, E0 est la tension constante de batterie,
Rb est la résistance interne, Ib est le courant de sortie, K est la tension de polarisation, Q
est la capacité de batterie (Ah), ∫ 𝐼𝑏 𝑑𝑡 est la charge actuelle de la batterie, A est l’
amplitude de la zone exponentielle, B est la constante inverse de la zone exponentielle
de temps. La batterie est dit complètement déchargée ou chargée, si la valeur de SOC
est, respectivement, 0 ou 100%. En conséquence, les processus de charge/décharge
d’une batterie doivent être soigneusement contrôlé afin d’éviter la surcharge et la
décharge profonde.
Rb Ib
+
E + Source de
- tension Vb
-
𝒕
𝑸 it
𝑬 = 𝑬𝟎 − 𝑲 + 𝑨𝒆𝒙𝒑(−𝑩. 𝒊𝒕) ⬚
𝑸 − 𝒊𝒕
𝟎
35
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
V1 C1 Q1 C2 V2
V1 C1 C2 V2 V1 C1 C2 V2
(a) (b)
Figure II.16 : Circuit équivalent du convertisseur boost avec : a) Q1 passant, b) Q1 bloquant [79].
Temps(s)
36
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
37
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
i2
+
Q1
L
i1 C2 V2
V1
Q2 -
Figure II.18 : Circuit de puissance d’un convertisseur réversible abaisseur-élévateur [60].
S1 S2 S3
𝑖𝑔𝑎
a vga
VDC CDC b vgb
n
vao c
vbo 𝑅 𝐿 vga
S4 S5 S6
vco
o
Figure II.19 : Circuit de puissance de l’onduleur à deux niveaux couplé au réseau [72].
38
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
39
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
40
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
Tableau II.1 : Etats de commutation de chaque phase (y=a,b,c) d’onduleur 3L-NPC [72].
La tension de sortie de chaque bras d’onduleur, 𝑣𝑦𝑁 (y = a, b, c), peut être calculée
en fonction de l’état de commutation Sy d’onduleur comme suit [85]:
𝑉𝐷𝐶 𝑉𝐷𝐶 𝑉𝐷𝐶
𝑣𝑎𝑁 = 2
𝑆𝑎 , 𝑣𝑏𝑁 = 2
𝑆𝑏 , 𝑣𝑐𝑁 = 2
𝑆𝑐 (II.64)
Donc, les composantes αβ du vecteur de tension de sortie de l’onduleur 3L-NPC
sont exprimées comme suit :
𝑉𝐷𝐶
𝑣𝛼 = (2𝑆𝑎 − 𝑆𝑏 − 𝑆𝑐 )
6
{ 𝑉𝐷𝐶 (II.65)
𝑣𝛽 = (𝑆𝑏 − 𝑆𝑐 )
2√3
Da Db Dc
a ica
VDC N b c
icb icc
Sa1’ Sb1’ Sc1’
C2 Da ’ Db’ Dc’
Sa2’ Sb2’ Sc2’
41
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
En effet, la commande d’un système a pour but de lui rendre capable de répondre
dynamiquement à une commande. Une tâche de commande d’un système peut être
de régulation ou de poursuite de la trajectoire d’une consigne imposée. Pour achever
la régulation d’un système, on doit lui forcer de rester dans une situation bien précise
même en présence des perturbations. Tandis que, la tâche de poursuite consiste à
contraindre le système de suivre fidèlement la trajectoire d’une référence donnée quel
que soit le type de perturbation [86].
D’ailleurs, la commande d’un SEH dynamique n’est pas une tâche triviale. Les
schémas de commande classiques peuvent mener à des mauvaises performances de
régulation et de poursuite en présence des non-linéarités, des incertitudes et des
variations paramétriques. C’est dans ce contexte que les TIAs ont été reconnues
comme des outils très adéquats pour la commande des SEHs. Dans ce qui suit, on
donne un bref aperçu sur les TIAs qu’on va les appliquer dans cette thèse.
II.3.1. Logique floue
Les systèmes d’inférence floue permettent d’imiter le raisonnement chez l’être
humain qui est capable de prendre des bonnes décisions même en présence de
l’incertitude. Les définitions et les concepts de base de la LF introduits dans cette
section sont compilés des références [87-94].
II.3.1.1. Ensemble flou
En 1965, le Professeur Lotfi A. Zadeh [87] a introduit pour la première fois le
concept d’ensemble flou qui présente une extension de la notion d’ensemble classique.
Dans la théorie d’ensemble flou, un ensemble flou A, dans l’univers de discours B, est
caractérisé par une fonction dite Fonction d’Appartenance (FA) qui assigne à chaque
élément b Є B une valeur dans l’intervalle [0,1] et non pas dans l’ensemble de deux
éléments {0,1} (voir figure (II.21)). L’ensemble A défini à l’aide d’une FA est alors
appelé ensemble flou.
A B
B
b3 b1
b4 b2
42
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
(a) (b)
(c) (d)
Figure II.22 : Exemples de types de FA : (a) Fonction-L; (b) Fonction-Λ (triangulaire); (c) Fonction
Gaussienne; (d) Fonction-Π (trapézoïdale) .
L M R
Degré d’appartenance
1
0.8
0.6
0.4
0.2
0
10 20 30 40 50
Vitesse de rotation (rad/s)
43
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
Moteur
d’inférence
Base
des règles
Un CLF est un système qui permet de faire un passage non-linéaire des grandeurs
de contrôle en entrée vers les variables de sortie en se basant sur les expériences et les
perceptions humaines.
Module de Fuzzification : Cette unité permet de convertir à l’aide des FAs les
valeurs numériques, appartenant aux univers de discours des entrées du CLF,
à des valeurs floues compatibles avec les règles d’inférence floue. Le degré
d’appartenance à chaque ensemble flou associé à chacune des entrées est
déterminé dans cette étape.
Base des règles : Cette base détermine la stratégie de commande adoptée par
un CLF permettant de relier les entrées aux sorties. Les règles floues ayant la
forme d’implication « Si-Alors » sont établies en se basant sur les connaissances
des experts sur le système commandé. Ces règles « Si-Alors » sont interprétées
à l’aide d’un opérateur d’implication flou tel que, l’implication de Mamdani.
A titre d’exemple, le tableau (II.2) donne la base des règles d’un régulateur
PI flou proposée par MacVicar-Whelan en 1976. Ce CLF prend en entrée l’erreur
(e) et la variation d’erreur (Δe) et renvoie en sortie l’action de commande (Δu).
Les termes linguistiques M+(Moyen positif), S+(Petit positif), M-(Moyen
négatif), S-(Petit négatif) et O (Zéro) sont utilisés pour qualifier les ensembles
floues des variables (e, Δe et Δu). Pour chacune des valeurs floues des entrées e
et Δe, l’incrément de commande en sortie (Δu) trouve sa valeur dans la base
des règles du tableau (II.2). Par exemple, la première règle floue est donnée
comme suit :
R1 : Si e est M+ et Δe est M- , Alors Δu est O
Partie prémisse ou antécédent Partie conséquent ou conclusion
Δe
M- S- O S+ M+
e
M+ O S- M- M- L-
S+ S+ O S- S- M-
O M+ S+ O S- M-
S- M+ S+ S+ O S-
M- L+ M+ M+ S+ O
44
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
45
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
𝑧0 = 1 𝑤0
𝑧1 𝑤1
𝑦
𝑧2
𝑤2 ∑ 𝑓
⋮ 𝑤𝑁
𝑧𝑁
46
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
Tableau II.4 : Quelques exemples de fonction d’activation d’un neurone formel [95, 97].
47
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
même couche d’un perceptron ne sont pas interconnectés entre eux. En outre, le
neurone de chaque couche possède des entrées provenant de la couche en amont et
une sortie renvoyée à la couche en aval. Si chaque neurone est interconnecté à chaque
neurone de la couche adjacente, il est dit totalement connecté, sinon il est dénommé
partiellement connecté [92, 100].
La sortie d’un neurone j dans la couche L d’un perceptron multicouche se calcule
comme suit [96]:
(𝐿) (𝐿) (𝐿) (𝐿) 𝑝 (𝐿) (𝐿−1) (𝐿)
𝑦𝑗 = 𝑓𝑗 (𝑛𝑒𝑡𝑗 ) = 𝑓𝑗 ((∑𝑖=1
𝐿−1
𝑤𝑖𝑗 𝑦𝑖 ) + 𝑤0𝑗 ) 𝐿 = 1, … , 𝑀 (II.68)
Sous forme matricielle :
(𝐿)
𝒚(𝐿) = 𝑓 (𝐿) (𝑛𝑒𝑡 (𝐿) ) = 𝑓 (𝐿) (𝑊 (𝐿) 𝒚(𝐿−1) + 𝑊𝑜 ) (II.69)
Notons que : 𝑊𝑜(𝐿) = {𝑤0𝑗
(𝐿)
}, 𝑊 (𝐿) = {𝑤𝑖𝑗(𝐿) }, 𝑖 = 1, … , 𝑝𝐿−1 , 𝑗 = 1, … , 𝑝𝐿 , 𝑓 (𝐿) (𝑛𝑒𝑡 (𝐿) ) =
𝑇
[𝑓1(𝐿) (𝑛𝑒𝑡1(𝐿) ), 𝑓2(𝐿) (𝑛𝑒𝑡2(𝐿) ), … , 𝑓𝑗(𝐿) (𝑛𝑒𝑡𝑝(𝐿)
𝐿
)]
(𝐿) (𝐿) (𝐿)
Tel que𝑓𝑗 , 𝑛𝑒𝑡𝑗 et 𝑤0𝑗 désignent, respectivement, la fonction d’activation, le
(𝐿)
niveau d’activation et le biais du neurone j de la Lème couche. 𝑤𝑖𝑗 représente le poids
de connexion entre le neurone i de la couche (L-1) et le neurone j de la couche (L). 𝑝𝐿
est le nombre des neurones dans la Léme couche. M est le nombre des couches du RNA.
(2)
𝑤01
(2) (3) (3)
𝑤11 ∑𝒇 𝑤11 𝑤01
𝑧1 𝑦1
∑𝒇
⋮ ⋮
𝑧𝑝1 ∑𝒇 𝑦𝑝𝑀
⋮
Couche d’entrée (3)
𝑤0𝑝
∑𝒇
Couche de sortie
(2)
𝑤0𝒑𝟐
Couche cachée
48
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
⋮ ⋮
𝑧𝑁 𝒉𝟑(𝒌) 𝑦𝑝
𝒉𝒎 (𝒌)
49
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
(a)
(b)
Figure II.29 : Types d’apprentissage : a) Supervisé, b) Non-supervisé [105].
50
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
51
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
𝑦(𝑘) = 𝑤0 + 𝑤𝑖 𝑧𝑖
𝑧1 𝑤1 𝑂𝑑 (𝑘)
∑
𝑖=1
𝑧2
𝑤2 ∑
- +
⋮ 𝑤𝑁
𝑧𝑁
𝑒𝑜 (𝑘)
Mise à jour des
𝒛(𝑘) poids (LMS)
52
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
53
Chapitre II : Modélisation des systèmes de conversion d’énergie et
aperçu sur les TIAs
La figure (II.31) illustre l’architecture d’un ANFIS à deux entrées (y1 et y2) et une
seule sortie (f). L’ANFIS est composé de cinq couches, tel que les couches deux, trois
et cinq sont constituées de nœuds fixes qui effectuent des opérations mathématiques.
Alors que, la couche d’entrée et la quatrième couche sont composées de nœuds
adaptatifs qui doivent être ajustés durant la phase d’apprentissage. [89]. Ce système
flou se base sur le modèle d’inférence de TSK d’ordre un. Plus de détails sur
l’architecture de l’ANFIS est donné dans l’annexe C.
y1 y1 y2
μA1
A1
w1 w1
y1 П N
y
μA12 w1 . f1
A2
f
y
∑
B1 μB21
y2 П N w 2 . f2
B2 y w2 w2
μB22 Couche 5
y1 y2
Couche 1 Couche 2 Couche 3
Couche 4
Figure II.31 : Architecture d’un ANFIS [89].
Le rôle fondamental d’un algorithme d’apprentissage d’un ANFIS est d’ajuster les
paramètres non-linéaires de la partie prémisse et les paramètres linéaires de la partie
conclusion. Ceci est achevé en utilisant les informations provenant d’une base de
données d’apprentissage [105,89].
D’après [89], l’apprentissage d’un ANFIS peut être accompli soit à l’aide de
l’algorithme de rétropropagation d’erreur, soit via l’algorithme hybride qui combine
la règle de la descente du gradient et la méthode des moindres carrés. En utilisant
l’algorithme hybride, l’identification des paramètres de l’ANFIS se fait en deux étapes
à chaque itération : D’abords, les paramètres conséquents sont déterminés à l’aide de
la méthode des moindres carrés. Puis, les paramètres non linéaires des prémisses sont
ajustés via la règle de la descente du gradient [89].
II.4. Conclusion
Dans ce chapitre, nous nous sommes intéressés au développement des modèles
dynamiques des sources d’énergie (PV, éolienne et PAC) et de différents types de
convertisseurs de puissance. Ces modèles facilitent la tâche de simulation du
comportement des sources énergétiques et permettent d’étudier leurs structures de
commande d’une manière appropriée. En outre, nous avons dressé un état de l’art sur
les techniques d’intelligence artificielle (LF, RNAs, ANFIS). Les notions de base sur les
TIAs ont été introduites afin de les appliquer aux schémas de commande étudiés dans
cette thèse.
Le chapitre suivant est consacré à l’étude et à la simulation de trois systèmes de
production PV, éolienne et à base de PAC-OS.
54
Chapitre III
III.1. Introduction
Les panneaux photovoltaïques et les piles à combustibles sont des sources de
courant électrique continu tandis qu’un générateur éolien fournit un courant alternatif
de fréquence incompatible avec celle du réseau public. De ce fait, l’intégration de ce
type de sources d’énergie dans le réseau requière la conversion et l’adaptation de leur
puissance produite afin d’avoir une puissance cohérente avec le réseau électrique. En
effet, le système de conditionnement de puissance et sa structure de commande
associée représentent les éléments clés dans la phase de conception des systèmes de
production décentralisée.
Dans ce chapitre, une investigation sur les structures de commande dédiées aux
certains systèmes de production connectés au réseau électrique est achevée. Nous
examinons certaines configurations des chaines de conversion d’énergie à base des
sources PV, éolienne et PAC. Des notions de base sur les stratégies de commande
conventionnelles destinées au pilotage des convertisseurs de puissance sont également
introduites. De plus, des TIAs sont appliquées afin d’optimiser le rendement
énergétique et d’améliorer les réponses des systèmes commandés en régimes
transitoire et permanent. Le comportement dynamique des systèmes de génération
étudiés est examiné par simulation via le logiciel de Matlab/Simulink. Les résultats de
simulation obtenus sont présentés afin de vérifier l’efficacité des schémas de
commande et d’évaluer les performances dynamiques des systèmes étudiés.
55
Chapitre III : Commande des systèmes de production connectés au réseau
Réseau
GPV Convertisseur Convertisseur Filtre
VPV VDC
CC-CC CC-CA RL
- -
Spv Sa,b,c
MPPT DPC vga,b,c
IncCond
Figure III.1 : Structure du SPVCR étudié.
56
Chapitre III : Commande des systèmes de production connectés au réseau
𝒅𝑰𝑷𝑽 𝑰𝑷𝑽
=− PPM
𝒅𝑽𝑷𝑽 𝑽𝑷𝑽
𝒅𝑰𝑷𝑽 𝑰𝑷𝑽
<−
𝒅𝑽𝑷𝑽 𝑽𝑷𝑽
PPV (W)
𝒅𝑰𝑷𝑽 𝑰𝑷𝑽
>−
𝒅𝑽𝑷𝑽 𝑽𝑷𝑽
VPV (V)
Figure III.2 : Position du point de fonctionnement sur la courbe PPV(VPV).
Début
Non Oui
ΔVPV=0
Vref (k)=Vref (k-1)+ΔVref Vref (k)=Vref (k-1)-ΔVref Vref (k)=Vref (k-1)-ΔVref Vref (k)=Vref (k-1)+ΔVref
Retour
57
Chapitre III : Commande des systèmes de production connectés au réseau
La structure de DPC [118,122], adoptée dans cette étude, est illustrée par la figure
(III.6). Les comparateurs à hystérésis du schéma DPC classique sont remplacés par des
contrôleurs linéaires de type PI. De plus, un bloc de modulation MLI est utilisé au lieu
du tableau de commutation afin d’assurer une fréquence de commutation constante.
L’incorporation d’une boucle à verrouillage de phase (PLL : Phase Locked Loop) dans
cette structure est nécessaire afin d’assurer la synchronisation avec le réseau. Le PLL
calcule l’angle du vecteur de tension du réseau (g) qui est utilisé pour la transformation
58
Chapitre III : Commande des systèmes de production connectés au réseau
R L
Réseau
VDC CDC
VDC
VDCref - + sa,b,c iga,b,c vga,b,c
SPWM
Calcul des puissances
algorithme
PI dq/abc g vga,b,c
KI Qref PLL
vd* vq*
+ Q
PI PI
-
Pref P
+ -
Figure III.6. Schéma de DPC adopté pour l’onduleur côté réseau.
59
Chapitre III : Commande des systèmes de production connectés au réseau
Comme (𝑣𝑔𝑞 = 0 ), les puissances active et réactive varient donc en fonction des
composantes directe (igd) et en quadrature (igq) du courant, respectivement, selon les
équations suivantes :
3
𝑃 = 𝑣𝑔𝑑 𝑖𝑔𝑑
2
{ 3 (III.9)
𝑄 = − 𝑣𝑔𝑑 𝑖𝑔𝑞
2
Le schéma simplifié de la boucle de régulation de la puissance active est illustré par
la figure (III.7).
𝑃𝑟𝑒𝑓 𝑲𝒊 𝟏 𝑖𝑔𝑑 𝟑
+ 𝒗 𝑃
+ 𝑲𝒑 + 𝟐 𝒈𝒅
- 𝒔 𝑳. 𝒔 + 𝑹
- 𝑉
𝑔
Tel que, la constante de temps du système en boucle fermée Tbf est définie par :
𝐿
𝑇𝑏𝑓 = 𝑉 (III.15)
𝑀 𝐾𝑝
60
Chapitre III : Commande des systèmes de production connectés au réseau
Nouvelle position 𝑥𝑖 (𝑘 + 1)
Position courante
𝑥𝑖 (𝑘) Meilleur position globale 𝐺𝑏 (𝑘)
Figure III.8 : Déplacement d’une particule de l’essaim dans l’espace de recherche [125].
Ces dernières années, l’algorithme de PSO est largement utilisé dans divers champs
d’application [125-128], vue les avantages qu’il présente à savoir la simplicité de
l’algorithme, l’efficacité de calcul et la facilité d’implémentation. En automatique,
l’algorithme de PSO est fréquemment utilisé pour l’optimisation des paramètres des
contrôleurs classiques ou même pour la substitution des algorithmes de contrôle
traditionnels permettant d’améliorer les performances dynamiques et statiques des
systèmes commandés.
A. Formulation
Etant donné un espace de recherche de dimension D, chaque particule i d’une
population de N individus est caractérisée par son vecteur position 𝑥𝑖 =
(𝑥𝑖1 , 𝑥𝑖2 , … , 𝑥𝑖𝐷 ) qui représente une solution candidate. Le taux de déplacement de la
particule i est défini par un vecteur vitesse 𝑣𝑖 = (𝑣𝑖1 , 𝑣𝑖2 , … , 𝑣𝑖𝐷 ) [125,129]. Par ailleurs,
61
Chapitre III : Commande des systèmes de production connectés au réseau
62
Chapitre III : Commande des systèmes de production connectés au réseau
Début
Non
Fonction coût actuelle < Meilleur coût local ?
Oui
Mettre meilleur coût local = Coût actuel
Non
Coût actuel < Meilleur coût global ?
Oui
Mettre meilleur coût global = Coût actuel
63
Chapitre III : Commande des systèmes de production connectés au réseau
64
Chapitre III : Commande des systèmes de production connectés au réseau
4
x 10
active power(W)
1
-1
0 0.2 0.4 0.6 0.8 1 1.2
t(s)
4
reactive power(var)
x 10
1
-1
0 0.2 0.4 0.6 0.8 1 1.2
t(s)
Figure III.10 : Variation des puissances active et réactive sous un niveau d’éclairement variable.
500
400
DC-link voltage (V)
300
200
100
0
0 0.2 0.4 0.6 0.8 1 1.2
t(s)
10
Inverter output current (A)
-5 10
Inverter output current (A)
-10
0 0.2 0.4 0.6 -5 0.8 1 1.2
t (s)
-10
0.25 0.26 0.27 0.28 0.29 0.3 0.31 0.32 0.33 0.34 0.35
t(s)
65
Chapitre III : Commande des systèmes de production connectés au réseau
R L vga
CDC vgb
GSAP
vgc
Sa,b,c Sa,b,c
isabc , vsabc VDC
iga,b,c
MPPT + Schéma de
Pt , Te , ωt DTC VOC vga,b,c
Figure III.13 : Topologie du système de production éolien étudié.
66
Chapitre III : Commande des systèmes de production connectés au réseau
𝝓* + 𝑆𝝓
Sa
𝝓 - Tableau de Sb
commutation Sc
𝝎∗𝒕 + Te*+ 𝑆𝑇
PI abc i
- - sabc
Te s isα
𝝎𝒕
Estimation du αβ
flux et du couple isβ GSAP
VDC
Dans cette section, le schéma modifié de DTC-SVM, basé sur des contrôleurs PI, est
adopté pour la commande du redresseur interfaçant la GSAP. Où, l’algorithme de PSO
est utilisé pour déterminer hors-ligne les paramètres des régulateurs PI du flux et du
couple électromagnétique. La figure (III.15) montre le schéma de DTC adopté.
Estimation du couple et du flux
Dans le repère de stationnaire αβ, les composantes α et β du flux statorique de la
GSAP sont données comme suit [139]:
𝑡
ϕ𝑠𝛼 = ∫0 (𝑣𝑠𝛼 − 𝑅𝑠 𝑖𝑠𝛼 )𝑑𝑡 + 𝜙𝑟
{ 𝑡 (III.24)
ϕ𝑠𝛽 = ∫0 (𝑣𝑠𝛽 − 𝑅𝑠 𝑖𝑠𝛽 ) 𝑑𝑡
67
Chapitre III : Commande des systèmes de production connectés au réseau
Comme le montre le schéma de la figure (III.15), les erreurs des modules du flux et
du couple sont régulées à l’aide des contrôleurs PI. Ensuite, les composantes vsd* et vsq*
du vecteur de tension de référence sont transformées en coordonnées αβ. Ces
composantes sont envoyées au modulateur SVM afin de générer les signaux de
commande du redresseur. On note que la structure de DTC adoptée performe un
contrôle découplé du flux et du couple sans utilisation de termes de découplage.
III.3.2.2. Poursuite du PPM du GE
Suite à la nature fluctuante du vent, la puissance mécanique générée par la TE varie
non-linéairement dans le temps. La figure (III.16) montre les caractéristiques
puissance-vitesse de rotation d’une éolienne à axe horizontal pour différentes vitesses
de vent (où β=0°). On peut observer que pour chaque vitesse de vent, la courbe de
puissance possède un unique PPM. Le maximum de puissance est obtenu lorsque le
GE tourne avec une certaine vitesse de rotation.
Comme déjà vu, l’implémentation d’un contrôleur MPPT pour la commande d’une
éolienne à vitesse variable (dans la gamme des vitesses moyennes) est indispensable.
Il permet de forcer le GE de fonctionner en permanence au PPM en agissant souvent
sur sa vitesse de rotation. Sans tenir compte de la vitesse de vent, chaque vitesse de
rotation optimale (𝜔𝑡−𝑜𝑝𝑡 ) de l’éolienne permet de maintenir le rapport de vitesse à sa
valeur optimale unique ( 𝜆𝑜𝑝𝑡 ), qui correspond ainsi au coefficient de puissance
maximal ( 𝐶𝑃𝑚𝑎𝑥 ) [25,69] (voir figure (III.16)). Les techniques de MPPT les plus
68
Chapitre III : Commande des systèmes de production connectés au réseau
9000 13m/s
8000
12m/s
Puissance Pt(W)
7000
6000
11m/s
5000
10m/s
4000
3000
8m/s
2000
1000
0
0 10 20 30 40 50 60 70 80 90 100
Vitesse de rotation (rad/s)
Figure III.16 : Courbes de puissance Pt(ωt) d’une TE pour différentes vitesses de vent.
69
Chapitre III : Commande des systèmes de production connectés au réseau
𝝈𝝎 (𝒌) Règle
d’apprentissage
Couche cachée
Couche d’entrée (𝟐)
𝒈𝟏 Couche de sortie
∆𝑷𝒕 (𝒌) (𝟐)
𝒈𝟐
𝑻∗𝒆
∆𝝎𝒕 (𝒌) 𝑤0𝑛
𝑤𝑗𝑛
(𝟐)
𝒈𝟑
𝑤𝑖𝑗𝑛 1
Tel que 𝑑𝑒𝑤 (𝑘) = 𝑒𝑤 (𝑘) − 𝑒𝑤 (𝑘 − 1) et 𝜆𝜔 est une constante positive. Les dérivées
partielles de la fonction coût par rapport aux poids du perceptron sont calculées
comme suit (voir annexe B) :
2
𝜕𝜎𝜔 (𝑘) 𝜕𝑆𝜔 (𝑘) (1)
𝑛
𝜕𝑤𝑖𝑗
= 𝜕𝑦 (3)
𝑆𝜔 (𝑘)𝑤𝑗𝑛 [1 − (𝑦(𝑗 2) ) ] 𝑦𝑖
𝜕𝜎𝜔 (𝑘) 𝜕𝑆𝜔 (𝑘)
𝜕𝑤0𝑛
= 𝜕𝑦 (3)
𝑆𝜔 (𝑘) j=1,2,3, i=1,2 (III.31)
𝜕𝜎𝜔 (𝑘) 𝜕𝑆𝜔 (𝑘) (2)
= 𝑆𝜔 (𝑘)𝑦𝑗
{ 𝜕𝑤𝑗𝑛 𝜕𝑦 (3)
70
Chapitre III : Commande des systèmes de production connectés au réseau
On remplace (III.31) dans (III.29) en approximant le terme 𝜕𝑆𝜔 (𝑘)/𝜕𝑦 (3)à un. La
mise à jour des paramètres du FFNN se fait alors selon les règles suivantes :
2
𝑛( 𝑛( ) (1)
𝑤𝑖𝑗 𝑘 + 1) = 𝑤𝑖𝑗 𝑘 − 𝜂𝑆𝜔 (𝑘)𝑤𝑗𝑛 [1 − (𝑦(𝑗 2) ) ] 𝑦𝑖
𝑤0𝑛 (𝑘 + 1) = 𝑤0𝑛 (𝑘) − 𝜂𝑆𝜔 (𝑘) j=1,2,3, i=1,2 (III.32)
(2)
{ 𝑤0𝑛 (𝑘 + 1) = 𝑤0𝑛 (𝑘) − 𝜂𝑆𝜔 (𝑘)𝑦𝑗
Par rapport à la technique de contrôle à TSR optimal, cette méthode de MPPT à
base de FFNN adaptatif ne nécessite pas la connaissance des courbes caractéristiques
de l’éolienne. Elle requiert seulement le calcul des erreurs de puissance et de vitesse
de la TE.
III.3.3. Commande à tension orientée de l’onduleur
Le schéma de commande vectorielle VOC, adopté pour le contrôle de l’onduleur
connecté au réseau, est illustré par la figure (III.19). Ce schéma est établi dans le
référentiel synchrone dq, tournant avec une vitesse 𝜔, et orienté de sorte que l’axe d est
aligné avec le vecteur de tension du réseau (voir figure (III.20))[141]. Le principe de la
technique de VOC est basé sur la transformation des grandeurs contrôlées entre le
repère naturel abc et le repère synchrone dq. En régime permanent, les variables de
contrôle deviennent des grandeurs continues dans le référentiel tournant [9,69]. On
note ici que la vitesse de rotation du repère dq est égale à la pulsation de la tension du
réseau 𝜔 = 2𝜋𝑓𝑟é𝑠 , où 𝑓𝑟é𝑠 est la fréquence du réseau.
Comme le montre la figure (III.19), le schéma de VOC est constitué de deux boucles
en cascade de contrôle de tension du bus CC et du courant de réseau. L’objectif de la
commande vectorielle de l’onduleur est de contrôler le courant triphasé de sortie dans
le but de générer des formes d’ondes sinusoïdales (à faible THD ) ayant une fréquence
égale à celle de la tension du réseau. En outre, elle doit maintenir la tension du bus
continu constante en dépit du changement des flux des puissances active et réactive
[9,21].
VDC vgd
PSO
VDCref + - igd* +
PI PI vd*
-
igd + + dq
vga,b,c g ωL Sa
PLL SPWM
Sb
igq ωL
abc Sc
g g igq*=0 + -
PI vq*
abc vgd abc igd
vgabc + +
PSO
dq vgq igabc dq vgq
igq
71
Chapitre III : Commande des systèmes de production connectés au réseau
β
q
𝑖𝑔
𝑖𝑔𝛽 d
ω
𝑣𝑔
𝑖𝑔𝑞 𝑖𝑔𝑑
𝜃𝑔 α
𝑖𝑔𝛼
Figure III.20 : Orientation du repère dq selon le vecteur de tension du réseau [141].
En effet, les boucles de contrôle des composantes d et q du courant réseau (igd et igq)
sont responsables de la régulation indirecte des puissances active et réactive et de
l’amélioration de la qualité du courant du réseau. Alors que, la boucle externe de
contrôle de tension (VDC) se charge de la stabilisation de la tension du bus CC en
ajustant la référence du courant direct (igd*). Pour assurer le fonctionnement à un
facteur de puissance unitaire, la consigne du courant en quadrature est fixée à zéro
(igq*=0) [9,69].
Comme le montre la figure (III.19), les tensions des trois phases (a,b,c) du réseau
sont mesurées et envoyées à la PLL afin de déterminer l’angle du vecteur de tension
du réseau g utilisé par les transformées de Park directe et inverse. En outre, le vecteur
du courant réseau dans le repère fixe (abc) est transformé à un vecteur biphasé à deux
composantes igd et igq dans le repère tournant dq. Des contrôleurs de type PI s’occupent
par la suite de la régulation des composantes du courant mesuré (igd et igq) pour qu’elles
suivent les trajectoires de leurs références correspondantes (igd* et igq*).
D’après l’équation (II.62), les dérivées des courants direct et en quadrature peuvent
être exprimées comme suit :
𝑑𝑖𝑔𝑑 1
= 𝐿 (−𝑅𝑖𝑔𝑑 − 𝐿𝜔𝑖𝑔𝑞 − 𝑣𝑔𝑑 + 𝑣𝑑 )
{ 𝑑𝑖𝑑𝑡 1
(III.33)
𝑔𝑞
𝑑𝑡
= 𝐿 (−𝑅𝑖𝑔𝑞 + 𝐿𝜔𝑖𝑔𝑑 − 𝑣𝑔𝑞 + 𝑣𝑞 )
Afin d’achever un contrôle découplé des courants igd et igq et d’améliorer les
performances dynamiques du schéma VOC, les termes de découplage (Lωigd et Lωigq)
et les tensions de compensation 𝑣𝑔𝑑 et 𝑣𝑔𝑞 sont pris en compte comme le montre la
figure (III.19). Les tensions de référence 𝑣𝑑∗ et 𝑣𝑞∗ à l’entrée du générateur SPWM sont
calculées comme suit [69]:
𝑣𝑑∗ = 𝑣𝑖𝑑 + 𝐿𝜔𝑖𝑔𝑞 + 𝑣𝑔𝑑
{ ∗ (III.34)
𝑣𝑞 = 𝑣𝑖𝑞 − 𝐿𝜔𝑖𝑔𝑑 + 𝑣𝑔𝑞
Tel que 𝑣𝑖𝑑 et 𝑣𝑖𝑞 sont les tensions de sortie des contrôleurs des courants selon les
axes d et q, respectivement. Le modulateur SPWM s’occupe ensuite de la génération
des signaux de commande des interrupteurs de puissance de l’onduleur.
III.3.3.1. Contrôle de la tension du bus continu
Le bus continu subit des variations aléatoires de l’énergie issue du redresseur. En
outre, les perturbations du réseau peuvent affectées considérablement la tension du
bus CC. Donc, elle doit être maintenue constante afin d’assurer l’équilibre d’énergie
72
Chapitre III : Commande des systèmes de production connectés au réseau
CDC représente la capacité du bus CC. Cette fonction de transfert a la même forme
canonique suivante [61]:
1+𝑇𝑛 𝑠
𝑇𝑉 (𝑠) = 1 2 2𝜉 (III.36)
𝑠 + 𝑠+1
𝜔2𝑛 𝜔𝑛
73
Chapitre III : Commande des systèmes de production connectés au réseau
74
Chapitre III : Commande des systèmes de production connectés au réseau
Paramètre Valeur
Nombre maximal d’itération Kmax=35.
Nombre des particules d’essaim N=6.
Les valeurs maximale et minimale des C1min=C2min=1.
coefficients cognitif et social C1max=C2max=2
Coefficients d’inertie maximal et minimal wmax=0.9, wmin=0.4
Paramètres des contrôleurs du couple et 𝐾𝑝𝑇 = 15.5, 𝐾𝑖𝑇 = 600, 𝐾𝑝𝑑 =
du courant direct. 690.2872, 𝐾𝑖𝑑 = 2
Dans le premier cas, nous avons testé le SGEE sous la vitesse nominale du vent
(VW=14m/s). La figure (III.22) montre la variation de la puissance mécanique et de la
vitesse rotationnelle de la TE. L’évolution du couple de la TE et du flux statorique est
décrite par la figure (III.23). Les composantes d et q du courant du réseau, la tension
du bus continu et les puissances active et réactive sont représentés dans la figure
(III.24). Il est clair de la figure (III.22) que le GE a atteint rapidement le point
d’opération optimal après un temps de réponse de l’ordre de 0.02s. En comparaison
avec le contrôleur par TSR optimal, le FFNN offre un dépassement minimal et une
réponse transitoire plus rapide. Les régulateurs PI ajustés à l’aide de l’algorithme de
PSO ont garanti des réponses transitoire et statique satisfaisantes. Nous pouvons
constater de la figure (III.23.a) que l’indice de performance choisi (ITSE) permet
d’avoir des meilleures performances de contrôle que le critère ISE (Integral Square
Error). En outre, la figure (III.24) montre que les objectifs de contrôle de l’onduleur
triphasé sont achevés à l’aide du schéma de VOC. Les composantes d et q du courant
réseau suivent parfaitement leurs consignes qui sont atteintes après un temps de
réponse d’environ 0.032s. Le fonctionnement à un facteur de puissance unitaire est
donc assuré en maintenant la puissance réactive constamment nulle. De plus, la
puissance générée par le GE est totalement injectée au réseau électrique. Le régulateur
PI de tension du bus CC offre ainsi des performances de contrôle satisfaisantes : un
dépassement maximal de l’ordre de 1031V et un temps de réponse d’environ 0.35s.
10000
Pussance de la TE (W)
8000
6000
4000
2000
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
t(s)
(a)
100
vitesse de rotation de la TE (rad/s)
PI
80
FFNN
60
40
20
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
t(s)
(b)
75
Chapitre III : Commande des systèmes de production connectés au réseau
200
Couple de la TE (N.m)
150
100 ISE
ITSE
50
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
t(s)
(a)
0.09
0.085
Module du flux (Wb)
0.08
flux
0.075 flux*
0.07
0.065
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
t(s)
(b)
Figure III.23 : Evolution du : a) Couple de la TE, b) Module du flux dans le premier cas de
simulation.
50
1000
composante d du courant (A)
0 800
-50 Id 600
Idref
400 Vdcref
-100
Vdc
200
-150
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
t(s) t(s)
(a) x 10
4 (c)
30 3
Iq
composante q du courant (A)
25 2
Puissance P (W) e Q (VAr)
Iqref
20 1
15
0
10
-1
P
5
Q
-2
0
-3
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
t(s) t(s)
(b) (d)
76
Chapitre III : Commande des systèmes de production connectés au réseau
150
couple Te et Tt (Nm)
100
Te
Te*
50 Tt
0
0 0.5 1 1.5 2
t(s)
(a)
vitesse de rotation wt(rad/s)
60
50
40
30 TSR-optimal
20 FFNN
10
0
0 0.5 1 1.5 2
t(s)
(b)
Idref
40 Id
Iq
Iqref
20
-20
0 0.5 1 1.5 2
t(s)
(a)
1000 Vdcref
Vdc
Tension du bus CC (V)
800
600
400
200
0
0 0.5 1 1.5 2
t(s)
(b)
77
Chapitre III : Commande des systèmes de production connectés au réseau
78
Chapitre III : Commande des systèmes de production connectés au réseau
la tension du bus continu est basée sur un contrôleur flou. Tandis que, la boucle interne
de régulation du courant de sortie de la PAC est basée sur un contrôleur PI classique.
Ce dernier génère le rapport cyclique (D) utilisé pour la génération du signal de
commande MLI du convertisseur boost.
VDCref + Ifcref + D Générateur
Sfc
CLF PI MLI
-
- 𝑑
𝑑𝑡 Ifc
VDC
Figure III.28 : Schéma de contrôle du convertisseur boost.
0.8
0.6
0.4
0.2
0
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1
Figure III.29 : FAs des entrées et de la sortie du CLF de la tension du bus CC.
𝒆𝑽𝒅𝒄
NB NM NS ZE PS PM PB
∆𝒆𝑽𝒅𝒄
NB NB NB NB NB NM NS ZE
NM NB NM NB NM NS ZE PS
NS NB NB NS ZE ZE PS PM
ZE NB NM ZE ZE ZE PM PB
PS NM NS ZE ZE PS PB PB
PM NS ZE PS PM PB PM PB
PB ZE PS PM PB PB PB PB
79
Chapitre III : Commande des systèmes de production connectés au réseau
Pref + idref + + +
CLF PI
- - +
P id vcd *
ωL dq /
abc PWM
ωL
iq vcq*
Qref + - + -
CLF PI 𝜽𝒈
- iqref + +
Q vgq
𝜽𝒈 vgd
vga,b,c dq
PLL /abc vgq
Les puissances active et réactive instantanées sont calculées en fonction des tensions
des phases du réseau et des courants mesurés au PCC (𝑖𝑐𝑎 ,𝑖𝑐𝑏 ,𝑖𝑐𝑐 ) comme suit :
𝑃 = 𝑣𝑔𝑎 𝑖𝑐𝑎 + 𝑣𝑔𝑏 𝑖𝑐𝑏 + 𝑣𝑔𝑐 𝑖𝑐𝑐
{ 1 (III.44)
𝑄= [(𝑣𝑔𝑏 − 𝑣𝑔𝑐 )𝑖𝑎 + (𝑣𝑔𝑐 − 𝑣𝑔𝑎 )𝑖𝑏 + (𝑣𝑔𝑎 − 𝑣𝑔𝑏 )𝑖𝑐 ]
√3
Dans cette étude, les régulateurs traditionnels PI des puissances active et réactive
sont remplacés par des CLFs en raison de leurs propriétés avantageuses. L’objectif de
commande ici est de maintenir des bonnes performances de contrôle et d’assurer une
haute qualité d’énergie sous des conditions variables du réseau et de demande
d’énergie. Les CLFs proposés sont basés sur le modèle d’inférence de TSK d’ordre zéro
(voir section II.3.1.4.). Dans ce cas, l’inférence par Somme-Produit et le calcul de sortie
par la méthode de moyenne pondérée (wtaver dans Matlab) sont adoptés. Les entrées
de chaque régulateur flou sont l’erreur (e) et la variation d’erreur (Δe) de la variable
contrôlée (P ou Q). Alors que, la sortie de chaque CLF (P ou Q) est l’incrément de la
consigne de la composante d ou q du courant de référence (Δidref ou Δiqref). Les lois de
commande en sortie des CLFs sont donc calculées comme suit :
80
Chapitre III : Commande des systèmes de production connectés au réseau
Degree of membership
1
Degree of membership
0.8
0.8
0.6
0.6
0.4 0.4
0.2 0.2
0 0
-1 -0.5 0 0.5 1 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1
error and change of error (e and de) control output
Δe
NB NM NS ZE PS PM PB
e
NB NB NB NB NB NM NS ZE
NM NB NM NM NM NS ZE PS
NS NB NM NS NS ZE PS PM
ZE NM NM ZE ZE ZE PM PM
PS NM NS ZE PS PS PM PB
PM NS ZE PS PM PM PM PB
PB ZE PS PM PB PB PB PB
Tableau III.4 : Base des règles du CLF des puissances active et réactive.
Par ailleurs, les erreurs des courants selon les axes d et q sont contrôlés séparément
à l’aide des correcteurs PI. Les gains de ces derniers sont déterminés à l’aide de la
méthode de placement des pôles. Comme l’indique la figure (III.30), les tensions de
référence fournies à l’entrée du modulateur SPWM sont utilisées pour générer les
signaux de commande de l’onduleur 3L-NPC.
III.4.3.Résultats de simulation
Afin de valider l’efficacité du schéma de commande adopté deux tests de
simulation sont effectués :
Un test pour un changement des consignes des puissances.
Un test en présence d’un creux de tension du réseau.
Le comportement du système de GD est examiné pour une variation de la consigne
de la puissance active (Pref) à l’instant t=0.3s de 100 à 150kW, puis un changement à
l’instant t=0.5s vers la valeur de 100kW. La référence de la puissance réactive croît de
zéro à la valeur de 30kVAr durant la période de temps [0.3s-0.5s].
La figure (III.32) illustre la tension et le courant de sortie de la PAC-OS en régimes
transitoire et statique. L’évolution de la tension du bus continu est décrite par la figure
(III.33). Il est clair que la PAC-OS fonctionne au point de puissance correspondant à la
valeur nominale (100kW) durant toute la période de simulation. La figure (III.33)
démontre que le CLF de tension a maintenu la tension du bus CC constamment fixe à
la valeur désirée de 700V. De plus, la figure (III.34) montre que les puissances active et
réactive suivent correctement leurs consignes correspondantes grâce aux CLFs
adoptés. Le contrôle découplé des puissances active et réactive est assuré avec des
81
Chapitre III : Commande des systèmes de production connectés au réseau
300 360
200
340
100
0 320
0 0.05 0.1 0.15 0.2 0.25 0 0.05 0.1 0.15 0.2
t(s) t(s)
(a) (a)
400 400
FC current (A)
FC current (A)
200 300
200
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0 0.05 0.1 0.15 0.2
t(s) t(s)
(b) (b)
Figure III.32 : Evolution de la tension et du courant de sortie de la PAC-OS durant le premier cas
de simulation.
900 900
800 850
800
700
750
DC-link voltage (V)
600
DC-link voltage (V)
700
500
650
400
600
300
550
200
500
100
450
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0 0.05 0.1 0.15 0.2
t(s) t(s)
82
Chapitre III : Commande des systèmes de production connectés au réseau
5
x 10
Pgrid
-1
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
t(s)
4 (a)
x 10
Reactive power (VAr)
4 Qgrid
Qload
Qinverter
2
-2
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
t(s)
(b)
Figure III.34 : Variation des puissances active et réactive de l’onduleur, du réseau et de la
charge CA durant le premier cas de simulation.
5
Pinverter(W) and Qinverter(VAr)
x 10
Pinv
Qinv
1
-1
0 0.1 0.2 0.3 0.4 0.5 0.6
t(s)
(a)
800
DC-link voltage (V)
600
400
200
0
0 0.1 0.2 0.3 0.4 0.5 0.6
t(s)
(b)
Figure III.35 : Evolution de : a) Puissances active et réactive d’onduleur, b) Tension du bus-CC
durant le deuxième cas de simulation.
III.5. Conclusion
Dans ce chapitre, trois types de systèmes de production connectés au réseau
électrique sont étudiés. L’objectif visé était d’étudier le comportement dynamique des
sources d’énergie (PV, éolienne et PAC-OS) et d’examiner les schémas de commande
les plus populaires dans la littérature. Des TIAs sont adoptées dans le but d’optimiser
les méthodes classiques de régulation des correcteurs PI ainsi que pour remplacer les
contrôleurs PI classiques. L’algorithme de PSO présente une alternative très
prometteuse pour la synthèse des régulateurs PI classiques. Son avantage principal
réside dans la simplicité et la non exigence de modélisation mathématique du système
commandé. Par ailleurs, les contrôleurs flous peuvent efficacement substituer les
contrôleurs classiques et assurer des bonnes performances dynamiques et statiques.
83
Chapitre III : Commande des systèmes de production connectés au réseau
84
Chapitre IV
IV.1. Introduction
Ces dernières décennies, les RNAs reçoivent un intérêt majeur de la part de la
communauté de recherche dans différents domaines. Dans le domaine d’électronique
de puissance, les RNAs présentent un outil très efficace pour la commande non linéaire
des systèmes énergétiques qui subissent divers types de perturbation et des variations
paramétriques. Ceci se justifie par les propriétés pertinentes des RNAs à savoir la
capacité d’entraînement, d’approximation et de généralisation.
Comme nous l’avons évoqué précédemment, on discrimine selon la manière avec
laquelle l’apprentissage est réalisé deux catégories [12]: les RNAs entraînés hors-ligne
et les RNAs entraînés en-ligne. L’apprentissage hors-ligne des réseaux neuronaux est
couramment accompli en utilisant une base de données collectées auparavant via des
expériences ou des modèles mathématiques du système à commander. Néanmoins, la
tâche d’apprentissage hors-ligne d’un RNA est souvent lente. Dans certaines
situations, la collecte des patrons d’apprentissage qui couvrent une large gamme de
fonctionnement du système commandé est laborieuse. Lorsque le système est difficile
à modéliser et/ou les données expérimentales sont insuffisantes, l’utilisation des
RNAs qui s’entraînent en-ligne s’avère plus adéquate. Ces derniers possèdent
l’aptitude de s’adapter dynamiquement aux incertitudes des systèmes contrôlés [12].
Dans ce chapitre, nous nous intéressons à la problématique de la commande
intelligente des micro-réseaux hybrides. D’abord, on commence par un état de l’art sur
les études antérieures réalisées sur les MRs hybrides. Ensuite, nous allons présenter
un bref aperçu sur les techniques de commande neuronale déjà proposées dans la
littérature pour la MPPT des sources PV et éolienne. Après, nous allons étudier le
comportement d’un microréseau hybride (CA/CC) en adoptant une stratégie de
commande intelligente basée sur les RNAs adaptatifs. Enfin, les résultats de
simulation sont fournis afin d’affirmer l’efficacité et la robustesse de la structure de
supervision/commande proposée face au changement des conditions de
fonctionnement.
85
Chapitre IV : Commande neuronale d’un micro-réseau (CA/CC) hybride
86
Chapitre IV : Commande neuronale d’un micro-réseau (CA/CC) hybride
Dans [46], l’étude est faîte sur un MR connecté au réseau qui se compose d’une
source PV, d’une PAC-OS et d’un banc de batteries. Les auteurs ont proposé des
contrôleurs flous par mode glissant pour assurer le contrôle découplé des puissances
active et réactive. Des résultats de simulation dans Matlab/Simulink sont présentés
afin de valider le schéma de contrôle suggéré.
Un schéma de contrôle basé sur un superviseur ANFIS central et des contrôleurs
flous par mode glissant est proposé dans [60] pour un SEH (PAC/super-condensateur)
connecté au réseau. Leur objectif de commande était de garantir un fonctionnement
correct du MR même lors l’occurrence du creux de tension du réseau. L’efficacité de la
méthode proposée est approuvée par simulation dans Matlab/Simulink.
Une structure de commande en puissance basée sur des ANFISs est proposée dans
[150] pour le pilotage d’un onduleur au sein d’un SEH (PV/éolien/PAC/
électrolyseur/réservoir d’hydrogène/SSEB) couplé au réseau. En tenant compte le
SOC du SSEB, un superviseur d’énergie centralisé à base d’ANFIS est également
développé pour déterminer les puissances qui doivent être générées ou stockées dans
le système d’hydrogène et le SSEB. Les résultats de simulation sont fournis afin de
démontrer la faisabilité de ce système de commande.
Un SGE centralisé établi à base de la LF est proposé dans [151] afin d’assurer la
supervision des flux énergétiques dans un SEH (PV/TE/PAC/SSEB/électrolyseur)
autonome. Ce superviseur assigne les consignes des puissances à l’électrolyseur et à la
PAC suivant la puissance demandée par la charge et le SOC du SSEB. Les
performances du SGE flou sont évaluées par simulation.
La tendance actuelle de recherche est l’utilisation des TIAs pour l’ajustement des
paramètres ou la substitution des contrôleurs classiques insérés dans les structures de
commande traditionnelles. Néanmoins, la stratégie de supervision/commande d’un
MR doit être choisie de sorte qu’un bon compromis entre l’efficacité et la simplicité
soit assuré.
IV.3. Aperçu sur les techniques de MPPT des SERs utilisant des RNAs
Comme déjà mentionné, la production d’énergie PV et éolienne est très fluctuante
en raison de leur forte dépendance des conditions météorologiques. En conséquence,
il est indispensable d’associer aux SERs des organes de commande qui se chargent de
la recherche des points de fonctionnement optimal afin d’extraire en permanence le
maximum d’énergie disponible. D’ailleurs, un algorithme de MPPT doit être capable
de s’adapter dynamiquement aux changements climatiques en assurant des bonnes
performances de poursuite du PPM.
De nombreux travaux antérieurs ont focalisé leurs études sur la thématique de
MPPT des SERs dont l’objectif de commande était d’améliorer le rendement
énergétique des SERs. Récemment, un grand intérêt est accordé aux contrôleurs MPPT
basés sur des RNAs, notamment, ceux qui qui s’entraînent en ligne. Grâce à la capacité
d’entraînement et de traitement des non-linéarités, les suiveurs neuronaux peuvent
assurer des bonnes performances statiques et dynamiques qui surpassent celles des
algorithmes MPPT classiques.
87
Chapitre IV : Commande neuronale d’un micro-réseau (CA/CC) hybride
Le Tableau (IV.1) (voir page 89) donne un bref aperçu sur les RNAs déjà proposés
dans la littérature pour la commande MPPT des sources PV et éolienne. On mentionne
pour chaque RNA : le type, la structure, les entrées/sorties, l’algorithme
d’apprentissage ainsi que la fonction coût associée.
Les remarques qui peuvent être faites à propos ces travaux sont les suivantes :
L’algorithme de RP basé sur la règle de la descente du gradient est le plus utilisé
pour l’entraînement des RNAs.
Tous les RNAs évoqués utilisent le type d’apprentissage supervisé.
Le type des fonctions d’activation, le taux d’apprentissage, les termes de biais
sont les facteurs qui influent d’une manière significative sur les performances
des RNAs.
A l’opposé du RNA entrainé hors ligne, le RNA qui apprend en ligne ne requiert
pas une connaissance a priori du modèle et des caractéristiques de la SER (PV ou
éolienne). En outre, il n’exige que peu d’information provenant du système afin de
décider l’action de commande adéquate [12,157]. Vue ces avantages, nous nous
focalisons dans cette étude sur les RNAs entrainés en ligne comme des moyens de
contrôle intelligent et adaptatif.
Dans ce chapitre, une étude sur un MR hybride (CA/CC) fonctionnant en mode
couplé au réseau électrique est achevée. Ce système consiste en une agrégation d’une
source PV, d’un GE, d’une PAC-OS, d’un SSEB et des charges locaux à CC et à CA. Un
SGE centralisé basé sur la LF est proposé pour la supervision des flux énergétiques
dans le MR. Des contrôleurs neuronaux avec une adaptation en-ligne de leurs
paramètres sont proposés pour la poursuite des PPMs des SERs. De plus, des
contrôleurs linéaires adaptatifs basés sur des ADALINEs sont introduits pour la
commande du SSEB et de la PAC-OS. En outre, un schéma de commande prédictive
est appliqué au SGEE afin d’achever un contrôle découplé du couple
électromagnétique et du flux de la GSAP. D’autre part, une stratégie de contrôle direct
de puissance avec orientation du flux virtuel (VF-DPC : Virtuel Flux Direct Power
Control) est adoptée pour le pilotage de l’onduleur triphasé. Un FFNN est proposé
pour la régulation de la tension du bus continu. Alors que, deux contrôleurs basés sur
des ENNs sont adoptés pour assurer le contrôle du transfert bidirectionnel de la
puissance active ainsi que la compensation de la puissance réactive de charge CA.
L’algorithme d’interaction adaptative est utilisé pour l’ajustement en temps réel des
poids des connexions des réseaux FFNN et ENN. Le MR étudié est simulé dans
l’environnement Matlab/Simulink pour différents scénarios afin de vérifier et valider
l’efficacité des contrôleurs neuronaux proposés [158].
88
Type de Référence et Algorithme
Type de RNA Entrées Sortie Fonction coût
source année d’apprentissage
FFNN à trois - Vitesse du générateur. Erreur quadratique
Incrément de l’angle de - En ligne.
[152], 2012. couches - Consigne de vitesse. de vitesse de
phase du GE. - Règle de Widrow–Hoff.
(3×4×1) - Sortie du RNA. rotation.
- Supervisé, en ligne.
- Erreur de vitesse de
- RP basé sur la descente du Erreur quadratique
ENN rotation de la TE. Consigne du courant en
[153], 2012. gradient. de la puissance de la
- Variation d’erreur de quadrature de la GSAP.
(2×9×1) - Ajustement des taux TE.
vitesse.
d’apprentissage via PSO.
- Supervisé, en ligne.
- Erreur de vitesse de
Éolienne - RP basé sur la descente du Erreur quadratique
RBFN rotation. Consigne du courant en
[154], 2013. gradient. de la vitesse de
- Courant en quadrature de quadrature de la GSAP.
(2×9×1) - Adaptation des paramètres rotation.
la GSAP.
d’apprentissage par PSO.
- Erreur de puissance de la
- Supervisé, en ligne. Erreur quadratique
génératrice. Référence de l’angle de
[54], 2013. ENN modifié - RP basé sur la descente du de la puissance de la
- Dérivée d’erreur de calage.
gradient. génératrice.
puissance.
- Supervisé, en ligne. Produit de la surface
- Surface de glissement Consigne de l’angle de
[155],2014. RBFN - RP basé sur la descente du de glissement par sa
(S(t)). calage.
gradient. dérivée (𝑆. 𝑆̇).
RBFN - Tension PV. - Supervisé.
Erreur quadratique
[54], 2013. - Courant PV. Tension de référence. - RP basé sur la descente du
(3×9×1) de tension PV.
- Température. gradient.
- Tension PV.
GRNN - Algorithme modifié de Définie en fonction
[155],2014. - Courant PV. Tension de référence.
PSO. de l’erreur de tension.
Quatre couches. - Température.
PV - Supervisé, hors ligne.
- Tension PV.
[156],2015. FFNN Tension de référence. - RP basé sur la méthode de Ce n’est pas précisé.
- Courant PV.
Levenberg-Marquardt.
- Erreur de la pente de la Erreur quadratique
Un seul neurone - Supervisé, en ligne.
courbe PPV(VPV) (𝑒 = 0 − Rapport cyclique du de la pente de la
[157], 2015. caché à fonction - RP basé sur la descente du
𝑑𝑃𝑃𝑉 ⁄𝑑𝑉𝑃𝑉 ). hacheur buck. courbe PPV(VPV) (𝑒 2 /
sigmoïdale. gradient.
- Dérivée d’erreur. 2).
Tableau IV.1 : Aperçu sur les RNAs utilisés pour la MPPT des sources PV et éolienne.
89
Chapitre IV : Commande neuronale d’un micro-réseau (CA/CC) hybride
Charge
CC
TE+ Convertisseur
GSAP CA-CC
Charge
isabc
Contrôleur du CA
Vw , Te , ωe GE
+ +
Convertisseur Convertisseur Filtre
PAC-OS Vfc VDC CDC
CC-CC CC-CA LC
- -
Ifc ica,b,c Réseau
Contrôleur de
VDC
Pfc*, Vfc PAC VFDPC-SVM Qref
SOC Pnet
+
Convertisseur Convertisseur
PV VPV SSEB SGE-LF
CC-CC CC-CC
- Pb* Pfc*
IPV Ib
Contrôleur du
Pb*
VPV Contrôleur du
GPV SSEB Vb
90
Chapitre IV : Commande neuronale d’un micro-réseau (CA/CC) hybride
basée sur des RNAs adaptatifs est développée afin d’assurer une gestion efficace
d’énergie dans le MR hybride. Dans ce qui suit, les organes de commande du système
de supervision/commande proposé seront amplement présentés.
IV.5.1. Commande neuronale de la source PV
Comme déjà mentionné, l’algorithme d’IncCond est largement utilisé pour la
poursuite des PPMs des sources PV. Cette technique de MPPT assure la régulation de
la tension de sortie du GPV selon la valeur de la consigne correspondante au point de
fonctionnement optimal (voir section (III.2.2.1)). Les avantages principaux de
l’algorithme d’IncCond résident dans sa simplicité et sa facilité de mise en œuvre.
Toutefois, sa vitesse de convergence et leurs oscillations statiques autour du PPM
dépendent notamment de la largeur du pas de perturbation de la tension de référence :
Un incrément (ΔVref) assez grand de la consigne (Vref) augmente simultanément
l’erreur statique et la vitesse de convergence du contrôleur d’IncCond et vice versa.
Dans cette étude, un neuro-contrôleur basé sur un RBFN à un seul neurone caché
(SN-RBFN: Single Neuron Radial Basis Function Network) est adopté pour accomplir
la tâche de MPPT du GPV. Le but visé ici est d’améliorer la précision de poursuite du
PPM et d’optimiser les performances dynamiques de l’algorithme d’IncCond. La
régulation du point d’opération de la source PV est donc réalisée à l’aide du
convertisseur survolteur. La figure (IV.2) illustre le schéma de commande du hacheur
PV établi dans l’environnement Matlab/Simulink.
Le suiveur du PPM proposé est conçu en se basant sur le principe de la méthode
d’IncCond [113,115]. Où, la capacité d’apprentissage du SN-RBFN garantit son auto-
adaptation aux changements des conditions climatiques. Donc, le contrôleur SN-RBFN
adopté est composé d’un seul neurone dans la couche cachée dont la fonction
d’activation gaussienne est définie par [101,159]:
‖𝒙−𝒄‖2
𝑓 (𝒙) = 𝑒𝑥𝑝 (− 2𝑏 2
)=ℎ (IV.1)
Avec c=[c11,c21,c31] est le vecteur centre de la fonction gaussienne f(x), b est la largeur
de f(x) et x=[x1,x2,x3] est le vecteur des entrées du SN-RBFN. ‖ ‖ désigne la norme
Euclidienne. La sortie du réseau SN-RBFN est calculée comme suit [159]:
𝑦 = 𝑉𝑝𝑣𝑟𝑒𝑓 = 𝑎0 + 𝑎1 𝑓(𝒙) (IV.2)
a0 et a1 représentent, respectivement, le biais et le poids du SN-RBFN et Vpvref est la
tension de référence du GPV. Les entrées choisies du régulateur SN-RBFN sont :
La conductance instantanée (x1=IPV/VPV).
La conductance incrémentale (x2=ΔIPV/ΔVPV).
L’erreur de la tension de référence (x3=ΔVpvref(k)=Vpvref(k)-Vpvref(k-1)).
Comme l’indique la figure (IV.2), un contrôleur PI se charge de la régulation de la
tension du GPV (VPV) pour qu’elle suive sa consigne issue du SN-RBFN. Le rapport
cyclique en sortie de ce correcteur est ensuite comparé avec un signal triangulaire afin
de générer le signal de commande du hacheur.
91
Chapitre IV : Commande neuronale d’un micro-réseau (CA/CC) hybride
92
Chapitre IV : Commande neuronale d’un micro-réseau (CA/CC) hybride
Une fois Gin(k) s’annule, sous certaines conditions climatiques, le contrôleur SN-
RBFN stabilise au point de fonctionnement atteint qui correspond au PPM.
IV.5.2. Commande prédictive du SGEE
Actuellement, la commande prédictive apparaît comme une alternative très efficace
et adéquate pour le pilotage et le contrôle des convertisseurs de puissance et des
machines électriques. De manière générale, le principe d’une méthode de contrôle
prédictif consiste à utiliser le modèle du système à commander pour la prédiction du
comportement futur des variables contrôlées. La décision sur l’action de commande
adéquate est faîte par l’intermédiaire d’un critère de performance [72]. Le schéma de
commande prédictive à réponse pile (en anglais: deadbeat control) est le plus répandu.
La structure de commande du redresseur triphasé, reliant la GSAP au bus CC, est
illustrée par la figure (IV.3). Un organe de commande MPPT basé sur un réseau
ADALINE est implémenté afin d’assurer un asservissement précis de la vitesse de
rotation du GE. De plus, une méthode de contrôle prédictif à réponse pile (DPCM :
Deadbeat Predictive Control Method) [161-163] est appliquée pour le pilotage du
convertisseur côté génératrice. Le rôle principal du contrôleur local du SGEE est de
poursuivre en permanence le PPM du GE ainsi que de contrôler séparément le couple
électromagnétique et le flux de la GSAP.
93
Chapitre IV : Commande neuronale d’un micro-réseau (CA/CC) hybride
94
Chapitre IV : Commande neuronale d’un micro-réseau (CA/CC) hybride
Te*(k+1)
Te*(k)
Te*(k-1)
Ts Ts
Figure IV.4 : Prédiction de la consigne future du couple Te* [161].
En remplaçant par (IV.15) et (IV.16) dans (IV.14), les composantes selon les axes d
et q du vecteur de tension statorique de référence sont finalement calculées par :
2 𝐿
𝑣𝑠𝑑 (𝑘) = −𝑅𝑠 𝑖𝑠𝑑 (𝑘) + 3𝑝𝜙 𝐿𝑠 𝜔𝑒 𝑇𝑒 (𝑘) − 𝑇𝑠 ∆𝑖𝑠𝑑 (𝑘)
𝑟 𝑠
{ 2𝑅𝑠 2 𝐿𝑠 (IV.17)
𝑣𝑠𝑞 (𝑘) = − 3𝑝𝜙 𝑇𝑒 (𝑘) − 𝐿𝑠 𝜔𝑒 𝑖𝑠𝑑 (𝑘) + 𝜔𝑒 𝜙𝑟 − 3𝑝𝜙 (∆𝑇𝑒 (𝑘) + 𝑑𝑇𝑒∗ (𝑘))
𝑟 𝑟 𝑇𝑠
Avec ∆𝑇𝑒 (𝑘) et ∆𝑖𝑠𝑑 (𝑘) sont, respectivement, les erreurs instantanées de poursuite
des consignes du couple et du courant direct. 𝑑𝑇𝑒∗ (𝑘) représente la variation actuelle de
95
Chapitre IV : Commande neuronale d’un micro-réseau (CA/CC) hybride
la consigne du couple :
∗ (𝑘)
∆𝑖𝑠𝑑 (𝑘) = 𝑖𝑠𝑑 − 𝑖𝑠𝑑 (𝑘)
{ ∆𝑇𝑒 (𝑘) = 𝑇𝑒∗ (𝑘) − 𝑇𝑒 (𝑘) (IV.18)
𝑑𝑇𝑒∗ (𝑘) = 𝑇𝑒∗ (𝑘) − 𝑇𝑒∗ (𝑘 − 1)
96
Chapitre IV : Commande neuronale d’un micro-réseau (CA/CC) hybride
source secondaire qui permet d’améliorer la réponse du MR aux besoins des charges
électriques. Comme le montre la figure (IV.6), un contrôleur basé sur un réseau
ADALINE est adopté pour le contrôle du courant du banc de batteries (Ib). Ce
régulateur doit garantir la minimisation de l’écart entre la consigne et la mesure du
courant du SSEB pour différents modes d’opération.
Les deux entrées du contrôleur du courant du SSEB sont choisies comme : l’erreur
actuelle du courant 𝑒𝑏 (𝑘) et l’erreur précédente 𝑒𝑏 (𝑘 − 1). Tel que, l’erreur du courant
(𝑒𝑏 (𝑘)) est égale à la différence entre le courant de référence (Ib*(k)) et le courant
mesuré en sortie du SSEB (Ib(k)) (𝑒𝑏 (𝑘) =Ib*(k)-Ib(k)). La référence du courant du SSEB
(Ib*) est déterminée en divisant la consigne de puissance (Pb*) issue du SGE par la
tension de sortie du SSEB. Le rapport cyclique Db(k) en sortie d’ADALINE est comparé
avec un signal triangulaire afin de générer les signaux de commutation du
convertisseur bidirectionnel.
L’adaptation des poids du réseau ADALINE est réalisée en utilisant la surface
d’erreur suivante :
𝐸𝑏 (𝑘) = 0 − 𝑆𝑏 (𝑘) = −[𝜆2 . 𝑒𝑏 (𝑘) + 𝑑𝑒𝑏 (𝑘)] (IV.21)
Tel que 𝑆𝑏 (𝑘) est la surface de glissement du courant de SSEB. 𝜆2 est une constante
positive et 𝑑𝑒𝑏 (𝑘) = 𝑒𝑏 (𝑘) − 𝑒𝑏 (𝑘 − 1) . Le vecteur des poids d’ADALINE (𝑊𝑏 ) est
modifié en-ligne selon la règle d’entraînement de LMS comme suit :
𝑊𝑏 (𝑘 + 1) = 𝑊𝑏 (𝑘) + 2𝛼3 𝐸𝑏 (𝑘)𝑋𝑏 (𝑘) (IV.22)
𝑋𝑏 et 𝛼3 présentent, respectivement, le vecteur d’entrée et le taux d’apprentissage du
contrôleur de SSEB. La vérification de la condition 𝑆𝑏 (𝑘) = 0, en régime permanent,
assure l’élimination de l’erreur statique du courant du SSEB.
97
Chapitre IV : Commande neuronale d’un micro-réseau (CA/CC) hybride
q β
Ψ𝐿
Ψ𝑖𝑛𝑣
𝑣𝑔 d
𝑖𝑐
Ψ𝑔
𝑖𝑔𝑞 𝑖𝑔𝑑 α
𝜃Ψ
98
Chapitre IV : Commande neuronale d’un micro-réseau (CA/CC) hybride
Ψ𝑔𝛼 = ∫ 𝑣𝛼 𝑑𝑡 − 𝐿 𝑖𝑐𝛼
{ (IV.25)
Ψ𝑔𝛽 = ∫ 𝑣𝛽 𝑑𝑡 − 𝐿 𝑖𝑐𝛽
(𝑣𝛼 , 𝑣𝛽 ) et (𝑖𝑐𝛼 , 𝑖𝑐𝛽 ) représentent, respectivement, les composantes αβ des vecteurs
de tension et du courant de sortie d’onduleur. Tenant compte l’équation (II.60), les
composantes αβ du FV se calculent en fonction des états de commutation de l’onduleur
(Sa,Sb,Sc), des composantes αβ du courant et de la tension mesurée du bus continu (VDC)
comme suit :
𝑉𝐷𝐶
Ψ𝑔𝛼 = ∫ (2𝑆𝑎 − 𝑆𝑏 − 𝑆𝑐 )𝑑𝑡 − 𝐿 𝑖𝑐𝛼
3
{ 𝑉𝐷𝐶 (IV.26)
Ψ𝑔𝛽 = ∫ (𝑆𝑏 − 𝑆𝑐 )𝑑𝑡 − 𝐿 𝑖𝑐𝛽
√3
Dans cette étude, nous proposons de remplacer les contrôleurs PI compris dans le
schéma de VFDPC-SVM par des contrôleurs neuronaux adaptatifs afin d’améliorer les
performances dynamiques et statiques. Le but est de rendre l’onduleur capable de
réagir dynamiquement aux changements des conditions opérationnelles du MR.
Comme le montre la Figure (IV.8), un réseau FFNN est utilisé dans la boucle externe de
contrôle de la tension du bus CC. En outre, deux réseaux récurrents d’Elman (ENNs)
99
Chapitre IV : Commande neuronale d’un micro-réseau (CA/CC) hybride
100
Chapitre IV : Commande neuronale d’un micro-réseau (CA/CC) hybride
(𝑙) (𝑙) 2
(𝑙) (𝑙) (𝑙) 1−𝑒 −𝑛𝑒𝑡𝑛 (𝑙) ′ (𝑙) 𝜕𝑥𝑛 1 (𝑙)
𝑥𝑛 = 𝑓𝑛 (𝑛𝑒𝑡𝑛 ) = (𝑙) ⟹ 𝑓𝑛 (𝑛𝑒𝑡𝑛 ) = (𝑙) = 2 (1 − (𝑥𝑛 ) ) (IV.32)
1+𝑒 −𝑛𝑒𝑡𝑛 𝜕𝑛𝑒𝑡𝑛
Comme déjà vu dans le chapitre II, le réseau d’Elman est constitué de quatre couches
de neurones. On décrit dans ce qui suit la fonction principale de chaque couche du
contrôleur ENN :
Les neurones de la CE utilisent des fonctions de transfert de type identité. La
sortie de chaque neurone d’entrée est donnée par :
(1) (1) (1) (1)
𝑥𝑖 (𝑘) = 𝑓𝑖 (𝑛𝑒𝑡𝑖 ) = 𝑛𝑒𝑡𝑖 , 𝑖 = 1,2 (IV.33)
(1)
𝑛𝑒𝑡𝑖 est la iéme entrée du ENN.
La fonction de chaque neurone de la CH peut être décrite par l’équation
suivante :
(2) (2) (2) (2) (3) (1)
𝑥𝑗 (𝑘) = 𝑓𝑗 (𝑛𝑒𝑡𝑗 ) = 𝑓𝑗 (∑𝑟 𝑤𝑟𝑗 𝑥𝑟 (𝑘) + ∑𝑖 𝑤𝑖𝑗 𝑥𝑖 (𝑘)) , 𝑗 = 1, . . ,5 (IV.34)
(2)
Avec 𝑥𝑗 est la sortie du neurone j de la CH, 𝑤𝑖𝑗 est le poids de connexion entre le
(3)
neurone i de la CE et le neurone j de la CH, 𝑥𝑟 (𝑘) est la sortie du riéme neurone de
la couche contexte, 𝑤𝑟𝑗 est le poids reliant le neurone contexte r au neurone caché j,
(2)
et 𝑓𝑗 est la fonction d’activation du jiéme neurone de la CH.
La sortie du neurone r de la Couche ConteXte (CX) est définie comme :
(3) (2)
𝑥𝑟 (𝑘) = 𝑥𝑗 (𝑘 − 1), 𝑟 = 1, . . ,5 (IV.35)
La sortie du neurone unique de la CS est calculée comme suit :
(4) (4) (4) (4) (2)
𝑥𝑜 (𝑘) = 𝑓𝑜 (𝑛𝑒𝑡𝑜 ) = 𝑓𝑜 (∑𝑗 𝑤𝑜𝑗 𝑛𝑒𝑡𝑗 (𝑘)) (IV.36)
(4) (4)
𝑥𝑜 est le signal de sortie du ENN, 𝑓𝑜 est la fonction tangente-hyperbolique et
𝑤𝑗𝑜 est le poids reliant le neurone j de la CH au neurone de sortie.
En tenant compte l’équation (IV.32), les poids des connexions (𝑤𝑖𝑗 ,𝑤𝑟𝑗 , 𝑤𝑜𝑗 ) du
réseau ENN sont ajustés en-ligne selon la règle d’AI (équation (IV.31)) comme suit:
(1)
1 (2) 2 𝑥𝑖
∆𝑤𝑖𝑗 (𝑘) = 2 [1 − (𝑥𝑗 ) ] (2) 𝑤𝑜𝑗 ∆𝑤𝑜𝑗 𝑃𝑜𝑖𝑑𝑠 𝑒𝑛𝑡𝑟𝑒 𝑙𝑎 𝐶𝐸 𝑒𝑡 𝑙𝑎 𝐶𝐻
𝑥𝑗
(3)
1 (2) 2 𝑥𝑟 (IV.37)
∆𝑤𝑟𝑗 (𝑘) = 2 [1 − (𝑥𝑗 ) ] (2) 𝑤𝑜𝑗 ∆𝑤𝑜𝑗 𝑃𝑜𝑖𝑑𝑠 𝑒𝑛𝑡𝑟𝑒 𝑙𝑎 𝐶𝑋 𝑒𝑡 𝑙𝑎 𝐶𝐻
𝑥𝑗
𝛾 (4) 2 (2)
{ ∆𝑤𝑜𝑗 (𝑘) = − 2 [1 − (𝑥𝑜 ) ] 𝑥𝑗 𝑒𝑛 𝑃𝑜𝑖𝑑𝑠 𝑒𝑛𝑡𝑟𝑒 𝑙𝑎 𝐶𝐻 𝑒𝑡 𝑙𝑎 𝐶𝑆
deP,Q (k)
.
Vcd,q
Z-1 +
Z-1 .
Couche de sortie
Z-1
Couche contexte 𝑤𝑟𝑗 5 neurones 𝑊𝑜𝑗
Comme déjà suggéré dans [155], la condition d’invariance 𝑆𝑃,𝑄 (𝑘). 𝑑𝑆𝑃,𝑄 (𝑘) = 0 qui
doit être vérifiée en mode de glissement est adoptée comme objectif d’apprentissage
101
Chapitre IV : Commande neuronale d’un micro-réseau (CA/CC) hybride
d’ENN. Donc, le terme 𝑆𝑃,𝑄 (𝑘). 𝑑𝑆𝑃,𝑄 (𝑘) est utilisé au lieu de l’erreur de sortie (𝑒𝑛 ) dans
les règles d’adaptation de l’équation (IV.37), tel que :
𝑆𝑃,𝑄 (𝑘) = 𝜆3 . 𝑒𝑃,𝑄 (𝑘) + 𝑑𝑒𝑃,𝑄 (𝑘)
{ (IV.38)
𝑑𝑆𝑃,𝑄 (𝑘) = 𝑆𝑃,𝑄 (𝑘) − 𝑆𝑃,𝑄 (𝑘 − 1)
Avec 𝜆3 est une constante positive. 𝑆𝑃,𝑄 et 𝑑𝑆𝑃,𝑄 représentent la surface de glissement et
sa variation pour la puissance active et réactive (P et Q). Donc, l’objectif de commande
des ENNs proposés est de diriger dans un temps fini puis maintenir les variables d’état
P et Q, respectivement, sur leurs surfaces de commutation correspondantes (𝑆𝑃 et 𝑆𝑄 ).
Avec 𝑁𝐸𝑇𝑖(1) est la ième entrée du FFNN, 𝑋𝑖(1) (𝑘) est la sortie i de la CE, 𝑋𝑗(2) (𝑘) est la
sortie du neurone j de la CH, 𝑁𝐸𝑇𝑗(2)et 𝑁𝐸𝑇𝑜(3)sont les niveaux d’activation du neurone
caché j et du neurone de sortie. 𝑋𝑜(3) (𝑘) est la sortie du FFNN, 𝜔𝑖𝑗 désigne le poids de
connexion entre le neurone i de la CE et le neurone j de la CH, 𝜔𝑗𝑜 est le poids reliant
(2) (3)
le jième neurone de la CH au neurone de sortie, 𝐹𝑗 , 𝐹𝑜 sont, respectivement, les
fonctions tangente-hyperbolique de la CH et de la CS. L’erreur actuelle de la tension du
bus CC (𝑒𝑉 (𝑘)= VDC*(k)-VDC(k)) ainsi que sa valeur précédente (eV(k-1)) représentent les
entrées du contrôleur FFNN proposé. Le signal de commande en sortie du FFNN est
multiplié par la tension mesurée (VDC) afin de générer la consigne actuelle de la
puissance active.
L’entraînement du réseau FFNN est réalisé à l’aide de la règle d’AI selon les
expressions suivantes :
(2) 2
1 [1−(𝑋𝑗 ) ] (1)
∆𝜔𝑖𝑗 (𝑘) = (2) 𝑋𝑖 𝜔𝑗𝑜 ∆𝜔𝑗𝑜 𝑃𝑜𝑖𝑑𝑠 𝑒𝑛𝑡𝑟𝑒 𝑙𝑎 𝐶𝐸 𝑒𝑡 𝑙𝑎 𝐶𝐻
2 𝑋𝑗 (IV.40)
𝛾1 (3) 2 (2)
∆𝜔𝑗𝑜 (𝑘) = − [1 − (𝑋𝑜 ) ] 𝑋𝑗 𝑒𝑛 𝑃𝑜𝑖𝑑𝑠 𝑒𝑛𝑡𝑟𝑒 𝑙𝑎 𝐶𝐻 𝑒𝑡 𝑙𝑎 𝐶𝑆
{ 2
𝛾1 est le taux d’apprentissage. L’erreur de sortie (𝑒𝑛 ) dans (IV.40) est remplacée par le
terme 𝑆𝑉 (𝑘) + 𝑑𝑆𝑉 (𝑘), tel que :
𝑆 (𝑘) = 𝜆4 . 𝑒𝑉 (𝑘) + 𝑑𝑒𝑉 (𝑘)
{ 𝑉 (IV.41)
𝑑𝑆𝑉 (𝑘) = 𝑆𝑉 (𝑘) − 𝑆𝑉 (𝑘 − 1)
102
Chapitre IV : Commande neuronale d’un micro-réseau (CA/CC) hybride
X(1) 1
dW11 K Ts W11 net(2) 1
1 1
Xi(1) X(1) 2 z-1
Discrete-T ime x(2) 1
Integrator2 tansig
2 1/u
2
X(2) 1
4 dW21 K T s W21 net(2) 1
dW01 3 z-1
W01 Discrete-T ime
Integrator1
1/2
1
Figure IV.10 : Exemple d’ajustement en-ligne des poids du FFNN par AI dans Matlab/Simulink .
103
Chapitre IV : Commande neuronale d’un micro-réseau (CA/CC) hybride
0.8
0.6
0.4
0.2
0
-3 -2 -1 0 1 2 3 4
Pnet(W) x 10
4
(a)
L M H
1
Degree of membership
0.8
0.6
0.4
0.2
(b)
NB NS ZE PS PB
1
Degree of membership
0.8
0.6
0.4
0.2
(c)
104
Chapitre IV : Commande neuronale d’un micro-réseau (CA/CC) hybride
PS PM PB
1
Degree of membership
0.8
0.6
0.4
0.2
(d)
Figure IV.11 : FAs des variables a) Pnet, b) SOC, c) Pb*, d) Pfc*.
Pnet
NM NS PS PM PB PB+
SOC
H - - - PM PB PB
M - - PS PB PB PB
L - - PM PB PB PB
Tableau IV.2 : Tableau des règles floues pour la gestion de la PAC-OS.
Pnet
NM NS PS PM PB PB+
SOC
H ZE ZE PB PS PB PB
M NB NS ZE ZE PS PB
L NB NS NB NS NS NS
Tableau IV.3 : Tableau des règles floues pour la gestion du SSEB.
105
Chapitre IV : Commande neuronale d’un micro-réseau (CA/CC) hybride
irradiance (W/m2)
1000
800
600
400
40
20
0 0.1 0.2 0.3 t(s) 0.4 0.5 0.6 0.7
(b)
400
PV voltage (V)
380
360
340
320
0 0.1 0.2 0.3 t(s) 0.4 0.5 0.6 0.7
(c)
400
voltage Vref (V)
300
200
0 0.1 0.2 0.3 t(s) 0.4 0.5 0.6 0.7
4
x 10 (d)
4 2.1129
x 10
2.1128
2.1128
2
2.1128
1.5
PV output power (W)
0.5
IncCond method
proposed SN-RBFN
0
0 0.1 0.2 0.3 t(s) 0.4 0.5 0.6 0.7
(e)
Figure IV.12 : Variation de : a) Irradiation solaire, b) Courant PV, c) Tension PV, d) Tension de
référence en sortie du SN-RBFN, e) Puissance PV de sortie.
Il est clair de ces résultats que le neuro-contrôleur MPPT du GPV s’adapte d’une
manière automatique aux changements d’ensoleillement. A chaque niveau
d’irradiation (G), le suiveur SN-RBFN atteint rapidement le PPM correcte grâce au
processus d’entraînement en temps réel. Au niveau G=1000W/m2, le temps de réponse
du contrôleur SN-RBFN est d’environ 38ms qui est plus court que le temps transitoire
du régulateur d’IncCond (de l’ordre de 92ms). En régime permanent, le SN-RBFN
montre des oscillations autour du PPM considérablement atténuées en comparaison
avec le contrôleur d’IncCond (voir figure (IV.12.e)). L’erreur statique de puissance
obtenue avec le réseau SN-RBFN est de l’ordre de 0.305W qui est plus réduite que
l’erreur commise par le contrôleur d’IncCond (de 0.887W environ). En outre, les
106
Chapitre IV : Commande neuronale d’un micro-réseau (CA/CC) hybride
40
39.5
39
70 0.25 0.3 0.35 t(s) 0.4 0.45 0.5
(a)
60
Rotating speed (rad/s)
50
40
30
20 wr with PI controller
wr with ADALINE controller
10
wref
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7
t(s)
(a)
50
0
Te(N.m) and Isd(A)
-50
Te
Teref
-100
Isd
-150
-200
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7
t(s)
(b)
0.11
stator flux amplitude (wb)
0.1
0.09
0.08
0.07
0.06
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7
t(s)
(c)
880
860
DC-link voltage (V)
840
820
800
780
107
Chapitre IV : Commande neuronale d’un micro-réseau (CA/CC) hybride
108
Chapitre IV : Commande neuronale d’un micro-réseau (CA/CC) hybride
4
x 10
15000 Pfc*
10000 Pbess*
5000 Pbess
0 Pfc
-5000
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
t(s)
Qgrid, Qinv and QLac (Var)
4
x 10 (c)
4
Qinverter
Qgrid
2
QLac
-2
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
t(s)
(d)
Figure IV.14 : Variation de : a) Demande des charges CA et CC, b) Puissances actives de l’onduleur
et du réseau, c) Puissances de sortie de PAC-OS et du SSEB, d) Puissances réactives de l’onduleur,
de charge CA et du réseau principal.
10000
Pfc* (W)
5000 PMS-FL
PMS-states
0
0 0.2 0.4 t(s) 0.6 0.8 1
(a)
10000
Pbess*(W)
5000
0
PMS-FL
-5000 PMS-states
Nous avons accompli un autre test de simulation pour un niveau bas d’irradiation
solaire (G= 0.1kW/m2) et une faible vitesse de vent (VW=8m/s). Le niveau de SOC initial
est ainsi considéré bas (SOCinitial=30%). Les exigences des charges CC et CA sont fixées
comme suit : 𝑃𝐿𝑑𝑐 = 20kW et (PLac=5kW, QLac=0VAr). C’est le cas où la demande de la
charge CC excède la puissance totale disponible du MR. La figure (IV.16) illustre les
résultats obtenus.
La figure (IV.16.a) montre que le régulateur FFNN a réussi à la stabilisation de la
tension du bus CC autour de la valeur de consigne. Nous pouvons constater ainsi que
109
Chapitre IV : Commande neuronale d’un micro-réseau (CA/CC) hybride
800
750
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
t(s)
4
x 10 (a)
Pmg
Pmg, Pgrid and PL (W)
4
PL
2 Pg
-2
-4
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
t(s)
(b)
110
Chapitre IV : Commande neuronale d’un micro-réseau (CA/CC) hybride
4
x 10
0
0 0.1 0.2 0.3 0.4 0.5 0.6
t(s)
(a)
10000
SOFC and BESS powers (W)
Pbess
Pfc
5000
-5000
0 0.1 0.2 0.3 0.4 0.5 0.6
t(s)
(b)
900
DC-link voltage (V)
800
700
600
0 0.1 0.2 0.3 0.4 0.5 0.6
t(s)
(c)
100
inverter current (A)
50
-50
-100
0.25 0.3 0.35 0.4 0.45 0.5
t(s)
(d)
15
alpha-negative component (A) beta-negative component (A)
10
-5
-10
-15
0 0.1 0.2 0.3 0.4 0.5 0.6
t(s)
15
10
5
0
-5
-10
-15
0 0.1 0.2 0.3 0.4 0.5 0.6
t(s)
(e)
111
Chapitre IV : Commande neuronale d’un micro-réseau (CA/CC) hybride
=90%). La figure (IV.18.a) montre que le contrôleur SN-RBFN a assuré des bonnes
performances de poursuite du PPM avec des oscillations statiques de tension plus
réduites que celles de l’algorithme d’IncCond. Le temps de convergence de ce dernier
(de l’ordre de 70ms) est lent par rapport à celui du suiveur SN-RBFN (d’environ 35ms).
D’un autre côté, la figure (IV.18.c) confirme les bonnes performances statiques et
dynamiques du contrôleur ADALINE de vitesse pour un profil aléatoire du vent.
Comme l’indique la figure (IV.18.d), il offre une erreur de poursuite de consigne
considérablement réduite en comparaison avec le contrôleur PI classique. On constate
ainsi de la figure (IV.18.e) que le couple électromagnétique suit parfaitement la
trajectoire de sa référence avec une erreur statique quasi-nulle.
376
374
PV voltage (V)
372
370
368
366
IncCond method
364 proposed SN-RBFN
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7
t(s)
(a)
13
Wind speed (m/s)
12.5
12
11.5
11
10.5
10
0 0.5 1 t(s) 1.5 2 2.5
(b)
Rotating speed wr(rad/s)
wrref
wr-PI
50
Wr-ADALINE
45
40
0 0.5 1 1.5 2 2.5
t(s)
(c)
speed tracking error (rad/s)
PI controller
-1 ADALINE controller
-100
-120
Te
Te*
-140
112
Chapitre IV : Commande neuronale d’un micro-réseau (CA/CC) hybride
4
x 10
5
Pinverter with PI controller
PLac
3.5
2.5
0 0.2 0.4 0.6 0.8 1 1.2
4 t(s)
x 10
0
Reactive power (Var)
-1
-2
-3
Qinverter with PI controller
Qinverter with ENN
-4 QLac
Figure IV.18 : Evolution de : a) Tension PV, b) Profil du vent, c) Vitesse de rotation du GE, d)
Erreur statique de vitesse, e) Couple de la GSAP, f) Puissances active et reactive de l’onduleur.
[Link]
La contribution principale de l’étude, réalisée dans ce chapitre, porte sur la
conception d’un système de commande/contrôle d’un MR hybride en se basant sur
des réseaux de neurones adaptatifs. Le choix des RNAs a été motivé par le souhait de
construire une structure de commande capable de réagir rapidement et
dynamiquement au changement du comportement du MR. Donc, des contrôleurs
neuronaux ont été adoptés pour assurer la poursuite des PPMs des générateurs
renouvelables ainsi que le contrôle de l’échange d’énergie entre le MR, les charges et
le réseau électrique. Un système flou de gestion d’énergie a été également proposé
dans le but de minimiser l’énergie soutirée du réseau principal. Les avantages
principaux de la structure de commande proposée résident dans l’adaptativité, la non
nécessité de connaissance du modèle du système commandé, le faible coût de calcul
et la facilité d’implémentation. Des tests de simulation pour différents scénarios ont
été effectués dans l’environnement Matlab/Simulink. Les résultats obtenus ont
affirmé le gain apporté par les neuro-contrôleurs en termes de performances
dynamiques et statiques par rapport aux contrôleurs classiques.
Dans le chapitre suivant, nous traitons la problématique de réglage des gains des
régulateurs PI, utilisés au sein des schémas de contrôle classiques, en faisant recours à
la LF et aux ANFISs.
113
Chapitre V
V.1. Introduction
Comme déjà mentionné, les convertisseurs de puissance CC-CC et CC-CA jouent un
rôle clé dans les SEHs. Ils permettent d’incorporer des algorithmes de commande afin
d’assurer le contrôle des dynamiques non linéaires des ressources énergétiques de
différents types. Généralement, les schémas de commande classiques dédiés au pilotage
des convertisseurs statiques (onduleurs, redresseurs et hacheurs) utilisent des
régulateurs PI linéaires. Ces derniers peuvent garantir des bonnes performances sous
certaines conditions de fonctionnement. Néanmoins, en présence des non linéarités et de
différents types de perturbations dans le SEH, ces contrôleurs linéaires deviennent
moins efficaces. Les alternatives les plus populaires [3,11,83,175-177], pour s’affranchir
de la linéarité des contrôleurs PI, consistent à ajuster leurs paramètres à l’aide des SIFs,
des ANFISs , de l’algorithme PSO , etc.
Ce chapitre est divisé en deux parties : dans la première nous nous intéressons à
l’étude d’un SEH en appliquant une stratégie de commande intelligente. Notre choix
s’est porté sur la régulation floue des paramètres proportionnel et intégral des
correcteurs PI utilisés dans la structure de commande adoptée. L’étude des
performances du schéma de commande est réalisée dans l’environnement
Matlab/Simulink. La deuxième partie de ce chapitre est consacrée à l’étude d’un
système PV-éolien hybride connecté au réseau électrique. L’objectif est d’établir un
schéma de commande numérique pour l’implémenter sur une carte FPGA. Les résultats
de simulation sont fournis afin de valider le circuit de commande développé.
114
Chapitre V : Commande à FV orienté des systèmes PV hybrides
115
Chapitre V : Commande à FV orienté des systèmes PV hybrides
référence Vref du GPV (VPV) sera perturbée. Si la puissance à l’itération actuelle (PPV(k)) a
augmentée par rapport à celle de l’instant précédent (PPV (k-1)), alors on conserve la
même direction de perturbation de Vref avec un pas constant (ΔVref). Dans le cas
contraire (i.e. ΔPPV négative), on inverse le sens de perturbation de Vref pour converger
vers le PPM. Le processus de perturbation de la tension VPV se répète périodiquement
jusqu’à l’atteinte du PPM.
Début
ΔPPV=PPV(k)-PPV(k-1)
ΔVPV=VPV(k)-VPV(k-1)
ΔPPV=0
Non Oui
ΔPPV>0
ΔVPV>0 ΔVPV>0
Non Oui Oui Non
PPV(k-1)=PPV(k)
VPV(k-1)=VPV(k)
Retour
116
Chapitre V : Commande à FV orienté des systèmes PV hybrides
Z-1
𝑤𝑗 𝑤0
CX 1
Figure V.3 : Schéma bloc de commande à base d’ENN du convertisseur interfaçant le GPV.
On donne dans ce qui suit les équations décrivant la tâche de chaque couche du
contrôleur ENN adopté :
Les neurones de la CE permettent de faire passer les signaux d’entrée vers la CH
:
(1) (1) (1) (1)
𝑧𝑖 (𝑘) = 𝑆𝑖 (𝑛𝑖 ) = 𝑛𝑖 (𝑘) = 𝐸𝑖 (𝑘) , 𝑖 = 1,2
{ (V.1)
𝐸1 (𝑘) = ∆𝐼𝑃𝑉 (𝑘) , 𝐸2 (𝑘) = ∆𝑃𝑃𝑉 (𝑘)
(1) (1)
Avec 𝐸𝑖 , 𝑆𝑖 et 𝑧𝑖 présentent, respectivement, l’entrée, la fonction d’activation et la
(1)
sortie du iième neurone de la CE. 𝑛𝑖 est le niveau d’activation du nœud i.
La sortie d’un neurone caché j se calcule suivant l’équation suivante :
(2) (2) (2) (2) (3) (1)
𝑧𝑗 (𝑘) = 𝑆𝑗 (𝑛𝑗 ) = 𝑆𝑗 (∑𝑟 𝑧𝑟 (𝑘) + ∑𝑖 𝑤𝑖𝑗𝑐 𝑧𝑖 (𝑘)) , 𝑖 = 1,2 , 𝑗 𝑒𝑡 𝑟 = 1,2,3 (V.2)
(2) (2) (2)
𝑧𝑗 , 𝑆𝑗 et 𝑛𝑗 dénotent, respectivement, la sortie, la fonction de transfert et le niveau
d’activation du neurone j de la CH, 𝑤𝑖𝑗𝑐 est le poids de connexion entre le neurone i de
la CE et le neurone j de la CH et 𝑧𝑟(3) (𝑘) est la sortie du neurone r de la CX. On note que
nous avons choisi d’utiliser une fonction d’activation non linéaire de type tangente
hyperbolique pour les neurones de la CH.
Les sorties des neurones de la CH sont renvoyées aux entrées des neurones de la
CX, tel que :
(3) (2)
𝑧𝑟 (𝑘) = 𝑧𝑗 (𝑘 − 1) , 𝑗 𝑒𝑡 𝑟 = 1,2,3 (V.3)
En fin, la sortie d’ENN se calcule comme suit :
(4) (4) (4) (4) (2)
𝑧𝑜 (𝑘) = 𝑆𝑜 (𝑛 ) = 𝑆𝑜 (𝑤0 + ∑𝑗 𝑤𝑗 𝑧𝑗 (𝑘)) = 𝑉𝑃𝑉𝑟𝑒𝑓 , 𝑗 = 1,2,3 (V.4)
(4)
𝑆𝑜 estla fonction de transfert sigmoïdale du neurone de sortie, 𝑤𝑗 représente le poids
synaptique reliant le neurone j de la CH au neurone de sortie et 𝑤0 est le biais du
neurone de sortie. La dérivée de la fonction d’activation de la CS est donnée par :
117
Chapitre V : Commande à FV orienté des systèmes PV hybrides
(4)
(4) −𝛽.𝑛
(4) (4) (4) 1 𝜕𝑧𝑜 𝛽𝑒 (4) (4)
𝑧𝑜 = 𝑆𝑜 (𝑛 )= (4) ⇒ (4) = (4) 2
= 𝛽𝑧𝑜 [1 − 𝑧𝑜 ] (V.5)
−𝛽.𝑛 𝜕𝑛 −𝛽.𝑛
1+𝑒 (1+𝑒 )
Tel que β est une constante positive qui représente la pente de la courbe sigmoïdale.
Dans la présente étude, l’algorithme de RP avec momentum [18,20,101] est adopté
pour accomplir l’apprentissage en-ligne du ENN. Les lois d’adaptation des poids des
connexions selon la règle de la méthode de la descente du gradient sont exprimées
comme suit :
𝜕𝜎𝑒 (𝑘)
𝑤𝑖𝑗𝑐 (𝑘 + 1) = 𝑤𝑖𝑗𝑐 (𝑘) − 𝜇 𝑐 + 𝛼 (𝑤𝑖𝑗𝑐 (𝑘) − 𝑤𝑖𝑗𝑐 (𝑘 − 1))
𝜕𝑤𝑖𝑗
𝜕𝜎𝑒 (𝑘)
𝑤0 (𝑘 + 1) = 𝑤0 (𝑘) − 𝜇
𝜕𝑤0
+ 𝛼(𝑤0 (𝑘) − 𝑤0 (𝑘 − 1)) (V.6)
𝜕𝜎𝑒 (𝑘)
𝑤𝑗 (𝑘 + 1) = 𝑤𝑗 (𝑘) − 𝜇 + 𝛼 (𝑤𝑗 (𝑘) − 𝑤𝑗 (𝑘 − 1))
{ 𝜕𝑤𝑗
tel que, 𝜔𝑐 est la fréquence de coupure du FPB. En calculant les gradients de la fonction
coût dans le système d’équation (V.6) (voir annexe B), on obtient :
𝛽 𝜕𝑒𝑆 (4) (4) (2) 2 (1)
𝑤𝑖𝑗𝑐 (𝑘 + 1) = 𝑤𝑖𝑗𝑐 (𝑘) − 𝜇. 𝑤𝑗 . 𝑧𝑜 (1 − 𝑧𝑜 ) [1 − (𝑧𝑗 ) ] 𝑧𝑖 𝑒𝑆 (𝑘) + 𝛼 (𝑤𝑖𝑗𝑐 (𝑘) − 𝑤𝑖𝑗𝑐 (𝑘 − 1))
2 𝜕𝑧𝑜(4)
𝜕𝑒𝑆 (4) (4)
𝑤0 (𝑘 + 1) = 𝑤0 (𝑘) − 𝜇. 𝛽 (4) 𝑧𝑜 (1 − 𝑧𝑜 )𝑒𝑆 (𝑘) + 𝛼(𝑤0 (𝑘) − 𝑤0 (𝑘 − 1))
𝜕𝑧𝑜
𝜕𝑒𝑆 (4) (4) (2)
𝑤𝑗 (𝑘 + 1) = 𝑤𝑗 (𝑘) − 𝜇. 𝛽 (4) 𝑧𝑜 (1 − 𝑧𝑜 )𝑧𝑗 𝑒𝑆 (𝑘) + 𝛼 (𝑤𝑗 (𝑘) − 𝑤𝑗 (𝑘 − 1))
{ 𝜕𝑧𝑜
(V.10)
(4)
En approximant 𝜕𝑒𝑆 ⁄𝜕𝑧𝑜
à 1, les règles finales d’ajustement en-ligne des paramètres
du suiveur ENN sont données comme suit :
118
Chapitre V : Commande à FV orienté des systèmes PV hybrides
Product1 -1 W1(4)
1 Ei(1) X(2) 1 Add2
Xi(1) 1
2 Ei(1) Ei(2)
-K-
Xi(1) X(2) 2 z
Add8 momentium3
W2(4)
1 Unit Delay7
Ei(2) input layer X(2) 1
X(3)r
X(2) z 1
Unit Delay1 W3(4)
X(2) 3
Wjo(3) x(2) X(2) 1 Wjo(3)
Xo(4) 1
-2
en en Wjo(3) Wo io Xo(4) Product2
Add1
X(2) x(3) r 1 -K-
X(4) X(4) z W0(4)
output layer Unit Delay6
Add7 momentium2
1
OL Wjo
context layer hidden layer z
Unit Delay2
1
2
en Product4 -3
Add3
1 1 -K-
Constant -K- z
Add6 momentium1
num(s) Add4 Unit Delay5
eu
1
3 0.8 1 Product1 Product6 Gain
z
den(s) f(u) f(u) 3
en u learning Unit Delay8
FPB Gain1 Exp1 1 X(4) 1
Sum1 Reciprocal1 Fcn8 Fcn9
Product5 -5
one1 Add10
1 -K-
z
Add11 momentium6
Unit Delay9
(a) (b)
Figure V.4 : Implémentation dans Matlab/Simulink de : a) Contrôleur ENN, b) La règle d’adaptation
des poids de la CS.
V.2.2.2. Contrôleur flou de la PAC-OS
Le hacheur survolteur associé à la PAC-OS est utilisé pour élever sa tension de sortie
ainsi que contrôler l’énergie générée par le stack. En fonction des conditions de
génération et de demande d’énergie, la puissance produite par la PAC est ajustée en
régulant son courant de sortie (Ifc) comme le décrit le schéma de la figure (V.5). La
consigne du courant de la PAC-OS (Ifcref) est calculée en divisant la puissance de
référence (Pfcref), reçue du superviseur central, par la tension mesurée aux bornes du
stack (Vfc). La référence obtenue est ensuite limitée dans un intervalle qui correspond à
la gamme désirée du facteur d’utilisation de combustible Uf (Uf ∈[0.8,0.9] ) [31]:
𝑖𝑛 𝑖𝑛
0.8 𝑞𝐻 0.9 𝑞𝐻
𝐼𝑓𝑐_𝑚𝑖𝑛 =
2𝐾𝑟
2
≤ 𝐼𝑓𝑐𝑟𝑒𝑓 ≤ 𝐼𝑓𝑐_𝑚𝑎𝑥 =
2𝐾𝑟
2
(V.12)
Ensuite, l’écart entre le courant mesuré et sa référence est calculé puis, envoyé au
contrôleur PI-flou adaptatif. La sortie du correcteur du courant est comparée avec un
signal triangulaire de fréquence haute pour générer le signal de contrôle du dispositif de
commutation du convertisseur élévateur (voir figure (V.5)).
119
Chapitre V : Commande à FV orienté des systèmes PV hybrides
Ifmin
3 fuzzy_gain_tuner
Umax
Pfc_ref
Satur1 Ifc_ref 1 Kp
f(u)
max error
z >= 1
Product3 limit Unit Delay Ki Product2
1 1/u Saturation PWM
qH2 Add2 Relational
Vfc Operator1
Satur Fcn1 4 -K- -K- 2 K Ts
Ifcmax Ifc
1/2Kr min z-1 triangle
Product1
120
Chapitre V : Commande à FV orienté des systèmes PV hybrides
Avec (KPf, KIf) représentent les valeurs finales des paramètres PI. ΔKPi et ΔKIi dénotent
les valeurs initiales des gains du contrôleur PI. Ces derniers sont déterminés en utilisant
la méthode classique de placement de pôles.
NB NMNB NS NM
ZE PS NS
PM PB B
1 1 S MS MB
Degree of membership
Degree of membership
0.8 0.8
0.6 0.6
0.4 0.4
0.2 0.2
0 0
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Error and change of error (A) Kp and Ki
(a) (b)
Figure V.6: FAs assignées aux : (a) Variables d’entrée (efc et Δefc), (b) Variables de sortie (KP,KI) du RGF
du contrôleur de PAC-OS.
efc Δefc NB NM NS ZE PS PM PB
NB B B B B B B B
NM MB MB MB MS MB MB MB
NS MB MS MS MS MS MS MB
ZE MS MS S S S MS MS
PS MB MS MS MS MS MS MB
PM MB MB MB MS MB MB MB
PB B B B B B B B
Tableau V.1 : Bases de règles floues pour la régulation du gain KP du contrôleur du courant de PAC.
efc Δefc NB NM NS ZE PS PM PB
NB S MB MB MB MB MB S
NM MB MS MS MB MS MS MB
NS MB MB B B B MB MB
ZE B B B B B B B
PS MB MB B B B MB MB
PM MB MS MS MB MS MS MB
PB B MB MB MB MB MB B
Tableau V.2 : Bases de règles floues pour la régulation du gain KI du contrôleur du courant de PAC.
121
Chapitre V : Commande à FV orienté des systèmes PV hybrides
NOT 2
2
Kp PWM1
ibatt Product2 Zero-Order Logical
Error Hold Operator2
3 >=
Ki
1
Pbref Saturation3
Saturation2 PWM
Add3 K Ts Relational
Product3 FGT
Operator1
1 1/u z-1
Product1
Vb
Saturation1 Fcn1
Triangle
122
Chapitre V : Commande à FV orienté des systèmes PV hybrides
Le calcul des puissances instantanées (P,Q) en fonction de (Ψgα, Ψgβ) est alors accompli
selon les formules suivantes [171]:
3
𝑃 = 𝜔 ∙ (Ψ𝑔𝛼 𝑖𝑐𝛽 − Ψ𝑔𝛽 𝑖𝑐𝛼 )
2
{ 3 (V.15)
𝑄 = 𝜔 ∙ (Ψ𝑔𝛼 𝑖𝑐𝛼 + Ψ𝑔𝛽 𝑖𝑐𝛽 )
2
L’équation modélisant l’ensemble « onduleur 3L-NPC - filtre RL – réseau » dans le
repère dq (en négligeant la chute de tension aux bornes de la résistance du filtre) est
donnée par :
𝑑𝑖𝑐𝑑
𝑣𝑑 = 𝐿 + 𝐿𝜔𝑖𝑐𝑞 + 𝑣𝑔𝑑
𝑑𝑡
{ 𝑑𝑖𝑐𝑞 (V.16)
𝑣𝑞 = 𝐿 − 𝐿𝜔𝑖𝑐𝑑 + 𝑣𝑔𝑞
𝑑𝑡
La commande vectorielle par orientation du FV est achevée en alignant l’axe d du
repère tournant dq avec le vecteur du FV du réseau (voir figure (IV.7)). En conséquence,
le vecteur de tension du réseau est aligné avec l’axe q en quadrature avec le vecteur du
FV [141]. Donc, dans le schéma de VFOC, la puissance active de l’onduleur est contrôlée
en agissant sur la composante en quadrature du courant (𝑖𝑐𝑞 ). Alors que, la puissance
réactive est régulée en altérant la composante directe du courant d’onduleur (𝑖𝑐𝑑 )
[141,167].
Comme pour la structure de VOC, les composantes d et q de la tension de référence
de l‘onduleur (𝑣∗𝑐𝑑 et 𝑣∗𝑐𝑞 ), dans le schéma de VFOC, sont déterminées en fonction des
termes de découplage comme suit [141,167]:
𝑣∗𝑐𝑑 = 𝑣𝑟𝑑 + 𝐿𝜔𝑖𝑐𝑞 + 𝑣𝑔𝑑
{ ∗ (V.17)
𝑣𝑐𝑞 = 𝑣𝑟𝑞 − 𝐿𝜔𝑖𝑐𝑑 + 𝑣𝑔𝑞
Tel que 𝑣𝑟𝑑 et 𝑣𝑟𝑞 désignent les signaux de sortie des contrôleurs des courants (𝑖𝑐𝑑 et
𝑖𝑐𝑞 ) dans le repère tournant dq. En utilisant la transformée de Park, les composantes du
vecteur de tension du réseau 𝑣𝑔𝑑 et 𝑣𝑔𝑞 sont estimées, à chaque période d’échantillonnage
Ts, en fonction des composantes du FV comme suit [168] :
Ψ𝑔𝛼 (𝑘+1)−Ψ𝑔𝛼 (𝑘)
𝑣𝑔𝛼 (𝑘) = 𝑣𝑔𝑑 𝑣𝑔𝛼
𝑇𝑠
{ Ψ𝑔𝛽 (𝑘+1)−Ψ𝑔𝛽 (𝑘)
⟹ [𝑣 ] = [𝑇𝛼𝛽→𝑑𝑞 (𝜃𝑔 )] [𝑣 ] (V.18)
𝑔𝑞 𝑔𝛽
𝑣𝑔𝛽 (𝑘) =
𝑇𝑠
123
Chapitre V : Commande à FV orienté des systèmes PV hybrides
En fin, les références des courants selon les axes d et q (𝑖𝑐𝑑,𝑟𝑒𝑓 , 𝑖𝑐𝑞,𝑟𝑒𝑓 ) sont obtenues en
appliquant la transformée de Park (𝑇𝛼𝛽→𝑑𝑞 ) comme suit :
𝑖𝑐𝑑𝑟𝑒𝑓 𝑖𝑐𝛼,𝑟𝑒𝑓
[
𝑖𝑐𝑞,𝑟𝑒𝑓
] = [𝑇𝛼𝛽→𝑑𝑞 (𝜃𝑔 )] [
𝑖𝑐𝛽,𝑟𝑒𝑓
] (V.21)
La figure (V.8) indique que des signaux de correction issus du CFMG de la tension du
bus CC (noté ∆𝑖𝑞𝑟𝑒𝑓 ) et du FPB (noté 𝑖𝑞ℎ ) sont utilisés pour corriger la valeur de la
consigne du courant en quadrature ( 𝑖𝑐𝑞,𝑟𝑒𝑓 ). Le FPB d’ordre un, rajouté à la boucle de
régulation du courant en quadrature, a pour but de réduire les oscillations statiques du
courant [170,181]. D’autre part, la sortie du régulateur PI de la puissance réactive ∆𝑖𝑑𝑟𝑒𝑓
est rajoutée à la composante d du courant de référence ( 𝑖𝑐𝑑,𝑟𝑒𝑓 ) dans le but d’améliorer
les performances de contrôle. Les composantes d et q du courant mesuré sont comparées
avec leurs consignes correspondantes et les erreurs sont ensuite fournies aux entrées des
Régulateurs de Gain Neuro-Flou (RGNF) des contrôleurs PI des courants d et q. Enfin,
les tensions de référence obtenues (𝑣∗𝑐𝑑 et 𝑣∗𝑐𝑞 ) sont utilisées pour générer les signaux de
contrôle du convertisseur 3L-NPC.
gα,gβ Estimation 𝒗𝒈𝒅
Q
θg de vgd et vgq 𝒗𝒈𝒒
Qref - +
PI
icα , icβ 𝒗𝒈𝒅
icd Δidref KP et KI
Sa Calcul des θg αβ icd*
Sb composantes 𝒗𝒓𝒅 + 𝒗∗𝒄𝒅
Sc αβ du FV
+
icdref + RGNFs
dq icq PI
VDC gαgβ - + -
∫e
Pref , Qref icd
ωL Sa
icα gα,gβ
Calcul de +
icabc abc SPWM Sb
θg icdref et icqref icq
+
ωL Sc
αβ icβ + - icq*
s + +
+
VDCref + Surface de icqref
- RGNFs PI
glissement
CFMG - 𝒗∗𝒄𝒒
- Δiqref + + + 𝒗𝒓𝒅 +
+ ∫e
VDC ds FPB iq,h KP et KI
𝒗𝒈𝒒
- +
+
Figure V.8 : Schéma de VFOC adopté pour la commande de l’onduleur.
124
Sous forme discrète, la surface de commutation S(k) et sa variation dS(k) sont
calculées comme suit :
𝑆(𝑘) = 𝜆𝑖 ∑𝑘𝑖=1 𝑒𝑑𝑐 (𝑖) + 𝜆𝑒 𝑒𝑑𝑐 (𝑘) + 𝑑𝑒𝑑𝑐 (𝑘)
{ (V.23)
𝑑𝑆(𝑘) = 𝑆(𝑘) − 𝑆(𝑘 − 1)
𝑒𝑑𝑐 (𝑘) = 𝑉𝐷𝐶𝑟𝑒𝑓 (𝑘) − 𝑉𝐷𝐶 (𝑘)
Tel que { (V.24)
𝑑𝑒𝑑𝑐 (𝑘) = 𝑒𝑑𝑐 (𝑘) − 𝑒𝑑𝑐 (𝑘 − 1)
Avec edc(k) est l’erreur actuelle de la tension du bus CC, dedc(k) est la dérivée discrète de
l’erreur de tension et VDCref est la valeur de référence de la tension du bus CC. 𝜆𝑖 et
𝜆𝑒 sont des constantes positives non nulles. Comme le montre la figure (V.9), la loi de
commande en sortie du CFMG est calculée selon la formule suivante [183] :
𝑈(𝑘) = 𝑑𝑈(𝑘) + 𝑈(𝑘 − 1) (V.25)
dU est l’incrément de l’action de commande. Le CFMG développé est donc basé sur le
modèle flou de TSK d’ordre-zéro. Les FAs trapézoïdale et gaussienne assignées aux
entrées du CFMG sont illustrées par la figure (V.10). La distribution des sorties
singletons du CFMG sur l’univers de discours est comme suit : NB=-1.1, NM=-0.75,
NS=-0.4, ZE=0, PS=0.4, PM=0.75 and PB=1.1. L’ensemble des règles floues du CFMG est
donné par le tableau (V.3). La méthode d’inférence de « Max-Min » est également
adoptée. En outre, la sortie du CFMG est calculée par la méthode de moyenne pondérée.
Donc, l’objectif de commande du CFMG est de générer le signal de contrôle correct
qui garantit la convergence rapide vers l’état d’équilibre et la stabilité asymptotique du
système. Quand 𝑆 > 0 et 𝑑𝑆 > 0 (voir tableau (V.3)), l’incrément du signal de contrôle
𝑑𝑈 doit être négatif afin de conduire la trajectoire de la variable d’état (VDC) vers la
surface de glissement dans un temps fini. Toutefois, en régime de glissement 𝑆 ≈ 0 et
𝑑𝑆 ≈ 0, la perturbation de la loi de commande (𝑑𝑈) du CFMG doit être quasiment nulle
afin de maintenir la trajectoire d’état sur la surface de glissement et ne la quitte jamais.
Memory
1 -K- -K- dU
S
edc(t) -K- 1
c1 normaliser Satur1
1 U(k)
Add normaliser2
1
-K- Fuzzy Logic z
z dS Controller Unit Delay1
1
Unit Delay normaliser1 Satur
K Ts z
-K-
Unit Delay2
z-1
Discrete integrator c2
0.8
Degree of membership
0.8
0.6 0.6
0.4 0.4
0.2 0.2
0 0
(a) (b)
Figure V.10 : FAs des entrées du CFMG : a) Surface de glissement S, b) Sa variation dS.
Chapitre V : Commande vectorielle à FV orienté des systèmes PV hybrides
S dS NB NM NS ZE PS PM PB
NB PB PB PB PM PM PS PS
NM PB PM PM PM PS PS NS
NS PM PS PS PS PS NS NM
ZE PM ZE ZE ZE ZE ZE NM
PS PM PS NS NS NS NS NM
PM PS NS NS NM NM NM NB
PB NS NS NM NM NB NB NB
126
Chapitre V : Commande vectorielle à FV orienté des systèmes PV hybrides
construire une base de données d’apprentissage. Les paramètres des FAs des entrées et
de sortie de chaque ANFIS sont modifiés hors-ligne en utilisant l’algorithme hybride
d’entraînement (voir section II.3.3). La figure (V.11) illustre les FAs des entrées des
RGNFs de KP et KI obtenues après la phase d’apprentissage.
(a)
(b)
Figure V.11 : FAs des entrées des régulateurs des gains KP et KI de contrôleur de courant direct.
127
Chapitre V : Commande vectorielle à FV orienté des systèmes PV hybrides
Tant que le SOC est inférieur au niveau minimal (𝑆𝑂𝐶 < 𝑆𝑂𝐶𝑚𝑖𝑛 ), le SSEB doit
opérer en mode de charge sans tenir compte les conditions de génération et de
demande d’énergie.
Quand la puissance générée par la source PV est supérieure à la puissance requise
par la charge CA (Pnet<0), l’excès d’énergie est stocké dans le SSEB tant que le 𝑆𝑂𝐶 ≤
𝑆𝑂𝐶𝑚𝑎𝑥 pour éviter la surcharge des batteries.
Dans le cas où le SSEB est complètement chargé (𝑆𝑂𝐶 > 𝑆𝑂𝐶𝑚𝑎𝑥 ), dans une période
de surproduction, la puissance PV excédentaire est injectée totalement au réseau
principal.
Si l’extra de puissance demandée (Pnet) surpassant la puissance disponible du GPV
est moyen ou élevé, la PAC-OS doit répondre aux exigences des charges.
D’autre part, si le 𝑆𝑂𝐶 < 𝑆𝑂𝐶𝑚𝑖𝑛 durant la période d’insuffisance de génération (0
<Pnet ≤ Pfcmax), la PAC-OS génère plus d’énergie afin de charger le SSEB tout en
satisfaisant la demande de charge CA.
En revanche, si la demande énergétique surpasse la puissance maximale de la PAC
(Pnet>Pfcmax), le banc de batteries intervient afin de restituer le manque d’énergie (tel
que 𝑆𝑂𝐶 ≥ 𝑆𝑂𝐶𝑚𝑖𝑛 ).
Lorsque la PAC-OS alimente les charges locales, le SSEB est utilisé pour lisser la
puissance produite.
Par ailleurs, le réseau principal doit fournir l’énergie requise, si l’énergie totale
produite par le MR ne peut pas satisfaire la demande des charges CA (PMR<Plaod).
Le schéma bloc du SGE implémenté dans l’environnement Matlab/Simulink est
illustré par la figure (V.12.b).
128
Pload PPV
Pnet= Pload-PPV
Pfcref et Pbref représentent les puissances de références pour la PAC-OS et le SSEB, respectivement. Pfc et Pb sont les puissances estimées en sortie de PAC-OS et du SSEB,
respectivement. PMR et PG désignent les puissances du MR et du réseau principal, respectivement. Pbmax et Pbmin représentent, respectivement, les valeurs maximale et minimale
de puissance de charge du SSEB. Pfcmax est la puissance maximale de la PAC-OS.
(a)
129
Chapitre V : Commande à FV orienté des systèmes PV hybrides
meet
Pf cref
Hsoc -C-
Msoc
Pbmax 1
< 40 NOT Lsoc
Pbref
Product1 Pfcref
Compare
AND Pf c Subtract3
Logical1 5 -K-
To Constant NOT
Logical MET DEMAND
1 > 80 Operator2 1
Pfc Gain1 2
SOC
Logical Pnet 4 Product3
Pf cref Pfcref
C excess Pbref
[Pnet]
Subtract7 Lsoc AND 0 Subtract4
Hsoc 1
== 0
Msoc
Goto2 Pbref 2 Pbre
C1 Lsoc
meet Logical Product4
Pbref
Subtract4 Operator1
2 4 excess
Pload
<0
Pb
Pbest OR
3 Pf cest AND
5 Subtract2 2
Co
[PL] Pfc Pload
Pf cref Product2
Ppv [PL]
Logical 0
>0 def icit2 Hsoc Logical Operator8
Goto21 Pnet Operator2
Com zero2
Hsoc 3
Msoc AND
Pbref Msoc
Lsoc
Ppv Logical
deficit Operator6
(b)
Figure V.12: SGE adopté: a) Organigramme de sa stratégie de supervision, b) Son implémentation dans
Matlab/Simulink.
V.2.3. Résultats de simulation
Afin de prouver l’efficacité de la stratégie de commande proposée dans cette partie,
différents tests de simulation sont effectués à l’aide de l’outil Matlab/Simulink. Le MR est
étudié pour des variations des conditions climatiques et des changements des exigences
énergétiques. Dans cette étude, les valeurs minimale et maximale du SOC de SSEB sont
fixées, respectivement, à 40%et 80%. Les paramètres de simulation sont donnés dans
l’annexe. Les cas suivants sont donc considérés :
A. Conditions variables de demande de charge
Dans le premier cas de simulation, les conditions climatiques sont stables (une
température de T=25°C et un ensoleillement de G=1000W/m2) et le niveau de SOC est
initialement de l’ordre de SOCinitial=70%. La demande de la charge ohmique-inductive
varie à l’instant t=5s de (PL,QL) = (100kW, 10kVAr) à (PL,QL) = (95kW, 20kVAr) puis change
à t=10s vers (PL,QL) =(120kW, 5kVAr) et finalement stabilise dans l‘intervalle de temps
[15s,20s] à (PL,QL)=(78kW,5kVAr). Les résultats obtenus dans ce cas sont illustrés par les
figures (V.13) , (V.14) et (V.15).
La figure (V.13.a) montre qu’en adoptant le contrôleur MPPT proposé, le point de
fonctionnement optimal de la source PV est atteint après un temps de réponse si rapide
d’environ 0.0335s qui est plus court que celui du régulateur P&O (environ 0.138s). De plus,
les oscillations statiques autour du PPM, obtenues avec le contrôleur ENN, sont
considérablement réduites par rapport à celles du suiveur P&O (voir figure (V.13.d)).
La figure (14.b) montre que durant la période [0,5s], le SSEB est inactif puisque la
puissance produite par le GPV et la PAC-OS satisfait la demande totale de la charge CA.
A l’instant t=5s, la demande d’énergie reste plus élevée que la puissance crête du GPV
(Pload> PPV). En conséquence, le stack de PAC-OS continue à fournir l’énergie requise (voir
figure (V.14.b)). À l’instant t=10s, le SSEB commence à fonctionner en mode de décharge
130
Chapitre V : Commande à FV orienté des systèmes PV hybrides
afin de satisfaire le besoin d’énergie (voir figure (V.14.a) et (V.14.b)). Durant la période de
temps [15s,20s], le GPV et la PAC-OS s’occupent de la réponse à la demande de charge
CA. Alors que, le SSEB opère en mode de charge. La puissance active soutirée du réseau
principal durant la période de test [0,20s] est constamment nulle puisque la production
d’énergie par le MR satisfait entièrement les exigences d’énergie. En outre, la figure
(V.14.c) montre que l’onduleur a compensé la puissance réactive de la charge CA.
500
PV voltage (V)
400
300
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9
t(s)
(a)
PV current (A)
120
115
110
105
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9
t(s)
(b)
500
Vref (V)
450
400
4
4
P&O
x 10
2
5.06
ENN
5.05
0 5.04
t(s) 5.02
(d) 5.01
P&O
5
ENN
4.99
4.98
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
t(s)
Figure V.13 : Variation de : a) Tension PV, b) Courant PV, c) Tension de référence imposée par le ENN,
d) Puissance PV de sortie.
131
Chapitre V : Commande à FV orienté des systèmes PV hybrides
4
x 10
10
5
Pinv
0 Pload
Pgrid
0 2 4 6 8 10 12 14 16 18 20
t(s)
x 10
4 (a)
6
Pfc et Psseb (W)
4
Pb
2 Pbref
Pfc
0
Pfcref
-2
0 2 4 6 8 10 12 14 16 18 20
t(s)
4
(b)
Qinv, Qgrid, Qlaod (VAr)
x 10
Qgrid
4
Qinv
2 Qload
-2
0 2 4 6 8 10 12 14 16 18 20
t(s)
(c)
Iq-PI
300 Iqref
Iq-RGNF
200
100
0 2 4 6 8 10 12 14 16 18 20
t(s)
(a)
450 Kp
Kp and Ki
Ki
400
350
300
0 2 4 6 8 10 12 14 16 18 20
t(s)
(b)
DC-link -voltage (V)DC-link -voltage (V)
1500
1000
Vdc-PI
500
1600
Vdc-CFMG
1500
Vdc-PI
1400
Vdc-CFMG
1300
0
1200
0 2 4 6 8 10 12 14 16 18 20
t(s)
1100
1000
900 (c)
800
700
600
0 0.5 1 1.5
t(s)
(c)
132
Chapitre V : Commande à FV orienté des systèmes PV hybrides
D’autre part, les résultats représentés dans la figure (V.15) démontrent que le schéma
de VFOC a garanti des bonnes performances de suivi de trajectoire de référence malgré la
variation continue des consignes des puissances active et réactive. La figure (V.15.a)
montre qu’en utilisant les RGNFs du contrôleur de courant q, des meilleures réponses
transitoire et statique sont obtenues en comparaison avec le contrôleur PI. Après chaque
changement de consignes, les paramètres des régulateurs PI des courants d et q sont
adaptés en temps réel comme l’illustre la figure (V.15.b). En outre, il est clair de la figure
(V.15.c) que le CFMG a surperformé le contrôleur PI de tension en assurant un
dépassement minimal et une réponse transitoire plus rapide.
B. Conditions variables d’éclairement et de demande d’énergie
Dans le deuxième test de simulation, la demande de la charge ohmique décroît de
PL=110kW à PL=70kW à l’instant t=2.5s puis accroît à PL=135kW à t=4s et change à l’instant
t=10s à PL=120kW (tel que SOCinitial=90%). D’autre part, le niveau d’irradiation solaire (G)
varie à une température constante (T=25°C) comme l’illustre la figure (V.16.a). Les
résultats trouvés dans ce cas sont représentés dans les figures (V.16) et (V.17).
Ces résultats démontrent que le système de commande adopté a réagi rapidement aux
changements d’éclairement et de demande de charge. Le contrôleur ENN adopté assure
des bonnes performances de poursuite pour un éclairement variable comme l’illustre la
figure (V.16). Pour chaque niveau d’irradiation, la tension de référence (Vref) du GPV est
déterminée précisément par le suiveur ENN. Les amplitudes des oscillations de la
puissance PV en régime permanent sont considérablement atténuées en comparaison avec
l’algorithme de P&O (voir figure (V.16.e)).
D’autre part, la figure (V.17) montre que le SGE centralisé détermine correctement les
consignes des puissances de la PAC-OS et du SSEB conformément aux exigences
d’énergie. Durant la période [0s,2.5s], la PAC et le SSEB sont activés afin de répondre à la
demande de charge CA comme le montre les figures (V.17.a) et (V.17.b). Pendant
l’intervalle du temps [4s,10s], l’exigence de charge CA dépasse l’énergie totale produite
par le MR comme l’indique la figure (V.17.c). En conséquence, l’extra de l’énergie requise
est soutiré du réseau électrique. Les sources PAC et SSEB sont donc commandées, durant
cette période, pour fournir le maximum d’énergie. Par ailleurs, on constate des figures
(V.17.a) et (V.17.b) que les courants de sortie de la PAC-OS et du SSEB (IFC, Ib) suivent
précisément les consignes imposées par le SGE (IFC*, Ib*). Les contrôleurs PI basés sur les
RGFs offrent des temps de réponse rapide et des erreurs statiques minimales en
comparaison avec les régulateurs PI à gains fixes (voir (V.17.b)).
133
Chapitre V : Commande à FV orienté des systèmes PV hybrides
1000
G (W/m2)
800
600
400
0 1 2 3 4 5 6 7 8 9 10
t(s)
PV voltage (V) (a)
500
450
400
0 1 2 3 4 5 6 7 8 9 10
t(s)
(b)
500
Vref (V)
450
400
0 1 2 3 4 5 6 7 8 9 10
t(s)
(c)
PV current (A)
150
100
50
0 1 2 3 4 5 6 7 8 9 10
t(s)
(d)
4
x 10
6
P&O
ENN
5
4
PV Power (W)
3 4
x 10
4.5
PV Power (W)
2 4.45
4.4
4.35
1 P&O
4.3 ENN
4.25
7 7.5 8 8.5 9 9.5 10
0 t(s)
0 1 2 3 4 5 6 7 8 9 10
t(s)
(e)
Figure V.16 : Variation de : a) Irradiation solaire, b) Tension PV, c) Tension de référence imposée par le
ENN, d) Courant PV, e) Puissance PV de sortie.
134
Chapitre V : Commande à FV orienté des systèmes PV hybrides
130
Ifc
120
115
0 2 4 6 8 10 12 14 16
t(s)
(a)
Ibref et Ib (A)
80
90
60 Ib-PI 85
40 Ib-RGF
Ibref et Ib (A)
80
Ib-ref 75
20 Ib-PI
Ib-RGF
70 Ibref
0 65
0 2 4 6 8 10 12 14 16
t(s) 60
0 2 4 6 8 10 12
t(s)
4
(b)
Pinv , Pgrid et Pload (W)
x 10 (b)
20
15
10
5 Pinv
0 Pgrid
Pload
-5
0 2 4 6 8 10 12 14 16
t(s)
(c)
135
Chapitre V : Commande à FV orienté des systèmes PV hybrides
sortie de l’onduleur au PCC est synchronisé avec la tension du réseau grâce à l’estimation
correcte de l’angle de phase du FV. Le taux de distorsion harmonique du courant
d’onduleur est acceptable qui est de l’ordre de 2.03%.
36
34
Température (°C)
32
30
28
26
24
22
0 2 4 6 8 10 12 14 16 18 20
t(s)
(a)
4
x 10
4.5
4
5
x 10 P&O
ENN
PV Power (W)
4
4 P&O
PV Power (W)
3 ENN
2
3.5 1
0
0 0.05 0.1 0.15 0.2 0.25
4
t(s) x 10
3 5.05
P&O
PV Power (W)
5 ENN
2.5
4.95
4.9
2 5 6 7 8 9 10 11 12
0 2 4 6 8 10 12 14 t(s) 16 18
t(s)
(b)
Figure V.18. Variation de la : a) Température, b) Puissance de sortie du GPV.
136
Chapitre V : Commande à FV orienté des systèmes PV hybrides
4
x 10
10
6
Pinv
4 Pload
Pgrid
2
-2
0 5 10 15 20 25
t(s)
x 10
4 (a)
Qinv, Qload et Qgrid (VAr)
-5 Qinv
Qload
Qgrid
-10
0 5 10 15 20 25
t(s)
(b)
Figure V.19 : Variation de : a) Puissance active d’onduleur, de charge et du réseau, b) Puissance réactive
d’onduleur, de charge et du réseau.
100
d-axis current (A)
50 Id
Idref
0
-50
-100
0 2 4 6 8 10 12 14 16 18 20
t(s)
(a)
q-axis current (A)
200
150 Iq
Iqref
100
0 2 4 6 8 10 12 14 16 18 20
t(s)
(b)
420
400
Kp and Ki
Kp
380
Ki
360
340
0 2 4 6 8 10 12 14 16 18 20
t(s)
(c)
Figure V.20 : Variation de : a) Courant direct d’onduleur, b) Courant en quadrature d’onduleur, c) Gains
du contrôleur PI de courant.
137
Chapitre V : Commande à FV orienté des systèmes PV hybrides
950
900
850
800
0 2 4 6 8 10 12 14 16 18 20 22
t(s)
4
Pfc, Pfcref, Pb, Pbref (W) x 10 (a)
6
4
Pbref
2 Pb
0 Pfc
Pfcref
-2
0 5 10 15 20 25
t(s)
(b)
Igrid (A) and Vgrd (V)
400
Vgrid
200 Igrid
-200
-400
3.9 4 4.1 4.2 4.3 4.4 4.5
t(s)
(c)
Figure V.21 : Evolution de a) Tension du bus CC, b) Puissance de PAC-OS et du SSEB, c) Courant
d’onduleur et tension du réseau (phase a).
V.2.4. Conclusion 1
Dans cette partie du chapitre, une stratégie de commande intelligente d’un SEH a été
proposée et évaluée dans l’environnement Matlab/Simulink. Des contrôleurs locaux
basés sur les TIAs ont été développés pour le pilotage des convertisseurs CC-CC et CC-
CA. En outre, un superviseur central de puissance a été implémenté afin de gérer
efficacement les flux énergétiques à travers le MR. Les résultats de simulation obtenus ont
confirmé l’efficacité de la technique de commande adoptée pour différents scénarios. Le
SGE centralisé a déterminé précisément les puissances de référence pour la PAC-OS et le
SSEB sous des conditions variables de génération et de demande d’énergie. En
comparaison avec le suiveur P&O classique, le contrôleur ENN adopté a donné des
meilleures performances de poursuite du PPM en régimes transitoire et permanent. De
plus, le schéma de VFOC, basé sur des régulateurs neuro-flous, a assuré des bonnes
performances de régulation et de suivi de consigne. Les puissances fournies par les
sources d’appoints ont été contrôlées dynamiquement avec une adaptation en ligne des
paramètres de leurs contrôleurs PI. Les avantages principaux du schéma de commande
proposé résident dans la non-exigence de modélisation mathématique du système, la
capacité de traitement des non-linéarités et la robustesse vis-à-vis les perturbations.
Comme perspective à ce travail, le système de commande développé doit être implémenté
sur des plateformes FPGA/DSP afin de vérifier expérimentalement leur efficacité et
faisabilité. En outre, une étude du MR en mode d’islanding peut être réalisée où, le SGE
proposé doit être modifié pour pouvoir réagir correctement dans ce mode d’opération.
138
Chapitre V : Commande à FV orienté des systèmes PV hybrides
139
Chapitre V : Commande à FV orienté des systèmes PV hybrides
R2011b ISE14.1
Développement du système de commande dans Matlab/Simulink
Synthèse
Placement et routage
Figure V.22 : Etapes de développement et d’implémentation d’un contrôleur sur une carte FPGA à
l’aide de XSG [185,190].
V.3.2. Topologie du SEH étudié
La configuration du système de production étudié dans cette partie est illustrée par la
figure (V.23). La TE et le GPV sont interfacés au bus CC par le biais d’un redresseur
triphasé et d’un hacheur élévateur, respectivement. L’onduleur triphasé à deux niveaux
s’occupe de la conversion et l’injection de l’énergie totale générée au réseau électrique.
CA-CC CC-CA
R L vga
TE+
GSAP CDC vgb
vgc
Sa,b,c Sa,b,c
isd , isq VDC
iga,b,c
MPPT + FOC VFOC à base de
Vw , Te , ωt Backstepping vga,b,c
+
Convertisseur
GPV VPV
élévateur
-
IPV Spv
VPV P&O-MPPT
140
Chapitre V : Commande à FV orienté des systèmes PV hybrides
141
Chapitre V : Commande à FV orienté des systèmes PV hybrides
0.0140380859375 a
z-3
1 In a -1 ab a
z a pas
Ipv ab Sy stem b a+b a
Gateway In1 b a
a-b Generator b
b not cast a+b
a>b Mult5
Mult1 -1 b
z b AddSub2
AddSub5 Convert2
Inverter
Delay3 Relational AddSub1 Out 1
0 xor
cast a z-1 cast
-3
z Gateway Out1 Vref
2 In a
cts a Logical ab Convert
Vpv a-b Convert3 b Delay1
Gateway In3 -1 b
z a>b -0.0140380859375
0 b Mult3
Delay4 AddSub6
Relational1 pas1
cts1
(a)
a
Delay1 sel
a>b
1 b
-1 d0
Sy stem z Relational1
Generator MAX
a cast a a d1
x2
a a
a+b Convert2 a+b a>b sel Mux1
1 In a Kp b cast b b ab Out 1
a-b b
vpv In3 a-b 0 Gateway Out1 PWM
AddSub4 b Relational2 d0
b x 0.0002441 Convert1 AddSub2 Relational
2 In MIN
-1 AddSub3 d1
vref AddSub1 z
In1 Ki*Ts
Bloc de génération du signal ++ reinterpret
Delay Mux2
MLI de commande Reinterpret
2KHZ
(b)
Figure V.25: Développement à base de XSG de : a) Contrôleur P&O, b) L’ensemble contrôleur PI-
générateur du signal MLI.
𝑻∗𝒆 𝟐 𝒊∗𝒔𝒒 𝝎𝒕
𝑲𝒐𝒑𝒕 × +- PI SGEE
𝟑𝒑𝜳𝒓
𝒊𝒔𝒒
𝝎𝒕 𝟐
D’après l’équation (II.10), la puissance de la turbine Darrieus peut être calculée comme
suit :
1 𝜔𝑡 3
𝑃𝑡 = 𝜌𝐻𝑅4 𝐶𝑃 (𝜆) (V.29)
2 𝜆3
Tel que 𝑆 = 𝑅𝐻 . Pour une vitesse de vent donnée, le rapport de vitesse est optimal au
PPM (i.e 𝜆 = 𝜆𝑜𝑝𝑡 ) et la valeur correspondante du coefficient de puissance est maximale
(𝐶𝑃 = 𝐶𝑃−𝑚𝑎𝑥 𝑒𝑡 𝑃𝑡 = 𝑃𝑡−𝑚𝑎𝑥 ). Dans ces conditions, on obtient un couple éolien optimal
(𝑇𝑡 = 𝑇𝑡−𝑜𝑝𝑡 ) qui s’exprime comme suit :
142
Chapitre V : Commande à FV orienté des systèmes PV hybrides
𝑃𝑡−𝑚𝑎𝑥
𝑇𝑡−𝑜𝑝𝑡 = 𝜔𝑡
= 𝐾𝑜𝑝𝑡 . 𝜔𝑡 2
{ 1 𝐶𝑃−𝑚𝑎𝑥 (𝜆) (V.30)
Tel que 𝐾𝑜𝑝𝑡 = 2 𝜌𝐻𝑅4 𝜆3𝑜𝑝𝑡
Cette valeur du couple optimal 𝑇𝑡−𝑜𝑝𝑡 est donc définie comme la consigne du couple
électromagnétique de la génératrice synchrone dans le schéma de CCO i.e : 𝑇𝑒∗ = 𝐾𝑜𝑝𝑡 𝜔𝑡 2
𝑑𝜔𝑡
(voir figure (V.26)). En régime permanent (où 𝑑𝑡
≈ 0), sous une vitesse de vent constante,
∗
la technique de CCO assure que 𝑇𝑡 = 𝑇𝑒 = 𝑇𝑒 ce
qui signifie que le système a atteint le
PPM.
Alors, nous avons choisi d’appliquer cette méthode de MPPT en raison de sa simplicité
et la facilité de son implémentation. Des blocs de multiplication de XSG suffisent pour
implémenter ce contrôleur.
V.3.3.3. Schéma de contrôle à flux orienté
La technique de commande à flux orienté (FOC), connue ainsi sous le nom de
commande vectorielle, a été proposée en 1972 par F. Blaschke pour la commande des
machines électriques à CA [197]. Par rapport au schéma de DTC, elle permet d’obtenir des
meilleures performances dynamiques en accomplissant un contrôle découplé du flux et
du couple électromagnétique d’une machine à CA [69]. Dans cette étude, on choisit
d’appliquer la technique de commande par orientation indirecte du flux à la GSAP de
l’éolienne [67-69 ,50].
Comme le montre la figure (V.27), le schéma de FOC est établi dans le repère tournant
(dq) tel que, les courants des phases statoriques (abc) sont mesurés et transformés à des
composantes d et q dans le repère de Park (dq). Le contrôle séparé du flux et du couple
peut être achevé en agissant sur les composantes d et q du courant statorique. Le flux est
∗
orienté en maintenant la composante d du courant statorique constamment nulle ( 𝑖𝑠𝑑 =
0 ). Alors que, le couple électromagnétique est contrôlé en altérant indépendamment la
composante en quadrature du courant. Tel que, la consigne du courant q est calculée en
fonction de la référence du couple comme suit :
∗ 2 𝑇∗
𝑖𝑠𝑞 = 3𝑝𝜙𝑒 (V.31)
𝑟
Comme l’indique la figure (V.27), des contrôleurs PI s’occupent de la régulation des
erreurs des courants statoriques (𝑖𝑠𝑑 et 𝑖𝑠𝑞 ). On note que les paramètres de ces régulateurs
sont déterminés par la méthode de compensation des pôles. En tenant compte les termes
′ ′
de découplage, les tensions statoriques de référence (𝑣𝑠𝑑 et 𝑣𝑠𝑞 ) sont calculées comme suit :
′
𝑣 = 𝑣𝑠𝑑 + 𝐿𝑠 𝜔𝑒 𝑖𝑠𝑞
{ ′ 𝑠𝑑 (V.32)
𝑣𝑠𝑞 = 𝑣𝑠𝑞 − 𝐿𝑠 𝜔𝑒 𝑖𝑠𝑑 − 𝜔𝑟 ϕ𝑟
′ ′
Enfin, la transformée inverse de Park est appliquée aux composantes 𝑣𝑠𝑑 et 𝑣𝑠𝑞 pour
générer les signaux de commutation du redresseur à l’aide du modulateur SPWM.
143
Chapitre V : Commande à FV orienté des systèmes PV hybrides
UDC
isd
isabc abc
𝑳𝒔 𝝎𝒓 𝒊𝒔𝒒 r
dq isq
isd
isd* + - vsd + + 𝒗′𝒔𝒅
VW r PI dq
SPWM
Te* 𝟐 isq* + abc
CCO PI
𝟑𝒑𝝓𝒓 - vsq + - 𝒗′
ωt 𝒔𝒒
isabc
d/dt r isq 𝑳𝒔 𝝎𝒓 𝒊𝒔𝒅 + 𝝎𝒓 𝛟𝒓
r
La figure (V.28.a) illustre le schéma de FOC développé dans Simulink à l’aide de XSG.
Comme le montre la figure (V.28.b), chacun des contrôleurs PI des courants (𝑖𝑠𝑑 et 𝑖𝑠𝑞 )
est établi en utilisant une formule récurrente semblable à celle de l’équation (V.27).
La figure (V.28.c) représente le bloc XSG du modulateur SPWM qui génère les signaux
de commande des interrupteurs du redresseur. Ce modulateur réalise une
comparaison entre un signal triangulaire symétrique HF (porteuse) et les tensions de
′ ′ ′
référence du redresseur (𝑣𝑠𝑎 , 𝑣𝑠𝑏 , 𝑣𝑠𝑐 ). Pour la génération de la porteuse à base de XSG,
nous avons adopté la même méthode présentée dans [185] :
Comme l’illustre la figure (V.28.c), on utilise cinq blocs XSG pour générer une
porteuse triangulaire de fréquence ftri =5kHz. Le décompteur (de 7bits) répète une
séquence de décomptage, à chaque période Ttri=(1/ftri)=0.2ms, ce qui permet de générer
un signal triangulaire asymétrique dont l’amplitude varie entre 127 et 0. On note ici
que la limite de décomptage est déterminée suivant la relation (V.28). En outre, le bloc
compteur/décompteur (de 9bits) génère un signal triangulaire symétrique
d’amplitude variable entre -32 et 32. Le bloc «Reinterpret» permet d’obtenir finalement
un signal triangulaire de fréquence 5kHz variant dans l’intervalle [-1,1].
1 In a z-3
cast a
ba b x (-0.001953) In1
wt Gateway In2 Out1 a
x 0.007813 z
-3 a-b
Mult Convert1 In2 ab cast b
b
Kopt/(1,[Link]) PI controller CMult1 AddSub1
Convert3
Mult1 a
3 In a-b Out In2
MPPT par CCO iq Gateway In1 b
x 1.328 Gateway Out2
In1 Out1 1
AddSub2
5 In In1 CMult2 PWM
Out1 cast 2 In3
idref Gateway In4
In2 teta_m
Convert4 PWM generator
4 In PI controller1
Sy stem id Gateway In3 a
Generator a -3 a+b Out
z
ab cast b
x 0.007813 b Gateway Out1
Convert2 AddSub
Mult2
CMult3
Contrôleur PI
(a)
144
Chapitre V : Commande à FV orienté des systèmes PV hybrides
Delay1
-1 a
z
a>b sel
b
Sy stem x 230 a cast a d0 1
Generator Relational
a+b a a+b Out1
Convert2
Kp b a-b cast b d1
1000000
1 a sel
b a
AddSub4 Convert1 AddSub2 Uup1 Mux
In1 a-b x 0.01367
2 b a>b
-1 AddSub3 -1000000 b d0
z
In2
AddSub1 Ki*Ts Uup Relational1
Delay
d1
Mux1
(b)
1 a
a > b sel
Va 1 b
d0
U2 Rel ati onal 3
Sy s t em
a d1 Generat or Gateway Out2
sel
a > b a
b d0 M ux Out
-1 a b
b
Rel ati onal 4
d1 Rel ati onal not Out
U1
M ux1 a Inverter2
Gateway Out1
a a b Out
2 b 1
a > b sel Gateway Out3
Vb 1 b Rel ati onal 1 not Out Sabc
d0 Inverter Gateway Out5
Rel ati onal 5 a
U
d1 a b
b Out
a
a > b sel Rel ati onal 2 Gateway Out6
b M ux3 not Out
-1 d0
Rel ati onal 6 Inverter1Gateway Out4
-2 d1
a M ux4
3
a > b sel -- [a:b] not
up rei nterpret
Vc 1 b
Sl i ce1 Inverter3 Rei nterpret
d0 Counter1 Counter2
1 Rel ati onal 7
a sel d1
b
a > b d0
M ux5
Bloc de génération du signal
-1
Rel ati onal 8 d1 triangulaire symétrique
-1
M ux6
(c)
Figure V.28: Blocs XSG du a) Schéma de FOC de la GSAP, b) Contrôleur PI du courant isd c)
Générateur des signaux MLI.
145
Chapitre V : Commande à FV orienté des systèmes PV hybrides
146
Chapitre V : Commande à FV orienté des systèmes PV hybrides
alors l’équation (III.34) en fonction des variables de contrôle (𝑖𝑔𝑞 = 𝑦1) et (𝑖𝑔𝑑 = 𝑦2) comme
suit :
1
𝑦̇ 1 = (−𝑅𝑦1 + 𝐿𝜔𝑦2 − 𝑣𝑔𝑞 + 𝑣𝑞 )
𝐿
{ 1 (V.33)
𝑦̇ 2 = 𝐿 (−𝑅𝑦2 − 𝐿𝜔𝑦1 − 𝑣𝑔𝑑 + 𝑣𝑑 )
La loi de commande de type backstepping est donc conçue étape par étape en utilisant
des fonctions candidates de Lyapunov qui assurent que les objectifs de contrôle
mentionnés ci-dessus soit atteints :
Etape 1 :
On définit l’erreur de régulation de la composante q du courant du réseau comme suit:
∗
𝑒𝑔𝑞 = 𝑖𝑔𝑞 − 𝑖𝑔𝑞 = 𝑦1∗ − 𝑦1 (V.34)
Sa dérivée temporelle est calculée, en tenant compte l’équation (III.34), comme suit :
𝑣𝑞 𝑅 𝑣𝑔𝑞
𝑒̇𝑔𝑞 = 𝑦̇ 1∗ − 𝑦̇ 1 = −𝑦̇ 1 = − 𝐿
+ 𝐿 𝑦1 − 𝜔𝑦2 + 𝐿
(V.35)
On définit la fonction de Lyapunov suivante :
2 1
𝑉1 = 𝑒𝑔𝑞 (V.36)
2
Sa dérivée est donc calculée comme suit :
𝑣 𝑅 𝑣𝑔𝑞
𝑉̇1 = 𝑒𝑔𝑞 . 𝑒̇𝑔𝑞 = 𝑒𝑔𝑞 (− 𝑞 + 𝑦1 − 𝜔𝑦2 + )
𝐿 𝐿 𝐿
(V.37)
Selon la condition de stabilité de Lyapunov [205], On doit avoir :
𝑣𝑞 𝑅 𝑣𝑔𝑞
− + 𝑦1 − 𝜔𝑦2 + = −𝑔1 . 𝑒𝑔𝑞 avec 𝑔1 > 0 (V.38)
𝐿 𝐿 𝐿
Pour que : 𝑉1̇ = −𝑔1 𝑒𝑔𝑞
2
soit négative. Ceci assure l’élimination de l’erreur du courant
en quadrature. Selon la procédure de backstepping [205], on définit vq comme entrée de
commande virtuelle qui est exprimée comme suit :
𝑣𝑞 = 𝑅𝑦1 − 𝜔𝐿𝑦2 + 𝑣𝑔𝑞 + 𝑔1 𝐿 𝑒𝑔𝑞 (V.39)
Etape 2 :
L’erreur du courant du réseau selon l’axe d est définie comme suit :
∗
𝑒𝑔𝑑 = 𝑖𝑔𝑑 − 𝑖𝑔𝑑 = 𝑦2∗ − 𝑦2 (V.40)
Sa dérivée est donc donnée comme suit :
𝑣𝑑 𝑅 𝑣𝑔𝑑
𝑒̇𝑔𝑑 = 𝑦̇ 2∗ − 𝑦̇ 2 = −𝑦̇ 2 = − 𝐿
+ 𝐿 𝑦2 + 𝜔𝑦1 + 𝐿
(V.41)
Pour ce sous système, on définit une deuxième fonction de Lyapunov comme :
2 1
𝑉2 = 𝑉1 + 2 𝑒𝑔𝑑 (V.42)
Sa dérivée est calculée comme suit :
𝑣 𝑅 𝑣𝑔𝑑
𝑉̇2 = −𝑔1 𝑒𝑔𝑞
2 2
+ 𝑒𝑔𝑑 . 𝑒̇𝑔𝑑 = −𝑔1 𝑒𝑔𝑞 + 𝑒𝑔𝑑 (− 𝑑 + 𝑦2 + 𝜔𝑦1 + )
𝐿 𝐿 𝐿
(V.43)
On doit avoir :
𝑣𝑑 𝑅 𝑣𝑔𝑑
−𝑔2 𝑒𝑔𝑑 = − 𝐿
+ 𝐿 𝑦2 + 𝜔𝑦1 + 𝐿
avec 𝑔2 > 0 (V.44)
Pour que la dérivée 𝑉̇2 soit négative : 𝑉̇2 = 2
−𝑔1 𝑒𝑔𝑞 2
− 𝑔2 𝑒𝑔𝑑 < 0.
La loi de commande 𝑣𝑑 est donc définie comme suit :
𝑣𝑑 = 𝑅𝑦2 + 𝐿𝜔𝑦1 + 𝑣𝑔𝑑 + 𝑔2 𝐿 𝑒𝑔𝑑 (V.45)
147
Chapitre V : Commande à FV orienté des systèmes PV hybrides
Cette loi de commande permet d’avoir 𝑉2 > 0 et 𝑉̇2 < 0 ce qui assure la convergence
asymptotique des erreurs des courants 𝑒𝑔𝑑 et 𝑒𝑔𝑞 vers le zéro.
Les équations (V.39) et (V.45) sont alors utilisées pour calculer les tensions de référence
(𝑣𝑑∗et 𝑣𝑞∗ ) fournies à l’entrée du modulateur SPWM (voir figure(V.29)).
1 a -3
z
x2 Ial pha a b a
1 b
a -b 1
a
ia 3 Mul t2 b id
CMult3 a-b x 0.3359 1 si n 4 a -3
a z AddSub3
2 b ialpha cos a b
a+b b
ib CMult1 2
b AddSub1 Mul t1
3 Ibeta a -3
z
ic AddSub a b a
b 2
a + b
a Mul t4 b iq
2 a
a-b x 0.5781 z-3 AddSub1
a b
b ibeta b
(a)
1 a 1 1
z-3
Id a b ialpha ia
b a
3 a + b 1
Mult2
2 x 1.732 a
sin b ialpha
a -3
ibeta a-b x 0.5 2
z a b AddSub1
b ib
b sqrt(3)
4
cos Mult1 AddSub1 CMult
a -3
z a b a
2 a
b a -b 2
Iq a -3 b ibeta a+b x (-0.5) 3
z a b Mult4
AddSub3
b ic
b
(b)
148
Chapitre V : Commande à FV orienté des systèmes PV hybrides
2 Out 3
4 a fluxi-alpha
z
-3 -0.572265625 a Gateway Out5 teta
1 x2 -1 cast a
vdc ab x 0.332 z x 2e-006 a
b a+b
Sa Constant b
a+b a -b
CMult a Convert2 cast b x_in x_out
Delay1 Delay2 b
a -b a Mult2 CMult3 Ts AddSub5
2 b Convert3 AddSub2 phase_out phase_in 4
a -b AddSub6
Sb b z-1 5 x 0.008057 Cos-Sin
AddSub y_in y_out
ilpha
3 AddSub3 a
z-3 CMult2
-1 CORDIC 4.0 CORDIC 4.0 1
Sc a ab x 0.5771 z x 2e-006 cast a a
a -b b
b a+b
Convert cast b a -b
Delay3 b
Mult1 CMult1 Ts1 1
AddSub1
Convert1 AddSub4
AddSub7 fluxi-beta
-1 CMult4
z System
6 x 0.008057 Generator
Delay4
ibeta
(c)
z-1 cast
Delay2 Convert3
a
a -3
1 a z a -b
x 0.5 b
a b
Vdcref a -b b
gain1 b AddSub6
x 0.5 a -3 2 Mult2
z x 0.5
a b AddSub4
gain2 b Vdc
gain
Mult3
w1 a
a x 0.4
-1 a
z a + b
a -b a + b b learning
b b
Delay w2 a -3
z AddSub3 AddSub5
AddSub2 a b
b 1 -1 cast
a z
Iqref
a - bx (-0.002502) z-1 Mult1 Convert1
Delay1
b a
Kp Delay3 a -3
z a -b
AddSub7 a b b
b
AddSub1
Mult4
(d)
ia
x 0.1875 a In ialpha
a+b a Gateway In1 ib
1 In a b a 3 In
Rg1 a+b
Id_ref Gateway In9 a -b b Iabc_inv Gateway In5 ic
ibeta
b AddSub4 a+b
b PMW
In
AddSub5 q VA Va
AddSub3 x 2.514 6 In Sabc 1 Gateway In14 abc/alpha-beta
[Id] In AddSub7
Vq Gateway In2 x8 d VB Vb
Gateway In8 Sa
Lg In flux-alpha 3
Sabcvf
2 In a In3 VC Vc Gateway In11 Sb
a g2*Lg flux-beta
Iq_ref Gateway In10 a -b x 0.1875 In
b a dq/abc control PWM
a -b Gateway In12 Sc
[Iq] b a In
In AddSub1 Rg a+b flux-beta 2
b Gateway In13 VDC
Gateway In4 AddSub2 a+b
x 2.514 b flux-alpha
AddSub6 5 In
4 In Ig alpha
AddSub8 Iabc Id [Id] Vdc Gateway In7
Lg1 Vd Gateway In3
x8 Sin_Cos
Ig beta
Backstepping
System Sin-cos Iq [Iq]
Generator g1*Lg estimation du FV
abc/dq
(e)
Figure V.30: Blocs XSG de : a) Transformation des coordonnées (abc/αβ)et (αβ/dq), b)
Transformation des coordonnées (dq/αβ) et (αβ/abc), c) Estimateur du FV du réseau, d) Contrôleur
ADALINE ,e) Schéma global de VFOC par Backstepping.
149
Chapitre V : Commande à FV orienté des systèmes PV hybrides
(V.32), (V.33) et (V.34) illustrent les résultats obtenus. On constate que les blocs de
contrôle de XSG ont réussi à commander correctement les convertisseurs de
puissance du SEH. Les contrôleurs P&O et CCO déterminent correctement les
points de fonctionnement optimal des SERs. La figure (V.31) montre qu’il existe une
grande ressemblance entre les résultats obtenus avec l’algorithme P&O établi à
l’aide des blocs de Simulink et celui établi à base de XSG. Néanmoins, la figure
(V.32.c) montre que les amplitudes des oscillations statiques du couple
électromagnétique en utilisant les modules de XSG sont plus importantes que celles
obtenues avec les blocs de Simulink. Ceci est dû à l’erreur de quantification
introduite par la représentation numérique à point fixe. D’autre part, la figure
(V.33.c) montre que le régulateur ADALINE adopté a parfaitement maintenu la
tension du bus continu constante à la valeur désirée (400V). La puissance active
générée par le SEH est totalement injectée au réseau électrique comme l’indique la
figure (V.34.a). Tandis que, la puissance réactive est gardée quasiment nulle ce qui
permet d’obtenir un facteur de puissance unitaire. Le taux de distorsion harmonique
du courant d’onduleur est d’environ 2.71%, ce qui répond bien aux exigences de
connexion au réseau (𝑇𝐻𝐷𝐼 < 5%). Ce résultat affirme que la commande
Backstepping a atteint les objectifs visés. La figure (V.34.c) confirme ainsi que
l’estimation du FV du réseau est réalisée correctement.
Current(A)
3.8
3.7
3.6
150 205
PV voltage (V)
100 200
Vpv-Simulink
50
Vpv-XSG 195
Vpv-Simulink
0 190
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 Vpv-XSG
t(s) 185
0 0.05 0.1 0.15 0.2 0.25
(b) t(s)
Voltage reference (V)
200 205
Voltage reference (V)
150
200
100 Vpvref-Simulink
Vpvref-XSG 195 Vpv-Simulink
50
Vpv-XSG
0 190 0.3
0 0.05 0.1 0.15 0.2 0.25 0.04 0.06 0.35
0.08 0.4
0.1 0.12 0.14 0.16 0.18
t(s) t(s)
(c)
Figure V.31: Performances de poursuite du PPM du GPV dans le premier cas de simulation.
150
Chapitre V : Commande à FV orienté des systèmes PV hybrides
805 30
Pt(W) 804.5
804 Te-XSG
25
Te-Simulink
803.5 Tt
Te et Tt (Nm)
803
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 20
t(s)
(a)
13.6
15
13.5
wt(rad/s)
13.4
10
13.3
13.2
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45
t(s) t(s)
(b) (c)
500
q-axis stator current (A)
0
Isq-XSG 450
Isq-Simulink
-5 400
-10 350
DClink voltage (V)
300
-15
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45
t(s) 250
(a)
d-axis stator current (A)
10 200
0 150
-30 0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0 0.1 0.2 0.3 0.4 0.5 0.6
t(s) t(s)
(b) (c)
4000
3000 0.4
2000
1000
0.2
0
beta-axis VF
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
t(s) 0
(a)
Reactive power (VAr)
0 -0.2
-5000 -0.4
-10000
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 -0.6 -0.4 -0.2 0 0.2 0.4 0.6
t(s) alpha-axis VF
(b) (c)
151
Chapitre V : Commande à FV orienté des systèmes PV hybrides
Dans le deuxième test, le SEH subit une variation graduelle du niveau d’irradiation
de 0.6kW/m2 vers 0.9kW/m2 puis, un changement soudain à l’instant t=0.7s au
niveau de 0.5kW/m2. La vitesse de vent est maintenue constante à VW=10m/s. Les
figures (V.35) et (V.36) montrent les résultats de simulation obtenus à l’aide de XSG.
On constate que le système de commande adopté a assuré des performances
satisfaisantes pour un niveau d’irradiation variable. Le contrôleur P&O converge
rapidement au PPM correspondant à chaque niveau d’irradiation comme l’illustre
la figure (V.35). Les erreurs statiques des courants direct et en quadrature sont
éliminées grâce à l’action de la commande Backstepping (voir figure (V.36.a) et
(V.36.b)). L’estimateur du FV a garanti la synchronisation du courant d’onduleur
(phase a) avec la tension du réseau comme le montre la figure (V.36.c). En outre,
l’ADALINE numérique a stabilisé la tension du bus CC à la valeur de référence avec
un dépassement atténué (d’environ 479.2V) (voir figure (V.36.d)).
210
Vref
PV voltage (V)
205 Vpv
200
195
190
0 0.2 0.4 0.6 0.8 1
t(s)
(a)
PV current (A)
2.5
2
60 1
0
40
-1
20 -2
2
(a) 1
80 0
q-axis current (A)
-1
60
-2
0
0 0.2 0.4 0.6 0.8 1
t(s)
(b)
152
Chapitre V : Commande à FV orienté des systèmes PV hybrides
-100
-200
0.45 0.5 0.55 0.6
t(s)
(c)
500
DC-link voltage (V)
400
300
200
100
0
0 0.2 0.4 0.6 0.8 1
t(s)
(d)
V.3.5. Conclusion 2
Dans cette partie du chapitre, l’objectif a été d’examiner la procédure de conception
d’un circuit de commande d’un SEH, à l’aide de XSG, dans le but de l’implémenter sur
une carte FPGA. Les contrôleurs locaux des SERs ont été développés en adoptant des
algorithmes classiques de MPPT. En outre, le contrôle de l’onduleur a été assuré via un
schéma de VFOC établi à base de la commande Backstepping. Le choix porté à ces
techniques de contrôle se justifie par la simplicité de leurs algorithmes ce qui assure une
faible consommation des ressources du circuit FPGA. La simulation XSG/Simulink a été
accomplie afin de valider le schéma de commande construit à base de l’outil XSG. Les
résultats obtenus ont démontré la précision et la faisabilité du circuit numérique
développé. Cette étape nous permettra alors de passer directement à l’étape de réalisation
pratique et de validation expérimentale du système de commande étudié.
153
Conclusion générale
Les travaux scientifiques qui ont été présentés dans cette thèse sont le fruit d’un long travail
de recherche réalisé au niveau du laboratoire des énergies renouvelables de l’université de
Jijel. Ils peuvent être positionnés à l’intersection de deux principaux axes : la commande des
systèmes PV hybrides et l’intelligence artificielle. En effet, la stratégie de gestion/commande
est l’élément essentiel d’un système énergétique multi-sources sans lequel un fonctionnement
correct à haut rendement ne peut avoir lieu. La présence de sources énergétiques de nature
différente au sein du même système exige leur pilotage d’une manière adéquate en assurant
une bonne coordination entre elles. En outre, les systèmes de stockage d’énergie à base de
batteries doivent être soigneusement contrôlés afin de respecter certaines contraintes de
fonctionnement en vue de les protéger et de prolonger leur durée de vie. Alors, le rôle attribué
à un tel système de commande consiste à assurer la supervision du flux énergétique et à traiter
de manière optimale les non-linéarités du SEH pour différentes conditions de fonctionnement.
D’ailleurs, les algorithmes de commande doivent se baser sur des mécanismes susceptibles
d’améliorer leur efficacité.
Dans nos travaux, différentes structures de systèmes énergétiques interconnectés au réseau
électrique ont été étudiées. L’accent a été mis sur les TIAs qui présentent un moyen pertinent
pour l’optimisation des schémas de commande classiques en vue de leur permettre d’agir
intelligemment sur les SEHs dynamiques. Plusieurs aspects de la problématique de
commande des systèmes de production d’énergie électrique ont été traités :
La poursuite du PPM : Un régulateur MPPT d’une SER (PV ou éolienne) a pour rôle de
suivre le point d’opération optimal qui change sa position conformément aux conditions
climatiques. Dans cette thèse, la tâche de contrôle MPPT, pour chacune des sources PV et
éolienne, a été affectée à des RNAs adaptatifs qui utilisent des algorithmes d’apprentissage en
temps réel. À l’inverse du RNA entrainé hors-ligne, le RNA qui apprend en-ligne ne requiert
pas une connaissance à priori du système commandé. En outre, il n’exige qu’un minimum
d’information provenant du système commandé contrairement au CLF. Les résultats de
simulation obtenus démontrent la supériorité des suiveurs neuronaux, adoptés pour les GPV
et GE, par rapport aux régulateurs MPPT classiques. Les oscillations statiques sont
Conclusion générale
considérablement réduites et les vitesses de convergence vers les PPMs sont améliorées d’une
manière satisfaisante.
Le réglage des paramètres des régulateurs PI classiques : L’algorithme de PSO a été
appliqué dans le chapitre III pour l’ajustement hors-ligne des gains proportionnel et intégral
des contrôleurs PI utilisés au sein des schémas de commande conventionnels. Dans un SPVCR,
la synthèse du contrôleur PI de la tension du bus CC est accomplie à l’aide de la technique
PSO. En outre, les régulateurs PI inclus dans les schémas de contrôle DTC et VOC, qui sont
adoptés pour un SGEE, sont ajustés hors-ligne via l’algorithme PSO. Cette technique de
réglage de paramètres est simple et ne requiert pas un modèle mathématique du système
contrôlé. Elle permet de trouver les gains PI optimaux en vérifiant hors-ligne un ensemble de
solutions éventuelles. Des bonnes réponses en régimes transitoire et statique ont été obtenues
en exprimant le coût de l’algorithme de PSO en fonction des performances de contrôle qu’on
désire les optimiser.
D’autre part, dans le chapitre V, des RGNFs ont été introduits afin d’assurer l’adaptation
des gains des contrôleurs PI des composantes d et q du courant d’onduleur dans un SEH.
L’apprentissage de l’ANFIS est accompli hors-ligne à l’aide d’une base de données, ce qui lui
permet d’acquérir des connaissances à priori sur le système commandé. En outre, un
régulateur flou a été utilisé pour décider, en se basant sur les entrées d’erreur et de variation
d’erreur, les valeurs des gains des contrôleurs de la PAC-OS et du SSEB. Les superviseurs flous
adoptés performent alors une auto-adaptation en temps réel des paramètres PI permettant de
s’affranchir de la linéarité des correcteurs PI à gains fixes.
La substitution des contrôleurs classiques : En vue d’améliorer les performances des
schémas de commande traditionnels, des régulateurs PI ont été remplacés par des contrôleurs
intelligents basés sur la LF et les RNAs. Un CFMG a été utilisé pour la régulation de la tension
bus CC du système PV hybride adopté dans le chapitre V. Dans le chapitre IV, un schéma de
commande basé sur des RNAs adaptatifs a été élaboré afin d’assurer la supervision d’un MR
hybride. Où, des contrôleurs neuronaux ajustés en-ligne, à l’aide de l’algorithme d’AI, sont
proposés pour optimiser les performances du schéma de DPC de l’onduleur. En outre, des
régulateurs locaux de type ADALINE ont été adoptés pour assurer le contrôle du GE et des
sources d’appoint du MR. Un superviseur central du flux d’énergie a été établi à base de la LF
afin d’assurer l’alimentation continue des charges de consommation. Les simulations
accomplies prouvent la supériorité et l’efficacité des contrôleurs intelligents adoptés en
comparaison avec les contrôleurs classiques.
Conception d’un circuit de commande numérique pour une future implémentation sur
FPGA: Dans la deuxième partie du dernier chapitre, un schéma de commande simplifié d’un
système PV-éolien a été établi à l’aide de l’outil de développement XSG. Les étapes de
conception ont été présentées et les résultats de simulation ont été fournis afin de valider le
circuit numérique ainsi développé.
De nombreuses perspectives de recherche sont envisagées afin d’améliorer et compléter les
travaux menés dans cette thèse. Nous suggérons particulièrement de :
156
Conclusion générale
157
Annexe A
Les caractéristiques électriques du module BPMSX 120, fournies par le datasheet, sont
récapitulées dans tableau suivant :
Caractéristiques électriques BPMSX 120
Puissance maximale (Pmax) 120 W
Nombre de cellules en série (NS) 72
Tension au point de puissance maximale (VPPM) 33.7 V
Courant au point de puissance maximale (IPPM) 3.56 A
Courant de court-circuit (ICC) 3.87 A
Tension en circuit ouvert (VCO) 42.1 V
Coefficient de température de la puissance -(0. 5±0.05) %/°C
Coefficient de température du courant de court-circuit (0.065±0.015)%/°C
Coefficient de température de la tension de circuit ouvert -(80±10) mV/°C
NOCT 47±2 °C
𝐶𝑃 (𝜆) = ∑ 𝐶𝑘 𝜆𝑘
𝑘=0
Tel que: 𝐶0 = 0.110898,𝐶1 = −0.02493,𝐶2 = 0.057456, 𝐶3 = −0.01098, 𝐶4 = 0.00054.
158
J Coefficient d’inertie totale 5kg.m2
F Coefficient de frottement total 9.08 10-3Nm.s
Pnom Puissance nominale 900W
Paramètres de la turbine Darrieus
R Rayon de la turbine 1m
H Hauteur 2m
Q2
D
V1 C1 C2 V2
159
Paramètres du MR étudié dans le chapitre 5
Q Capacité nominale 50 Ah
E0 Tension nominale 288 V
Rb Résistance interne 0.0576 Ω
VDC* Consigne de tension du bus CC 900V
r Résistance de la PAC-OS de 50kW 3.2813e-4 Ω
VRMS RMS de la tension du réseau 380V
CDC Capacité du bus CC 1200µF
Annexe B
Calcul des dérivées partielles d’équation (III.29) :
(2) (2)
𝜕𝜎𝜔 (𝑘) 𝜕𝜎𝜔 (𝑘) 𝜕𝑆𝜔 𝜕𝑦 (3) 𝜕𝑦𝑗 𝜕𝑁𝑗
𝑛
𝜕𝑤𝑖𝑗
= 𝜕𝑆𝜔 𝜕𝑦 (3) 𝜕𝑦 (2) 𝜕𝑁 (2) 𝜕𝑤𝑖𝑗
𝑛
𝑗 𝑗
𝜕𝜎𝜔 (𝑘) 𝜕𝜎𝜔 (𝑘) 𝜕𝑆𝜔 𝜕𝑦 (3) j=1,2,3, i=1,2
𝜕𝑤0𝑛
= 𝜕𝑆𝜔 𝜕𝑦 (3) 𝜕𝑤0𝑛
𝜕𝜎𝜔 (𝑘) 𝜕𝜎𝜔 (𝑘) 𝜕𝑆𝜔 𝜕𝑦 (3)
𝜕𝑤𝑗𝑛
= 𝜕𝑦 (3) 𝜕𝑤𝑗𝑛
{ 𝜕𝑆𝜔
(2) (2)
−𝑁 2
(2) (2) (2) 1−𝑒 𝑗 𝜕𝑦𝑗 1 (2)
Sachant que : 𝑦𝑗 (𝑘) = 𝑔𝑗 (𝑁𝑗 ) = (2) ⟹ (2) = 2 [1 − (𝑦𝑗 ) ],
−𝑁 𝜕𝑁𝑗
1+𝑒 𝑗
(2)
𝜕𝜎𝜔 𝜕𝑦 (3) 𝜕𝑁𝑗 (1)
𝜕𝑆𝜔
= 𝑆𝜔 , (2) = 𝑤𝑗𝑛 et 𝑛
𝜕𝑤𝑖𝑗
= 𝑦𝑖 .
𝜕𝑦𝑗
160
Calcul des dérivées partielles d’équation (V.6) :
(4) (4) (2) (2)
𝜕𝜎𝑒 𝜕𝜎𝑒 𝜕𝑒𝑆 𝜕𝑧𝑜 𝜕𝑛 𝜕𝑧𝑗 𝜕𝑛𝑗
𝑐 = . . . . (2) .
(4)
𝜕𝑤𝑖𝑗 𝜕𝑒𝑆 𝜕𝑧𝑜 𝜕𝑛 (4)
𝜕𝑧
(2)
𝜕𝑛 𝜕𝑤𝑖𝑗𝑐
𝑗 𝑗
(4) (4)
𝜕𝜎𝑒 (𝑘) 𝜕𝜎𝑒 𝜕𝑒𝑆 𝜕𝑧𝑜 𝜕𝑛
= . (4) . (4) .
𝜕𝑤0 𝜕𝑒𝑆 𝜕𝑧𝑜 𝜕𝑛 𝜕𝑤0
(4) (4)
𝜕𝜎𝑒 (𝑘) 𝜕𝜎𝑒 𝜕𝑒𝑆 𝜕𝑧𝑜 𝜕𝑛
= . (4) . (4) .
{ 𝜕𝑤𝑗 𝜕𝑒𝑆 𝜕𝑧𝑜 𝜕𝑛 𝜕𝑤𝑗
(2) (2)
−𝑛 2
(2) (2) (2) 1−𝑒 𝑗 𝜕𝑧𝑗 1 (2)
Sachant que : 𝑧𝑗 (𝑘) = 𝑆𝑗 (𝑛𝑗 ) = (2) ⟹ (2) = 2 [1 − (𝑧𝑗 ) ] ,
−𝑛 𝜕𝑛𝑗
1+𝑒 𝑗
(4) (4) (2)
𝜕𝑧𝑜 (4) (4) 𝜕𝜎 𝜕𝑛 𝜕𝑛𝑗 ( )
(4) = 𝛽𝑧𝑜 [1 − 𝑧𝑜 ] , 𝑒 = 𝑒𝑆 (𝑘), (2) = 𝑤𝑗 , 𝑐 = 𝑧𝑖 1
𝜕𝑛 𝜕𝑒𝑆 𝜕𝑧𝑗 𝜕𝑤𝑖𝑗
Annexe C
Architecture d’ANFIS : On considère une base de règles floues constituée de deux règles Si-
Alors qui sont exprimées comme suit [89] :
Règle 1: Si y1 est A1 et y2 est B1, Alors f1=p1 y1+q1 y2 +r1.
Règle 2: Si y1 est A2 et y2 est B2, Alors f2=p2 y1+q2 y2 +r2.
Avec (A1, A2) et (B1, B2) sont, respectivement, les labels linguistiques associés à chacune des
entrées y1 et y2. {𝑝𝑗 , 𝑞𝑗 , 𝑟𝑗 } (j=1,2) est l’ensemble des coefficients linéaires du polynôme de
conclusion de la règle j. Les nœuds de la même couche sont définis par des fonctions
semblables. Les informations suivantes sont compilées des références [89,95] :
Couche 1 : La sortie de chaque nœud i dans la couche d’entrée, notée 𝑂𝑖1 est donnée par :
𝑂𝑖1 = 𝜇𝐴𝑖 (𝑦1 ) 𝑝𝑜𝑢𝑟 𝑖 = 1,2.
Ou 𝑂𝑖1 = 𝜇𝐵𝑖−2 (𝑦2 ) 𝑝𝑜𝑢𝑟 𝑖 = 3,4.
Avec 𝑂𝑖1
représente le degré d’appartenance d’une entrée (y1 ou y2) à un ensemble flou (A1 ,
A2 ou B1, B2). La fonction d’appartenance 𝜇𝐴𝑖 (𝑦𝑖 ) est souvent choisie d’être une fonction en
forme de cloche (Generalized–bell function) qui est définie comme suit :
1
𝜇𝐴𝑖 (𝑦𝑖 ) =
𝑦𝑖 − 𝑐𝑖 2𝑏𝑖
1+( 𝑎 )
𝑖
{𝑎𝑖 , 𝑏𝑖 , 𝑐𝑖 } est l’ensemble des paramètres non-linéaires qui altère la forme de la fonction
d’appartenance à l‘ensemble flou Ai. Les paramètres de cette couche sont dits paramètres de la
partie prémisse.
Couche 2 : Chaque nœud i de cette couche effectue le produit des signaux d’entrée comme:
𝑂𝑖2 = 𝑤𝑖 = 𝜇𝐴𝑖 (𝑦1 ) × 𝜇𝐵𝑖 (𝑦2 ) , 𝑖 = 1,2
La sortie de chaque nœud de cette couche représente le degré d’activation d’une règle.
Couche 3 : Le nœud i dans cette couche calcule le rapport entre le degré d’activation de la
règle i et la somme des degrés d’activation de toutes les règles :
𝑤𝑖
𝑂𝑖3 = 𝑤
̅𝑖 = , 𝑖 = 1,2.
𝑤1 + 𝑤2
Cette couche renvoie alors en sortie des degrés d’activation normalisés.
Couche 4 : La fonction linéaire de chaque nœud i adaptatif, dans cette couche, est définie par:
𝑂𝑖4 = 𝑤
̅ 𝑖 𝑓𝑖 = 𝑤
̅ 𝑖 (𝑝𝑖 𝑦1 + 𝑞𝑖 𝑦2 + 𝑟𝑖 )
̅ 𝑖 est la sortie de chaque nœuds i de la couche 3, et {𝑝𝑖 , 𝑞𝑖 , 𝑟𝑖 }, 𝑖 = 1,2 présente
Avec 𝑤
l’ensemble des paramètres de la partie conclusion.
161
Couche 5 : Le nœud fixe de la couche de sortie calcule la somme de tous les signaux d’entrée
comme suit :
∑𝑖 𝑤𝑖 𝑓𝑖
𝑂15 = ∑ 𝑤
̅ 𝑖 𝑓𝑖 =
∑𝑖 𝑤𝑖
𝑖
Principe du contrôle par mode glissant
La commande par mode glissant est une technique de commande non linéaire qui permet
de forcer l’état d’un système de suivre une trajectoire qui conduit, dans un temps fini, à une
surface de glissement. On considère un système non linéaire ayant la forme suivante [199]:
𝑥̇ = 𝑓(𝑥, 𝑡) + 𝑔(𝑥). 𝑢(𝑡)
x est le vecteur d'état et f (x,t) est la fonction décrivant l'évolution du système au cours du temps.
Les étapes principales pour la conception d’une loi de commande par mode glissant peuvent
être résumées comme suit [199]:
La première étape est la sélection de la surface de commutation 𝜎 (𝑥, 𝑡) = 0 qui fournie le
comportement asymptotique désiré.
Le calcul de la commande équivalente 𝑢𝑒𝑞 (𝑡) en prenant en considération les conditions
d’invariance suivantes :
𝜎 (𝑥, 𝑡) = 0 𝑒𝑡 𝜎̇ (𝑥, 𝑡) = 0 , avec 𝑢(𝑡) = 𝑢𝑒𝑞 (𝑡)
L’existence de la commande équivalente assure la faisabilité du mouvement glissant sur la
surface de commutation 𝜎 (𝑥, 𝑡) = 0 .
Finalement, nous choisissons une entrée de commande non linéaire 𝑢𝑛 (𝑡) qui doit vérifier
le critère de stabilité de Lyapunov c.à.d. 𝜎̇ 𝜎 < [Link] contrôle 𝑢𝑛 (𝑡) , dit attractif, permet de
contrôler le système en régime transitoire de telle sorte que les grandeurs contrôlées
tendent vers leurs références.
162
Références Bibliographiques
[1] N. Bizon , H. Shayeghi, N. M. Tabatabaei, “Analysis, Control and Optimal Operations in Hybrid Power
Systems Advanced Techniques and Applications for Linear and Nonlinear Systems” Green Energy and
Technology, Springer, London, 2013, ISBN 978-1-4471-5537-9.
[2] S. Upadhyay , M.P. Sharma, “A review on configurations, control and sizing methodologies of hybrid
energy systems”, Renewable and Sustainable Energy Reviews , vol. 38, pp.47–63, 2014.
[3] P. [Link]ٌno , A. J. Gil-Mena , F. [Link] , C. A. [Link], L. M. Fernandez-Ramirez , F. Jurado,
“Power control based on particle swarm optimization of grid-connected inverter for hybrid renewable
energy system”, Energy Conversion and Management, vol.91 , pp.83–92,2015.
[4] Nicola Femia, Giovanni Petrone, Giovanni Spagnuolo, and Massimo Vitelli, “Power Electronics and
Control Techniques for Maximum Energy Harvesting in Photovoltaic Systems”, Taylor & Francis
Group, LLC, 2013, ISBN 978-1-4665-0691-6.
[5] G. Dileep, [Link], “ Maximum power point tracking of solar photovoltaic system using modified
perturbation and observation method , Renewable and Sustainable Energy Reviews, vol. 50, pp.109–
129, 2015.
[6] A. [Link] , [Link] , A. Cheriti , [Link], P. Sicard , M. Benghanem, A survey on control
of electric power distributed generation systems for microgrid applications, Renew. Sustain. Energy
Review, vol. 44 , pp.751–766, 2015.
[7] M. H. Nehrir, C. Wang, K. Strunz, H. Aki, R. Ramakumar, J. Bing, Z. Miao, and Z. Salameh, “A Review
of Hybrid Renewable/Alternative Energy Systems for Electric Power Generation: Configurations,
Control, and Applications”, IEEE Transactions on sustainable energy, vol. 2, no. 4, pp. 392-403, October
2011.
[8] M. S. Carmeli , F. [Link] , M. Mauri , G. Marchegiani , D. Rosati, « Control strategies and
configurations of hybrid distributed generation systems”, Renewable Energy, vol. 41, pp. 294-305, 2012.
[9] F. Blaabjerg, R. Teodorescu, M. Liserre, A.V. Timbus, “Overview of Control and Grid Synchronization
for Distributed Power Generation Systems”, IEEE Transactions on industrial electronics, vol. 53, no. 5,
pp.1398- 1409, 2006.
[10] L. Hassaine , [Link] , [Link] , [Link], “Overview of power inverter topologies and control
structures for grid connected photovoltaic systems “, Renewable and Sustainable Energy Reviews,
vol.30, pp.796–807, 2014.
[11] H. M. Jabr, D. Lu, N. C. Kar, Design and Implementation of Neuro-Fuzzy Vector Control for Wind-
Driven Doubly-Fed Induction Generator, IEEE Trans. Sustain. Energy, vol.2, pp.404-413, 2011.
[12] T.O. Kowalska, F. Blaabjerg, J. Rodríguez, “Advanced and Intelligent Control in Power Electronics and
Drives,” Studies in Computational Intelligence, Vol.531, Springer [Link] Switzerland, 2014,
ISBN 978-3-319-03400-3.
[13] G. Bartolini ; A. Ferrara ; A.A. Stotsky, “Robustness and performance of an indirect adaptive control
scheme in presence of bounded disturbances”, IEEE Transactions on Automatic Control, Vol.44, Issue.
4, Apr 1999.
[14] A. Labouret , M. Villoz, “Solar Photovoltaic Energy”, IET renewable energy series 9, London, United
Kingdom,2010, ISBN 978-1-84919-154-8.
[15] M. R. Patel, «Wind and Solar Power Systems: Design, Analysis, and Operation» 2ème edition, Taylor
& Francis Group, 2006, ISBN 978-0-8493-1570-1.
[16] A. Labouret , P. Cumunel, J.P. Braun , B. Faraggi, « Cellules solaires : Les bases de l’énergie
photovoltaïque » 5ème édition ETSF , Dunod, Paris,2010, ISBN 978-2-10-055598-7.
[17] A. Cid Pastor, “Conception et réalisation de modules photovoltaïques électroniques», Thèse de Doctorat,
Institut des Sciences Appliquées de Toulouse, 2006.
[18] S. Chakraborty, M. G. Simُ oes, W. E. Kramer , « Power Electronics for Renewable and Distributed
Energy Systems : A Sourcebook of Topologies, Control and Integration», Green Energy and
Technology, Springer, 2013, ISBN 978-1-4471-5103-6.
[19] M. H. Quang, »Optimisation de la production de l'électricité renouvelable pour site isolé» Thèse de
Doctorat, Université de Reims Champagne-Ardenne,2013.
[20] Q.C. Zhong, T. Hornik, « Control of power inverters in renewable energy and smart grid integration»
IEEE Press, John Wiley & Sons, Ltd, 2013, ISBN: 978-0-470-66709-5.
[21] H. Abu-Rub, M. Malinowski, K. Al-Haddad, « Power Electronics For Renewable Energy Systems,
Transportation And Industrial Applications» IEEE Press and John Wiley & Sons Ltd, 2014, ISBN 978-
1-118-63403-5.
163
[22] M. Nadjah, M. Khechana, L. Laiche, T. Ouksel, C. Mahfoudi, «Etude de l’hélice d’une éolienne de 5
kW » Revue des Energies Renouvelables CISM’08 Oum El Bouaghi , pp. 257 – 264, 2008.
[23] D. Rekioua, «Wind Power Electric Systems: Modeling, Simulation and Control», Green Energy and
Technology, Springer-Verlag, London, 2014, ISBN 978-1-4471-6424-1.
[24] Siegfried Heier, «Grid integration of wind energy: onshore and offshore conversion systems», Third
Edition, John Wiley & Sons, Ltd, 2014, ISBN 978-1-119-96294-6.
[25] D. Kumar, K. Chatterjee, «A review of conventional and advanced MPPT algorithms for wind energy
systems», Renewable and Sustainable Energy Reviews vol.55, pp. 957–970,2016.
[26] M.A. Abdullah, A.H.M. Yatim, C.W. Tan, R. Saidur, «A review of maximum power point tracking
algorithms for wind energy systems», Renewable and Sustainable Energy Reviews vol.16, pp.3220–
3227, 2012.
[27] Z. Zeng, [Link], [Link], [Link], «Topologies and control strategies of multi-functional grid-
connected inverters for power quality enhancement: A comprehensive review «, Renewable and
Sustainable Energy Reviews vol. 24, pp.223–270, 2013.
[28] G. Hoogers, «fuel cell technology handbook», CRC Press LLC, 2003, ISBN 0-8493-0877-1.
[29] M. H. Nehrir, C. Wang, «Modeling and control of fuel cells: distributed generation applications», IEEE
Press, A John Wiley & Sons, Inc., Publication, 2009, ISBN 978 0 470 23328 3.
[30] Pierre Mayé, “Générateurs électrochimiques : piles, accumulateurs et piles à combustibles” Dunod, Paris,
2010, ISBN 978-2-10-055564-2.
[31] Y. Zhu, K. Tomsovic, “Development of models for analyzing the load-following performance of
microturbines and fuel cells”, Electric Power Systems Research Vol.62, pp.1-11, 2002.
[32] B. Huang, Y. Qi, A. M. Murshed,, “ dynamic modelling and predictive control in solid oxide fuel cells:
first principle and data-based approaches”, John Wiley & Sons Ltd, 2013, ISBN 978-0-470-97391-2.
[33] E. F. MOGOŞ, « Production Décentralisée dans les Réseaux De Distribution. Etude Pluridisciplinaire
De La Modélisation Pour Le Contrôle Des Sources »,Thèse de doctorat ,Ecole Supérieure d’Arts et
Métiers Centre de Lille,2005.
[34] Y. H. Li, S. Rajakaruna, S. S. Choi, « Control of a Solid Oxide Fuel Cell Power Plant in a Grid-Connected
System”, IEEE Transactions on Energy Conversion, Vol. 22, No. 2, June 2007, pp. 405-413.
[35] A. Chauhan , [Link], “A review on Integrated Renewable Energy System based power generation for
stand-alone applications : Configurations ,storage options, sizing methodologies and control”,
Renewable and Sustainable Energy Reviews, vol.38, pp 99–120, 2014.
[36] S. Sumathi, L. Ashok Kumar , P. Surekha, “Solar PV and Wind Energy Conversion Systems An
Introduction to Theory, Modeling with MATLAB/SIMULINK, and the Role of Soft Computing
Techniques”, Green Energy and Technology, Springer International Publishing Switzerland, 2015, ISBN
978-3-319-14940-0.
[37] A. Dhaker, « Contribution au dimensionnement et à l'optimisation des systèmes hybrides éoliens-
photovoltaïques avec batteries pour l'habitat résidentiel autonome », Thèse de doctorat, université de
Poitiers, 2012.
[38] D. Rekioua , E. Matagne, “ Optimization of Photovoltaic Power Systems: Modelization, Simulation and
Control”, Green Energy and Technology, Springer, 2012, ISBN 978-1-4471-2348-4.
[39] N. Karami , N. Moubayed , R. Outbib, “Energy management for a PEMFC–PV hybrid system”, Energy
Conversion and Management, vol.82 , pp.154–168, 2014.
[40] [Link], [Link], [Link] “Control strategies for a hybrid renewable energy
system: A review”, Renewable and Sustainable Energy Reviews, vol. 42, pp.597–608, 2015.
[41] L. Stoyanov, « Structures de systèmes hybrides à sources d’énergie renouvelables », Revue de Génie
Industriel, vol 5, pp.1-19, 2010, ISSN 1313-8871.
[42] L. Stoyanov « Etude de Différentes Structures de Systèmes Hybrides à Sources D’énergie Renouvelables
«, Thèse de Doctorat, Université Technique de Sofia ,2011.
[43] Y.M. Chen, C.S. Cheng, H.C. Wu, “Grid-Connected Hybrid PV/Wind Power Generation System with
Improved DC Bus Voltage Regulation Strategy”, IEEE Twenty first annual applied power electronics
conference and exposition, APEC’06, pp.1088-1094, 2006.
[44] N.A. Ahmed, A.K. Al-Othman, M.R. AlRashidi, “Development of an efficient utility interactive
combined wind/photovoltaic/fuel cell power system with MPPT and DC bus voltage regulation”, Electric
Power Systems Research, vol. 81, pp. 1096–1106, 2011.
[45] J.E. Paiva , A.S. Carvalho, “Controllable hybrid power system based on renewable energy sources for
modern electrical grids”, Renewable Energy, vol.53, pp.271-279, 2013.
[46] M. Mohammadi, M. Nafar, Fuzzy sliding-mode based control (FSMC) approach of hybrid micro-grid in
power distribution systems, Electrical Power and Energy Systems, Vol. 51 , pp. 232–242, 2013.
164
[47] [Link], M. Cebeci, “Grid connected fuel cell and PV hybrid power generating system design with
Matlab Simulink”, international journal of hydrogen energy, Vol. 39, pp.8803-8812,2014.
[48] T. Vigneysh, N. Kumarappan, Autonomous operation and control of photovoltaic/solid oxide fuel cell/
battery energy storage based microgrid using fuzzy logic controller, Hydrogen Energy,vol.41, No.3,
pp.1877-1891, 2016.
[49] P. Thounthong, A. Luksanasakul, P. Koseeyaporn, B. Davat, “Intelligent Model-Based Control of a
Standalone Photovoltaic/Fuel Cell Power Plant With Supercapacitor Energy Storage”, IEEE transactions
on sustainable energy, vol. 4, no. 1,pp.240-249 , january 2013.
[50] G. Boukettaya, L. Krichen, “A dynamic power management strategy of a grid connected hybrid
generation system using wind, photovoltaic and Flywheel Energy Storage System in residential
applications”, Energy, vol. 71, pp.148-159, 2014.
[51] N. Mezzai , D. Rekioua , T. Rekioua , A. Mohammedi , K. Idjdarane , S. Bacha, “Modeling of hybrid
photovoltaic/wind/fuel cells power system”, international journal of hydrogen energy vol.39, pp.15158-
15168,2014.
[52] V. Dash , P. Bajpai, “Power management control strategy for a stand-alone solar photovoltaic-fuel cell.
battery hybrid system”, Sustainable Energy Technologies and Assessments, vol.9, pp.68–80,2015.
[53] Tareq Alnejaili , Said Drid , Driss Mehdi , Larbi Chrifi-Alaoui , Rafik Belarbi , Aziz Hamdouni,
“Dynamic control and advanced load management of a stand-alone hybrid renewable power system for
remote housing”, Energy Conversion and Management vol.105 , pp.377–392, 2015.
[54] C.M. Hong , T.C. Ou , K.H. Lu, « Development of intelligent MPPT (maximum power point tracking)
control for a grid-connected hybrid power generation system”, Energy, vol.50, pp.270-279, 2013.
[55] F. Baghdadi , K. Mohammedi , S. Diaf , O. Behar, “Feasibility study and energy conversion analysis of
stand-alone hybrid renewable energy system”, Energy Conversion and Management, vol.105 pp. 471-
479, 2015.
[56] M. Kalantar , S.M. [Link] , “Dynamic behavior of a stand-alone hybrid power generation system of
wind turbine, microturbine, solar array and battery storage”, Applied Energy, vol. 87, pp.3051–
3064,2010.
[57] M. Parvez , M.F.M. Elias , N.A. Rahim , N. Osman, « Current control techniques for three-phase grid
interconnection of renewable power generation systems: A review”, Solar Energy , vol.135, pp.29–42,
2016.
[58] E. Félice, P. Révilla, “ Qualité des réseaux électriques et efficacité énergétique », Dunod, Paris, ISBN
978-2-10-054618-3, 2009.
[59] F. Spertino , G. Graditi, « Power conditioning units in grid-connected photovoltaic systems: A
comparison with different technologies and wide range of power ratings”, Solar Energy, vol.108, pp.
219–229, 2014.
[60] A. Hajizadeh , M. [Link], ‘Intelligent robust control of hybrid distributed generation system under
voltage sag”, Expert Systems with Applications, vol.37 pp.7627–7638, 2010.
[61] L. Belhadji , « Optimisation du contrôle commande d'un système hydraulique réversible à vitesse
variable », Engineering Sciences, Université de Grenoble, 2013. French. <tel-00881633v1>.
[62] R. Kadri, H. Andrei, J.P. Gaubert, T. Ivanovici, G. Champenois, P. Andrei, « Modeling of the
photovoltaic cell circuit parameters for optimum connection model and real-time emulator with partial
shadow conditions, Energy , Vol.42, pp. 57-67, 2012.
[63] K. Ishaque , Z. Salam, H. Taheri, Syafaruddin, “Modeling and simulation of photovoltaic (PV) system
during partial shading based on a two-diode model”, Simulation Modelling Practice and Theory, vol.19,
pp.1613–1626, 2011.
[64] G. Walker, “Evaluating MPPT converter topologies using a Matlab PV Model”, Journal of Electrical
and Electronics Engineering, Australia, January 2001.
[65] S. Roy Chowdhury, [Link], “Maximum power point tracking of partially shaded solar photovoltaic
arrays”, Solar Energy Materials & Solar Cells , Vol.94, pp.1441-1447, 20.
[66] D. Lalili, A. Mellit, N. Lourci, B. Medjahed, C. Boubakir, « State feedback control and variable step size
MPPT algorithm of three-level grid-connected photovoltaic inverter”, Solar Energy, Vol. 98, pp.561–
571, 2013.
[67] [Link]ée , F. Auger , M. Macé , L. Loron,” Design of an efficient small wind-energy conversion system
with an adaptive sensorless MPPT strategy”, Renewable Energy, vol.86 pp. 280-291,2016.
[68] J. Brahmi, L. Krichen , A. Ouali, « A comparative study between three sensorless control strategies for
PMSG in wind energy conversion system”, Applied Energy, vol.86,pp. 1565–1573, 2009.
[69] B. Wu, Y. Lang, N. Zargari, S. Kouro ,“Power conversion and control of wind energy systems,” John
Wiley & Sons Inc. and IEEE Press; 2011, ISBN: 978-1-118-02900-8.
[70] T. GHENNAM, « Supervision d’une ferme éolienne pour son intégration dans la gestion d’un réseau
électrique, Apports des convertisseurs multi niveaux au réglage des éoliennes à base de machine
asynchrone à double alimentation », Thèse de doctorat, Ecole Militaire Polytechnique d’Alger et Ecole
Centrale de Lille, 2011.
[71] [Link]é, « Aide mémoire Electrotechnique », DUNOD, Paris, 2006, ISBN 2 10 049578 X.
165
[72] J. Rodriguez , P. Cortes, “Predictive control of power converters and electrical drives”,IEEE A John
Wiley & Sons, Ltd.,2012, ISBN 978-1-119-96398-1.
[73] M. M. M. EZZAT , “Commande non lineaire sans capteur de la machine synchrone a aimants
permanents”, Thèse de Doctorat, école centrale de nantes,2011.
[74] Jean-Paul Louis,” Modélisation des machines électriques en vue de leur commande : concepts
généraux », Lavoisier, 2004, ISBN 2-7462-0916-0.
[75] J. Padulles, G.W. Ault, J.R. McDonald, “An integrated SOFC plant dynamic model for power systems
simulation,” Journal of Power Sources vol.86, pp.495–500, 2000.
[76] Y. Zhu, K. Tomsovic, “Development of models for analyzing the load-following performance of
microturbines and fuel cells”, Electric Power Systems Research, vol. 62, pp.1-11,2002.
[77] O. Tremblay, L.A. Dessaint, A.I. Dekkiche, A Generic Battery Model for the Dynamic Simulation of
Hybrid Electric Vehicles, In: IEEE Proc. Vehicle Power and Propulsion Conf.2007, pp.284-289.
[78] H. S.-Ramirez , R.S.-Ortigoza, “Control Design Techniques in Power Electronics Devices”, Springer,
2006, ISBN-13: 9781846284588.
[79] F. L. Luo , [Link], « Renewable energy systems advanced conversion technologies and applications”, by
Taylor & Francis Group, LLC, 2013, ISBN: 978-1-4398-9110-0.
[80] C. Kong Tse, “Complex behavior of switching power converters”, CRC Press, 2004, ISBN 0-8493-
1862-9.
[81] H. R. Karshenas, H. Daneshpajooh, A. Safaee, P. Jain, A. Bakhshai, “Bidirectional DC-DC Converters
for Energy Storage Systems”, Energy Storage in the Emerging Era of Smart Grids, pp161-178, Prof.
Rosario Carbone (Ed.), ISBN: 978-953-307-269-2, InTech.
[82] El Fadil, H., Giri, F., & Guerrero, J. M. ,” Lyapunov based control of hybrid energy storage system in
electric vehicles”, In Proceedings of the American Control Conference (ACC), 2012, pp. 5005-5010,
Montreal, QC: IEEE Press.
[83] [Link], “ Développment d’un superviseur flou appliqué à la commande d’un moteur à induction “,
These de doctorat, université M’HAMED bougara-Boumerdes, 2013-2014.
[84] G. Abad, [Link], M. A. Rodrıguez, L. Marroyo, G. Iwanski, “Doubly Fed Induction Machine Modeling
And Control For Wind Energy Generation”, IEEE press, John Wiley & Sons, Inc.,2011, ISBN: 978-1-
118-10494-1.
[85] J. Rodríguez, S. Bernet, [Link], J. O. Pontt, [Link],” Multilevel Voltage-Source-Converter Topologies
for Industrial Medium-Voltage Drives”, IEEE Transactions on Industrial Electronics, Vol. 54, No. 6,
pp2930-2945 , 2007.
[86] G. Dreyfus,G. Dreyfus, J.-M. Martinez, M. Samuelides, M. B. Gordon, F. Badran, S.
Thiria, « Apprentissage statistique : Réseaux de neurones, Cartes topologiques, Machines à vecteurs
supports», Éditions Eyrolles, ISBN : 978-2-212-11464-5.
[87] L. Zadeh (1965) Fuzzy sets. Information and Control, Vol. 8, pp.338–353.
[88] C. Pappis, C. Siettos (2005) Fuzzy reasoning. In: Burke, E.K., Kendall, G. (eds.) Search Methodologies:
Introductory Tutorials in Optimization and Decision Support Techniques, chapter. 15, pp.437–474.
Springer.
[89] [Link], [Link], [Link] “Neuro-Fuzzy and Soft Computing : a computational approach to
learning and machine intelligence”,Prentice Hall, 1997,ISBN 0-13-261066-3.
[90] N. Siddique “Intelligent Control:A Hybrid Approach Based on Fuzzy Logic, Neural Networks and
Genetic Algorithms”,Springer , 2014, ISBN 978-3-319-02134-8.
[91] [Link], “Fuzzy logic in control systems : Fuzzy logic controller –Part I”, IEEE transactions on systems
, man, and cybernetics , vol 20, no.2, 1990.
[92] M.N. Cirstea, A. Dinu, J.G. Khor, M. McCormick, “Neural and Fuzzy Logic Control of Drives and
Power Systems”, Elsevier, 2002, ISBN 0 7506 55585.
[93] T. Takagi, M. Sugeno (1985) Fuzzy Identification of Systems and Its Applications to Modeling and
Control. IEEE Transactions on Systems Man and Cybernetics, Vol. 15, N.1, pp. 116-132.
[94] P. P.-Cruz, F. D. R.-Figueroa, “Intelligent Control Systems with LabVIEW”, Springer,2010, ISBN 978-
1-84882-683-0.
[95] E. Czogalat . J. Leski, “Fuzzy and Neuro-Fuzzy Intelligent Systems”, A Springer-Verlag Company,
2000, ISBN 978-3-662-00389-3.
[96] [Link], [Link], [Link]ège, « Les réseaux de neurones : Présentation et applications », Edition
Technip, 2007, ISBN :978-2-7108-0896-1.
[97] [Link] “Neural Networks : A comprehensive Foundation”, deuxième edition, Pearson Education,
Inc.,1999, ISBN:81-7808-300-0.
[98] O. GUENOUNOU, « Méthodologie de conception de contrôleurs intelligents par l'approche génétique-
application à un bioprocédé », Thèse de doctorat, Université De Toulouse, 2009.
[99] A. S. Poznyak, E. N. Sanchez, Wen Yu, “Differential neural networks for robust nonlinear control”,
2001, ISBN 981-02-4624-2.
[100] K. Karabacak, [Link], “Artificial neural networks for controlling wind–PV power systems: A
review”, Renewable and Sustainable Energy Reviews, Vol.29, pp.804–827, 2014.
166
[101] J. Liu, « Radial Basis Function (RBF) Neural Network Control for Mechanical Systems : Design,
Analysis and Matlab Simulation »,Springer, 2013, ISBN 978-7-302-30255-1.
[102] M. Ünal, A. Ak, V. Topuz, H. Erdal, “Optimization of PID Controllers Using Ant Colony and
Genetic Algorithms”, Springer, 2013, ISBN 978-3-642-32899-2.
[103] F.-J. Lin, Y.C. Hung, S.-Yi Chen, “FPGA-Based Computed Force Control System Using Elman
Neural Network for Linear Ultrasonic Motor”, IEEE Transactions On Industrial Electronics, Vol. 56,
No. 4, pp.1238 -1253, 2009.
[104] [Link], “ RESEAUX DE NEURONES” , Université de LAVAL,2004.
[105] [Link] Nguyen, “Approche neuromimétique pour l’identification et la commande des systèmes
électriques : application au filtrage actif et aux actionneurs Synchrones”, These de doctorat, Université
De Haute-Alsace, Mulhouse, 2010.
[106] Krose.B., P.van_der_Smagt, “An introduction to Neural Networks », Eighth edition , November
1996.
[107] W.Y. Han, S. M. Kim, S.J. Kim, C.G. Lee, “Sensorless Vector Control of Induction Motor using
Improved Self-Tuning Fuzzy PID Controller”, SICE Annual Conf. in Fukui-Japan, August 4-6 (2003)
3112-3117.
[108] F. BLAYO,M. VERLEYSEN, « Les réseaux de neurones artificiels », Presses Universitaires de
France, 1996, isbn 2 13 047355 5.
[109] M._Gupta,_L._Jin,_N._Homma, “Static and dynamic neaural networks: from fundamentals to
advanced theory”, IEEE Press , A JOHN WILEY SONS, INC.2003,ISBN 0-471-21948-7.
[110] T.O. Kowalska, F. Blaabjerg, J. Rodríguez, “Advanced and Intelligent Control in Power
Electronics and Drives,” Studies in Computational Intelligence, Vol. 531, Springer, 2014.
[111] M. Parvez , M.F.M. Elias , N.A. Rahim , N. Osman, « Current control techniques for three-phase
grid interconnection of renewable power generation systems: A review”, Solar Energy, vol.135, pp.29–
42, 2016.
[112] N. Chettibi, A. Mellit, “PSO based Direct Power Control of grid connected photovoltaic
system“, 4th International Conference on Electrical Engineering (ICEE), PP.1-6, DOI:
10.1109/INTEE.2015.7416732 , 2015.
[113] T. Esram, P. L. Chapman, “Comparison of Photovoltaic Array Maximum Power Point Tracking
Techniques”, IEEE Transactions on energy conversion, vol. 22, no. 2, 439-449, June 2007.
[114] K.H. Hussein, I. Muta , T. Hoshino, M. Osakada «Maximum photovoltaic power tracking: an
algorithm for rapidly changing atmospheric conditions » IEE Proc.-Gener. Transm Distrib, Vol. 142, No.
1 , pp.59-,64, January I995.
[115] A. Safari, S. Mekhilef,” Simulation and Hardware Implementation of Incremental Conductance
MPPT With Direct Control Method Using Cuk Converter”, IEEE Transactions On Industrial Electronics,
Vol. 58, No. 4, pp.1154-1161,April 2011.
[116] T. Noguchi, H. Tomiki, S. Kondo, T. Takahashi, "Direct Power Control of PWM Converter
Without Power-Source Voltage Sensors", IEEE Transactions on industry applications, vol. 34, No. 3,
pp.473-479, 1998.
[117] A. L. de Heredia Bermeo, “ Commandes avancées des systèmes dédiés a l'amélioration de la qualité
de l'énergie : de la basse tension a la montée en tension », Thèse de doctorat, Institut polytechnique de
Grenoble, 2006.
[118] M. Malinowski, [Link], M. P. Kazmierkowski, “Simple Direct Power Control of Three-Phase
PWM Rectifier Using Space-Vector Modulation (DPC-SVM)”, IEEE Transactions on industrial
electronics, vol. 51, no. 2, pp.447-454, 2004.
[119] M. Cichowlas, M. Malinowsk, M. Jasinsk ,M. P. Kazmierkowsk,” DSP Based Direct Power
Control for three-phase PWM Rectifier with Active Filtering Function”, IEEE International Symposium
on Industrial Electronics ISIE, 2003.
[120] M. Malinowski, M. P. Kazmierkowski, , S. Hansen, F. Blaabjerg, G. D. Marques, "Virtual-Flux-
Based Direct Power Control of Three-Phase PWM RectifIers", IEEE transactions on industry
applications, vol.37, No. 4, pp. 10 19-1027 , 200 1.
[121] A. Bouafla, J-P. Gaubert, F. Krim, "Predictive Direct Power Control of Three-Phase Pulsewidth
Modulation (PWM) RectifIer Using SpaceVector Modulation (SVM)", IEEE Transactions on Power
Electronics, Vol. 25, No. I ,pp. 228- 236, January 20 10.
[122] L. Zhang , G. Zheng , L. Jishun, "Direct Power Control Strategy of Three-phase PWM RectifIer",
Journal ofTnformation & Computational Science 10:9 , pp 2807-2813 ,10 June 20 13.
[123] M. Malinowski, S. Bernet, “ Simple Control Scheme of Three-Level PWM Converter Connecting
Wind Turbine with Grid”, [Link]
[124] 1. Kennedy, [Link], "Particle swarm optimization", IEEE International Conference on Neural
Networks, Perth, W A, vol.4, pp. 1942-1948, 1995.
[125] A. EL DOR, "Perfectionnement des algorithmes d'Optimisation par Essaim Particulaire :
Applications en segmentation d' images et en electronique ",These de doctorat, UNIVERSITÉ PARIS-
EST, 2012.
167
[126] V. Rajinikanth , K. Latha, “Tuning and Retuning of PID Controller for Unstable Systems Using
Evolutionary Algorithm”, International Scholarly Research Network, ISRN Chemical Engineering, Vol.
2012, ID 693545, doi:10.5402/2012/693545.
[127] S. S. Patnaik , A.K. Panda, “Particle SwarmOptimization and Bacterial Foraging Optimization
Techniques for Optimal Current Harmonic Mitigation by Employing Active Power Filter”, Applied
Computational Intelligence and Soft Computing, Hindawi Publishing Corporation, Vo1.20 12, Article
ID 897127.
[128] Y. H. Liu, S. C. Huang, 1. W. Huang, W. C. Liang, "A Particle Swarm Optimization-Based
Maximum Power Point Tracking Algorithm for PV Systems Operating Under Partially Shaded
Conditions", IEEE Transactions On Energy Conversion, Vol. 27, No. 4, pp. 1027-1035, 2012.
[129] C.C. Kao, C. W. Chuang, R. F. Fung , "The self-tuning PID control in a slider-crank mechanism
system by applying particle swarm optimization approach", Journal of Mechatronics, Vol. 16, No.8, pp.
5 13-522, 2006.
[130] Y. Shi, R.C. Eberhart, "Empirical study of particle swarm optimization", in Proceeding of the 1999
Congress on Evolutionary Computation, Washington, Vol.3, pp. 1945-1950, 1999.
[131] I. Takahashi and T. Noguchi, “A new quick response and high efficiency strategy of an induction
motor,” in Conf. Rec. IEEE-IAS Annu. Meeting, 1985, pp. 495–502.
[132] L. Zhong, M. F. Rahman, W. Y. Hu, and K. W. Lim,” Analysis of Direct Torque Control in
Permanent Magnet Synchronous Motor Drives”, IEEE Transactions on Power Electronics, Vol. 12,
No.3,pp. 528-536, May 1997.
[133] B. Bossoufi, M. Karim, A. Lagrioui, B. Bossoufi, S. Ioniţă “ Performance Analysis of Direct
Torque Control (DTC) for Synchronous Machine Permanent Magnet (PMSM) », 2010 IEEE 16th
International Symposium for Design and Technology in Electronic Packaging (SIITME), 23-26 Sep
2010, Pitesti, Romania.
[134] [Link], L. Zhong, M. F. Rahman, Y. Hu, “A Novel Direct Torque Control for Interior Permanent-
Magnet Synchronous Machine Drive With Low Ripple in Torque and Flux—A Speed-Sensorless
Approach”, IEEE Transactions On Industry Applications, Vol. 39, No. 6, pp. 1748 -1756, 2003.
[135] N. Khalil, “ Stratégies de commande et techniques intelligentes appliquées aux machines de type
synchrone », Thèse de doctorat,Universite Mentouri, Constantine,2010.
[136] H. Merabet Boulouiha , A. Allali , M. Laouer , A. Tahri , M. Denaï, A. Draou, « Direct torque
control of multilevel SVPWM inverter in variable speed SCIG-based wind energy conversion system”,
Renewable Energy, Vol. 80, pp.140-152, 2015.
[137] J. Zhang, M. Faz Rahman, [Link] “A New Scheme to Direct Torque Control of Interior
Permanent Magnet Synchronous Machine Drives for Constant Inverter Switching Frequency and Low
Torque Ripple”, IEEE, IPEMC 2006.
[138] C.-Yi Chen, C.-H. Hsu, S.H. Yu, C.F. Yang, H.H. Huang,” Cascade PI Controller Designs for
Speed Control of Permanent Magnet Synchronous Motor Drive Using direct torque approach”, IEEE
Fourth International Conference on Innovative Computing, Information and Control, pp. 938-941,2009.
[139] K. Chikh, M. Khafallah , A. Saâd, “Improved DTC Algorithms for Reducing Torque and Flux
Ripples of PMSM Based on Fuzzy Logic and PWM Techniques”, MATLAB – A Fundamental Tool for
Scientific Computing and Engineering Applications, Vol. 1. Chapitre 8, INTECH.
[140] A. Boubakir, F. Boudjema, [Link], “A Neuro-fuzzy-sliding Mode Controller Using Nonlinear
Sliding Surface Applied to the Coupled Tanks System”, International Journal of Automation and
Computing,Vol.6, no.1, pp.72-80, February 2009.
[141] R. Teodorescu, M. Liserre, P. Rodriguez, Grid converters for photovoltaic and wind power systems,
IEEE John Wiley & Sons, Ltd , 2011, ISBN 978-0-470-05751-3.
[142] M. A. GOLKAR, A. HAJIZADEH, “Control strategy of hybrid fuel cell/battery distributed
generation system for grid-connected operation”, Journal of Zhejiang University SCIENCE A, Vol.10,
no.4, pp. 488-496, 2009.
[143] N. Chettibi, A. Mellit, G. Sulligoi, A. Massi Pavan, “Fuzzy-based power control for
distributed generators based on solid oxide fuel cells”, International Conference on Clean
Electrical Power (ICCEP), PP. 580 - 585, DOI: 10.1109/ICCEP.2015.7177550, 2015.
[144] E. Unamuno , J. A. Barrena, Hybrid ac/dc microgrids—Part II: Review and classification of control
strategies, Renew. Sustain. Energy Rev. 52 (2015) 1123–1134.
[145] F. Nejabatkhah, Y.W. Li, “Overview of Power Management Strategies of Hybrid AC/DC
Microgrid”, IEEE Trans. Power Electronics, Vol 30, no.12,pp.7072-7089, 2015.
[146] X. Liu, P. Wang, P. C. Loh, “A Hybrid AC/DC Microgrid and Its Coordination Control,” IEEE
Trans. Smart Grid, vol. 2, no. 2, pp. 278- 286 , June 2011.
[147] A. Choudar , D. Boukhetala , S. Barkat , J.-M. Brucker, “A local energy management of a hybrid
PV-storage based distributed generation for microgrids”, Energy Conversion and Management, vol. 90,
pp. 21–33, 2015.
168
[148] [Link] , F. Pourfayaz, Optimal sizing of autonomous hybrid photovoltaic/wind/battery power
system with LPSP technology by using evolutionary algorithms, Solar Energy, vol. 115, pp. 471–483,
2015.
[149] M.S. Ismail, M. Moghavvemi, T.M.I. Mahlia, Genetic algorithm based optimization on modeling
and design of hybrid renewable energy systems, Energy Conversion Management , vol.85, pp.120-130,
2014.
[150] [Link] , [Link], L.M. Fernandez, F. Llorens, [Link] , ANFIS-Based Control of a Grid-
Connected Hybrid System Integrating Renewable Energies, Hydrogen and Batteries, IEEE Transaction
on Industrial Information, vol 10,no. 2, pp. 1107 – 1117, 2014.
[151] O. Erdinc, M. Uzunoglu, ‘The importance of detailed data utilization on the performance evaluation
of a grid-independent hybrid renewable energy system’, Hydrogen energy, vol.36, pp.12664-12677,
2011.
[152] H. M. Hasanien, S.M. Muyeen, “Speed control of grid-connected switched reluctance generator
driven by variable speed wind turbine using adaptive neural network controller”, Electric Power Systems
Research, vol. 84, pp. 206– 213,2012.
[153] C. H. Chen, C.-M. Hong, F.S. Cheng,” Intelligent speed sensorless maximum power point tracking
control for wind generation system”, Electrical Power and Energy Systems, vol. 42, pp. 399–407,2012.
[154] C.-M. Hong, C.-H. Chen , C.S. Tu, « Maximum power point tracking-based control algorithm for
PMSG wind generation system without mechanical sensors”, Energy Conversion and Management, vol.
69, pp. 58–67, 2013.
[155] T.C. Ou , C.M. Hong, “ Dynamic operation and control of microgrid hybrid power systems”,
Energy, Vol.66, pp.314-323, 2014.
[156] S. [Link], G. Scelba, “ANN based MPPT method for rapidly variable shading conditions”,
Applied Energy, vol. 145 , pp. 124–132, 2015.
[157] A. I. Dounis , P. Kofinas, G. Papadakis, C. Alafodimos, “A direct adaptive neural control for
maximum power point tracking of photovoltaic system”, Solar Energy, vol. 115, pp.145–165,2015.
[158] N. Chettibi , [Link] , G. Sulligoi , A. M. Pavan “Adaptive Neural Network-Based Control
of a Hybrid AC/DC Microgrid”, IEEE Transactions on Smart Grid, PP.1-13 DOI:
10.1109/TSG.2016.2597006, 2016.
[159] P. K. Dash, S. Mishra, and G. Panda, “A Radial Basis Function Neural Network Controller for
UPFC,” IEEE Trans. Power Syst., vol. 15, no. 4, November 2000.
[160] W.M. Lin, C.M. Hong, C.H. Chen, “Neural-Network-Based MPPT Control of a Stand-Alone
Hybrid Power Generation System,” IEEE Trans. Power Electron., Vol. 26, No. 12, pp. 3571-3581,
December 2011.
[161] A. Bouafia, J.P. Gaubert, F. Krim, “Predictive Direct Power Control of Three-Phase Pulsewidth
Modulation (PWM) Rectifier Using Space-Vector Modulation (SVM), ” IEEE Trans. Power Electron.,
Vol. 25, No. 1, pp.228-236, January 2010.
[162] P. Cortés, M. P. Kazmierkowski, R. M. Kennel,D. E. Quevedo, José Rodríguez, “Predictive Control
in Power Electronics and Drives”, IEEE Trans. Indus. Electron., Vol. 55, No. 12, pp. 4312-4324, 2008.
[163] N. Li ,Y. Ming, X. Dian-guo, “Predictive Current Control for Permanent Magnet Synchronous
Motor Based on deadbeat Control”, 7th IEEE Conf. on Industrial Electronics and Applications (ICIEA),
2012.
[164] D. Ould Abdeslam, P. Wira, D. Flieller, J. Mercklé, "Artificial Neural Networks to Control an
Inverter in a Harmonic Distortion Compensation Scheme," in ISIE'2008, Cambridge, UK, 2008,
pp.1873-1878, 2008.
[165] D. Biel, F. Guinjoan, E. Fossas, J. Chavarria, “Sliding-Mode Control Design of a Boost–Buck
Switching Converter for AC Signal Generation”, IEEE Trans. circuits and systems—I, Vol. 51, No. 8,
pp. 1539-1551, August 2004.
[166] M. Malinowski, M.P. Kazmierkowski, S. Hansen, F. Blaabjerg, G.D. Marques, Virtual-Flux-Based
Direct Power Control of Three-Phase PWM Rectifiers, IEEE Trans. Indus. Appl. 37 (2001)1019-1027.
[167] M. Malinowski, M. P. Kazmierkowski, A. M. Trzynadlowski, A Comparative Study of Control
Techniques for PWM Rectifiers in AC Adjustable Speed Drives, IEEE Trans. Power Electron. 18 No. 6
(2003) 1390- 1396.
[168] H. [Link], D. A. Khaburi, V. E. Torbati, Dead Beat Control of Three Phase PWM Rectifier
Using Virtual Flux Based Direct Power Control (DPC) and with no Line Voltage Measurements, 6th
conf. Power Electronics, Drives systems & Technologies PEDSTC (2015).
[169] J.A. Suul, A. Luna, P. Rodriguez, T. Undeland, Virtual-Flux-Based Voltage-Sensor-Less Power
Control for Unbalanced Grid Conditions, IEEE Transactions Power Electronics, Vol.27, No.9, pp.4071-
4087, 2012.
[170] R. Teodorescu, M. Liserre, P. Rodriguez, Grid converters for photovoltaic and wind power systems,
IEEE John Wiley & Sons, Ltd., 2011, ISBN:978-0-470-05751-3.
[171] L. A. Serpa, S. D. Round, J. W. Kolar, A Virtual-Flux Decoupling Hysteresis Current Controller
for Mains Connected Inverter Systems, IEEE Transactions Power Electronics, 22-5 (2007) 1766-1777.
169
[172] R.D. Brandt, F. Lin, “Adaptive interaction and its application to neural networks,” Information
Science, vol. 121, pp.201–215, 1999.
[173] X. Li, Z.Q. Chen, Z.Z. Yuan, “Generating chaos by an Elman network,” IEEE Transaction on
Circuit System I: Fundam. Theory Application, Vol. 48, pp.1126–1131, 2001.
[174] J. G. de Matos, F.S. F. e Silva, L.A. de S. Ribeiro, “Power Control in AC Isolated Microgrids With
Renewable Energy Sources and Energy Storage Systems”, IEEE Trans. Indus. [Link]. 62,
No.6,pp.3490-3498, June 2015.
[175] [Link],[Link],.[Link], “Fuzzy gain scheduling of PID controllers”, IEEE transactions on
systems man and cybernetics, vol.23,[Link].1392-1398,1993.
[176] A. I. Dounis, P. Kofinas, C. Alafodimos, D. Tseles, Adaptive fuzzy gain scheduling PID controller
for maximum power point tracking of photovoltaic system, Renewable Energy, vol. 60, pp. 202-214,
2013.
[177] [Link], Y. Yali, Z. Guosheng, S. Xiaoliang, Fuzzy Auto-adjust PID Controller Design of
Brushless DC Motor, Inter. Conf. Medical Physics and Biomedical Engineering, Physics Procedia, vol.
33 ,pp. 1533 – 1539, 2012.
[178] L. Gil-Antonio, M. B. Saldivar-Marquez, [Link]-Rodriguez, “ Maximum power point tracking
techniques in photovoltaic systems: A brief review”,13th International Conference on Power Electronics
(CIEP), 2016.
[179] M. [Link], L. Mathew, “Different Kinds of Maximum Power Point Tracking Control Method for
Photovoltaic Systems: A Review”, Springer, Arch Computat Methods Eng.,2016.
[180] [Link], Z. Salam, A critical evaluation on maximum power point tracking methods for partial
shading in PV systems, Renew. Sustain. Energy Rev. 47 (2015) 933–953.
[181] G. Tsengenes, G. Adamidis”, “A multi-function grid connected PV system with three level NPC
inverter and voltage oriented control”, Solar Energy, vol. 85, pp.2595–2610,2011.
[182] S.W. Kim, J.J. Lee, Design of a fuzzy controller with fuzzy sliding surface, Fuzzy Sets Syst. 71
(1995) 359-367.
[183] B. Allaoua, A. Laoufi, “A Novel Sliding Mode Fuzzy Control based on SVM for Electric Vehicles
Propulsion System” TerraGreen 13 International Conference - Advancements in Renewable Energy and
Clean Environment, Energy Procedia, vol.36 , pp.120 – 129, 2013.
[184] Y. Liu, Q. Zhang, C. Wang, N. Wang, “A control strategy for microgrid inverters based on adaptive
three-order sliding mode and optimized droop controls”, Electric Power Systems Research vol.117 ,
pp.192–201, 2014.
[185] R. Selvamuthukumaran, R. Gupta, “Rapid prototyping of power electronics converters for
photovoltaic system application using Xilinx System Generator”, IET Power Electron., 2014, Vol. 7, Iss.
9, pp. 2269–2278.
[186] Y. S. K. BABU, G. T. R. DAS, “The use of field programmable gate array (FPGA) in direct torque
control of induction motor”J Control Theory Appl 2013 11 (4) 642–650.
[187] S. Krim, S. Gdaim, A. Mtibaa, M. F. Mimouni, “Design and Implementation of Direct Torque
Control Based on an Intelligent Technique of Induction Motor on FPGA”, J Electr Eng Technol.2015;
10: 30-40.
[188] N. SGHAIER, [Link], S. GDAIM, B. BOSSOUFI, M. F. MIMOUNI, “Vector control for
induction motor with rotor flux oriented using xilinx system generator”, Journal of Theoretical and
Applied Information Technology, Vol. 54 No.2, August 2013.
[189] R. Faraji, A. Rouholamini, H. R. Naji, R. Fadaeinedjad, M. R. Chavoshian, “FPGA-based real
time incremental conductance maximum power point tracking controller for photovoltaic systems”, IET
Power Electron., 2014, Vol. 7, Iss. 5, pp. 1294–1304.
[190] S. Krim, S. Gdaim, A. Mtibaa, M. F. Mimouni “FPGA implementation of the Direct Torque
Control with constant switching frequency of induction motor”, 12th International Multi-Conference on
Systems, Signals & Devices, 2015.
[191] K. Murugesan, R. Muthu, S. Vijayenthiran, J.B. Mervin, “Prototype hardware realization of the
DSTATCOM for reactive power compensation”, Electrical Power and Energy Systems, vol.65, pp. 169–
178, 2015.
[192] B. Bossoufi , [Link] , A. Lagrioui , M. Taoussi , A. Derouich, « Observer backstepping control
of DFIG-Generators for wind turbines variable-speed: FPGA-based implementation,” Renewable
Energy, vol.81, 903-917, 2015.
[193] M. Knezevic, “An introduction to Xilinx System Generator”, http:/[Link].
[194] “Xilinx System Generator v2.1 Reference Guide”,
[Link]
[195] [Link], Z. Miao, L. Piyasinghe, L. Fan, “Real-time digital simulation-based modeling of a
single-phasesingle-stage PV system”, Electric Power Systems Research, vol.123, pp.85–91, 2015.
[196] K. Yenduri , [Link],” Maximum Power Point Tracking of Variable Speed Wind Turbines
With Flexible Shaft”, IEEE Transactions on sustainable energy, vol. 7, no. 3, pp.956-965, JULY 2016.
170
[197] [Link] , «The principal of field orientation on applied to the transvector closed loop control
system for rotating field machines”, Siemens , review, vol.34, pp.217-220, may 1972.
[198] M. Parvez, M.F.M. Elias , N.A. Rahim, N. Osman, “Current control techniques for three-phase grid
interconnection of renewable power generation systems: A review”, Solar Energy, vol. 135, pp. 29–42,
2016.
[199] Il-Song Kim, “Robust maximum power point tracker using sliding mode controller for the three-
phase grid-connected photovoltaic system”, Solar Energy, vol. 81, pp. 405–414, 2007.
[200] M.W. Naouar , B. Ben Hania , I. Slama-Belkhodja , E. Monmasson, A.A. Naassani, “FPGA-based
sliding mode direct control of single phase PWM boost rectifier”, Mathematics and Computers in
Simulation, vol.91, pp.249–261, 2013.
[201] G.-D. Wang, R.-J. Wai, Y. Liao, “Design of backstepping power control for grid-side converter of
voltage source converter-based high-voltage dc wind power generation system”,IET Renew. Power
Gener., Vol.7, Issue 2, pp. 118–133, 2013.
[202] C. Aouadi, A. Abouloifa, « Backstepping Based Control of PV system Connected to the Grid”,
International Journal of Computer and Information Technology (ISSN: 2279 – 0764) ,Vol. 03, Issue 05,
2014.
[203] K. Naoufel, M. Zazi, H. Mahmoudi, “Grid-Connected Photovoltaic System Using An Advanced
Backstepping Approach”, Journal of Theoretical and Applied Information Technology, Vol.79. No.2,pp.
330-337 , September 2015.
[204] J. Zhou, C. Wen, « Adaptive Backstepping Control of Uncertain Systems: Nonsmooth
Nonlinearities, Interactions or Time-Variations”, Springer-Verlag Berlin Heidelberg, 2008, ISBN 978-
3-540-77806-6.
[205] C.-X. Chen; Y.-X. Xie, Y.-H. Lan, “Backstepping Control of Speed Sensorless Permanent Magnet
Synchronous Motor Based on Slide Model Observer”, International Journal of Automation and
Computing, Vol. 12, Issue.2, pp.149-155, April 2015.
[206] B. Bossoufi , [Link], A. Lagrioui, M. Taoussi, A. Derouich, « Observer backstepping control
of DFIG-Generators for wind turbines variable-speed: FPGA-based implementation”, Renewable
Energy, vol. 81, pp. 903-917,2015.
171