3ème conférence Internationale des énergies renouvelables CIER-2015

3ème
International conférence
Journal Internationale
of Scientific Research desTechnology
& Engineering énergies(IJSET)
renouvelables CIER-2015
International Journal of Scientific Research & Engineering Technology (IJSET)
Vol.4 pp.135-139

Modélisation et Simulation de Poursuite Aveugle

d’un Héliostat pour Centrale à Tour
A. Takilalte, O. Mahfoud, M. Debbache, H. Karoua, S. Bouaichaoui
Centre de Développement des Energies Renouvelables, CDER
BP 62 Route de l'Observatoire, Bouzaréah, 16340, Algiers, Algeria
[email protected]
[email protected]

Abstract— Le présent travail est une synthèse et simulation des Avant de dériver les équations de poursuite du soleil, les
équations de poursuite Azimut-Elévation (AE) d’un héliostat angles indiquant la position du soleil doivent être présentés en
appartient à un champ d’une centrale à tour. La simulation premier qui peuvent être décrites mathématiquement par
d’évolution diurne des angles AE nous a permet d’analyser la rapport à n'importe quel endroit sur la surface de la terre, par
plage de mouvement en fonction des saisons de l’année et en
l’angle d'altitude (α) du soleil et l’angle d’azimut (A). La
fonction des positions d’Héliostat dans le champ définies par
l’angle focale et l’angle de face. Les simulations établies seront figure 1 représente la composition du système de coordonnées
menées par des discutions et des interprétations en tirant des attaché à la fois au centre de la terre et à la surface de la terre.
conclusions.

Mots clés — Héliostat, poursuite AE, centrale à tour,

modélisation, plage de mouvement.

I. INTRODUCTION
Les centrales solaires à tour (CSP) offrent une alternative
sérieuse pour la production d’électricité d’origine
renouvelable [1]. Le champ optique des centrales solaire,
composé exclusivement des Héliostat, représente près de 40 %
du cout d’investissement de la centrale et est le maillon faible
de la chaine de conversion énergétique [2]. Lorsque
l’Héliostat est équipé d’un système suiveur, la production
électrique atteint alors son maximum [3],[4],[5].
Bien que les modèles des Héliostat contrôlés par capteurs
(poursuite sensible) existent, la plupart des Héliostat en
opération maintenant, sont contrôlés par des ordinateurs ou
microcontrôleurs à partir de la théorie astronomique Fig. 1 Repérage du soleil lié au système de coordonnées centre de la terre et
(poursuite aveugle) [4], [5]. surface de la terre [1].
Dans ce travail, la poursuite aveugle Alta-Azimutale de
Dans le système de coordonnées centre de la terre, (S’) est
l’Héliostat a été adoptée (Fig. 1). Ceci est doté de deux
le vecteur pointé vers le soleil décrit en termes de l’angle
moteurs d’Azimuth-Elevation (AE) respectivement, lui
permettant de corriger à tout instant (en temps réel), les rayons horaire   et l’angle de déclinaison   .
incidents du soleil au cadre de récepteur situé au sommet de la L'angle de déclinaison   , peut être calculé en utilisant la
tour [6], [7], [8] (Fig. 3). formule suivante [9] :
La modélisation et la simulation du mouvement d’un
Héliostat est indispensable afin d’avoir une idée sur la vitesse (1)
et la plage parcourue dans les sens Azimut et Elévation quel
que soit les paramètres d’entrées notamment les paramètres
géographiques et la position d’Héliostat dans le champ. Aussi, Où N est le nombre de jours dans une année, N=1 pour 1
cette étude peut être considérée comme une base, nous sert à Janvier et N=42 pour le 11 Février,
valider le système de la poursuite qui sera réalisé L'angle horaire   , peut être calculé comme suit [10] :
ultérieurement.

II. REPERAGE DE LA POSITION DU SOLEIL   15TSV  12 (2)

Copyright IPCO-2015
ISSN 2356-5608
3ème conférence Internationale des énergies renouvelables CIER-2015
International Journal of Scientific Research & Engineering Technology (IJSET)
Vol.4 pp.135-139

