Université Mohamed Khider de Biskra

Faculté des Sciences et de la Technologie

Département de génie électrique

MÉMOIRE DE MASTER
Sciences et Technologies
ELECTROTECHNIQUE
Réseau électrique

Réf. : ……………………….

Présenté et soutenu par :

MENZER Meryem

Le : lundi 25 juin 2018

Analyse des performances d’un réseau

électrique en présence des sources à
énergie renouvelable

Jury :

M. MAHDAD Belkacem MCA Université Mohamed Khider Président

M. SALHI Ahmed MCB Université Mohamed Khider Encadreur

M. ROUINA Abd elhafid MCA Université Mohamed Khider Examinateur

Année universitaire : 2017 – 2018

 ‫اﻟﺟﻣﮭورﯾﺔ اﻟﺟزاﺋرﯾﺔ اﻟدﯾﻣﻘراطﯾﺔ اﻟﺷﻌﺑﯾﺔ‬
République Algérienne Démocratique et Populaire
‫وزارة اﻟﺗﻌﻠﯾم اﻟﻌﺎﻟﻲ واﻟﺑﺣث اﻟﻌﻠﻣﻲ‬
Ministère de l’enseignement Supérieur et de la recherche scientifique

Université Mohamed Khider Biskra

Faculté des Sciences et de la Technologie

Département de Génie Electrique
Filière : ELECTROTECHNIQUE
Option : Réseau électrique
Mémoire de Fin d'Etudes
En vue de l’obtention du diplôme :

MASTER

Analyse des performances d’un réseau électrique en présence

des sources à énergie renouvelable
Présenté par : Avis favorable de l’encadreur :
MENZER Meryem SALHI Ahmed signature

Avis favorable du Président du Jury

Nom Prénom Signature

Cachet et signature
 ‫اﻟﺟﻣﮭورﯾﺔ اﻟﺟزاﺋرﯾﺔ اﻟدﯾﻣﻘراطﯾﺔ اﻟﺷﻌﺑﯾﺔ‬
République Algérienne Démocratique et Populaire
‫وزارة اﻟﺗﻌﻠﯾم اﻟﻌﺎﻟﻲ و اﻟﺑﺣث اﻟﻌﻠﻣﻲ‬
Ministère de l’enseignement Supérieur et de la Recherche scientifique

Université Mohamed Khider Biskra

Faculté des Sciences et de la Technologie

Département de Génie Electrique
Filière : Electrotechnique
Option : Réseau électrique

Analyse des performances d’un réseau électrique en présence des sources à

énergie renouvelable

Proposé et dirigé par : SALHI Ahmed

RESUMES (Français et Arabe)

Ce mémoire traite l’analyse de l’intégration des sources à énergie renouvelable dans le réseau électrique, plus
particulièrement l’énergie éolienne et l’énergie solaire. On étudie l’impact de ces dernières sur la stabilité de
réseau. L’étude doit être faite avec CPF pour l’analyse de la stabilité statique considérant l’emplacement des
sources renouvelables, puis la création d’un défaut et déterminer le TCID pour voir le fonctionnement de réseau
avec cette intégration pour l’analyse dynamique. Cette étude est appliquée sur les réseaux test IEEE 14 JB et IEEE
30 JB.

Mots clés : énergie renouvelable, éolienne, solaire, CPF, la stabilité, TCID

‫اﻟﻣﻠﺧص‬
‫ ﺳﻧدرس ﺗﺄﺛﯾر اﻷﺧﯾرﯾن ﻋﻠﻰ ﺛﺑﺎت واﺳﺗﻘرار‬،‫ھذه اﻟﻣذﻛرة ﺗﻌﺎﻟﺞ ﺗﺣﻠﯾل دﻣﺞ ﻣﺻﺎدر اﻟطﺎﻗﺎت اﻟﻣﺗﺟددة ﻓﻲ اﻟﺷﺑﻛﺔ اﻟﻛﮭرﺑﺎﺋﯾﺔ ﺑﺎﻟﺗﺣدﯾد طﺎﻗﺔ اﻟرﯾﺎح واﻟطﺎﻗﺔ اﻟﺷﻣﺳﯾﺔ‬
‫ ﻟﻧﻼﺣظ‬TCID ‫ ﺛم ﺧﻠق ﺧﻠل ﻓﻲ ھذه اﻟﺷﺑﻛﺔ واﺳﺗﺧراج‬،‫ ﻓﻲ اﻟﺗﺣﻠﯾل اﻟﺳﺎﻛن ﺑﺎﻋﺗﺑﺎر ﻣﻛﺎن ﺗواﺟد ﻣﺻﺎدر اﻟطﺎﻗﺎت اﻟﻣﺗﺟددة‬CPF ‫اﻟدراﺳﺔ ﺳﺗﻛون ﺑواﺳطﺔ‬. ‫اﻟﺷﺑﻛﺔ‬
.‫ ﻗﺿﯾب ﺗﺟﻣﯾﻊ‬IEEE 30 ‫ و‬IEEE14 ‫ ھذه اﻟدراﺳﺔ طﺑﻘت ﻋﻠﻰ ﺷﺑﻛﺔ اﻹﺧﺗﺑﺎر‬.‫ﻋﻣل ھذه اﻟﺷﺑﻛﺔ ﻣﻊ ھﺎذ اﻹدﻣﺎج ﻓﻲ اﻟﺗﺣﻠﯾل اﻟدﯾﻧﺎﻣﯾﻛﻲ‬

،‫ طﺎﻗﺔ اﻟﺮﯾﺎح‬،‫ اﻟﻄﺎﻗﺔ اﻟﻤﺘﺠﺪدة‬،‫اﻟﻄﺎﻗﺔ اﻟﺸﻤﺴﯿﺔ‬,CPF ،‫ اﻻﺳﺘﻘﺮار‬،TCID : ‫اﻟﻜﻠﻤﺎت اﻟﻤﻔﺘﺎﺣﯿﺔ‬

 TABLE DE MATIERE
Titres Pages

Remerciement…………………………………………………………………………….. I

Liste des tableaux………………………………………………………………………… II

Liste des figures ………………………………………………………………………….. III

Liste des Symboles et Acronymes ………………………………………………………. VI

Introduction générale …………………………………………………………………… 1

Chapitre I : Ecoulement de puissance

1. Introduction ………………………………………………………………………... 3
2. C’est quoi l’écoulement de puissance ?.................................................................. 4
3. Le but et l’objectif de l’écoulement de puissance………………………………… 4
4. Les équations de l’écoulement de puissance ……………………………………... 5
4.1. Equation de l’écoulement de puissance dans les lignes………………………... 7
4.2. Quelques observations sur l’écoulement de puissance…………….................. 7
5. Classification des variables de l’écoulement de puissance……………………… 9
5.1. Variable de perturbation (P)……………………………………………………. 9
5.2. Variable d’état (X)……………………………………………………………… 9
5.3. Variable de contrôle (U)………………………………………………………... 10
6. Classification des jeux de barres………………………………………………….. 10
6.1. Jeu de barre de charge (PQ)……………………………………………………. 10
6.2. Jeu de barre de contrôle (PV)………………………………………………….. 10
6.3. Jeu de barre de référence (| | )………………………………………………. 10
6.4. Jeu de barre de contrôle (PQ)………………………………………………….. 11
7. Résoudre le problème de l’écoulement de puissance……………………………. 11
7.1. La méthode itérative Newton-Raphson………………………………………… 11
7.2. Méthode de NR- modifié (fast-decoupled)…………………………………….. 13
8. Ecoulement de puissance continue………………………………………………... 13
8.1. Facteur de charge……………………………………………………………….. 13
8.2. Calcul de continu de l’écoulement de puissance……………………………….. 14
a. Pas prédicteur……………………………………………….......................... 15
b. Pas correcteur……………………………………………….......................... 16
9. Conclusion …………………………………………………………………………... 17

Chapitre II : Energies renouvelables

1. Introduction ………………………………………………………………………... 18
2. Energie renouvelable………………………………………………………………. 18
2.1. Evolution de la production en énergies renouvelables…………………………. 18
3. Différente types des énergies renouvelables……………………………………… 19
3.1. Hydroélectrique………………………………………………………………… 19
3.2. Solaire………………………………………………………………………….. 20
3.3. Eolienne ………………………………………………………………………... 21
3.4. Biomasse……………………………………………………………………….. 22
3.5. Géothermie …………………………………………………………………….. 22
3.6. Energie marine …………………………………………………………………. 23
4. Principaux sources à énergie renouvelables……………………………………… 24
4.1. Energie éolienne ……………………………………………………………….. 24
4.1.1. Evolution de l’énergie éolienne dans le monde…………………………. 25
4.1.2. Types d’éoliennes………………………………………………………... 26
4.1.2.1. Eolienne à axe horizontale……………………………………….. 26
4.1.2.2. Eolienne à axe verticale …………………………………………. 26
4.2. Energie solaire………………………………………………………………….. 27
4.2.1. Différents types d’énergie solaire……………………………………….. 27
4.2.1.1. L’énergie solaire thermique……………………………………… 27
a. Capteur parabolique…………………………………………….. 27
b. Capteur cylindro-parabolique…………………………………... 28
4.2.1.2. L’énergie solaire photovoltaïque………………………………… 28
5. Conclusion………………………………………………………………………….. 29
 Chapitre III : Intégration des sources à énergie renouvelable
1. Introduction ………………………………………………………………………... 30
2. Le besoin de l’intégration dans le réseau électrique …………………………….. 30
3. Qualité de la tension et perturbations du réseau électrique ……………………. 30
3.1. Qualité de la tension……………………………………………………………. 30
3.2. Différents perturbations dans le réseau électrique …………………………….. 31
3.2.1. Creux de tension ………………………………………………………… 31
a. Harmonique……………………………………………………....... 31
b. Inters harmoniques………………………………………………… 31
c. Infra-harmoniques…………………………………………………. 31
d. Conséquences des harmoniques…………………………………… 31
e. Les principales sources d’harmonique…………………………….. 31
3.2.2. Surtension ……………………………………………………………….. 31
3.2.3. Variation et fluctuation de tension ……………………………………… 31
3.2.4. Déséquilibre……………………………………………………………... 31
4. Impact de la production d’énergie renouvelable sur le réseau électrique …….. 31
4.1. Impacte sur la qualité de service………………………………………………... 33
4.2. Impact sur la continuité de service ……………………………………………. 33
4.3. Impact sur le plan de protection………………………………………………... 33
4.3.1. Courant de court-circuit…………………………………………………. 33
4.3.2. Influence sur la sélectivité des protections………………………………. 34
4.4. Impact sur la tenue de fréquence et de tension…………………………………. 34
4.4.1. Impact sur la tenue de fréquence………………………………………… 34
4.4.2. Impact sur la tenue de tension…………………………………………… 35
5. Intégration en Algérie……………………………………………………………… 35
5.1. Energie éolienne dans l’Algérie ……………………………………………….. 36
5.1.1. Les ressources en énergie éolienne ……………………………………... 36
5.1.2. Le vent en Algérie……………………………………………………….. 36
5.1.3. La ferme éolienne d’Adrar………………………………………………. 37
5.2. Energie solaire en l’Algérie…………………………………………………….. 38
5.2.1. Atlas solaire de l’Algérie ……………………………………………….. 38
5.2.1.1. Projet « Atlas 1 »………………………………………………… 38
5.2.2. Potentiel solaire………………………………………………………….. 39
6. Conclusion …………………………………………………………………………. 40

Chapitre IV : Analyse des résultats de stabilité

1. Introduction……………………………………………………………………... 41
2. La stabilité statique de tension ………………………………………………… 42
3. PSAT (Power Système Analysis Toolbox)…………………………………….. 42
3.1.C’est quoi le PSAT ?..................................................................................... 43
3.2.Utilisation de PSAT…………………………………………………………. 43
4. Ecoulement de puissance continue (continuation power flow « CPF »)……. 45
5. Analyse de la stabilité statique de réseau électrique …………………………. 45
5.1. Analyse de RE de 14 JB……………………………………………………... 45
5.1.1. Analyse de réseau avec la méthode de CPF………………………….. 48
5.2. Analyse de RE de 30 JB…………………………………………………… 56
5.2.1. Analyse de RE avec la méthode de CPF……………………………... 58
6. Analyse de la stabilité dynamique du réseau électrique …………………….. 65
6.1. Le temps critique d’isolement le défaut (TCID)……………………………. 65
6.2. Réseau teste (Réseau d’étude)………………………………………………. 65
6.3.Représentation la stabilité dynamique du RE de 30JB (en cas de défaut)…… 65
6.3.1. Etapes d’analyse……………………………………………………… 65
7. Intégration des sources à énergie solaire dans le réseau de 30JB…………. 68
8. Conclusion …………………………………………………………………….. 72
Conclusion générale …………………………………………………………………… 73
Bibliographie …………………………………………………………………………... 74
 Remerciements

Tout d’abord, nous remercions Dieu le tout puissant de nous avoir donné le courage

et la patience durant toutes ces années d’études.

Nous tenons à exprimer notre profonde gratitude et nos reconnaissances envers

Mr. Ahmed SALHI de nous avoir d'abord proposé ce thème, pour le suivi continuel tout le

long de la réalisation de ce mémoire et pour les conseils qu’il n’a cessé de nous donner et les

remarques aussi nous remerciement Mr NAIMI Djemai pour son remarques, informations et

son conseilles.

Nos remercîments vont aussi au Mr. MAHDAD Belkacem pour avoir accepté de

présider ce jury de soutenance et au Mr. ROUINA Abd-elhafid d’avoir accepté d’examiner

ce travail.A travers ce mémoire, j’adresse mes reconnaissances à tous mes enseignants

qui ont contribué à ma formation depuis la première classe de primaire jusqu’à

aujourd’hui, ainsi qu’à tous la promotion 2018 de la spécialité réseaux électrique.

 Liste des tableaux

Tableau pages

Tableau 3.1 Potentiel solaire dans l’Algérie 40

Tableau 4.1 Caractéristiques du réseau 14 JB 46
Tableau 4.2 les résultats de l’écoulement de puissance à l’état initial du modèle 46
de 14JB
Tableau 4.3 Les résultats de l’écoulement de puissance l’état à la marge de 49
stabilité-modèle 14JB.
Tableau 4.4 Les réserves de la puissance réactive des générateurs 50
Tableau 4.5 La différence entre l’état initial et l’état à la marge de stabilité en ce 51
qui concerne le module de tension dans les JB les plus sensibles du
modèle de 14JB.
Tableau 4.6 Résultats de simulation du premier scénario du réseau 14JB à la 53
marge de stabilité
Tableau 4.7 Résultats de simulation de deuxième scénario du réseau 14JB à la 54
marge de stabilité.
Tableau 4.8 Caractéristiques du réseau 30JB. 56
Tableau 4.9 les résultats de l’écoulement de puissance à l’état initial du modèle 57
de 30JB.
Tableau 4.10 Résultats de l’écoulement de puissance à la marge de stabilité du 58
modèle de 30JB
Tableau 4.11 La déférence entre l’état initial et l’état à la marge de stabilité du 60
module de tension dans les JB le plus sensibles du modèle de 30JB.
Tableau 4.12 résultats de simulation du premier scénario du réseau 30JB 61
Tableau 4.13 Résultats de simulation de deuxième scénario du réseau 30JB 63
Tableau 4.14 Classification des scénarios optimaux d’analyse statique. 64
Tableau 4.15 Classification le TCID sans intégration des DG. 66
Tableau 4.16 TCID de JB06 après l’intégration 67
Tableau 4.17 Classification de TCID après l’intégration de Photovoltaïque avec 71
l’éolienne

II
 Liste des figures

