Hindawi Publishing Corporation

International Journal of Photoenergy

Volume 2006, Article ID 73650, Pages 1–8
DOI 10.1155/IJP/2006/73650

Sizing Stand-Alone Photovoltaic Systems

A. Balouktsis,1 T. D. Karapantsios,2 A. Antoniadis,3 D. Paschaloudis,1 A. Bezergiannidou,1 and N. Bilalis4

1 Department of Informatics & Communications, Technological Educational Institution of Serres, Terma Magnesias,
P.O. Box 62124, Serres, Greece
2 Department of Chemistry, Aristotle University of Thessaloniki, University Box 116, 541 24 Thessaloniki, Greece
3 Department of Natural Resources & Environment, Technological Educational Institution of Crete, Romanou 3, 731 33 Chania, Greece
4 Department of Production Engineering & Management, Technical University of Crete, University Campus, 73100 Chania, Greece

Received 25 February 2005; Revised 28 August 2005; Accepted 5 September 2005

A method of sizing stand-alone photovoltaic systems regarding the reliability to satisfy the load demand, economy of components,
and discharge depth exploited by the batteries is presented in this work. Solar radiation data simulated by an appropriate stochastic
time series model, and not actual measurements, are used in the sizing procedure. This offers two distinct advantages: (a) sizing
can be performed even for locations where no actual data exist, (b) the influence of the variation of the statistical parameters of
solar radiation in sizing can be examined. The method has been applied and tested for several representative locations all over
Greece for which monthly daily average values of solar radiation are given by ELOT (Hellenic Organization of Standardization).

Copyright © 2006 A. Balouktsis et al. This is an open access article distributed under the Creative Commons Attribution License,
which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. INTRODUCTION (radiation distributions, correlations, etc.) of the actual data.

Such stochastic models have been sporadically communi-
An essential first step when designing a stand-alone photo- cated in the past [6–9].
voltaic (PV) system is to determine the total power of the The present study is concerned with sizing stand-alone
PV system and the energy of the batteries in such a way photovoltaic systems for different locations in Greece. To
that the system becomes capable of satisfying the load (de- produce hourly solar radiation data, the stochastic model of
mand) at a given reliability level and minimum cost. Sizing [9] is employed (with some refinements), for it involves the
a PV system is, apparently, a quite complex issue because fewer and, virtually, the more salient statistical parameters
several stochastic parameters are involved having a signifi- of solar radiation (tabulated monthly average daily values
cant contribution, e.g., meteorological data, variation of de- and the respective standard deviations). If tabulated stan-
mand on daily and seasonal bases, and economical consider- dard deviations values do not exist, they can still be es-
ations. timated by the stochastic model of Bent et al. [10]. Cer-
So far, several methods have been proposed in the lit- tain key simplifications in the theoretical problem formu-
erature for sizing stand-alone photovoltaic systems [1–5]. lation in that study allowed the development of a realis-
These methods either rely on tabulated seasonal-monthly av- tic, yet sufficiently simple, model with broad applicability
erage values of solar radiation and employ simplified cal- which can be used as a convenient diagnostic tool for in-
culations, or rely on daily/hourly measured values of solar vestigating the sensitivity of the design of a stand-alone
radiation and incorporate simulation techniques. The sim- photovoltaic system to the variation of the aforementioned
plified methods are generally not satisfactory; especially for statistical parameters. Tabulated monthly average values of
cases where increased reliability is required as regards the daily solar radiation data, needed as input to the simula-
response to the electric consumptions. On the other hand, tion code, are taken by ELOT [11], whereas standard devia-
the simulation methods require detailed information on a tions are estimated from measured daily values of solar radi-
daily or hourly basis for prolonged periods of time (above 10 ation for the Greek cities of Serres, Kavala, Athens, and Ha-
years), which is often not available. An alternative and per- nia. This is done in the absence of original standard devia-
haps more efficient approach when using simulation meth- tion data for most of the locations considered herein. How-
ods is to employ stochastic models for generating solar ra- ever, to partly overcome this uncertainty, the possible ef-
diation data which retain all the essential stochastic features fect of the standard deviation variation to the system design
2 International Journal of Photoenergy

50

40

Frequency
30

20

10

0
−3 −2 −1 0 1
εd

Figure 1: Histogram of εd distribution.

100

80

60
Percent

40

20

0
−4 −3 −2 −1 0 1
εd

Shape 47.50 Thresh −38.22

Scale 38.17 N 100

Figure 2: 3-parameter Weibull distribution approximation of original εd distribution.