TSV est le temps solaire en 24 heures plutôt que AM / PM .

Il peut s’écrit de la forme suivante [11], [12]:

(3)

TL : est le temps donné par la montre.

DE : décalage horaire par rapport au méridien de Greenwich
(Égale 1 pour l'Algérie).
L : longitude de lieu.
Et : la correction de l’équation des temps.

(4)
Avec:
) (5)
Fig.2 Rayons incident et réfléchi [6].
Etant donné que l'observateur Q  (Fig. 1) est situé sur la
surface de la terre à la latitude   , la position du vecteur
solaire S  est définie dans le système de coordonnées
surface de la terre en fonction de l'angle d'altitude du soleil
  et un angle d'azimut du soleil  A .
Les angles d'altitude et d'azimut du soleil sont définis
respectivement comme suite [1] :*

  sin 1 sin  sin   cos cos cos  (6)

 sin  cos   cos cos sin   (7)

A  cos1  
 cos 

Si sin   0 donc : A  2  A Fig. 3 Définition des deux angles de poursuite AE

III. LA METHODE DE POURSUITE TYPE AZIMUTH-ELEVATION

(AE) Pour la méthode de poursuite du soleil type Azimuth-
En fait, les expressions mathématiques des formules de Elevation  AE  , le système de coordonnées est lié au repère
poursuite pour la méthode AE ont été déjà obtenues par (surface-terre) où le vecteur normal de l'héliostat H AE  est
« Stine » et « Harrigan » en 1985 [13] et «Chen et al» en 2001
[14]. défini en fonction de l'angle d'azimut et de l'angle
Les positions des Héliostats par rapport à la tour peuvent d'élévation (voir fig. 3). Ces derniers peuvent être
être définies en utilisant l'angle de face ( et l'angle focal   . dérivés exclusivement à partir de la loi de réflexion
Ce dernier est une relation directe avec la distance (d) concernons le vecteur de position du soleil (S), le vecteur de
Héliostat-Tour [1], [15] (voir la fig. 2). position de la cible (T) et le vecteur normal de
l’Héliostat H AE  . Ce qui donne [16] :

(8)
(9)

Où H t : Hauteur de la tour
r : Hauteur d’Héliostat (10)

La position de repos de l’Héliostat est telle qu’il est dirigé

vers le zénith ( = +π/2) et orienté vers le sud ( = 0). Le
sens positif d’azimut est dirigé vers l’Ouest.

Copyright IPCO-2015
ISSN 2356-5608
3ème conférence Internationale des énergies renouvelables CIER-2015
International Journal of Scientific Research & Engineering Technology (IJSET)
Vol.4 pp.135-139

