Ingeniare. Revista chilena de ingeniería, vol.

 30 Nº 1, 2022, pp. 69-81

Determinación y análisis temporal de la radiación solar

global en el Altiplano de Puno

Determination and temporal analysis of global solar

radiation in the Puno Altiplano

Lelia Quispe Huamán1*  Guina Sotomayor Alzamora2

Recibido 14 de septiembre de 2020, aceptado 23 de diciembre de 2021

Received: September 14, 2020   Accepted: December 23, 2021

RESUMEN

El sol proporciona una fuente estable y confortable de calor, luz y es la fuente central de energía renovable.
Su estudio es vital para evaluar su disponibilidad y utilidad en diversas aplicaciones para establecer
recursos de energía solar utilizando energía limpia y gratuita. Este trabajo analiza la disponibilidad de
radiación solar global en tres zonas ecológicas del altiplano de Puno como Circunlacustre, Puna Húmeda
y Clima de Altura; involucrando catorce estaciones meteorológicas del SENAMHI desde 2007 al 2013.
El modelo utilizado es el de Bristow y Campbell, que es simple, confiable y efectivo, y utiliza datos
que se encuentran fácilmente en las estaciones meteorológicas. Como resultado, se ha encontrado una
correlación de 0,938, entre la radiación solar global estimada y medida, además, se ha observado que la
radiación solar global tiene comportamientos estacionales durante los años analizados. Posteriormente,
se obtuvieron los promedios anuales máximos y mínimos de las zonas evaluadas. Los valores son 6,342
y 5,481; 6,783 y 5,908; 7,298 y 5,945 kWh.m–2.dia–1, respectivamente. Se concluye que los resultados
en el área de estudio son apropiados y muy rentables para aplicaciones de energía solar debido a que son
superiores a 4,000 kWh.m–2.dia–1 en el período analizado. Finalmente, se busca que la radiación solar
global se transforme en energía térmica y eléctrica, para ayudar a solucionar problemas como las heladas,
que afronta la población de la zona altoandina de la región Puno.

Palabras clave: Altiplano de Puno, Bristow y Campbell, variabilidad de la energía solar, radiación solar
global, amplitud térmica, transmitancia atmosférica máxima.

ABSTRACT

The sun provides a stable and comfortable source of heat, light and is the central source of renewable
energy. Its study is vital to evaluate its availability and usefulness in various applications to establish
solar energy resources using clean and free energy. This work analyzes the availability of global solar
radiation in three ecological zones of the high plateau of Puno as Circunlacustre, Puna Húmeda and
Climate of Height; involving fourteen SENAMHI meteorological stations from 2007 to 2013. The model
used is that of Bristow and Campbell, which is simple, reliable and effective, and uses data that are easily
found in meteorological stations. As a result, a correlation of 0,938 was found between the estimated and

1 Universidad Nacional del Altiplano. Escuela Profesional de Ciencias Físico Matemáticas. Puno, Perú.
E-mail: [email protected]
2 Universidad Nacional del Altiplano. Escuela Profesional de Ingeniería de Sistemas. Puno, Perú.

E-mail: [email protected]
* Autor de correspondencia: [email protected]
Ingeniare. Revista chilena de ingeniería, vol. 30 Nº 1, 2022

measured global solar radiation, observing the global solar radiation has seasonal behaviors during the
years analyzed. Subsequently, the maximum and minimum annual averages of the evaluated areas were
obtained. The values are 6,342 and 5,481; 6,783 and 5,908; 7,298 and 5,945 kWh.m–2.day–1, respectively.
It is concluded that the results in the study area are appropriate and very profitable for solar energy
applications because they are higher than 4,000 kWh.m–2.day–1 in the period analyzed. Finally, it is sought
that global solar radiation is transformed into thermal and electrical energy, to help solve problems such
as frost, faced by the population of the high Andean area of the Puno region.

Keywords: Altiplano de Puno, Bristow and Campbell, variability of solar energy, global solar radiation,
thermal amplitude, maximum atmospheric transmittance.

INTRODUCCIÓN En zonas altoandinas del Perú, no existen muchas

estaciones meteorológicas del SENAMHI que
La fuente de energía para la vida proviene del sol, recolectan datos para que se puedan realizar mediciones
que proporciona una fuente estable y confortable de radiación solar, por lo que la información existente
de calor y luz [1]. La energía recibida por el puede no ser de fácil acceso o puede estar incompleta
sistema de la tierra y atmósfera es también crucial [13]. En este escenario, surgen los modelos empíricos
para equilibrar la energía en nuestro planeta [2]; [14, 15, 16], que sirven como herramientas esenciales y
además, es la principal fuente de energía renovable económicas, pues dependen de datos básicos, como el
que soporta la biósfera y estimula los procesos geoposicionamiento y las variables meteorológicas [17,
físicos, químicos y biológicos en la tierra, tales 18, 19]. Por ejemplo, los modelos Bristow y Campbell
como los procesos ecológicos y biofísicos [3], y Hargreaves-Samani [20] utilizan la temperatura del
como la fotosíntesis y los procesos ambientales aire, ubicación, el día del año, amplitud térmica y las
relacionados a la temperatura del aire y suelo, condiciones atmosféricas para estimar la radiación
afectando la transferencia de calor por medio de solar global [21], así estos modelos empíricos pueden
la evaporación y transpiración [4]. Sin embargo, ser usados como predictores [22].
no toda la radiación solar consigue incidir a la
superficie terrestre, esto se debe a que interactúa Dada la importancia que tiene el estudio y aplicación de
con componentes atmosféricos al atravesar la la radiación solar en proyectos energéticos, se presenta
masa atmosférica y sufre la acción de procesos el Modelo Bristow y Campbell para determinar la
de reflexión, dispersión y absorción [5]. Así, el radiación solar global y analizar la disponibilidad
51% de la radiación extraterrestre está disponible en tres zonas altoandinas, ubicadas dentro del
para procesos en medios físicos y biológicos, la departamento de Puno: La zona Circunlacustre, la
fracción que efectivamente llega a la superficie zona de Puna Húmeda y la zona de Clima de Altura,
terrestre se denomina radiación solar global [6]. durante el periodo de 2007 a 2013.

El conocimiento preciso de la radiación solar MATERIALES Y MÉTODOS

