Impact de l’injection des productions PV décentralisées

dans le réseau BT de Lubumbashi

Yves Heri, Moise Nguz

To cite this version:

Yves Heri, Moise Nguz. Impact de l’injection des productions PV décentralisées dans le réseau BT de
Lubumbashi. 2020. �hal-02510716�

HAL Id: hal-02510716

[Link]
Preprint submitted on 18 Mar 2020

HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est

archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents
entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non,
lished or not. The documents may come from émanant des établissements d’enseignement et de
teaching and research institutions in France or recherche français ou étrangers, des laboratoires
abroad, or from public or private research centers. publics ou privés.
Impact of the injection of decentralized PV productions into the low
voltage network of the city of Lubumbashi

Impact de l’injection des productions PV décentralisées dans le réseau

BT de Lubumbashi

Yves HERI*1, Moise NGUZ2

Abstract

The injection of decentralized photovoltaic production into the low-voltage distribution network can contribute to improving
the reliability of the services of public distribution system operators (DSO). This article is based on the study conducted on the
residential block LA MISSION in the Kasungami district located on the outskirts of the city of Lubumbashi in the Democratic
Republic of Congo.
Through the analysis by simulation of the tension profiles on the whole of the network before and after PV injections on
different nodes of the network, we define the optimal powers and the optimal point for the integration of these productions by
basing ourselves on the power flow on the lines.

Résumé

L’injection des productions photovoltaïques (PV) décentralisées dans les réseaux de distribution basse tension peut contribuer
à l’amélioration de la fiabilité des services des opérateurs publics. Cet article a pour principal objectif l’étude des faisabilités de
l’injection de l’énergie photovoltaïque dans le réseau basse tension (BT) de la cellule DE LA MISSION dans le quartier
KASUNGAMI se trouvant à la périphérie de Lubumbashi.
En analysant par simulation les profils de tension sur l’ensemble du réseau sans injections et avec injections sur différents
nœuds, nous définissons les puissances optimales et le point optimal pour l’intégration de ces productions en nous basant sur
le flux de puissance sur les lignes.

Key Words / Mots-clés

Modeling /Modélisation; Integration/ Injection; LV Grid /Réseau BT

Contenu
Abréviations
1. Introduction et contexte

2. Subdivision de l’article Grandeur Unité

Symbole Désignation Nom Symbole
3. Présentation du site PV Photovoltaïque - -
3.1. Présentation du réseau test BT Basse tension - -
P Puissance active Watt W
3.2. Les données de charges
Q Puissance Volt-Ampere VAr
4. Modélisation et simulation réactive réactif
5. Présentation et analyse des résultats SNEL Société Nationale - -
6. Conclusion et interprétation des résultats d’électricité
Références

1 Département de Génie Electrique : Ecole des Ingénieurs Industriels / Université de Lubumbashi, yvesheri@[Link]

2
Ecole des Ingénieurs Industriels / Université de Lubumbashi
 International Journal of Applied Research Vol (N°...),

dans le réseau de Lubumbashi. Les auteurs [12] en étudiant

1. Introduction et contexte
un réseau urbain en Australie montre que l’utilisation des
Suite à une extension géographique et une croissance systèmes de stockage (BESS : Batteries Energy Storage
démographique effrénées au cours des deux dernières System) permet de limiter les contraintes thermiques sur
décennies, la ville de Lubumbashi, en République les lignes car les batteries aident à injecter une puissance
Démocratique du Congo, doit faire face à des coupures active quand le système de génération PV est en sous-
d’électricité intempestives. On observe notamment une production et permettent d’emmagasiner de l’énergie
augmentation sensible de la consommation d’énergie lorsque le système de génération PV est en surproduction.
électrique en basse tension (BT) entre 2007 et 2016, allant
de 13.041.000 kWh à environ 62.848.364 kWh [1, 2]. Cette En Italie, [13] estime que la facture d’électricité des
augmentation a provoqué des fortes perturbations dans le ménages peut-être réduite de 49,5 % en cas d’utilisation
réseau de distribution conduisant ainsi la Société Nationale d’un système PV avec BESS.
d’Electricité (SNEL) à opérer des délestages de longue
durée pour un rééquilibrage de la consommation.
2. Subdivision de l’article

