Next Article in Journal
An Optimized Sensorless Charge Balance Controller Based on a Damped Current Model for Flyback Converter Operating in DCM
Next Article in Special Issue
An LM-BP Neural Network Approach to Estimate Monthly-Mean Daily Global Solar Radiation Using MODIS Atmospheric Products
Previous Article in Journal
Comparison of Cassava Starch with Corn as a Feedstock for Bioethanol Production
Previous Article in Special Issue
Global Solar Radiation Prediction Using Hybrid Online Sequential Extreme Learning Machine Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multiple-Regression Method for Fast Estimation of Solar Irradiation and Photovoltaic Energy Potentials over Europe and Africa

1
Department of Control and Computer Engineering, Politecnico di Torino, 10129 Turin, Italy
2
Department of Energy, Politecnico di Torino, 10129 Turin, Italy
*
Author to whom correspondence should be addressed.
Energies 2018, 11(12), 3477; https://doi.org/10.3390/en11123477
Submission received: 15 November 2018 / Revised: 3 December 2018 / Accepted: 10 December 2018 / Published: 13 December 2018
(This article belongs to the Special Issue Solar and Wind Energy Forecasting)

Abstract

:
In recent years, various online tools and databases have been developed to assess the potential energy output of photovoltaic (PV) installations in different geographical areas. However, these tools generally provide a spatial resolution of a few kilometers and, for a systematic analysis at large scale, they require continuous querying of their online databases. In this article, we present a methodology for fast estimation of the yearly sum of global solar irradiation and PV energy yield over large-scale territories. The proposed method relies on a multiple-regression model including only well-known geodata, such as latitude, altitude above sea level and average ambient temperature. Therefore, it is particularly suitable for a fast, preliminary, offline estimation of solar PV output and to analyze possible investments in new installations. Application of the method to a random set of 80 geographical locations throughout Europe and Africa yields a mean absolute percent error of 4.4% for the estimate of solar irradiation (13.6% maximum percent error) and of 4.3% for the prediction of photovoltaic electricity production (14.8% maximum percent error for free-standing installations; 15.4% for building-integrated ones), which are consistent with the general accuracy provided by the reference tools for this application. Besides photovoltaic potentials, the proposed method could also find application in a wider range of installation assessments, such as in solar thermal energy or desalination plants.

1. Introduction

The world population is estimated to reach nearly 8.1 billion in 2040, with an average global economic growth of 3.4% per year. As a result, energy scenarios forecast a remarkable increase in global energy demand of about 30% by 2040 [1]. In line with the ratified agreements for worldwide sustainable development [2], a significant step up of renewables over conventional fossil fuels is expected in the global energy mix. The average net capacity addition per year of renewable energy is foreseen to grow steadily in the near future, led mainly by electricity production via solar photovoltaic (PV) conversion (see Figure 1a).
In recent years, the PV market in Europe has experienced a tremendous expansion, owing to the EU directive 2009/28/EC on the promotion of energy from renewable sources. Such initiative aims at fulfilling 20 % of the overall energy needs of Europe with renewables by 2020 [4], and 27 % by 2030 [5]. Considering that building-integrated installations have the advantage to preserve natural landscape with respect to large free-standing installations (PV farms), they have been promoted by local governments via specific policies, which led to a continuously increasing number of installations [6]. On the other hand, despite the outstanding potential in terms of solar energy and thus of sustainable electricity production [7], investments are still at the beginning in Africa [8]. In this sense, proper economic arguments on this potential would help to deploy investments, with beneficial effects particularly in rural areas, where no connection to the electric power grid is available. Photovoltaic systems could also favor the democratization of the energy availability [9], especially in Sub-Saharan Africa [10]. To promote and drive investments in Europe and Africa in this rapidly expanding market, one key aspect is the availability of fast and reliable screening tools for estimating: (i) the intensity of the available solar resource, which depends on the geographical location and climate; (ii) the potential energy output, which depends on the PV module technology, installation and local ambient conditions.
Methods for estimating the annual energy yield of PV systems can be classified into direct and indirect approaches [11]. The former evaluates the electrical energy output directly from the solar irradiation (i.e., insolation), whereas the latter obtains the energy output from the solar irradiation, ambient temperature, and some additional ancillary parameters [3,12]. The common starting point of these methods is a reliable measure or estimate of the available solar potential for specific geographical locations. In general, solar irradiation is best obtained via experimental measures by transducers (e.g., pyranometers and pyrheliometers in dedicated stations). However, this implies installation and monitoring difficulties in remote areas, especially those characterized by poorly developed technological access. Furthermore, such measurement stations are not suitable for collecting data in large-scale territories, because of their high capital and operating costs.
For this reason, a growing number of (online) tools and databases with solar irradiation data have been developed lately, for instance PVWatts [13], PVGIS [14], Global Atlas [15] or Solargis [16]. In particular, PVWatts is a web application implemented by the National Renewable Energy Laboratory (NREL) that can be used to predict the electricity generation from PV systems given the geographical position of the installation site and some technical details of the plant (e.g., size, module type, array type, losses, tilt, azimuth). This tool has been developed to be easily accessible and usable by both non-experts and advanced users. The prediction errors of the PVWatts model are claimed to be in the range ± 10 % for annual energy totals [13]. PVGIS, instead, is an online calculator of potential electricity production by PV systems developed by the Joint Research Center (JRC) [17]. In addition to the annual electricity production, PVGIS is also capable to provide estimations of the monthly and hourly ones. The required inputs are again the geographical location of the PV system and its technological characteristics. The database used by PVGIS for solar radiation in Europe and Africa (see e.g., the total in-plane solar irradiation in Figure 1b) demonstrated a mean bias error of about 2 % against satellite data [17].
However, such online tools sometimes provide low-resolution data (i.e., in the range of only a few kilometers), which may be scattered and require multiple queries to the online databases in case of large-scale analysis. Thus, various models have been developed for forecasting the local solar resource using only the most widespread climatic data [18]. Both empirical and non-empirical models have been implemented for estimating monthly, daily, and hourly global solar irradiation [19,20]. While short-term estimations methods deserve greater attention [21,22] for dealing with the conversion efficiency [3,23], long-term estimations are essential to analyze possible investment scenarios with respect to the PV solar energy potential [24,25]. These can be obtained using either data pattern analysis [26,27] or physical modeling approaches [28,29].
In this work, we propose a simple yet effective model to estimate the yearly solar irradiation per unit area, only taking into account the latitude, altitude, and average temperature of a certain location as input parameters. The yearly PV potential output is then obtained considering the different climatic and technological factors affecting electricity production, namely the temperature, reflection, module, and installation efficiencies, as well as the non-optimal azimuthal angle of the module. This procedure is first formulated and then validated against well-established databases for both free-standing and building-integrated systems. This approach may be particularly beneficial when only scattered data of solar irradiation are available over large territories but high-resolution analysis would be required. Furthermore, the presented methodology—which is now implemented for Europe and Africa—does not require querying online databases and thus has the capability to operate offline.
The layout of this work is as follows. The proposed models for the yearly sum of global irradiation and PV energy potential are introduced in Section 2. The parametrization of these models for Europe and Africa is reported in Section 3, together with a comparison of the estimated solar irradiation and electricity generation by PV systems with those of a well-established online database, using standard error metrics. In Section 4, the final conclusions are drawn and an outlook on the perspective developments and applications of the present work are given.

2. Methodology