permite evaluar su disponibilidad y utilidad en
aplicaciones de diversas áreas, como la meteorología, Área de estudio y datos utilizados
climatología, hidrología, agricultura, entre otros La presente investigación se realizó en el departamento
[7, 8, 9]. Además, los datos de la radiación solar de Puno, específicamente en el altiplano Puneño,
global son útiles para poder predecir la eficiencia ubicado a más de 3810 m.s.n.m. y abarcó tres zonas
y optimización del diseño de dispositivos solares específicas: La zona Circunlacustre, la zona de Puna
como los colectores solares térmicos y sistemas Húmeda y la zona de Clima de Altura, considerando
fotovoltaicos [10, 11]. Predecir la radiación solar las estaciones que se encuentran dentro de cada
global es fundamental para establecer la eficacia zona [23]. Algunas especificaciones se detallan a
de los recursos de la energía solar, haciendo uso de continuación:
energía libre y limpia, identificando y evaluando el
uso de energía solar en diferentes lugares alrededor • Zona Circunlacustre, corresponde al área
del mundo [12]. que bordea al lago Titicaca, comprende a las

70
 Quispe y Sotomayor: Determinación y análisis temporal de la radiación solar global en el Altiplano de Puno

estaciones de Puno, Juli, Capachica, Moho y efectivo [3]. Su uso se simplifica pues utiliza
y Yunguyo. En esta zona ocurre un efecto datos comunes, los cuales pueden ser fácilmente
termorregulador, que consiste en que el lago encontrados en las estaciones meteorológicas del
absorbe el calor del sol durante el día y lo pierde SENAMHI.
de forma lenta durante la noche, debido a este
efecto, disminuye fuertemente la amplitud Cuando la radiación solar llega a la superficie
térmica. terrestre, una parte de ella se refleja y la otra se
• Zona Puna Húmeda, comprende a las estaciones absorbe. Además, al calentarse el suelo, se emite
de Azángaro, Crucero, Desaguadero, Huancané, una radiación de onda larga, que calienta el aire
Ilave y Putina. En esta zona la temperatura sufre adyacente, contribuyendo a cambiar su temperatura.
fuertes descensos, debido a la presencia de Basados en ese argumento, Bristow y Campbell
heladas más o menos intensas, específicamente presentaron un modelo empírico para estimar
en estaciones de invierno. la radiación solar global diaria, en función de la
• Zona Clima de Altura, comprende a las radiación extraterrestre y la diferencia diaria entre
estaciones de Ananea, Capazo y Mazocruz. En las temperaturas máximas y mínimas [26], por esta
esta zona de altura, la topografía es sumamente razón, este modelo depende de la relación entre el
accidentada, con suelos erosionados y vegetación calor sensible y el calor latente. El calor sensible
natural muy pobre. Existe poca información depende de la radiación solar diaria entrante y
meteorológica; sin embargo, se estima que es responsable de las temperaturas máximas del
la temperatura promedio es sumamente baja, aire. Por la noche, el calor sensible se pierde en el
clara indicación de heladas intensas frecuentes espacio, en forma de radiación de onda larga, junto
durante todo el año. con los flujos radiantes, esto da lugar a un descenso
de la temperatura del aire, hasta que se alcanza la
Los datos utilizados para la determinación de la temperatura mínima diaria, normalmente justo
radiación solar global han sido medidos por las antes del amanecer [27].
estaciones meteorológicas del Servicio Nacional
de Meteorología e Hidrología (SENAMHI). El La radiación solar global diaria es estimada como
periodo de estudio considerado fue de enero del una función de radiación extraterrestre y la diferencia
2007 a diciembre de 2013, con intervalos de 15 entre las temperaturas máxima y mínima diarias,
minutos. Como excepción a la toma de datos, los como se muestra en la ecuación (1), dónde Rg es la
correspondientes a la estación de Puno corresponden radiación solar global (kWh.m–2), Re la radiación solar
al periodo de enero del 2007 a diciembre del 2017. extraterrestre (kWh.m–2), Tmax la temperatura máxima
Todos los datos fueron almacenados en la base de (ºC), Tmin la temperatura mínima (ºC). Además, se
datos que proporciona el software Weatherlink. utilizan los coeficientes empíricos con significado
Además, para la validación, se consideraron los datos físico A, B y C, siendo A el valor máximo de la
de radiación solar global medidos con el Piranómetro transmitancia atmosférica característica a cada área
de la marca Kipp-Zonen, específicamente utilizando de estudio, depende de la contaminación y altitud
el modelo de clase CMP11, que es un piranómetro de la zona, los coeficientes B (ºC –1) y C determinan
estándar secundario ISO, el cual posee un nivel de el efecto del incremento en la amplitud térmica
sensibilidad de 7,000 a 14,000 µV.W–1.m –2. (representada también como ΔT) en los máximos
valores de transmitancia atmosférica (At) [26], vale
Modelo Bristow y Campbell decir que estos últimos coeficientes dependen de
El costo de realizar mediciones de radiación solar la amplitud térmica y latitud correspondiente a la
global es alto, además se requiere de equipos zona evaluada.
altamente especializados y un tiempo extenso de
procesamiento, por esta razón datos de la radiación Rg ⎡ −B T −T
c
⎤
solar global no están disponibles en la mayoría de las = A ⎢1− e ( max min ) ⎥ (1)
Re ⎣ ⎦
regiones del mundo [24, 25]. Sin embargo, se puede
estimar la radiación solar global utilizando métodos
empíricos como el modelo de Bristow y Campbell, Los coeficientes B y C propuestos, pueden ser
el cual es un método simple, rápido, confiable aplicados solamente en áreas con similares

71
Ingeniare. Revista chilena de ingeniería, vol. 30 Nº 1, 2022