Les habitants de Lubumbashi recourent ainsi à l’utilisation Après la présentation du contexte de l’étude, nous
d’autres sources d’énergie pour palier à ce manque, d’après présenterons premièrement le site de l’étude dans la
l’étude menée par [3] en cas de coupure d’électricité, les section suivante. Cela sera suivi la section 4 par une
groupes électrogènes restent la principale source d’énergie présentation de méthode de modélisation choisie, en
à Lubumbashi. Ces travaux révèlent que 29,2 % de la précisant les hypothèses considérées et nous procéderons
population utilisent les groupes électrogènes, 13,1 % sont à la section 5 à la présentation et l’analyse des résultats des
raccordés à des mini-réseaux alimentés par les groupes simulations de ce modèle.
diesels communément appelés dynamo et seulement 1,9 %
3. Présentation du site
de la population utilisent l’énergie photovoltaïque. A ces
sources d’énergie s’ajoutent l’utilisation du charbon de bois Cette étude sera menée dans la ville de Lubumbashi. Nous
qui cause une forte déforestation des zones environnant nous baserons toute de même plus sur les nouveaux
Lubumbashi. quartiers de la ville car ces quartiers possèdent le plus
faible taux d’électrification et on remarque que ce sont ces
Notre étude consistera à analyser les impacts qu’aurait une quartiers qui ont une plus grande proportion des ménages
injection d’énergie solaire photovoltaïque (PV) possédant des installations PV [14]. Ainsi, nous avons opté
décentralisée dans le réseau de distribution de la ville de pour la cellule DE LA MISSION dans le quartier Kasungami
Lubumbashi avec comme cas d’étude la cellule « De la dans la commune Annexe.
mission » dans le quartier Kasungami. Des études similaires
ont déjà été menées dans plusieurs pays et les résultats Située à 11° 40’ de latitude et 27° 29’ de longitude,
obtenus, nous encouragent dans la démarche. Lubumbashi est la capitale de la province du Haut-Katanga.
La ville a 413.000 habitants en 1973 et 700.000 en 1988.
Au Brésil par exemple, il en ressort que l’injection de PV Récemment la ville comptait 2.786.397 habitants en 2015.
reste techniquement faisable et économiquement Une carte de la ville est donnée à la figure 1.
avantageux pour tous les réseaux de distribution [4 – 8].
Tandis qu’en Allemagne certaines études montrent la
corrélation entre le besoin d’autosuffisance énergétique et
les paramètres techno-économiques devant promouvoir
l’achat d’installation PV pour l’injection dans les zones
résidentielles. Ces études montrent que le désir
d’autosuffisance énergétique des ménages, associé à des
législations incitatives (telle que le Feed-In-Tariff) accroit le
taux de pénétration de production PV dans une zone
résidentielle [9, 10].

Ces mesures d’incitation peuvent découler d’un ajustement

de la législation, pour permettre l’intégration à grande
échelle des productions PV des ménages. Pour l’Espagne
Figure 1. Carte de la ville de Lubumbashi
[11] recommande l’adoption d’un mécanisme de
rétribution afin de rémunérer les ménages ayant un surplus
de production PV. 3.1. Présentation du réseau test

La société nationale d’électricité assure la distribution

Dans notre cas, les mesures d’incitations doivent être
électrique dans notre milieu d’étude via un feeder
accompagnées des procédures et des techniques qui
provenant du poste Katuba. C’est un réseau du type radial
permettent de faciliter l’intégration des productions PV
ayant un transformateur de 500 kVA.
3 International Journal of Applied Research and Technology, Vol (N°), pp……

La canalisation est du type aérien, constituée des poteaux - 𝑝𝑘 , 𝑘𝑢𝑘 , 𝑘𝑠𝑘 respectivement la puissance, le facteur
(nœuds) séparés l’un par rapport à l’autre de 40 m. Sur ce d’utilisation et le facteur de simultanéité du
réseau on retrouve des tronçons avec câbles en aluminium ménage k ;
et des tronçons en cuivre. Les caractéristiques sont reprises - 𝑚 nombre de parcelles connectées au nœud ;
dans le tableau 1. - 𝑛 nombre des ménages dans la parcelle j.

Le tableau 2 présente les résultats de cette collecte.

