République Algérienne Démocratique et Populaire

Ministère de l'Enseignement Supérieur et de la Recherche Scientifique

Université 20 Août 1955-Skikda

Mémoire de fin D’études en vue de l’obtention du diplôme

MASTER 2 en Réseau Electrique

Thème:
Dispatching Economique Environnemental
d’énergie électrique par Essaim de particule

Soutenu le: 25/06/2018

Réaliser par: Président: LAOUAFI Farida M.A.A

ASLOUNI Sofiane Encadreur: MALLEM Aicha M.A.A
KOUICEM Messaoud Examinateur: LEKHECHINE Salima M.C.B

Année : 2017/2018
 PLAN DE TRAVAIL

Introduction
Objectif de travail
Dispatching Economique

Simulation et résultats
Conclusion

2
 Introduction
Toute entreprise électrique a pour mission technique de la
couverture de la charge de son réseau avec un minimum de
défaillances, avec une qualité acceptable de l'énergie livrée.
Cette mission doit être menée à bien à chaque instant présent et
futur sous la qualité et du coût. Les problèmes qui peut
rencontrer dans l’entreprise c’est la distribution de l'économie
d'énergie, pour réduire ce problème, les unités de production
doivent être augmentées mais ce n'est pas la meilleure solution.
Le problème du Dispatching Economique Environnemental
(DEE) est un problème d’optimisation qui consiste à minimiser
les fonctions des coûts de production des unités de production et
l’émission.

3
 Objectif
Ce présent travail nous a permis la minimisation du coût
de production de l’énergie électrique et l’émission qui
est connue du point de vue économique sous le nom
«Dispatching Economique Environnemental».
L’objectif principal de ce travail était de résoudre ce
problème avec et sans considération des pertes de
transmission par la méthode des coûts marginaux et la
PSO pour l’amélioration de la performance des systèmes
d’énergie.

4
 Dispatching Economique
Le Dispatching Economique est l'un des problèmes
d'optimisation dans les systèmes d'alimentation qui ont
l'objectif de diviser la demande de puissance entre les
générateurs en ligne économiquement tout en satisfaisant
diverses contraintes.
Depuis le coût de la production d'électricité est exorbitante,
une répartition optimale économise une quantité
considérable d'argent.

5
Puissance
à l’entrée
Fi Pgi   ai Pgi2  bi Pgi  ci
(DA/h)

Pgi min Pgi max Puissance à la sortie du

générateur (MW)

Figure III.1. Courbe de coût typique (entrée-sortie)

d’un générateur.

6
Minimisation des coûts de génération

Le but principal du Dispatching Economique est la minimisation du coût

de production de chaque générateur.
Cette minimisation peut être traduite par la condition suivante:

Avec:

7
 Dispatching Economique
Environnemental

DISPATCHING ECONOMIQUE ENVIRONNEMENTAL

Le problème du Dispatching Environnemental consiste à minimiser la

fonction des émission des gaz toxiques dans l'atmosphère provenant des
centrales de production, comme les centrale nucléaire, les centrales
thermiques à combustibles fossiles, ...etc.

8
 FORMULATION MATHEMATIQUE

La formulation mathématique peut être décrit de

la manière suivante:

Avec:

Ou:

9
 Méthode des coûts marginaux
Méthode des coûts marginaux sans considération des
pertes
La méthode des coûts marginaux consiste à construire le
Lagrangien qui tient compte des contraintes d'égalité et des
contraintes d'inégalité:

10
Méthode des coûts marginaux avec considération des
pertes
Le traitement mathématique dans ce cas est le même que dans
un réseau sans pertes, sauf que nous devons maintenant
intervenir les pertes de transmission comme une contrainte
additionnelle dans l’équation de H, Pour résoudre le problème
du dispatching économique avec pertes, on peut formulée le
lagrangien avec pertes en utilisant l'équation comme suit:

11
Organigramme de la méthode
des coûts marginaux

12
 Début

Initialisation de , et