Fig.1.1 : Vue d’ensemble des composantes du réseau électrique 3

Fig.1.2 : jeu de barres quelconque avec génération, charge et ligne de 5
transmission
Fig.1.3 : calcul continu de l’écoulement de puissance: Pas prédicteur par les 16
moyens de la tangente
Fig1.4 : Calcul continu de l’écoulement de puissance: pas correcteur obtenu 16
par les moyens d'intersection perpendiculaire.
Fig.1.5: Calcul continu de l’écoulement de puissance: pas correcteur obtenu 17
par les moyens de la paramétrisation locale.
Fig.2.1 : Evolution de la consommation d’énergie mondiale. 19
Fig.2.2 : Centrale hydroélectrique 20
Fig.2.3 : une centrale solaire 21
Fig.2.4 : Différents types des énergies renouvelables 24
Fig.2.5 : Une ferme d’éoliens. 25
Fig.2.6 : Capacité mondiale installé d’éolienne 2001-2017. 26
Fig.2.7 : éolienne type horizontale 26
Fig.2.8 : éolienne type verticale 26
Fig.2.9 : le fonctionnement d’un capteur parabolique. 28
Fig.2.10 : Capteur cylindro parabolique 28
Fig.2.11 : Cellule PV 29
Fig.2.12 : Panneau PV 29
Fig.2.13 : Champ PV 29
Fig.3.1 : Flux d'énergie sur un réseau en présence de GED 32
Fig.3.2 : Un défaut sur la ligne 2 provoque la coupure inutile de la ligne 1 34
Fig.3.3 : Atlas de la vitesse moyenne du vent de l'Algérie estimée à 10 m du 37
sol.
Fig.3.4 : La ferme éolienne de Kabertène à 75 km au nord est d’Adrar 38
Fig4.1 : Classification des types de stabilité des réseaux électriques. 42
Fig4.2 : page d’accueil de PSAT 43
Fig4.3 : Bibliothèque du Simulink. 44
Fig. 4.4: Le schéma de topologie du modèle de 14-JB 45

III
 Liste des figures
Fig.4.5 : Profile de module de tension à l’état initial 47
Fig.4.6 : La modification du paramètre de CPF 48
Fig.4.7 : La comparaison entre V init et V cpf de 14JB 50
Fig.4.8 : Histogramme de variation du module de tension entre l état initiale 51
et l état à la marge de stabilité de14JB.
Fig.4.9 : Le module de tension au niveau des JB sensibles en fonction de  - 51
modèle 14 JB
Fig4.10 : différents scénarios d’étude de réseau 14JB 52
Fig.4.11: paramètre de source auxiliaire renouvelable de forme PV du modèle 53
14JB.
Fig4.12 : profile de tension des trois premiers scénarios de 20% du réseau 54
14JB
Fig4.13 : profile de tension des quatre derniers scénarios de 20% du réseau 54
Fig4.14 : profile de tension des trois premiers scénarios de 25% du 14JB 55
Fig4.15 : profile de tension des quatre derniers scénarios de 25% du 14JB 55
Fig.4.16: Schéma topologique du modèle de 30JB 56
Fig.4.17 : La tension de modèle de 30JB à l’état initial. 58
Fig.4.18 : La comparaison entre le module de tension initial et le module de 59
tension à l’état de la marge de stabilité de réseau 30JB.
Fig.4.19 : Histogramme de variation du module de tension entre l état initiale 60
et l état à la marge de stabilité du réseau 30JB.
Fig.4.20 : Le module de tension au niveau les JB sensibles en fonction de  - 60
modèle 30JB
Fig4.21 : différents scénarios d’intégration d’énergie renouvelable du réseau 61
30JB.
Fig4.22 : profile de tension des trois premiers scénarios de 20% du réseau 62
30JB
Fig4.23 : profile de tension des quatre derniers scénarios de 20% du réseau 62
30JB
Fig4.24 : profile de tension des trois premiers scénarios de 25% du réseau 63
30JB
Fig4.25 : profile de tension des quatre derniers scénarios de 25% du réseau 64
30JB

IV
 Liste des figures
Fig4.26 : modèle de la machine à double alimentation sous PSAT 66
Fig4.27 : variation de vitesse angulaire du rotor avec le temps de défaut du 67
20% de Pd tot
Fig4.28 : variation de vitesse angulaire du rotor avec le temps de défaut du 68
25% de Pd tot.

Fig.4.29 : La vitesse angulaire de synchronisme avec Photovoltaïque sans 69

défaut
Fig 4.30 : La vitesse angulaire de synchronisme en présence de PV du 20% de 70
Pd tot.
Fig.2.31 : La vitesse angulaire de synchronisme en présence de PV du 25% de 71
Pd tot.

V
 Liste d’acronymes et symboles