Tableau 1. Caractéristiques des tronçons du réseau étudié
Section
Tronçon Délimitation Matériaux Tableau 2. Charges du réseau
(mm2)
Puissances
A N0 à N9 Al 35
Nœud Réactives
B N9 à N13 Cu 16 Actives (W)
(VAr)
C N0 à N15 Al 50
D N15 à N18 Cu 16 1 1700 559
2 4240 1393
N = Nœud 3 2700 888
4 1440 474
La figure 2 montre une vue aérienne de la zone d’étude et le
5 2425 798
tracé de la partie considérée du réseau.
6 1040 342
7 1840 604
8 1700 559
9 2600 855
10 1200 395
11 1770 582
12 1660 546
13 780 256
14 1810 595
15 960 317
16 1680 552
17 1672 550
18 1660 546
19 1700 559
20 440 145
Figure 2. Carte de la ville de Lubumbashi 21 2500 822
22 780 256
3.2. Les données de charges 23 1810 595
24 1770 582
Une récolte de données sur l’ensemble du réseau test nous 25 1000 329
a permis de trouver les puissances raccordées sur chaque 26 1672 550
nœud du réseau. La méthodologie considérée pour évaluer 27 3500 1150
ces puissances est la suivante : 28 1660 546
- Un inventaire des équipements des ménages. 29 3100 1019
- Analyse des tableaux d’installations électriques 30 7000 2300
afin de déterminer les coefficients de 31 7800 2565
simultanéités et d’utilisation. 32 1700 559
- Agrégation des puissances des différents ménages 33 6000 1972
connectés au nœud. 34 950 312
- Considérer un facteur de puissance (cos 𝜑) de 35 1400 460
0,95 pour l’ensemble du réseau. 36 2540 835
37 2500 8217
𝑃𝑖 = ∑𝑚 𝑛
𝑗=1((∑𝑘=1 𝑝𝑘 ) × 𝑘𝑢𝑘 × 𝑘𝑠𝑘 ) (1)
4. Modélisation et simulation
𝑄𝑖 = 𝑃𝑖 × tan 𝜑 (2)
Le modèle du réseau a été élaboré partant d’un modèle
existant, IEEE à 15 bus, facile à développer dans
Avec MATLAB/Simulink. Ce modèle de référence a permis de
- 𝑃𝑖 , 𝑄𝑖 puissance active et réactive connectées au modifier les paramètres du réseau ainsi que la structure
nœud i ; afin de retrouver la topologie et les caractéristiques de
notre réseau test. Ainsi deux hypothèses simplificatrices
sont considérées :
4 International Journal of Applied Research and Technology, Vol (N°), pp……

- Le réseau est équilibré (permettant une analyse 5. Présentation et analyse des résultats
simple d’écoulement de puissance).
- La puissance de court-circuit du réseau amont est Nous présentons premièrement les graphiques de
de 1000 MW (valeur prescrite par la SNEL). l’écoulement de puissance dans les lignes (figure 3a) et du
profil de tension (figure 3b) sans injection. On constate
La simulation a été faite en deux étapes, en premier lieu, sans surprise que les deux caractéristiques sont tombantes,
nous avons lancé l’analyse d’écoulement des puissances c’est-à-dire que leurs valeurs diminuent au fur et à mesure
(load-flow) sur notre réseau sans qu’aucune injection n’y qu’on s’éloigne du poste source.
soit raccordée nous avons obtenu deux courbes l’une du
profil de tension en fonction des nœuds et l’autre de La tension va jusqu’à une tension de 0.63 p.u au nœud 13
l’écoulement de puissance active par nœuds. En second qui est déjà de loin inferieure des valeurs recommandées
faire progressivement des injections de trois paliers de dans la norme CEI-60038 qui est de plus ou moins 10% soit
puissance à savoir ; 60, 100 et 150 kW aux quatre nœuds 0.9 p.u minimum et 1.1 p.u maximum. Notons aussi que
ayant les chutes de tension les plus élevées pour essayer de l’augmentation de tension que l’on constate entre les nœuds
voir les impacts que chacune des injections peut apporter 13 et 14 n’a pas lieu réellement, les deux n’étant pas sur le
sur le flux de puissance ainsi que le profil de tension du même départ et par conséquent pas sur la même ligne.
réseau.