Calculer les perte initiales et  initial

Calculer les pertes et les puissances générées

Si: Pour:
Si: Pour:

Calculer l’erreur E:

Non Oui
<ɛ

Calculer 
Afficher les puissances et le coût total

FIN

13
PSO

14
 Essaim de particules (PSO)
L’optimisation par l’Essaim de particules (PSO) repose sur un
ensemble d’individus originellement disposés de façon
aléatoire et homogène, que nous appellerons dès lors des
particules.

L’Essaim de particules est constitué de n particules et la

position de chaque particule représente une solution dans
l'espace de recherche.
Les particules changent d’état selon les trois principes
suivants:

15
Garder son inertie.

Changer d'état en fonction de sa position la plus

optimiste.

Changer d’état selon la position la plus optimiste du groupe.

16
La formule classique peut être décrite de la façon suivante:

: Position de la particule
: Vitesse de la particule

: Meilleure fitness obtenue pour la particule

ω : Coefficient d'inertie,
c1 et c2: Coefficients d'accélération,
r1 et r2: Sont deux nombres aléatoires tirés uniformément de
l’intervalle [0,1] à chaque itération et pour chaque dimension.

17
Organigramme de la
méthode PSO

18
 Début

Initialisation de population

Calcul de la fitness de chaque particule

Calcul de la meilleure fitness de la population actuel

Calcul de la meilleure fitness de toutes les populations

Calcul de la vitesse
Calcul de la position
T=T+1
Calcul de la fonction sélective de la nouvelle
population

Test

Résultat

Fin

19
SEMULATION ET
RESULTATS

20
 INTRODCTION

Le but de ce chapitre est de présenter un exemple de

simulation de deux réseaux électriques, pour démontrer
l’utilité et l’efficacité de la méthode PSO dans la résolution du
problème de Dispatching Economique Environnemental dans
les réseaux électriques.
Le programme à été développé sous l’environnement
(Matlab R2017b) pour optimiser la fonction du coût et la
fonction émission.

21
 Application sur un réseau de 14nœuds
Ce réseau électrique IEEE 14 nœuds est constitué de 5
centrales électriques de production et 20 lignes.

Figure 1. Schéma unifilaire du réseau électrique de14 nœuds.

22
Paramètres des générateurs du réseau

Les données des générateurs et les limites minimales et

maximales des puissances actives des cinq générateurs et les
coefficients caractéristiques des émissions sont représentées
respectivement dans les tableaux 1 et 2.

23
 Paramètres des générateurs du réseau

Tableau 1. Coefficients de la fonction du coût des générateurs.

N° du P min P max
générateur ($/MWh2) ($/MWh) ($/h) [MW] [MW]
1 0.00315 2.00 0.00 50 200
2 0.0175 1.75 0.00 20 80
3 0.0625 1.00 0.00 15 50
6 0.00834 3.25 0.00 10 35
8 0.025 3.00 0.00 10 30

24
 Tableau 2. Coefficients caractéristiques des émissions.

N° du P min P max
générateur (Kg/MWh2) (Kg/MWh) (Kg/h) [MW] [MW]

1 0.0126 -0.90 22.983 50 200

2 0.0200 -0.10 25.313 20 80

3 0.0270 -0.01 25.505 15 50

6 0.0291 -0.005 24.900 10 35

8 0.0290 -0.004 24.700 10 30

25
Application de la méthode des coûts marginaux

Cette partie et consacré à l'étude du Dispatching Economique

Environnemental avec et sans considération des pertes.

Application sans considération des pertes

Le tableau 3 présente les résultats qu'on a obtenus par la

méthode des coûts marginaux sans considération des pertes.

26
Tableau 3. Résultats obtenus sans considération des pertes par
la méthodes des coût marginaux.
Résultats Méthode Coût marginaux
PG1 [MW] 174.3172

PG2 [MW] 42.2351

PG3 [MW] 16.8127

PG6 [MW] 17.0534

PG8 [MW] 10.0000