Liste d’acronymes
N-R Newton-Raphson Irriguée
KWh Kilo-Watt-heurs CDER Centre de Développement des
CPF Continuation Power Flow Energies Renouvelables
TW Térawatts JB Jeu de barre
GED Génération d’Energie TCID Temps Critique d’Isolement le
Dispersée(Distribuée Défaut
ASA l’Agence Spatiale Allemande IEEE Institute of Electric al and
TWh/an Térawatts par an Electronics Engineers
SAT Superficie Agricole Totale PSAT Power System Analysis Toolbox
SAU Superficie Agricole Utilisée RE Réseau Electrique
SAUI Superficie Agricole Utilisée MVA Méga volt ampère

Liste des symboles

f Fréquence
∆ Pertes de puissance actives
∆ Pertes de puissance réactive
Puissance transitée
Puissance transitée maximale
, Puissance générée
, Puissance demandée
Admittance série
Admittance parallèle ou shunt
Vi Tension de jeu de barre i
Tension de jeu de barre k
L’admittance de la ligne entre les deux jeux de barres i et k
La contribution shunt du courant au jeu de barres i

Angle de phase

VI
 Liste d’acronymes et symboles

module de la tension du jeu de barres i

│ │
puissance active injectée au jeu de barres i

puissance réactive injectée au jeu de barres i

puissance réactive demandée au jeu de barres i

P Variable de perturbation
U Variable de contrôle
X Variable d’état
La différence des phases entre les jeux de barres i et j
+ La partie réelle et imaginaire de la matrice
J matrice de jacobienne
Equation de l’écoulement de puissance
Variables dépendant
λ Facteur de charge
Puissance active du générateur
Puissance active de charge
Fourniture de l’offre (Supply bids)
Demande d’offre (demand bids)
λ max Le facteur de charge maximal que le réseau peut supporter sans perdre sa stabilité
∆ Le différent du module de tension entre l’état initial et l’état à la marge de stabilité
Puissance intégrée
Puissance demandée totale

VII
 INTRODUCTION GENERALE

A cause de l’augmentation de la population et le développement technologique, la

consommation d’énergie électrique n’a pas cessé d’accroitre nécessitant un développement sur le
plan de la production de l’énergie électrique. Malheureusement, ce développement nécessite la
construction de nouvelles centrales électriques de plus avec des dépenses de plus en plus
importantes sur l’énergie électrique produite, accompagné par un dégagement important des gaz
toxiques issue de la combustion de la matière première pour libérer l’énergie calorifique
entrainant une pollution de l’atmosphère, laissant des traces sur le plan de l’épuisement des
matières premières énergétiques (pétrole, charbon,…etc.).

Dans ces dernières années, le monde à connu un développement des énergies renouvelables
plus particulièrement l’énergie éolienne et l’énergie solaire. Les recherches scientifiques
élaborées pour exploiter ces énergies ont connu une révolution sur le plan énergétique dans le
monde entier. La production d’une énergie propre présente plusieurs avantages comme la
propreté de l’environnement, l’augmentation de la durée de vie des réserves fossiles en
minimisant la pollution atmosphérique causée par le dégagement des gaz toxique à partir des
centrales classiques. D’autre part, ces énergies renouvelables présentes des inconvénients en ce
qui concerne leur intégration. Le pourcentage d’intégration dans le réseau de cette énergie est
limité pour préserver la stabilité du système énergétique et le prix de KWh reste relativement
élevé.

Pour étudier les effets de l'intégration des énergies renouvelables dans le réseau sur la
stabilité et la qualité de l’énergie électrique du système énergétique, nous proposons dans ce
mémoire l’analyse des performances des réseaux électriques en présences des sources à énergies
renouvelables. Pour cet objectif, le mémoire sera organisé en quatre chapitres selon les points
traités comme suit :

Dans le premier chapitre .le problème du l’écoulement de puissance sera détaillé ainsi que les
méthodes de résolution notamment la méthode de Newton-Raphson qui nous aide à résoudre le
problème de l’écoulement de puissance pour qu’on puisse analyser notre réseau électrique. Une
partie sera consacré à la stabilité statique.

Par contre la totalité du deuxième chapitre est consacrée aux l’énergie renouvelables et les
différents ressources.

Le troisième chapitre vient pour l’intégration ces énergies renouvelables au réseau électrique, le
besoin d’intégration et ces différents impacts de ce dernier.

1
 INTRODUCTION GENERALE

Le dernier chapitre traite l’application que nous présentons et les modèles d’étude des réseaux
comme 14-30JB. Ensuite, on fait l’analyse statique puis dynamique pour retirer les résultats à
partir la simulation avec PSAT qui fonctionne sous MATLAB ver. R2009b, ainsi que leurs
interprétations.

Finalement une conclusion générale sera établir pour mettre fin à ce mémoire.

2
Chapitre I [ECOULEMENT DE PUISSANCE]

1. Introduction

Un réseau électrique compose d’un part des générateurs (la production) et d’autre part
des charges électriques à différent niveau de tension (les consommateurs d’énergie
électrique).

Les capacités de production des différents générateurs étant connues, comme calculer
l’état du réseau électrique complet, c.à.d. les différents paramètres et grandeurs de réseau
électrique courant, tension et puissance.

Ce problème général connu sous le nom de l’écoulement de puissance (Power Flow)

c’est un problème d’équilibre entre la production et la consommation. Pour résoudre le
problème de l’écoulement de puissance on été fait par des méthodes approfondies, utilisant les
méthodes numériques appliquées en ce domaine.

Dans ce chapitre, on va utiliser la méthode de Newton-Raphson (NR) pour résoudre ce

problème.

Fig 1.1 : Vue d’ensemble des composantes du réseau électrique [MOS 08]

3
Chapitre I [ECOULEMENT DE PUISSANCE]

2. C’est quoi l’écoulement de puissance ?

L’étude de l'écoulement de puissance consiste à calculer les grandeurs électriques

(tensions) du réseau électrique pour des extrémités spécifiées et des conditions données
aux jeux de barres.

3. Le but et l’objectif de l’étude de l’écoulement de puissance

L'étude de l'écoulement de puissance permet d'avoir la solution des grandeurs

d'un réseau électrique en fonctionnement normal, équilibré en régime permanent. Ces
grandeurs sont les tensions aux nœuds, les puissances injectées aux nœuds et celles qui
transitent dans les lignes. Les pertes et les courants s'en déduisent [GUP 05].

Le but de calcul des écoulements de puissance permet d’investiguer les points suivants :

 la détermination des tensions en tout point du réseau

 la détermination des puissances actives et réactives
 l’effet de la modification de la topologie du réseau
 l’étude les pertes (perte d’un générateur, d’une ligne ou autre)
 l’optimisation du fonctionnement du réseau
 l’optimisation des pertes

L’objectif fondamental des opérateurs des réseaux électriques est d’assurer la continuité de
service (alimentation de leurs consommateurs), tout en respectant plusieurs exigences :

 Les tensions et la fréquence de réseau doivent d’être dans les limites acceptables
Vmin< V< Vmax = = fmin <f <f max
 La forme d’onde du courant et de la tension doit être sinusoïdale et périodique
(éviter les harmoniques).
 Avoir le minimum des pertes de puissance P min et Qmin
 Eviter le dépassement du Imax ou (Sij max) puissance transitée maximale Sij< Sijmax
 Les lignes de transmission doivent fonctionner dans leurs limites thermiques et de
stabilité. (Ptr < Pmax pour la dégagement totale de chaleur et plus durée de vie des
câbles).
 La durée de coupure de l’alimentation doit être ou minimum possible.

4
Chapitre I [ECOULEMENT DE PUISSANCE]

 En outre, à cause de la dérégulation (ouverture) du marché d’électricité, le courant de

kWh doit être le plus bas possible.

Pour évaluer toute ces conditions, on doit procéder à des méthodes qui déterminent
l’état des grandeurs électriques.

4. Les équations de l’écoulement de puissance [MOS 08]

Dans le cas général, pour déterminer les équations de l’écoulement de puissance, on

considère que les puissances au jeu de barres i sont équilibrées, donc le bilan des puissances
électriques à un jeu de barres i d’un réseau électrique de n jeu de barres est la différence entre
la puissance générée et la puissance demandée ; au niveau du même jeu de barres.

= − = − + ( − ) (1.1)

Les points de connexion des branches des réseaux sont codés par des numéros désignant
chaque jeu de barres. Ces nombres spécifient les arrivées des lignes de transmission et des
transformateurs. Les numéros sont utilisés pour identifier les types de jeu de barres et
l’emplacement des éléments du réseau (condensateurs, inductances shunts et impédances).

Les lignes de transmission connectent le jeu de barres i aux autres jeux de barres k dans le
réseau électrique. Un jeu de barres peut être connecté au maximum à (n-1) jeu de barres. On
peut représenter chacune de ces lignes par un schéma équivalent avec une admittance série
et une admittance parallèle « Si une ligne n’existe pas l’admittance sera égale à
zéro » (fig I.2).

Fig.1.2 : jeu de barres quelconque avec génération, charge et ligne de transmission [MOS 08]

5
Chapitre I [ECOULEMENT DE PUISSANCE]

La construction du modèle mathématique, par lequel on peut prédire la forme de

l’écoulement de puissance, est essentielle pour la compréhension du mécanisme de
l’écoulement de puissance. Dans la formulation des équations de l’écoulement de puissance
l’une des deux matrices (admittance ou impédance) est nécessaire.

La puissance apparente injectée au jeu de barres i est donnée par l’équation

∗
= où est le courant qui entre dans le jeu de barres i. ce courant est composé de deux
composantes, une composante qui s’écoule à travers l’admittance shunt, et une
composante

( − ) qui s’écoule à travers l’admittance série de la ligne équivalente du réseau. D’après

la loi de Kirchhoff appliquée à un jeu de barres, pour une seule phase, nous avons l’équation
du courant :

∗
= ∗ = ∑ +∑ ( + )= ∑ + +∑ (− )

i=1,…n (1.2)

On peut écrire l’équation (1.2) sous la forme suivante :

∗
= ∗ = + + ⋯+ + ⋯+ i=1,…n (1.3)

Où =∑ ( + ) = =−

Du fait que les puissances sont connues et les courants sont inconnus, donc les
équations prendront des forme non linéaires :

∗ ∗
= − = ∑ = 1, … (1.4)

Séparons la partie réelle et imaginaire de l’équation (1.4) on aura les équations (1.5) et
(1.6) qui expriment l’équation des puissances actives et réactives au jeu de barres i :

=∑ | | ( + sin )= = 1, … (1.5)

=∑ | | − = = 1, … (1.6)

6
Chapitre I [ECOULEMENT DE PUISSANCE]

4.1. Equation de l’écoulement de puissance dans les lignes [MOS 08]

La détermination des valeurs des puissances dans les lignes de transport est
indispensable afin de localiser les lignes électriques surchargées, de calculer la valeur de
pertes de puissance et de déterminer les contraintes qui forcées.

Le courant de la branche entre les deux jeux de barres i et k, qui a le sens positif de i vers k est
donné par :

= ( − )+ (1.7)

Avec : : L’admittance de la ligne entre les deux jeux de barres i et k

: L’admittance shunt au jeu de barres i

: La contribution shunt du courant au jeu de barres i

∗
− = (1.8)
∗ ∗
− = ( − ) + (1.9)
∗( ) ∗
− = − + (1.10)

La valeur des pertes de puissance dans la ligne entre les deux jeux de barres i et k est
la somme algébrique de la répartition des puissances déterminées à partir des relations (1.9) et
(1.10)

4.2. Quelques observations sur les équations de l’écoulement de puissance [MOS 08]

Les observations concernant les équations du système (1.5) et (1.6) sont :

1. Les équations sont algébriques car elles représentent un système qui travaille en
régime permanent.
2. Les équations son non linéaires, donc il est très difficile d’avoir des solutions
analytiques, mais on peut obtenir facilement une solution numérique.
3. L’équilibre des puissances actives donne :

∑ =∑ +∑ (1.11)

Où : les pertes de puissances actives

4. De même l’équilibre de la puissance réactive

7
Chapitre I [ECOULEMENT DE PUISSANCE]

= + (1.12)

Où : les pertes de puissances réactives

5. : sont des équations en fonction des tensions

6. Dans les équations du système (1.5) et (1.6) les angles de phase , apparaissent
sous forme de différence − et non séparément.
7. D’après les équations du système (1.5) et (1.6) à chaque jeu de barres i nous avons à
considérer six variables : , , , . deux de ces variables sont toujours
spécifiées, , . pour résoudre les équations (1.5) et (1.6), nous devons préciser
deux autres variables supplémentaire. Du point de vue mathématique, le choix est
indifférent ; pourtant d’un point de vue physique nous ne pouvons spécifier que deux
variables sur lesquelles nous avons un moyen d’action. Nos choix sont influences
jusqu’à un certain degré par le type de dispositifs connectés à chaque jeu de barres.
Nous distinguons alors trois types de jeux de barres : jeu de barres du type P, Q ou
charge qui ne comprenant pas de générateur, jeu de barres du type P, V que comporte
certainement au moins un générateur et jeu de barres du type V, qui est toujours
connecté un générateur considéré comme source de tension.
8. Les réseaux électriques réels peuvent contenir des certains de jeu de barres, des
générateurs et des milliers des branches, donc pour un système de n jeux de barres on
aura :
 n variables :
 n variables : (angle de phase de la tension du jeu de barres i)
 n variables : (puissance active injectée au jeu de barres i)
 n variables : (puissance réactive injectée au jeu de barres i)
 n variables : (puissance active demandée au jeu de barres i)
 n variables : (puissance réactive demandée au jeu de barres i)

8
Chapitre I [ECOULEMENT DE PUISSANCE]

5. Classification des variables de l’écoulement de puissance

Les (6n) variables de l’écoulement de puissance sont regroupées en trois vecteurs, le

vecteur de perturbation constitue de et car elles ne sont pas contrôlables et dépendent
seulement des consommateurs, le vecteur de contrôle constitue de la puissance active et
réactive de chaque générateur en service ( ) et enfin le vecteur d’état constitue des
tensions des jeux de barres (| | ). [MOS 08].

5.1.Variables de perturbation (P)

D’après les variables des équations de l’écoulement de puissance, les variables de

perturbations ce sont les variables qu’on ne peut pas contrôler, dans les réseaux électriques, ce
sont des puissances demandées représentés par le vecteur d’état P.

⎡ ⎤
⎢ ⋮ ⎥
⎢ ⎥
⎢ ⎥
(P)=⎢ ⋯ ⎥ (1.13)
⎢ ⎥
⎢ ⎥
⎢ ⋮ ⎥
⎣ ⎦

5.2.Variables d’état (X)

Ce sont l’amplitude et l’angle de la tension au niveau de jeux de barres, ils

sont qualifiés d’état suite à leur valeur qui peut déterminer l’état de réseau .Elles sont
représentés par le vecteur d’état X.

⎡ ⎤
⎢ ⎥
⋮
⎢ ⎥
⎢ ⎥
=⎢ ⋯ ⎥ ,| | : tension de charge (1.14)
⎢| | ⎥
⎢| |⎥
⎢ ⋮ ⎥
⎣| |⎦

9
Chapitre I [ECOULEMENT DE PUISSANCE]

5.3.Variables de contrôle (U)

Les puissances générée, à savoir ou (actives ou réactives) sont des

variables de contrôles représentées par le vecteur de contrôle U.

⎡ ⎤
⎢ ⋮ ⎥
⎢ ⎥
⎢ … ⎥
=⎢ ⎥ (1.15)
⎢ ⎥
⎢ ⎥
⎢ ⋮ ⎥
⎣ ⎦

n: nombre de jeux de barres de tout le réseau.

6. Classification des jeux de barres

Les jeux de barres sont classés suivant les deux variables connus parmi les quatre
variables citées ultérieurement. Généralement en les classes sons trois types [ART 00], mais
l’actualité du réseau oblige un quatrième type.

6.1.Jeu de barres de charges (PQ)

Ce sont les jeux de barres de la demande clientèles, pour ce genre, les puissances
actives et réactives sont données, mais l’amplitude et l’argument de la tension sont à
déterminer [MER 14].ou bien c’est le JB où il n’y a pas de générateur ou un connait .

6.2.Jeu de barres de contrôle (PV)

Ce sont les jeux de barres de générateurs du réseau où la puissance active et le module

de la tension sont connus, puisqu’elles sont contrôlables ; par contre la puissance réactive et
l’argument de la tension sont inconnus. [MER 14]

6.3.Jeu de barres de référence (| | )

C’est un jeu de barres connecté a une source de tension considérée constante, et son
angle de phase est pris comme référence de calcul. Ses puissances actives et réactives, doivent
être calculées d’après le calcule de la puissance pour s’adapter aux exigences de jeu de barres
de contrôle. [MER 14]

10
Chapitre I [ECOULEMENT DE PUISSANCE]

6.4.Jeu de barre de contrôle (PQ)

Dans le calcule des variables inconnus (Q et ) d’un jeu de barre de contrôle on peut
trouver la valeur de l’énergie réactive Q dépasse les limites d’égalité ou d’inégalité de cette
dernier ; on est obligés donc de fixer la valeur de Q, ce jeu de barre reçoit un jeu de barre de
contrôle (P,Q). [MER 14]

7. Résoudre le problème de l’écoulement de puissance

7.1.La méthode itérative de Newton-Raphson

Le problème de l’écoulement de puissance peut être résoudre par la méthode de

Newton-Raphson. En réalité, parmi les nombreuses méthodes disponibles pour l’analyse de
l’écoulement de puissance, la méthode de Newton-Raphson est considérée comme la plus
raffinée et la plus importante. Elle n’est pas plus simple que la méthode de Gauss-Seidel,
mais elle est plus rapide (convergence) surtout pour les réseaux de tailles importantes.
Equation de puissance réelle [ART 00]

= , = 1,2,3, … , (1.16)

∗ ∗
= ( )∗ = ( (1.17)

Alors : = | | (1.18)

On définit :

= − (1.19)

= + (1.20)

= | | cos − + jsin( − ( − )) (1.21)

: La différence des phases entre les jeux de barres i et j

+ : La partie réelle et imaginaire de la matrice .

11
Chapitre I [ECOULEMENT DE PUISSANCE]

Etape 1 : on définit le type de jeu de barres i=1, un jeu de barre référence ( , ) connues
( , ) a déterminer [ART 00].

= ⋮ (1.22)

= ⋮ (1.23)

= (1.24)

Etape 2 : On définit les fonctions de ( ) et ( ):

= | | cos + sin (1.25)

= | | sin + cos (1.26)

Etape 3 : On forme l’équation ( )=0

( )− =0
( )− =0

 Ou , sont des puissances spécifies.

2( ) − 2
)
⎛ 3( − 3⎞
⎜ ( )− ⎟ ∆
( )=⎜ … ⎟ = (0) (1.27)
⎜ 2( ) − ∆
⎜ 2⎟⎟
3( ) − 3
⎝ ( )− ⎠

Etape 4 :

( ) ( )
11 = , 12 =
(1.28)
( ) ( )
21 = , 22 =

Alors :

12
Chapitre I [ECOULEMENT DE PUISSANCE]

∆ ∆
= ∗ (1.29)
∆ ∆

J=matrice de jacobienne.

7.2.La méthode de NR modifié (fast-decoupled):

En se basant sur les cas des différents réseaux étudiés à constater les valeurs des Sub-
matrices J12 et J21 que sont très petites. |MER 14]

0
=
0

L’explication de ces faibles valeurs revient au fait que la puissance active P

dépend nécessairement du θi; et non du |Vi | par contre Qi dépend essentiellement du |Vi | et
non θi. Généralement Bij≈ 0 et θij ≈ 0 d’où θi– θj = 0 alors sinθij = 0. La conductance est
presque nulle car les éléments connectés aux réseaux sont d’origine réactive.

8. Ecoulement de puissance continue

8.1.Facteur de charge [BEK 13]

L'outil analytique le plus accepté pour l’étude des phénomènes de

l’effondrement de tension est la théorie de la bifurcation, qui est une théorie mathématique
générale capable de classer les instabilités, d’étudier le comportement du système dans
le voisinage d’effondrement ou des points instables et de donner une information
quantitative sur les actions correctives pour éviter les conditions critiques.

Dans la théorie de la bifurcation, il est supposé que les équations du système dépendent d'un
ensemble de paramètres avec des variables d'état, comme suit:

0= ( , ) (1.30)

Avec :

: Equation de l’écoulement de puissance

: Variables dépendant

λ: Facteur de charge

13
Chapitre I [ECOULEMENT DE PUISSANCE]

Alors, du point de vue de stabilité des réseaux électriques, les propriétés de la

stabilité/instabilité sont évaluées en variant « lentement » les paramètres. Dans ce mémoire, le
paramètre utilisé pour étudier la proximité du système de l’effondrement de tension est le soi-
disant facteur de charge ( ∈ ℛ), ce qui modifie les puissances du générateur et la charge
comme suit:

= (1 + )( + ), = (1 + )( + ) (1.31)

Avec :

: Puissance active du générateur

: Puissance active de charge

: Fourniture de l’offre (Supply bids)

: Demande d’offre (demand bids)

Les puissances qui multiplient λ sont appelées les puissances de direction. Les
équations (1.25) différent du modèle généralement utilisé dans l’analyse du calcul continu de
l’écoulement de puissance (CPF), c.à.d.

=( + ), =( + ) (1.32)

Où le facteur de charge touche uniquement les variables puissances PS et PD. Dans

les diagrammes typiques de bifurcation, les tensions sont tracées en fonction de λ, c'est à dire
la mesure de la capacité de charge du système, obtenant ainsi les courbes V(p).

Les indices 0, 1 et 2 dénotent le cas de base, le premier point et le deuxième point de

directions des puissances respectivement.

8.2.Calcul Continu de l’écoulement de puissance [BEK 13]

Les techniques du calcul continu de l’écoulement de puissance sont largement

reconnues comme un outil précieux pour déterminer les courbes V(P) du système de
puissance et permettent d'estimer les conditions de chargement maximal et des
solutions «critiques». Malgré que les grands systèmes nécessitent un calcul numérique
exigeant, le CPF n'est pas affectée par des instabilités numériques. En fait, il est capable de
déterminer la partie stable et instable des courbes V(p) et peut fournir des informations
14
Chapitre I [ECOULEMENT DE PUISSANCE]

supplémentaires, telles que les facteurs de sensibilité de la solution actuelle à l'égard des
paramètres pertinents.

D'un point de vue mathématique, le CPF est une technique d'homotopie et permet
d'explorer la stabilité des équations du système électrique en faisant varier un paramètre du
système, qui, dans les études de stabilité statique et dynamique typique de tension, est
le paramètre de charge λ. De manière générale, le CPF consiste en un pas prédicteur réalisé
par le calcul du vecteur tangent et un pas correcteur qui peut être obtenu soit par
l'intermédiaire d'un paramétrage local ou à une intersection perpendiculaire.

a. Pas prédicteur

Au point d'équilibre générique p, la relation suivante s'applique:

, =0⇒ = | + =0 (1.33)

Et le vecteur tangent peut être approximé par :

∆
= ≈ (1.34)
∆

A partir de l’équation (1.27) et (1.28) on a :

=− | , ∆ = ∆ (1.35)

A ce point un pas de grandeur k de contrôle doit être choisi pour déterminer la

quantité∆ ∆ , avec une normalisation afin d’éviter de grands pas quand est grand.

∆ ≜ ∆ ≜ (1.36)

Où ‖.‖est la norme euclidienne et = ±1. Le signe de k détermine l’augmentation ou la

diminution de λ. La figure (1.3) présente une représentation graphique du pas prédicteur.

15
Chapitre I [ECOULEMENT DE PUISSANCE]

Fig.1.3 : calcul continu de l’écoulement de puissance: Pas prédicteur par les moyens de la
tangente

b. Pas correcteur

Pour le pas correcteur, l’ensemble des équations n+1 est résolu :

( , )=0 ( , )=0 (1.37)

Où la solution de doit être dans le collecteur de bifurcation et  est une équation

additionnelle pour garantir un ensemble non singulier au point de bifurcation. Pour le choix
de  il y’a deux options: l’intersection perpendiculaire et la paramétrisation locale.

Dans le cas de l’intersection perpendiculaire, dont la représentation est donnée par le logiciel
PSAT du MATLAB, l’expression de  devient (éq: 1.38)

Fig1.4 : Calcul continu de l’écoulement de puissance: pas correcteur obtenu par les moyens
d'intersection perpendiculaire.

16
Chapitre I [ECOULEMENT DE PUISSANCE]

∆ −( +∆ (1.38)
( , )= =0
∆ −( −∆ )

Tandis que pour la paramétrisation locale, soit, le paramètre soit la variable est forcée
d’être une valeur fixe.

( , )= − −∆ (1.39)

Ou

( , )= − −∆ (1.40)

Le choix pour que la variable soit fixe dépend du collecteur de bifurcation de , comme c’est
présenté à la figure (fig1.5)

Fig.1.5: Calcul continu de l’écoulement de puissance: pas correcteur obtenu par les moyens
de la paramétrisation locale.

9. Conclusion

Dans ce chapitre, afin d’analyser le réseau électrique, il nous faut déterminer

l’écoulement de puissance dont le calcul doit passer par la résolution des équation non
linéaires, où le recours aux méthodes numériques est inévitable. Quelques méthodes sont
étudiées et nous avons détaillé que la méthode de Newton-Raphson représente des avantages
énormes par rapport aux autres méthodes, c’est pour cette raison que nous allons
l’utiliser dans la partie du calcul de l’écoulement de puissance.

17
 Chapitre II [ENERGIES RENOUVELABLES]

1. Introduction

La plus grande part de la consommation énergétique mondiale est dédiée actuellement à

la production d’électricité.

Généralement, les sources à énergie fossile (classique) sont exploitées pour la production
de l’énergie électrique, mais elles présentent des inconvénients plus que des avantages
(pollution, l’effet de serre, trou d’Ozone, épuisement des gisements d’origine fossiles).
Aujourd’hui, l’exploitation d’énergie de forme propre a connu une grande préoccupation avec le
développement de la technologie dans le domaine énergétique. Cette forme d’énergie propre n’a
pas des effets nuisibles sur l’environnement et elle est inépuisable. Les sources des énergies
propres sont appelées les sources à énergie renouvelable.

2. Energie renouvelable

Une énergie renouvelable est une source d’énergie qui se renouvèle assez rapidement pour
être considérée comme inépuisable à l’échelle naturelle. Les énergies renouvelables sont
issues de phénomènes naturels réguliers, constants ou aléatoires provoqués par les astres,
principalement le Soleil (photovoltaïque), le vent (éolien), les marées (marémotrice) et la
chaleur du sous-sol (géothermique).

Aujourd'hui, nous assimilons souvent par abus de langage les énergies renouvelables
aux énergies propres.

2.1. Evolution de la production mondiale en énergies renouvelables.

Le recours aux énergies renouvelables en matière de la production de l’énergie électrique

a connu des aléas conformément aux contraintes politiques ; économiques et écologiques. La
production et la consommation sont deux faces d’une seule pièce de monnaie. La figure 2.1
illustre l’évolution de la consommation mondiale en énergie renouvelable dans l’année
2015.

18
 Chapitre II [ENERGIES RENOUVELABLES]

Fig.2.1 : Evolution de la consommation d’énergie mondiale.

Source : U.S. energy consumption by energy source, 2015
3. Différents types des énergies renouvelables
L’énergie renouvelable a plusieurs ressources, on cite quelque ressource connue comme :
3.1. Hydroélectrique
Les centrales hydroélectriques convertissent l’énergie de l’eau en mouvement en
énergie électrique. L’énergie provenant de la chute d’une masse d’eau est tout d’abord
transformée dans une turbine hydraulique en énergie mécanique. Cette turbine entraîne
un alternateur dans le lequel l’énergie mécanique est transformée en énergie électrique
(figure 2.2).

Avantage :

 Faible coût du KWh pour le gros hydraulique

 Souplesse d’utilisation
 Possibilité de « stocker » de l’électricité pour le turbinage
 L'hydraulique est une énergie renouvelable mise en valeur mais pourtant
rarement comptabilisée, elle fournit 13 % de la production et contribue également
à l'écrêtement des crues. [ZER 10]

Inconvénients :

 Contrainte géographique
 Aléas climatiques (sécheresse)

19
 Chapitre II [ENERGIES RENOUVELABLES]

 Nécessité de concilier les différents usages de l’eau

 Impacts sur l’environnement à la construction des centrales
 L'inondation de certains reliefs lors de la construction de barrages.

Fig.2.2 : Centrale hydroélectrique

3.2. Solaire

Le photovoltaïque, est un autre moyen de fabriquer de l’électricité avec l’énergie

solaire utilisant les rayonnements lumineux du soleil, qui sont directement transformés en un
courant électrique par des cellules à base de matériaux semi-conducteurs ayant des
propriétés de conversion lumière/électricité. Chaque cellule délivrant une faible tension, les
cellules sont assemblées en modules et en panneaux.

Avantages

 Acceptation du public
 Adapté aux pays chauds et aux sites isolés
 L'énergie solaire présente, outre tous les avantages des énergies renouvelables,
ses propres avantages, à savoir : énergie maîtrisable et adaptable aux situations de
toutes les régions.
 Frais de maintenance relativement faibles.

Inconvénients

 Coût encore élevé du KWh

 Puissance assez faible

20
 Chapitre II [ENERGIES RENOUVELABLES]

 Occupation des espaces agricoles et rurales

 Incapacité de fonctionnement du système photovoltaïque sans système de stockage.

Fig.2.3 : une centrale solaire

3.3. Eolienne
L'énergie éolienne est produite sous forme d'électricité par une éolienne en exploitant le
vent pour entrainer les hélices de la turbine éolienne et en convertissant l’énergie issue de la
turbine en énergie électrique. Cela se fait par l’intermédiaire d’un générateur (machine
asynchrone).
Avantage
 L'énergie éolienne est propre et peut être associée à des panneaux de photopiles et à
une batterie d'accumulateur afin d'optimiser la fourniture d'électricité.
 L’énergie éolienne est une énergie renouvelable ç.a.d. que contrairement aux énergies
fossiles, les générations futures pouvant toujours en bénéficier
 L’énergie éolienne n’est jamais être de risque comme l’énergie nucléaire
 L’exploitation de l’énergie éolienne n’est pas en procédé continu puisque les
éoliennes en fonctionnement peuvent facilement être arrêtées.

Désavantage

 L’impact visuel, qui reste néanmoins un thème subjectif.

21
 Chapitre II [ENERGIES RENOUVELABLES]

 La source d’énergie éolienne étant aléatoire donc la puissance électrique produite

n’est pas constante.
 Le coût de l’énergie éolienne par rapport aux sources d’énergie classique reste
encoure plus élevé.
3.4. Biomasse
L’énergie de la biomasse est une bioénergie qui utilise le combustible de la matière
végétale ou organique, appelée biomasse. En brûlant la source d’énergie, la centrale
thermique dédiée va alors produire de l’électricité afin d’alimenter le réseau électrique.
L’énergie de la biomasse peut aussi être utilisée pour le chauffage. Faisons un point sur les
avantages et les inconvénients de l’énergie de la biomasse

Avantages

 L’énergie de la biomasse est inépuisable. Le déchet étant produit par chaque espèce
vivante.
 L’énergie de la biomasse est une énergie qui peut préserver les réserves fossiles.
 L’énergie de la biomasse est un type d’énergie fait appel à des ressources peu
coûteuses, qui sont disponibles en grandes quantités.
 En développant l’énergie de la biomasse, il serait également possible de réduire la
production de méthane, qui est un autre gaz à effet de serre.

Désavantage

 L’utilisation du biogaz généré par la fermentation des matières de la biomasse, pour

produire de l’électricité permettrait en plus d’utiliser d’autre déchet.
3.5. Géothermie
L’énergie de la terre, la chaleur naturellement présente dans le sous-sol de planète
représente une formidable source d’énergie. Plus on creuse profondément, plus on atteint des
températures élevées.
La géothermie, utilise cette chaleur pour le chauffage et pour production l’électricité.
Avantages
 Ressource énergétique locale qui ne nécessite pas de transport
 Ressource énergétique renouvelable préservant l’environnement
 Capacités de production importantes comparées aux autres énergies renouvelables
 Énergie de base, indépendant des conditions climatiques

22
 Chapitre II [ENERGIES RENOUVELABLES]

Désavantages
 Ressource naturelle du sous-sol dont l’extraction nécessite la réalisation de forages.
 Délai important dans la mise en valeur d’une ressource géothermique.
 Coût d’investissement important comparé aux autres énergies renouvelables (travaux
d'exploration, réalisation de forages).
3.6. Energie marine
L'énergie marine ou marémotrice, est une source d’énergie renouvelable très importante
qui exploite le mouvement des marées. Ce phénomène est dû aux interactions
gravitationnelles entre la Terre et des autres astres comme la Lune et le Soleil (produite par la
gravité lunaire plutôt que par le rayonnement solaire).
Avantages
 Une des avantages de cette énergie, c’est qu’elle ne représente qu’une petite partie des
besoins de l’humanité. Et au contraire des autres énergies renouvelables, celle-ci est
totalement prédictible, c’est-à-dire qu’on peut la prévoir facilement.
 Aussi c’est une énergie à la fois propre, naturelle et renouvelable.
 Les marées qui arrivent des mers et océans constituent une source d'énergie durable et qui
n'est pas soumise aux aléas climatiques contrairement à l'énergie solaire ou éolienne. Les
usines marémotrices peuvent donc être des solutions alternatives contre le danger de
radioactivité des centrales nucléaires.
Désavantages
 Dans le monde il y a peu de pays que utilisent cette énergie parce qu’il a besoin d’être
équipe de centrales de production marémotrices car ils nécessitent des conditions
hydrodynamiques particulière. La composition des sols doit permettre à l'usine de rester
en équilibre.
 Aussi l’énergie marémotrice a besoin des aménagements importants, qui perturbent les
équilibres écologiques des sites d’implantation.

23
 Chapitre II [ENERGIES RENOUVELABLES]

Fig.2.4 : Différents types des énergies renouvelables

4. Principaux sources à énergie renouvelables

4.1. Energie éolienne

Parmi toutes les énergies renouvelables, à part l'énergie du bois, c'est l'énergie du vent qui a
été exploitée en premier lieu par l'homme. Depuis l'antiquité, elle fut utilisée pour la propulsion
des navires et ensuite les moulins à blé et les constructions permettant le pompage d'eau.

En Europe, les premiers moulins à vent ont fait leur apparition au début du Moyen Age.
Utilisés tout d'abord pour moudre le grain, d'où leur nom de " moulins ", ils furent
aussi utilisés aux Pays-Bas pour assécher des lacs ou des terrains inondés. Dès le XIV siècle, les
moulins à vent sont visibles partout en Europe et deviennent la principale source
d'énergie [HLA 84].

A l'arrivée de la machine à vapeur, les moulins à vent commencent leur disparition

progressive. L'arrivée de l'électricité donne l'idée à Poul La Cour en 1891 d'associer à une
turbine éolienne une génératrice. Ainsi, l'énergie en provenance du vent a pût être
«redécouverte » et de nouveau utilisé (dans les années 40 au Danemark 1300 éoliennes). Au
début du siècle dernier, les aérogénérateurs ont fait une apparition massive (6 millions de
pièces fabriquées) aux Etats-Unis où ils étaient le seul moyen d'obtenir de l'énergie électrique
dans les campagnes isolées. Dans les années 60, fonctionnait dans le monde environ 1 million

24
 Chapitre II [ENERGIES RENOUVELABLES]

d'aérogénérateurs. La crise pétrolière de 1973 a relancé de nouveau la recherche et les

réalisations éoliennes dans le monde [RAI 02].

Cette énergie est produite par la force exercée par le vent sur les pales d'une hélice. Cette
hélice est montée sur un arbre qui est relié à des systèmes mécaniques qui servent à moudre le
grain ou à pomper de l'eau, soit à un générateur qui transforme l'énergie mécanique en énergie
électrique.

Les éoliennes sont une forme très ancienne d’exploitation du vent. Une importante
évolution 1970-2000.Une éolienne qui est un dispositif pour extraire de l'énergie cinétique
du vent en la transformant en énergie mécanique ou en énergie électrique. En enlevant
une certaine partie de son énergie cinétique le vent doit ralentir mais seulement cette masse d'air
qui traverse le disque du rotor sera affectée [PER 11].

Fig.2.5 : Une ferme d’éoliens.

4.1.1. Evolution de l’énergie éolienne dans le monde

En 2017, l’énergie éolienne a couvert plus de 3% de la consommation d’électricité
mondiale totale grâce à une puissance de 540’000 MW. [source : association mondiale de
l’éolien: GWEC]

25
 Chapitre II [ENERGIES RENOUVELABLES]

Fig.2.6 : Capacité mondiale installé d’éolienne 2001-2017.

[Source : association mondiale de l’éolien: GWEC].

4.1.2. Types d’éolienne

Les éoliennes se divisent en deux grandes familles : celles à axe vertical et celles à axe
Horizontal

4.1.2.1. Eoliennes à axe horizontal

Les éoliennes les plus répandues sont à axe horizontal, on peut distinguer :
- les moulins à vent classique,
- les éoliennes lentes,
- les éoliennes rapides.
4.1.2.2. Eoliennes à axe verticale

Les éoliennes à axe verticale peut classer on deux catégories

- Les éoliennes lentes

- Les éoliennes rapides

Fig.2.7 : éolienne type horizontale Fig.2.8 : éolienne type verticale

26
 Chapitre II [ENERGIES RENOUVELABLES]

4.2. Energie solaire

Le soleil, bien que distant de plus de 150 millions de kilomètres de nous, demeure
notre plus grande source d'énergie même si elle est intermittente avec l'alternance jour / nuit...

La source d’énergie renouvelable la plus abondante est sans aucun doute l’énergie
solaire. Le soleil envoie sur Terre une puissance de 90.000 térawatts (TW), à comparer aux 14
TW consommés par la population mondiale, et va continuer à le faire pendant quelques milliards
d’années. Cependant, en raison de son caractère intermittent et de sa faible densité, il est
capital de pouvoir stocker et concentrer cette énergie solaire sous forme de carburant. [PER
11]

4.2.1. Différents types d’énergie solaire

On peut classer l’énergie solaire on deux catégories :

4.2.1.1. L’énergie solaire thermique

Le principe de l’énergie thermique consiste à transformer la chaleur du rayonnement solaire en

énergie thermique grâce à un fluide qui circule dans des panneaux exposés au soleil.

Il y a plusieurs manière pour capturer et exploiter ce chaleur du rayonnement solaire on cite :

a. Capteur parabolique

Les capteurs paraboliques fonctionnent d'une manière autonome. Ils sont constitués d’une
grande parabole de révolution réfléchissante et d’un moteur « Stirling » au foyer (on dite le foyer
ou bien le récepteur) de la parabole. Le tout pivote sur deux axes pour suivre le déplacement du
soleil afin de concentrer son rayonnement sur le foyer de la parabole réfléchissante. Le rapport
de concentration est généralement d’environ 4000 et la température obtenue entre 500 et 1000°C.

Ainsi, les rayons solaires chauffés le fluide comprimé afin de générer de l'électricité.

27
 Chapitre II [ENERGIES RENOUVELABLES]

Fig.2.9 : le fonctionnement d’un capteur parabolique.

b. Capteur cylindro-parabolique

Ce type de centrale se compose d’alignements parallèles de longs miroirs hémicylindriques,

qui tournent autour d’un axe horizontal pour suivre la course du soleil. Les rayons solaires
sont concentrés sur un tube horizontal dans lequel circule un fluide caloporteur
(généralement une huile synthétique). Les tuyaux étant noirs, ils absorbent toute la chaleur
du soleil et permettent à la température du fluide de monter jusqu’à 500° C.

La chaleur ainsi récupérée produit de la vapeur via un échangeur, cette vapeur a son tour
va actionner des turbines pour produire de l’électricité.

Fig.2.10 : Capteur cylindro parabolique

4.2.1.2. Energie solaire photovoltaïque

Contrairement à l’énergie solaire thermique qui utilise le soleil comme source de chaleur,
l’énergie solaire photovoltaïque utilise le soleil comme source de lumière en transformant en

28
 Chapitre II [ENERGIES RENOUVELABLES]

électricité l’énergie des photons arrivant à la surface de la terre. La lumière solaire (photons)
transmet son énergie aux électrons continus dans un semi-conducteur (qui constitue une cellule
photovoltaïque).

Fig.2.11 : Cellule PV Fig.2.12 : Panneau PV Fig.2.13 : Champ PV

5. Conclusion

Nous avons exposé dans ce chapitre les principaux types d’énergies renouvelables, plus
particulièrement l’énergie éolienne et l’énergie solaire, en donnant un aperçu sur leur
exploitation dans le monde et leur importance énergétique sur le plan économique,
environnemental et sociale. D’où l’importance de l’intégration des énergies renouvelables dans
les réseaux électriques, cela constitue le thème du chapitre suivant qui va exposer l’impact de ces
énergies sur le réseau électrique.

29
 Chapitre III [INTEGRATION DES SOURCES RENOUVELABLES]

1. Introduction

Les énergies renouvelables sont nécessaires pour réduire la pollution des gaz toxique émis à
partir des centrales classiques et préserver les réserves fossiles. C’est dans ce but que réside
l’importance d’intégrer les sources renouvelables dans le réseau, mais jusqu’à maintenant cette
intégration a des impacts sur la stabilité du réseau électrique et sur le coût d’investissement des
centrales à énergie renouvelable.

Dans ce chapitre, nous traitons le besoin de l’intégration des énergies renouvelables dans le
réseau, ainsi que les différents impacts de leur insertion dans le réseau électrique.

2. Le besoin de l’intégration dans le réseau électrique

Une des conséquences de l’ouverture du marché de l’électricité est le développement

d’une production décentralisée, de sources d’énergie renouvelable.

L’intégration dans les réseaux électriques des sources à énergie renouvelables, et plus
particulièrement celles soumises aux aléas du climat telles que les énergies éolienne et solaire, et
d’une manière plus générale de la production décentralisée va nécessiter d’importants
aménagements de ces réseaux, ainsi que la mise en œuvre de nouveaux équipements et de
nouvelles méthodes de gestion.

Le défi est de maintenir la fiabilité et la qualité de l’alimentation des particuliers et des

entreprises en énergie électrique malgré la libéralisation du marché de l’électricité et l’utilisation
croissante de sources à énergies renouvelables [HEI 05].

L’intégration nécessite de garder une bonne qualité d’énergie électrique sous différentes
conditions de fonctionnement du système énergétique.

3. Qualité de la tension et perturbations du réseau électrique

3.1. Qualité de la tension

La qualité de tension est le concept d'efficacité de classer les équipements sensibles d'une
manière qui convient à l'opération de l'équipement.
Pour rappel, la tension possède quatre caractéristiques principales :

 Amplitude.
 Fréquence.
 Forme d'onde.

30
 Chapitre III [INTEGRATION DES SOURCES RENOUVELABLES]

 Symétrie.
3.2. Différents perturbations dans le réseau électrique

Le réseau électrique peut être soumis à de multiples perturbations qui peuvent être mesurés
grâce à l’utilisation de centrales de mesure, peuvent endommager les machines et fragiliser
toutes les installations électriques.

Il existe différents types de perturbation électrique, on cite :

3.2.1. Creux de tension et coupures

3.2.2. Harmoniques
a) Inters harmoniques
b) Infra-harmoniques
c) Conséquences des harmoniques
d) Les principales sources d’harmoniques
3.2.3. Surtension
3.2.4. Fluctuation de la tension (Flicker)
3.2.5. Déséquilibre
4. Impact de la production d’énergie renouvelable sur le réseau électrique

