➢ Chamawa Larry

➢ Nemadeu Luc Jordan

Modélisation Dynamique et Analyse d'un Système Hybride PV-Éolien-Diesel pour

l'Alimentation d'une Infrastructure Hospitalière

Mots-clés : Résumé

Micro-réseau, Ce travail présente la conception et la simulation sous MATLAB/Simulink

Énergie Photovoltaïque, d'un système énergétique hybride destiné à secourir une charge hospitalière
Énergie Éolienne, critique. Le système intègre des sources renouvelables, notamment des
Groupe Électrogène, panneaux photovoltaïques (PV) et des turbines éoliennes, couplées à un
MATLAB/Simulink, groupe électrogène diesel pour garantir la continuité de service. Un bus
Bus DC. continu (DC) centralise la production renouvelable via des convertisseurs
Boost, tandis qu'un onduleur triphasé assure la conversion vers le bus
alternatif (AC) de 400 V. L'étude analyse les performances de chaque sous-
système face à des variations de charge dynamiques. Les résultats de
simulation démontrent la capacité des contrôleurs (MPPT pour le PV et PI
numérique pour l'éolien) à stabiliser les tensions de sortie, malgré
l'intermittence des ressources naturelles.

Figure 1: Architecture générale d'un micro-réseau hybride PV-Éolien-Diesel

1. Introduction
La sécurité énergétique des infrastructures critiques, telles que les centres solaire et de l'éolien avec la fiabilité des sources conventionnelles comme le
hospitaliers, constitue un défi majeur face à la croissance de la demande et à diesel.
la nécessité de réduire l'empreinte carbone. Les énergies renouvelables, bien Ce projet s'inscrit dans cette démarche en proposant un modèle dynamique
que compétitives, présentent une intermittence intrinsèque liée aux d'un système multi-sources. L'objectif est de maintenir une alimentation
conditions météorologiques. Pour pallier cette instabilité, l'hybridation des stable de 18 kW (puissance crête de l'hôpital) en optimisant l'usage des
sources est une solution viable, permettant de combiner la durabilité du ressources renouvelables. Le rapport détaille la modélisation de chaque
composant, les stratégies de contrôle implémentées et discute des défis
rencontrés lors de la gestion de puissance.

Figure 2: schéma logique montrant comment l'énergie circule des sources vers la charge.

Figure 3:Architecture globale de ton modèle Simulink

2. Données Météorologiques et Environnement de Simulation
La performance d'un système hybride dépend directement des conditions
climatiques. Pour cette étude, les données ont été extraites de la base de
données NASA POWER pour la journée du 2 Février 2024.
2.1. Analyse des gisements (NASA)
Les données récupérées présentent les caractéristiques suivantes :

• Irradiance Solaire (𝑮) : Mesurée en 𝑊/𝑚2 . On observe une

courbe en cloche typique, débutant à 06ℎ45 et s'achevant à
18ℎ15, avec un pic d'intensité à midi solaire.

2.2. Justification des valeurs constantes de simulation

