MÁSTER EN ENERGÍAS RENOVABLES

Ingeniería de los Sistemas Fotovoltaicos

Radiación Solar

Antonio Urbina
 Sumario

1. El Sol y el espectro solar

2. Irradiancia e irradiación solar
3. Geometría solar
4. EU - JRC- PVGIS
 El Sol
El Sol es una estrella
El Sol es un reactor termonuclear de fusión
Composición: por # átomos por masa
Hidrogeno 91% 70.6%
Helio 8.9% 27.4%
Otros 0.1% 2%

Masa: 1.9891×1030 kg; 333,000 × Tierra

Radio: 695,700 km; 110 x Tierra

Temperatura: Centro (modelada): 1.57×107 K

Fotosfera (efectiva): 5,772 K
Corona: ≈ 5×106 K

Luminosidad: 3.828×1026 W Williams, D.R (2013) "Sun Fact Sheet“, NASA

Goddard Space Flight Center.
 La “constante” solar
La constante solar es la cantidad de potencia que el Sol deposita por
unidad de área expuesta a la luz solar. Su valor es aproximadamente
1367 W/m2 a la distancia de una Unidad Astronómica (AU*), que es
la distancia media Sol-Tierra:

𝐺𝑠𝑐 = 1367 𝑊/𝑚2

Variaciones del ciclo solar Variaciones estacionales

3.3%

360𝑛
(*) 1AU = 1.495978707×1011 m 𝐺𝑜𝑛 = 𝐺𝑠𝑐 1 + 0.033𝑐𝑜𝑠
365
(aprox. 150 millones de km)
𝑛 = día del año
 La “radiación de cuerpo negro”

Max Planck (1900): “lichtquantum”

Un cuerpo a temperatura T emite radiación

electromagnética con la siguiente distribución
espectral:

Otras expresiones útiles usando la frecuencia

angular y el momento (dualidad onda-partícula)

Donde: h = 6.62607015×10−34 J⋅s constante de Planck

kB = 1.380649×10−23 J⋅K−1 constante de Boltzmann p = momentum
c = 299792458 m⋅s-1 velocidad de la luz k = wave vector
 El espectro solar: irradiancia espectral
Aplicando el concepto de radiación de cuerpo negro al Sol (a T = 5772 K)

Considerando el tamaño angular del Radiación fuera de la atmósfera (WRC:

Sol visto desde la Tierra e integrando: World Radiation Centre standard, a 1UA)
 El espectro solar “AM1.5”
Cuando la radiación solar penetra en la atmósfera es atenuada. Se define la
“masa de aire” (Air Mass, AM) para calcular la atenuación:

AM0: fuera de la atmósfera

AM1: para z = 0

AM1.5: para z = 48.19° (STC)

 Sumario

1. El Sol y el espectro solar

2. Irradiancia e irradiación solar
3. Geometría solar
4. EU - JRC- PVGIS
 Irradiancia solar
Irradiancia (Irradiance, G): potencia por unidad de superficie (W/m2), se
obtiene integrando la irradiancia espectral a todas las longitudes de onda.
La irradiancia se define en función de la orientación del plano sobre el que
cae la radiación solar. Las dos definiciones más utilizadas son en el plano
perpendicular a la dirección del Sol o en plano horizontal (en la superficie
del suelo). Su cálculo dependerá de la geometría solar.

Irradiancia en Murcia 𝐺 = 𝐺𝑏 + 𝐺𝑑 + 𝐺𝑟
Consta de tres componentes principales:
• Irradiancia directa (“beam”, Gb), que incide
directamente desde el Sol
• Irradiancia difusa (“diffuse”, Gd) que incide
desde la bóveda celeste debido al scattering
de la luz solar por la atmósfera
• Albedo o irradiancia reflejada (“reflected”, Gr)
que incide desde la reflexión en el paisaje,
edificios, suelo, etc…
 Irradiación solar
Irradiación (Irradiation, H, I, E): energía por unidad de superficie (J/m2,
kWh/m2). Normalmente se expresa en valores medios para un período de
tiempo determinado: hora, día, mes o año. Se obtiene integrando la curva
de irradiancia sobre dicho periodo de tiempo.
Insolación (Insolation, I) es el término que se dá a la irradiación proveniente del Sol.
 Sumario

1. El Sol y el espectro solar

2. Irradiancia e irradiación solar
3. Geometría solar
4. EU - JRC- PVGIS
 Coordenadas celestes: “horizontales” alt-az
El horizonte se establece como “círculo fundamental”,
y su polo se establece con la vertical local: el zenith (y
el nadir). Las coordenadas son:

altitud (a) del objeto X es la distancia angular a lo

largo del círculo vertical trazado desde el
horizonte hasta X (+ 90° a - 90°).

azimuth (A) del objeto X es la distancia angular a