Dans cette partie, nous allons citer les différents impacts de l’insertion des énergies
renouvelables dans le réseau électrique.

Les études ont montré qu’un taux de pénétration croissant des énergies d’origine
renouvelable n'est pas sans impacts prévisibles sur l'exploitation future des réseaux.

En particulier, le plan de tension peut être gravement modifié par la présence de GED
(Génération d’Energie Dispersée(Distribuée)), au point que la tension risque d’avoir des
fluctuations en certains nœuds du réseau alors qu'elle est maintenue à une valeur normale au
poste source. Le plan de protection risque également d'être affecté par un fort taux de
pénétration des GED du fait de la puissance de court circuit qu'elles apportent en aval
des protections, et de l'inversion possible des flux de puissance active sur certaines lignes,
ainsi que de la diminution du temps d'élimination critique de défauts. [SAD 14]

Une partie de ces GED à, de plus, des sources d'énergie primaire intermittentes
(éolien, solaire) qui ne permettent pas de prévoir aisément la production disponible à court
terme. Elles ne peuvent donc pas garantir une puissance de sortie et proposer toute la puissance

31
 Chapitre III [INTEGRATION DES SOURCES RENOUVELABLES]

disponible sur le marché. D'autre part, ces nombreuses sources sont trop petites pour être
observables et dispatchables par les gestionnaires de réseaux et ne participent donc pas,
aujourd'hui, aux services système. Cela peut poser des problèmes en cas de fort taux de
pénétration si les moyens de réglage classiques de la distribution deviennent inaptes à assurer la
tenue en tension. Cela peut contraindre par exemple les gestionnaires de réseaux à engager
des investissements non prévus initialement. [ARB 06]