is thoroughly examined. In addition, it is investigated how clearness index Kd can be described by a time series:
the batteries’ discharge depth affects the optimum design
considering that it is directly associated with their durabil- Kd = µd + σd · εd , (2)
ity.
It is of great importance to realize prior to PV sizing
where µd = E(Kd ) is the expected Kd value at day d, σd =
whether the system design is accurate enough when it is 
based on simulated data alone. In pursuing this goal, a com- E[(Kd − µd )2 ] is the standard deviation of Kd at day d, and
parison is made between designs based on actual measure- εd is a static (error) time series with E(εd ) = 0, E(εd 2 ) = 1
ments of solar radiation in Serres and simulated data com- with a correlation coefficient among its successive elements
puted from the average values of the actual measurements of ρ = 0.23 and a distribution as displayed in Figure 1.
for a period of a year. The above distribution can be approximated by a 3-
parameter Weibull distribution as it is shown in Figure 2.
2. STOCHASTIC MODEL OF SOLAR RADIATION This is a refined approach with respect to what was used
in [9].
If the daily clearness index is defined as The daily values of µd for the whole year are derived
from tabulated average monthly values (12 values per year)
Gd
Kd = (1) by simple linear interpolation. The same applies also for the
Go,d daily σd .
(Gd is the daily solar radiation and Go,d the extraterrestrial In cases that there are no available σd data, these val-
solar radiation for the same day), then according to [9] the ues can be obtained from Figure 3 which was deduced from
A. Balouktsis et al. 3

0.2

Standard deviation of mean KT

0.18

0.16

0.14

0.12

0.1
0.3 0.4 0.5 0.6 0.7
Mean KT

Figure 3: Estimated standard deviations, σd , for different mean clearness indexes.

Ed bus
Gd
Controller
Eload
PV
units nc ni
ndeg nw1 nw2
nPV Inverter
nw3
+
nb
−

Figure 4: Schematic representation of a stand-alone PV system. nPV : efficient factor of PV units, nw1 : line efficient factor between PV-
converter, nw2 : line efficient factor between inverter-load, ni : efficient factor of inverter, nc : efficient factor of converter, nw3 : line efficient
factor between converter-batteries, nb : efficient factor of batteries, ndeg : degradation coefficient of PV units.

the clearness index distributions in [10]. The corresponding units, Gt . Total irradiance on the PV array plane is calculated
standard deviations are estimated from the relation: using an isotropic model for both the diffuse irradiation and

 Kd,max the ground reflected irradiation. Calculations by this model
take into account the location latitude. Details on the above

σd =  (Kd − µd )2 pdfµd (Kd )dKd (3)
Kd,min can be found in classic solar engineering books, e.g., [12];
(b) The energy of the photovoltaic, Ed , is estimated for
3. SIMULATION METHODOLOGY the period of a day by the relation

A typical stand-alone photovoltaic system consists of a num- Ed = ndeg nPV Gt . (4)

ber of PV units connected in series and in parallel, batteries, a Then, the energy at the level of the alternating current bal-
converter or controller for supervising the charging and dis- ance (bus point) is calculated from the formula
charging of the batteries, and, finally, one or more inverters
connected in parallel, for the case that an alternating current Ed,load = Ed nw1 nw2 ni . (5)
balance is included in the system. Figure 4 illustrates better
the energy approach that is incorporated in the sizing pro- Next, this value is compared to the daylight energy demand
cedure where, in addition to the individual units, the corre- Elight,load where two different cases are considered:
sponding efficiency coefficients are also depicted. (b1) if Ed,load ≥ Elight,load , the batteries are charged with en-
In the simulation, the energy balance between produc- ergy equal to
tion and demand is inspected during the daytime in two
nw3 nc nb
stages: one from sunrise to sunset and another for the rest of (Ed,load − Elight,load ) (6)
nw2 ni
the day. The most important inspection steps are as follows.
(a) The daily solar irradiation is converted from the with the simulation checking at the same time for the
horizontal, Gh , to the inclined plane of the photovoltaic batteries not to be overcharged,
4 International Journal of Photoenergy

Interpolation
Gt
Simulations for R11· · · R1m Cb,11· · ·Cb,1m
PPV Results . . Contours
PPV . PPV . Cbat
different pairs . .
Cbat Rn1· · ·Rnm Cb,n1· · ·Cb,nm
PPV -Cbat Cbat R PPV

Figure 5: Block diagram of simulation methodology to produce Cbat versus PPV curves.