(a) (b)
Figure 3. Résultats de simulation sans injection

Les résultats des premières simulations nous ont permis dire qu’on se retrouve avec des puissances « négatives ».
d’identifier quatre nœuds ayant des fortes chutes de Cela est remarqué lorsqu’on injecte une puissance de 100
tensions à savoir les nœuds 10, 11, 12 ,13. C’est sur quatre kW ou de 150 kW.
nœuds que nous procéderons aux simulations des
injections des différentes puissances. Les résultats pour une injection au nœud 11 sont présentés
par la figure 5. On constate comme pour le cas précèdent
Les résultats de la simulation (profil de tension et
que les profils de tension restent dans une marge
écoulement de puissance) pour une injection faite aux
acceptable, cependant ces caractéristiques ont un peu
nœuds 10 sont présentés à la figure 4.
régressé par rapport à celles ou l’injection est faite au nœud
On remarque que par rapport au scenario sans injection le 10 c’est-à-dire que les valeurs de tension se sont un peu
profil de tension (figure 4a) est largement amélioré car abaissée.
aucune valeur de tension n’est en dessous de 0.96 p.u pour Le constat et le même pour l’écoulement de puissance que
toutes les puissances injectées. En ce qui concerne dans le cas précèdent, plus la puissance injectée est
l’écoulement de puissance (b) on voit bien que plus la importante plus elle permet de décharger le poste source,
puissance injectée est grande, plus les puissances transitées ici pour une injection de 150 kW la puissance active
sur les lignes deviennent réduites, pour une puissance de transitée au bus 0 passe d’environ 270 kW à 100 kW
150 kW injectée par exemple, on passe de 0,275 MW à l’atténuation dans ce cas est moins importante que celle de
0,080 le réseau devient déchargé. Cependant on constate l’injection réalisée au nœud 10.
aussi que pour une puissance injectée supérieure à la
La figure 6 présente respectivement le profil de tension
puissance demandée en aval du point d’injection, il y a un
ainsi que l’écoulement de puissance pour chaque nœud du
changement de la direction du flux de puissance c'est-à-
réseau pour des injections faites au nœud 12.
5 International Journal of Applied Research and Technology, Vol (N°), pp……

On remarque encore une légère régression des profils de La figure 7 présente le profil de tension et l’écoulement de
tension par rapport aux cas précèdent ici ils tombent puissance pour chaque nœud du réseau pour des injections
jusqu’à 0.95 p.u au nœud 9, ils demeurent toutefois dans faites au nœud 13. On voit que les profils de tension
une marge acceptable. continuent à se détériorer, cette fois, on va même en
dessous de la barre de 0.95 p.u pour une injection de 150
En analysant les résultats présentés par la figure 6b le kW. En ce qui concerne le flux de puissance, il demeure
constat est tel qu’il se diffère des autres, en effet le positif jusqu’au 4ième nœud pour toutes les puissances
changement du sens du flux de puissance en amont du injectées il est aussi négatif dans les trois nœuds précédant
nœud où se fait l’injection que l’on voyait dès le premier le nœud 13 où les injections sont réalisées.
nœud a disparu, on se retrouve avec un écoulement de On constate aussi une très grande atténuation de la
puissance dans le sens positif jusqu’au deuxième nœud puissance transitée au bus 0 pour une puissance injectée de
pour une injection de 150 kW et jusqu’au sixième nœud 150 kW ce qui va se traduire par un soulagement du poste
pour une injection de 100 kW. source.

(a) (b)
Figure 4. Résultats de simulation sans injection 10

(a) (b)
Figure 5. Résultats de simulation sans injection 11
6 International Journal of Applied Research and Technology, Vol (N°), pp……

(a) (b)
Figure 6. Résultats de simulation sans injection 12

(a) (b)
Figure 7. Résultats de simulation sans injection 13

6. Conclusion et interprétation des résultats peut s’attendre à des déclenchements intempestifs de