Une partie de ces GED produisant par construction du courant continu (pile a
combustible, panneau solaire) doit être raccordée au réseau par l'intermédiaire d'interfaces
d'électronique de puissance, injectant ainsi des harmoniques nuisibles à la qualité de la tension
délivrée. [SAD 14]

Enfin, la présence de GED en aval d'un poste source dont le transformateur est équipé
d'un régleur en charge régulé par compoundage perturbe son fonctionnement basé sur la
mesure du courant absorbé. [SAD 14]

En effet, la production de puissance par les GED réduit le courant traversant le

transformateur, provoquant une action du régleur en charge et diminuant ainsi la tension au
niveau du poste source.

Fig.3.1 : Flux d'énergie sur un réseau en présence de GED

Nous trouvons dans la littérature plusieurs études qui ont été faites afin d’identifier
et d’évaluer les impacts de la production décentralisée sur le réseau, puis, de proposer des
solutions adaptées.

L’opérateur du réseau se préoccupe de ce type d’étude lors de la phase de planification ou de

l’approbation d’un projet de raccordement des GED. Les principaux impacts à rappeler
sont les suivants : [CAM 03]

32
 Chapitre III [INTEGRATION DES SOURCES RENOUVELABLES]

4.1. Impact sur la qualité de service

Les GED de type asynchrone consomment de la puissance réactive afin de magnétiser leur
circuit magnétique. Lors de la connexion au réseau, elles appellent un courant fort, ce qui
contribue au creux de tension (en profondeur).

D’ailleurs, la présence d’interfaces d’électronique de puissance peut faire augmenter le taux

des harmoniques qui nuisent gravement à la qualité de service fournie. [PER 11]

4.2. Impact sur la continuité de service

Pour la même raison concernant la caractéristique d’intermittence, l’indisponibilité des GED

lors que le système les sollicite peut occasionner la rupture d’électricité par manque de
puissance.

4.3. Impacts sur le plan de protection

La connexion de nouveaux générateurs au réseau modifiera également les courants de courts

circuits et donc le réglage des seuils des protections.

D’ailleurs, ce changement peut provoquer une mauvaise sélectivité, des déclenchements

intempestifs ou l’aveuglement du système de protection. [GIP 07]

4.3.1. Courant de court-circuit

En matière de protection, il faut tenir compte de ce que les générateurs contribuent à une
augmentation du courant de court-circuit dans le réseau. Dans le cas de générateurs, directement
raccordés au réseau, cette contribution est de l’ordre de 4 à 8 fois le courant nominal de la
machine.

Dans le cas d’un raccordement au travers d’un convertisseur électronique cette

contribution peut être ramenée au courant nominal par contrôle rapide du convertisseur.
Signalons toutefois qu’avec ce dernier système en cas de fonctionnement sur charge isolée
du réseau suite à un incident, les protections par surintensité ne fonctionneraient pas pour
éliminer un court-circuit dans le réseau. [SAD 14]

Ce problème pourrait systématiquement se rencontrer avec des systèmes photovoltaïques

ou à piles à combustibles sur réseau isolé. Des précautions particulières sont à prévoir.

33
 Chapitre III [INTEGRATION DES SOURCES RENOUVELABLES]

4.3.2. Influence sur la sélectivité des protections

La sélectivité des protections peut être affectée par l’introduction d’unités de production
dans le réseau de transport dont la structure est de type maillé. Pour illustrer cette influence
considérons le schéma de la figure (3.2) ci-dessous d’alimentation par une sous-station du réseau
d’une charge par la liaison 2, et la connexion à la sous-station par la liaison 1 d’une unité de
production. Chacune de ces deux liaisons est protégée par un disjoncteur contre les surintensités
comme c'est l'usage.

Dans cet exemple extrêmement simple la protection de la ligne 1 déconnecte

inutilement l’unité de production lors d’un défaut sur la ligne 2, car le courant du générateur lors
de ce défaut est supérieur au seuil de protection en cas de défaut sur la ligne 2. Les seuils de
protection doivent donc être revus pour que le seule de la ligne en défaut soit déconnectée.

Fig.3.2 : Un défaut sur la ligne 2 provoque la coupure inutile de la ligne 1

La sélectivité de la protection est modifiée car les unités de production peuvent modifier
l’amplitude ; la durée et la direction des courants de défaut

Pour toute implantation d’une unité de production dans le réseau, il faut impérativement
vérifier la sélectivité des protections et le cas échéant les ajuster.

4.4. Impact sur la tenue de fréquence et de tension

4.4.1. Impact sur la tenue de fréquence

La fréquence doit impérativement être maintenue dans les limites contractuelles.

Les principaux problèmes apparaissant en cas d'excursion non maîtrisée de la fréquence
sont les suivants:

 Rupture de synchronisme des alternateurs pouvant entraîner la perte de groupes en

cascade.
 En cas de baisse de fréquence :

34
 Chapitre III [INTEGRATION DES SOURCES RENOUVELABLES]

Saturation des circuits magnétiques de transformateurs ou compensateurs à bobine d'énergie

réactive entraînant un fonctionnement dégradé et des pertes fer. [WAG 10]

4.4.2. Impact sur la tenue de tension

Nous savons que les nœuds proches d’une source ont une tension plus élevée, la
tension diminuant au fur et à mesure que l’on s’en éloigne. La connexion de sources de GED à
proximité des charges modifiera donc naturellement le profil de tension sur le réseau [GOU 12]

5. Intégration en Algérie

L’Algérie en particulier et les pays du Maghreb ont un potentiel solaire élevé. Les taux
d’irradiation solaire effectués par satellites par l’Agence Spatiale Allemande DLR, montrent
des niveaux d’ensoleillement exceptionnels de l’ordre de 1200 kWh/m²/an dans le Nord du
Grand Sahara.

Par contre, les meilleurs taux d’irradiation solaire en Europe sont de l’ordre de 800
kWh/m²/an limités à la partie sud de l’Europe. Suite à une évaluation par satellites, l’Agence
Spatiale Allemande (ASA) a conclu, que l’Algérie représente le potentiel solaire le plus
important de tout le bassin méditerranéen, soit: 169.000 TWh/an pour le solaire thermique,
13,9TWh/an pour le solaire photovoltaïque et 35 TWh/an pour l’éolien [PER 11]

Dans cette partie, nous décrivons le potentiel Algérien en matière d’énergie

renouvelable surtout solaire et éolien, ainsi que leurs intégrations dans une vision de
développement a travers plusieurs scénarios énergétiques, en mettant l’Algérie comme un noyau
de production d’énergie renouvelable par rapport à son environnement surtout européen, dont
les études de la banque mondiale indique qui il sera obligé dans les années à venir d’ importer
son énergie du grand Sahara. D’ou vient l’initiative allemande pour exporter l’énergie
solaire du grand Sahara à l’Europe. [SAD 14]

35
 Chapitre III [INTEGRATION DES SOURCES RENOUVELABLES]

5.1. L’énergie éolienne en l’Algérie

L'énergie cinétique du vent constitue une ressource énergétique inépuisable mais

fluctuante. Pour toute installation d'une éolienne dans un endroit, il faut d'abord une étude des
caractéristiques du site considéré. Dans cette partie, on s'intéresse au potentiel éolien en
Algérie.

5.1.1. Les ressources en énergie éolienne

L'estimation des ressources en énergie éolienne présente une difficulté majeure. La

quantité d'énergie disponible varie avec la saison et l'heure du jour.

La quantité totale d'énergie éolienne produite est en fonction de deux facteurs, le premier
facteur est celui de la disponibilité de l’énergie éolienne produite en fonction du site
d’emplacement du parc éolien sur un territoire, c'est-à-dire la quantité d'énergie qui peut être
réellement produite par la mise en œuvre à l'échelle d'une région (énergie éolienne disponible), le
deuxième facteur est défini par les caractéristiques techniques des éoliennes placées dans le parc
éolien, c’est au niveau d’éolienne que l’énergie récupérable dépend des caractéristiques du
système de conversion utilisé. Chacun de ces facteurs présente des avantages et des
inconvénients suivant son application. L’Algérie est un pays qui est placé dans un lieu trop
riche avec des sources des énergies développées. Le vent est très fort dans quelques régions
(Adrar, Tiaret et Telemcen...ect).

5.1.2. Vents en Algérie

Les vents « GISEMENT EOLIEN » de l'Algérie commencent à être timidement étudiés.

L'évaluation globale du gisement éolien dans notre pays se fait en première phase à partir des
données des stations météorologiques classiques, les vents sont mesurés à un ou deux mètres du
sol. Lorsque les vents dans une zone ou une région sont définis comme source d'énergie, des
études approfondies doivent être faite selon les aérogénérateurs que l'on désire installer [PER
11].

Toutefois, la vitesse du vent subit des variations en fonction des saisons qu'on ne doit pas
négliger, en particulier, lorsqu'il s'agit d'installer des systèmes de conversion de l'énergie
éolienne.

Néanmoins, nous remarquons que la région de TIARET fait l'exception avec une
vitesse plus faible en été qu'en hiver. D'un autre côté, des régions telles que Biskra, Adrar et

36
 Chapitre III [INTEGRATION DES SOURCES RENOUVELABLES]

Annaba sont caractérisés par une vitesse relativement constante tout au long de l'année. [SAD
14]

Fig.3.3 : Atlas de la vitesse moyenne du vent de l'Algérie estimée à 10 m du sol.

5.1.3. La ferme éolienne d’Adrar en Algérie

La Wilaya d’ADRAR s’étend sur une superficie de 427.968 km², soit près d’un
cinquième (1/5) du territoire national. L’espace utile demeurant toutefois réduit, la wilaya
d’ADRAR située à l’extrême sud ouest du pays, compte onze (11) daïra et vingt huit (28)
communes. Sa population estimée à 336 046 habitants dont 76% de population rurale, soit 254
991habitants, repartie à travers 294 Ksour (localités) éparpillés au niveau de quatre régions

• Le Touat.
• Le Tidikelt.
• Le Gourara.
• Le Tanezrouft.

La Superficie Agricole Totale (SAT) est de 337 650 HR, la Superficie Agricole utilisée
(SAU) couvre une superficie de 31 657 HR dont 24 675 HR de Superficie Agricole Utilisée
Irriguée (SAUI) repartie comme suit : [SAD 14]

• 19 303 HR irriguée par système gravitaire.

• 2 331 HR irriguée par système localisée.
37
 Chapitre III [INTEGRATION DES SOURCES RENOUVELABLES]

• 3 041 HR irriguée sous pivots (irrigation par aspersion).

Fig.3.4 : La ferme éolienne de Kabertène à 75 km au nord est d’Adrar

5.2. L’énergie solaire en Algérie

L'Algérie, de part sa situation géographique, bénéficie des conditions favorables à

l'utilisation des énergies renouvelables, en particulier l'énergie solaire, l'une des possibilités de
l'exploitation de l'énergie solaire est sa transformation directe en une énergie électrique,
au moyen des convertisseurs photovoltaïques. Plusieurs projets et programmes solaires ont
été initiés et développé notamment par le centre de développement des énergies
renouvelables (CDER).

L'intensité de l'énergie solaire reçue à la terre est variable et très souvent les besoins en
électricité ne correspondent pas aux heures d'ensoleillement. Conséquemment, il est nécessaire
d'équiper le système par des batteries d'accumulateurs qui permettent de stocker l'électricité et
de la restituer en temps voulu.

Généralement, les travaux effectués sont axés et orientés vers l'optimisation de la puissance
du système photovoltaïque:

 Par l'orientation du panneau photovoltaïque.

 Par la recherche du point de la puissance maximale [MPPT].
5.2.1. Atlas solaire de l'Algérie
5.2.1.1. Projet « Atlas 1 »

Le plan Algérien en matière d’énergies renouvelables passe la vitesse supérieure. Pour

preuve, la création d’un ministère consacré aux énergies renouvelables constitue la dernière
initiative de développement de la production d’énergie solaire et éolienne menée par le
gouvernement.

38
 Chapitre III [INTEGRATION DES SOURCES RENOUVELABLES]