Bien que les données réelles soient indispensables pour une étude annuelle,
nous avons opté pour des valeurs constantes lors des tests de validation
des sous-systèmes sous Simulink :
1. Stabilité de l'analyse : L'utilisation de constantes (𝐺 =
1000 𝑊/𝑚2 , 𝑇 = 25 °C pour le PV et 𝑣 = 12 𝑚/𝑠 pour
l'éolien) permet de vérifier la précision des régulateurs (MPPT et
PI) sans subir les fluctuations chaotiques du climat.
2. Régime établi : Cela permet de s'assurer que le bus DC atteint
bien sa valeur de consigne (400 𝑉) et de mesurer le temps de
• Vitesse du Vent (𝒗) : Mesurée à 10𝑚 de hauteur en 𝑚/𝑠. Le
réponse des convertisseurs de manière rigoureuse.
profil montre une variabilité importante, avec des rafales
intermittentes dépassant le seuil de démarrage de la turbine.
3. Modélisation du Sous-système Photovoltaïque
Cette section décrit la conversion de l'énergie solaire en énergie électrique.
3.1. Le Générateur PV
Le modèle mathématique utilisé est celui à une diode, défini par l'équation
du courant de sortie 𝐼𝑝𝑣 :
𝑉 + 𝑅𝑠 𝐼 𝑉 + 𝑅𝑠 𝐼
𝐼 = 𝐼𝑝ℎ − 𝐼0 [exp ( ) − 1] −
𝑛𝑉𝑡 𝑅𝑠ℎ
Où 𝐼𝑝ℎ est le courant photogénéré (A), 𝑅𝑠 la résistance série (Ω) et 𝑅𝑠ℎ la
résistance shunt (Ω).
3.2. Convertisseur Boost et MPPT
Pour extraire la puissance maximale (𝑃𝑚𝑎𝑥 ), un algorithme Perturb and
Observe (P&O) a été implémenté. Il ajuste le rapport cyclique (𝐷) du
convertisseur Boost selon la logique suivante :

• Si Δ𝑃 / Δ𝑉 > 0, alors on augmente 𝑉.

• Température (𝑻) : Mesurée à 2𝑚 en °C. Elle influence • Si Δ𝑃 / Δ𝑉 < 0, alors on diminue 𝑉.
directement le rendement des cellules photovoltaïques (une
hausse de température diminue la tension en circuit ouvert 𝑉𝑜𝑐 ).
 Figure 4: Bloc PV connecté au Boost piloté par le MPPT

3.3. Dimensionnement du convertisseur Boost PV

Le convertisseur élévateur (Boost) est dimensionné pour assurer le passage