condiciones de régimen térmico, esto debido a Otro trabajo más reciente, utilizando el modelo
que los valores absolutos y las diferencias entre las Bristow y Campbell fue presentado por Camayo
temperaturas máximas y mínimas son fuertemente et al. [28], quienes comprobaron la validación del
influenciados por la topografía, latitud y altitud, modelo, a través de la radiación solar global media
entre otros factores. Las ecuaciones empíricas, mensual en Junín, cuyos resultados indican que
ecuación (2) y ecuación (3) determinan el valor de no existen mayores diferencias entre los valores
tales coeficientes [28], donde φ representa a la latitud. registrados y los calculados por el modelo aplicado.
Así mismo, concluyeron que para que una región
sea económicamente rentable, la energía solar debe
C = 2,116 − 0,072 ( Tmax − Tmin ) + 57,574e φ (2)
ser mayor que 4,000 kWh.m –2 por día. Los autores
consideraron el valor empírico A de la estación
B = 0,107C−2,6485 (3)
Huayao del Atlas de Energía Solar [34], cabe
mencionar que este Atlas informa también sobre
Para la aplicación de este modelo, es necesario el fenómeno del efecto termorregulador de Puno,
contar con las coordenadas geográficas de las zonas por estar a orillas del lago, evitando que valores de
de estudio, representadas por la latitud φ, que es temperatura mínima sean extremadamente bajos,
la distancia angular entre la línea ecuatorial a un como debería corresponder a zonas de gran altura,
punto de la tierra, al norte o al sur del ecuador, como es la zona evaluada en el presente estudio.
entre los –90º≤ φ ≤90º [29]; la longitud es el ángulo
formado por el meridiano de Greenwich y un punto Factores astronómicos
de la superficie de la tierra, el ángulo es negativo Para comprender el comportamiento de la radiación
en el oeste y positivo hacia el este [30]; y la altitud, solar extraterrestre (Re) se debe considerar que el sol
que es la cantidad de metros sobre el nivel del mar es una estrella variable debido a fenómenos como las
del lugar evaluado, en la Tabla 2 se muestran tales manchas y fáculas solares, presentes principalmente
valores para la zona de estudio, localizada en el en las zonas de convección, superficie y atmósfera,
Altiplano Puneño. provocando que el sol emita radiación variable a lo
largo del año [35]. Se conoce que estas variaciones
El modelo de Bristow y Campbell ha sido utilizado son menores a ± 1,5%; sin embargo, la distancia
en diversos estudios, debido a la simplicidad de tierra sol varía un 6,7% a lo largo del año [36].
datos utilizados [31]. En estados Unidos se demostró
que su desempeño pudo explicar entre un 70% y Para determinar la radiación solar global diaria,
90% de la variación de la radiación solar de la zona utilizando un modelo determinado, es indispensable
evaluada, esto se dio de manera similar en China [3], calcular la radiación solar extraterrestre, como se
27], en Bolivia se obtuvo una validez al 89% [21] y muestra en la ecuación (4), donde se consideran
al sureste de Brasil se obtuvo un mejor desempeño factores astronómicos como el factor de corrección
de la estimación de radiación solar global [32, 33]. E0, la declinación solar δ, el ángulo horario ω,
y la constante solar ICS; valores que pueden ser
En el Perú se han desarrollado diversos trabajos encontrados con la ecuación (5), ecuación (6) y
sobre Energía Solar, entre ellos se destaca la ecuación (7) respectivamente. Recordando además
generación de un Atlas de Energía solar [34], que φ representa la latitud. Además, se ha utilizado
trabajo multisectorial que buscaba establecer las ICS = 1367 W.m –2 valor indicado en [29].
zonas con elevado potencial de energías renovables,
cuantificando su disponibilidad a fin de conocer 24 ( 60 )
su distribución temporal a lo largo del territorio Re = E 0 Ics (ωsenφsenδ + cos φ cosδsenω) (4)
π
peruano, utilizando además diferentes metodologías
para probar la efectividad de diferentes modelos para El factor de corrección (E0) está relacionado con
la costa, sierra y selva del Perú, donde, al igual que la variación de la distancia entre la tierra y el sol,
en el trabajo de Baigorria et al. [25] concluyen que representado por el valor r, que varía a lo largo del
el modelo Bristow y Campbell es el que se adecua año con respecto a su valor medio r0 = 1,496 x 108
mejor a las condiciones de Perú, esto debido a que km, distancia denominada Unidad Astronómica [37].
la topografía peruana es compleja. Aplicaciones tecnológicas y de ingenierías aplican la

72
 Quispe y Sotomayor: Determinación y análisis temporal de la radiación solar global en el Altiplano de Puno

ecuación (5), que se da en función del día del año dn Tabla 1. Escala de transmitancia atmosférica.
[29], que corresponde al número del día Juliano del año
Transmitancia
(1≤dn ≤365), que varía desde el día 1 para el primero Tipo de día
Atmosférica
de enero y hasta el día 365 para el 31 de diciembre.
Nublado 0 < At ≤ 0,2
⎛ r ⎞2 ⎡ 2πd n ⎤
E 0 = ⎜ 0 ⎟ = 1+ 0, 33cos ⎢ (5)
Parcialmente nublado 0,2 < At ≤ 0,6
⎝ r ⎠ ⎣ 365 ⎥⎦ Soleado 0,6 < At ≤ 0,75
Muy soleado 0,75 < At <1
La declinación solar (δ) corresponde a la posición Fuente: [40].
angular del sol al mediodía solar respecto al plano
del ecuador terrestre, positivo hacia el norte [38], y
se puede calcular mediante la ecuación (6), obtenida dispersión por aerosoles, CO2, ozono, vapor de agua,
por Perrín de Brichambaut indicada en [29]. etc.; causados por diferentes agentes atmosféricos
y contaminantes [39, 5]. La interpretación física de
⎡ 2πd n ⎤ la transmitancia atmosférica se refiere al grado de
δ = 0, 409.sen ⎢ −1, 39⎥ (6)
⎣ 365 ⎦ claridad del cielo, considerando bandas que clasifican
el tipo de día, según la escala de transmitancia
El ángulo horario (ω) representa el desplazamiento atmosférica que se muestra en la Tabla  1 [40].
angular del sol al este o al oeste del meridiano local Además, la distribución de probabilidad indicada,
debido a la rotación de la tierra [35]. La ecuación se considera esencial para la predicción de la
(7) corresponde al ángulo horario de salida del sol. energía de la radiación solar global [41]. Cuando
el día es soleado, la radiación solar se incrementa y
ω = cos−1 (−tgφtgδ) (7) en presencia de nubosidad disminuye [42]; así, las
nubes que cubren de 40% a 60% de la superficie
Transmitancia atmosférica de la tierra, desempeñan un papel fundamental en
La transmitancia de la atmósfera At representa la el balance de energía del planeta, así como en la
fracción de energía solar que incide en la superficie variación de la radiación solar global [43].
terrestre, con respecto a la radiación solar que
llega hasta el tope de la atmósfera [33]. Al llegar a RESULTADOS Y DISCUSIÓN
la atmósfera, esta no se comporta como un filtro
homogéneo, sino que ocurren múltiples fenómenos Los resultados presentados en esta sección fueron
como la difracción, reflexión y absorción, así como obtenidos utilizando los datos de las estaciones de