Nul ne peut nier aussi que le gisement d’utilisation des énergies renouvelables est
potentiellement important dans notre pays. Selon les professionnels, pour réussir un tel défi, une
stratégie bien définie, voir bien pensée doit être mise en place. Plusieurs projets sont lancés dans
le domaine des énergies renouvelables et nous ne pouvons pas parler de ce secteur névralgique
sans évoquer le projet « Atlas 1 ». Les dernières informations à ce sujet parlent du lancement
d’un appel d’offres national. Le projet exclusivement Algérien dénommé « Atlas1 » est d’une
capacité globale de 4050 mégawatts réparti en trois lots de 1350 mégawatts chacun en énergie
d’origine solaire de type photovoltaïque. [ACT 18]

5.2.2. Potentiel solaire

L’Algérie possède un gisement solaire parmi les plus élevé dans le monde ,la
durée moyenne d’ensoleillement dans le Sahara algérien est de 3500 heures ,ce potentiel
peut constituer un facteur important de développement durable dans cette région, s’il est exploité
de manière économique ,le tableau suivant indique le taux d’ensoleillement pour chaque
région de l’Algérie De par sa situation géographique, l’Algérie dispose d’un des
gisements solaires les plus importants du monde. La durée d’insolation sur la quasi totalité du
territoire national dépasse les 2000 heures annuellement et atteint les 3900 heures (hauts plateaux
et sahara).

De par sa situation géographique, l'Algérie dispose d'un gisement solaire énorme surtout
dans le Sud du pays (Tableau 3.1). La durée d'insolation sur la quasi totalité du territoire national
dépasse les 2000 heures annuellement et peut atteindre les 3900 heures

Tableau (3.1) : Potentiel solaire dans l’Algérie.

Région Région côtières Hauts plateaux Sahara

Superficie % 4% 10% 86%
Durée moyenne d’ensoleillement 2650 3000 3500
(Heures/an)
Energie moyenne reçue 1700 1900 2650
(KWh/m²/an)

39
 Chapitre III [INTEGRATION DES SOURCES RENOUVELABLES]

6. Conclusion

Dans ce chapitre, nous avons montré que l’intégration des énergies renouvelables est
nécessaire pour réduire la pollution des gaz toxiques, mais elles ont différents impacts
d’insertion dans le réseau électrique.

Ainsi, nous allons voir que l’Algérie, de part sa situation géographique, elle est bénéficié par
des conditions favorables pour l’utilisation des énergies renouvelables, comme l’énergie éolienne
et l’énergie solaire.

40
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE]

1. Introduction

La stabilité occupe une grande part dans les travaux de recherche sur l’analyse de réseau
électrique.

Les perturbations pouvant entraîner de sérieux dégâts sur le fonctionnement normal du

réseau électrique, notamment les générateurs et les lignes de transmission. Ces perturbations à
cause de plusieurs défauts soit externe (comme la foudre) ou bien interne (comme la variation
graduelle ou brusque des charges), ces derniers affectent le bon fonctionnement du système
énergétique en aboutissant à une dégradation de profile de la tension (entrainent un effondrement
de la tension) ou déséquilibre de production-consommation de l’énergie électrique (peut
conduire à la perte de synchronisme) par conséquent la perte de la stabilité du réseau.

La stabilité d’un réseau électrique est la propriété qui lui permet de rester dans un état
d’équilibre, pour des conditions de fonctionnement normales, et de retrouver un état d’équilibre
acceptable, suite à une perturbation [BHL 91].

Dans ce chapitre on étude et analyse la stabilité statique des réseaux de 14 JB et de 30

JB de norme IEEE par le « CPF » de PSAT qui est un logiciel fonctionne sous MATLAB, à
travers ces études, on confirme l’emplacement optimale pour l’intégration d’une source à énergie
renouvelable pour améliorer la marge de stabilité de ces réseaux, puis on analyse la stabilité
dynamique le réseau de 30JB.

On peut classer les types de la stabilité comme suit :

41
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE
[ANLYSE STABILITE]

Fig4.1 : Classification des types de stabilité des réseaux électriques.

2. La stabilité statique de tension

L’analyse statique de la stabilité implique seulement la solution des équations algébriques

du problème d’écoulement de puissance. De ce fait, du point et du vue calcul, elle est beaucoup
plus efficace que l’analyse dynamique. L’analyse statique est idéale pour la majeure partie des
études dans lesquelles des limites de stabilité de tension pour des pré-contingences
pré contingences et des post
post-
contingences doivent être déterminées.

Généralement, l’analyse de la stabilité statique de tension est utilisée afin de vérifier si un

point de fonctionnement est stable ou non, évaluer la marge de stabilité d’un point de
fonctionnement et identifier le point d’effondrement de
d tension.

Plusieurs outils de simulation on été utilisée pour l’analyse de la stabilité tels que
(MATLAB, Power World, PSAT,…etc). Ces logiciels permettent de nous renseigner sur les
différentes grandeurs d’un réseau électrique et à simuler les différents cas pour avoir les
paramètres du réseau. Ces derniers utilisent les mêmes concepts. Dans le présent travail, on
choisi PSAT comme outil de simulation.

42
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE]

3.1. C’est quoi le PSAT ?

C’est un logiciel développé par Dr. Federico Milano en 2001. Le PSAT est une boîte à
outils qui s’exécute sous environnement MATLAB. Ce logiciel à développer pour l’analyse
statique et dynamique et le contrôle des réseaux électriques. Il utilise l’écoulement de puissance,
l’écoulement de puissance optimal, la stabilité de tension… . L’écoulement de puissance à
utilisée la méthode de Newton-Raphson (NR)
Le PSAT à plusieurs versions apparut, nous exploitons dans notre présent mémoire la
version 2.1.3 datée en February 14,2009 (sous MATLAB 2009b) dont la page principale est
illustrée ci-dessous :

Fig4.2 : page d’accueil de PSAT

3.2. Utilisation de PSAT

Le choix du PSAT est justifié par la richesse de sa bibliothèque, pour l’analyser, contrôler,
régler, simuler… des systèmes électriques. Comme suit :

43
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE]

Fig4.3 : Bibliothèque du Simulink.

Toutes les opérations de PSAT sont réparties en deux types d’analyses :
 La première analyse et le calcul de l’écoulement de puissance. Cette application
s’effectue dans la page principale de PSAT comme la figure 4.2
 La seconde opération consiste à construire le réseau à étudier en utilisant la bibliothèque
de Simulink.
On utilise le PSAT :
 Ouvrir MATLAB
 Ecrit « psat » dans la « command window » de MATLAB
 Ouvrir un nouveau fichier
 Créer du réseau électrique à étudier
 Lancer le modèle de la bibliothèque de simulation
 Introduire les données numériques des éléments (générateur, lignes, charges,
transformateurs,….).
 Enregistrer le modèle.
 Charger le modèle :
 Click deux fois sur « data file » de page d’accueil de PSAT.
 Choisir « .mdl » de « Filters » puis « tests » de ‘Folders in currentpath’ et a
partir de « Files in currentpath » on charger le modèle d’étude et click sur
« load » pour la terminaison de charger.
Pour l’étude à l’état initial :
44
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE]

 Exécuter l’écoulement de puissance (Power Flow).

Puis pour l’étude à l’état statique :

 Exécuter l’écoulement de puissance continue « CPF »

4. Ecoulement de puissance continue (continuation power flow « CPF »)

La définition et les types de l’écoulement de puissance continue sont représentés dans le

chapitre I.

5. Analyse de la stabilité statique de réseau électrique

5.1. Analyse de RE de 14 JB

Bus 13

Bus 14
Bus 10
Bus 12

Bus 09
Bus 11

Bus 07

Bus 06
Bus 08
Bus 04

Bus 05
Bus 01

Bus 02

Bus 03

Fig. 4.4: Le schéma de topologie du modèle de 14-JB

45
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE]

Remarque : la fréquence de ce réseau est 60Hz et la simulation aussi avec 60Hz

On analyse ce réseau à l’état initial et on obtient les résultats suivants de « static report » :

Tableau 4.1.Caractéristiques du réseau 14 JB

Nombre de Nombre Nombre de ligne Nombre de Nombre de

jeux de barres de charges de transmission générateurs Transformateur
14 11 16 5 4

Tableau4.2. les résultats de l’écoulement de puissance à l’état initial du modèle de 14JB

Jeu de La tension Angle Pgen Qgen Pcharge Qcharge

Barre (JB) de JB [pu] [rad] [pu] [pu] [pu] [pu]
01 1.06 0 3.5206 -0.27943 0 0
02 1.045 -0.13564 0.4 0.9532 0.3038 0.1778
03 1.01 -0.33196 0 0.59933 1.3188 0.266
04 0.9975 -0.26398 0 0 0.6692 0.056
05 1.0024 -0.22713 0 0 0.1064 0.022 [1]4
06 1.07 -0.3962 0 0.45686 0.1568 0.105
07 1.0349 -0.36155 0 0 0 0
08 1.09 -0.36155 0 0.34081 0 0
09 1.0117 -0.41098 0 0 0.413 0.2324
10 1.0112 -0.41542 0 0 0.126 0.0812
11 1.0351 -0.40862 0 0 0.049 0.0252
12 1.0461 -0.41759 0 0 0.0854 0.0224
13 1.0364 -0.41882 0 0 0.189 0.0812
14 0.99617 -0.4405 0 0 0.2086 0.07

Les pertes totales : Puissance active : 0.29455 [pu] = =29.455MW

46
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE]

profile de module de tension

1.1
Vinit

1.08

1.06
V [p.u.]

1.04

1.02

0.98
1 2 3 4 5 6 7 8 9 10 11 12 13 14
Les jeux de barre
Fig.4.5 : Profile de module de tension à l’état initial

On observe de résultats et des courbes précédents que :

 Les tensions à des valeurs diminuées quand on a loin des générateurs c.à.d. chutes de
tension (pertes au niveau des lignes de transmissions).
Maintenant on applique la méthode CPF

Etapes d’étude

A partir la page d’accueil de PSAT

 Cliquer sur « Edit »

 Choisir « CPF Setting »
 Poser les limites de puissance réactive Q (l’étude est la stabilité statique de tension au
niveau de tous les JB c.à.d. on intéresse à puissance réactive qui influer sur la module de
la tension si pour ça on poser sur la limite de Q dans les paramètres).

47
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE
[ANLYSE STABILITE]

Fig.4.6 : La modification du paramètre de CPF

Pour examiner l’influence de la surcharge dans le réseau électrique sur la stabilité de ce dernier
et identifier la limite de cette surcharge on a adopté sur :
 Augmenter graduellement les valeurs de charge au niveau de toutes les charges (avec la
méthode de CPF)
 Remarquer l’influence de cette variation sur les résultats du PSAT ((static report)
 Comparer avec le résultat initial
 C’est la limite de surcharge obtenir extrait la profile de module de tension et puissance
réactive.

5.1.1. Analyse de réseau avec la méthode

méthod de CPF
Après la simulation par la méthode de CPF on obtient les résultats suivants
*De « command Window » :
 On remarque que les JB n°3, 6,8 (de compensations) sont fonctionnent le maximum de
puissance réactive Q.
 Le générateur PV au JB n°02 à fonctionne à maximum de puissance réactive générée à
lambda = 0.70035.
= 1.1896.

 λ max : Le facteur de charge maximal que le réseau peut supporter sans perdre sa
stabilité c’est à dire après cette valeur le réseau sera instable (physiquement
(physiquement pas de

48
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE]

fonctionnement, mathématiquement la divergence de programme de l’écoulement de

puissance).
*De « static report »

Tableau 4.3. Les résultats de l’écoulement de puissance l’état à la marge de stabilité-modèle

14JB.

Jeu de La tension Angle Pgen Qgen Pcharge Qcharge

Barre (JB) de JB [pu] [rad] [pu] [pu] [pu] [pu]
01 1.06 0 4.6705 3.5422 0 0
02 0.88497 -0.15939 0.47274 0.5 0.35905 0.21013
03 0.71859 -0.53518 0 0.4 1.5586 0.31437
04 0.71178 -0.38894 0 0 0.79089 0.06618
05 0.74426 -0.31289 0 0 0.12575 0.02647
06 0.62884 -0.75991 0 0.24 0.18531 0.12409
07 0.63902 -0.64085 0 0 0 0
08 0.69947 -0.64085 0 0.24 0 0
09 0.57545 -0.80014 0 0 0.4881 0.27466
10 0.56193 -0.82024 0 0 0.14891 0.09597
11 0.58444 -0.80054 0 0 0.05791 0.02978
12 0.58014 -0.83634 0 0 0.10093 0.02647
13 0.56312 -0.84235 0 0 0.22337 0.09597
14 0.50628 -0.92768 0 0 0.24653 0.08273

49
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE]

profile de module de tension

1.2
Vinit
Vcpf
1.1

0.9
V [p.u.]

0.8

0.7

0.6

0.5
1 2 3 4 5 6 7 8 9 10 11 12 13 14
Les jeux de barre

Fig.4.7 : La comparaison entre V init et V cpf de 14JB

On remarque que : à l’état d’étude de CPF il y a des chutes de tension plus que l’état initial
parce que la méthode de CPF fonctionne de augmente graduellement les valeurs de charges au
niveau de chaque JB.

Tableau 4.4 : Les réserves de la puissance réactive des générateurs

∆ = −

Générateurs JB n°02
Puissance Q
0.5
0.5
∆ 0

50
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE]

profile de module de tension

1.4
Vinit
Vcpf
1.2

0.8
V [p.u.]

0.6

0.4

0.2

0
1 2 3 4 5 6 7 8 9 10 11 12 13 14
Les jeux de barre

Fig.4.8 : Histogramme de variation du module de tension entre l état initiale et l état à la marge
de stabilité de14JB.

A partir de cette simulation, on peut constater les JB les plus sensibles (les chutes de tension plus
basses).

Tableau 4.5 : La différence entre l’état initial et l’état à la marge de stabilité en ce qui concerne
le module de tension dans les JB les plus sensibles du modèle de 14JB.

Nbr de 6 7 8 9 10 11 12 13 14
JB
∆ % 44.12% 39.59% 39.05% 43.63% 44.9% 45.07% 46.6% 47.3% 48.9%
∆ : le défferent de la module de tension entre l état initial et l état à la marge de stabilité

1.1

0.9

0.8
X: 1.19
V [pu]

VBus 11 Y: 0.682
0.7
VBus 12
0.6 VBus 13
VBus 14
0.5

0.4

0.2 0.4 0.6 0.8 1

Loading Parameter  (p.u.)

Fig.4.9 : Le module de tension au niveau des JB sensibles en fonction de  -modèle 14 JB

51
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE]

Après cette étude, on constate que ce réseau contient des nœuds (JB) critiques en ce qui concerne
la stabilité statique de la tension. Pour l’intégration des sources à énergies renouvelables, on
cherche l’emplacement optimal pour améliorer la marge de stabilité statique.

 La puissance active demandée totale de réseau=4.2854 pu

On choisit la puissance active totale générée à partir des sources à énergie renouvelable égale à
un certain pourcentage de la puissance demandée totale du réseau. De cette façon, on définit le
pourcentage de pénétration de l’énergie renouvelable en fonction de la puissance demandée
totale.
 Si on a un taux de pénétration de 20% de la puissance demandée totale, la puissance totale
intégrée à partir des sources renouvelables Pintg dans le réseau est :
Pintg=Pg tot-Pd tot -Pertes
 20%( )= =0.85708pu
 25%( )= = 1.07135 pu
La puissance ce décompose sur les quatre JB les plus sensibles (11, 12, 13, 14) par
différents pourcentages ou bien scénarios comme suit :

Fig4.10 : différents scénarios d’étude de réseau 14JB