IV. RESULTATS ET ANALYSE de l'année. Ces jours ont été choisis librement à 21 Décembre
(solstice d'hiver), 21 Mars (équinoxe de printemps) et 21 Juin
A. Validation des Equations Héliostatiques (solstice d'été). De telles variations peuvent être accomplies
Notre modèle de simulation a été validé par rapport une pour un jour de l'année en utilisant la procédure numérique
étude intéressante menée pour Targastonne par MONTENON développée dans cette recherche. Les données sur la base
[16]. desquelles la Fig. 5 est obtenue, sont liées aux positions sur le
côté sud de la tour avec un angle de face =0°) ainsi qu’un
45
Simulation de la position d,Héliostat à Targasonne
angle focal   15 qui est équivalent d’un (d=35m).
La Fig. 5 (haut), représente l’évolution de l’angle d’azimut
40 de l’héliostat en fonction du temps, nous pouvons remarquer
en premier lieu que durant les trois saisons, l’Héliostat est
Angle d,élévation (Degré)

35 dirigé vers l’Est au matin. Ainsi, sa direction varie le matin de

l’est vers l’ouest sur la plage de [-20 0], [-38 0], [-60 0] durant
30
les trois saisons successivement [hiver], [printemps], [été]. Par
contre, l’après-midi sa course varie sur la plage de [0 +20], [0
25
+22], [0 +40] durant les trois saisons successivement [hiver,
équinoxe, été].
20
Notre simulation
Les trajectoires d’élévation d’Héliostat sont symétriques par
Similation de MONTENON rapport à midi comme cela apparaît sur la Fig. 5(bas). Le point
15
-60 -40 -20 0 20 40 60 H s’élève en début de journée, atteint son point le plus haut
Angle d,azimut (Degré)
vers 13 heures et décroît jusqu’à la fin de journée.
Fig. 4 Validation de notre modèle de poursuite de la méthode (AE) avec Les courses sont plus hautes vers le solstice d’été et plus
ceux utilisé par MONTENON [16] à 21 Juin pour Targastonne
( La légère
basses vers celui d’hiver. Donc, en été, elle varie sur la plage
différence apparue est due au faite qu’on n’a
, L=1.44°),
[17 50], en printemps sur la plage [0 +38] et en hiver sur la
plage [0 +25].
La légère différence apparue est due au faite qu’on n’a
pas utilisé les mêmes modèles inclus dans la détermination de
Simulation de l angle d azimut d Heliostat
course de soleil. Ainsi, la négligence du terme de la correction 60

de la hauteur de soleil dû à la réfraction atmosphérique, car la N=81 (equinoxe)

zone d’implantation dans notre cas est dans un climat 40

N=173 (solstice d été)
N=356 (solstice d hiver)
favorable et stable (Ghardaïa), et ce paramètre ne perturbe pas
Angle d azimut d Heliostat (Degré)

l’estimation de la hauteur [16]. 20

B. Simulation des Equations Héliostatiques 0

Une simulation des angles de poursuite est faite du 6 H au
18H pour les paramètres d’entrés cités ci-dessous : -20
Dans cette étude, nous avons choisi la zone de Ghardaïa qui
est un site, par sa situation géographique au sud de l’Algérie, -40
aride, de gisement solaire promoteurs (spécialement
rayonnement direct) et il est prévu de recevoir prochainement, -60
6 8 10 12 14 16 18
une implantation de plusieurs projets de CSP [17], [18]. Temps Légal (Heure)

Simulation de l angle d élévation d Heliostat

50
TABLEAU I
DONNEES GEOGRAPHIQUE DU SITE DE GHARDAIA, PARAMETRE
HELIOSTAT-TOUR 40
Angle d élévation d Heliostat (Degré)

Localisation du Site Hauteur de la Hauteur de

30
Tour (Ht) l’Héliostat (r)
Latitude ( Longitude (L) 11m 1m
20
32.4° 3.8°

10

1) Etude du Mouvement Diurne de l’Héliostat Durant 0

N=81 (equinoxe)
N=173 (solstice d été)
Certaines Saisons de l’Année (Impact de Nombre du Jour N: N=356 (solstice d hiver)

Les variations horaires des angles caractéristiques d'un -10

6 8 10 12 14 16 18
Héliostat typique ont été étudiées numériquement. Les Temps Légal (Heure)

résultats sont montrés sur la Fig. 5, pour trois différents jours

Fig. 5 Evolution des angles d’Azimut (en haut) et d’Elévation (en
bas) de l’héliostat à Ghardaïa en fonction du nombre du jour pour
=0°,

Copyright IPCO-2015
ISSN 2356-5608
3ème conférence Internationale des énergies renouvelables CIER-2015
International Journal of Scientific Research & Engineering Technology (IJSET)
Vol.4 pp.135-139

2) Etude du Mouvement Diurne de l’Héliostat en  L’impact de l’angle de face : les graphes de la Fig. 7,
Fonction de sa Position sur le Champ : illustrent l’évolution des deux angles de poursuite de
l’Héliostat pour N=81 (équinoxe), de distance d=35m,
 L’impact de l’Angle Focal : Dans ce cas, nous avons =15°). Dans tous les cas étudiés les Héliostats sont situés au