Tabla 2. Estaciones meteorológicas SENAMHI Zonal Puno.

Estación Provincia Distrito Latitud Longitud Altitud m.s.n.m. Piso Ecológico
Ananea San Antonio de Putina Ananea –14º40’43’’ –69º32’04’’ 4660 Clima de altura
Capazo El Collao Ilave –17º11’16’’ –69º44’08’’ 4530 Clima de altura
Mazocruz Chucuito Mazocruz –16º44’25’’ –69º42’22’’ 4100 Clima de altura
Crucero Carabaya Crucero –14º21’44’’ –70º01’28’’ 4100 Puna húmeda
Huancané Huancané Huancané –15º12’05’’ –69º45’13’’ 3890 Puna húmeda
Ilave El Collao Ilave –16º05’18’’ –69º38’42’’ 3880 Puna húmeda
Putina San Antonio de Putina Putina –14º54’53’’ –69º52’04’’ 3878 Puna húmeda
Azángaro Azángaro Azángaro –14º54’52’’ –70º11’27’’ 3863 Puna húmeda
Desaguadero Chucuito Desaguadero –16º34’07’’ –69º02’25’’ 3860 Puna húmeda
Capachica Puno Capachica –15º36’56’’ –69º50’38’’ 3933 Circunlacustre
Moho Moho Moho –15º23’18’’ –69º29’04’’ 3890 Circunlacustre
Yunguyo Yunguyo Yunguyo –16º18’28’’ –69º04’29’’ 3890 Circunlacustre
Puno Puno Puno –15º49’35’’ –70º00’44’’ 3825 Circunlacustre
Juli Chucuito Juli –16º12’14’’ –69º27’36’’ 3812 Circunlacustre

73
Ingeniare. Revista chilena de ingeniería, vol. 30 Nº 1, 2022

Tabla 3. Radiación solar global diaria durante el mes de diciembre del 2017 de la estación Puno.
Radiación Radiación
Día Factor de Ángulo Radiación
Coeficiente Tmax Tmin Coeficiente Coeficiente solar global solar global
Fecha juliano corrección horario extraterrestre
A [–] [ºC] [ºC] C [–] B [–] estimada Rg medida Rg
dn [–] E0 [–] ω [rad.] Re (kWh.m–2)
(kWh.m–2 ) (kWh.m–2 )
01-dic 0,797 335 18,8 6,1 1,029 1,686 11,369 1,202 0,066 6,815 6,521
02-dic 0,819 336 18,9 5,5 1,029 1,687 11,374 1,151 0,074 7,159 7,369
03-dic 0,814 337 18,9 5,5 1,029 1,688 11,378 1,151 0,074 7,112 7,048
04-dic 0,812 338 17,9 6,3 1,029 1,688 11,382 1,281 0,056 6,679 6,328
05-dic 0,825 339 17,9 5,9 1,030 1,689 11,386 1,252 0,059 6,895 6,423
06-dic 0,824 340 17,6 6,4 1,030 1,690 11,390 1,310 0,052 6,665 6,377
07-dic 0,933 341 17,7 6,4 1,030 1,690 11,394 1,302 0,053 7,583 7,501
08-dic 0,745 342 17,7 5,1 1,030 1,691 11,397 1,209 0,065 6,365 6,317
09-dic 0,754 343 17,4 5,1 1,031 1,691 11,400 1,230 0,062 6,378 6,498
10-dic 0,844 344 19 3,2 1,031 1,692 11,403 0,978 0,113 7,844 7,295
11-dic 0,788 345 19 3,2 1,031 1,692 11,406 0,978 0,113 7,325 7,901
12-dic 0,815 346 18,4 6,3 1,031 1,693 11,409 1,245 0,060 6,853 6,610
13-dic 0,915 347 19 6,5 1,031 1,693 11,411 1,216 0,064 7,798 7,259
14-dic 0,947 348 19 6,5 1,032 1,693 11,413 1,216 0,064 8,078 8,926
15-dic 0,820 349 18,4 7,9 1,032 1,693 11,415 1,360 0,047 6,427 6,483
16-dic 0,923 350 17,8 6,8 1,032 1,694 11,417 1,324 0,051 7,421 7,664
17-dic 0,787 351 17,8 6,8 1,032 1,694 11,419 1,324 0,051 6,324 6,098
18-dic 0,693 352 16,4 5,7 1,032 1,694 11,421 1,346 0,049 5,493 5,938
19-dic 0,731 353 16,4 5,7 1,032 1,694 11,422 1,346 0,049 5,790 5,396
20-dic 0,660 354 17,6 6,9 1,032 1,694 11,424 1,346 0,049 5,233 5,697
21-dic 0,642 355 17,6 6,7 1,033 1,694 11,425 1,331 0,050 5,136 5,315
22-dic 0,660 356 17,7 6,6 1,033 1,694 11,426 1,317 0,052 5,330 5,596
23-dic 0,809 357 19,1 6,6 1,033 1,694 11,426 1,216 0,064 6,907 6,437
24-dic 0,759 358 19,1 6,6 1,033 1,694 11,427 1,216 0,064 6,480 6,307
25-dic 0,795 359 19 6,1 1,033 1,694 11,428 1,187 0,068 6,879 6,865
26-dic 0,720 360 17,7 5,7 1,033 1,693 11,428 1,252 0,059 6,037 5,801
27-dic 0,772 361 17,7 5,7 1,033 1,693 11,428 1,252 0,059 6,473 6,845
28-dic 0,872 362 16,3 6,3 1,033 1,693 11,428 1,396 0,044 6,650 6,287
29-dic 0,775 363 16,3 5,4 1,033 1,693 11,428 1,331 0,050 6,205 6,414
30-dic 0,849 364 12,4 4,7 1,033 1,692 11,427 1,562 0,033 5,328 5,411
31-dic 0,867 365 17,8 5,1 1,033 1,692 11,427 1,202 0,066 7,451 7,346