On commence par 20% du
Premièrement, on règle la source d’énergie renouvelable en vérifiant les grandeurs de base.

52
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE]

Fig.4.11: paramètre de source auxiliaire renouvelable de forme PV du modèle 14JB.

Remarque : les limites du Q sont [0.2 -0.2] à 20% et [0.4 -0.2] à 25%

Tableau 4.6 : Résultats de simulation du premier scénario du réseau 14JB à la marge de

stabilité

scénarios (pu) ∆ ∆ ∆ ∆
1 1.7092 1.4507 0 0 0 0
2 1.7064 1.0675 0 0 0 0
3 1.7115 1.4469 0 0 0 0
4 1.7063 1.4353 0 0 0 0
5 1.7006 1.0732 0 0 0 0
6 1.7048 1.4425 0 0 0 0
7 1.7146 1.4478 0 0 0 0

53
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE]

Fig4.12 : profile de tension des trois premiers scénarios de 20% du réseau 14JB
profile de tension
1.1

0.9
V [p.u.]

0.8

0.7

0.6 V dcmp4
V dcmp5
0.5 V dcmp6
V dcmp7

0.4
1 2 3 4 5 6 7 8 9 10 11 12 13 14
Les jeux de barre

Fig4.13 : profile de tension des quatre derniers scénarios de 20% du réseau 14JB
Maintenant, on étudie le même réseau pour 25% du avec la même chose par des « DG »
Tableau 4.7 : Résultats de simulation de deuxième scénario du réseau 14JB à la marge de
stabilité.
scénarios ∆ ∆ ∆ ∆
1 1.7599 1.0858 0 0 0 0
2 1.7547 1.0918 0 0 0 0
3 1.7634 1.0637 0 0 0 0
4 1.8491 1.3004 0 0.0919 0 0
5 1.8385 1.2094 0 0.23744 0 0
6 1.8083 1.4819 0 0.03909 0 0
7 1.8243 1.4135 0 0.0488 0 0

54
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE]

profile de module de tension

1.1
V dcmp 2.1
1.05 V dcmp2.2
V dcmp2.3
1

0.95

0.9
V [p.u.]

0.85

0.8

0.75

0.7

0.65

1 2 3 4 5 6 7 8 9 10 11 12 13 14
Les jeux de barre

Fig4.14 : profile de tension des trois premiers scénarios de 25% du 14JB

Voltage Magnitude Profile
1.2

V
dcmp2.4
1.1
V
dcmp2.5
V
dcmp2.6
1
V dcmp2.7

0.9
V [p.u.]

0.8

0.7

0.6

0.5
1 2 3 4 5 6 7 8 9 10 11 12 13 14
Bus #

Fig4.15 : profile de tension des quatre derniers scénarios de 25% du 14JB

A partir ces études on constate
 La chute de tension est minimisée (les modules des tensions ont augmenté par rapport
l’état sans intégration).
 Le facteur de charge a augmenté par rapport à l’état sans intégration c'est-à-dire que la
marge de stabilité c’est augmenté et le réseau peut supporter l’augmentation de la charge.
 Quand le taux d’intégration augmente, la marge de stabilité statique augmente.

55
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE]

5.2. Analyse de RE de 30 JB

Fig.4.16: Schéma topologique du modèle de 30JB

On analyse ce réseau à l’état initial et on obtient les résultats suivants de « static report » :

Tableau 4.8.Caractéristiques du réseau 30JB.

Nombre de Nombre Nombre de ligne Nombre de Nombre de

jeux de barres de charges de transmission générateurs Transformateur

31 21 34 6 8

56
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE]

Tableau4.9. les résultats de l’écoulement de puissance à l’état initial du modèle de 30JB.

Jeu de La tension Angle Pgen Qgen Pcharge Qcharge

Barre (JB) de JB [pu] [rad] [pu] [pu] [pu] [pu]

01 1.06 0 2.1157 -0.16401 0 0

02 1.045 -0.07733 0.4 0.37437 0.217 0.127
03 1.0317 -0.11334 0 0 0.024 0.012
04 1.0248 -0.13639 0 0 0.076 0.016
05 1.01 -0.20995 0.1 0.27702 0.942 0.19
06 1.0183 -0.15961 0 0 0 0
07 1.0072 -0.18943 0 0 0.228 0.109
08 1.01 -0.16632 0.1 0.11332 0.3 0.3
09 1.0511 -0.19709 0 0 0 0
10 1.0384 -0.22862 0 0 0.058 -0.18489
11 1.082 -0.18549 0.12 0.305 0 0
12 1.0375 -0.21556 0 0 0.112 0.075
13 1.071 -0.20044 0.12 0.25724 0 0
14 1.0275 -0.2304 0 0 0.062 0.016
15 1.0282 -0.23275 0 0 0.082 0.025
16 1.0348 -0.22272 0 0 0.035 0.018
17 1.0323 -0.23055 0 0 0.09 0.058
18 1.0196 -0.24362 0 0 0.032 0.009
19 1.0176 -0.24668 0 0 0.095 0.034
20 1.022 -0.24319 0 0 0.022 0.007
21 1.0251 -0.23677 0 0 0.0175 0.112
22 1.0254 -0.23661 0 0 0 0
23 1.0172 -0.24045 0 0 0.032 0.016
24 1.0101 -0.24403 0 0 0.087 0.02313
25 0.99912 -0.23974 0 0 0 0
26 0.98111 -0.24734 0 0 0.035 0.023
27 1.0011 -0.2325 0 0 0 0
28 1.0138 -0.16841 0 0 0 0
29 0.98077 -0.25495 0 0 0.024 0.009
30 0.96901 -0.27109 0 0 0.106 0.019

57
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE]

profile de module de tension

1.15
V init30

1.1

1.05
V [p.u.]

0.95

0.9
0 1 2 3 4 5 6 7 8 9 101112131415161718192021222324252627282930
Les jeux de barre

Fig.4.17 : La tension de modèle de 30JB à l’état initial.

 Les pertes active totale= 0.12167 [pu]

Maintenant on applique la méthode de CPF
5.2.1. Analyse de RE avec la méthode de CPF
Tableau4.10 : Résultats de l’écoulement de puissance à la marge de stabilité du modèle de 30JB
Jeu de La tension Angle Pgen Qgen Pcharge Qcharge
Barre (JB) de JB [pu] [rad] [pu] [pu] [pu] [pu]

01 1.06 0 3.9265 3.0181 0 0

02 0.91179 -0.1237 0.6436 0.5 0.34915 0.20434
03 0.8291 -0.19679 0 0 0.03862 0.01931
04 0.78245 -0.2447 0 0 0.12228 0.02574
05 0.75293 -0.45528 0.1609 0.4 1.5157 0.30571
06 0.74431 -0.30342 0 0 0 0
07 0.72767 -0.39298 0 0 0.36685 0.17538
08 0.73191 -0.32699 0.1609 0.4 0.4827 0.4827
09 0.6839 -0.433 0 0 0 0
10 0.62528 -0.56439 0 0 0.09332 0.03218
11 0.71994 -0.38985 0.19308 0.24 0 0
12 0.65716 -0.52339 0 0 0.18021 0.12067
13 0.70378 -0.46491 0.19308 0.24 0 0

58
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE]

14 0.62503 -0.58737 0 0 0.09976 0.02574

15 0.622 -0.59675 0 0 0.13194 0.04022
16 0.64069 -0.55148 0 0 0.05631 0.02896
17 0.61645 -0.57783 0 0 0.14481 0.09332
18 0.58757 -0.64256 0 0 0.05149 0.01448
19 0.57635 -0.65482 0 0 0.15285 0.05471
20 0.58601 -0.63599 0 0 0.0354 0.01126
21 0.58765 -0.60192 0 0 0.28157 0.18021
22 0.58844 -0.601 0 0 0 0
23 0.58047 -0.6286 0 0 0.05149 0.02574
24 0.54601 -0.63739 0 0 0.13998 0.1078
25 0.54658 -0.64356 0 0 0 0
26 0.48783 -0.68852 0 0 0.05631 0.03701
27 0.57628 -0.62199 0 0 0 0
28 0.71731 -0.3279 0 0 0 0
29 0.49992 -0.74524 0 0 0.03862 0.01448
30 0.45569 -0.85356 0 0 0.17055 0.03057

En appliquant la méthode de CPF, on obtient le facteur de charge maximal =1.6593

profile de module de tension
1.1

V init30
0.9
V cpf30
V [p.u.]

0.8

0.7

0.6

0.5

0.4
0 1 2 3 4 5 6 7 8 9 101112131415161718192021222324252627282930
Les jeux de barre

Fig.4.18 : La comparaison entre le module de tension initial et le module de tension à l’état de la

marge de stabilité de réseau 30JB.

59
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE]

profile de module de tension

1.4
V init30
V cpf30
1.2

1
V [p.u.]

0.8

0.6

0.4

0.2

0
0 1 2 3 4 5 6 7 8 9 10 11 1213 14 15 1617 18 1920 21 2223 24 2526 27 2829 30
Les jeux de barre

Fig.4.19 : Histogramme de variation du module de tension entre l état initiale et l état à la marge
de stabilité du réseau 30JB.
A partir de cette simulation, on peut constater les JB les plus sensibles.

Tableau 4.11 : La déférence entre l’état initial et l’état à la marge de stabilité du module de
tension dans les JB le plus sensibles du modèle de 30JB.

Nbr de 18 19 20 21 22 24 25 26 29 30
JB
∆ % 43.2 44.13 43.6 43.74 43.7 46.41 45.25 49.33 48.65 51.33

1.4

1.2

0.8
V [pu]

VBus 24
0.6 X: 1.612
VBus 25 Y : 0.6339

0.4 VBus 26
VBus 29
0.2
VBus 30

0
0.2 0.4 0.6 0.8 1 1.2 1.4 1.6
Loading Parameter  (p.u.)

Fig.4.20 : Le module de tension au niveau les JB sensibles en fonction de  -modèle 30JB

60
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE]

D’après les résultats de simulation, on constate que ce réseau présent des chutes de tension
énormes dans plusieurs nœuds (JB) à la marge de stabilité.
Les générateurs fonctionnent à leurs limites maximales (JB n°2, 5, 8, 11, 13).
Pour ce réseau du 30JB, la puissance demandée totale =4.6655 (pu)
 20%( )= =0.9331 pu et = [0.2 -0.2]
 25%( )= =1.1664 pu et = [0.4 -0.2]
 ∑ , , , , = 0.84 =84MW
La somme des puissances des sources classique est inférieur par rapport les puissances de
à20% et 25%.
Alors, il faut faire la même chose comme le 14JB pour ce réseau, ces puissance intégrer ce
décompose sur des scénarios comme suit :

Fig4.21 : différents scénarios d’intégration d’énergie renouvelable du réseau 30JB.

On commence par 20% du
Tableau 4.12 : résultats de simulation du premier scénario du réseau 30JB
scénarios ∆ ∆ ∆ ∆
1 2.1707 0.77418 0 0 0.0121 0
2 2.1976 0.85573 0 0 0 0
3 2.1321 0.79044 0 0 0 0
4 2.2438 0.79801 0 0 0.00192 0
5 2.1303 0.8697 0 0 0 0
6 2.129 0.79307 0 0 0 0
7 2.1298 0.83373 0 0 0 0

61
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE]

profile de tension
1.1

1.05
V30.1
1 V30.2
V
0.95 30.3

0.9
V [p.u.]

0.85

0.8

0.75

0.7

0.65

0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
Les jeux de barre

Fig4.22 : profile de tension des trois premiers scénarios de 20% du réseau 30JB

profile de tension
1.1

1.05

1 V
30.4
0.95 V 30.5
V
0.9 30.6
V [p.u.]

V 30.7
0.85

0.8

0.75

0.7

0.65

0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
Les jeux de barre
Fig4.23 : profile de tension des quatre derniers scénarios de 20% du réseau 30JB

Maintenant, on étudue le même réseau pour 25% du .

62
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE]

Tableau 4.13 : Résultats de simulation de deuxième scénario du réseau 30JB

scénarios ∆ ∆ ∆ ∆
1 2.2199 1.4046 0 0.02394 0.02584 0.11613
2 2.2367 1.1853 0 0.16026 0.07585 0.13021
3 2.2011 1.217 0 0.00464 0.05838 0.15611
4 2.2786 1.2283 0 0.10902 0.1021 0.14755
5 2.1812 1.2866 0 0.12924 0.2169 0.09382
6 2.2033 1.2249 0 0.00671 0.16434 0.03436
7 2.1988 1.22 0 0.00165 0 0.21463

profile de tension
1.1

V30.2.1
1
V30.2.2
V30.2.3
0.9
V [p.u.]

0.8

0.7

0.6

0.5

0.4
0 1 2 3 4 5 6 7 8 9 10 1112 13 14 1516 17 1819 20 21 2223 24 25 2627 28 29 30
Les jeux de barre
Fig4.24 : profile de tension des trois premiers scénarios de 25% du réseau 30JB

63
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE]

profile de tension
1.3
V 30.2.4
1.2
V 30.2.5
V 30.2.6
1.1
V 30.2.7
1
V [p.u.]

0.9

0.8

0.7

0.6

0.5
0 1 2 3 4 5 6 7 8 9 101112131415161718192021222324252627282930
Les jeux de barre

Fig4.25 : profile de tension des quatre derniers scénarios de 25% du réseau 30JB

On constate a partir ce simulation de 20% et de 25% du réseau 30JB que les modules de tension
à augmenté conséquence d’intégration des sources auxiliaires distribuées
La marge de stabilité a augmenté ç.à.d. le réseau peut supporter une augmentation auxiliaire des
charges.

Tableau 4.14 : Classification des scénarios optimaux d’analyse statique.

JB Pourcentage Scénarios
20% 7 1.7146
14
25% 4 1.8491
20% 4 2.2438
30
25% 4 2.2786

64
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE]

6. Analyse de la stabilité dynamique du réseau électrique

Dans cette partie consacrée à l’analyse dynamique de réseau du 30JB, on calculer le temps
critique d’isolement le défaut(TCID) de tous les JB sauf les JB des générateurs et JB
d’intégration puis on fait l’extraction du JB qui a la plus petite valeur de TCID, après on intègre
les sources à énergies renouvelable (appliquer les scénarios optimales d’analyse statique) et
calculer le TCID a chaque fois.
L’analyse dynamique contre l’analyse statique il faut existe les types de générateurs et le
générateurs pour le réseau peut fonctionne à l’état dynamique.
6.1. Le temps critique d’isolement le défaut (TCID)
C’est le paramètre le plus décisif dans l’analyse de la stabilité transitoire d’un réseau
électrique.
TCID est le temps maximal pendent lequel le réseau peut supporter le défaut sans perdre
sa stabilité.
6.2. Réseau teste (Réseau d’étude)
Dans cette étude d’analyse de stabilité dynamique, on prend le même réseau de 30JB (norme :
IEEE) comme on déjà fait à l’analyse statique mais dans cette partie on analyse ce réseau en cas
de défaut.
6.3. Représentation la stabilité dynamique du RE de 30JB (en cas de défaut)
L’analyse de ce réseau dans ce cas est pour observe les grandeurs électriques de ce réseau et
comparer les résultats avec l’état d’analyse statique.
6.3.1. Etapes d’analyse
On analyse le réseau à l’état dynamique c.à.d. le cas en présence d’un défaut (court-circuit
triphasé à la terre car est le plus sévère qui affecter le réseau électrique) et on calcule le TCID.
Pour réaliser ces études on procède comme suit
1. Charger le réseau d’étude (30JB IEEE) sur le PSAT
2. Pour crée le défaut ont besoin d’un disjoncteur (Breaker) et le défaut de block
« Fault&operations » de bibliothèque de PSAT comme montre dans la figure
« Fig.4.3 »
3. Crée un défaut au niveau de tous les JB sauf les générateurs et JB d’intégration et en
change le temps d’isolement (fin de défaut) avec le temps de déclanchement de
disjoncteur à chaque fois (on commence le calcule de 1 s).
4. Dans chaque changement on observe les courbes de vitesse angulaire des générateurs de
ce réseau
5. A chaque observation on extraire deux cas :
65
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE]