lo largo del horizonte desde el punto cardinal N (a
veces S) hasta el círculo vertical que pasa por X
(0°a 360°).
Relación con las coordenadas terrestres de la localidad geográfica del observador:
La altitud (a) del polo norte celestial coincide con la latitud () de la localidad geográfica.
Para obtener la longitud, es necesario tener en cuenta la rotación angular (en horas)
de la Tierra y las definiciones del tiempo “sideral” del meridiano de Greenwich (GST),
que se toma como referencia y el tiempo sideral local (LST, t) del observador:
LST = GST – longitud (Oeste)
Fuente de los gráficos: Univ. St. Andrews (Scotland) Prof. Fiona Vincent
Para detalles sobre LST y rotación terrestre, ver USNO Circular No.179
 Coordenadas celestes: ecuatoriales I y II
I) HA-dec II) RA-dec

I y II) Declinación (dec, d): es la distancia angular desde el ecuador celeste hasta el objeto X
I) Angulo horario (HA, H): distancia angular II) Ascensión recta (RA, a): es el ángulo
entre “el” meridiano celeste y el meridiano medido a lo largo del ecuador celeste entre el
de X; se mide hacia el oeste en tiempo de 0h meridiano de X y el meridiano que pasa por el
a 24h (la tierra rota 360o en 24h) punto del equinocio de primavera ()
Un objeto en “el” meridiano tiene H = 0; y en el tránsito (RA se mide también en horas, pero en sentido contrario a
inferior (posiblemente bajo el horizonte) tiene H = 12 HA, y tiene la ventaja de no depender del tiempo local
[1 hora = 15o] sideral LST, t)
 Posición solar en la esfera celeste
Posición del Sol en la esfera celeste (en coordenadas celestiales) depende del tiempo
(fecha y hora) y la localidad geográfica que se considera
Declinación (d): altitud del Sol cuando está en el meridiano (a mediodía)
Ascensión recta (RA): tiempo local sideral (LST) cuando el Sol cruza el meridiano.

El camino “aparente” que sigue el Sol es la eclíptica

Septiembre Junio Marzo Diciembre Septiembre

Septiembre (E)
Julio (A) +23°26’ > d > -23°26'
E: Equinocio
Enero (P) S: Solsticio
A: Afelio
Diciembre (S) Junio (S) P: Perihelio
Marzo (E)
Conversión entre coordenadas ecuatoriales y horizontales

Utilizando el triángulo XPZ, llamado

“triángulo astronómico” y aplicando
trigonometría esférica se obtienen las
ecuaciones de conversión.

Además tenemos, por definición:

LST (t) = GST - longitud oeste
HA (H) = LST (t) - RA (a)

H=t-α
ecuatorial → horizontal sin(a) = sin(δ) sin(φ) + cos(δ) cos(φ) cos(H)
(H, d) → (A, a) sin(A) = - sin(H) cos(δ) / cos(a)
cos(A) = { sin(δ) - sin(φ) sin(a) } / cos(φ) cos(a)

sin(δ) = sin(a)sin(φ) + cos(a) cos(φ) cos(A)

horizontal → ecuatorial sin(H) = - sin(A) cos(a) / cos(δ)
(A, a) → (H, d) cos(H) = { sin(a) - sin(δ) sin(φ)} / cos(δ) cos(φ)
α=t–H
 Relación entre coordenadas eclípticas y ecuatoriales

Las coordenadas eclípticas son muy

útiles para los objetos del sistema solar.
Para el resto de objetos se utilizan las
coordenadas ecuatoriales.

longitud eclíptica (o celeste) de X: b

latitud eclíptica (o celeste) de X: l

Aplicando trigonometría al triángulo KPX

se convierten las coordenadas eclípticas
en ecuatoriales:

sin(δ) = sin(β) cos(ε) + cos(β) sin(ε) sin(λ)

ecuatorial  eclíptica
sin(β) = sin(δ) cos(ε) - cos(δ) sin(ε) sin(α)
(H, d)  (b, l) cos(λ) cos(β) = cos(α) cos(δ)

e es la “oblicuidad” de la elíptica: e = 23°26’ que es la inclinación del eje de rotación

terrestre (diario) respecto al plano de rotación de la Tierra alrededor del Sol.
 Proyección de la irradiancia sobre superficie horizontal
Máxima irradiancia solar sobre una superficie horizontal (en función de la posición
solar en la esfera celeste y las coordenadas geográficas de la localidad :

Siento w el ángulo horario solar. Debido a las leyes de Kepler (que describen las
órbitas planetarias), la hora solar varía a lo largo del año: el Sol se mueve más len-
tamente en el perihelio y más rápido en el afelio debido a la excentricidad de la órbita
elíptica. Además se debe considerar la oblicuidad de la órbita (n es el día del año).

La solución del sistema de ecuaciones

para obtener w no es analítica, y un
cálculo numérico nos dá el siguiente
resultado si se compara la hora solar
real (“true”) y una hora solar media
(calculada como si el Sol se moviera
uniformemente). Esta diferencia es la
“ecuación del tiempo”:
 Cartas de posición solar
La solución numérica de las ecuaciones se suministra en tablas de datos o cartas
que permiten obtener la posición del Sol en cualquier fecha/hora y en cualquier
lugar geográfico (cada latitud tendrá su propia carta)

La proyección de la topografía o edificios permite calcular el efecto de las sombras

 Orientación del panel solar

Además de los ángulos que definen:

Sol
1) La posición del Sol en los distintos
sistemas de coordenadas (horizontal,
ecuatorial, eclíptico) al que se añade la
definición del “ángulo zenital”: z = 90o – a