(b2) if Ed,load < Elight,load , the system does not satisfy the de- where PPV is the peak PV power in kWp and Cbat the capacity
mand for a period of time equal to of batteries in kWh.
τd The whole procedure is presented in Figure 5.
(Elight,load − Ed,load ) , (7)
Elight,load
where τd is the duration of the daylight which has a 4. DESIGN CONSIDERATIONS: PARAMETRIC
different value for every day of the year; INVESTIGATION

(c) For the night hours, the energy demand Enight,load is In the present study, the proposed stochastic time series
satisfied exclusively by the batteries and, therefore, they are model has been utilized to generate solar radiation data be-
discharged by an amount of energy equal to cause original measurements are not available for most of
Enight,load the locations encountered in the calculations. The success
, (8) of this effort relies largely on whether the theoretically gen-
nw1 nw3 ni
erated data resemble sufficiently the lacking measurements
where Eload = Elight,load + Enight,load . or, in other words, whether their statistical parameters ade-
The simulation checks at the same time for the batter- quately approximate the respective parameters of the missing
ies not to be discharged below the minimum depletion level actual data. Therefore, a comparison is made first between
which is determined by the discharge depth. For the case that the simulated solar radiation values and measured data for
the balance of the energy at bus point requires discharge of a period of a year. The measured time series data and those
the batteries below their minimum level, then the system is determined from the simulation are displayed in Figures 6a
considered as not satisfying the demand for a period of time and 6b, respectively. Table 1 presents the corresponding av-
proportional to the energy not covered. erage values and standard deviations.
Following the above steps, the behavior of the system is As can be seen, the maximum error in the average value
examined separately for every different location and for an is less than 9%. Contrasting the sizing results obtained using
extended period of 10 years. Specifically, for every pair of PV the actual and the simulated data for Serres, at two differ-
power-batteries energy, the reliability index of the system is ent reliability levels (0.9, 0.99), the following are observed:
calculated, defined as the ratio of the hours satisfying the de- at reliability 0.9 the design is almost identical (deviation less
mand over the total hours of simulation: than 1%) and at reliability 0.99 the deviation is less than 8%,
hours satisfying the load the data from the simulation producing always higher sizing
R= . (9) values.
total hours of simulation
Next, the possible effect of the variation of standard de-
Running the simulation code for an individual pair of PV viation (σ) of the radiation on the optimum design of the
power-batteries energy, a single value for the reliability in- system is investigated. Figure 7 demonstrates the solar radia-
dex is calculated. For a series of different pairs, a table is tion time series produced by simulation with 0.2σ, σ, and 2σ
constructed of reliability indexes versus pairs of PV power- standard deviation for the city of Athens.
batteries energy. The resolution of this table can be made as
fine as desired in order to attain the required accuracy. Any The mean percentage difference between the cases of 0.2σ
PV power value is then interpolated among the elements of and σ is: (a) for the economical components −30%, (b) for
this table to obtain the corresponding value of the batter- the capacity of the batteries −50%, (c) for the power of the
ies, capacity for the different reliability levels of the system PV −15%. When comparing the cases of 2σ and σ the respec-
specified by the user. In other words, points of the function tive values are: (a) for the economical components 16%, (b)
f (PPV , C) = R are evaluated. Subsequently, using nonlin- for the capacity of the batteries 50%, (c) for the power of the
ear least squares numerical methods, the above functions are PV 1%.
materialized into the general form: Finally, the effect of the discharge depth of the batteries
is examined as regards the design procedure. For this pur-
a1 (R) + b1 (R) · (PPV )2 pose, some necessary information from the literature [14] is
= Cbat , (10)
1 + b2 (R) · (PPV )2 utilized as reproduced in Table 2.
A. Balouktsis et al. 5

8 8

6 6

4 4

2 2

50 100 150 200 250 300 350 200 400 600 800 1000 1200 1400
(a) (b)

Figure 6: (a) Solar radiation data (kWh) measured in Serres (Greece) for one year, (b) solar radiation values simulated for Serres (Greece)
for a period of 4 years.

Table 1: Comparison between actual and simulated monthly average solar radiation data for Serres (Greece).
Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec
Gd , actual,
1.25 1.95 2.8 5.02 6.33 7.14 6.12 5.82 4.6 3.63 2.05 1.03
kWh
Gd , simul,
1.21 1.8 2.78 4.6 5.92 6.6 5.91 5.5 4.42 3.32 1.9 0.99
kWh
Rel. error
3.2 7.7 0.7 8.4 6.5 7.6 3.4 5.5 3.9 8.5 7.3 3.9
%

8
8

6 6