simulé la position d’héliostat pour différente distance par Nord de la tour.
rapport à la tour, ce qui traduit par la variation de l’angle Pour =45° (Fig. 7 en haut), la variation de l’angle
focal   , nous avons considéré un Héliostat positionné au d’azimut au matin est très lente car les Héliostats sont au côté
nord de la tour avec un angle de face =0°) et nous avons Est, cette vitesse va s’accélérer plus en plus au-delà de midi
pris comme saison le printemps (Equinoxe du printemps), dû à l’éloignement du soleil à l’autre côté d’Héliostat. Le
c’est-à-dire (N = 81). Les courbes de la Fig. 6, représentent phénomène est réciproquement induit à ceux situés au côté
clairement l’effet de la variation de l’angle focal sur Ouest =-45°).
l’évolution des angles de poursuite de l’héliostat considéré. Contrairement aux Héliostats en face de la tour ( =0°), une
Nous pouvons remarquer en premier lieu que les courses en variation linéaire de la vitesse près de 5°/h revient au faite que
Azimut et en Elévation sont d’autant plus importantes que la ceci garde la symétrie entre le coté d’aube et coucher de soleil
hauteur de cible est plus grande, c’est- -dire, d’autant plus que par rapport au point focalisant.
l’héliostat est proche de la tour (d=10m) et on constate aussi Pour les angles d’Elévation (Fig.7 en bas), cette symétrie de
que cette variation n’est pas linéaire au fil de la distance =0° est évidement apparue de forme Gaussien atteinte sa
focale. crête environ midi. L’angle d’élévation démarre par un angle
faible orientant les miroirs vers le bas et continue de s’élever
par une vitesse de 5°/h. En poursuivant la hauteur de soleil, les
Simulation de l angle d azimut d Heliostat miroirs vont se baisser vers le soir.
Si = 45°, les mêmes formes sont obtenues sauf qu’elles
40
D=10m(Lambda=45°)
30 D=35m(Lambda=15.94°) semblent décalées, où les miroirs ont une hauteur plus large
D=70m(Lambda=8.13°)
20 gamme et à =+45° ils prennent beaucoup de temps de
Angle d azimut d Heliostat (Degré)

monter à l’inverse de =-45°.

10

0
Simulation de l,angle d,azimut d,Heliostat
-10 30

-20 20
Angle d,azimut d,Heliostat (Degré)

-30 10

-40
0

-50
6 8 10 12 14 16 18 -10
Temps Légal (Heure)

-20
Simulation de l angle d élévation d Heliostat
60 -30

angle de face=-45°
50 -40
angle de face= 0°
Angle d élévation d Heliostat (Degré)

angle de face=+45°
40 -50
6 8 10 12 14 16 18
Temps Légal (Heure)
30
Simulation de l,angle d,élévation d,Heliostat
45
20

40
10
Angle d,élévation d,Heliostat (Degré)

35

0 D=10m(Lambda=45°)
30
D=35m(Lambda=15.94°)
D=70m(Lambda=8.13°)
-10 25
6 8 10 12 14 16 18
Temps Légal (Heure) 20

15
Fig. 6 Evolution des angles d’Azimut (en haut) et d’Elévation
(en bas) de l’Héliostat à Ghardaïa en fonction l’angle focale 10
pour =0°, N=81 (équinoxe) angle de face=-45°
5 angle de face= 0°
angle de face=+45°
0
6 8 V. 10CONCLUSION
12 14 16 18
Temps Légal (Heure)

Fig.7 Evolution des angles d’Azimut (en haut) et d’Elévation (en bas)
de l’Héliostat à Ghardaïa en fonction l’angle de face pour =15°, N=81
(équinoxe)

Copyright IPCO-2015
ISSN 2356-5608
3ème conférence Internationale des énergies renouvelables CIER-2015
International Journal of Scientific Research & Engineering Technology (IJSET)
Vol.4 pp.135-139

V. CONCLUSION [18] [18] A. Takilalte, La production de l’électricité à base des CSP en