SENAMHI Zonal Puno listadas en la Tabla 2, las diaria (At), a fin de contrastar su comportamiento
cuales están dentro del área de estudio seleccionado con el propuesto por otros trabajos de investigación
en el presente trabajo. Tales datos permitieron analizar [25, 34]. Finalmente, se verifica si existen resultados
el comportamiento de la radiación solar global diaria, favorables para el aprovechamiento de la energía
mensual y anual a lo largo de los años seleccionados solar en el altiplano de Puno.
de 2007 a 2013, involucrando además datos que
pudieron ser recuperados de la Estación de Puno, Como ejemplo de todas las mediciones efectuadas y
correspondientes a los años 2014, 2015 y 2017. Los los datos estimados con el modelo propuesto, en la
datos estimados, obtenidos por el modelo de Bristow y Tabla 3 se muestran los valores correspondientes a la
Campbell, fueron contrastados con los datos medidos estación Puno, con latitud φ = 15º49’35’’, durante el
por las Estaciones, a fin de encontrar el nivel de mes de diciembre de 2017, donde se observa que la
correlación existente. Así mismo, se ha evaluado radiación solar global, durante todos los días supera al
el valor máximo de la transmitancia atmosférica valor de 4,000 kWh.m–2.dia–1, indicado por Camayo

74
 Quispe y Sotomayor: Determinación y análisis temporal de la radiación solar global en el Altiplano de Puno

et al. [44] como el mínimo valor aceptable para el de setiembre, octubre, noviembre y diciembre; sin
aprovechamiento de la energía solar. Así mismo, embargo, durante los meses de mayo, junio, julio y
viendo la columna de A, se observa que los valores agosto, la radiación solar global disminuye, debido
diarios superan 0,6, indicando que durante todo al solsticio de invierno en el hemisferio sur.
ese mes, el tipo de día fue soleado o muy soleado,
como se indica en la Tabla 1. Además, se observa que la radiación solar global
se comporta de forma estacional por ejemplo, los
La Figura 1 muestra el resumen de la variación entre valores más altos fueron alcanzados durante el mes
la Radiación solar global medida por el SENAMHI de noviembre, en la estación de Capazo, donde se
en la estación Puno (círculos rojos) y la estimada obtuvieron los valores máximos de 9,488; 9,290;
(cuadrados negros) usando el modelo Bristow y 9,290; 9,168; 9,174 y 9,092 kWh.m –2 .dia –1 en el
Campbell correspondiente, para los años 2014, 2015 periodo 2008 a 2013, respectivamente; y los valores
y 2017. Se observa la similitud entre ambos datos, más bajos fueron alcanzados durante el mes de
donde las oscilaciones estacionales de radiación solar junio, en la estación de Ananea, donde se obtuvieron
global se deben a los solsticios en el hemisferio sur valores de 4,700; 4,858 y 4,523 kWh.m –2.dia –1 en
del planeta, obteniéndose la correlación de Pearson los años 2007, 2008 y 2013 respectivamente.
R = 0,938, como se muestra en la Figura 2, que
indica un alto valor de confianza, al 93,8%. La Figura 4 muestra el valor máximo de transmitancia
atmosférica promedio, de los años 2015, 2016 y
En la Figura 3 se muestran los valores del promedio 2017, para ello se utilizaron datos medidos de la
mensual de la radiación solar global, estimadas con el estación meteorológica Davis del Laboratorio de
modelo Bristow y Campbell, correspondientes a los Física de la Universidad Nacional del Altiplano
años 2007 a 2013, de todas las estaciones involucradas Puno. Se observa que existe una mayor dispersión
en el presente estudio. Como era esperado, se ve en los extremos, durante los meses enero, febrero,
que, para los meses de enero, febrero, marzo y abril, marzo y abril, así como en los meses de setiembre,
los valores son altos; al igual que para los meses octubre, noviembre y diciembre, esto debido a la

Figura 1. Radiación solar global diaria estimada y medida en la estación

de Puno durante los años 2014, 2015 y 2017.

75
Ingeniare. Revista chilena de ingeniería, vol. 30 Nº 1, 2022

Figura 2. Correlación de Pearson de la radiación solar global

estimada y medida de la estación Puno durante los
años 2014, 2015 y 2017.

Figura 3. Promedio mensual de la Radiación solar global estimada en el altiplano puneño de

2007 a 2013.

presencia de días parcialmente nublados o nublados, Así mismo, la Figura 4 muestra que los valores de
algunas veces acompañados por precipitaciones, la transmitancia atmosférica máxima diaria (A) son
esto en contraste con las zonas menos dispersas, mayores a 0,6, lo cual indica que la mayoría de días
durante los meses de mayo, junio, julio y agosto, durante el año son soleados o muy soleados [40].
debido a la temporada de sequía y heladas, donde Además, se obtiene el promedio anual A = 0,84, que
los días son soleados o muy soleados. difiere en dos décimas el promedio anual sugerido

76
 Quispe y Sotomayor: Determinación y análisis temporal de la radiación solar global en el Altiplano de Puno

Figura 4. Promedio de la transmitancia atmosférica máxima diaria de la estación

Davis Pro Plus 2, años 2015, 2016 y 2017.

para Puno (A = 0,82) el 2003 y 2004 por [25] y [34], La Figura 5 muestra que el 2012 se obtuvo la radiación
estudios que utilizaron el modelo Angstrom-Prescott, solar global más baja, aunque con valores superiores
basados en insolación solar. Así, al compararlo con a 5,000 kWh.m –2.dia –1, lo que representa valores
los resultados de esta investigación, se encontró un alentadores para la aplicación de energías limpias
error relativo porcentual de 2,38%. Por esta razón, se y verdes. Además, en los resultados de radiación
considera el resultado confiable para realizar futuras solar global por zonas de estudio se observa en la
investigaciones, relacionadas con la determinación zona de Clima de Altura dos de las tres estaciones
de la radiación solar global en zonas altoandinas obtuvieron resultados superiores a 7,000 kWh.m –2.
del Perú, tales como: Puno, Junín, Cajamarca, dia –1, Capazo y Mazocruz, que obtuvieron
Ayacucho, Cusco, Cerro de Pasco, Huancayo, promedios anuales de 7,060 y 7,298 kWh.m –2.dia –1
Arequipa, Huancavelica, Huánuco, entre otros. respectivamente, a excepción de Ananea, que obtuvo
5,945 kWh.m –2.dia –1; en la zona de Puna húmeda,
La Figura 5 muestra el valor promedio anual de donde se encuentran las estaciones de Azángaro,
la radiación solar global estimada de 2007 a 2013 Crucero, Desaguadero, Huancané, Ilave y Putina
de las zonas de estudio, donde el valor más alto se obtuvieron los promedios anuales de 6,551;
corresponde a la Zona de Clima de Altura, estación 6,744; 5,908; 6,340; 6,228 y 6,783 kWh.m –2.dia –1
de Mazocruz, a 4100 m.s.n.m., con un promedio respectivamente; y la zona Circunlacustre, donde se
de 7,298 kWh.m –2.dia –1 entre los años analizados; encuentran las estaciones de Capachica, Juli, Moho,
así mismo, el valor mínimo corresponde a la Zona Puno y Yunguyo se obtuvieron promedios anuales
Circunlacustre, estación de Juli, a 3812 m.s.n.m., de 5,983; 5,481; 6,071; 5,997 y 6,342 kWh.m –2.dia–1
con un promedio de 5,481  kWh.m –2 .dia –1 entre respectivamente.
los años analizados. Tales resultados superan las
expectativas mínimas esperadas para el uso del Según Villicaña et al. [45], para el desarrollo de
recurso solar, pues se sabe por la Física Atmosférica instalaciones solares fotovoltaicas de tipo aislado
que, a medida que aumente la altura, la atmósfera se se requiere un nivel de radiación solar global entre
hace más tenue, por lo que el valor de la radiación 3,000 y 4,000 kWh.m –2.dia –1; para instalaciones
solar debe ser mayor en sitios de mayor altitud [35], solares con conexión a red se requiere un nivel de
esto se cumple en el altiplano de Puno. radiación solar de 4,000 kWh.m–2.dia–1, en cuanto a

77
Ingeniare. Revista chilena de ingeniería, vol. 30 Nº 1, 2022

Figura 5. Promedio anual de la radiación solar global estimada de la zona

de estudio, periodo de 2007 a 2013.

la producción de agua caliente sanitaria se requiere del altiplano puneño y el desarrollo de instalaciones
menores niveles de radiación solar diaria, entre 2,500 solares del tipo eléctrico y térmico.
y 3,000 kWh.m –2.dia –1. Por los valores obtenidos
en este estudio, se puede afirmar que los resultados En la literatura, otros estudios han considerado
son importantes para el uso del recurso solar a lo también el modelo Bristow y Campbell, obteniendo
largo del año, principalmente en zonas aisladas resultados particulares, como el estudio en Cúcuta,
de las redes del Servicio Eléctrico Interconectado Colombia, que posee disponibilidad de energía solar
Nacional (SEIN), más aún en zonas de mayor altitud constante durante todo el año, con un promedio
y de clima continental con gran amplitud térmica; anual de 5,335 kWh.m–2.dia–1 [47] y Nepal, donde la
sin embargo, cabe mencionar que costos, como el radiación solar global estimada y medida fue 3,210
de los paneles fotovoltaicos, no fueron considerados y 3,440 kWh.m –2.dia –1 respectivamente, resultado
en el presente estudio. poco alentador para la zona de estudio, ubicada
a 1350 m.s.n.m. [48], altitud muy por debajo del
Camayo et al. [28] realizaron un estudio aplicando área de estudio del presente trabajo, cuya altitud se
el modelo Bristow y Campbell para estimación de encuentra entre los 3812 y 4660 m.s.n.m.
la radiación solar global, obteniendo un promedio
anual de 6,000  kWh.m –2 .dia –1 en la sierra y CONCLUSIONES
4,000 kWh.m–2.dia–1 en la selva de la región de Junín,
donde concluyen que existe potencial en el recurso El modelo empírico Bristow y Campbell permite
solar de dicha región. Así mismo, Vanegas et al. determinar la radiación solar global diaria de
[46] realizaron la cuantificación y caracterización forma confiable en el altiplano de Puno en las
de la radiación solar en La Guajirá, Colombia, zonas Circunlacustre con 5,975 kWh.m –2 .dia –1;
encontrando valores en el intervalo de 6,000 y Puna Húmeda con 6,429 kWh.m –2.dia –1 y Clima
6,800 kWh.m –2.dia –1, indicando que los resultados de Altura con 6,768 kWh.m –2.dia –1 con datos de
son adecuados para el desarrollo de instalaciones entrada de temperaturas máximas, temperaturas
solares del tipo eléctrico y térmico. En el presente mínimas y latitud; así mismo, se ha obtenido el
trabajo, se ha obtenido los valores promedios anuales coeficiente de correlación R = 0,938 en relación a la
por piso ecológico, con 5,975 kWh.m –2.dia –1 en la ciudad de Puno tal como se muestra en la Figura 2.
zona Circunlacustre, 6,429 kWh.m –2 .dia –1 en la También, se demuestra que existe una influencia
zona Puna Húmeda y 6,768 kWh.m –2.dia –1 en la directa de los elementos climáticos de amplitud
zona Clima de Altura, tales resultados superan las térmica y latitud de la zona de estudio en relación
expectativas en la potencialidad del recurso solar a la radiación solar global.

78
 Quispe y Sotomayor: Determinación y análisis temporal de la radiación solar global en el Altiplano de Puno

Se han realizado estimaciones de radiación solar solar, según la Organización Latinoamericano de

global diaria en el Altiplano de Puno durante los años Energía Solar OLADE; como son el uso de bombas
2007 al 2013 con datos de entrada de SENAMHI hidráulicas solares, sistemas de transferencia de calor,
Zonal Puno obteniéndose un máximo valor promedio colectores solares térmicos, módulos fotovoltaicos,
de los años estudiados de 7,298 kWh.m–2.dia–1 en viviendas bioclimatizadas, entre otros. Además, se
la estación de Mazocruz ubicada a una altitud de espera que se pueda transformar la energía solar
4100 m.s.n.m. en la zona Clima de Altura y un valor en energía térmica y eléctrica, lo que ayudaría a
mínimo promedio de 5,481 kWh.m –2.dia –1 en la resolver problemas como heladas que son afrontadas
estación de Juli ubicada a una altitud de 3812 m.s.n.m. por pobladores, animales y plantas en las zonas
en la zona Circunlacustre, resultados propicios para altoandinas de la Región Puno.
la instalación de paneles fotovoltaicos.
AGRADECIMIENTOS
Finalmente, los resultados obtenidos de los niveles
de radiación solar global en las zonas del altiplano Al Servicio Nacional de Meteorología e Hidrología
de Puno, son mayores a 4,000 kWh.m –2 .dia –1 y (SENAMHI) Zonal Puno por los datos cedidos para
son muy rentables para las aplicaciones en energía la realización de esta investigación.

REFRENCIAS

[1] D.L. Hartmann. “Global Physical Climatology”. Radiation”. Appl. Energy. Vol. 179, pp. 437-
Elsevier. Second edition. 2016. 450. 2016.
[2] A.M. Manzano, L. Martín, F. Valero and C. [9] S. Steinmetz, A. Nunes Deibler e J. Baptista
Armenta. “A single method to estimate the da Silva. “Estimativa da produtividade de
daily global solar radiation from monthly arroz irrigado em função da radiação solar
data”. Atmos. Res. Vol. 166, pp. 70-82. 2015. global e da temperatura mínima do ar”.
[3] T. Pan, S. Wu, E. Dai and Y. Liu. “Estimating Ciência Rural. Vol. 43 Nº 2, pp. 206-211.
the daily global solar radiation spatial 2013.
distribution from diurnal temperature ranges [10] B. Jamil and E. Bellos. “Development of
over the tibetan plateau in China”. Appl. empirical models for estimation of global
Energy. Vol. 107, pp. 384-393. 2013. solar radiation exergy in India”. J. Clean.
[4] V.J. da Silva, C.R. da Silva, R.R. Finzi and N. Prod. Vol. 207, pp. 1-16. 2019.
da S. Dias. “Métodos para estimar radiação [11] Y. Feng, N. Cui, Y. Chen, D. Gong and X.
solar na região noroeste de Minas G erais”. Hu. “Development of data-driven models
Ciência Rural. Vol. 42 Nº 2, pp. 276-282. for prediction of daily global horizontal
2012. irradiance in Northwest China”. J. Clean.
[5] Q. Fu. “Radiación transfer in the Atmosphere”. Prod. Vol. 223, pp. 136-146. 2019.
Encyclopedia of Atmospheric Sciences. [12] J. Fan, L. Wu, F. Zhang, H.Cai, W. Zeng, X.
Vol. 5, pp. 37-44. 2015. Wang and H. Zou.“Empirical and machine
[6] G.C. Beruski, A.B. Pereira and P.C. Sentelhas. learning models for predicting daily global
“Desempenho de diferentes modelos de solar radiation from sunshine duration: A
estimativa da radiação solar global em ponta review and case study in China”. Renew.
grossa, PR”. Rev. Bras. Meteorol. Vol. 30 Sustain. Energy Rev. Vol. 100, pp. 186-212.
Nº 2, pp. 205-213, 2015. 2019. DOI: 10.1016/j.rser.2018.10.018.
[7] S.G. Gouda, Z. Hussein, S. Luo and Q. [13] J. Almorox, M. Benito and C. Hontoria.
Yuan. “Model selection for accurate daily “Estimation of global solar radiation in
global solar radiation prediction in China”. Venezuela”. Interciencia. Vol.  33 Nº  4,
J. Clean. Prod. 2019. pp. 280-283. 2008.
[8] G.E. Hassan, M.E. Youssef, Z.E. Mohamed, [14] D. Ocampo y R. Rivas. “Estimación de la
M.A. Ali and A.A. Hanafy. “New Temperature- radiación neta diaria a partir de Modelos de
based Models for Predicting Global Solar Regresión Lineal Múltiple”. Rev. Chapingo,

79
Ingeniare. Revista chilena de ingeniería, vol. 30 Nº 1, 2022