Puissance actives générée [MW] 260.4185
Fonction du coût de génération [$/h] 674.3083
Emission [Kg/h] 399.5393
27
La représentation graphique des résultats est illustrée par la
figure 2
Pgmin
Pg
200
Pgmax
Puissance active [MW]

150

100

50

0
1 2 3 4 5
Nombre générateur

Figure 2. Représentation graphique des puissances générées sans

considération des pertes par la méthode des coûts marginaux.

28
Application avec considération des pertes
Le tableau 4 présente les résultats qu'on a obtenus par la méthode des
coûts marginaux avec considération des pertes.
Tableau 4. Résultats obtenus avec considération des pertes par la méthodes des
coûts marginaux.

Méthode Cout Marginaux

PG1 [MW] 178.5145
PG2 [MW] 43.0945
PG3 [MW] 17.1541
PG6 [MW] 19.3814
PG8 [MW] 10.0000
Puissance actives générée [MW] 268.1446
Pertes actives [MW] 8.7522
Fonction du coût de génération [$/h] 699.4948
Emission [Kg/h] 418.5686
29
La représentation graphique des résultats est illustrée par la
figure 3
Pgmin
Pg
200 Pgmax
Puissance active [MW]

150

100

50

0
1 2 3 4 5
Nombre générateur

Figure 3. Représentation graphique des puissances générées avec

considération des pertes par la méthode coût marginaux.

30
Application de PSO
Dans cette partie d'étude de Dispatching Economique Environnemental
par PSO. Pour l'algorithme PSO il faut faire un choix des valeurs basé
sur l'expérience de l'application de cet algorithme introduits dans le
calcul, on sélectionne le choix suivant:

Nombre d'itération= 100.

Taille de population= 50.

Coefficients d'inertie (w)= 0.5 < 0.5 < 0.9

Constants de la vitesse: C1+C2 ≤ 4

31
Application sans considération des pertes
Le tableau 5 présente les résultats qu'on obtenus par l'application de PSO
sans considération des pertes.
Tableau 5. Résultats obtenus sans considération des pertes par PSO.
Méthode PSO

C1 = 2 C1 = 1 C1 = 1.8 C1 = 1.2 C1 = 2.5

Paramètre
C2 = 2 C2 = 3 C2 = 2.2 C2 = 2.8 C2 = 1.5

PG1 130.2061 130.1616 127.3730 129.9361 130.9762

PG2 52.4663 52.7470 55.6626 52.5372 52.2611
PG3 32.6303 32.4669 31.2439 32.9245 32.7908
PG6 23.6281 23.1409 24.1266 23.6226 23.4163
PG8 20.0692 20.4837 20.5939 19.9797 19.5557

Puissance actives générée

259.0000 259.0000 259.0000 259.0000 259.0000
[MW]

Fonction du coût de
704.7053 704.6677 705.3873 705.3194 704.1399
génération [$/h]

Emission [Kg/h] 325.7886 325.8031 324.6661 325.6949 326.6255

La représentation graphique des résultats est illustrée par la
figure 4

200
Puissance active [MW]

150

100

50

0
1 2 3 4 5
Nombre de générateur

Figure 4. Représentation graphique des puissances générées sans

considération des pertes par PSO.

33
Application avec considération des pertes
Le tableau 6 présente les résultats qu'on obtenus par l'application de PSO
avec considération des pertes.
Tableau 6. Résultats avec considération des pertes par la méthode PSO.
Méthode PSO
C1 = 2 C1 = 1 C1 = 1.8 C1 = 1.2 C1 = 2.5
Paramètre
C2 = 2 C2 = 3 C2 = 2.2 C2 = 2.8 C2 = 1.5
PG1 132.5308 132.1324 134.1388 133.2302 133.3929
PG2 48.1965 53.3648 54.2775 54.3530 52.7741
PG3 34.8954 34.4382 33.2697 33.3082 34.1903
PG6 32.2817 24.6832 23.9462 23.8316 24.8335
PG8 17.5724 21.4357 20.5256 21.3202 20.9067

