Economic Dispatch and Cost Analysis On A Power System Network Interconnected With Solar Farm
Economic Dispatch and Cost Analysis On A Power System Network Interconnected With Solar Farm
Abstract- Frequent escalation of fuel prices, concerns on environment and diminution of fossil fuel reserves have forced the
incorporation of Renewable Energy resources with the existing generation. This paper presents an Economic Dispatch model
developed for a system consisting of thermal units and Photo Voltaic (PV) plants. Since solar power is intermittent in nature,
Beta distribution function is utilized for estimating the power output of PV panels. The seasonal changes of desired location
are also included for effective analysis of the output. Case studies are carried out for various generation levels and seasonal
changes. The possible choices of installing solar plants in the test system are analysed. The proposed methodology is tested and
validated on a standard IEEE 30 bus test system..
Keywords Economic dispatch, solar irradiance, Beta distribution function, thermal generator, PV array.
balance etc. are important considerations in this problem. Ø Cost comparison by replacing some thermal units
Many researchers concentrated on the improvement in with Solar farm
general economic dispatch problem where as the research on
The rest of the paper is organised as follows: Section 2
dispatch considering renewables is limited [12].
describes the general economic dispatch problem and also
Economic dispatch using equal incremental cost was the changes to be adopted with the integration of renewable
used initially and transmission loss and penalty factor were energy in economic dispatch problem. Section 3 gives the
incorporated in that later. A review on ED methods from information on solar irradiance modelling and estimation of
1977 to 1988 are given in [19]. To evaluate the economic power output from a PV array. In section 4 integration of
dispatch problem thermal generators were used extensively. solar energy in economic dispatch and the cost optimisation
Intelligent methods like Particle swarm optimization (PSO), for various scenarios is described.
Differential Evolution (DE), Genetic Algorithm (GA),
Evolutionary Programming (EP) are used for complex 2. Economic Dispatch problem
dispatch problem considering prohibited operating zones,
valve point and quadratic cost functions. ED with cost and The main objective of economic load dispatch problem is to
emission is considered as a multi objective problem which minimise the fuel cost of the thermal generators [13]. The
requires optimization techniques to solve this.. In [20] PSO objective function can be formulated as
based ED with valve point effect is presented, which depicts
n
the effects of valve point loading of thermal generation in
generation cost. In [21] multi objective combined economic
Fi ( Pi ) = ∑ (ai Pi 2 + bi Pi + ci ) (1)
i =1
emission dispatch is converted to single objective economic
dispatch using penalty factor and ABC algorithm is used to where Fi ( Pi ) : Fuel cost of all the generators ($/hr)
solve the problem. In [22] economic dispatch considering
single solar and a wind farm is carried out on standard IEEE ai,bi,ci : Cost coefficients of ith generator.
30 bus system and SPEA algorithm is used for solving the
problem. In [23] DED is carried out with the incorporation Pi : Power generated by ith generator.
of two solar and two wind farms in addition to three thermal n : No of generators
units using a PSO algorithm and shows that the emission is
minimised with inclusion of these renewable energy sources. subjected to:
Equality constraints:
1.1. Proposed Work:
The power balance of the system is given by
m
Considering the literature survey on ED and renewable
energy sources, the solar plants are considered lossless and ∑P −P
i =1
i D − PL = 0 (2)
the climatic effects of their power outputs are not considered.
Since, the same transmission line is utilized for the where PD : Load demand (MW)
transmission of solar power also; considerable losses will
occur because of the resistance and reactance of the PL : Transmission losses (MW)
transmission lines. The solar power depends upon the Transmission losses can be represented as
irradiation which in turn depends on the climatic conditions.
n n
Considering all these factors, in this work, a fixed percentage
of loss is considered while transmitting solar power through PL = ∑∑ PB
i ij Pj
(3)
the existing transmission lines. Since the output of solar i =1 j =1
power depends in climate, the case studies are considered for where Bij : Transmission loss coefficient
different climatic conditions. In this model solar power
plants are assumed to be installed by the public utility and Inequality constraints:
the cost for installation is not considered. The power generation of all the generators has maximum and
IEEE 30 bus test system consists of 6 thermal generators, minimum limits
41 transmission lines interconnected to the 30 buses. The
thermal generator cost coefficients, Power generation limits Pi min ≤ Pi ≤ Pi max (4)
and Schematic diagram are given in appendix A and B. The
following case studies are conducted on the above said test where Pi min , Pi max are the minimum and maximum
system. generation limits
Ø Cost comparison considering solar farm output as
negative load 2.1. Coordinated Economic dispatch model considering solar
energy:
Ø Cost comparison including losses from Solar farm
Ø Cost comparison by incorporating Solar farm at The output power from a solar panel depends mainly on
desired location the irradiance. So, the power output for various irradiance
values is the be estimated which requires proper functional
1099
INTERNATIONAL JOURNAL of RENEWABLE ENERGY RESEARCH
V. Suresh and S.Suresh, Vol.5, No.4, 2015
model. The the best adopted model is beta distribution The average value of variation in solar irradiance is shown in
function. The historical data of solar irradiance is processed Fig. 1. It is observed that the irradiance is maximum in
according to seasons and then it is utilised for modelling the summer season (for the months March to June) and
beta distribution function. Using this function the output of a minimum for the months November to January.
solar panel is estimated and then the total output obtained for
the entire solar farm is calculated. This power generated by 3.2. Solar irradiance modelling
the solar farm is considered as negative demand and is
incorporated at the specific point. Then economic dispatch is A stochastic model [16] of Solar panel is constructed based
carried out using this model and the results are compared. on Beta distribution function. Beta distribution is considered
to be the most suitable model for statistical representation of
For incorporating the solar energy into the exiting
the probability density function. The Solar irradiance
generation, the power generated by PV arrays is considered
distribution of the panel is given by
as a negative load and equation (2) is updated as follows
n Γ(α + β ) (α −1) ( β −1)
PD' = PD − ∑ PiS (5)
fb (s) = s (1 − s ) ,0 ≤ s ≤ 1;α , β ≥ 0 (6)
Γ(α )Γ( β )
iS =1
n ⎛ µ (1 + µ ) ⎞ µβ
where PD' is new power demand and ∑ PiS is the sum β = (1 − µ ) ⎜ 2
− 1⎟ , α =
⎝ σ ⎠ 1− µ
iS =1
of solar power generators. where fb(s) is Beta distribution function, Γ(°) is the
gamma function and s is the random variable of solar
3. Solar Power generation model
irradiance (kw/m2), α and β are the parameters of the Beta
distribution function. µ and σ are the mean and standard
Solar power depends on meteorological conditions such as
deviation of s for the corresponding time segment.
irradiance, ambient temperature which are directly related to
geographical location [8,14]. For effective utilisation of PV 3.3. Power generation from PV array
arrays the characteristics should be desperately analysed. For
The expected output of PV is given by
the above case study, Vellore region which is located in
Tamilnadu, India is considered. P(s) = Po ( s)* fb ( s) (7)
is divided into 24 segments which refers to hours of a day. where power generation of panel at solar irradiance s is
Mean and standard deviation of solar irradiance are given by
calculated from the historical data [15]. For this study, mean
and standard deviation at 12o clock of the day for the three Po ( s ) = N * FF * Vy * I y (9)
seasons are calculated and furnished in Table 1.
where N is the total number of PV modules.
Table 1. Mean and standard deviation for three seasons
The voltage - current characteristics of a PV module for
Parameter Summer Spring Winter a given radiation level and ambient temperature are
Mean (µ) 0.9402 0.7834 0.7152 determined using the following relations
Standard Deviation (σ) 0.1088 0.1772 0.2402 VMPPT * I MPPT
FF = (10)
VOC * I SC
Vy = Voc − K v * Tcy (11)
⎛ N − 20 ⎞
Tcy = TA + s ⎜ OT ⎟ (13)
⎝ 0.8 ⎠
where FF is the fill factor, VMPPT,IMPPT are the voltage
and current maximum power point, VOC, Ioc are the open
circuit voltage and short circuit current of PV module, Kv and
Fig. 1. Monthly variation of solar irradiation (kw/m2)
1100
INTERNATIONAL JOURNAL of RENEWABLE ENERGY RESEARCH
V. Suresh and S.Suresh, Vol.5, No.4, 2015
Ki are the voltage temperature coefficient and current The discrete solar distribution of the solar panel during
temperature coefficient, TA, TCY,NOT are the ambient 12o clock of a day in winter is shown in Fig. 2 and the power
temperature, PV cell temperature and Normal operating variation during summer and spring is shown in Fig. 3.
temperature respectively. The specifications of a 220 W PV
module[17] are given in Table 2. 4. Results and Discussion
Table 2. Specifications of 220W solar panel
Proposed methodology has been simulated on standard IEEE
Parameter 30 bus system. The details of the test system are given in
Value [18]. Newton- Raphson power flow is used to evaluate the
Maximum Power point Voltage, VMPPT 28.36 V test system. A constant load model is assumed throughout
the year. It is also assumed that the geographical conditions
Maximum Power point Current, IMPPT 7.76 A analysed for considered location are similar for the test
Open circuit Voltage, VOC 36.96 V system. Based on the solar power generation model shown in
Section 2 various cases are considered for the study.
Short Circuit current, ISC 8.38 A
4.1. Scenario 1: Combined Economic dispatch with Thermal
Nominal Operating Temperature, NOT 43 °C
and Solar generation
Ambient Temperature, TA 30.76 °C
The solar farm is designed with 20000 PV modules which
Voltage Temperature Coefficient, Kv 0.1278 V/°C
comprises an installed capacity of 44 MW (20000*220W).
Current Temperature Coefficient, Ki 0.00545 A/°C The power output for each season is estimated using Beta
distribution function. The contribution in each season for
meeting the demand is plotted in Figure 4. Here four possible
combinations of utilising the solar power are discussed.
Fig. 2. Discrete Power distribution (W) of a PV module in Fig. 4. Contribution of solar energy during various seasons
Winter season
4.1.1 Case 1:
1101
INTERNATIONAL JOURNAL of RENEWABLE ENERGY RESEARCH
V. Suresh and S.Suresh, Vol.5, No.4, 2015
Table 2. Generating capacities and Fuel costs of the combined solar thermal systems
Table 4. Generating capacities and fuel costs after replacing some thermal generators by solar
1102
INTERNATIONAL JOURNAL of RENEWABLE ENERGY RESEARCH
V. Suresh and S.Suresh, Vol.5, No.4, 2015
Table 5. Cost comparison by incorporating Solar farm at specific period such that there is negligible start up cost. But
desired location practically all the generators have finite start up time. Hence,
in this case analysis is made considering the start up
Gen No Winter ($/h) Spring($/h) Summer($/h)
constraint. Here output from solar generator is sufficient to
Pg 1 - Pg S 707.41 700.10 664.00 turn off the units 4,5 and 6 for the particular time interval.
Pg 2 - Pg S 709.32 703.09 674.72 But it is not feasible to turn off the units as it imposes
additional penalty costs for the upcoming intervals if the
4.2 Scenario 2: Combined Economic Dispatch considering thermal units are to be immediately started. So for this
start up constraints situation we will run the thermal units 4,5 and 6 at minimum
MW limit and the available solar power is utilised by
The solar farm is designed with 10000 PV modules which reducing the generation of unit 2. The cost analysis is
comprises an installed capacity of 22 MW (10000*220W). presented in Table 7.
Here two case studies are done for various seasons
considering the start up constraints for the thermal From the case studies above in Scenario 2, it is evident that
generators. the scheduling of generation is an essential part. The power
generation from thermal units cannot be turned off as
4.2.1. Case 1: required depending on the solar power available. In this
This case follows the equation (5) given in section 2 and the condition, these generators are allowed to run in its minimum
economic dispatch is obtained. Here the cost comparison for generation value and the dispatch is carried out, which aviods
various seasons is presented in Table 6. the additional penalty costs. However, this may result in
increased generation cost if the costly generators are put into
4.2.2. Case 2: service.
1103
INTERNATIONAL JOURNAL of RENEWABLE ENERGY RESEARCH
V. Suresh and S.Suresh, Vol.5, No.4, 2015
given load demand results in reduced fuel cost of the loading effect. " International Journal of Electrical
convectional generators. The future work is aimed at Power & Energy Systems 73 (2015): 456-464.
investigating the real time economic dispatch considering
[14] Kayal, Partha, and C. K. Chanda. "Optimal mix of solar
multiple fuels, cost recovery of newly established plants in a
and wind distributed generations considering
scheduled period
performance improvement of electrical distribution
network." Renewable Energy 75 (2015): 173-186.
References
[15] http://rredc.nrel.gov/solar/new_data/India/nearestcell.cgi
[1] Sangroya, Deepak, and Jogendra Kumar Nayak.
[16] Teng, Jen-Hao, et al. "Optimal charging/discharging
"Development of Wind Energy in India." International
scheduling of battery storage systems for distribution
Journal of Renewable Energy Research (IJRER) 5.1
systems interconnected with sizeable PV generation
(2015): 1-13.
systems." Power Systems, IEEE Transactions on 28.2
[2] Khan, Naveed Ahmed, et al. "Combined emission (2013): 1425-1433.
economic dispatch of power system including solar photo
[17] 220 W PV Module, Solar Gate Technology Corporation,
voltaic generation." Energy Conversion and
August 2015.[Online]. Available:
Management 92 (2015): 82-91.
http://de.enfsolar.com/pv/panel-
[3] Tigas, K., et al. "Wide scale penetration of renewable datasheet/Polycrystalline/12187.
electricity in the Greek energy system in view of the
[18] Lee, K. Y., Y. M. Park, and J. L. Ortiz. "Fuel-cost
European decarbonisation targets for 2050."Renewable
minimisation for both real-and reactive-power
and Sustainable Energy Reviews 42 (2015): 158-169.
dispatches." Generation, Transmission and Distribution,
[4] Kalantari, Amir, and Francisco D. Galiana. "Generalized IEE Proceedings C. Vol. 131. No. 3. IET, 1984.
Sigma approach to unit commitment with uncertain wind
[19] Chowdhury, Badrul H., and Saifur Rahman. "A review of
power generation." International Journal of Electrical
recent advances in economic dispatch." Institute of
Power & Energy Systems 65 (2015): 367-374.
Electrical and Electronics Engineers, 1990.
[5] http://mnre.gov.in/file-manager/annual-report/2014-
[20] Lim, Shi Yao, Mohammad Montakhab, and Hassan
2015/EN/Chapter%204/chapter_4.htm
Nouri. "Economic dispatch of power system using
[6] Ela, Erik, and Mark O'Malley. "Studying the variability particle swarm optimization with constriction
and uncertainty impacts of variable generation at multiple factor."International Journal of Innovations in Energy
timescales." Power Systems, IEEE Transactions on 27.3 Systems and Power 4.2 (2009): 29-34.
(2012): 1324-1333.
[21] Sonmez, Y. "Multi-objective environmental/economic
[7] Fan, Miao, et al. "Probabilistic power flow analysis with dispatch solution with penalty factor using Artificial Bee
generation dispatch including photovoltaic Colony algorithm." Scientific Research and Essays 6.13
resources." Power Systems, IEEE Transactions on 28.2 (2011): 2824-2831.
(2013): 1797-1805.
[22] Brini, Saoussen, Hsan Hadj Abdallah, and Abderrazak
[8] Bou-Rabee, Mohammed, and Shaharin Sulaiman. "On Ouali. "Economic dispatch for power system included
Seasonal Variation of Solar Irradiation in wind and solar thermal energy." Leonardo Journal of
Kuwait." International Journal of Renewable Energy Sciences 14 (2009): 204-220.
Research (IJRER) 5.2 (2015): 367-372.
[23] ElDesouky, AZZA Ahmed. "Security and stochastic
[9] Atwa, Y. M., et al. "Optimal renewable resources mix for economic dispatch of power system including wind and
distribution system energy loss minimization." Power solar resources with environmental
Systems, IEEE Transactions on 25.1 (2010): 360-370. consideration." International Journal of Renewable
Energy Research (IJRER)3.4 (2013): 951-958.
[10] Salameh, Ziyad M., Bogdan S. Borowy, and Atia RA
Amin. "Photovoltaic module-site matching based on the Appendix A: Cost coefficients and generation limits for an
capacity factors." Energy Conversion, IEEE Transactions IEEE 30 bus system
on 10.2 (1995): 326-332. Cost Coefficients Generation
Generator
($/hr) Limits(MW)
[11] Chowdhury, Badrul H., and Saifur Rahman. "A review of Bus No
ai bi ci Pgmin Pgmax
recent advances in economic dispatch." Institute of
1 0.00375 2 0 50 200
Electrical and Electronics Engineers, 1990.
2 0.0175 1.75 0 20 80
[12] Wang, Jianhui, et al. "Impact of wind power forecasting 5 0.0625 1 0 15 50
on unit commitment and dispatch." Proc.8th Int. 8 0.0083 3.25 0 10 55
Workshop Large-Scale Integration of Wind Power into 11 0.025 3 0 10 30
Power Systems. 2009. 13 0.025 3 0 12 40
[13] Banerjee, Sumit, Deblina Maity, and Chandan Kumar
Chanda. "Teaching learning based optimization for
economic load dispatch problem considering valve point
1104
INTERNATIONAL JOURNAL of RENEWABLE ENERGY RESEARCH
V. Suresh and S.Suresh, Vol.5, No.4, 2015
1105