4 4

2 2

50 100 150 200 250 300 350 50 100 150 200 250 300 350
(a) (b)

50 100 150 200 250 300 350

(c)

Figure 7: Yearly solar radiation (kWh) for (a) 0.2σ, (b) σ, and (c) 2σ standard deviation for the city of Athens.
6 International Journal of Photoenergy

40 40
36 36
32 32
28 28
24 24
20 20
16 16
12 12
8 8
4 4
0 0
0.5 0.7 0.9 1.1 1.3 1.5 1.7 1.9 2.1 2.3 2.5 0.5 0.67 0.84 1.01 1.18 1.35 1.52 1.69 1.86 2.03 2.2
(a) (b)

40
36
32
28
24
20
16
12
8
4
0
0.5 0.63 0.76 0.89 1.02 1.15 1.28 1.41 1.54 1.67 2.8
(c)

Figure 8: PV power-batteries energy curves at two different reliability levels (0.9, 0.99) for (a) Serres, (b) Athens, and (c) Hania.

Table 2: Battery life cycle (years) versus discharge depth (%) for Table 3: Economical parameters [13].
two different battery types.
PV Battery
Life cycle, years Discharge depth, % Purchase cost 12 $/Wp 120 $/kWh
For battery For battery Maintance 0$ 0$
24–48 V > 120 V
1 50 40
5 35 25 economical considerations, costs are taken from Table 3 ac-
10 20 15 cording to the work of Lasnier and Ang [13].
15 15 10 Furthermore, the PV cost recently announced by the
American Department of Energy, 5 $/Wp, is also employed
for the sake of comparison. The discharge depth is assumed
The sizing examination is conducted for PV systems for to be constant at 35% which results in a batteries’ life cycle of
some Greek cities for which measured solar radiation data 5 years.
could be found. In all cases, the most economical design (for Figure 8 displays PV batteries curves of equal reliability as
a 30-year life cycle of the total PV system) is attained by using well as optimum points for the Greek cities of Serres, Athens,
batteries with discharge depths <= 35. This value is eventu- and Hania.
ally selected for the calculations of the present study. The PV sizing method and the stochastic radiation time
series model (both described above) combined with tab-
5. RESULTS ulated monthly radiation values were successfully applied
for sizing stand-alone PV systems at different locations in
For the course of the present calculations, the daily electric Greece.
demand is taken as 1 kWh, constant throughout the year. The Table 4 summarizes (a) information as provided by
tilt of the photovoltaic units is set at a value equal to the ELOT [11] on the latitude and minimum monthly aver-
latitude of each examined location plus 20 degrees. For all age daily solar radiation on the tilted plane Gd,t , (b) battery
A. Balouktsis et al. 7

Table 4: Sizing PV systems for different locations in Greece.

12 $/Wp 5 $/Wp
Reliability Min Gd,t Cbat PVpower Cbat PVpower
Location Latitude Zones
index (kWh) (kWh) (Wp) (kWh) (Wp)
0.9 5.6 803 4.7 887
Serres 0.99 41◦ 07′ 1.405 17.7 1601 14.7 1881
0.9 4.4 754 3.63 824 G
Thessaloniki 0.99 40◦ 33′ 2.331 16 1272 13.56 1496
0.9 5.24 787 4.24 878
Ioannina 0.99 39◦ 42′ 2.015 16.34 1193 13.8 1424
0.9 3.89 675 3.21 738
Lamia 0.99 38◦ 54′ 2.535 9.35 983 7.74 1130 F
0.9 5.54 790 4.46 890
Corfu 0.99 39◦ 47′ 2.378 15.56 1240 12.9 1490
0.9 5.26 820 4.2 920
Aliartos 0.99 38◦ 23′ 2.14 11.3 1150 9.15 1350
E
0.9 5.3 780 4.25 880
Limnos 0.99 39◦ 55′ 1.768 9.17 1260 7.13 1450
0.9 4 665 3.2 743
Athens 0.99 37◦ 58′ 2.304 9.4 1142 7.5 1324
0.9 5.08 660 4 760 D
Araxos 0.99 38◦ 10′ 2.73 8.5 1050 6.7 1220
0.9 4.62 760 3.4 840
Mytilini 0.99 39◦ 06′ 2.24 11.46 1190 9.34 1410
0.9 5 620 4.2 720
Kalamata 0.99 37◦ 04′ 2.625 9.5 1070 7.56 1250
C
0.9 4.77 760 3.6 870
Hios 0.99 38◦ 22′ 2.27 10.65 1150 8.5 1350
0.9 4.08 617 3.2 694
Hania 0.99 35◦ 30′ 2.746 11.8 1114 9.7 1310
B
0.9 5.5 600 4.4 700
Naxos 0.99 37◦ 06′ 2.29 9.8 1040 7.8 1230
0.9 3.5 620 2.8 690
A
Rhode 0.99 36◦ 23′ 2.841 11.13 1030 8.98 1230