The aim of the proposed methodology is to provide: (i) a fast means for evaluating the global energy irradiation in locations where no or scattered data are available; (ii) a fast and easy-to-use screening tool for the preliminary assessment of the perspective PV energy yield over large territories. As sketched in Figure 2, the approach relies on three basic and easily accessible data types, namely the latitude ( ϕ ), altitude above sea level (h), and average daily ambient temperature ( T 24 ) for a given site. We remark that, for the purpose of the analysis, the average daytime temperature would be a more suitable parameter as it is most closely related to sunlight; however, this information is generally less accessible. Therefore, we use T 24 values because they are more easily retrievable for a large number of sites. Clearly, this methodology takes into account local weather conditions only implicitly, since wind, clouds, humidity, fog, pollution, or other environmental conditions indirectly affect the global energy irradiation estimate through the average daytime temperature. Additionally, system-dependent factors such as installation type (free-standing or building-integrated), azimuth of the module respect to the optimal orientation ( α z ), solar exposed surface ( S inst ), PV module ( η mod ), and installation ( η inst ) efficiencies are also required for the final estimation of the yearly PV energy output. The multiple-regression method allowed the correlation of these inputs with the yearly solar irradiation ( H y ) and PV energy output ( Π y ) for the considered installation.
First, based on the three selected geodata types, the following general equation is found to provide an accurate estimation of the average yearly sum of global solar irradiation:
H y = w 1 | ϕ | + w 2 h + w 3 T 24 2 + w 4 | ϕ | T 24 2 + w 5 .
In the above equation, the latitude ϕ is considered in absolute value, meaning that the dependence of H y on this parameter is symmetrical with respect to the equator. The model also reports that the yearly sum of global solar irradiation depends linearly on the altitude of the selected location and quadratically on the average temperature. This evidence is physically meaningful, since insolation has been repeatedly found to rely on latitude [31], altitude [32], and climate conditions (e.g., average temperature) [33]; whereas, no effect of longitude ( λ ) is typically noticed. Finally, w i ( i = 1 , , 5 ) are the weights (or coefficients) of the polynomial model, which will be fitted in Section 3 over an extensive set of geographic locations throughout Europe and Africa (see Figure 3).
Second, to predict the specific electrical output by a PV system at a given location, one must consider the various climatic and technological losses ψ and the corresponding efficiencies η = 1 ψ . In particular, we consider: (i) the efficiency resulting from losses due to temperature and low-irradiance effects η temp ; (ii) the efficiency resulting from losses due to angular reflectance effects η refl ; (iii) the efficiency of the module technology itself η mod , and (iv) the additional efficiency resulting from losses due to system installation η inst (e.g., inverter and cables). The final correction factor applied to the available solar insolation is the product of the aforementioned efficiencies, namely
η tot = η temp η refl η mod η inst .
Among the various correlations available in the literature [34,35], here the temperature and low-irradiance efficiency is considered to have the following form [24]:
η temp = p 1 T 24 2 + p 2 T 24 + p 3 ,
where the model coefficients p i depend on the installation type (free-standing or building-integrated) and PV technology. Similarly, one of the possible correlations [36] to model the reflectance efficiency is a fourth-order polynomial using the azimuthal angle between the module orientation and either the South (northern hemisphere) or North (southern hemisphere) direction as independent variable [24]:
η refl = q 1 α z 4 + q 2 α z 3 + q 3 α z 2 + q 4 α z + q 5 .
Note that the model coefficients p i in Equation (3) and q i in Equation (4) will be obtained in Section 3 by fitting the data of the geographic locations reported in Figure 3. The remaining coefficients, namely η mod and η inst , depend on user choice for the PV technology and installation, respectively.
One further correction of the PV energy output must be considered when non-optimal azimuthal angles of the PV installation are imposed. In fact, H y refers to the yearly sum of global solar irradiation under optimal orientation of the PV module, namely α z = 0 . Usually, this mounting angle is not always possible (especially in building-integrated installations), thus a correction factor F ( α z ) must be included (note that the tilt angle of the modules may also be considered; however, in this work we assume optimal tilt). This scaling factor for the non-optimal azimuth can be obtained as a function of the azimuthal angle only as [24]
F ( α z ) = r 1 α z 4 + r 2 | α z 3 | + r 3 α z 2 + r 4 | α z | + r 5 ,
being r i the model coefficients. In the next Section, we will show that proper parametrizations of this model allow to obtain well-representative functions even for large geographical extents.
Finally, the electrical energy output of the PV system is computed as
Π y = η tot H y F S inst ,
being S inst the net exposed surface of the installed modules.
In summary, as schematically depicted in Figure 2, the proposed methodology for fast estimation of solar irradiation (Equation (1)) depends directly on geographical data (latitude, altitude) and indirectly on climatic data (average daytime temperature, which is implicitly influenced by wind, clouds, humidity, fog, pollution or other local environmental conditions). Instead, the estimation of photovoltaic energy potential (Equations (2) and (6)) relies on both climatic data (average daytime temperature, Equation (3)) and installation characteristics of the considered PV system (photovoltaic technology, Equation (2); azimuth, Equations (4) and (5); installation type—that is building-integrated or free-standing).

3. Results

In this Section, the proposed methodology is tuned, validated, and then adopted to provide possible scenarios of the yearly PV production over large territories in Europe and Africa. First, the yearly sum of global solar irradiation is obtained according to Equation (1), which requires as inputs the latitude ϕ , the altitude h and the mean temperature T 24 for a given geographical location. Then, the temperature and reflectance efficiencies are computed according to Equations (3) and (4) for optimal tilt and azimuthal angles of the modules. These efficiencies are then combined with the user-defined efficiencies η mod and η inst , to compute the total efficiency η tot according to Equation (2). The yearly electrical energy output is eventually estimated via Equation (6), based on the net surface of the installed PV modules. The analysis is finally carried out also for module installations that present non-optimal azimuthal angles, by considering the correction factor F ( α z ) .

3.1. Solar Irradiation Model