5.1.Un système stable : quand les valeurs des vitesses angulaire tendent vers une valeur
stable et fixe.
5.2.Un système instable : quand les vitesses angulaire est diverge (un seul générateur
capable faire l’écroulement de réseau) alors système instable.
6. Quand le système perd sa stabilité en calcule le temps avant ce diverge pour préciser le
TCID.
Tableau 4.15 : Classification le TCID sans intégration des DG.
JB 06 07 14 15 17 20 25 27 28
TICD 0.0201 0.034 0.167 0.124 0.122 0.159 0.222 0.216 0.054

D’après ce tableau, le JBn° 06 présente le temps plus bas par rapport aux autres.
On adopte le scénario le plus optimale (de l’analyse de la stabilité statique en présence
des DG et prendre 20% de Pd tot puis 25% de Pd tot) et on utilise la machine asynchrone à
double alimentation d’éolienne comme la montre le figure ci-dessous :
L’avantage de cette machine est :
 Peut absorber et injecter la puissance réactive par comparaison avec d’autre types
des machines d’éolienne.

Fig4.26 : modèle de la machine à double alimentation sous PSAT

On régler cette source éolienne comme suit :
 De 20% :
 Le JB 24= Pintg=40% Pd tot= 0.37324 pu
- la puissance de base (Pb)=38 MVA
- Park éolienne= 0.37324/1.5 = 25 éoliennes
 Le reste des JB Pintg =20% Pd tot =0.1866 pu
- La puissance de base (Pb)=19 MVA
- Park éolienne =13 éoliennes
66
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE]

 De 25% :
 Le JB 24= Pintg= 40%de Pd tot= 0.46656 pu
- La puissance de base =47 MVA
- Park d’éolienne =31 éoliennes
 Le reste des JB Pintg=20% Pd tot =0.2333 pu
- La puissance de base =24 MVA
- Park d’éolienne =16 éoliennes

Tableau 4.16 : TCID de JB06 après l’intégration

Pintg Initial 20%Pd tot 25%Pd tot
(s) 0.0201 0.355 0.49

De 20% Pd tot

1.4 1.5
Vitesse angulaire de synchronisme [rad/s]

Vitesse angulaire de synchronisme [rad/s]

Syn 1 Syn 1
1.2 Syn 2 Syn 2
Syn 3 Syn 3
1
Syn 4 1 Syn 4
0.8 Syn 5 Syn 5
T=1.355
Syn 6 T=1.356 Syn 6
0.6
mDfig 1 0.5 mDfig 1
0.4
mDfig 2 mDfig 2

0.2 mDfig 3 mDfig 3

mDfig 4 mDfig 4
0 0
0 5 10 15 20 0 5 10 15 20
time (s) time (s)

Fig4.27 : variation de vitesse angulaire du rotor avec le temps de défaut du 20% de Pd tot

67
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE]

De 25% Pd tot

1.4 1.5
Vitesse angulaire de synchronisme

Vitesse angulaire de synchronisme

Syn 1 Syn 1
1.2 Syn 2 Syn 2
Syn 3 Syn 3
1
Syn 4 1 Syn 4
0.8 Syn 5 Syn 5
T=1.50 s
0.6 Syn 6 Syn 6
T=1.49 s mDfig 1 0.5 mDfig 1
0.4
mDfig 2 mDfig 2

0.2 mDfig 3 mDfig 3

mDfig 4 mDfig 4
0 0
0 5 10 15 20 0 5 10 15 20
time (s) time (s)
Fig4.28 : variation de vitesse angulaire du rotor avec le temps de défaut du 25% de Pd tot.

A partir les deux cas on remarque que :

 Les sources classiques (le générateur n°1 et le générateur n°3) sont les plus sensibles pour
ce cas de défaut.
 Le TCID du JB06 a augmenté quand la puissance intégrer augmente c.à.d. le réseau et
peut supporter le défaut et permette de faire les réparations pour éliminer le défaut.
 L’intégration des énergies éoliennes avec le réseau classique améliorer la marge de
stabilité statique aussi amélioré les caractéristiques du réseau à la stabilité dynamique.

A partir de ces études, on fait l’intégration des sources photovoltaïque avec les sources éoliennes
et on le cherche est que le réseau plus améliorer ou non ? .

7. Intégration des sources à énergie solaire dans le réseau de 30JB

Dans cette partie d’étude on utilise le MATLAB 2009b avec le PSAT 2.1.9 pour l’intégrer le
photovoltaïque dans le réseau électrique, il faut noter que le modèle de photovoltaïque ne
fonctionne pas à l’analyse statique, juste à l’analyse dynamique et avec ce type de forme PV

On choisit les JB n° (19, 21 et 25) pour intégrer les sources de photovoltaïques

68
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE]

 Pintg =20% Pd tot =0.9331 pu

0.9331/7=0.1333 pu (à chaque JB)

 Pour l’éolienne :
- La puissance de base =14MVA
- Park d’éolienne =9 éoliennes
 Pour photovoltaïque
- Puissance de base =13.33 MW

Sans défaut

avec photovoltaïque/ sans défaut

1.4
Syn 1
vitesse angulaire de synchronisme

1.2 Syn 2
Syn 3
1
Syn 4

0.8 Syn 5
Syn 6
0.6
mDfig 1

0.4 mDfig 2
mDfig 3
0.2 mDfig 4

0
0 5 10 15 20
time (s)

Fig.4.29 : La vitesse angulaire de synchronisme avec Photovoltaïque sans défaut

Dans le cas de défaut on à des courbes suivants :

69
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE]

avec photovoltaïque/ T=1.135 avec photovoltaïque / T=1.136

1.4 1.4

Vtesse angulaire de synchronisme

Syn 1 Syn 1
Vitesse angulaire de synchronisme

1.2 Syn 2 1.2 Syn 2

Syn 3 Syn 3
1 1
Syn 4 Syn 4
0.8 Syn 5 0.8 Syn 5
Syn 6
Syn 6
0.6 0.6
mDfig 1
mDfig 1
0.4 mDfig 2
0.4 mDfig 2
mDfig 3
0.2
mDfig 3
mDfig 4 0.2
mDfig 4
0 0
0 5 10 15 20 0 0.5 1 1.5 2
time (s) time (s)
Fig 4.30 : La vitesse angulaire de synchronisme en présence de PV du 20% de Pd tot.

 Pintg=25%Pd tot=1.16644 pu

1.16644/7=0.1666 pu (à chaque JB)

 pour l’éolienne :
- La puissance de base =17 MVA
- Park d’éolienne =11 éolienne
 Pour photovoltaïque
- Puissance de base =16 .66 MW

70
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE]

avec photovoltaïque /T=1.148 avec photovoltaïque /T=1.149

1.4 1.4

Vitesse angulaire de synchronisme

Syn 1 Syn 1
Vitesse angulaire de synchronsme

1.2 Syn 2 1.2 Syn 2

Syn 3 Syn 3
1 1
Syn 4 Syn 4
0.8 Syn 5 0.8 Syn 5
Syn 6
Syn 6
0.6 0.6
mDfig 1
mDfig 1
0.4 mDfig 2
0.4 mDfig 2
mDfig 3
0.2 mDfig 3
mDfig 4 0.2
mDfig 4
0 0
0 5 10 15 20 0 0.5 1 1.5 2 2.5 3 3.5
time (s) time (s)
Fig.2.31 : La vitesse angulaire de synchronisme en présence de PV du 25% de Pd tot.

Tableau 4.17 : Classification de TCID après l’intégration de Photovoltaïque avec l’éolienne

P intg initial 20% Pd tot 25% Pd tot

éolienne Eolienne+PV éolienne Eolienne+PV
0.0201 0.355 0.135 0.49 0.148

On remarque que :
 Quand l’intégration de l’éolienne seul le TCID va augmente ç.à.d. le réseau capable du
supporte le défaut et les organes de protections peut fonctionne pour éliminer le défaut
 L’intégration de PV avec l’éolienne fait une diminution de TCID car le système
photovoltaïque n’est pas supporter l’énergie réactive qui influence sur la module de
tension ce pour cela les performances de réseau à dégradé.

71
 Chapitre IV [ANLYSE DES RESULTATS DE STABILITE]

8. Conclusion

Nous avons présenté dans ce chapitre, une étude de l’amélioration de stabilité des réseaux du
14JB et du 30JB à l’état statique, puis l’amélioration de 30JB à l’état dynamique.

Dans la première partie, on commence par l’analyse statique de 14JB à l’état normal et pour
améliorer la marge de stabilité statique on cherche l’emplacement optimal pour intégrer des
sources à énergie éolienne à différents scénarios. Même étude sera présentée pour le réseau
30JB.

Après l’analyse statique, on passe à deuxième partie c’est l’analyse dynamique, on

commence cette partie avec la détermination de TCID dans le réseau de 30JB, puis l’intégration
aussi des sources éoliennes et déterminer le TCID. Cette intégration a donné des bons résultats
par rapport à l’état normal concernant l’amélioration du TCID. Comme un choix, on intègre des
sources à énergie solaires de forme photovoltaïque avec les sources d’énergie éolienne pour voir
l’influence de ces sources solaires sur la stabilité dynamique, ces dernières ont données des
résultats insatisfaisants concernant à la stabilité dynamique.

72
 CONCLUSION GENERALE

Dans ce travail, nous avons présenté une approche de formulation de la méthode itérative
de Newton-Raphson pour résoudre le problème de l’écoulement de puissance nécessaire pour
l’analyse des performances des réseaux électriques.

L’analyse des performances du réseau électrique en présence de deux types de source

d’énergie renouvelable a été accomplit. Deux types de sources d’énergie renouvelable sont à
considérer, la 1ere c’est la source à énergie éolienne qui à donnée des bons résultats d’analyse
statique et d’analyse dynamique, la 2eme c’est source d’énergie solaire photovoltaïque, mais
cette dernière ont données des résultats insatisfaisant concernant à la stabilité dynamique.

D’après toutes ces analyses on a soulevé les points suivant :

 Le choix de l’emplacement d’intégration est important.

 Comprendre et maitriser les impacts de l’intégration sur la stabilité (stabilité statique et
stabilité dynamique).
 Les avantages et des impacts d’intégration sur la stabilité de réseau électrique.
 L’emplacement de défaut influe sur l’intégration.
 Découvrir les impacts et l’influence de sources photovoltaïque sur la stabilité dynamique
du réseau.

Dans cette étude, on vient d’analyser ces impacts et on à proposé des solutions :

 L’intégration des sources d’énergie renouvelable doit être accompagnée des dispositifs
de compensation de l’énergie réactive pour garder un bon profil de tension.
 L’intégration de ces énergies peut contribuer à la congestion de quelque ligne de réseau
électrique.

Enfin, Comme perspective, des études plus avancées à envisager en introduisant des
autres types d’énergie renouvelable et utilise le système FACTS pour améliorer la stabilité avec
des applications sur le réseau algérien et des autres réseaux plus large.

73
 Bibliographie

A [ACK 05] : T. Ackermann, “Wind Power in Power Systems”, John Wiley & Sons, Ltd,
ISBN 0-470-85508-8, 2005.
[ACT 18] : A.Fatiha, « une production globale de 4.050 mégawatts visée »,
INTERNATIONAL L’ACTUEL, N°192 Janvier 2018, pp06.
[ARB 06] : Abdelaziz Arbaoui « Aide a la décision pour la définition d’un
système éolien adéquation au site et à un réseau faible ». Thèse de doctorat
de l’Ecole Nationale Supérieure d’Arts et Métiers 2006.
[ART 00] : Arthur R.Bergen et Vijay Vihal, « Power Systeme Analysis », 2000
B [BEK 13] : Bekri oum elfadhel loubaba, « Combution à l’étude des systèmes FACTS et
leurs emplacements optimaux dans le réseau électrique », Thèse de doctorat
en Electrotechnique, Université Djilali Laibes de Sidi Bel.Abbès, 2013.
[BHL 91] : B.Ha Lee and K.Y.Lee, « Astudy on Voltage Collapse Mechanisme in
Electric Power Système », IEEE Tren.PS.Vol.6N°3,PP.966-974, Augest
1991.
[BOU 07] : Boukadoum Aziz, « Contribution à l’analyse et la réduction de la pollution
harmonique dans le réseau électrique », Mémoire de magister, Université
Badji Mokhtar, Annaba, 21octobre2007.
C [CAM 03] : Haritza CAMBLONG « minimisation de l’impact des perturbations
d’origine éolien dans la génération d’électricité par des aérogénérateurs a
vitesse variable ». Thèse de doctorat de l’Ecole Nationale Supérieure
d’Arts et Métiers 2003.
G [GIP 07] : Paul GIPE « Le Grand livre de l'éolien », édition du Moniteur, 27 février
2007
[GOU 12] : D. LE GOURIERES. « Eolienne : Théorie conception et calculs pratiques
des installations », Edition Evrolles, Energie, 2012
[GUP 05] : J.B, GUPTA, « A cours in power systems », 14 juillet 2005
H [HEI 05] : S. Heier, “ Grid Integration of Wind Energy Conversion Systems”, 2nd ed.
Wiley, ISBN: 0-470-86899-6, 2005.
[HLA 84] : J. Hladik , “ Énergétique éolienne “, Éd. Masson et Presses de l’université
du Québec,1984
K [KER 13] : Kerbaa amel, « etude de l’influence des systèmes FACTS sur la qualité de

74
 Bibliographie
l’énergie électrique », Mémoire de fin d’étude Master, Université Mohamed
Khider, Biskra, 02 Juin 2013.
M [MER 14] : Merzougui Nassima, « gestion de congestion dans les réseaux électriques »,
Mémoire de fin d’étude Master, Université Mohamed Khider, Biskra 03
Juin 2014.
[MOS 08] : Mostefaoui Abdelkader, « etude la répartition optimale des puissances dans
un réseau électrique par la programmation non linéaire», Mémoire de
Magister, Université Hassiba Benbouali, Chlef. 2008.
P [PER 11] : PROGRAMME DES ENERGIES RENOUVELABLES ET DE
L’EFFICACITE ENERGETIQUE .mars 2011 Ce document a été élaboré
par le ministère de l’énergie et des mines Conception et réalisation
SATINFO Société du Groupe Sonelgaz.
R [RAI 02] : B. Raison et al. “L’énergie électrique éolienne, Partie I : Présentation
générale et approche probabilistique,” Revue Internationale de Génie
Electrique, vol. 5/3-4, pp.405-484, 2002.
S [SAD 14] : Sadi Mohammed el-amine et Hamlat Zohra, « L’impact Des Energies
Renouvelables Sur Le Réseau Electrique Ouest Algérien 220 KV »,
Mémoire de Fin d’Etudes Master (LMD), Université D r. Tahar Moulay
de Saïda, 26-06-2014.
W [WAG 10] : H.J. WAGNER. Livre Energy Technologies 2010, Sub volume C :
Renewable Energy
Z [ZER 10] : Zerguine Bilal, « MODELISATION D’UN SYSTEME DE
PRODUCTION ELECTRIQUE PAR LA CELLULE
PHOTOVOLTAIQUE », Mémoire de Magister, Université BADJI
MOKHTAR, Annaba, 2010

75