Puissance actives générée [MW] 265.4768 266.0544 266.1578 266.0431 266.0974

Pertes actives [MW] 6.4533 7.0536 7.1578 7.0431 7.0974

Fonction du coût de génération
730.5760 733.3886 728.6701 729.3549 730.7004
[$/h]
La représentation graphique des résultats est illustrée par la
figure 5
Pgmin
Pg
200 Pgmax
Puissance active [MW]

150

100

50

0
1 2 3 4 5
Nombre générateur

Figure 5. Représentation graphique des puissances générées avec

considération des pertes par PSO.

35
Convergence de l’algorithme PSO

La figure 6 montre l’évolution de la fonction du coût de génération par

PSO.
1000

950

900
Fonction coût [$/h]

850

800

750

700

650
0 10 20 30 40 50 60 70 80 90 100
Itérations

Figure 6. Evolution de la fonction du coût.

36
La figure 7 montre l’évolution de la fonction émission par
PSO.
600

550

500
Fonction émission [Kg/h]

450

400

350

300

250
0 10 20 30 40 50 60 70 80 90 100
Itérations

Figure V.7. Evolution de la fonction émission.

37
La figure 8 montre l’évolution de la fonction multi-objective
par PSO.
104
6

5.5

5
Fonction multi-objective [$/h]

4.5

3.5

2.5

1.5

0 10 20 30 40 50 60 70 80 90 100
Itérations

Figure 8. Evolution de la fonction multi-objective.

38
Comparaison des résultats entre la méthode des coût marginaux et
la méthode PSO
Les tableaux 7 et 8 présente respectivement une comparaison entre les
résultats des deux méthodes sans considération des pertes et avec
considération des pertes
Tableau 7. Comparaison entre les résultats des deux méthodes sans considération des
pertes.

Méthodes PSO Coût marginaux

Fonction du coût de génération[$/h] 704.7053 674.3083

Emission [Kg/h] 326.6255 399.5393

Tableau 8. Comparaison entre les résultats des deux méthodes avec considération des
pertes.

Méthodes PSO Coût marginaux

Fonction du coût de génération [$/h] 730.5760 699.4948

Emission [Kg/h] 338.7041 418.5686

 Application sur un réseau de 30 nœuds
Ce réseau est constitué de 6 nœuds générateurs et 24 nœuds de
charge.

Figure 9. Schéma unifilaire du réseau IEEE 30 nœuds.

40
Paramètres des générateurs du réseau

Les coefficients de la fonction quadratique du coût, les limites

minimales et maximales des puissances des six générateurs et
les coefficients de la fonction émission sont présentées
respectivement dans les tableaux 9 et 10.

41
 Paramètres des générateurs du réseau

Tableau 9. Coefficients de la fonction du coût des générateurs.

N° du P min P max
générateur ($/MWh2) ($/MWh) ($/h) [MW] [MW]

1 0.00315 2.00 0.00 50 200

2 0.0175 1.75 0.00 20 80
5 0.0625 1.00 0.00 15 50
8 0.00834 3.25 0.00 10 35
11 0.0250 3.00 0.00 10 30
13 0.02500 3.00 0.00 12 40
 Tableau V.10. Coefficients caractéristiques des émissions.

N° du P min P max
générateur (Kg/MWh2) (Kg/MWh) (Kg/h) [MW] [MW]

1 0.0126 -1.100 22.983 50 200

0.0200 -0.1000 22.313
2 20 80

0.0270 -0.1000 25.505

5 15 50

0.0291 -0.0050 24.900

8 10 35

0.0291 -0.0400 24.700

11 10 30

0.0271 -0.0055 25.300

13 12 40
Application de la méthode des coûts marginaux sur le
réseau
Cette partie et consacré à l'étude du Dispatching Economique
Environnemental avec et sans considération des pertes.

Application sans considération des pertes