The model coefficients for Equation (1) are fitted to a random set of 80 different geographical locations in Europe and Africa, see Figure 3 and Table A1 and Table A2, with latitudes covering a range from about 60 to −30 . Note that the exact geographical coordinates have been chosen based on the availability of the required data, so they do not generally coincide with the center of the indicated cities. These cities are indeed indicated as the most significant nearby location to the considered geographical point. The altitude, average temperature in the 24 h and the yearly sum of global solar irradiation for each location are extracted from the online database PVGIS [37] of the Joint Research Centre [17,38,39]. The solar irradiation database of PVGIS is based on data records of the Satellite Application Facility on Climate Monitoring (CM-SAF) [40]. The average local temperature ( T 24 ) is provided by the interactive map, which covers Europe and part of the Northern Africa. All the temperature data for the locations in Africa with latitude lower than 32.5 are obtained from Berkeley Earth™ [41]. The H y values are extracted from the PVGIS database considering both optimal azimuthal angle and optimal tilt angle.
The full set of data for these 80 random locations is split in two different subsets, namely a training and a validation set. The model coefficients of Equation (1) are initially fitted on the training data set; then, the best-fitted model is used to predict the responses for the observations in the second data set (the validation one), which provides an unbiased evaluation of the model fitting. Here, the training data set is taken as 70% (i.e., 56 locations) of the whole data set; whereas, the validation one includes 24 locations. The training/validation process is iteratively repeated 10,000 times, with a random distribution of locations in the two subsets, and standard error metrics are computed with respect to the H y values provided by the PVGIS tool [42]. The distribution of the mean absolute percent error (MAPE) for the 10,000 repetitions is reported in Figure 4a, for both the training and the validation steps. The MAPE distribution shows a Gaussian shape, with average values equal to 4.3% (training) and 4.9% (validation).
The fitted coefficients of Equation (1) per each training step are averaged over the 10,000 repetitions, and the resulting mean values and standard deviations are reported in Table 1, where units for w 1 and w 4 are indeed given per degree of latitude (the unit for latitude is the decimal degree, as recommended by ISO 31, instead of using degrees, minutes and seconds [43]). These optimal coefficients are then used to predict H y for the 80 considered locations in Europe and Africa. The predicted global sum of solar irradiation is compared to that provided by the PVGIS tool, and the relative percent error is reported in the last column of Table A1 and Table A2. Figure 4b shows a graphical comparison between the H y values obtained by the proposed model (Equation (1) with the optimal coefficients in Table 1) and those given by the PVGIS tool for all the 80 considered locations. The difference in the global irradiation for Europe (from about 60 to 37 latitude) and Africa (from about 37 to −30 latitude) can be clearly appreciated; nevertheless, no systematic errors with the latitude are noticeable. Results show that the maximum relative error of the model prediction is 13.6% in absolute value, which is within the standard range of estimation errors for the yearly global irradiation. For instance, Solargis [44] reports a data accuracy in the range of ±4% to ±8% for the global horizontal irradiation, and ±8% to ±15% for the direct normal irradiation [45]. The inset of Figure 4b, presents the cumulative distribution function of the percent error between the predicted and the PVGIS values of H y . The MAPE is 4.4% and the normalized root mean square error (NRMSE) in percentage is 5.5%. These errors are similar to the ones typically obtained in other solar irradiation models from the literature. For example, methodologies based on artificial neural networks or inverse distance weighting algorithm have shown mean absolute percent errors equal to 5.9 % , 3.4 % and 4.3 % at a nationwide level (Malaysia [46], Indonesia [47], and Taiwan [48], respectively). Notably, such discrepancies are all consistent with investors requests, namely average errors within 5% in prediction accuracy [49].
This evidence shows that the number of chosen locations is sufficient to train adequate values of w i ( i = 1 , , 5 ) coefficients and, therefore, that an accurately representative model for the yearly sum of global solar irradiation in Europe and Africa has been obtained.

3.2. PV Energy Model (Optimal Azimuth)

The proposed methodology is now applied to estimate the potential electrical energy output by PV system installations for the whole set of considered locations. Initially, H y and thus Π y are computed considering installations with optimal tilt and azimuth, namely α z = 0 and thus F = 1 , and crystalline silicon solar cells. The methodology is implemented in MATLAB®, resulting in a simple and fast algorithm that can be straightforwardly applied for the analysis of large data sets.
First, the model coefficients for the temperature and low-irradiance efficiency in Equation (3) are obtained by fitting the values extracted from PVGIS tool over the 80 locations in Europe and Africa. Best-fitted values for free-standing and building-integrated installations are listed in Table 2. The coefficients of determination of these fittings, that is R2 = 0.92 (free-standing) and R2 = 0.90 (building-integrated), demonstrate that Equation (3) is an accurate correlation between η temp and T 24 , as also evident in Figure 5a. Second, the reflectance efficiency of Equation (4) is fitted on the 80 European and African locations as well. Since, in this case, we consider optimal azimuthal angle, Equation (4) simplifies to η refl = q 5 , being q 5 = 0.972 the best-fitted value. As reported in Figure 5b, the η refl values from the PVGIS tool range from 0.967 to 0.978, namely within a ± 0.6 % from η refl = 0.972 . In fact, η refl demonstrates only a slight dependence on latitude.
The Π y values are finally computed by Equation (6), considering an indicative efficiency for the crystalline silicon solar cells equal to η mod = 0.25 [50], a system installation efficiency equal to η inst = 0.84 , and a unitary surface exposed to the sun. The predicted Π y values for the 80 locations in Europe and Africa are compared to those provided by PVGIS tool, and the relative percent error is reported in Figure 6a,b (free-standing installations), and in Figure 6c,d (BI installations), respectively. Results show that the maximum relative error for FS installations is 14.8% in absolute value, MAPE is 4.3% and NRMSE is 5.5%. Similar errors are obtained for BI solutions, namely a maximum relative error equal to 15.4% in absolute value, MAPE equal to 4.3% and NRMSE to 5.5%. These values are essentially coherent with the prediction errors of H y discussed in the previous Section, which therefore represent the main source of uncertainty in the Π y estimation. These uncertainties are consistent with the general accuracy provided by other tools for estimating the PV potential throughout large territories. For example, PVWatts tool shows errors as high as ± 10 % for annual PV electricity production, with actual performance in a specific year up to ± 20 % respect to long-term average [13]. These uncertainties are also similar to those indicated by JRC (i.e., from 20 % to + 5 % ) when the effect of irradiation and temperature on PV module performance are both considered [3].

3.3. Effect of Non-Optimal Azimuth

In this Section we analyze the effect of non-optimal azimuthal angles on the reflectance efficiency η refl of Equation (4) and on the scaling factor F ( α z ) of Equation (5) for various geographical locations in Europe and Africa. The locations are specifically selected to be representative of the whole extent of each continent, namely maximum, minimum, and intermediate latitudes. Using this criterion, the following spots are chosen from the data set: St. Petersburg (Russia), Edinburgh (UK), Rome (Italy) and Sevilla (Spain) for Europe; Rabat (Morocco), Aswan (Egypt), Karthoum (Sudan) and Kampala (Uganda) for Africa. Note that the selected spots also differ significantly in longitude; however, no appreciable effect of this latter parameter has been noticed in the analysis. The raw data for the reflectance efficiency and insolation is obtained for each location using standard queries to the online PVGIS tool, for azimuthal angles varying in the range 90 α z 90 .
The results are shown in Figure 7a,b, where we also compare the solution given by Equations (4) and (5), respectively. The parametrization of these two latter equations was previously obtained for Italy, using the raw data of seven cities throughout the peninsula for the reflectance efficiency and the only data of Rome for the scaling factor [24]. We note that the parametrization for Italy gives well-representative curves also for the other locations in Europe: the maximum relative error is 2.1% on the reflectance efficiency and 13.8% on the scaling factor at the maximum azimuthal angles considered ( ± 90 ). We can then assert that, given the geographical location of Italy at intermediate continental latitudes, this parametrization can be assumed to be representative also for Europe—and the coefficients are reported in Table 3. On the other hand, the results for Africa require a specific treatment, as the decrease in latitude towards the equator significantly modifies both the reflectance efficiency and the scaling factor with respect to those of the European spots. We then proceed with a new parametrization of the models by regression, using the data of the four considered locations in Africa. We obtain the coefficients reported in Table 3. In this case, the maximum relative error is 0.3% on the reflectance efficiency and 6.4% on the scaling factor at the maximum azimuthal angles. Considering that α z = ± 90 is the most unfavorable and limiting situation for an installation, the proposed models provide acceptable approximations for the reflectance efficiency and scaling factor. In particular, according to the results obtained, the parametrization for Europe can be adopted in the range of latitudes 37 ϕ 60 , while that for Africa from ϕ = 37 to the equator ( ϕ = 0 ). Symmetry applies for locations in Africa below the equator (negative latitudes).
An example application of the proposed models and related analysis is reported in Figure 8a,b for a BI installation in the city of Rome (Italy), which has azimuthal angle α z = 35 , η mod = 0.25 for crystalline silicon solar cells [50], η inst = 0.84 and unitary exposed surface. The main contribution to the total losses is given by the module efficiency, then by temperature and installation losses. Note that the losses due to reflectance effects and azimuthal angle contribute only for a small share on the total losses. For Europe, a similarly small contribution of the reflectance losses would be obtained for any azimuthal angle, as they consist in any case of a few percent (see Figure 7a). The contribution of losses due to the azimuthal angle can be expected to be more consistent for large azimuthal angles, e.g., for α z = ± 60 , where losses could reach nearly 10% (see Figure 7b). For Africa, both contributions to the total losses tend to be smaller for any azimuth.

