Optimisation de L'irradiation Solaire Sur Le Plan Des Capteurs Photovoltaïques À Structure Porteuse Statique
Optimisation de L'irradiation Solaire Sur Le Plan Des Capteurs Photovoltaïques À Structure Porteuse Statique
MEMOIRE DE MASTER
DOMAINE : SCIENCES ET TECHNOLOGIE
FILIERE : ELECTRONIQUE
OPTION : TECHNOLOGIE DES SYSTEMES ELECTR-ENERGETIQUES DE
SOURCES D’ENERGIES RENOUVELABLES
THEME
N° D’ordre:2016/S3ER05/89/385
PROMOTION: JUIN 2016
SOMMAIRE
Sommaire
Introduction Générale
2.1 Introduction................................................................................................................................ 51
... Irradiation globale sur un plan incline ..................................................................................... 52
2.2.1 Rayonnement extraterrestre . .................................................................................... 53
2.2.2 Irradiation sur un plan horizontal ............................................................................ 54
. a) Fraction diffuse ................................................................................................ 55
b) Indice de clarté ................................................................................................. 56
iii
SOMMAIRE
2.3 Estimation des moyennes mensuelles de l’irradiation globale journalière.......................57
2.3.1 Irradiation journalière globale ........................................................................... 58
2.4 Estimation des moyennes mensuelles des irradiations journalières ............................... 59
2.4.1 Modèle de Liu et Jordan. .................................................................................... 60
iv
SOMMAIRE
3.2.10 Irradiation solaire de la ville de Tindouf ............................................................ 86
3.2.11 Irradiation solaire de la ville de Skikda ............................................................... 87
3.2.12 Irradiation solaire de la ville d'Adrar ................................................................... 88
3.2.13 Irradiation solaire de la ville de Sétif ................................................................... 89
3.2.14 Irradiation solaire de la ville d’El-tarf ................................................................ 90
3.2.15 Irradiation solaire de la ville de Bejaia ................................................................ 91
3.2.16 Irradiation solaire de la ville de Ghardaïa ...................................................... 92
3.2.17 Irradiation solaire de la ville de Boumardes .................................................. 93
3.2.18 Irradiation solaire de la ville d’Illizi .............................................................. 94
3.2.19 Irradiation solaire de la ville de Tlemcen ...................................................... 95
3.2.20 Irradiation solaire de la ville de Mila ............................................................. 96
3.3 Pertes d’irradiation solaire annuelle ................................................................................ 97
3.3.1 Pertes d’irradiation solaire annuelle de la ville de M'sila ................................. 98
3.3.2 Pertes d’irradiation solaire annuelle de la ville d’Ouargla ................................ 99
3.3.3 Pertes d’irradiation solaire annuelle de la ville d’Adrar .................................. 100
v
Listes des figures et tableaux
vi
Listes des figures et tableaux
vii
Introduction générale
L'énergie solaire occupe une des places les plus importantes parmi les diverses sources
d'énergies possibles. Une connaissance précise de distribution du rayonnement solaire à un
endroit géographique donné est d'une importance majeure pour le développement des
dispositifs d'énergie solaires. Malheureusement, pour beaucoup de pays en voie de
développement les mesures du rayonnement solaire ne sont pas facilement disponibles.
Donc, il est important d'élaborer des modèles d'estimation du rayonnement sur la base des
données astronomiques et météorologiques. Au cours des années, beaucoup des modèles ont
étés proposés pour prévoir la valeur du rayonnement solaire en utilisant plusieurs
paramètres, à savoir : le rayonnement extraterrestre, la durée d'insolation, l'altitude, la
latitude, la longitude et l’albédo. Notons également que les composantes du rayonnement
solaire direct et diffus sont mesurées essentiellement sur un plan horizontal. Or, pour capter
maximum, les capteurs solaires que ce soient photovoltaïque ou thermique doivent être
minutieusement orientés et inclinés.
1
radiation solaire globale sur une surface inclinée et nous donnerons une présentation de
modèle utilisés dans notre étude pour l'optimisation de l'irradiation solaire globale.
Le dernier chapitre est réservé aux résultats obtenus par l’application du modèle de Perez
pour déterminer les angles optimaux de l’orientation et l’inclinaison (αopt et βopt)des capteurs
PV installés dans différentes sites algérien.
2
Chapitre 1
Rayonnement solaire
1.1. Introduction
Plusieurs paramètres peuvent influencer l‘énergie captée par un système photovoltaïque. le
mouvement apparent de la terre et son influence sur la définition du jour et du temps, la
position du soleil, l‘orientation et l’inclinaison du capteur, la situation géographique,
quelques grandeurs qui sont en relation avec l'évaluation de l'irradiation solaire.
Dans ce chapitre on va parler sur L'énergie solaireetprésenterle rayonnement solaire hors
atmosphère et le rayonnement solaire au niveau du sol .Et aussi nous présentons les
paramètres essentiels pour l‘estimation des énergies captées par les systèmes PV.
3
CHAPITRE 1 RAYONNEMENT SOLAIRE
4
CHAPITRE 1 RAYONNEMENT SOLAIRE
1.4.1. Troposphère
C'est la couche atmosphérique la plus proche du sol, dans laquelle nous vivons. Son
épaisseur varie de 7 kilomètres de hauteur au-dessus des pôles, 18 kilomètres au-dessus de
l'équateur et environ 13 kilomètres selon les saisons, dans la zone tempérée. La température
diminue à mesure qu'on s'élève dans la troposphère, c'est-à-dire à mesure qu'on s'éloigne du
sol, atteignant -56C° à la tropopause (zone séparant la troposphère de la stratosphère); l’air
près du sol est plus chaud qu'en altitude car la surface le réchauffe.
5
CHAPITRE 1 RAYONNEMENT SOLAIRE
1.4.2. Stratosphère
C'est la couche atmosphérique au-dessus de la troposphère, dans la stratosphère on trouve la
couche d'ozone. Sa formation il y a 3 milliards d’années, a permis le développement de la
vie sur notre planète. Car elle filtre une partie des rayons ultraviolets du soleil qui sont
extrêmement nocifs pour tout être vivant. Cette absorption provoque un dégagement
d'énergie sous forme de chaleur. C'est pourquoi la température augmente lorsqu'on s'élève
dans la stratosphère. Les équations chimiques suivantes décrivent la formation de la
molécule d’ozone.
1.4.3. Mésosphère
Elle se trouve juste au-dessus de la stratosphère. Dans cette couche, la température
recommence à décroître avec l'altitude pour atteindre un minimum de -80 degrés Celsius à
une altitude d'environ 80 kilomètres. Les poussières et particules qui proviennent de l'espace
(les météores) s'enflamment lorsqu'elles entrent dans la mésosphère à cause de la friction de
l'air. Ce phénomène nous apparaît sous la forme « d'étoiles filantes ».
1.4.4. Thermosphère
C’est la couche la plus haute. Dans cette couche, la température augmente avec l'altitude et
peut atteindre environ 100 degrés Celsius. La thermosphère atteint des milliers de kilomètres
d'altitude et disparaît graduellement dans l'espace. La partie inférieure de la thermosphère
réfléchit les ondes courtes (ondes radio), elle s'appelle l'ionosphère.
6
CHAPITRE 1 RAYONNEMENT SOLAIRE
𝐡·𝐂
𝐄=𝐡·𝐯= (1.1)
𝛌
Où
h :est la constante de Planck.
𝑣 : est la fréquence.
C : est la vitesse de la lumière.
Le spectre du rayonnement extraterrestre correspond environ à l’émission d’un corps noir
porté à 5800°K. Une courbe standard, compilée selon les données recueillies par les
satellites, est désignée sous le nom de AM0. Sa distribution en énergie est répartie en :
ultraviolet UV 0.20 < l < 0.38 mm 6.4%.
visible 0.38 < l < 0.78 mm 48.0%.
7
CHAPITRE 1 RAYONNEMENT SOLAIRE
Plus le soleil est bas sur l’horizon, plus il va traverser une épaisseur importante
d’atmosphère et plus il va subir de transformation [10].
On appelle <masse d’air>,(m) la masse d’atmosphère traversée par le rayonnement direct
pour atteindre le sol (par ciel clair), par rapport à une traversée verticale au niveau de lamer.
Cela dépend surtout de la hauteur angulaire du soleil.
8
CHAPITRE 1 RAYONNEMENT SOLAIRE
Pour le soleil au zénith (ℎ =90°), m* =1 (AM1 = Air Mass One en anglais) et pour le soleil
vu sous un angle (ℎ =30°), m* =2(AM2). Le cas m* =0 correspond à l’étude du
rayonnement solaire à la limite de l’atmosphère terrestre (référence AM0).
𝐂𝐬
𝐇𝐨 = [𝐜𝐨𝐬(ф) 𝐬𝐢𝐧(𝛚) + 𝛚𝐬𝐬 𝐬𝐢𝐧(𝛅) 𝐬𝐢𝐧(ф)] (1.4)
𝛑
Où
Nj: Nombre de jours dans l’année.
CS0: Constante solaire. (Elle est égale à 1367W/m2).
ф: Latitude.
δ: Déclinaison solaire.
9
CHAPITRE 1 RAYONNEMENT SOLAIRE
Au sol le rayonnement global reçu par une surface plane d’inclinaison quelconque est
constitué de plusieurs composantes principales [14] ;[15] :
10
CHAPITRE 1 RAYONNEMENT SOLAIRE
On en distingue :
Le rayonnement diffus isotropique provenant de toute la voute céleste.
Le rayonnement diffus circonférence provenant de la couronne solaire.
Le rayonnement diffus provenant de l’horizon.
11
CHAPITRE 1 RAYONNEMENT SOLAIRE
Après son voyage dans l’espace le rayonnement solaire atteint l’atmosphère, une certaine
quantité de ce rayonnement est absorbée par l’atmosphère, les nuages et les particules en
suspension dans l’atmosphère, une autre partie est réfléchie dans l’espace et une certaine
quantité est absorbée par la surface de la terre. La surface de la terre renvoie aussi une
certaine partie sous forme d’un rayonnement de différentes longueurs d’ondes en raison de
la température superficielle de la terre. Environ 70% du rayonnement solaire total reste dans
la surface de la terre et dans l’atmosphère.
La rotation de la terre autour de son axe produit des variations horaires dans l’intensité du
rayonnement reçue par un emplacement donné pendant le jour et aboutit même à l’extinction
complète pendant la nuit. En passant par l’atmosphère, le rayonnement solaire est
partiellement dispersé dans toutes les directions, partiellement absorbé et partiellement
réfléchie, par les gaz atmosphériques, les aérosols et les gouttelettes de nuage. Le
rayonnement dispersée atteignant la surface de la terre constitue ce qu’on appelle le
rayonnement diffus. À nouveau, une partie de ce rayonnement diffus revient en arrière dans
l’espace, alors qu’une partie atteinte la surface de la terre sous forme d’un rayonnement,
direct, diffus et réfléchi (figure 1.5).
Le rayonnement direct comme son nom l’indique est le rayonnement solaire reçu à
n’importe quel point sur terre directement du soleil sans aucune perturbation. Le
rayonnement réfléchi est principalement renvoyé par la surface réceptrice, cette composante
est plus importante dans des régions montagneuses [16].
1.11.1. Pyranomètre
C'est le capteur le plus fréquemment utilisé pour mesurer le rayonnement global ou diffus. Il
se compose d'une thermopile composé d'une série de thermocouples avec alternance des
soudures chaudes et soudures froides [17].
Les soudures chaudes sont en contact thermique avec la surface absorbante
une «noire» avec une capacité thermique faible.
12
CHAPITRE 1 RAYONNEMENT SOLAIRE
Les soudures froides sont en contact avec un matériau ayant une capacité
thermique relativement élevée et reste à la température ambiante.
Lorsque le rayonnement solaire tombe sur la surface absorbante (noire) la
quasi-totalité du rayonnement incident est absorbé.
1.11.3. Pyrhéliomètre
C’est est un instrument destiné à la mesure de l'intensité du rayonnement solaire direct. Il est
généralement constitué par une thermopile peinte en noir, donc très absorbante, montée sur
un système de suivi. Ainsi cet instrument ne mesure que le rayonnement directement issu du
soleil.
1.11.4. Spectro-radiométre
13
CHAPITRE 1 RAYONNEMENT SOLAIRE
l’éclairement énergétique.
1.11.5. Albedométre
Il mesure l’albédo qui est définie par le rapport du rayonnement réfléchi au rayonnement
global on utilise deux pyranomètres opposés à l’horizontal.
14
CHAPITRE 1 RAYONNEMENT SOLAIRE
a) Hauteur du soleil
C’est l’angle formé par la direction du soleil et sa projection sur le plan horizontal. ℎ = 0 :
correspond au lever et au coucher du soleil, la hauteur du soleil varie entre +90 et −90.Il est
donné par la relation suivante :
b) Zénith
C’est le point le plus élevé de l'hémisphère céleste, se trouvant directement à la verticale de
l'observateur. Il est donné par l’équation suivante :
Alors
𝐬𝐢𝐧(𝜽𝒛) = 𝐬𝐢𝐧(𝒉) (1.7)
c) Azimut du soleil
C'est l'angle que fait sur le plan horizontal et la projection de la direction du soleil avec la
direction du sud.
L’azimut du soleil varie à chaque instant de la journée .α : varie entre -180 et +180.
On prend : α>0 vers l’ouest.
α<0 vers l’est.
α =0 vers le sud
15
CHAPITRE 1 RAYONNEMENT SOLAIRE
𝐜𝐨𝐬(𝛅) . 𝐬𝐢𝐧(𝛚)
𝐬𝐢𝐧(𝛂) = (1.8)
𝐜𝐨𝐬(𝐡)
a)Déclinaison(δ)
C’est l’angle entre la direction terre-soleil et le plan équatorial. La déclinaison varie de (-
23°45’) au solstice d’hiver à (+23°45’) au solstice d’été et elle est nulle aux équinoxes. Elle
est bien représentée par la formule suivante :
𝟐𝛑
𝛅(°) = 𝟐𝟑. 𝟒𝟓. 𝐬𝐢𝐧 (( ) . (𝟐𝟖𝟒 + 𝐍𝐣)) (1.9)
𝟑𝟔𝟓
16
CHAPITRE 1 RAYONNEMENT SOLAIRE
TSV : est le temps solaire vrai (TSV en heures et ω en degrés), tel que le soleil se trouve au
zénith à midi. A chaque heure qui s'écoule correspond une augmentation de l'angle horaire
de 15°.
ωss=-ωsr=cos⁻¹(-tan(Ф)×tan(δ)) (1.11)
Où
δ : est la déclinaison
Ф: est la latitude du lieu
17
CHAPITRE 1 RAYONNEMENT SOLAIRE
a)Latitude du lieu Ф:
L’angle de latitude est l’angle entre une ligne tracée d’un point sur la surface de la terre avec
le centre de cette dernière, et le cercle équatorial de la terre. L’intersection du cercle
équatorial avec la surface de la terre forme l’équateur, est indiquée en tant que latitude de 0°,
le pôle nord par latitude +90° et le pôle sud par la latitude -90°.Ф varie de -90° à +90°, tel
que :Ф> 0 vers le nord et Ф< 0 vers le sud
b)Longitude λ :
Permet de localiser un point à l’Est ou à l’Ouest d’une ligne Nord- Sud de référence appelée
le méridien Greenwich. λ varie de -180° à +180°.
tel que : λ > 0 à l’est du méridien de Greenwich et λ < 0 à l’ouest du méridien de Greenwich.
18
CHAPITRE 1 RAYONNEMENT SOLAIRE
19
CHAPITRE 1 RAYONNEMENT SOLAIRE
en coordonnées horizontale du lieu, cet angle dépend des paramètres de la position du soleil
et de l’angle d’inclinaison du capteur. Il est défini par :
𝒄𝒐𝒔(𝜽𝒊) = 𝒔𝒊𝒏(𝜹) . 𝒔𝒊𝒏(Ф) 𝒔𝒊𝒏(𝜷) − 𝒔𝒊𝒏(𝜹) 𝒄𝒐𝒔(Ф) 𝒔𝒊𝒏(𝜷) 𝒄𝒐𝒔(𝜶) 𝒄𝒐𝒔(Ф) 𝒄𝒐𝒔(𝜷) 𝒄𝒐𝒔(𝝎)
+ 𝒄𝒐𝒔(𝜹) 𝒔𝒊𝒏(Ф) 𝒔𝒊𝒏(𝜷) 𝒄𝒐𝒔(𝜶) 𝒄𝒐𝒔(𝝎) + 𝒄𝒐𝒔(𝜹) 𝒔𝒊𝒏(𝜶) 𝒔𝒊𝒏(𝜷) 𝒔𝒊𝒏(𝝎)
(1.12)
CONCLUSION
L'irradiation solaire globale pour un site donné est fonction de plusieurs facteurs;
astronomiques et météorologiques. Dans ce chapitre on a parlé de ces facteurs
astronomiques et de leurs définitions. Elle est influencée par deux principaux paramètres, en
premier lieu on trouve les constituants de l'atmosphère, qu'ils soient permanents ou ceux
dont la composition varie largement avec les activités de l'homme, et en second c'est le
mouvement de la terre, alors il important de faire une étude préliminaire du climat pour un
site donné avant tout expose d'exploitation de l'énergie solaire.
20
Chapitre 2
Modélisation du rayonnement
solaire
2.1. Introduction
Beaucoup de modèles pour estimer le rayonnement solaire global, diffus et direct à partir de la
durée d’insolation des valeurs moyennes mensuelles des valeurs moyennes journalières ou
horaires d’irradiation moyennant. Dans ce chapitre, sont présentés les différents modèles
d’évaluation de l’éclairement en un site donnée. Ces modèles sont basé sur les différentes
contributions du ciel au rayonnement diffus, à savoir l’anisotropie et l’isotropie du ciel et la
valeur estimée albédo pour le rayonnement réfléchi du sol.
Pour une surface placée parallèlement à la terre, le taux de la radiation solaire extraterrestre
Gh0, incidente sur cette surface horizontale à un moment donné de l'année est donné par [25]
𝟑𝟔𝟎. 𝐍𝐣
𝐆𝐞𝐱𝐭 = 𝟏𝟑𝟔𝟕. [𝟏 + 𝟎. 𝟎𝟑𝟑. 𝐜𝐨𝐬 ( )] . 𝐜𝐨𝐬(Ф) . 𝐜𝐨𝐬(𝛅) . + 𝐬𝐢𝐧(Ф) 𝐜𝐨𝐬(𝛅) (2.1)
𝟑𝟔𝟓
L’irradiation globale incidente sur une surface horizontale extraterrestre durant la journée peut
être obtenue par l'intégration d'équation (2.1) sur la période du lever au coucher du soleil.
L'équation résultante est :
𝟑𝟔𝟎. 𝐍𝐣
𝐇𝐨 = 𝟖𝟔𝟒𝟎𝟎. [𝟏 + 𝟎. 𝟎𝟑𝟑. 𝐜𝐨𝐬 ( )] . 𝐜𝐨𝐬(Ф) 𝐜𝐨𝐬(𝛅) 𝐜𝐨𝐬(𝛚) + 𝛚𝐬𝐬 (𝛑/𝟏𝟖𝟎). 𝐬𝐢𝐧(Ф) 𝐬𝐢𝐧(𝛅)
𝟑𝟔𝟓
(2.2)
21
CHAPITRE 2 MODELISATION DU RAYONNEMENT SOLAIRE
a) Fraction diffuse
L'indice de diffuse est défini comme le rapport entre la moyenne mensuelle journalière de la
radiation diffuse sur une surface horizontale et la valeur correspondante de l’irradiation
globale journalière Gh(0).est définit par la relation suivante :
𝐃𝐡 (𝟎)
𝐤𝐃 =
𝐆𝐡 (𝟎) (2.3)
L’rayonnement global sur une surface horizontale, Gh est la somme de rayonnement direct
horizontal et de rayonnement diffus, c'est-à-dire :
𝐆𝐡 = 𝐁𝐡 + 𝐃𝐡 (2.4)
b) Indice de clarté
22
CHAPITRE 2 MODELISATION DU RAYONNEMENT SOLAIRE
L'indice de clarté, kT, est défini comme le rapport entre la moyenne mensuelle de la radiation
journalière sur une surface horizontale et la valeur correspondante du hors atmosphère: clarté
de l’atmosphère.
𝐆𝐡
𝐤𝐓 =
𝐆𝐞𝐱𝐭
(2.6)
Où :
Gext: Irradiation globale extraterrestre journalière
Gh(0) : Irradiation globale horizontale journalière
Dh(0) : Irradiation diffuse horizontale journalière
𝝎𝒔𝒓
∫ 𝐆𝐞(𝝎) 𝐝𝝎 (2.7)
𝝎𝒔𝒔
𝝎𝒔𝒓
∫ 𝐆𝐦(𝝎) 𝐝𝝎 (2.8)
𝝎𝒔𝒔
Où :
Ge : irradiation journalière globale estimée.
Gm : irradiation journalière globale mesuré.
Ils ont établi une relation empirique entre le rapport de la moyenne mensuelle de l’irradiation
diffuse journalière Dh(0)à Gextet la moyenne mensuelle de la fraction d’irradiation KT.[26]
23
CHAPITRE 2 MODELISATION DU RAYONNEMENT SOLAIRE
𝐃𝐡(𝟎)
= 𝟏. 𝟑𝟗𝟎 − 𝟒. 𝟎𝟐𝟕𝐊𝐓 + 𝟓. 𝟓𝟑𝟏𝐊𝐓² − 𝟑. 𝟎𝟏𝟖 𝐊𝐓 𝟑 (2.9)
𝐆𝐞𝐱𝐭
Ce résultat a été très fréquemment utilisé mais il a suscité plusieurs critiques : Il a été
obtenu avec une valeur erronée de la constante solaire (1367W/m2).
Les mesures du rayonnement diffus qui ont été utilisées ont été effectuées avec un
pyranomètre à anneau sous-estime ce rayonnement.
Il faudrait tenir compte d’un effet saisonnier. En effet, l’intensité du rayonnement
diffus dépend de la masse atmosphérique qui en moyenne est beaucoup plus grande en
hiver qu’en été.
Kd=1.188-2.272KT+9.473KT²-21.856KT³+14.648KT⁴ (2.11)
En tenant compte de ces divers facteurs Collares-Pereira et Rabl qui dépendra, cette fois-ci,
de l’angle horaire du lever du soleil ωss, paramètre caractérisant la saison considérée. Cette
relation est donnée comme suit:
𝟎. 𝟐𝟔𝟏
𝐊𝐝 = 𝟎. 𝟕𝟕 + (𝟎. 𝟑𝟒𝟕/𝟏𝟖𝟎)(𝛚𝐬𝐬 − 𝟗𝟎) − 𝟎. 𝟓𝟎𝟓 + [ (𝛚𝐬𝐬 − 𝟗𝟎) . 𝐜𝐨𝐬(𝟐(𝐊𝐓 − 𝟎. 𝟗))] (2.12)
𝟏𝟖𝟎
Où
ωss :est l’angle horaire moyen du lever du soleil pendant la période considérée.
24
CHAPITRE 2 MODELISATION DU RAYONNEMENT SOLAIRE
L’algorithme de calcul sur une surface inclinée nécessite la connaissance des ensoleillements
direct et diffus pour chaque heure d’une <<journée moyenne>>. En premier lieu, il convient
donc de calculer la moyenne mensuelle du rayonnement diffus journalier Dh(0) à partir de la
moyenne mensuelle du rayonnement global journalier Gh et en utilisant la corrélation d’Erbs
[30]
Quand l’angle horaire du soleil à son coucher pour le jour moyen du mois est inférieur à 81.4°
Kd=1.311-3.0.22KT+3.1427KT²-1.821KT³ (2.15)
Notons que les deux formules sont valables uniquement si est compris entre [0.3 ; 0.8].
Quand l’angle horaire du soleil à son coucher est supérieur à 81,4° (la moyenne mensuelle de
l’indice de clarté). Ensuite, l’ensoleillement quotidien moyen est réparti en valeurs horaires.
Cela est réalisé grâce aux formules de ≪ Collares-Pereira et Rabl≫ pour le rayonnement
global :
Où
rg : est le rapport de la valeur horaire sur le total quotidien de l’ensoleillement global.
ωss : est l’angle horaire du soleil à son lever exprimé en radians.
ω :est l’angle horaire du soleil pour le milieu de l’heure pour laquelle le calcul est fait, exprimé
aussi en radians ; et rd peut être donné par la de Liu et Jordan :
𝛑 𝐜𝐨𝐬(𝛚) − 𝐜𝐨𝐬(𝛚𝐬𝐬)
𝐫𝐝 = ( (2.18)
𝟐𝟒 𝐬𝐢𝐧(𝛚𝐬𝐬) − 𝛚𝐬𝐬𝐜𝐨𝐬(𝛚𝐬𝐬)
25
CHAPITRE 2 MODELISATION DU RAYONNEMENT SOLAIRE
Où
rd: est le rapport de la valeur horaire et le total quotidien de l’ensoleillement diffus. Pour
chaque heure de la « journée type »
Gh :est l’éclairement global horizontal.-
Dh(0) et Bh(0) : sont les composantes diffuse et directe. Elles sont données par les trois
formules suivantes :
Gh=rg.Gh(0) (2.19)
GT(Β.Α)=BT(Β.Α)+DT(Β.Α)+RT(Β.Α) (2.22)
Où
Gt :est le rayonnement global.
Bt :est le rayonnement direct.
Dt :est le rayonnement diffus.
Rt :est le rayonnement réfléchi par le sol.
Bt(β.α)=B . cos(θi)
(2.23)
26
CHAPITRE 2 MODELISATION DU RAYONNEMENT SOLAIRE
L’éclairement du au rayonnement direct sur un plan horizontal, il est calculé aussi par relation
suivante :
Bh(0)=B .cos(θz) (2.24)
Avec
θz :Angle zénithal du soleil
B :est l’éclairement du au rayonnement solaire direct sur un plan normal et θi l’angle
d’incidence à la surface inclinée. De cela on obtient le rapport suivant :
𝐜𝐨𝐬(𝛉𝐢)
𝐑𝐛 = (2.25)
𝐜𝐨𝐬(𝛉𝐳)
B(β.α)=Rb(Gh(0)-Dh(0) (2.27)
Où
Dh(0), Gh(0), Bh(0) sont respectivement les éclairements dus aux rayonnements global, diffus
direct sur un plan horizontal.
La majorité des mesures des radiations solaires sont pour les surfaces horizontales, en
connaissant les valeurs des radiations directes et diffuses ainsi que les angles solaires. Les
modèles des radiations du ciel, dans ce contexte sont des représentations mathématiques des
radiations diffuses, quand les radiations directe et diffuse sont combinées elles fournissent la
moyenne calculée des radiations solaires sur une surface inclinée à partir des mesures
effectuées sur une surface horizontale. Un des modèles est le modèle des radiations diffuses
isotropes, proposé par Liu et Jordan. [32]
27
CHAPITRE 2 MODELISATION DU RAYONNEMENT SOLAIRE
Dans ce modèle, le rayonnement solaire global incident sur une surface inclinée est supposé
être composer de trois éléments: le des radiation directs, des radiations diffusées du ciel
isotopique et des radiations réfléchie par le sol sou forme de radiations diffusées due a l‘albédo
du sol.(celui qui arrive uniformément de la voûte céleste indépendamment de sa direction).Le
rayonnement diffus sur surface inclinée𝑫(𝜷. 𝜶) peut être obtenu en multipliant le rayonnement
diffus sur surface horizontale par un facteur de visibilité entre la surface et le ciel peut être une
estimer par la relation suivante:
Le rayonnement réfléchi du ciel sur un plan incliné peut être une estimer par l'expression
suivante: 𝐑(𝛃. 𝛂) = 𝛒. 𝐆𝐡(𝟎). [((𝟏 − 𝐜𝐨𝐬(𝛃) ))/𝟐] (2.29)
Le rayonnement global sur une surface inclinée pour un ciel isotrope se résume par :
(𝟏 + 𝐜𝐨𝐬(𝛃))
𝐆𝐭(𝛃. 𝛂) = 𝐑𝐛(𝐆𝐡(𝟎) − 𝐃𝐡(𝟎)) + 𝐃𝐡(𝟎). [ ] + 𝐃𝐡(𝟎). [((𝟏 + 𝐜𝐨𝐬(𝛃)))/𝟐]
𝟐 (2.30)
28
CHAPITRE 2 MODELISATION DU RAYONNEMENT SOLAIRE
(𝟏 + 𝐜𝐨𝐬(𝛃)) 𝛃
𝐆𝐭(𝛃. 𝛂) = 𝐃𝐡(𝟎). ( ) (𝟏 + 𝐅𝐬𝐢𝐧 ( )) . [𝟏
𝟐 𝟐
𝟏 − 𝐜𝐨𝐬(𝛃) (2.31)
+ 𝐜𝐨𝐬²(𝛃)𝐬𝐢𝐧³(𝛉𝐳 + 𝛒𝐆𝐡(𝟎)) ( )]
𝟐
𝐃𝐡(𝟎) 𝟐
𝐅 = [𝟏 − ( ) ] (2.32)
𝐆𝐡
Sa valeur est nulle lorsque le ciel est couvert et l'équation devient pour l’isotrope.
(𝟏 + 𝐜𝐨𝐬(𝛃))
𝐆𝐭(𝛃, 𝛂) = 𝐑𝐛. (𝐁𝐡(𝟎) + 𝐅𝐢. 𝐃𝐡(𝟎)) + 𝐃𝐡(𝟎)(𝟏 − 𝐅𝐢) ( )
𝟐
(𝟏 − 𝐜𝐨𝐬(𝛃)) (𝟏 − 𝐜𝐨𝐬(𝛃))
+ 𝛒𝐆𝐡(𝟎) [ ]( ) (2.33)
𝟐 𝟐
29
CHAPITRE 2 MODELISATION DU RAYONNEMENT SOLAIRE
(𝟏 + 𝐜𝐨𝐬(𝛃)) 𝛃 𝐁𝐡(𝟎)
𝐆𝐭(𝛃, 𝛂) = 𝐃𝐡(𝟎. )(𝟏 − 𝐅𝐢). ( ) + [𝟏 + 𝐅𝐢. 𝐑𝐛 + 𝐬𝐢𝐧³( ). √ ] (2.35)
𝟐 𝟐 𝐆𝐡(𝟎)
L’éclairement global reçu sur une surface inclinée est donnée par la relation suivante:
(𝟏 + 𝐜𝐨𝐬(𝛃))
𝐆𝐭(𝛃, 𝛂) = 𝐑𝐛. 𝐁𝐡 + 𝐃𝐡(𝟎. )(𝟏 − 𝐅𝐢) ( )
𝟐
𝐁𝐡(𝟎) (𝟏 − 𝐜𝐨𝐬(𝛃))
+ [𝟏 + 𝐅𝐢. 𝐑𝐛 + 𝐬𝐢𝐧³(𝛃/𝟐). √ + 𝛒𝐆𝐡(𝟎) ( )]
𝐆𝐡(𝟎) 𝟐
(2.36)
Dans ce modèle, la composante directe dans un plan d'inclinaison (β )et d'orientation (α) et
pour une hauteur (h)est donnée par la relation suivante [27] :
𝐌𝐚𝐱(𝟎, 𝐜𝐨𝐬(𝛉𝐢))
𝐆𝐁(𝛃. 𝛂) = [𝐅𝟏 ] . 𝐅𝐓(𝛉𝐢). 𝐁(𝟎)
𝐜𝐨𝐬(𝛉𝐳) (2.37)
Où:
𝛉𝐢: est l’angle d’incidence
30
CHAPITRE 2 MODELISATION DU RAYONNEMENT SOLAIRE
(𝟏+𝐜𝐨𝐬(𝛃) (2.41)
𝐆𝐝𝐑(𝛃. 𝛂) = 𝐆𝐝(𝟎). [(𝟏 − 𝐅𝟏) ]
𝟐
Les expressions de l’indice de clarté (ε) et la luminosité Δ sont données respectivement par:
𝐁(𝟎) (2.44)
𝐆𝐝(𝟎) +
𝐜𝐨𝐬(𝛉𝐳)
𝛆=
𝐆𝐝(𝟎)
𝐆𝐝(𝟎) (2.45)
𝚫=
𝐆𝐞𝐱𝐭. 𝐜𝐨𝐬(𝛉𝐳)
Perez classifie les divers ciels en huit catégories basées sur l’indice de clarté ε de 1.000 à 1.065
pour ciel couvert (rayonnement direct nul); au-delà de 10.08 pour un ciel sans nuages et entre
les valeurs extrêmes on a des états de ciels nuageux (voir tableau 2.1):
Tableau(2.1):Coefficients du modèle de Perez
31
CHAPITRE 2 MODELISATION DU RAYONNEMENT SOLAIRE
GDt(β.α)=FT.GdC(α,β)+0.856.GdH(α,β)+0.856.GdR(α,β) (2.46)
L’éclairement réfléchi reçu sur une surface inclinée se donne par l’équation suivante:
32
CHAPITRE 2 MODELISATION DU RAYONNEMENT SOLAIRE
ETAPE 8 : calcul des valeurs des pertes de l’irradiation global sur la surface inclinée.
33
CHAPITRE 2 MODELISATION DU RAYONNEMENT SOLAIRE
CONCLUSION
La production d’énergie des systèmes photovoltaïque dépend de l’énergie solaire incidente sur
leurs surfaces captatrices donc il est nécessaire de développer des modèles de prédiction de
l’ensoleillement qui engendre les paramètres du lieu (site).Dans ce chapitre nous distinguons le
modèle de ciel isotropique, les modèles du ciel anisotrope qui tiennent en compte de la
présence de nuages dans le ciel et leur effet sur la diffusion du rayonnement et le modèle
(Perez) qui tiennent en compte de tous les aspects de la diffusion du rayonnement.
34
CHAPITRE 3
OPTIMISATION DE L’IRRADIATION
SOLAIRE - RESULTATS
ET DISCUSSIONS
3.1. Introduction
Dans ce chapitre nous allons présenter les résultats de simulation de calcul de l’irradiation
solaire que peut capter un générateur photovoltaïque à structure porteuse arbitrairement
incliné et incliné. Il est question également de déterminer les deux angles optimaux
(inclinaison, orientation) donnant les maximumes d’énergie reçue sur le plan d’un générateur
photovoltaïque. L’étude de simulation a été réalisée grâce au logiciel Matlab et concerne une
vingtaine de stations en Algérie (20 villes.
Le tableau 3.1 ci-après montre les données de l’irradiation journalière moyenne mensuelle sur
le plan horizontale pour déférentes régions Algérien utilise pour le calcule de l’irradiation
solaire globale sur le plan inclinée. Pour identifier les angles d’inclinaison et d’orientation
optimaux pour lesquels l’énergie reçue attient sa valeur maximale pour différentes sites
algérien, nous procédons à la variation de l’angle d’inclinaison de 0 ° à 90 °et l’angle
d’orientation de -90 à 90.
35
CHAPITRE 3 RESULTAT ET DISCUSSIONS
Ouargla 304 4350 5620 706 7250 769 7910 7420 597 433 3270 2660
0 0 0 0 0
Oran 272 3640 4740 594 6550 708 6920 6260 525 394 2800 2400
0 0 0 0 0
Alger 248 3380 4590 569 6490 720 7130 6440 528 382 2630 2150
0 0 0 0 0
Tamanrasse 442 5480 6060 671 7130 709 7180 6660 602 540 4570 4040
t 0 0 0 0 0
Becher 420 5330 6220 727 7600 754 7420 7080 647 537 4380 3880
0 0 0 0 0
36
CHAPITRE 3 RESULTAT ET DISCUSSIONS
Tableau (3.2) : Irradiation solaire journalière annuelle en fonction azimut et l’inclinaison (M’sila )
Compte tenu de les composantes diffuse et globale du rayonnement solaire sur un plan incliné
relativement importante dans la région M’sila, l’irradiation solaire annuelle (Wh/m2), rapporté
à sa valeur maximale varie avec l’orientation(α) et avec l’inclinaison (β) comme l’indique la
figure :
Figure (3.1) : Variation de l’irradiation journalière annuelle en fonction (α, β) pour le site M’sila
37
CHAPITRE 3 RESULTAT ET DISCUSSIONS
Par la lecture du tableau (3.2) et figure (3.1) on peut déduire que le maximum de valeur de
l’irradiation journalière annuelle pour l’orientation (α) et l’inclinaison (β) varier
respectivement (-20° à 20°) et (30° à 60°) elle est égale à (10400 à 10680) Wh/m². Alors que
les deux angles optimaux pour le site de M’sila sont l’orientation Sud (0°) et inclinaison (40°)
qui donnent l’irradiation maximale journalière annuel. Elle est égale à 10672 Wh/m².
3.2.2. Irradiation solaire pour le site d’Ouargla
Le tableau (3.3) suivant présente la variation de l’irradiation solaire journalière annuelle en
fonction de l’orientation (α) et l’inclinaison(β). pour le site d’Ouargla les données
géographiques utilisées sont latitude (31°) et longitude (5°.24 E)
Tableau (3.3) : Valeurs de l’irradiation solaire annuelle en fonction des angles (α, β) pour le site
d’Ouargla
α/β -60 -50 -40 -30 -20 -10 0 10 20 30 40 50 60
0 8663 8663 8663 8663 8663 8663 8663 8663 8663 8663 8663 8663 8663
10 9177 9315 9435 9531 9603 9646 9661 9646 9603 9531 9435 9315 9177
20 9449 9721 9956 10147 10287 10373 10402 10373 10287 10147 9956 9721 9449
30 9470 9868 10211 10490 10695 10821 10863 10821 10695 10490 10211 9868 9470
40 9240 9751 10193 10551 10815 10977 11031 10977 10815 10551 10193 9751 9240
50 8766 9375 9901 10328 10643 10835 10900 10835 10643 10328 9901 9375 8766
60 8062 8750 9345 9828 10183 10401 10475 10401 10183 9828 9345 8750 8062
70 7149 7897 8542 9065 9451 9687 9767 9687 9451 9065 8542 7897 7149
80 6057 6840 7516 8064 8469 8716 8800 8716 8469 8064 7516 6840 6057
90 4816 5612 6298 6855 7265 7517 7601 7517 7265 6855 6298 5612 4816
La figure (3.2) suivante indique la variation journalière annuelle de l’irradiation pour les
différents orientations et inclinaisons pour la ville d’Ouargla:
Figure (3.2) : Variation de l’irradiation journalier annuelle en fonction (α, β) pour le site Ouargla
38
CHAPITRE 3 RESULTAT ET DISCUSSIONS
D’après la figure (3. 2)et le tableau (3.3), la valeur de l’irradiation solaire annuelle pour le site
Ouargla est maximale pour une orientation (αopt=0°) et une inclinaison (βopt=40°)et La valeur
maximale est de 1031 Wh/m.
Tableau (3.4) :Valeurs de l’irradiation solaire annuelle en fonction des angles (α ,β) pour le site
d’Oran
α/β -60 -50 -40 -30 -20 -10 0 10 20 30 40 50 60
0 7596 7596 7596 7596 7596 7596 7596 7596 7596 7596 7596 7596 7596
10 8067 8187 8290 8374 8436 8473 8486 8473 8436 8374 8290 8187 8067
20 8327 8563 8766 8931 9053 9127 9152 9127 9053 8931 8766 8563 8327
30 8368 8713 9010 9251 9429 9538 9575 9538 9429 9251 9010 8713 8368
40 8189 8631 9014 9324 9552 9692 9740 9692 9552 9324 9014 8631 8189
50 7794 8322 8777 9147 9419 9586 9643 9586 9419 9147 8777 8322 7794
60 7197 7794 8309 8726 9034 9223 9286 9223 9034 8726 8309 7794 7197
70 6415 7063 7621 8074 8408 8613 8682 8613 8408 8074 7621 7063 6415
80 5473 6151 6736 7211 7561 7775 7848 7775 7561 7211 6736 6151 5473
90 4398 5086 5681 6163 6518 6736 6809 6736 6518 6163 5681 5086 4398
A partir le tableau (3.4) et la figure (3.3) ci-dessous, nous obtenons une valeur maximale de
l’irradiation annuelle se trouvant dans la zone (8800 à 9600) Wh/m² pour une orientation (α)
et inclinaison(β) du GPV allant respectivement de -10°SE à 10°OS et de 40° à 50°. Les deux
angles optimaux sont l’orientation α=0° et β =40 donnent une irradiation maximale de 9740
Wh/m².
39
CHAPITRE 3 RESULTAT ET DISCUSSIONS
Figure (3.3): Variation de l’irradiation solaire annuelle en fonction de (α, β) pour le site d’Oran
Tableau (3.5) : Variations de l’irradiation solaire journalier annuelle en fonction (α) et (β) site
d’Alger.
α/β -60 -50 -40 -30 -20 -10 0 10 20 30 40 50 60
0 7592 7592 7592 7592 7592 7592 7592 7592 7592 7592 7592 7592 7592
10 8214 8327 8425 8505 8563 8599 8611 8599 8563 8505 8425 8327 8214
20 8633 8857 9050 9207 9322 9393 9417 9393 9322 9207 9050 8857 8633
30 8838 9165 9448 9676 9845 9949 9983 9949 9845 9676 9448 9165 8838
40 8823 9243 9606 9900 10117 10250 10340 10250 10117 9900 9606 9243 8823
50 8587 9087 9520 9870 10129 10287 10294 10287 10129 9870 9520 9087 8587
60 8137 8703 9192 9588 9881 10060 10120 10060 9881 9588 9192 8703 8137
70 7489 8103 8633 9063 9380 9574 9640 9574 9380 9063 8633 8103 7489
80 6660 7304 7859 8310 8642 8846 8914 8846 8642 8310 7859 7304 6660
90 5678 6331 6895 7353 7690 7896 7966 7896 7690 7353 6895 6331 5678
40
CHAPITRE 3 RESULTAT ET DISCUSSIONS
La variation de l’irradiation solaire annuelle sur un plan incliné et arbitrairement orienté pour
le site d’Alger est illustrée par la figure ci-dessous:
Figure (3.4): Irradiation solaire journalier annuelle en fonction de (α, β) pour site d’Alger
Au vu des résultats donnés au tableau (3.5) et la figure (3.4), on peut déduire que la valeur de
l’irradiation solaire journalière moyenne annuelle pour le site d’Alger est maximale (10340
Wh/m²) pour une orientation (αopt=0°)et une inclinaison (βopt=40°)se trouvant dans la zone
d’orientation qui varier de -40° à 40°et d’inclinaison qui varie de 30° à 60°.
41
CHAPITRE 3 RESULTAT ET DISCUSSIONS
Figure(3. 5): Variation de l’irradiation solaire annuelle en fonction des angles (α, β) pour le site
de Tamanrasset
A partir du tableau (3.6 et la figure (3.5), la valeur de l’irradiation solaire journalière moyenne
annuelle pour le site de Tamanrasset est maximale (10751 Wh/m²) pour une orientation
(αopt=0°) et une inclinaison (βopt=30°).
42
CHAPITRE 3 RESULTAT ET DISCUSSIONS
Tableau (3.7) : Valeurs de l’irradiation solaire journalière en fonction (β et α)pour le site de Bechar
α/β -60 -50 -40 -30 -20 -10 0 10 20 30 40 50 60
0 9564 9564 9564 9564 9564 9564 9564 9564 9564 9564 9564 9564 9564
10 9978 10091 10189 10268 10326 10362 10374 10362 10326 10268 10189 10091 9978
20 10125 10348 10540 10696 10811 10882 10905 10882 10811 10696 10540 10348 10125
30 10001 10327 10608 10836 11004 11107 11142 11107 11004 10836 10608 10327 10001
40 9610 10029 10390 10683 10899 11031 11076 11031 10899 10683 10390 10029 9610
50 8963 9462 9893 10242 10499 10657 10710 10657 10499 10242 9893 9462 8963
60 8081 8645 9131 9526 9817 9995 10055 9995 9817 9526 9131 8645 8081
70 6990 7601 8129 8557 8873 9066 9131 9066 8873 8557 8129 7601 6990
80 5723 6363 6917 7365 7696 7899 7967 7899 7696 7365 6917 6363 5723
90 4318 4969 5531 5986 6322 6528 6597 6528 6322 5986 5531 4969 4318
Figure (3.7) : Variation de l’irradiation solaire journalière en fonction (β et α)pour le site de Bechar
Par la lecture du tableau (3.7) et la figure (3. 6) on peut déduire que l’irradiation solaire
journalier annuelle pour le site de Bechar est maximale pour une orientation (αopt=0°)et une
43
CHAPITRE 3 RESULTAT ET DISCUSSIONS
Tableau (3.8) : Valeurs de l’irradiation solaire journalière en fonction (β et α) pour le site d’El-Oued
α/β -60 -50 -40 -30 -20 -10 0 10 20 30 40 50 60
0 8203 8203 8203 8203 8203 8203 8203 8203 8203 8203 8203 8203 8203
10 8647 8764 8864 8946 9006 9043 9055 9043 9006 8946 8864 8764 8647
20 8862 9092 9290 9451 9569 9642 9666 9642 9569 9451 9290 9092 8862
30 8842 9178 9468 9703 9876 9982 10018 9982 9876 9703 9468 9178 8842
40 8588 9019 9392 9694 9916 10053 10099 10053 9916 9694 9392 9019 8588
50 8107 8621 9065 9425 9690 9852 9907 9852 9690 9425 9065 8621 8107
60 7414 7995 8497 8903 9203 9387 9449 9387 9203 8903 8497 7995 7414
70 6530 7161 7705 8146 8471 8671 8738 8671 8471 8146 7705 7161 6530
80 5483 6143 6713 7176 7516 7725 7795 7725 7516 7176 6713 6143 5483
90 4303 4973 5552 6022 6368 6579 6651 6579 6368 6022 5552 4973 4303
Figure (3.7) : Variation de l’irradiation solaire journalière en fonction (β et α)pour le site d’El-Oued
44
CHAPITRE 3 RESULTAT ET DISCUSSIONS
Au vu résultats donnés dans le tableau (3.8) et la figure (3.7), on peut déduire que l’irradiation
solaire journalière moyenne annuelle est maximale (10099 Wh/m²) reçue par une surface
incliné de βopt=40°et orientée αopt=0°.
Tableau (3.9) : Valeurs de l’irradiation solaire journalière en fonction (β et α) pour le site de Jijel
α/β -60 -50 -40 -30 -20 -10 0 10 20 30 40 50 60
0 7137 7137 7137 7137 7137 7137 7137 7137 7137 7137 7137 7137 7137
10 7604 7717 7816 7895 7954 7990 8002 7990 7954 7895 7816 7717 7604
20 7873 8097 8291 8448 8564 8635 8658 8635 8564 8448 8291 8097 7873
30 7937 8265 8548 8778 8947 9050 9085 9050 8947 8778 8548 8265 7937
40 7794 8216 8579 8874 9092 9225 9250 9225 9092 8874 8579 8216 7794
50 7448 7951 8384 8736 8995 9154 9207 9154 8995 8736 8384 7951 7448
60 6911 7478 7968 8366 8659 8838 8899 8838 8659 8366 7968 7478 6911
70 6197 6813 7344 7776 8093 8288 8354 8288 8093 7776 7344 6813 6197
80 5329 5974 6532 6983 7316 7520 7589 7520 7316 6983 6532 5974 5329
90 4334 4989 5555 6013 6351 6558 6628 6558 6351 6013 5555 4989 4334
Figure (3.8) : Variation de l’irradiation solaire journalière en fonction (β et α) pour le site de Jijel
45
CHAPITRE 3 RESULTAT ET DISCUSSIONS
A partir du tableau (3.9) et la figure (3.8), on peut déduire que l’irradiation journalière
moyenne annuelle atteint un maximum de 9250 Wh/m² pour βopt=40°et αopt=0°.
Au des résultats du tableau (3.10) et la figure (3.9), on peut déduire que la valeur l’irradiation
solaire journalière moyenne atteint un maximum de 9851 Wh/m²pour une orientation αopt=0°
et une inclinaison βopt=40°.
Figure (3.9) : Variation de l’irradiation solaire journalière en fonction (β et α)pour le site de Médéa
46
CHAPITRE 3 RESULTAT ET DISCUSSIONS
A partir du tableau (3.11) et la figure (3.10), on peut déduire que les angles d’orientation et
d’inclinaison correspondant à la valeur de l’irradiation solaire journalière moyenne annuelle
maximale de 10651 Wh/m². Sont respectivement (αopt=0°) et (βopt=30°).
Figure (3.10) : Variation de l’irradiation solaire journalière en fonction (β et α)pour le site de Tindouf
47
CHAPITRE 3 RESULTAT ET DISCUSSIONS
Tableau (3.12) : Valeurs de l’irradiation solaire journalière en fonction (β et α) pour le site de Skikda
α/β -60 -50 -40 -30 -20 -10 0 10 20 30 40 50 60
0 7138 7138 7138 7138 7138 7138 7138 7138 7138 7138 7138 7138 7138
10 7615 7731 7832 7913 7973 8010 8022 8010 7973 7913 7832 7731 7615
20 7895 8124 8322 8482 8600 8672 8697 8672 8600 8482 8322 8124 7895
30 7970 8304 8593 8827 9000 9106 9141 9106 9000 8827 8593 8304 7970
40 7836 8266 8638 8939 9161 9297 9343 9297 9161 8939 8638 8266 7836
50 7499 8012 8454 8813 9077 9240 9294 9240 9077 8813 8454 8012 7499
60 6969 7548 8048 8454 8753 8936 8997 8936 8753 8454 8048 7548 6969
70 6261 6890 7432 7872 8196 8395 8462 8395 8196 7872 7432 6890 6261
80 5398 6056 6624 7086 7425 7634 7704 7634 7425 7086 6624 6056 5398
90 4405 5073 5650 6118 6463 6675 6746 6675 6463 6118 5650 5073 4405
Au vu des résultats obtenu en tableau (3.12) et en figure (3.11),on peut déduire que les angles
(orientation, inclinaison) correspondant à la valeur de l’irradiation solaire annuelle maximale
(9343 Wh/m²)sont respectivement αopt=0° et βopt=40°
Figure (3.11) : Variation de l’irradiation solaire journalière en fonction (β et α) pour le site de Skikda
48
CHAPITRE 3 RESULTAT ET DISCUSSIONS
Tableau (3.13) : Valeurs de l’irradiation solaire journalière en fonction (β et α)pour le site d’Adrar
α/β -60 -50 -40 -30 -20 -10 0 10 20 30 40 50 60
0 9185 9185 9185 9185 9185 9185 9185 9185 9185 9185 9185 9185 9185
10 9573 9680 9772 9847 9902 9936 9948 9936 9902 9847 9772 9680 9573
20 9706 9916 10098 10245 10353 10420 10442 10420 10353 10245 10098 9916 9706
30 9578 9885 10151 10366 10525 10622 10655 10622 10525 10366 10151 9885 9578
40 9193 9589 9930 10207 10411 10536 10578 10536 10411 10207 9930 9589 9193
50 8565 9036 9442 9772 10015 10164 10214 10164 10015 9772 9442 9036 8565
60 7711 8243 8703 9076 9350 9519 9575 9519 9350 9076 8703 8243 7711
70 6658 7235 7734 8138 8436 8619 8680 8619 8436 8138 7734 7235 6658
80 5438 6043 6565 6989 7301 7492 7557 7492 7301 6989 6565 6043 5438
90 4087 4701 5232 5662 5979 6173 6239 6173 5979 5662 5232 4701 4087
Les résultats obtenus au tableau (3.13) et en figure (3.12) montrent que l’irradiation
journalière moyenne annuelle est maximale (10655 Wh/m²) pour les deux angles d’orientation
(αopt=0°) et inclinaison (βopt=30°) se trouvant respectivement dans les intervalles de -20° à
20° et de 20° à 40°.
Figure (3.12) : Variation de l’irradiation solaire journalière en fonction (β et α) pour le site d’Adrar
49
CHAPITRE 3 RESULTAT ET DISCUSSIONS
Tableau (3.14) : Valeurs de l’irradiation solaire journalière en fonction (β et α)pour le site de Sétif
α/β -60 -50 -40 -30 -20 -10 0 10 20 30 40 50 60
0 7269 7269 7269 7269 7269 7269 7269 7269 7269 7269 7269 7269 7269
10 7720 7832 7929 8007 8065 8100 8112 8100 8065 8007 7929 7832 7720
20 7970 8190 8381 8535 8649 8718 8742 8718 8649 8535 8381 8190 7970
30 8011 8333 8611 8837 9003 9105 9140 9105 9003 8837 8611 8333 8011
40 7842 8256 8614 8904 9118 9249 9197 9249 9118 8904 8614 8256 7842
50 7469 7962 8388 8734 8989 9145 9293 9145 8989 8734 8388 7962 7469
60 6902 7460 7941 8332 8620 8796 8856 8796 8620 8332 7941 7460 6902
70 6159 6764 7286 7710 8023 8214 8279 8214 8023 7710 7286 6764 6159
80 5262 5896 6444 6888 7215 7416 7483 7416 7215 6888 6444 5896 5262
90 4239 4883 5439 5890 6223 6426 6495 6426 6223 5890 5439 4883 4239
Au vu des résultats du tableau (3.14) et de la figure (3.13), on peut déduire que le maximum
de l’irradiation journalière moyenne annuelle se trouvant dans l’intervalle (de 8800 à 9600
Wh/m²) correspondant aux angles d’orientation (α) et inclinaison (β) variant respectivement
de -30°SE à 30°OS et de 30° à 50°. La valeur obtenue est de 9293 Wh/m² pour α =0° et
β =50°.
Figure (3.13) : Variation de l’irradiation solaire journalière en fonction (β et α) pour le site de Sétif
50
CHAPITRE 3 RESULTAT ET DISCUSSIONS
Tableau (3.15) : Valeurs de l’irradiation solaire journalière en fonction (β et α)pour le site d’El-Taref
α/β -60 -50 -40 -30 -20 -10 0 10 20 30 40 50 60
0 6509 6509 6509 6509 6509 6509 6509 6509 6509 6509 6508,8 6509 6509
10 6958 7065 7157 7231 7286 7320 7331 7320 7286 7231 7156,6 7065 6958
20 7228 7438 7619 7765 7874 7940 7962 7940 7874 7765 7618,7 7438 7228
30 7310 7616 7881 8096 8254 8351 8383 8351 8254 8096 7881 7616 7310
40 7202 7595 7935 8211 8415 8539 8558 8539 8415 8211 7935,5 7595 7202
50 6906 7375 7781 8109 8352 8500 8550 8500 8352 8109 7780,6 7375 6906
60 6433 6963 7421 7793 8066 8234 8291 8234 8066 7793 7421,1 6963 6433
70 5795 6371 6868 7271 7568 7750 7811 7750 7568 7271 6867,9 6371 5795
80 5014 5617 6138 6560 6872 7062 7126 7062 6872 6560 6137,9 5617 5014
90 4112 4725 5253 5682 5998 6192 6257 6192 5998 5682 5253,3 4725 4112
A partir des résultats donnés au tableau (3.15) et en (3.14), on peut déduire que l’irradiation
solaire journalière moyenne annuelle pour le site de EL-Tarf est maximale une orientation
(αopt=0°) et une inclinaison (βopt=40°) se trouvant dans les intervalles respectivement de -30°
à 30°et de 20° à 60°. La valeur maximale est de 8558 Wh/m².
Figure (3.14) : Variation de l’irradiation solaire journalière en fonction (β et α) pour le site d’El-Taref
51
CHAPITRE 3 RESULTAT ET DISCUSSIONS
Tableau (3.16) : Valeurs de l’irradiation solaire journalière en fonction (β et α)pour le site de Bejaia
Figure (3.15) : Variation de l’irradiation solaire journalière en fonction (β et α) pour le site de Bejaia
52
CHAPITRE 3 RESULTAT ET DISCUSSIONS
D’après le tableau (3.17) et la figure (3.16), l’irradiation solaire journalière moyenne annuelle
est maximale de 10926 Wh/m² pour une orientation αopt=0° et une inclinaison βopt=40°.
Figure (3.16) : Variation de l’irradiation solaire journalière en fonction (β et α) pour le site de Ghardaïa
53
CHAPITRE 3 RESULTAT ET DISCUSSIONS
Figure (3. 17) : Variation de l’irradiation solaire journalière en fonction (β et α) pour le site de Boumerdes
54
CHAPITRE 3 RESULTAT ET DISCUSSIONS
Tableau (3.19) : Valeurs de l’irradiation solaire journalière en fonction (β et α) pour le site d’Illizi
Les angles optimaux donnant l’irradiation solaire journalière moyenne annuelle maximale
10858 Wh/m² sont αopt=0° et βopt=30°.
Figure (3. 19) : Variation de l’irradiation solaire journalière en fonction (β et α) pour le site d’Illizi
55
CHAPITRE 3 RESULTAT ET DISCUSSIONS
Les angles optimaux donnant l’irradiation solaire journalière moyenne annuelle maximale de
9858Wh/m² sont αopt=0° et βopt=40°
Figure (3.19) : Variation de l’irradiation solaire journalière en fonction (β et α)pour le site de Tlemcen
56
CHAPITRE 3 RESULTAT ET DISCUSSIONS
Figure (3.20) : Variation de l’irradiation solaire journalière en fonction (β et α) pour le site de Milla
57
CHAPITRE 3 RESULTAT ET DISCUSSIONS
3 .2.21. Récapitulatif
Dans le tableau (3.22) ci-dessous, nous avons récapitulé les résultats obtenus auparavant
relatifs aux valeurs d’angles d’inclinaisons et d’orientations optimum et la valeur de
l’irradiation journalière moyenne annuelle maximale (Hmax).
Tableau (3.22) : Tableau récapitulant les valeurs de l’irradiation solaire journalière maximale
en fonction (β et α) pour tous les sites étudiés précédemment.
locale latitude longitude βopt αopt Hmax(Wh/m)
58
CHAPITRE 3 RESULTAT ET DISCUSSIONS
Au vu des résultats obtenus, il est à noter que l’orientation en Algérie optimale correspond au
plein sud (αopt=0°), or l’inclinaison optimale annuelle correspond à la latitude +10° (βopt=
latitude +10°).
ΔH=H(αopt;βopt)-H(α β)
Où
H(αopt;βopt) : est l’irradiation maximale
H(α;β) : est l’irradiation à α;β arbitraire
A titre illustratif, nous avons calculé les pertes d’irradiation solaire journalière moyenne
annuelle pour trois sites.
Tableau (2.23) : Valeurs des pertes d’irradiation solaire journalière en fonction de (β et α)pour
le site de M’sila
α/β -60 -50 -40 -30 -20 -10 0 10 20 30 40 50 60
0 2315 2315 2315 2315 2315 2315 2315 2315 2315 2315 2315 2315 2315
10 1805 1674 1562 1470 1403 1361 1347 1361 1403 1470 1562 1674 1805
20 1527 1270 1047 867,2 734 653 626 653 734 867,2 1047 1270 1527
30 1489 1113 789 525,1 331 212 172 212 331 525,1 789 1113 1489
40 1694 1210 793 454,1 204 51,5 0 51,5 204 454,1 793 1210 1694
50 2133 1557 1060 656,2 359 177 115 177 359 656,2 1060 1557 2133
60 2795 2144 1582 1125 789 583 514 583 789 1125 1582 2144 2795
70 3659 2952 2342 1847 1482 1259 1184 1259 1482 1847 2342 2952 3659
80 4698 3958 3318 2800 2418 2183 2105 2183 2418 2800 3318 3958 4698
90 5881 5130 4480 3954 3566 3328 3248 3328 3566 3954 4480 5130 5881
59
CHAPITRE 3 RESULTAT ET DISCUSSIONS
Tableau (2.24) : Valeurs des pertes d’irradiation solaire journalière en fonction de (β et α)pour
le site de Ouargla
α/β -60 -50 -40 -30 -20 -10 0 10 20 30 40 50 60
0 2368 2368 2368 2368 2368 2368 2368 2368 2368 2368 2368 2368 2368
10 1854 1716 1597 1500 1429 1385 1370 1385 1429 1500 1597 1716 1854
20 1582 1310 1075 884,7 744 658 629 658 744 884,7 1075 1310 1582
30 1561 1163 820 541,3 336 210 168 210 336 541,3 820 1163 1561
40 1791 1280 838 480,2 216 54,5 0 54,5 216 480,2 838 1280 1791
50 2266 1656 1130 703,2 389 196 131 196 389 703,2 1130 1656 2266
60 2970 2281 1686 1204 848 630 557 630 848 1204 1686 2281 2970
70 3882 3135 2489 1966 1580 1344 1264 1344 1580 1966 2489 3135 3882
80 4975 4192 3515 2967 2563 2315 2232 2315 2563 2967 3515 4192 4975
90 6215 5420 4733 4176 3766 3515 3430 3515 3766 4176 4733 5420 6215
Tableau (2.25) : Valeurs des pertes d’irradiation solaire journalière en fonction de (β et α)pour
le site d’Adrar
α/β -60 -50 -40 -30 -20 -10 0 10 20 30 40 50 60
0 1470 1470 1470 1470 1470 1470 1470 1470 1470 1470 1470 1470 1470
10 1081 974,4 882 807,4 752,3 718 707,1 718 752 807 882 974,4 1081
20 949 738,7 557 409,7 301,1 235 212,2 235 301 410 557 738,7 949
30 1077 769,6 504 288,8 130 32,8 0 32,8 130 289 504 769,6 1077
40 1461 1066 725 448 243,9 119 76,78 119 244 448 725 1066 1461
50 2090 1619 1212 882,5 639,4 491 440,3 491 639 882 1212 1619 2090
60 2944 2411 1952 1579 1304 1136 1079 1136 1304 1579 1952 2411 2944
70 3997 3419 2921 2516 2218 2036 1974 2036 2218 2516 2921 3419 3997
80 5217 4612 4090 3666 3353 3162 3098 3162 3353 3666 4090 4612 5217
90 6568 5953 5423 4993 4676 4481 4416 4481 4676 4993 5423 5953 6568
CONCLUSION
En utilisant les données de l’irradiation journalière moyenne mensuelle sur le plan
horizontale, nous avons développé un algorithme et un programme à l’aide du logiciel Matlab
calculant l’irradiation journalière moyenne annuelle maximale et les angles optimaux
d’orientation et d’inclinaison des capteurs photovoltaïques. Nous avons validé notre travail
pour plusieurs sites Algériens. Au des résultats obtenus, Il s’est avéré que, l’orientation
optimale annuelle correspond au plein sud quant à l’inclinaison est égale à la latitude+10°
Nous avons également procédé au calcul des pertes d’irradiation solaire pour les capteurs dont
l’orientation et l’inclinaison sont différentes des valeurs optimales.
60
Conclusion générale
L'énergie solaire est l’une des plus importantes parmi les diverses sources d'énergiespossibles.
Une connaissance précise de distribution de la radiation solaire à un endroitgéographique
donné est très indispensable pour l’évaluation de la production des capteurs photovoltaïques.
L'objectif principal de cette étude est de montrer l’importance des angles d'inclinaison et
orientation pour maximiser le rayonnement solaire global incident sur la surface des capteurs
photovoltaïques. Dans ce sens, le présent travail se veut une contribution au développement
d’une méthode d’optimisation de l’irradiation solaire journalière annuelle pour des capteurs
PV à structure porteuse statique.
La production d’énergie des systèmes photovoltaïque dépend de l’énergie solaire incidente sur
leurs surfaces captatrices donc il est nécessaire de développer des modèles de prédiction de
l’irradiation qui engendre les paramètres du lieu (site).Dans ce travail, nous avons développé
un algorithme de calcul en se basant sur le modèle anisotropique de Perez qui tient en compte
de 3 zones du ciel (circumsolaire, horizon, isotropie de la voûte céleste). Nous avons
également introduit dans le calcul l’effet de l’angle d’incidence qui peut engendrer des pertes
supplémentaires.
Afin de valider notre algorithme, nous avons développé un programme en utilisant le logiciel
Matlab. Nous avons effectué le calcul d’optimisation de l’irradiation solaire annuelle et les
angles (orientation et inclinaison) correspondant. Guidés par les résultats obtenus pour
plusieurs sites algériens, nous avons pu conclure que l’orientation et l’inclinaison optimales
correspondent respectivement au plein (α=0°) et β=la latitude du lieu +10°. Nous avons établi
une relation calculant les pertes d’irradiation pour des capteurs dont les deux angles sont
différents des angles optimaux
Comme perspective à notre étude nous préconisons l’étude des effets de masques solaires par
le développement de méthode de calcul de l’irradiation solaire en présence d’ombrage.
61
Références Bibliographie
[1] Le soleil pour tous - initiation à l'énergie solaire pratique, Raymond BRUCKERT, Au
signal,Lausanne, 1980.
[2] BENJAMIN.DESSUS, FRANCOIS PHARABOD – L’énergie Solaire- QUE SAIS-JE ?
presse Universitaire de France, VENDOME, OCT 2002
[3] Zekai ¸Sen. Solar Energy Fundamentals and Modeling Techniques, Springer 2008
[4] .Adolf Goetzberger, Joachim Knobloch Bernhard, “ Crystalline Silicon Solar Cells”,
JohnWiley & Sons, p 100 -135
[6] W. Palz, Solar (Electricity, UNESCO, Butterworths, 1978).
[7] P. Bessemoulin, J. Oliviéri. Le rayonnement solaire et sa composante ultraviolette. La
météorologie 8e série n° 31 septembre (2000), 42-59
[8] Milton Ohring, “The Materials Science of Thin Films”, Academic Press, p 234-242
[9] P, H, Communay. Héliothermique : le gisement solaire, méthodes et calculs, Groupe
de Recherche et d’Edition, France. (2002).
[10] Spectre solaire, courbe téléchargée ; http://rredc.nrel.gov/solar/spectra/am0/
[11] M. Iqbal .Une introduction au rayonnement solaire, Academic Press, New York (1983).
[12] N, Bellel. Contribution à la réalisation et simulation d’un capteur plan à eau type
collecteur; thèse de magistère, Université de Constantine. (1989).
[13] Capderou M. Atlas solaire de l'Algérie. Office des publications universitaires, T1-3,
1988.
[14] A. Duffle and W.A. Beckman, ‘Solar Engineering of Thermal Processes’, 3rd Ed
JohnWiley and Sons, 2006
[15] M. Iqbal .Une introduction au rayonnement solaire, Academic Press, New York (1983).
[16] Viorel Badescu. Modeling solar radiation at the earth’s surface, Springer 2008
[17] Tian Pau Chang . The Sun’s apparent position and the optimal tilt angle of asolar
collector in the northern hemisphere. Solar EnergyVolume 83, Issue 8, August 2009,
Pages 1274-1284.
[18] S, El Mokretar; R, Miri; M, Belhamel. Etude du bilan d’énergie et de masse d’un
séchoir de type serre, applications au séchage des produits agro-alimentaires, Revue des
Energies Renouvelables, Vol 7, p 109-123. (2004).
[19] John A. Duffie and William A. Beckman “Solar engineering of thermal processes
“Second ed 1991.New York. John Wiley.
[20] A. Labouret, M. Villoz. Energie solaire photo voltaique. Edition le Moniteur,
3emmeedition.(2006).
[21] chiba younes,mémoire de magister, «etude de l'utilisation d'un bassin d'eau comme
capteur solaire» universite mentouri constantine 2005.
[22] D.K Edwards, Capteurs solaire ; Edition SCM, Paris, 1979.
[23] A.Boyette ‘’ contrôle – commande d’un générateur asynchrone à double alimentation
avec système de stockage pour la production éolienne ‘’, Thèse de doctorat de l’université
Henry Poincaré, Nancy I,decembre 2006.
[24] M. Meinhardt and G. Cramer, “Past, present and future of grid connected photovoltaic-
and hybridpower-systems,” Proceedings of the IEEE Power Engineering Society Summer
Meeting, Vol. 2, pp. 1283–1288, 2000.
[25] M. Meinhardt and D. Wimmer, “ Multistring-converter: The next step in evolution of
string converter technology,” Proceedings of the EPE’01 Conference, Graz, Austria, 2001.
[26] Perez R, Ineichen P, Seals R, Michalsky J, Stewart R. Modeling daylight availability and
irradiance components from direct and global irradiance. Solar Energy 1990;44(5):pages271-
289.
[27] Gary L. Johnson “wind energy systems”, book, Electronic Edition Manhattan October
10, 2006, pp.61- 70-15733
[28] Efim G. Evseev, Avraham I. Kudish . The assessment of different models to predict the
global solar radiation on a surface tilted to the south. Solar Energy, March 2009, Pages
377-388.
[29] .T.M. Klucher,. Evaluation of models to predict insolation on tilted surfaces.
[30] Hay JE, Davies JA. Calculation of the solar radiation incident on an inclined surface.
First Canadian Solar Radiation Data Workshop 1980.
[31] T.Muneer. Solar Radiation and Daylight Models . Elsevier Butterworth-Hermann
photovoltaïques.
Résumé
Le présent travail de Master travail présente une étude de l’importance des angles
d'inclinaison et d’orientation optimaux pour maximiser l’irradiation solaire globale incidente
sur le plan des capteurs photovoltaïques.
En effet, nous avons développé un algorithme d’optimisation de l’irradiation solaire et le
calcul des pertes inhérents aux capteurs arbitrairement orientés et inclinés.
: مــلخــص
أهمية دراســة الزوايــا مــن اجــل تلقــي الخاليــا الشمسيــة للحـــد األقصى مـــن الشعـــــــــاع العمــل المقــدم في الماستر يقدم
الشمســـي علــــى سطــــح مائـــل ومــــن اجــــل هـــذا نـــدرس إحصــائيــة للشعـــاع الشــامــل و كذلك الضياع في الطاقة الشمسية
.علـــى بعـــض المنــاطــق الجــزائـريـة و ذلــك باختيـــار خوارزمية بسيطة
........اإلشعاع الشمسي ’ سطح مائل و موجه ’ا لواح شمسيــة ’ الضياع: الكلمات المفـتاحية
Abstract
The present work of Master work presents the importance study of optimum angles inclination
and direction allowing sensor the maximum global solar radiation incident in solar panel
photovoltaic.
An effect, we present the algorithm developed for estimating the optimization of radaition and
calcul of cost à in to solar panel oriented and inclined .