Le tableau 11 présente les résultats qu'on a obtenus par la
méthode des coûts marginaux sans considération des pertes.

44
 Tableau 11. Résultats sans considération des pertes par la
méthodes des coût marginaux
Méthode Coût marginaux

PG1 [MW] 189.3336

PG2 [MW] 47.7143

PG5 [MW] 19.3600

PG8 [MW] 10.1920

PG11 [MW] 10.0000

PG13 [MW] 12.0000

Puissance actives générée [MW] 288.5999

Fonction du coût de génération[$/h] 785.3117

Emission [Kg/h] 447.3727

La représentation graphique des résultats est illustrée par la
figure 10

Pgmin
Pg
200 Pgmax
Puissance active [MW]

150

100

50

0
1 2 3 4 5 6
Nombre générateur

Figure 10. Représentation graphique des puissances générées

sans considération des pertes par la méthode des coût marginaux.

46
Application avec considération des pertes
Le tableau 12 présente les résultats qu'on a obtenue par la méthode des
coûts marginaux avec considération des pertes.
Tableau 12. Résultats avec considération des pertes par la méthodes des
coût marginaux.
Méthode Cout Marginaux

PG1 [MW] 179.3821

PG2 [MW] 49.0640
PG5 [MW] 20.4121
PG8 [MW] 21.4779
PG11 [MW] 11.7247
PG13 [MW] 12.0000
Puissance actives générée [MW] 294.0608
Pertes actives [MW] 10.2705
Fonction du coût de génération[$/h] 805.7348
Emission [Kg/h] 426.9545
La représentation graphique des résultats est illustrée par la
figure 11
Pgmin
Pg
200 Pgmax
Puissance active [MW]

150

100

50

0
1 2 3 4 5 6
Nombre générateur

Figure 11. Représentation graphique des puissances générées

avec considération des pertes par la méthode des coût marginaux.

48
Application de PSO
Dans cette partie d'étude de Dispatching Economique Environnemental
par PSO. Pour l'algorithme PSO il faut faire un choix des valeurs basé
sur l'expérience de l'application de cet algorithme introduits dans le
calcul, on sélectionne le choix suivant:

Nombre d'itération= 100.

Taille de population= 50.

Coefficients d'inertie (w)= 0.5 < 0.5 < 0.9

Constants de la vitesse: C1+C2 ≤ 4

49
Application sans considération des pertes
Le tableau 13 présente les résultats qu'on obtenus par l'application de
PSO sans considération des pertes.
Tableau 13. Résultats obtenus sans considération des pertes par PSO.
Algorithme PSO
C1 = 2 C1 = 1 C1 = 1.2 C1 = 1.9 C1 = 2.5
Paramètre
C2 = 2 C2 = 3 C2 = 2.8 C2 = 2.1 C2 = 1.5

PG1 128.3335 129.7146 122.1653 129.8787 129.4109

PG2 52.2067 52.2853 48.4504 51.3407 51.5292

PG5 32.7064 32.5385 35.2934 33.2141 33.1879

PG8 23.3206 22.6559 25.0409 22.0544 23.3281

PG11 20.3753 20.0695 24.0208 20.4017 19.1508

PG13 26.4575 26.1362 28.4291 26.5104 26.1508

Puissance actives générée

283.4000 283.4000 283.4000 283.4000 283.4000
[MW]

Fonction du coût de
805.7547 808.5062 819.5855 806.3948 805.7117
génération [$/h]
La représentation graphique des résultats est illustrée par la
figure 12

Pgmin
200 Pg
Pgmax
Puissance active [MW]

150

100

50

0
1 2 3 4 5 6
Nombre générateur

Figure 12. Représentation graphique des puissances générées

sans considération des pertes par PSO.

51
Application avec considération des pertes
Le tableau 14 présent les résultats qu'on obtenus par l'application de la
méthode PSO avec considération des pertes.
Tableau 14. Résultats obtenus avec considération des pertes par PSO.