capacity Cbat , (kWh), (c) the power of the PV given in Wp. influences largely affects the optimum design considering
The above are provided for two PV costs, namely, 12 $/Wp that it is directly associated with their durability. The pro-
and 5 $/Wp. Classification of the various locations into dif- posed method is the first step in an effort to produce a sim-
ferent zones is performed according to the radiation levels as plified generalized procedure for sizing stand-alone PV sys-
reported in [15]. tems.

6. CONCLUSIONS ACKNOWLEDGMENT
The proposed method of sizing stand-alone photovoltaic sys- This work has been founded partially by the Greek Ministry
tems in combination with an appropriate stochastic time se- of Education (program ARCHIMIDES I).
ries model for solar radiation offers two distinct advantages
in sizing stand-alone PV systems: (a) sizing can be performed
even for locations where no actual data exist, (b) the influ- REFERENCES
ence of the variation of the statistical parameters of solar
[1] L. Barra, S. Catalanotti, F. Fontana, and F. Lavorante, “An ana-
radiation in sizing can be examined. The method has been lytical method to determine the optimal size of a photovoltaic
successfully applied for several representative locations all plant,” Solar Energy, vol. 33, no. 6, pp. 509–514, 1984.
over Greece. The capacities of the batteries incorporated in [2] H. Saha, “Design of a photovoltaic electric power system for
the system are seen. In addition, it is investigated how the an indian village,” Solar Energy, vol. 27, no. 2, pp. 103–107,
batteries’ discharge depth that the variance of solar radiation 1981.
8 International Journal of Photoenergy

[3] M. Buresch, Photovoltaic Energy Systems: Design and Installa-

tion, McGraw-Hill, New York, NY, USA, 1983.
[4] P. P. Groumpos and G. Papageorgiou, “An optimal sizing
method for stand-alone photovoltaic power systems,” Solar
Energy, vol. 38, no. 5, pp. 341–351, 1987.
[5] E. Ofry and A. Braunstein, “The loss of power supply proba-
bility as a technique for designing stand-alone solar electrical
(photovoltaic) systems,” IEEE Transactions on Power Appara-
tus and Systems, vol. 102, pp. 1171–1175, 1983.
[6] B. J. Brinkworth, “Autocorrelation and stochastic modelling of
insolation sequences,” Solar Energy, vol. 19, no. 4, pp. 343–347,
1977.
[7] C. Mustacchi, V. Cena, and M. Rocchi, “Stochastic simulation
of hourly global radiation sequences,” Solar Energy, vol. 23,
no. 1, pp. 47–51, 1979.
[8] L. Vergara-Dominguez, R. Garcia-Gomez, A. R. Figueiras-
Vidal, J. R. Casar-Corredera, and F. J. Casajus-Quiros, “Auto-
matic modelling and simulation of daily global solar radiation
series,” Solar Energy, vol. 35, no. 6, pp. 483–489, 1985.
[9] A. Balouktsis, D. Tsanakas, and G. Vachtsevanos, “Stochastic
modeling of daily global solar radiation,” International Journal
of Solar Energy, vol. 7, no. 1, pp. 1–10, 1989.
[10] P. Bendt, M. Collares-Pereira, and A. Rabl, “The frequency
distribution of daily insolation values,” Solar Energy, vol. 27,
no. 1, pp. 1–5, 1981.
[11] ELOT, Hellenic Organization of Standardization, No. 1291,
1991.
[12] J. A. Duffie and W. A. Beckman, Solar Engineering of Thermal
Processes, John Wiley & Sons, New York, NY, USA, 1991.
[13] F. Lasnier and T. G. Ang, Photovoltaic Engineering Handbook,
Adam Hilger, Bristol, UK, 1990.
[14] M. S. Imamura, P. Helm, and W. Palz, Photovoltaic System
Technology: A European Handbook, H.S. Stephens & Asso-
ciates, Bedford, UK, 1992.
[15] A. Hadj Arab, B. Ait Driss, R. Amimeur, and E. Lorenzo,
“Photovoltaic systems sizing for Algeria,” Solar Energy, vol. 54,
no. 2, pp. 99–104, 1995.