4. Concluding Remarks

In this work, we have proposed a simple methodology for a fast estimation of the yearly sum of solar irradiation in Europe and Africa and the resulting potential PV energy output. The model thus generated is intended to provide: (i) high spatial (continuous) resolution and (ii) a fast means for pre-screening of the PV potential over large-scale territories for any investor. In fact, the presented method relies on just a small set of input parameters. For instance, estimation of the yearly sum of global solar irradiation requires only the knowledge of a few basic geodata, such as latitude, altitude, and average daily temperature. One such tool operates offline and provides a much faster screening than systematic calls to online databases, and can be easily implemented with just a few lines of code.
Limitations of this work include the lack of information at local level, for instance, about the shading effect of buildings or other elements nearby, as well as horizon height. Thus, the reported tool should be considered as a fast (macroscopic) screening of the solar and PV potentials over large-scale territories with tolerable errors in the range ± 5 % , which should then be followed by a more accurate (microscopic) analysis considering also local elements for a limited set of locations. Furthermore, both the solar irradiation and the potential photovoltaic electricity generation are here estimated on an annual basis, namely they consider the cumulative value over the average year. The reason is that the yearly electricity output is a standard figure of merit needed to assess the techno-economic feasibility of a possible photovoltaic system in a given location [31]. Clearly, this information would not be enough to size the whole photovoltaic system, since detailed estimations of daily and hourly variability are required to design, for example, the energy storage system [51]. In perspective, our multiple-regression model could be further improved to consider daily variability as well, for instance by including well-established models from the literature [52].
This approach may find application in the design of micro-grids [53], polygeneration systems [54], BI PV solutions [55] or flexible storage photovoltaic systems [56]. Furthermore, this methodology could also support the large-scale potential assessment of sustainable technologies other than photovoltaics, for example, solar thermal systems [57], solar greenhouses or dryers for food production or conservation [58], and desalination plants driven by solar source [59].

Author Contributions

A.B. conceived the idea of this study. A.B., with help of L.B. (Luca Bergamasco) and M.F., performed the computations, analyzed the results and wrote the manuscript. L.B. (Lorenzo Bottaccioli) prepared the training and validation data sets. A.M., E.C. and P.A. provided methodological support and supervised the study. All authors have read and approved the final version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

