Ministère de l’Enseignement Supérieure et de la Recherche Scientifique

Université de Biskra
Faculté des Sciences et de la Technologie
Département de Génie Electrique
Filière de l’Energies Renouvelables

Rapport de TP01 :

Module : Gisements énergétiques renouvelables⎟ 1ère Master ER

De la préparation des étudiants: Sous la Supervision du:

Touaref Fares G01 Prof Saadi Aicha
Abderraouf Slimane G02
Boumedien Abderrahmane G01

Année Universitaire : 2021 – 2022

Introduction :
Dans ce TP, on précède sur quelques notions fondamentales d’astronomie
qui donnent des informations sur le rayonnement reçu par le sol (diffus et direct).
Les données nécessaires à l’établissement du bilan énergétique solaire (les durées
d’insolations, les mesures du rayonnement direct et diffus et les mesures du
rayonnement global) d’un lieu déterminé sont obtenues par des appareils de
mesures. Finalement, on va estimer le rayonnement solaire sur un plan horizontal
et incliné par modélisation selon de modèle approché de PERRIN DE
BRICHAMBEAUT dans le site de Biskra.

Le but de TP :
Etude les deux modèles :
- Modèle simplifiée de l’éclairement.
- Modèle de Perrin de Brichambeaut.

Pour faire connaissance son effet au gisement solaire + faire une comparaison
entre elle pour savoir le modèle la plus proche et la plus précise.

A - Partie théorique

Rappels fondamentaux d’astronomie :

1. Généralité sur la terre et le soleil :
Le soleil est une immense sphère de plasma de diamètre 109 fois
supérieur à celui de la terre, de densité 1,41 ; la pression qui y règne à sa surface
est de plusieurs milliards de fois plus forte que celle observée au niveau de la
surface terrestre, le constituant principal du soleil est l’hydrogène, des réactions
thermonucléaires se produisent alors en permanence provoquant la fusion de
l’hydrogène en Hélium. La température à la surface du soleil est de 5800°k et peut
atteindre 20. 106 à son centre. A chaque seconde, le soleil perd 4,5.106 tonnes de
matières sous forme de rayonnement électromagnétique se propageant à la
vitesse de 300.000 km/s. Notons toute fois que la perte de matière est infime vis-
à-vis de la masse totale du soleil. Ce rayonnement solaire est équivalent au
rayonnement qui serait émis par un corps noir à 5800°k.

Le gisement solaire est une science qui propose l’étude radiométrique et

météorologique de l’atmosphère. En effet, les aspects de l’environnement
relèvent immédiatement de l’objet précité. Les niveaux d’énergies au sol et leur

2
rythme de fluctuation saisonnière, les variations dans l’espace et dans le temps de
l’humidité et de la pression atmosphérique, le régime du vent avec d’autres
paramètres avec moins ressentis, reflètent le concept propre du climat. L’Algérie
offre des possibilités considérables pour le développement des applications solaire
et éolienne via l’exploitation du gisement solaire sous ses différentes formes selon
les diverses modes de captation.

2. Les émissions du Soleil :

L’énergie émise par le Soleil est d’abord sous forme de rayonnements
électromagnétiques dont l’ensemble forme le rayonnement solaire, qui constitue
la seule source externe notable d’énergie pour l’atmosphère.
L’œil humain perçoit une partie seulement du rayonnement solaire, celle
située dans le domaine dit visible, de longueurs d’onde comprises entre 0,40 et
0,70 μm, le rayonnement ultraviolet, est un rayonnement électromagnétique de
longueur d'onde plus courte que le visible (<0.4 μm ). Le rayonnement infrarouge
est un rayonnement électromagnétique de longueur d'onde >0.77 μm.
La plus grande partie de l’énergie solaire est cependant rayonnée dans les
domaines ultraviolet, visible et proche d’infrarouge : 99,2 % de l’énergie solaire
hors atmosphère se trouve entre 200 nm et 4 μm. Au sol, par suite de l’absorption
du rayonnement solaire par la vapeur d’eau, le spectre est limité vers le haut à 2,5
μm environ (figure suivant).

Densité spectrale de l’éclairement énergétique du rayonnement solaire

direct (a) Rayonnement solaire hors de l’atmosphère. (b) Rayonnement solaire direct
au sol pour une hauteur du soleil de 6 degrés et une atmosphère claire standard.

3. La constante solaire :

3
 Le flux solaire moyen sur un plan perpendiculaire à la direction des rayons
et placé en dehors de l’atmosphère terrestre est appelé ; la constante solaire (I).
La valeur moyenne actuellement retenue est 1390 W/m2. La variation de la
distance terre-soleil pendant l’année influence la valeur de la constante solaire. La
relation suivante donne la valeur corrigée de la constante solaire:
𝑰 = 1353(1 + 0.033𝑐𝑜𝑠 𝐽 )
Lorsqu’une plus grande précision est requise, on pourra utiliser le développement
en série de Fourrier :
𝑰 = 1353 + 45.326 𝑐𝑜𝑠 𝐽 + 0.88018 𝑐𝑜𝑠 2𝐽 – 0.0046 𝑐𝑜𝑠 3𝐽 + 1.8037 𝑠𝑖𝑛 𝐽
+ 0.09746 𝑠𝑖𝑛 2𝐽 + 0.18412 𝑠𝑖𝑛 3𝐽

Avec : J = 0.984j, j étant le nombre du jour de l’année à partir de janvier.

4. Déclinaison du soleil :

Le mouvement de la terre sur son axe et autour du soleil est schématisé

sur la figure suivant. Les quatre positions de la terre représentées sur cette figure
correspondent aux solstices et aux équinoxes. On notera en particulier que l’axe
de rotation de la terre est incliné par rapport à la normale du plan de l’écliptique
d’un angle de 23.27’, et la déclinaison du soleil (δ) varie donc de –23.27’ (solstice
d’hiver 21 Décembre) à +23.27’ (solstice d’été 21 Juin) en passant par 0 aux
équinoxes.

L’angle de la déclinaison δ est donné par l’équation suivante:

𝛿 = 23.45 𝑠𝑖𝑛 [0.980( 𝑗 + 284 )]
Lorsqu’une précision est requise, on pourra utiliser le développement en série de
Fourier :
𝛿 = 0.33281 − 22.984 cos 𝐽 − 0.3499 𝑐𝑜𝑠2𝐽 − 0.1398 𝑐𝑜𝑠3𝐽 + 3.7872𝑠𝑖𝑛𝐽
+ 0.03205 𝑠𝑖𝑛2𝐽 + 0.07187 𝑠𝑖𝑛3𝐽

4
La courbe suivante représente la variation de l’angle δ en fonction de nombre de
jours :

5. Angle horaire ω :
L’angle horaire est déterminé par la rotation régulière de terre autour de
son axe polaire, il est donc directement lié au temps solaire vrai TSV qui est un
paramètre essentiel permettant de calcul exact des coordonnées angulaire du
soleil.
360
𝜔= (𝑇𝑆𝑉 − 12)
24
TSV : temps solaire vrai donné en heures.

6. Les différent temps solaire:

 Le temps solaire vrai (TSV) : Directement lié à la rotation de terre sur elle
même ; il est défini par l’angle horaire ω c’est-à-dire angle formé par le plan
méridien passant par le soleil et le plan méridien du lieu.
𝜔
𝑇𝑆𝑉 = 12 +
15
 Le temps solaire moyen TSM : C’est le temps qui correspondrait à une rotation
uniforme de la terre autour du soleil. Il diffère peu du TSV (écart maximal ET=
16 minutes) donc :
𝑇𝑆𝑀 = 𝑇𝑆𝑉 + 𝐸𝑇
Où l’équation de temps (ET) est donnée par la relation suivante :
𝐸𝑇 = 0.0072𝑐𝑜𝑠𝐽 − 0.0528𝑐𝑜𝑠 2𝐽 − 0.0012𝑐𝑜𝑠3𝐽 − 0.1229𝑠𝑖𝑛 𝐽
− 0.1565𝑠𝑖𝑛 2𝐽 − 0.0041𝑠𝑖𝑛 3𝐽

5
 Equation de temps.

Pour repérer la position du soleil dans le ciel, il est d'usage d'utiliser un

système-de coordonnées locales dit de coordonnées azimutales, défini en un point
de la surface terrestre, la direction du soleil est repérée grâce à deux angles :
 Hauteur du soleil (h) : C’est l’angle que fait la direction du soleil avec sa
projection sur le sol, il varie de 0° à 90° dans l’hémisphère sud (Nadir) et
s’annule au lever et au coucher du soleil et est maximal au midi solaire.
 Azimut du soleil (a) : Ou que fait la projection de la direction du soleil avec la
direction du sud, cet angle étant mesuré positivement vers l’ouest.

Rayonnement solaire :
L’énergie émise par le soleil voyage jusqu’à la terre sous forme de
rayonnement électromagnétique. Ce rayonnement est semblable aux ondes
radioélectriques mais il possède une gamme de fréquence différente. La quantité
d’énergie disponible au cours du trajet entre le soleil et l’atmosphère extérieur de
la terre est égale à environ 1367 W/m2. Il en résulte qu’au cours d’une journée
ensoleillée, la quantité d’énergie solaire disponible à la surface de la terre dans la
direction du soleil est généralement plus ou moins de 1000 W/m2. Quel que soit le
moment, l’énergie solaire disponible est principalement fonction de la position du
soleil dans le ciel et de la nébulosité. Sur une base mensuelle ou annuelle, la
quantité d’énergie solaire disponible dépend également de l’emplacement
géographique. De plus, l’énergie solaire utilisable dépend de l’énergie solaire
disponible.
Le rayonnement solaire global qui atteint la surface du capteur comprend
un rayonnement direct, provenant en lignes droites venant directement du soleil,
et un rayonnement diffus, qui provient de manière non isotrope, de toutes les
directions de l’espace. Une partie du rayonnement diffusé est renvoyée vers
l'espace ; le reste est transmis jusqu’au sol, voir figure suivant.

6
 On définit l'albédo comme le rapport entre l'énergie réfléchie par
diffusion et l'énergie incidente tableau suivant :

Surface Réflectivité moyenne

Sol enneigé 0.7

Sol recouvert de feuilles mortes 0.30
Herbe verte 0.26
Forêt en automne ou champs dorés 0.26
Galets et pierres blanches 0.20
Herbe sèche 0.20
Sol argileux 0.17
Foret en hiver (arbres conifères sans neige) 0.07
Plan d’eau (soleil haut h > 030°) 0.07

7
 De part sa situation géographique, l’Algérie possède une source
énergétique importante. Elle reçoit annuellement une énergie de l’ordre de 2500
KWh /m2/ an. La durée d’insolation sur le quasi totalité du territoire national
dépasse les 2000 heures annuellement et peut atteindre les 3900 heures (hauts
plateaux et Sahara). L’énergie acquise quotidiennement sur une surface
horizontale de 1m2 est de l'ordre de 5 KWh sur la majeure partie du territoire
national, soit prés de 1700 KWh/m2/an au Nord et 2263 kWh/m2/an au Sud du
pays.
Régions Région côtière Hauts Plateaux Sahara

Superficie (%) 4 10 86

Durée moyenne
d’ensoleillement
2650 3000 3500
(Heures/an)

Energie moyenne
1700 1900 2650
reçue (KWh/m2/an)

Estimation du rayonnement solaire :

En fonction des besoins d’énergie, pour estimer le rayonnement solaire, il
existe plusieurs méthodes et des instruments classiques de mesure de
rayonnement solaire qui le permettent, les plus connus sont :
- Le Pyranomètre : C’est un appareil qui est utilisé pour la mesure du rayonnement
global reçu par une surface plane.
- L’héliographe : Cet instrument mesure la durée d’ensoleillement journalière.
- Pyrhéliomètre : C’est un radiomètre qui mesure le rayonnement direct reçu par
une surface normale aux rayons du soleil.

L’objet de cette TP, n’est pas de présenter ces instruments. Donc il existe
d’autres modèles basés sur le traitement d’images, permettant d’estimer le flux
de rayonnement solaire, les plus connus sont les approches physiques, analytiques
et statistiques. Pour cela on a choisi un modèle approché de PERRIN DE
BRICHAMBEAUT qui donne de bons résultats pour l’estimation d’irradiation sur un
plan horizontal et incliné respectivement (exemple d’application dans le site de
Biskra).

Modèle de PERRIN DE BRICHAMBEAUT (plan horizontal) :

Le rayonnement global sur un plan horizontal est la somme du
rayonnement direct et diffus, donné par:
𝐺ℎ = 𝐼ℎ . 𝑠𝑖𝑛 ℎ + 𝐷ℎ

8
Où :
Ih: le rayonnement direct reçu par une surface horizontale.
Dh: le rayonnement diffus provenant du ciel, reçu par une surface horizontale.
Gh: le rayonnement global reçu par une surface horizontale.

Le premier composant de rayonnement ‘direct’ est donné par l’expression

suivante:
−1
𝐼ℎ = 𝐴. 𝑒𝑥𝑝 ×
𝐵. sin⁡
(ℎ + 𝐶)
Avec A, B et C sont des constantes empiriques qui dépendent de la nature du ciel,
dont les valeurs sont:
A =1210 B=6 C=1 pour un ciel très clair (CTC)
A = 1230 B =3.8 C = 1.6 pour des conditions normales de ciel clair (CNCC)
A = 1260 B = 2.3 C=3 pour un ciel clair pollué. (CCP)
La deuxième composante du rayonnement ‘diffus’ est donnée par :
𝐷ℎ = 𝐸. (sin ℎ)0.4
La valeur de la constante empirique E est :
E = (3/4).125 pour un ciel très clair (CTC)
E = 125 pour des conditions normales de ciel clair (CNCC)
E = (4/3).125 pour un ciel clair pollué (CCP).
La formule globale de PERRIN DE BRICHAMBAUT est obtenue directement sous la
forme :
𝐺ℎ = 𝐴. (𝑠𝑖𝑛 ℎ)𝐷
Avec:
A = 1150 D = 1.15 pour un ciel très clair (CTC).
A = 1080 D = 1.22 pour des conditions normales de ciel clair (CNCC).
A = 990 D = 1.25 pour un ciel clair pollué (CCP).

B - Partie pratique

1. Traçage les deux modèles en utilisant Matlab/Simulink :

Le premier modèle : modèle simplifié :

9
Le deuxième modèle : modèle de Perrin Brichambeaut :

10
clc
clear all
em=1000;
i=1;
tss(i)=0;
ess(i)=0;
gt1(i)=0;
l=34.8;
l1=(l*pi)/180;
k=1;
if k==1;
a=1150;d=1.15;
end
if k==2;
a=1080;d=1.22;
end
if k==3;
a=990;d=1.25;
j=100
p=0.984*j;
end
for ts=4:1:20
hs=(pi/12)*(ts-6);
es=em*sin(hs)
tss(i)=ts;
ess(i)=es;
%module perrin
e2=0.33281-22.984*cos(j)-0.3499*cos(2*j)-
0.1398*cos(3*j)+3.7872*sin(j)+0.03205*sin(2*j)+0.07187*sin(3*j);
sig=(e2*pi)/180;
i1=(ts-12)*15;
w=(i1*pi)/180;
hr=asin(cos(l1)*cos(w)+sin(l1)*sin(sig));
g=a*(sin(hr))^d;
gt=real(g)
gt1(i)=gt;
i=i+1;
end
figure(1);
plot(tss,ess,'r')
figure(2);
plot(tss,gt1,'b')

11
Traçage les deux modèles en utilisant Matlab/Simulink :

12
13
 2. Exécution les modèles pour les modèles temps (4 saisons …) :

24/12/22 20/03/22

21/06/22 11/09/22

Interprétation les résultats obtenus :

A partir des résultats obtenus, nous constatons que pour le site de Biskra, le
modèle d’estimation de Perrin de Brichambaut donne bonne concordance avec les
valeurs mesurées. Par contre, le modèle simplifié pour un type de ciel les pluparts
des valeurs calculées par ce modèle sont supérieures aux valeurs mesurées sur ce
site. Les écarts observés entre les courbes expérimentales et les courbes
théoriques restent très réduits, sauf l’écart de modèle de Perrin de Brichambaut à
un type de ciel reste notablement élevé. Mais dans le modèle simplifié il n’y a pas
changement dans le type de ciel, c’est pourquoi il n’est pas pris en compte.

14
3. L’organigramme suivez pour les deux modèles :
Organigramme : Estimation du rayonnement solaire :

15
Pour le premier modèle : modèle simplifié :

Début

- Déclaration des paramètres :

em = 1000
- Initialisation :
ts = 4 ; tss = 0 ; ess =0
(ℎ𝑠)
𝑡𝑠 − 6
𝐸𝑞 (2): 𝑒𝑠 = 𝑒𝑚. 𝑠𝑖𝑛⁡

Calcule de hs
Eq (1)
12
𝜋
Eq 1 : ℎ𝑠 =

Calcule de es
Eq (2)

ts = ts + 1

Oui
if
ts ≤ 20

Non

Tracer es

Fin

16
Pour le deuxième modèle : modèle de Perrin Brichambeaut :

Début
eq(1) : hs.
eq(2) : es.
eq(3) : l1.
Initialisation ts=4, eq(4) : p.
tss=0, gt1=0, i=1, K=1.

Déclaration des paramètres :

Introduction des paramètres :
eq(1), eq(2), eq(3), eq(4).

eq(1),

eq(1), a,d
eq(1),

K=2
Type du ciel
a,d K=K+1

Caractéristique de la journée
Calculer hs, w, es, i1, sig, g.

Calculer e2/gr i=i+1

ts= ti+1

Oui
if ts ≤ 20

Non

Tracer gt1

Fin

17
4. Comparaison entre les deux modèles :
Il est vrai que des valeurs spéciale de modèle simplifié de l'éclairement bien
ensemble par rapport au la modèle de PERRIN de Brichambeaut mais n'oubliez pas
que le deuxième modal tient compte des conditions environnantes ( type du ciel,
du climat et de l'atmosphère), et il y a de nombreux éléments qu'il prend en
compte, comme le type de rayonnement, ce qui signifie qu'il ne calcule que la
valeur du rayonnement exploité. C'est pourquoi je préfère le deuxième modèle
car il est plus proche de la réalité, contrairement au premier modèle, que je
considère être pris comme étant dans des conditions proches de l'idéal.

Conclusion :
L’application des formules du rayonnement solaire pour la région de Biskra
donne de bons résultats en comparaison avec les valeurs expérimentales.
L’utilisation des formules simplifiées permet l’estimation rapide de la puissance
incidente et de l’énergie journalière et mensuelle .La détermination de l’angle
d’inclinaison optimal est très importante pour pouvoir obtenir le maximum
d’énergie. On constate d’après le calcul de l’énergie globale qu’on peut utiliser un
angle optimal mensuel ou bien plus facilement un angle saisonnier et cela en
changeant l’inclinaisons des plans receveurs du rayonnement solaire quatre fois
par an .Une autre méthode qui permet d'obtenir le maximum d'énergie c'est la
méthode des système de poursuite du soleil qui est très efficace surtout dans des
conditions idéales dépendant du climat.

18