Méthode de calcul du rayonnement solaire
Disponible avec une licence Spatial Analyst.
Les outils d’analyse du rayonnement solaire permettent de calculer l’insolation d’un paysage ou de localisations spécifiques, en s’appuyant sur les méthodes de l’algorithme du champ de vision hémisphérique développées par Rich et al. (Rich 1990, Rich et al. 1994) et reprises ultérieurement par Fu et Rich (2000, 2002).
La quantité totale de rayonnement calculée pour une localisation ou une zone en particulier est entrée sous la forme de rayonnement global. Le calcul de l’insolation directe, diffuse et globale est répété pour chaque localisation d’entité ou pour toutes les localisations de la surface topographique et génère des cartes d’insolation pour toute une zone géographique.
Les sections suivantes présentent les équations utilisées pour calculer le rayonnement solaire.
Calcul du rayonnement global
Le rayonnement global (\(Global_{tot}\)) est calculé comme étant la somme des rayonnements direct (\(Dir_{tot}\)) et diffus (\(Dif_{tot}\)), respectivement, de tous les secteurs de la carte d’ensoleillement et de la carte du ciel.
\(Global_{tot} = Dir_{tot} + Dif_{tot}\)
Rayonnement solaire direct
L’insolation directe totale (\(Dir_{tot}\)) d’une localisation donnée est la somme de l’insolation directe (\(Dir_{θ,α}\)) de tous les secteurs de la carte d’ensoleillement :
\(\large Dir_{tot} = Σ \space Dir_{θ,α} \space \space \space \normalsize (1)\)
L’insolation directe du secteur de la carte d’ensoleillement (\(Dir_{θ,α}\)) avec un centroïde à l’angle zénithal (\(θ\)) et à l’angle azimutal (\(α\)) est calculé à l’aide de l’équation suivante :
\(\large Dir_{θ,α} = S_{Const} * β ^{m(θ)} * SunDur_{θ,α} * SunGap_{θ,α} * cos(AngIn_{θ,α} ) \space \space \space \normalsize (2)\)
où :
\(S_{Const}\) représente le flux solaire à l’extérieur de l’atmosphère à la distance moyenne séparant la Terre du Soleil, appelé constante solaire. La constante solaire utilisée dans l’analyse est 1 367 W/m2. Cela est cohérent avec la constante solaire du World Radiation Center (WRC).
β représente la transmissivité de l’atmosphère (moyennée sur toutes les longueurs d’onde) pour le chemin le plus court (en direction du zénith).
\(m(θ)\) représente la longueur du chemin optique relative, mesurée comme étant une proportion relative à la longueur du chemin zénithal (voir l’équation 3 ci-dessous).
\(SunDur_{θ,α}\) correspond à la durée représentée par le secteur céleste. Pour la plupart des secteurs, cette valeur est égale à l’intervalle de jours (par exemple, un mois) multiplié par l’intervalle d’heures (par exemple, une demi-heure). En ce qui concerne les secteurs partiels (proches de l’horizon), la durée est calculée à l’aide de la géométrie sphérique.
\(SunGap_{θ,α}\) représente la fraction de trouée pour le secteur de la carte d’ensoleillement.
\(AngIn_{θ,α}\) représente l’angle d’incidence entre le centroïde du secteur céleste et l’axe normal par rapport à la surface (voir l’équation 4 ci-dessous).
La longueur optique relative, \(m(θ)\), est déterminée par l’angle zénithal solaire et l’élévation au-dessus du niveau de la mer. Pour les angles zénithaux inférieurs à 80°, elle peut être calculée à l’aide de l’équation suivante :
\(m(θ) = \Large \frac {EXP (-0.000118 \space * \space Elev \space - 1.638 * 10 ^{-9} \space * \space Elev ^{2} )}{ cos(θ)} \space \space \space \normalsize (3)\)
où :
\(θ\) représente l’angle zénithal solaire.
\(Elev\) représente l’élévation au-dessus du niveau de la mer, exprimée en mètres.
L’effet de l’orientation de la surface est pris en compte par la multiplication par le cosinus de l’angle d’incidence. L’angle d’incidence (\(AngInSky_{θ,α}\)) entre la surface d’interception et un secteur céleste donné ayant un centroïde à l’angle zénithal et à l’angle azimutal est calculé à l’aide de l’équation suivante :
\(\large AngIn_{θ,α} = acos( Cos(θ) * Cos(G_{z}) + Sin(θ) * Sin(G_{z}) * Cos(α-G_{a})) \space \space \space \normalsize (4)\)
où :
\(G_z\) représente l’angle zénithal de la surface. (Il convient de noter que pour les angles zénithaux supérieurs à 80°, la réfraction est importante.)
\(G_a\) représente l’angle azimutal de la surface.
Calcul du rayonnement diffus
Pour chaque secteur céleste, le rayonnement diffus à son centroïde (\(Dif\)) est calculé, intégré sur l’intervalle temporel, puis corrigé par la fraction de trouée et l’angle d’incidence à l’aide de l’équation suivante :
\(\large Dif_{θ,α} = R_{glb} * P_{dif} * Dur * SkyGap_{θ,α} * Weight_{θ,α} * cos(AngIn_{θ,α}) \space \space \space \normalsize (5)\)
où :
\(\large R_{glb}\) représente le rayonnement normal global (voir l’équation 6 ci-dessous).
\(\large P_{dif}\) représente la proportion du flux de rayonnement normal global qui est diffus. Elle est généralement égale à 0,2 pour des conditions de ciel dégagé et à 0,7 pour des conditions de ciel très couvert.
\(\large Dur\) représente l’intervalle temporel de l’analyse.
\(\large SkyGap_{θ,α}\) représente la fraction de trouée (proportion de ciel visible) pour le secteur céleste.
\(\large Weight_{θ,α}\) représente la proportion du rayonnement diffus émis par un secteur céleste donné par rapport à tous les secteurs (voir les équations 7 et 8 ci-dessous).
\(\large AngIn_{θ,α}\) représente l’angle d’incidence entre le centroïde du secteur céleste et la surface d’interception.
Le rayonnement normal global (\(R_{glb}\)) peut être calculé en additionnant le rayonnement direct de tous les secteurs (y compris les secteurs obstrués) sans correction de l’angle d’incidence, puis avec correction pour la proportion du rayonnement direct, qui est égale à \(1-P_{dif}\) :
\(\large R_{glb} = (S_{Const} \space Σ(β ^{m(θ)})) / (1 - P_{dif}) \space \space \space \normalsize (6)\)
Pour le modèle diffus de ciel uniforme, la valeur \(Weight_{θ,α}\) est calculée comme suit :
\(\large Weight_{θ,α} = (cosθ_{2} - cosθ_{1}) / Div_{azi} \space \space \space \normalsize (7)\)
où :
\(θ_1\) et \(θ_2\) représentent les angles zénithaux limites du secteur céleste.
\(Div_{azi}\) représente le nombre de divisions azimutales de la carte du ciel.
Pour le modèle de ciel couvert standard, la valeur \(Weight_{θ,α}\) est calculée comme suit :
\(\large Weight_{θ,α} = (2 cosθ_{2} + cos2θ_{2} - 2 cosθ_{1} - cos 2θ_{1}) / 4 \space Div_{azi} \space \space \space \normalsize (8)\)
Rayonnement solaire diffus total pour la localisation ($Dif_) is calculated as the sum of the diffuse solar radiation (\(Dif\)) de tous les secteurs de la carte du ciel :
\(\large Dif_{tot} = Σ \space Dif_{θ,α} \space \space \space \normalsize (9)\)
Bibliographie
Fu, P. 2000. A Geometric Solar Radiation Model with Applications in Landscape Ecology. Thèse de doctorat, Département de géographie, Université du Kansas, Lawrence, Kansas, États-Unis.
Fu, P. et P. M. Rich. 2000. The Solar Analyst 1.0 Manual. Helios Environmental Modeling Institute (HEMI), États-Unis.
Fu, P. et P. M. Rich. 2002. "A Geometric Solar Radiation Model with Applications in Agriculture and Forestry." Computers and Electronics in Agriculture 37:25–35.
Rich, P. M., R. Dubayah, W. A. Hetrick et S. C. Saving. 1994. "Using Viewshed Models to Calculate Intercepted Solar Radiation: Applications in Ecology. American Society for Photogrammetry and Remote Sensing Technical Papers, 524–529.
Rich, P. M. et P. Fu. 2000. "Topoclimatic Habitat Models." Proceedings of the Fourth International Conference on Integrating GIS and Environmental Modeling.