Las herramientas de análisis de radiación solar calculan la irradiación solar en un paisaje o correspondiente a ubicaciones específicas a partir de los métodos del algoritmo de cuenca visual hemisférica desarrollado por Rich et al. (Rich 1990, Rich et al. 1994) y posteriormente por Fu y Rich (2000, 2002).
La cantidad total de radiación calculada para un área o una ubicación se proporciona como radiación total. Los cálculos de irradiación solar directa, difusa y global se repiten para cada ubicación de entidad o ubicación en la superficie topográfica y así se producen mapas de insolación para un área geográfica completa.
Ecuaciones de radiación solar
Cálculo de radiación global
La radiación global (Globaltot) se calcula como la suma de la radiación directa (Dirtot) y difusa (Diftot) de todos los sectores del mapa solar y el mapa celeste, respectivamente.
Globaltot = Dirtot + DiftotRadiación solar directa
La insolación directa total (Dirtot) en una ubicación determinada es la suma de la insolación directa (Dirθ,α) desde todos los sectores del mapa solar:
Dirtot = Σ Dirθ,α (1)La insolación directa desde el sector del mapa solar (Dirθ,α) con un centroide en el ángulo cenital (θ) y el ángulo acimutal (α) se calcula con la siguiente ecuación:
Dirθ,α = SConst * βm(θ) * SunDurθ,α * SunGapθ,α * cos(AngInθ,α) (2)- donde:
- SConst: flujo solar en la parte externa de la atmósfera a la distancia media entre el sol y la tierra, conocida como constante solar. La constante solar utilizada en el análisis es 1367 W/m2. Coincide con la constante solar del Centro de Radiación Mundial (WRC).
- β: transmisividad de la atmósfera (como media de todas las longitudes de onda) correspondiente a la trayectoria más corta (en la dirección del cénit).
- m(θ): longitud de trayectoria óptica relativa, medida como una proporción relativa la longitud de trayectoria cenital (consulte la ecuación 3 siguiente).
- SunDurθ,α: duración en el tiempo representada por el sector celeste. Para la mayoría de sectores, es igual que el intervalo de días (por ejemplo, un mes) multiplicado por el intervalo de horas (por ejemplo, media hora). Para sectores parciales (cerca del horizonte), la duración se calcula utilizando la geometría esférica.
- SunGapθ,α: fracción de huecos correspondiente al sector del mapa solar.
- AngInθ,α: ángulo de incidencia entre el centroide del sector celeste y el eje normal a la superficie (consulte la ecuación 4 a continuación).
La longitud óptica relativa, m(θ), está determinada por el ángulo cenital del sol y la elevación sobre el nivel del mar. En el caso de los ángulos cenitales de menos de 80°, puede calcularse con la ecuación siguiente:
m(θ) = EXP(-0.000118 * Elev - 1.638*10-9 * Elev2) / cos(θ) (3)- donde:
- θ: ángulo cenital del sol.
- Elev: elevación sobre el nivel del mar en metros.
El efecto de la orientación de la superficie se tiene en cuenta multiplicando por el coseno del ángulo de incidencia. El ángulo de incidencia (AngInSkyθ,α) entre la superficie de captación y un sector celeste determinado con un centroide en el ángulo cenital y el ángulo acimutal se calcula con la ecuación siguiente:
AngInθ,α = acos( Cos(θ) * Cos(Gz) + Sin(θ) * Sin(Gz) * Cos(α-Ga) ) (4)- donde:
- Gz: ángulo cenital de la superficie.
La refracción es importante en el caso de los ángulos cenitales de más de 80°.
- Ga: ángulo acimutal de la superficie.
- Gz: ángulo cenital de la superficie.
Cálculo de radiación difusa
Para cada sector celeste se calcula la radiación difusa en su centroide (Dif), integrada a lo largo del intervalo de tiempo y corregida por medio de la fracción de huecos y el ángulo de incidencia mediante el uso de la ecuación siguiente:
Difθ,α = Rglb * Pdif * Dur * SkyGapθ,α * Weightθ,α * cos(AngInθ,α) (5)- donde:
- Rglb: radiación normal global (consulte la ecuación 6 a continuación).
- Pdif: proporción del flujo de radiación normal global que es difusa. Normalmente, es aproximadamente 0,2 para las condiciones de cielo muy despejado y 0,7 para las condiciones con cielo muy nublado.
- Dur: intervalo de tiempo para el análisis.
- SkyGapθ,α: fracción de huecos (proporción de cielo visible) del sector celeste.
- Weightθ,α: proporción de radiación difusa que se origina en un sector celeste determinado relativa a todos los sectores (consulte las ecuaciones 7 y 8 siguientes).
- AngInθ,α: ángulo de incidencia entre el centroide del sector celeste y la superficie de captación.
La radiación normal global (Rglb) puede calcularse sumando la radiación directa de cada sector (incluidos los sectores tapados) sin corrección del ángulo de incidencia y corrigiendo la proporción de radiación directa, lo que equivale a 1-Pdif:
Rglb = (SConst Σ(βm(θ))) / (1 - Pdif) (6)En el caso del modelo de difusión celeste uniforme, Weightθ,α se calcula como sigue:
Weightθ,α = (cosθ2- cosθ1) / Divazi (7)- donde:
- θ1 and θ2: ángulos cenitales delimitadores del sector celeste.
- Divazi: número de divisiones acimutales del mapa celeste.
En el caso del modelo de cielo cubierto estándar, Weightθ,α se calcula como sigue:
Weightθ,α = (2cosθ2 + cos2θ2 - 2cosθ1 - cos2θ1) / 4 * Divazi (8)La radiación solar difusa total de la ubicación (Diftot) se calcula como la suma de la radiación solar difusa (Dif) de todos los sectores del mapa celeste:
Diftot = Σ Difθ,α (9)Referencias
Fu, P. 2000. A Geometric Solar Radiation Model with Applications in Landscape Ecology. Ph.D. Thesis, Department of Geography, University of Kansas, Lawrence, Kansas, USA.
Fu, P., and P. M. Rich. 2000. The Solar Analyst 1.0 Manual. Helios Environmental Modeling Institute (HEMI), USA.
Fu, P., and 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, and 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., and P. Fu. 2000. "Topoclimatic Habitat Models." Proceedings of the Fourth International Conference on Integrating GIS and Environmental Modeling.