de la tension du champ PV (souvent variable selon l'ensoleillement) vers la
tension régulée du bus DC de 400 𝑉. Les paramètres critiques sont
calculés selon les formules suivantes :

• Rapport cyclique (𝑫) :

𝑉𝑖𝑛
𝐷 =1−
𝑉𝑜𝑢𝑡

• Inductance minimale (𝑳𝒎𝒊𝒏 ) : calculée pour garantir un mode

de conduction continue (CCM) :

𝑉𝑖𝑛 ⋅ 𝐷
𝐿𝑚𝑖𝑛 =
𝑓𝑠 ⋅ Δ𝐼𝐿

Où 𝑓𝑠 est la fréquence de découpage (10 𝑘𝐻𝑧) et Δ𝐼𝐿 l'ondulation du

courant. 1. Évolution de la tension du Bus DC (𝑽𝒅𝒄 )

• Capacité de sortie (𝑪𝒐𝒖𝒕) : dimensionnée pour limiter • Observation : À l'instant 𝑡 = 0 𝑠, la tension amorce une
l'ondulation de la tension sur le bus DC : montée rapide. On observe un léger dépassement (overshoot)
avant que le signal ne se stabilise précisément à 400 𝑉 aux
𝐼𝑜𝑢𝑡 ⋅ 𝐷
𝐶𝑜𝑢𝑡 = alentours de 𝑡 = 0.2 𝑠.
𝑓𝑠 ⋅ Δ𝑉𝑜𝑢𝑡

3.4. Résultats de simulation et interprétations • Interprétation : Ce délai de 0.2 𝑠 correspond au temps de

réponse du régulateur de tension. La stabilité maintenue à
Les courbes obtenues sous l'environnement Simulink permettent de valider 400 𝑉 après cette phase transitoire démontre l'efficacité du
le comportement dynamique du système PV-Boost. contrôle en boucle fermée pour rejeter les perturbations.

2. Suivi du Point de Puissance Maximale (MPPT)

• Observation : La puissance extraite du champ PV atteint sa

valeur nominale (𝑃𝑝𝑣 ≈ 10 𝑘𝑊 pour les conditions de test) de
manière stable. On remarque de très faibles oscillations autour
du point de fonctionnement.

• Interprétation : Ces oscillations sont caractéristiques de

l'algorithme "Perturb and Observe" (P&O) qui cherche
continuellement le sommet de la courbe de puissance. La faible
amplitude des oscillations confirme que le pas de perturbation
choisi est optimal pour minimiser les pertes énergétiques.
3. Courant d'inductance du Boost • Interprétation : Ce temps de réponse plus long est dû à l'inertie
mécanique de la turbine simulée et aux gains du correcteur
• Observation : Le courant présente une forme en dents de scie numérique (𝐾𝑝 = 0.01, 𝐾𝑖 = 0.5) choisis pour éviter les
stabilisée, sans jamais toucher la valeur zéro. oscillations sur le bus DC.

• Interprétation : Cela confirme que le convertisseur fonctionne 2. Comportement du correcteur numérique

bien en Mode de Conduction Continue (CCM), ce qui réduit
les contraintes sur les composants semi-conducteurs et assure un • Observation : On constate que le rapport cyclique 𝐷 se stabilise
transfert d'énergie plus fluide vers l'hôpital. autour d'une valeur moyenne sans sauts brusques, malgré le
bruit de mesure.
Valeurs des composants utilisés (PV Boost)
• Interprétation : L'utilisation de la fonction MATLAB permet
• Inductance (𝑳) : 5 𝑚𝐻
une saturation précise du signal (0.05 < 𝐷 < 0.95),
protégeant ainsi les composants de puissance contre les
• Capacité (𝑪) : 1000 𝜇𝐹 commandes aberrantes.

• Fréquence (𝒇𝒔) : 10 𝑘𝐻𝑧

4. Modélisation du Sous-système Éolien

Le système éolien convertit l'énergie cinétique du vent en énergie

électrique. Il se compose d'une turbine, d'une machine synchrone à aimants
permanents (PMSG) et d'un convertisseur Boost contrôlé numériquement.

4.1. Modèle de la turbine éolienne

La puissance mécanique 𝑃𝑚 extraite du vent est régie par l'équation suivante

:

1
𝑃𝑚 = 𝜌𝐴𝐶𝑝 (𝜆, 𝛽)𝑣 3
2

• 𝜌 : Masse volumique de l'air (1.225 𝑘𝑔/𝑚3).

• 𝐴 : Surface balayée par les pales (𝑚2 ).

• 𝑣 : Vitesse du vent (𝑚/𝑠).

• 𝐶𝑝 : Coefficient de puissance, représentant le rendement

aérodynamique de la turbine (limite de Betz).

4.2. Dimensionnement du Boost Éolien

Comme pour le PV, le Boost éolien adapte la tension redressée de la

génératrice vers les 400 𝑉 du bus DC.

• Inductance (𝑳𝒘𝒊𝒏𝒅 ) : 4 𝑚𝐻 (choisie pour lisser le courant

redressé).

• Capacité (𝑪𝒘𝒊𝒏𝒅 ) : 1500 𝜇𝐹 (plus élevée que le PV pour

compenser les fluctuations mécaniques de la turbine).

4.3. Contrôle par Correcteur Numérique (MATLAB Function)

L'innovation de ce sous-système réside dans l'utilisation d'un correcteur

Proportionnel-Intégral (PI) numérique implémenté via une fonction
MATLAB. Contrairement aux blocs classiques, ce code simule un
microcontrôleur traitant l'erreur de tension à chaque pas de temps 𝑇𝑠 .

L'algorithme calcule le rapport cyclique 𝐷(𝑘) selon l'équation de

récurrence :

𝐷(𝑘) = 𝐷(𝑘 − 1) + 𝐾𝑝 [𝑒(𝑘) − 𝑒(𝑘 − 1)] + 𝐾𝑖 𝑇𝑠 𝑒(𝑘)

Où 𝑒(𝑘) = 𝑉𝑟𝑒𝑓 − 𝑉𝑚𝑒𝑠 .

4.4. Résultats de simulation et interprétations

Les tests ont été effectués avec une vitesse de vent constante de 12 𝑚/𝑠 et
un angle de calage des pales 𝛽 = 45∘ .

1. Évolution de la tension de sortie

• Observation : La tension monte de manière plus amortie que

celle du PV. Elle atteint le régime permanent de 400 𝑉 en
environ 0.35 𝑠.
 Figure 5:Le bloc "Wind Turbine" connecté à la PMSG et au redresseur

5. Modélisation du Groupe Électrogène Diesel paramétré pour fournir la tension nominale requise par le bus
AC de l'hôpital.
Le groupe électrogène intervient comme source de secours. Il est composé
d'un moteur à combustion interne (Diesel) couplé à une machine synchrone 2. Réponse à la charge
à pôles saillants.
• Observation : Lors de l'enclenchement du breaker, on observe
5.1. Le Moteur Diesel et le Gouverneur une chute transitoire très brève de la vitesse 𝜔 avant un retour à
1 𝑝𝑢.
Le moteur Diesel convertit l'énergie chimique du carburant en énergie
mécanique. Le gouverneur (régulateur de vitesse) ajuste le débit de
• Interprétation : Le gouverneur réagit efficacement à l'appel de
carburant pour maintenir une fréquence constante de 50 𝐻𝑧 (1 𝑝𝑢).
charge en augmentant le couple mécanique, garantissant ainsi la
L'équation simplifiée du couple mécanique 𝑇𝑚 produit est : stabilité de la fréquence du micro-réseau.

𝐾𝑑
𝑇𝑚 (𝑠) = ⋅ 𝜙(𝑠)
1 + 𝜏𝑑 𝑠
Instructions pour tes captures d'écran (Section 5) :
• 𝐾𝑑 : Gain du moteur.
• Capture Simulink 3 : Le bloc "Diesel Engine" et le bloc
"Synchronous Machine".
• 𝜏𝑑 : Constante de temps de combustion.
• Courbe de Simulation 5 : Les trois tensions sinusoïdales
• 𝜙(𝑠) : Sortie du régulateur de vitesse.
(Vabc) du Diesel.
5.2. La Machine Synchrone (Alternateur)
• Courbe de Simulation 6 : La vitesse de rotation 𝜔 du groupe
La machine synchrone transforme le couple 𝑇𝑚 en puissance électrique. La (stabilisée à 1 pu).
tension de sortie dépend de la vitesse de rotation et du courant d'excitation
(𝐼𝑓 ).
6. Conversion AC et Interconnexion (L'Onduleur)
L'équation de la force électromotrice (FEM) induite $E$ est :
L'onduleur (VSC - Voltage Source Converter) est l'interface entre le bus DC
𝐸 = 𝐾 ⋅Φ⋅𝜔
stabilisé à 400 𝑉 et le bus AC de l'hôpital.
• 𝐾 : Constante de construction de la machine. 6.1. Modulation de Largeur d'Impulsion (PWM)

• Φ : Flux magnétique (lié à l'excitation 𝑉𝑓 ). Le pilotage des transistors (IGBT) utilise une stratégie de modulation
sinusoïdale asservie à la PLL (Phase Locked Loop) pour garantir le
• 𝜔 : Vitesse angulaire (𝑟𝑎𝑑/𝑠). synchronisme avec le Diesel ou le Réseau.

5.3. Résultats de simulation et interprétations L'indice de modulation 𝑚 est défini par :

Le Diesel a été configuré pour fournir une tension nominale de 400 V 𝑉𝑐𝑜𝑛𝑡𝑟𝑜𝑙
𝑚=
(RMS entre phases). 𝑉𝑡𝑟𝑖

1. Analyse de la tension de sortie 6.2. Filtre de sortie LC

• Observation : Sur le Scope, l'amplitude de la tension phase- Pour éliminer les harmoniques de découpage et obtenir une onde de tension
terre atteint 326 V (ce qui correspond bien à 400√2/√3). Le pure, un filtre LC est placé en sortie. Sa fréquence de coupure 𝑓𝑐 est
signal est parfaitement sinusoïdal à 50 𝐻𝑧. calculée par :

1
• Interprétation : La stabilisation de l'amplitude à 326 V crête 𝑓𝑐 =
2𝜋√𝐿𝑓 𝐶𝑓
confirme que le système d'excitation (AVR) est correctement
 contrôle supérieure. Le système fonctionne en mode "somme de
puissances" sans arbitrage économique ou prioritaire complexe.

3. Synchronisation Diesel-Onduleur : Caler précisément

l'excitation du diesel pour obtenir exactement la même
amplitude que l'onduleur (326 𝑉 crête) a nécessité de multiples
itérations sur les gains du régulateur de tension (AVR).

7.3. Suggestions d'amélioration

Pour transformer ce prototype en un système prêt pour une installation

hospitalière réelle, les pistes suivantes sont suggérées :

• Implémentation d'un stockage batterie : L'ajout d'un banc de

batteries permettrait de lisser l'intermittence des données NASA
(nuages, chutes de vent) sans solliciter immédiatement le groupe
Diesel.

• Finalisation de l'EMS : Développer un organigramme de

décision complet (Stateflow) pour gérer les priorités (Solaire >
Vent > Batterie > Réseau > Diesel).

• Analyse du THD : Effectuer une analyse harmonique

approfondie pour garantir un taux de distorsion inférieur à 5%,
conformément aux normes IEEE 519.

7. Intégration Globale, Analyse des Performances et Défis

Cette section présente les résultats du système complet face au profil de

charge réel de l'hôpital et discute des difficultés de mise en œuvre.

7.1. Comportement face au profil de charge dynamique (𝑷𝒍𝒐𝒂𝒅 )

La charge de l'hôpital a été simulée via le bloc Three-Phase Dynamic Load,

piloté par les données exportées du Workspace (profil 𝑃𝑙𝑜𝑎𝑑 ).

1. Suivi de la puissance consommée

• Observation : La puissance active appelée par la charge varie

entre 5 𝑘𝑊 et 18 𝑘𝑊 selon les paliers définis. On observe que
l'onduleur et les sources de secours (Diesel/Réseau) s'adaptent
pour fournir cette puissance.

• Interprétation : Le système répond à la demande, mais la

stabilité de la tension du bus AC dépend fortement de la rapidité
de réponse des sources d'appoint lors des échelons de charge.

2. Stabilité de la tension AC

• Observation : Malgré les variations brusques de $𝑃_{𝑙𝑜𝑎𝑑}$, la

tension du bus AC reste centrée sur une amplitude de 326 V
Figure 6 : Superposition de 𝑃𝑔𝑒𝑛é𝑟é𝑒 (PV+Vent) et 𝑃𝑑𝑒𝑚𝑎𝑛𝑑é𝑒
(crête). Des micro-creux de tension (sag) sont visibles
uniquement lors des commutations de breakers. (Pload)

Organigramme du Système de Gestion Énergétique (EMS)

• Interprétation : L'utilisation de la PLL (Phase Locked Loop)
est ici validée : elle permet à l'onduleur de rester synchronisé, 1. DÉBUT : Mesure des variables
évitant ainsi des déphasages qui seraient destructeurs pour les
équipements médicaux. o Lire 𝑃𝑝𝑣 (Solaire)

7.2. Difficultés rencontrées et problèmes non résolus o Lire 𝑃𝑤𝑖𝑛𝑑 (Éolien)

Le développement de ce micro-réseau hybride a fait face à plusieurs défis o Lire 𝑃𝑙𝑜𝑎𝑑 (Profil de l'hôpital)
techniques :
o Calculer 𝑃𝑒𝑛𝑟 = 𝑃𝑝𝑣 + 𝑃𝑤𝑖𝑛𝑑
1. Oscillations du Bus DC (Phénomène à 100 Hz) : L'une des
difficultés majeures fut l'apparition d'une ondulation sinusoïdale 2. TEST 1 : Les énergies renouvelables suffisent-elles ?
sur le bus DC lors de l'augmentation de la charge AC. Ce
problème a été atténué en augmentant la capacité de lissage o Si 𝑷𝒆𝒏𝒓 ≥ 𝑷𝒍𝒐𝒂𝒅 :
(𝐶𝑑𝑐 ) à 4700 𝜇𝐹, bien que cela ralentisse la simulation.
▪ Action : Maintenir le Breaker Diesel
2. Gestion de la puissance (EMS) - Problème non résolu : Par OUVERT (OFF).
manque de temps, l'algorithme de gestion intelligente des flux
(Energy Management System) n'a pas été totalement finalisé. ▪ Action : Maintenir le Breaker Réseau
Actuellement, la répartition de puissance entre le PV, l'Éolien et OUVERT (OFF).
le Diesel n'est pas optimisée en temps réel par une boucle de
 ▪ Optionnel (si surplus) : Envoyer • f (Fréquence nominale) : 50 𝐻𝑧
l'excédent vers le réseau ou les batteries.
• fs (Fréquence de découpage) : 10 000 𝐻𝑧
▪ Retour au début.
• G (Irradiance test) : 1000 𝑊/𝑚2
3. TEST 2 : Si 𝑷𝒆𝒏𝒓 < 𝑷𝒍𝒐𝒂𝒅 (Déficit énergétique)

o Action : Calculer le déficit 𝑃𝑔𝑎𝑝 = 𝑃𝑙𝑜𝑎𝑑 − 𝑃𝑒𝑛𝑟 .

• Kp (Gain Proportionnel) : 0.01 (Éolien)

o Vérifier la disponibilité du Réseau Public. • Ki (Gain Intégral) : 0.5 (Éolien)

o Si Réseau disponible : • L (Inductance Boost) : 5 × 10−3 𝐻

▪ Action : Fermer le Breaker Réseau (ON). • Pload (max) : 18 000 𝑊

▪ Action : Maintenir le Breaker Diesel • Ts (Pas de calcul) : 1 × 10−4 𝑠

OUVERT.
• Vdc (Consigne) : 400 𝑉
▪ L'onduleur reste en parallèle pour aider
le réseau. • Vrms (L-L) : 400 𝑉

o Si Réseau indisponible (Coupure/Black-out) :

▪ Action : Envoyer commande de

démarrage au Moteur Diesel.

▪ Action : Attendre stabilisation (𝑡 >

0.5𝑠).

▪ Action : Fermer le Breaker Diesel (ON).

4. FIN de boucle : Surveillance continue

8. Conclusion, Nomenclature et Constantes

8.1. Conclusion

Ce projet a permis de modéliser avec succès un micro-réseau hybride

complexe sous MATLAB/Simulink. Malgré les défis liés à la stabilisation
du bus DC et à la synchronisation des sources AC, les simulations
confirment qu'une combinaison intelligente de PV, d'éolien et de diesel peut
répondre aux exigences strictes d'un hôpital. L'utilisation d'un correcteur PI
numérique via une MATLAB Function a démontré la faisabilité d'un
contrôle précis et flexible, ouvrant la voie à une implémentation sur
microcontrôleur réel.

8.2. Nomenclature (Unités du Système International)

Symbole Grandeur Physique Unité (SI)

$V$ Tension électrique Volt (V)

$I$ Courant électrique Ampère (A)

$P$ Puissance active Watt (W)

$Q$ Puissance réactive Volt-Ampère Réactif (VAR)

$f$ Fréquence Hertz (Hz)

$G$ Irradiance solaire Watt par mètre carré (𝑊/𝑚2 )

$v$ Vitesse du vent Mètre par seconde (𝑚/𝑠)

$\omega$ Vitesse angulaire Radian par seconde (𝑟𝑎𝑑/𝑠)

$L$ Inductance Henry (H)

$C$ Capacité Farad (F)

$T$ Température Kelvin (K) ou Degré Celsius (°C)

8.3. Valeurs des Constantes et Paramètres

• C (Capacité Bus DC) : 4700 × 10−6 𝐹