Algorithme PSO
C1 = 2 C1 = 1 C1 = 1.2 C1 = 1.9 C1 = 2.5
Paramètre
C2 = 2 C2 = 3 C2 = 2.8 C2 = 2.1 C2 = 1.5
PG1 130.6467 151.4728 154.8514 16V.4235 162.4223
PG2 57.8116 63.4278 47.6309 56.2865 53.3502
PG5 23.5584 21.1887 20.0880 23.9615 17.5044
PG8 29.9021 27.4893 22.0965 12.3013 27.8824
PG11 20.2864 10.0000 27.5847 23.5600 11.6604
PG13 27.1336 18.5892 18.2539 12.0000 18.2763
Puissance actives
289.3389 292.1678 290.5053 293.5327 283.4000
générée [MW]
Pertes actives [MW] 5.9389 7.3555 5.9690 8.4692 7.6738
Fonction du coût de
818.8058 815.1670 819.0908 806.7408 802.3038
génération[$/h]
Emission [Kg/h] 347.5827 390.1945 372.2317 401.7208 391.6992
La représentation graphique des résultats est illustrée par la
figure 13
Pgmin
Pg
200 Pgmax
Puissance active [MW]

150

100

50

0
1 2 3 4 5 6
Nombre générateur

Figure 13. Représentation graphique des puissances générées

avec considération des pertes par PSO.

53
Convergence de PSO
La figure 14 montre l’évolution de la fonction du coût de génération par
PSO.
4000

3500

3000
Fonction cout [$/h]

2500

2000

1500

1000

500
0 10 20 30 40 50 60 70 80 90 100
Itérations

Figure 14. Evolution de la fonction du coût.

54
La figure 15 montre l’évolution de la fonction émission par
PSO.
1200

1100

1000

900
Fonction émission [Kg/h]

800

700

600

500

400

300

200
0 10 20 30 40 50 60 70 80 90 100
Itérations

Figure 15. Evolution de la fonction émission.

55
La figure 16 montre l’évolution de la fonction multi-objective
par PSO.
4000

3500
Fonction multi-objective [$/h]

3000

2500

2000

1500

0 10 20 30 40 50 60 70 80 90 100
Irérations

Figure 16. Evolution de la fonction multi-objective.

56
Comparaison des résultats entre la méthode des coût marginaux et
PSO
Les tableau 15 et 16 présentent une comparaison entre les résultats des
deux méthodes avec et sans considération des pertes.
Tableau 15. Comparaison entre les résultats des deux méthodes sans considération des
pertes.

Méthodes PSO Coût marginaux

Fonction du coût de génération[$/h] 805.7547 785.3117

Emission [Kg/h] 332.7515 447.3727

Tableau 16. Comparaison entre les résultats des deux méthodes avec considération des
pertes.

Méthodes PSO Coût marginaux

Fonction du coût de génération[$/h] 818.8058 805.7348

Emission [Kg/h] 347.5827 426.9545

Interprétation des résultats

Le problème de Dispatching Economique a été résolu par l'application de la méthode

des coût marginaux et PSO. Par comparaison des résultats obtenus on observe que la
méthode PSO et plus efficace par rapport à la méthode des coût marginaux.

Le coût total de production des unités des génération électriques est moins dans
l’application sans considération des pertes, ceci est plus économique mais n’est pas
pratique, pour cette raison l'étude de Dispatching Economique sans pertes reste
toujours théorique.

Les puissances générées sont dans les limites, c'est à dire les contraintes d'égalités et
inégalités sont respectées.

Les résultats et le temps de convergence de l’algorithme PSO dépend du poids

d’inertie, de la taille de population et du nombre de génération.

58
Conclusion

59
 Conclusion

Les résultats de la méthode PSO sont efficaces au niveau de

l’émission et au niveau du coût de production.
L’importance de l’étude de ce problème est l’économie de
l’énergie électrique qui nous ramène à un bilan de puissances
équilibre entre la génération et la consommation avec un coût
optimal et la diminution des émissions atmosphériques.

60