2) Las coordenadas geográficas del lugar:

latitud () y longitud.

Se añaden ahora los ángulos que definen la

orientación del panel solar:

Azimuth del panel (g): ángulo que forma la proyección sobre el plano horizontal de
la perpendicular a la superficie del panel con la dirección Sur. g = 0 para orientación
Sur, positivo hacia el Oeste y negativo hacia el Este (en el hemisferio Norte).

Inclinación del panel (b): ángulo que forma la superficie del panel con el plano
horizontal. b = 0 si el panel se coloca horizontal y b = 90o si se coloca vertical.
 Proyección de la irradiancia sobre el panel solar (I)

Cada componente de la irradiancia solar (directa, difusa, albedo) se proyecta de

manera diferente sobre la superficie del panel, según los factores Ri

𝐺 (𝛾, 𝛽) = 𝐺𝑏,0 ∙ 𝑅𝑏 + 𝐺𝑑,0 ∙ 𝑅𝑑 + 𝐺𝑟,0 ∙ 𝑅𝑟

Se incluye un nuevo ángulo i (*)

Coeficiente de radiación directa (beam, Rb): es la proyección sobre la

superficie inclinada
 Proyección de la irradiancia sobre el panel solar (II)

Coeficiente de radiación difusa (Rd):

Algunos modelos asumen que la intensidad de la variacion difusa del cielo es
uniforme, independientemente del sector del cielo del cual provenga (modelos
isotropos). Uno de los mas famosos es el modelo de Liu-Jordan (1960).

Otros modelos proponen que la radiacion difusa esta formada por 2 componentes:
una parte isotropica y una circunsolar (modelos anisótropos de 2 componentes):
Muneer (1990), Hay (1993).

Otros modelos además consideran una tercera componente, conocida como brillo del
horizonte (modelos anisótropos de 3 componentes): Pérez (1990).

El índice de claridad (“clearness index” KT , kT) es la relación entre la irradiancia

global horizontal G0 (kW/m2 calculando el valor medio por hora, día o mes) y la
irradiancia correspondiente disponible fuera de la atmósfera (varía entre 0 y 1, no
tiene unidades); se suele calcular para un promedio diario o mensual (K mayúscula)
u horario (k minúscula). Valores típicos entre 0.25 y 0.75.
 Proyección de la irradiancia sobre el panel solar (III)

Coeficiente de radiación reflejada (albedo, Rr):

Los modelos más sencillos asignan una dependencia de la inclinación del panel
(b) pues asumen que la radiación difusa es isotrópica. Además se incluye un valor
para la reflectancia o albedo del suelo ()
Material Albedo ()
Asfalto nuevo 0.04
Océano 0.06
Asfalto gastado 0.12
Bosque de coníferas 0.08 a 0.15
Bosque caduco 0.15 a 0.18
Suelo rural 0.17
Hierba verde 0.25
Desierto de arena 0.40
Cemento nuevo 0.55
Albedo superficial de la Tierra, datos del Hielo sobre el mar 0.5 a 0.7
Copernicus Global Land Service (EU) Nieve fresca 0.8
 Medida de la irradiancia
La medida de la irradiancia se realiza mediante un piranómetro. Son dispositivos
que convierten la irradiancia solar (fotones) en una señal eléctrica (electrones).
Dos tipos principales:
• Termopilas (basadas en el efecto termoeléctrico)
• Fotodiodos o células calibradas (basados en el efecto fotovoltaico)

Mediante su diseño permiten obtener datos de irradiancia

global, directa* y difusa (incluyendo albedo en este caso).
Algunos modelos avanzados permiten distinguir zonas del
espectro (fotodiodos)

Es importante realizar una calibración de forma periódica.

Su uso está regulado por el standard
IEC 61724-1:2017
(*) en este caso se llama pirheliómetro
 Sumario

1. El Sol y el espectro solar

2. Irradiancia e irradiación solar
3. Geometría solar
4. EU - JRC- PVGIS
 EU - JRC- PVGIS
Para obtener datos de irradiancia en cualquier localidad geográfica y para cualquier
orientación de panel se pueden realizar los cálculos descritos hasta ahora. Pero
también se puede acudir a base de datos basadas en observaciones de satélite,
complementadas con modelos y procedimientos de interpolación de datos. La más
utilizada es:
“Photovoltaic geographical information system” (PVGIS)
del Joint Research Centre (JCR) de la Comisión Europea (Unión Europea, EU)

Cálculo de producción
fotovoltaica (varias
tecnologías y tipos de
sistema)

Base de datos de
radiación solar (horario,
diario y mensual para
distintas orientaciones e
inclinaciones del panel)

Datos meteorológicos
¡Gracias por vuestra atención!