celles-ci à cause de ce changement de sens [15].
Les résultats que nous avons vu dans la section précédente
nous ont montré que l’injection de l’énergie PV dans le Se référant à ces idées, dans cette étude on considère qu’une
réseau BT aux nœuds ayant les tensions les plus faibles injection de 60 kW au nœud 10 est la meilleure option car
permet une amélioration du profil de tension sur la ligne, au-delà de ne pas créer des changements de sens dans
cela est vérifiable car, des quatre nœuds où nous avons l’écoulement de puissance donne un bon profil de tension.
choisi de faire nos injections les profils de tension
correspondant étaient de loin meilleurs que celui tracé sans Références
injection on en déduit que cela est conforme à ce qui se
trouve dans la littérature. [1] K. Ilunga, «Partenariats avec les exploitants miniers
pour réduire le déficit énergetique entre l'offre et la
Toutefois au niveau des puissances transitées sur les demande. Katanga Mining Briefing. Lubumbashi:
différentes lignes du réseau, le changement de sens Présentation DPK/SNEL.,» SNEL, Lubumbashi, 2016.
constitue un problème sérieux car le réseau n’a pas été
[2] DPK/SNEL, «Rapport Annuel SNEL,» Lumbubashi,
conçu au départ pour gérer ce flux bidirectionnel, les
2013.
protections déjà pas très performantes se trouvant dans la
cabine peuvent se retrouver aveuglées ou inversement on
7 International Journal of Applied Research and Technology, Vol (N°), pp……

[3] W. B. Banza, I. Sapembi, K. Koj et M. Kasumba, decentralised energy systems,» Solar Energy, pp.
«Identification and analysis of alternative energy 150-161, 2018.
sources usedin domestic scale lubumbashi, dr [10] J. Von Appen et M. Braun, «Interdependencies
congo.,» International Journal of Engineering and between self-sufficiency preferences, techno-
Scientific Research, vol. 15, 2017. economi cdrivers for investment decisions and grid
[4] C. Holderman, J. Kissel et J. Beigel, «Distributed integration of residential PV storage systems,»
photovoltaic generation in Brazil: An economic Applied Energy, pp. 1140-1151, 2018.
viability analysis of small-scale photovoltaic systems [11] A. Colmenar-Santos, S. Campınez-Romero, C. Perez-
in the residential and commercial sectors,» Energy Molina et M. Castro-Gil, «Profitability analysis of grid-
Policy, 2013. connected photovoltaic facilities for household
[5] P. V. Gomes, N. K. Neto, L. Carvalho, S. J, J. T. Saraivaa, electricity self-sufficiency,» Energy Policy, pp. 749-
B. H. Diasd, V. Miranda et S. M. Souza, «Technical- 764, 2012.
economic analysis for the integration of PV systems [12] E. L. Ratnam, S. R. Weller et C. M. Kellett, «Scheduling
in Brazil considering policy and regulatory issues,» residential battery storage withsolar PV: Assessing
Energy Policy, pp. 199-206, 2018. the benefits of net metering,» Applied Energy, pp.
[6] J. Haroldo de Faria, F. B. M. Trigoso et J. A. M. 881-891, 2015.
Cavalcanti, «Review of distributed generation with [13] G. Brusco, A. Burgio, D. Menniti, A. Pinnarelli et N.
photovoltaic grid connected systems in Brazil: Sorrentino, «The economic viability of a feed-in tariff
Challenges and prospects,» Renewable and scheme that solely rewards self-consumption to
Sustainable Energy Reviews, 2016. promote the use of integrated photovoltaic battery
[7] C. Lacchini et R. Rüther, «The influence of systems,» Applied Energy, p. 1075–1085, 2016.
government strategies on the financial return of [14] B. W. Banza, Croissance urbaine et gestion du service
capital invested in PV systems located in different d'électricité de Lubumbashi: Etat de lieu, impact de
climatic zones in Brazil,» Renewable Energy, pp. 786- l'étalement urbain et modes de gouvernance,
798, 2015. Bruxelles, Lubumbashi: Thèse de doctorat ULB, 2017.
[8] G. de Martino Jannuzzi et C. A. de Melo, «Grid- [15] R. Caire, Impact de la production décentralisée BT sur
connected photovoltaic in Brazil: Policies and la HTA - Réflexion sur les impacts envisagés., Notice
potential impacts for 2030,» Energy for Sustainable technique IDEA - Gie , Mai 2001.
Development, p. 40–46, 2013.
[9] H. Schwarz, H. Schermeyera, V. Bertschb et W.
Fichtnera, «Self-consumption through power-to-heat
and storage for enhanced PV integration in