Computational resources provided by HPC@POLITO, a project of Academic Computing of the Department of Control and Computer Engineering (DAUIN) at Politecnico di Torino (http://hpc.polito.it).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
PVPhotovoltaic
BIPVBuilding-Integrated Photovoltaic
PVGISPhotovoltaic Geographic Information System
JRCJoint Research Centre (European Commission)
CM-SAFSatellite Application Facility on Climate Monitoring
MAPEMean Absolute Percent Error
NRMSENormalized Root Mean Square Error

Appendix A

Table A1. Locations in Europe and Africa considered to train and validate the multiple-regression method. For each location, latitude ( ϕ ), longitude ( λ ), altitude (h), and mean temperature over 24 h ( T 24 ) are reported. The yearly sum of the global solar irradiation ( H y ) as obtained from the PVGIS tool [37] and the proposed estimator are compared. The last column reports the relative error of the estimates with respect to the PVGIS data. Note that the exact geographical coordinates have been chosen based on the availability of the required data, and the indicated cities/places represent the most significant locations nearby the considered geographical points.
Table A1. Locations in Europe and Africa considered to train and validate the multiple-regression method. For each location, latitude ( ϕ ), longitude ( λ ), altitude (h), and mean temperature over 24 h ( T 24 ) are reported. The yearly sum of the global solar irradiation ( H y ) as obtained from the PVGIS tool [37] and the proposed estimator are compared. The last column reports the relative error of the estimates with respect to the PVGIS data. Note that the exact geographical coordinates have been chosen based on the availability of the required data, and the indicated cities/places represent the most significant locations nearby the considered geographical points.
LocationCountry ϕ [ ] λ [ ]h [m] T 24 [ C] H y [kWh/m2] Δ [%]
PVGISComputed
EdinburghScotland55.94−3.30449.0114012086.0
VilniusLithuania54.6425.271867.1114011420.2
WarsawPoland52.2021.001108.9124012702.4
LondonEngland51.480.002810.2132013512.4
KievUkraine50.4530.461658.913401306−2.6
PragueCzech Republic50.1214.622809.3127013506.3
ViennaAustria48.1716.3922311.0141014764.7
BudapestHungary47.4719.1512311.515101505−0.3
VaduzLiechtenstein47.149.5045410.3142014774.0
BolzanoItaly46.4711.3223814.117401725−0.8
ZagrebCroatia45.8115.9712711.7150015392.6
BelgradeSerbia44.8020.388013.0159016332.7
BucharestRomania44.4326.009011.916401563−4.7
SofiaBulgaria42.6323.4157510.316401557−5.1
RomeItaly41.9712.535416.419401914−1.3
TiranaAlbania41.3619.8011116.4189019231.8
YerevanArmenia40.1644.52101113.619501844−5.4
LisbonPortugal38.75−9.158916.321701919−11.6
SevilleSpain37.38−5.951418.421802076−4.8
TunisiTunisia36.7410.241618.6209020910.0
GibraltarGibraltar36.15−5.35417.720502017−1.6
RabatMorocco33.96−6.877517.522002011−8.6
CairoEgypt29.7431.389622.723902368−0.9
AswanEgypt23.3132.3324027.4256025830.9
MoptiMali15.27−4.1726129.8227024206.6
KaolackSenegal13.66−15.69028.622802278−0.1
OuagadougouBurkina Faso12.05−1.6431929.0225022731.0
DjiboutiRepublic of Djibouti12.0542.7494230.6237024001.3
Dire DawaEthiopia10.4541.7367728.924902258−9.3
Addis AbabaEthiopia8.8438.11237919.0214023318.9
AccraGhana5.630.00927.121401985−7.2
BanguiCentral African Republic4.0218.4837625.520801997−4.0
DoualaCameroon4.0210.4532725.3188019925.9
KisanganiDR Congo0.8024.9144525.7185019233.9
MombasaKenya−4.0239.3821926.721301967−7.6
BrazzavilleRepublic of the Congo−4.0215.2739924.5185020078.5
HuamboAngola−13.6615.69156221.2226022881.3
HarareZimbabwe−18.4830.42128120.823402285−2.3
MaputoMozambique−26.5232.244821.91990226013.6
JohannesburgSouth Africa−26.5228.66160216.122502149−4.5
Table A2. Locations in Europe and Africa considered to train and validate the multiple-regression method. For each location, latitude ( ϕ ), longitude ( λ ), altitude (h), and mean temperature over 24 h ( T 24 ) are reported. The yearly sum of the global solar irradiation ( H y ) as obtained from the PVGIS tool [37] and the proposed estimator are compared. The last column reports the relative error of the estimates with respect to the PVGIS data. Note that the exact geographical coordinates have been chosen based on the availability of the required data, and the indicated cities/places represent the most significant locations nearby the considered geographical points.
Table A2. Locations in Europe and Africa considered to train and validate the multiple-regression method. For each location, latitude ( ϕ ), longitude ( λ ), altitude (h), and mean temperature over 24 h ( T 24 ) are reported. The yearly sum of the global solar irradiation ( H y ) as obtained from the PVGIS tool [37] and the proposed estimator are compared. The last column reports the relative error of the estimates with respect to the PVGIS data. Note that the exact geographical coordinates have been chosen based on the availability of the required data, and the indicated cities/places represent the most significant locations nearby the considered geographical points.
LocationCountry ϕ [ ] λ [ ]h [m] T 24 [ C] H y [kWh/m2] Δ [%]
PVGISComputed
St. PetersburgRussia59.9830.46185.41070941−12.1
MoscowRussia55.6537.631706.011601070−7.7
BerlinGermany52.5413.52559.7125013074.6
The HagueNetherlands52.064.36010.3132013462.0
BrusselsBelgium50.864.375410.41250137710.2
FrankfurtGermany50.138.7013310.5128014049.7
ParisFrance48.912.373811.4137014697.2
StuttgartGermany48.809.2024210.5132014388.9
BratislavaSlovakia48.1117.0613411.2146014781.2
BernSwitzerland46.967.435718.714401407−2.3
LjubljanaSlovenia46.0814.4831111.5146015476.0
TurinItaly45.117.7321012.917201641−4.6
BordeauxFrance44.79−0.53413.5160016603.7
SarajevoBosnia-Herzegovina43.8318.3451410.3150015322.1
PodgoricaMontenegro42.4219.264817.3188019915.9
IstanbulTurkey41.0728.779014.9180018030.2
MadridSpain40.35−3.7361514.520401851−9.3
AthensGreece37.9823.703018.421202080−1.9
SyracuseItaly37.2014.9534816.319901958−1.6
AlgiersAlgeria36.703.101219.121402132−0.4
NicosiaCyprus35.1433.3817120.3224022450.2
TripoliLibya32.7613.175121.122802273−0.3
NouakchottMauritania18.48−15.211228.2231024375.5
KarthoumSudan15.2732.5037930.3235024564.5
AsmaraEritrea15.2739.17133725.4225024056.9
BafatáGuinea Bissau12.05−14.79728.022302205−1.1
KanoNigeria12.058.2251427.222702255−0.7
KadunaNigeria10.457.3658926.8216022082.2
MoundouChad8.8415.4142027.822302149−3.6
BouafléCôte d’Ivoire7.23−5.6818526.521002056−2.1
JubaSouth Sudan4.0231.3480326.721802047−6.1
KampalaUganda0.8032.95107524.121102038−3.4
NgomaRwanda−2.4129.73169219.91970220111.7
Tarangire N.P.Tanzania−4.0236.16117522.1209021291.8
MbeyaTanzania−8.8433.24126120.922602194−2.9
LilongweMalawi−13.6633.85143820.7218022603.7
LusakaZambia−15.2727.50103021.423002238−2.7
AntananarivoMadagascar−18.4847.32139220.0217022714.7
WindhoekNamibia−23.3116.60179421.525602437−4.8
BloemfonteinSouth Africa−29.7425.85141116.824302152−11.4

References

  1. IEA World Energy Outlook 2017. Available online: http://www.iea.org/weo/ (accessed on 14 June 2018).
  2. Kyoto Protocol to the United Nations Framework Convention on Climate Change; U.N. Doc FCCC/CP/1997/7/ Add.1, 37 I.L.M. 22. Kyoto, Japan, 1998. Available online: http://www.globaldialoguefoundation.org/files/ENV.2009-jun.unframeworkconventionclimate.pdf (accessed on 11 December 2018).
  3. Huld, T.; Gracia Amillo, A.M. Estimating PV Module Performance over Large Geographical Regions: The Role of Irradiance, Air Temperature, Wind Speed and Solar Spectrum. Energies 2015, 8, 5159–5181. [Google Scholar] [CrossRef] [Green Version]
  4. Directive 2009/28/EC of the European Parliament and of the Council. 2009. Available online: http://data.europa.eu/eli/dir/2009/28/oj (accessed on 11 December 2018).
  5. Proposal for a Directive COM 2016/767/F2 of the European Parliament and of the Council. 2016. Available online: https://ec.europa.eu/transparency/regdoc/rep/1/2016/EN/COM-2016-767-F2-EN-MAIN-PART-1.PDF (accessed on 11 December 2018).
  6. Martínez-Rubio, A.; Sanz-Adan, F.; Santamaría-Peña, J.; Martínez, A. Evaluating solar irradiance over facades in high building cities, based on LiDAR technology. Appl. Energy 2016, 183, 133–147. [Google Scholar] [CrossRef]
  7. Numbi, B.P.; Malinga, S.J. Optimal energy cost and economic analysis of a residential grid-interactive solar PV system-case of eThekwini municipality in South Africa. Appl. Energy 2017, 186, 28–45. [Google Scholar] [CrossRef]
  8. Micangeli, A.; Del Citto, R.; Santori, S.G.; Gambino, V.; Kiplagat, J.; Viganò, D.; Poli, D.; Fioriti, D.; Del Citto, R.; Kiva, I.N. Energy Production Analysis and Optimization of Mini-Grid in Remote Areas: The Case Study of Habaswein, Kenya. Energies 2017, 10, 2041. [Google Scholar] [CrossRef]
  9. International Renewable Energy Agency (IRENA). Solar PV in Africa: Costs and Markets. 2016. Available online: https://www.irena.org/-/media/Files/IRENA/Agency/Publication/2016/IRENA_Solar_PV_Costs_Africa_2016.pdf (accessed on 11 December 2018).
  10. Quansah, D.A.; Adaramola, M.S.; Mensah, L.D. Solar Photovoltaics in sub-Saharan Africa—Addressing Barriers, Unlocking Potential. Energy Procedia 2016, 106, 97–110. [Google Scholar] [CrossRef]
  11. Rus-Casas, C.; Aguilar, J.D.; Rodrigo, P.; Almonacid, F.; Pérez-Higueras, P.J. Classification of methods for annual energy harvesting calculations of photovoltaic generators. Energy Convers. Manag. 2014, 78, 527–536. [Google Scholar] [CrossRef]
  12. Yousif, J.H.; Kazem, H.A.; Boland, J. Predictive models for photovoltaic electricity production in hot weather conditions. Energies 2017, 10, 971. [Google Scholar] [CrossRef]
  13. Dobos, A.P. PVWatts Version 5 Manual; Technical Report; National Renewable Energy Laboratory (NREL): Golden, CO, USA, 2014.
  14. Šúri, M.; Huld, T.A.; Dunlop, E.D. PV-GIS: A web-based solar radiation database for the calculation of PV potential in Europe. Int. J. Sustain. Energy 2005, 24, 55–67. [Google Scholar] [CrossRef]
  15. International Renewable Energy Agency (IRENA). Global Atlas for Renewable Energy: Overview of Solar and Wind Maps. 2014. Available online: https://www.irena.org/-/media/Files/IRENA/Agency/Publication/2014/GA_Booklet_Web.pdf (accessed on 11 December 2018).
  16. SOLARGIS. Available online: https://solargis.info/ (accessed on 14 June 2018).
  17. Huld, T.; Müller, R.; Gambardella, A. A new solar radiation database for estimating PV performance in Europe and Africa. Sol. Energy 2012, 86, 1803–1815. [Google Scholar] [CrossRef]
  18. Besharat, F.; Dehghan, A.A.; Faghih, A.R. Empirical models for estimating global solar radiation: A review and case study. Renew. Sustain. Energy Rev. 2013, 21, 798–821. [Google Scholar] [CrossRef]
  19. Zhang, J.; Zhao, L.; Deng, S.; Xu, W.; Zhang, Y. A critical review of the models used to estimate solar radiation. Renew. Sustain. Energy Rev. 2017, 70, 314–329. [Google Scholar] [CrossRef]
  20. Erbs, D.G.; Klein, S.A.; Duffie, J.A. Estimation of the diffuse radiation fraction for hourly, daily and monthly-average global radiation. Sol. Energy 1982, 28, 293–302. [Google Scholar] [CrossRef]
  21. Noorian, A.M.; Moradi, I.; Kamali, G.A. Evaluation of 12 models to estimate hourly diffuse irradiation on inclined surfaces. Renew. Energy 2008, 33, 1406–1412. [Google Scholar] [CrossRef]
  22. Perez, R.; Ineichen, P.; Seals, R.; Michalsky, J.; Stewart, R. Modeling daylight availability and irradiance components from direct and global irradiance. Sol. Energy 1990, 44, 271–289. [Google Scholar] [CrossRef]
  23. Skoplaki, E.; Palyvos, J.A. On the temperature dependence of photovoltaic module electrical performance: A review of efficiency/power correlations. Sol. Energy 2009, 83, 614–624. [Google Scholar] [CrossRef]
  24. Bocca, A.; Chiavazzo, E.; Macii, A.; Asinari, P. Solar energy potential assessment: An overview and a fast modeling approach with application to Italy. Renew. Sustain. Energy Rev. 2015, 49, 291–296. [Google Scholar] [CrossRef] [Green Version]
  25. Bocca, A.; Bottaccioli, L.; Chiavazzo, E.; Fasano, M.; Macii, A.; Asinari, P. Estimating photovoltaic energy potential from a minimal set of randomly sampled data. Renew. Energy 2016, 97, 457–467. [Google Scholar] [CrossRef]
  26. Pfenninger, S.; Staffell, I. Long-term patterns of European PV output using 30 years of validated hourly reanalysis and satellite data. Energy 2016, 114, 1251–1265. [Google Scholar] [CrossRef]
  27. Pfenninger, S. Dealing with multiple decades of hourly wind and PV time series in energy models: A comparison of methods to reduce time resolution and the planning implications of inter-annual variability. Appl. Energy 2017, 197, 1–13. [Google Scholar] [CrossRef]
  28. Bergamasco, L.; Asinari, P. Scalable methodology for the photovoltaic solar energy potential assessment based on available roof surface area: Application to Piedmont Region (Italy). Sol. Energy 2011, 85, 1041–1055. [Google Scholar] [CrossRef] [Green Version]
  29. Bergamasco, L.; Asinari, P. Scalable methodology for the photovoltaic solar energy potential assessment based on available roof surface area: Further improvements by ortho-image analysis and application to Turin (Italy). Sol. Energy 2011, 85, 2741–2756. [Google Scholar] [CrossRef] [Green Version]
  30. ESRI World Imagery. Available online: https://www.arcgis.com/ (accessed on 14 June 2018).
  31. Duffie, J.A.; Beckman, W.A. Solar Engineering of Thermal Processes; John Wiley & Sons: Hoboken, NJ, USA, 2013. [Google Scholar]
  32. Aldobhani, A.A.M.S. Effect of Altitude and Tilt Angle on Solar Radiation in Tropical Regions. J. Sci. Technol. 2014, 19, 96–109. [Google Scholar]
  33. Wong, L.; Chow, W. Solar radiation model. Appl. Energy 2001, 69, 191–224. [Google Scholar] [CrossRef]
  34. Dubey, S.; Sarvaiya, J.N.; Seshadri, B. Temperature Dependent Photovoltaic PV Efficiency and Its Effect on PV Production in the World—A Review. Energy Procedia 2013, 33, 311–321. [Google Scholar] [CrossRef]
  35. Bocca, A.; Bottaccioli, L.; Macii, A. Temperature Efficiency Analysis in Photovoltaics Using Basic Geodata: Application to Europe. In Proceedings of the 2018 International Symposium on Power Electronics, Automation and Motion (SPEEDAM), Amalfi, Italy, 20–22 June 2018; pp. 823–828. [Google Scholar]
  36. Martin, N.; Ruiz, J.M. Calculation of the PV modules angular losses under field conditions by means of an analytical model. Sol. Energy Mater. Sol. Cells 2001, 70, 25–38. [Google Scholar] [CrossRef]
  37. PVGIS, Online Tool. Available online: http://re.jrc.ec.europa.eu/pvgis (accessed on 14 June 2018).
  38. Šúri, M.; Huld, T.; Dunlop, E.; Ossenbrink, H. Potential of solar electricity generation in the European Union member states and candidate countries. Sol. Energy 2007, 81, 1295–1305. [Google Scholar] [CrossRef]
  39. Gracia Amillo, A.M.; Huld, T. Performance Comparison of Different Models for the Estimation of Global Irradiance on Inclined Surfaces. Validation of the Model Implemented in PVGIS; Technical Report EUR 26075 EN, JRC81902; 2013; ISBN 978-92-79-32507-6. Available online: https://publications.europa.eu/en/publication-detail/-/publication/4ef8c4e1-4397-4e27-8487-448786327f27 (accessed on 11 December 2018).
  40. Urraca, R.; Gracia-Amillo, A.M.; Koubli, E.; Huld, T.; Trentmann, J.; Riihelä, A.; Lindfors, A.V.; Palmer, D.; Gottschalg, R.; Antonanzas-Torres, F. Extensive validation of CM SAF surface radiation products over Europe. Remote Sens. Environ. 2017, 199, 171–186. [Google Scholar] [CrossRef] [Green Version]
  41. Berkeley Earth. Available online: http://berkeleyearth.lbl.gov (accessed on 12 September 2018).
  42. Zhang, J.; Hodge, B.M.; Florita, A.; Lu, S.; Hamann, H.F.; Banunarayanan, V. Metrics for Evaluating the Accuracy of Solar Power Forecasting. In Proceedings of the 3rd International Workshop on Integration of Solar Power into Power Systems, London, UK, 21–22 October 2013; pp. 1–19. [Google Scholar]
  43. Taylor, B.N.; Thompson, A. The International System of Units (SI)—NIST Special Publication 330; National Institute of Standards and Technology: Gaithersburg, MD, USA, 2008.
  44. Šúri, M.; Cebecauer, T.; Skoczek, A. SolarGIS: Solar data and online applications for PV planning and performance assessment. In Proceedings of the 26th European Photovoltaics Solar Energy Conference, Hamburg, Germany, 5–9 September 2011. [Google Scholar]
  45. Global Solar Atlas: Accuracy. Available online: http://globalsolaratlas.info/knowledge-base/accuracy? (accessed on 12 September 2018).
  46. Khatib, T.; Mohamed, A.; Sopian, K.; Mahmoud, M. Solar energy prediction for Malaysia using artificial neural networks. Int. J. Photoenergy 2012, 2012, 419504. [Google Scholar] [CrossRef]
  47. Rumbayan, M.; Abudureyimu, A.; Nagasaka, K. Mapping of solar energy potential in Indonesia using artificial neural network and geographical information system. Renew. Sustain. Energy Rev. 2012, 16, 1437–1449. [Google Scholar] [CrossRef]
  48. Kuo, P.H.; Chen, H.C.; Huang, C.J. Solar Radiation Estimation Algorithm and Field Verification in Taiwan. Energies 2018, 11, 1374. [Google Scholar] [CrossRef]
  49. Williams, S.R.; Gottschalg, R. Accuracy of European Energy Modelling Approaches. In Proceedings of the 21st European Photovoltaic Solar Energy Conference, Dresden, Germany, 4–8 September 2006; pp. 2452–2455. [Google Scholar]
  50. Green, M.A.; Hishikawa, Y.; Dunlop, E.D.; Levi, D.H.; Hohl-Ebinger, J.; Ho-Baillie, A.W.Y. Solar cell efficiency tables (version 51). Prog. Photovolt. Res. Appl. 2017, 26, 3–12. [Google Scholar] [CrossRef] [Green Version]
  51. Aderemi, B.; Chowdhury, S.; Olwal, T.; Abu-Mahfouz, A. Techno-Economic Feasibility of Hybrid Solar Photovoltaic and Battery Energy Storage Power System for a Mobile Cellular Base Station in Soshanguve, South Africa. Energies 2018, 11, 1572. [Google Scholar] [CrossRef]
  52. Gueymard, C. Prediction and performance assessment of mean hourly global radiation. Sol. Energy 2000, 68, 285–303. [Google Scholar] [CrossRef]
  53. Rose, A.; Stoner, R.; Pérez-Arriaga, I. Prospects for grid-connected solar PV in Kenya: A systems approach. Appl. Energy 2016, 161, 583–590. [Google Scholar] [CrossRef]
  54. Mohan, G.; Kumar, U.; Pokhrel, M.K.; Martin, A. A novel solar thermal polygeneration system for sustainable production of cooling, clean water and domestic hot water in United Arab Emirates: Dynamic simulation and economic evaluation. Appl. Energy 2016, 167, 173–188. [Google Scholar] [CrossRef]
  55. Pintér, G.; Baranyai, N.H.; Wiliams, A.; Zsiborács, H. Study of Photovoltaics and LED Energy Efficiency: Case Study in Hungary. Energies 2018, 11, 790. [Google Scholar] [CrossRef]
  56. Zsiborács, H.; Baranyai, N.H.; Vincze, A.; Háber, I.; Pintér, G. Economic and Technical Aspects of Flexible Storage Photovoltaic Systems in Europe. Energies 2018, 11, 1445. [Google Scholar] [CrossRef]
  57. Morciano, M.; Fasano, M.; Secreto, M.; Jamolov, U.; Chiavazzo, E.; Asinari, P. Installation of a concentrated solar power system for the thermal needs of buildings or industrial processes. Energy Procedia 2016, 101, 956–963. [Google Scholar] [CrossRef]
  58. Marucci, A.; Cappuccini, A. Dynamic photovoltaic greenhouse: Energy efficiency in clear sky conditions. Appl. Energy 2016, 170, 362–376. [Google Scholar] [CrossRef]
  59. Morciano, M.; Fasano, M.; Salomov, U.; Ventola, L.; Chiavazzo, E.; Asinari, P. Efficient steam generation by inexpensive narrow gap evaporation device for solar applications. Sci. Rep. 2017, 7, 11970. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (a) Average net capacity additions per year on the global power market by type (N: nuclear, G: gas, C: coal, R: renewables) according to the World Energy Outlook 2017 [1]. (b) Total in-plane solar irradiation for an equator-facing plane inclined 20 from the horizontal (color bar units are in kWh/m 2 ). Figure adapted (that is, cropped from original) from Reference [3], used under CC BY 4.0 license.
Figure 1. (a) Average net capacity additions per year on the global power market by type (N: nuclear, G: gas, C: coal, R: renewables) according to the World Energy Outlook 2017 [1]. (b) Total in-plane solar irradiation for an equator-facing plane inclined 20 from the horizontal (color bar units are in kWh/m 2 ). Figure adapted (that is, cropped from original) from Reference [3], used under CC BY 4.0 license.
Energies 11 03477 g001
Figure 2. Overview of the proposed methodology for fast estimation of solar irradiation and photovoltaic energy potential over large territories.
Figure 2. Overview of the proposed methodology for fast estimation of solar irradiation and photovoltaic energy potential over large territories.
Energies 11 03477 g002
Figure 3. Randomly chosen geographical locations for the analysis (the full list of coordinates is reported in the Appendix A). The data set consists of 80 locations, 40 in Europe and 40 in Africa. Map source: ESRI World Imagery [30].
Figure 3. Randomly chosen geographical locations for the analysis (the full list of coordinates is reported in the Appendix A). The data set consists of 80 locations, 40 in Europe and 40 in Africa. Map source: ESRI World Imagery [30].
Energies 11 03477 g003
Figure 4. (a) Distribution of the mean absolute percent error (MAPE) between the PVGIS and the predicted values of H y , for the 10,000 repetitions of the training/validation process for the coefficients of Equation (1). Red lines refer to the cumulative distribution functions of the MAPE values. (b) Distribution of H y for the 80 considered locations in Europe and Africa: comparison between PVGIS data (black dots) and estimations by Equation (1) with the optimal coefficients in Table 1 (red triangles). Note that the maximum average yearly sum of global solar irradiation is found at the Tropic of Cancer and Capricorn (±23 latitude). In the inset, the cumulative distribution function of the percent estimation errors is shown.
Figure 4. (a) Distribution of the mean absolute percent error (MAPE) between the PVGIS and the predicted values of H y , for the 10,000 repetitions of the training/validation process for the coefficients of Equation (1). Red lines refer to the cumulative distribution functions of the MAPE values. (b) Distribution of H y for the 80 considered locations in Europe and Africa: comparison between PVGIS data (black dots) and estimations by Equation (1) with the optimal coefficients in Table 1 (red triangles). Note that the maximum average yearly sum of global solar irradiation is found at the Tropic of Cancer and Capricorn (±23 latitude). In the inset, the cumulative distribution function of the percent estimation errors is shown.
Energies 11 03477 g004
Figure 5. (a) Temperature and low-irradiance efficiency ( η temp ) as a function of the average daily ambient temperature ( T 24 ) for the given locations: best-fitted models of Equation (3) (lines) are compared with the values taken from the PVGIS tool (dots) for the 80 locations in Europe and Africa shown in Figure 3. Crystalline silicon solar cells are considered. Free-standing (FS) installations are depicted in blue; whereas, building-integrated (BI) installations in red. (b) Reflectance efficiency ( η refl ) as a function of the latitude ϕ for the 80 locations in Europe and Africa: the values taken from the PVGIS tool (dots) and their average value (line) are shown.
Figure 5. (a) Temperature and low-irradiance efficiency ( η temp ) as a function of the average daily ambient temperature ( T 24 ) for the given locations: best-fitted models of Equation (3) (lines) are compared with the values taken from the PVGIS tool (dots) for the 80 locations in Europe and Africa shown in Figure 3. Crystalline silicon solar cells are considered. Free-standing (FS) installations are depicted in blue; whereas, building-integrated (BI) installations in red. (b) Reflectance efficiency ( η refl ) as a function of the latitude ϕ for the 80 locations in Europe and Africa: the values taken from the PVGIS tool (dots) and their average value (line) are shown.
Energies 11 03477 g005
Figure 6. (a) Percent error between the current model predictions of Equation (6) and the values extracted from the PVGIS tool for the yearly electrical energy output ( Π y ) of FS PV systems over the 80 considered locations in Europe and Africa. Results are reported as a function of the latitude ( ϕ ) and longitude ( λ ), as well as in the form of (b) error distribution. (c) Percent error between the current model predictions and the values extracted from the PVGIS tool for the yearly electrical energy output of BI PV systems. Results are reported as function of latitude ( ϕ ) and longitude ( λ ), as well as in the form of (d) error distribution.
Figure 6. (a) Percent error between the current model predictions of Equation (6) and the values extracted from the PVGIS tool for the yearly electrical energy output ( Π y ) of FS PV systems over the 80 considered locations in Europe and Africa. Results are reported as a function of the latitude ( ϕ ) and longitude ( λ ), as well as in the form of (b) error distribution. (c) Percent error between the current model predictions and the values extracted from the PVGIS tool for the yearly electrical energy output of BI PV systems. Results are reported as function of latitude ( ϕ ) and longitude ( λ ), as well as in the form of (d) error distribution.
Energies 11 03477 g006
Figure 7. Results obtained from the PVGIS tool for the reflectance efficiency (a) and scaling factor (b) for the selected locations in Europe (circles) and Africa (triangles). The solid lines correspond to the model in Equation (4) in (a) and Equation (5) in (b) for Europe (red) and Africa (yellow), with the model coefficients reported in Table 3 for the two cases.
Figure 7. Results obtained from the PVGIS tool for the reflectance efficiency (a) and scaling factor (b) for the selected locations in Europe (circles) and Africa (triangles). The solid lines correspond to the model in Equation (4) in (a) and Equation (5) in (b) for Europe (red) and Africa (yellow), with the model coefficients reported in Table 3 for the two cases.
Energies 11 03477 g007
Figure 8. (a) Sample analysis of the various contributions to the total losses ( η = 1 ψ ) for a building-integrated installation (BIPV) in the city of Rome (Italy) with unitary exposed surface. The azimuthal angle considered is α z = 35 ( ψ azim losses), the module efficiency is η mod = 0.25 ( ψ mod losses) for crystalline silicon solar cells [50], and the installation efficiency is η inst = 0.84 ( ψ inst losses). (b) The results obtained show that the main contribution to the total losses depends on the choice of the module technology.
Figure 8. (a) Sample analysis of the various contributions to the total losses ( η = 1 ψ ) for a building-integrated installation (BIPV) in the city of Rome (Italy) with unitary exposed surface. The azimuthal angle considered is α z = 35 ( ψ azim losses), the module efficiency is η mod = 0.25 ( ψ mod losses) for crystalline silicon solar cells [50], and the installation efficiency is η inst = 0.84 ( ψ inst losses). (b) The results obtained show that the main contribution to the total losses depends on the choice of the module technology.
Energies 11 03477 g008
Table 1. Model coefficients of Equation (1) obtained by the (training/validation) fitting procedure on the data sets in Table A1 and Table A2. The standard deviation indicates the variability of the fitted coefficients in the 10,000 different training sets. Note that, in Equation (1), the latitude is intended in decimal degrees.
Table 1. Model coefficients of Equation (1) obtained by the (training/validation) fitting procedure on the data sets in Table A1 and Table A2. The standard deviation indicates the variability of the fitted coefficients in the 10,000 different training sets. Note that, in Equation (1), the latitude is intended in decimal degrees.
ParameterAverage ValueStandard Deviation
w 1 [kWh·m 2 ]−21.569±2.073
w 2 [kWh·m 3 ]0.137±0.031
w 3 [kWh·m 2 · C 2 ]−0.421±0.133
w 4 [kWh·m 2 · C 2 ]0.071±0.003
w 5 [kWh·m 2 ]2119.345±108.680
Table 2. Coefficients of the model for the temperature efficiency in Equation (3) fitted over the 80 considered locations in Europe and Africa (see Figure 3). The reported values refer to crystalline silicon solar cells in either free-standing or building-integrated installations.
Table 2. Coefficients of the model for the temperature efficiency in Equation (3) fitted over the 80 considered locations in Europe and Africa (see Figure 3). The reported values refer to crystalline silicon solar cells in either free-standing or building-integrated installations.
ParameterFree-StandingBuilding-Integrated
p 1 [ C 2 ] 1 . 014 × 10 6 2 . 757 × 10 5
p 2 [ C 1 ] 3 . 430 × 10 3 4 . 598 × 10 3
p 3 [−] 9 . 484 × 10 1 9 . 114 × 10 1
Table 3. Coefficients of the model for the reflectance efficiency η refl in Equation (4) and for the scaling factor F ( α z ) in Equation (5) for Europe as originally reported for Italy in Ref. [24], and those obtained for Africa in the present work.
Table 3. Coefficients of the model for the reflectance efficiency η refl in Equation (4) and for the scaling factor F ( α z ) in Equation (5) for Europe as originally reported for Italy in Ref. [24], and those obtained for Africa in the present work.
η refl F ( α z )
EuropeAfrica EuropeAfrica
q 1 2 . 038 × 10 11 1 . 219 × 10 11 r 1 3 . 729 × 10 9 1 . 437 × 10 9
q 2 3 . 027 × 10 10 4 . 317 × 10 11 r 2 3 . 463 × 10 7 1 . 002 × 10 7
q 3 1 . 193 × 10 6 2 . 690 × 10 7 r 3 1 . 274 × 10 5 9 . 295 × 10 6
q 4 8 . 264 × 10 7 1 . 512 × 10 7 r 4 1 . 650 × 10 4 2 . 933 × 10 5
q 5 9 . 722 × 10 1 9 . 734 × 10 1 r 5 11

Share and Cite

MDPI and ACS Style

Bocca, A.; Bergamasco, L.; Fasano, M.; Bottaccioli, L.; Chiavazzo, E.; Macii, A.; Asinari, P. Multiple-Regression Method for Fast Estimation of Solar Irradiation and Photovoltaic Energy Potentials over Europe and Africa. Energies 2018, 11, 3477. https://doi.org/10.3390/en11123477

AMA Style

Bocca A, Bergamasco L, Fasano M, Bottaccioli L, Chiavazzo E, Macii A, Asinari P. Multiple-Regression Method for Fast Estimation of Solar Irradiation and Photovoltaic Energy Potentials over Europe and Africa. Energies. 2018; 11(12):3477. https://doi.org/10.3390/en11123477

Chicago/Turabian Style

Bocca, Alberto, Luca Bergamasco, Matteo Fasano, Lorenzo Bottaccioli, Eliodoro Chiavazzo, Alberto Macii, and Pietro Asinari. 2018. "Multiple-Regression Method for Fast Estimation of Solar Irradiation and Photovoltaic Energy Potentials over Europe and Africa" Energies 11, no. 12: 3477. https://doi.org/10.3390/en11123477

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop