Acessibilidade / Reportar erro

Determination of Diffused Irradiation from Horizontal Global Irradiation - Study for the City of Curitiba

Abstract

Abstract

Understanding the solar radiation values ​​in the surface is important for the development of solar energy projects, obtaining through radiometers installed in certain places or publications as in solar energy atlases. Typically, solar or weather stations do not have sensors to measure diffuse irradiation due to high investment and the need for constant maintenance. This paper presents methods for the determination of diffuse irradiation from surface measured data of global horizontal irradiation in the city of Curitiba, obtained through pyranometers installed at the Federal University of Technology – Paraná (UTFPR) headquarters and from Curitiba-A807 National Institute of Meteorology (INMET) station. Two models are presented for the estimation of the monthly mean daily diffuse irradiation, based on the studies of Liu and Jordan (1960) and Page (1961). Finally, these data were compared with those presented by the Solar Energy Atlas - Paraná, verifying the percentage differences presented. It was observed that the model elaborated by Page presented better results in the diffuse radiation estimates when compared to the values ​​presented by the Paraná Atlas, with an average variation of -1.39% and -1.55%, for data from INMET and UTFPR respectively.

Keywords:
Diffuse Irradiation; estimate; global irradiation; solar energy


INTRODUCTION

Solar energy is not a renewable source of energy but an inexhaustible source of energy [11 Pereira, E.; Martins, F.; Gonçalves, A.; Costa, R.; Lima F.; Rüther, R.; Abreu, S.; Tiepolo, G.; Pereira, S.; Souza, J. Brazilian Atlas of Solar Energy, 2nd ed.; National Space Research Institute: São José dos Campos, Brazil, 2017; pp. 09-68.]. The publication by INPE (National Institute for Space Research) of the Brazilian Solar Energy Atlas 2nd edition (2017) shows that Brazil has great solar potential, where even in the less sunny place it is viable to install photovoltaic solar energy systems.

In December 2017, the first Solar Energy Atlas - Paraná was published, whose project was the result of a partnership between UTFPR, INPE, Itaipu Binational and the Itaipu Technological Park (PTI). However, during the project, it was perceived a difficulty in validating the data estimated by the BRASIL-SR model with the surface’s collected data, since the state of Paraná does not have its own solarimetric grid, needing this information to improve it and develop new research in this sector, whether thermosolar or photovoltaic systems.

Each system has its way of operation: thermals depend exclusively on direct solar irradiance, whereas photovoltaics uses all available irradiance, either on a fixed inclined plane or through solar trackers [22 Souza, G. K.; Macagnan M. H. Comparison of hourly diffuse radiation models from data measured by SONDA station network in some cities in Brazil. VII CBENS 2018, 01-09.]. In Figure 1 is presented a basic representation of the solar radiation components, according to the Solar Energy Atlas - Paraná [33 Tiepolo, G; Pereira, E.; Urbanetz JR, J.; Pereira, S.; Gonçalves, A.; Lima, F.; Costa, R.; Alves, A. Solar Energy Atlas – Paraná, 1st ed.; National Space Research Institute: São José dos Campos, Brazil, 2017; pp. 10-84.].

Figure 1
Representation of Solar Radiation Components.

According to Pinho and Galdino [44 Pinho, J.; Galdino, M. Engineering Manual for Photovoltaic Systems. Rio de Janeiro: Cresesb Cepel; 2014. p. 66-101.] the radiation components are described as:

Extraterrestrial irradiance: Solar irradiance that reaches the top of the Earth's atmospheric layer.

Beam irradiation: Solar irradiation that focus directly on the surface without any influence.

Diffuse irradiation: Solar irradiation that reaches the surface after being scattered throughout the earth's atmosphere.

Reflected irradiation (albedo): irradiation reflected by the surrounding environment (soil, vegetation, obstacles, rocky terrain, etc.).

Global irradiation: amount of irradiation resulting from the sum of direct, diffuse and albedo solar irradiations.

For photovoltaic use, the most interesting is the global horizontal irradiation, which quantifies the radiation received by a horizontal flat surface. On cloudy days, the main portion is from the horizontal diffuse irradiation, whereas on clear days direct irradiation prevails [55 Energy Research Company. Technical Note: Analysis of Solar Generation Insertion in the Brazilian Electrical Matrix. Available online: http://www.cogen.com.br/content/upload/1/documentos/Solar/Solar_COGEN/NT_EnergiaSolar_2012.pdf (accessed on 14 July 2018).
http://www.cogen.com.br/content/upload/1...
].

Currently at INMET’s climatological stations, only the horizontal global component is collected. The accuracy and quality of solar radiation data collected by equipment is very important and depends on the manufacturer's recommended installation, operation, and maintenance, and should follow the International Organization for Standardization (ISO) and World Meteorological Organization (WMO) that establish ratings and specifications for solar irradiance measurement [11 Pereira, E.; Martins, F.; Gonçalves, A.; Costa, R.; Lima F.; Rüther, R.; Abreu, S.; Tiepolo, G.; Pereira, S.; Souza, J. Brazilian Atlas of Solar Energy, 2nd ed.; National Space Research Institute: São José dos Campos, Brazil, 2017; pp. 09-68.].

Due to the different characteristics of the components and the difficulty of obtaining the diffuse fraction out of the measured data, this work aims to determine the diffuse irradiation from global horizontal surface measured irradiation in the city of Curitiba, using the models studied by Liu and Jordan [66 Liu, B.Y.H.; Jordan, R.C. The interrelationship and characteristic distribution of direct, diffuse and total solar radiation. Solar Energy 1960, v.4, pp.01-19.] and Page [77 Page, J. K. The estimation of monthly mean values of daily total short-wave radiation on vertical and inclined surfaces from sunshine records for latitudes 40ºN- 40ºS. Proc. U.N. Conf. New Sources Energy 1961, Vol 4, pp. 378-390.], and compare the data estimated by these models with the data presented in the Solar Energy Atlas - Paraná.

MATERIAL AND METHODS

Global radiation is monitored by two KIPP & ZONEN - CMP3 Second Class thermoelectric pyranometers [88 KIPP&ZONEN. Kipp & Zonen. Available online: http://www.kippzonen.com/ (accessed on 14 July 2018).
http://www.kippzonen.com/...
], as shown in Figure 2 and positioned in a horizontal plane. The data obtained are from 12 months, from June 2017 to May 2018 and were acquired from INMET’s pyranometer installed at Curitiba-A807 Station [99 National Institute of Meteorology. Available online: http://www.inmet.gov.br (accessed on 16 July 2018).
http://www.inmet.gov.br...
], Curitiba - PR, with latitude of 25º26.922' S and longitude of 49º13.836' W. The other pyranometer is located at UTFPR, Curitiba central campus, located at 3165 Sete de Setembro Avenue, with latitude 25º26.352' S and longitude 49º16.176' W.

Figure 2
Pyranometer CMP3.

The data collected from the horizontal global component, both by the weather station and by the UTFPR, are in solar irradiance (W/m2), stored hourly, therefore, they are also hourly solar irradiation data (Wh/m2). To estimate diffuse irradiation, assumptions and calculations presented in this section are required.

Extraterrestrial solar radiation on a horizontal surface is obtained, according to Macagnan [1010 Macgnan, H.M. Introduction to solar radiation. Sinos River Valley University: São Leopoldo, Brazil, 2010; pp. 02-90.], by Equation 1:

Ion = Isc . Eo (1)

where, Ion is the extraterrestrial irradiance, Isc is the solar constant equal to 1367 W/m2 and Eo is the earth-orbit eccentricity correction factor, described in Equation 2:

Eo = Isc 1 + 0,033. cos 360. n 365 (2)

Where n is the day in the Julian calendar, ranging from 1 (January 1) to 365 (December 31), shown in Table 1.

Table 1
Standard average day of the month.

By obtaining the solar irradiation data on the terrestrial surface through the INMET and UTFPR pyranometers, and from the extraterrestrial irradiation calculation, it is obtained the values of the atmospheric transmissivity index (Kt), which is defined as the ratio between the solar radiation on the earth's surface (H) to the solar radiation reaching the top of the atmosphere (Ho). It varies according to the amount of clouds and aerosols in the atmosphere, and there may be increase or decrease of direct or diffuse components of solar radiation that reach the surface [1111 Querino, C. A.; Moura, M. A.; Silva Querino, J. K.; Filho, A. O. Study of Global Solar Radiation and External and Internal Transmissivity Index (Kt) in a mangrove forest in Alagoas - Brazil. RBMET 2011, v.26, pp.204 – 294.]. Kt is defined by Equation 3:

Kt = H Ho (3)

To calculate monthly mean daily horizontal extraterrestrial irradiation, the average day of the year is used, which is the day of the year when daily extraterrestrial radiation is almost equal to the monthly average value. It was noted that the 16th day of each month can lead to minor errors in Ho, especially for June and December, which are the solstices. The recommended days for each month are given in Table 1. It can be estimated Ho from Equation 4:

Ho = 24 π . Isc . Eo [ π 180 . ωs . sin δ . sin φ + cos δ . cos φ . sin ωs ] (4)

where, ωs is the angle of the sun's rise, in degrees, and can be obtained from Equation 5:

ωs = cos 1 tg φ . tg δ (5)

where φ is the local latitude, in degrees and δ is the solar declination, in degrees, and is calculated from Equation 6

δ = 23,45 ° . sin 360 284 + n 365 (6)

After determining the Kt index, it is estimated the variation of the daily diffuse fraction Kd, which is the ratio between diffuse irradiation Hd and global horizontal irradiation in the surface H, and is defined by the Equation 7:

Kd = Hd H (7)

The first method used to estimate the monthly mean daily diffuse solar irradiation on horizontal surface was Liu and Jordan’s [66 Liu, B.Y.H.; Jordan, R.C. The interrelationship and characteristic distribution of direct, diffuse and total solar radiation. Solar Energy 1960, v.4, pp.01-19.], is defined by Equation 8 and has as its parameter the Kt index.

Hd H = 1,39 4,027. Kt + 5,531. Kt ² 3,108. Kt ³ (8)

The second method used to estimate the monthly mean daily diffuse solar irradiation was Page’s [77 Page, J. K. The estimation of monthly mean values of daily total short-wave radiation on vertical and inclined surfaces from sunshine records for latitudes 40ºN- 40ºS. Proc. U.N. Conf. New Sources Energy 1961, Vol 4, pp. 378-390.] and is defined by Equation 9, also as a function of the Kt index:

Hd H = 1,00 1,13. Kt (9)

RESULTS

Through the data obtained by the INMET and UTFPR’s pyranometers, in the 12-month period (June 2017 to May 2018), it can be seen that these values, compared to the Solar Energy Atlas - Paraná, presented values close to the Atlas, which has estimates of solar radiation based on 17 years of satellite data (1999 to 2015).

Table 2 and Table 3 show the months of the year analyzed, the calculated values of transmissivity indices (Kt) and diffuse fraction (Kd) for INMET and UTFPR.

Table 2
Comparison of data with calculated Kd by the Liu and Jordan’s method.
Table 3
Comparison of data with calculated Kd by the Page’s method.

Measured horizontal global irradiation and estimated horizontal diffuse irradiation according to the Liu - Jordan and Page methods, respectively, were also presented. Finally, the percentage variation between the diffuse irradiation values estimated by the models and those presented by the Solar Energy Atlas - Paraná is presented.

The indices Kt and Kd are dimensionless values, since they represent the fraction of the components that arrive at surface, being the global and diffuse irradiations represented in W.h/m2.day.

The methods for determination of the diffuse irradiation component chosen from Liu and Jordan (Table 2) and Page (Table 3) presented expected diffuse irradiation data when compared to those presented by the Atlas, where it is assumed that the models are viable to estimate it.

DISCUSSION

It was found that the mean variation between diffuse irradiation obtained by the Liu and Jordan method for UTFPR in relation to Atlas is -16%, the same percentage found among the values obtained for INMET in relation to Atlas. This shows that the data obtained by the studied sites are very similar and have a considerable average variation in relation to the Atlas data, which was the reference used.

However, for the Page method, the average variation of diffuse irradiation presented between the values obtained for UTFPR in relation to Atlas was -1.55%, while for INMET the variation was -1.39%. These results show that the studied methods presented values very close to those found in the Atlas in the analyzed period.

Studies such as Crotti and Rampinelli [1313 Crotti, P.; Rampinelli, G. A. Estimation of the direct and diffuse solar components in horizontal surface for Araranguá/SC from meteorological networks. VII CBENS 2018, 01-09.] show that the lower the value of the Kt index, the higher the value of the Kd index, which can be observed in Figure 3 and Figure 4. This shows the convergence of the studies to different locations analyzed, proving the validation of the used methods.

Figure 3
Comparison between Kt and Kd (Liu e Jordan’s method).
Figure 4
Comparison between Kt and Kd (Page’s method).

CONCLUSION

This work presented analyzes for models in order to estimate the average daily diffuse irradiation by choosing Liu and Jordan’s and Page’s models. With the results presented in Tables 2 and 3, there was less variation between the data when calculated by the Page’s method. This does not mean that the Liu and Jordan’s method is less reliable, because it must be taken into account that the studied period is short compared to the values presented in the Solar Energy Atlas - Paraná in which a comprehensive historical data, between 1999 and 2015, was used.

It can also be seen from the comparisons that the data measured by the INMET-Curitiba Station and the UTFPR, of horizontal global irradiation are similar.

It was found that in the periods between April and August (autumn and winter seasons) the proportion of diffuse irradiation was smaller than in the other months of the year, characterized by less cloudiness and sunnier days compared to other periods of the year. This information is confirmed by the maps presented in the Solar Energy Atlas - Paraná, where seasonal daily averages of diffuse irradiation are shown.

Finally, to know which of these models presents the best estimate of diffuse irradiation, it would be ideal to obtain the measured values of this component, which is intended to be accomplished in a next step of this research.

Acknowledgments

The authors would like to thank UTFPR for the support and infrastructure available for the development of this research, and COPEL - Distribution for the support and funding of the resources for the realization of this R&D project “ANEEL RD 2866-0464 / 2017 - METHODOLOGY FOR ANALYSIS, MONITORING AND MANAGEMENT BY ENCOURAGED SOURCES”.

  • HIGHLIGHTS
  • • Review of solar radiation components.
  • • Acquisition of global irradiation data from INMET and UTFPR.
  • • Use of Liu and Jordan’s and Page’s method to calculate diffuse irradiation.
  • • Comparison with the results obtained in the Solar Energy Atlas - Paraná.

REFERENCES

  • 1
    Pereira, E.; Martins, F.; Gonçalves, A.; Costa, R.; Lima F.; Rüther, R.; Abreu, S.; Tiepolo, G.; Pereira, S.; Souza, J. Brazilian Atlas of Solar Energy, 2nd ed.; National Space Research Institute: São José dos Campos, Brazil, 2017; pp. 09-68.
  • 2
    Souza, G. K.; Macagnan M. H. Comparison of hourly diffuse radiation models from data measured by SONDA station network in some cities in Brazil. VII CBENS 2018, 01-09.
  • 3
    Tiepolo, G; Pereira, E.; Urbanetz JR, J.; Pereira, S.; Gonçalves, A.; Lima, F.; Costa, R.; Alves, A. Solar Energy Atlas – Paraná, 1st ed.; National Space Research Institute: São José dos Campos, Brazil, 2017; pp. 10-84.
  • 4
    Pinho, J.; Galdino, M. Engineering Manual for Photovoltaic Systems Rio de Janeiro: Cresesb Cepel; 2014. p. 66-101.
  • 5
    Energy Research Company. Technical Note: Analysis of Solar Generation Insertion in the Brazilian Electrical Matrix. Available online: http://www.cogen.com.br/content/upload/1/documentos/Solar/Solar_COGEN/NT_EnergiaSolar_2012.pdf (accessed on 14 July 2018).
    » http://www.cogen.com.br/content/upload/1/documentos/Solar/Solar_COGEN/NT_EnergiaSolar_2012.pdf
  • 6
    Liu, B.Y.H.; Jordan, R.C. The interrelationship and characteristic distribution of direct, diffuse and total solar radiation. Solar Energy 1960, v.4, pp.01-19.
  • 7
    Page, J. K. The estimation of monthly mean values of daily total short-wave radiation on vertical and inclined surfaces from sunshine records for latitudes 40ºN- 40ºS. Proc. U.N. Conf. New Sources Energy 1961, Vol 4, pp. 378-390.
  • 8
    KIPP&ZONEN. Kipp & Zonen. Available online: http://www.kippzonen.com/ (accessed on 14 July 2018).
    » http://www.kippzonen.com/
  • 9
    National Institute of Meteorology. Available online: http://www.inmet.gov.br (accessed on 16 July 2018).
    » http://www.inmet.gov.br
  • 10
    Macgnan, H.M. Introduction to solar radiation Sinos River Valley University: São Leopoldo, Brazil, 2010; pp. 02-90.
  • 11
    Querino, C. A.; Moura, M. A.; Silva Querino, J. K.; Filho, A. O. Study of Global Solar Radiation and External and Internal Transmissivity Index (Kt) in a mangrove forest in Alagoas - Brazil. RBMET 2011, v.26, pp.204 – 294.
  • 12
    Klein, S. A. Calculation of montlhly average insolation on tilted surfaces. Solar Energy 1977, Vol 19, pp. 325-329.
  • 13
    Crotti, P.; Rampinelli, G. A. Estimation of the direct and diffuse solar components in horizontal surface for Araranguá/SC from meteorological networks. VII CBENS 2018, 01-09.

Publication Dates

  • Publication in this collection
    25 Nov 2019
  • Date of issue
    2019

History

  • Received
    05 Nov 2018
  • Accepted
    26 July 2019
Instituto de Tecnologia do Paraná - Tecpar Rua Prof. Algacyr Munhoz Mader, 3775 - CIC, 81350-010 Curitiba PR Brazil, Tel.: +55 41 3316-3052/3054, Fax: +55 41 3346-2872 - Curitiba - PR - Brazil
E-mail: babt@tecpar.br