Ser. Ciencias For. y del Ambient. Vol.  19 transmissivity: Distribution and empirical
Nº 2, pp. 263-271. 2013. estimation around the central Andes”. Int. J.
[15] M. Bengulescu, P. Blanc and L. Wald. Climatol. Vol. 24 Nº 9, pp. 1121-1136. 2004.
“Characterizing Temporal Variability in [26] F. Meza and E. Varas. “Estimation of mean
Measurements of Surface Solar Radiation monthly solar global radiation as a function of
and its Dependence on Climate”. Energy temperature”. Agric. For. Meteorol. Vol. 100
Procedia. Vol. 97, pp. 164-171. 2016. Nº 2-3, pp. 231-241. 2000.
[16] J. Rojas Serrano, A. Vásquez Vega, F. García [27] K.L. Bristow and G.S. Campbell. “On the
Acevedo, D. Parra Peñaranda and E. Castro relationship between incoming solar radiation
Becerra. “Estimating missing data in historic and daily maximum an minimun temperatre”.
series of global radiation through neural Agric. For. Meteorol. Vol. 13 Nº 3, pp. 576. 1984.
network algorithms”. Sist. Telemática. Vol. 14 [28] B.F. Camayo Lapa, D.E. Condezo Hurtado,
Nº 37, pp. 9-22. 2016. A.Y. Ramos Cadillo, J.R. Massipe Hernández
[17] R. Torrez, A. Burgoa and E. Ricaldi. “Modelos and A.B. Camayo Vivas. “Estimation of global
de Estimación de la Radiación Solar para el solar radiation, using extreme temperatures,
Altiplano Central de Bolivia”. Rev. Boliv. applying the bristow-campbell model in
Física. Vol. 23, pp. 1-7. 2013. the junín region”. Ingeniare. Vol. 27 Nº 4,
[18] A.Dumas, A. Andrisani, M. Bonnici, G. Graditi, pp. 643-651. 2019.
G. Leanza, M. Madonia and M. Trancossi. “A [29] J.A. Duffie and W.A. Beckman. Solar
new correlation between global solar energy Engineering of Thermal Processes.Fourth
radiation and daily temperature variations”. edition. Canada. 2013.
Sol. Energy. Vol. 116, pp. 117-124. 2015. [30] A. Ben Othman, K. Belkilani and M. Besbes.
[19] A. Marzo, M. Trigo-Gonzalez, J. Alonso- “Global solar radiation on tilted surfaces in
Montesinos, M. Martínez-Durbán, G. López, Tunisia: Measurement, estimation and gained
P. Ferrada, E. Fuentealba, M. Cortés and energy assessments”. Energy Reports. Vol. 4,
F.J. Batlles. “Daily global solar radiation pp. 101-109. 2018.
estimation in desert areas using daily extreme [31] J. Almorox, M. Bocco and E. Willington.
temperatures and extraterrestrial radiation”. “Estimation of daily global solar radiation
Renew. Energy. Vol. 113, pp. 303-311. 2017. from measured temperatures at Cañada
[20] A.P. De Oliveira, A.J. Machado, J.F. Escobedo de Luque , Córdoba , Argentina”. Renew.
and J. Soares. “Diurnal eVolution of solar Energy. Vol. 60, pp. 382-387. 2013.
radiation at the surface in the city of São Paulo: [32] M. Conceição. “Evapotranspiração de
Seasonal variation and modeling”. Theor. Appl. referência com base na radiação solar global
Climatol. Vol. 71, Nº 3-4, pp. 231-249. 2003. estimada pelo método de Bristow-Campbell”.
[21] J. Álvarez, H. Mitasova and H.L. Allen. Eng. Agric. Jaboticabal. Vol. 30 Nº 4, pp. 619-
“Estimating Monthly Solar Radiation in 626. 2010.
South Central Chile”. Chil. J. Agric. Res. [33] A. Pacheco de Souza, A. Carvalho da Silva,
Vol. 71 Nº 4, pp. 601-609. 2012. A. Aki Tanaka, E. Morgan Uliana, F. Terra
[22] H. Oesterle. “Reconstruction of daily global de Almeida, A. Evaldo Klar e A. Almeida
radiation for past years for use in agricultural Gomes. “Global radiation by simplified
models”. Phys. Chem. Earth. Vol. 26, Nº 3, models for the state of Mato Grosso, Brazil”.
pp. 253-256. 2001. Pesqui. Agropecu. Bras. Vol. 52 Nº 4, pp. 215-
[23] C. Dejollx y A. litis. “El lago Titicaca: Sintesis 227. 2017.
del conocimiento limnologico actual”. 1991. [34] Servicio Nacional de Meteorología e
[24] Y. Feng, D. Gong, Q. Zhang, S. Jiang, L. Hidrología (SENAMHI). “Atlas de Energia
Zhao and N. Cui. “Evaluation of temperature- Solar del Perú”, pp. 72. 2003.
based machine learning and empirical models [35] C. Recalde, C. Cisneros, D. Vaca y C. Ramos.
for predicting daily global solar radiation”. “Relación de la transmitancia atmosférica con
Energy Convers. Manag. Vol. 198. 2019. la heliofanía y la diferencia de temperaturas
[25] G.A. Baigorria, E.B. Villegas, I. Trebejo, extremas diarias en la zona ecuatorial Andina”.
J.F. Carlos and R. Quiroz. “Atmospheric Inf. Tecnol. Vol. 26 Nº 1, pp. 143-150. 2015.

80
 Quispe y Sotomayor: Determinación y análisis temporal de la radiación solar global en el Altiplano de Puno

[36] L.J. Gray, J. Beer, M. Geller, J.D. Haigh, [43] E. Echer, M.P. Souza and N. J. Schuch.
M. Lockwood, K. Matthes, U. Cubasch, “A Lei de Beer Aplicada na Atmosfera
D. Fleitmann, G. Harrison, L. Hood, J. Terrestre”. Rev. Bras. Ensino Física. Vol. 23
Luterbacher, G.A. Meehl, D. Shindell, B. Geel Nº 3, pp. 276-283. 2001.
and W. White “Solar influences on climate”. [44] F.B. Camayo Lapa, R.J. Massipe Hernandez
Rev. Geophys. Vol. 50 Issue 1. 2012. E.J. Pomachagua Paucar, A. Torres-Ten
[37] M. Iqbal. “An introduction to solar radiation”. and M.O. Quispe-Flores. “Desarrollo del
Academic Press. First edition. Vancouver, modelo Bristow-Campbell para estimar la
Canada. 1983. radiación solar global de la Región de Junin,
[38] Y. El Mghouchi, A. El Bouardi, Z. Choulli Perú”. Tecnol. Química. Vol. XXXV Nº 2,
and T. Ajzoul. “Models for obtaining pp. 271-288. 2015.
the daily direct, diffuse and global solar [45] E. Villicaña-Ortiz, A.J. Gutiérrez-Trashorras,
radiations”. Renew. Sustain. Energy Rev. J.P. Paredes-Sánchez and J. Xiberta-Bernat.
Vol. 56, pp. 87-99. 2016. “Solar energy potential in the coastal zone of
[39] N. Falcón, F. Peña, H. Mavo y R. Muñoz. the gulf of Mexico”. Renew. Energy. Vol. 81,
“Irradiación solar global en la ciudad de pp. 534-542. 2015.
Valencia”. Rev. Ing. UC. Vol. 8 Nº 2, pp. 0. [46] M. Vanegas Chamorro, E. Villicaña Ortíz
2001. y L. Arrieta Viana. “Cuantificación y
[40] C. Yousif, G. Oña Quecedo and J. Bilbao Santos. caracterización de la radiación solar en
“Comparison of solar radiation in Marsaxlokk, el departamento de La Guajira-Colombia
Malta and Valladolid, Spain”. Renew. Energy. mediante el cálculo de transmisibilidad
Vol. 49 Nº June, pp. 203-206. 2013. atmosférica”. Prospect. Vol.  13 Nº  2,
[41] T.R. Ayodele and A.S.O. Ogunjuyigbe. pp. 54-63. 2015.
“Prediction of monthly average global solar [47] W. Contreras, M.G. Galban y S.B.Sepúlveda.
radiation based on statistical distribution of “Análisis estadístico de la radiación solar
clearness index”. Energy. Vol. 90, pp. 1733- en la ciudad de Cúcuta”. Entre Cienc. e Ing.
1742. 2015. Vol. 12 Nº 23, pp. 16. 2018.
[42] F.J. Lopes de Lima. “Previsão de irradiação [48] K.N. Poudyal, B.K. Bhattarai, B. Sapkota,
solar no nordeste do Brasil empregando o B. Kjeldstad and P. Daponte. “Estimation
modelo WRF ajustado por Redes Neurais of the daily global solar radiation; Nepal
Artificiais (RNAs)”. Tesis. São José dos experience”. Meas. J. Int. Meas. Confed.
Campos, Brazil. 2015. Vol. 46 Nº 6, pp. 1807-1817. 2013.

81