Algérie : réalité et perspectives, Bulletin des Energ. Ren., Algérie,
Dans ce travail, la simulation d’évolution diurne des angles N°30, 2014.
AE est réalisée en fonction des positions d’Héliostat dans le
champ, et en fonction des saisons de l’année.
En générale, la plage de mouvement obtenue de l’ensemble
des Héliostats dans le champ simulé peut varier en Azimut de
-100° à 50° (écart de 150°) et en Elévation de 0° à 50° (écart
de 50°) pour l’élévation avec une vitesse max près de 20°/h
(équivalent de 0.1 mrd/s) au cours de l’année.
Le résultat est une méthodologie simple et efficace qui peut
offrir des valeurs de la plage du mouvement et sa vitesse pour
les deux angles de poursuite et pour n’importe quelle saison,
quelque ce soit la position d’Héliostat par rapport à la tour
dans un site donné.

REFERENCES

[1] K.K. Chong et M.H. Tan, “Range of motion study for two different
sun-tracking methods in the application of heliostat field”, Sol. Energy,
2011.
[2] (2010) Ademe site web. [Online]. Available: http:// www.ademe.fr.
[3] T. Bendib et al, “Commande automatique d’un système de poursuite
solaire à deux axes à base d’un microcontrôleur PIC16F84A”, Revue
des Energies Renouvelables, Volume 11 N°4, Page 523 – 532, 2008.
[4] M. Haddadi, Élaboration d’une Commande pour Héliostat, Revue des
Energies Renouvelables, Volume 11 N°4, Page 153 – 157, 2003.
[5] A. Gama et al, Etude et réalisation d’un dispositif de poursuite solaire
sensible destiné à un prototype d’un concentrateur cylindro-
parabolique, Revue des Sciences Fondamentales et Appliquées,
Volume 2 N°.1, Pages 137-148, 2010
[6] Y. T. Chen et al, Comparison of Two Sun Tracking Methods in the
Application of a Heliostat Field, Journal of Solar Energy Engineering,
Voume 126, Février 2004.
[7] Minghuan Guo, Properties of a general azimuth–elevation tracking
angle formula for a heliostat with a mirror-pivot offset and other
angular errors, ScienceDirect - Solar Energy, Pages 159–167, 2013.
[8] A. Gama et al, Etude et réalisation d’un concentrateur cylindro
parabolique avec poursuite solaire aveugle, Revue des Energies
Renouvelables, Volume 11 N°3, Pages 437– 451, 2008.
[9] M. Iqbal, An Introduction to Solar Radiation, Academic Press,
Department of Mechanical Engineering, University of British
Columbia, Canada, 1983.
[10] W. Spencer, Fourier series representation of the position of the Sun.
Search, Volume 2, Pages 172, 1971.
[11] R. Mérat et R. Moreau, physique appliquée, livre édition NATHAN
TECHNIQUE, 1994.
[12] Oar Aliman et al, Simplification of Sun Tracking Mode to Gain High
Concentration Solar Energy, American Journal of Applied Sciences,
Pages 167-171, 2007.
[13] W.B. Stine et R.W. Harrigan, Solar Energy Fundamentals and Design
with Computer Applications. John Wiley, New York, Pages 135-262,
1985.
[14] Y.T. Chen et al, Non-imaging focusing heliostat, Solar Energy,
Volume 71(3), Pages 155-164, 2001.
[15] A. Gamil et al, Simulation of Incident Solar Power Input to Fixed
Target Central Receiver System in Malaysia‖, 2013 IEEE Conference
on Sustainable Utilization and Development in Engineering and
Technology
[16] A. Montenon, Analyse, mutualisation et optimisation par la commande
de la consommation énergétique des héliostats autonomes des centrales
à concentration solaire, Thèse Doctorat, Université de Toulouse,
France, 2013.
[17] K. Gairaa, Valorisation du gisement solaire de Ghardaïa, Bulletin des
Energ. Ren., Algérie, N°17, 2010.

Copyright IPCO-2015
ISSN 2356-5608