 Open Access
 Authors : Ashutosh Kumar Singh , Aditya Kumar Rawat , Dinesh Kumar Nishad , Ranitesh Gupta
 Paper ID : IJERTV10IS030230
 Volume & Issue : Volume 10, Issue 03 (March 2021)
 Published (First Online): 05042021
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Optimal Dispatch of Power System with Integrated System Constraints
Aditya Kumar Rawat1
Electrical Engineering Dept.
DSMNRU
Lucknow, India
Ashutosh Kumar Singp
Electrical Engineering Dept.
DSMNRU
Lucknow, India
Dinesh Kumar Nishad3
Electrical Engineering Dept.
DSMNRU
Lucknow, India
Ranitesh Gupta4
Electrical Engineering Dept.
DSMNRU
Lucknow, India
AbstractThis paper presents an economic load dispatch (ELD) model, incorporating wind power with thermal generator units.ELD is a fundamental issue that simulates optimal scheduling and dispatch of available generation by minimizing overall cost subjected to physical and functional constraints. Uncertain nature of wind results in the calculation of overestimation and underestimation. Forecasting of wind power is an issue in economic load dispatch, that is why, the wind turbine output prediction, is being done by scenario analysis method. In this, a mathematical model is used which quantifies the dispatch flexibility of thermal generation by integrating the renewable energy source (wind energy). To investigate system operations, one of the modern optimizations technique named as Particle Swarm Optimization (PSO) is being used. The results are demonstrated on the basis of calculation of overestimation and underestimation cost factor separately. The optimization problem is numerically solved for a scenario involving six thermal generators units and a wind farm.
KeywordsEconomic load dispatch, Scenario analysis, Wind power, Penalty cost, Reserve cost, Particle swarm optimization.

INTRODUCTION
In maximum costeffective power system services, ELD plays a pivotal role. The goal is to reduce energy costs by meeting all physical constraints together with the energy source [1]. The regular increase in fuel costs, environment concerns and the depletion of fossil fuels have prompted the inclusion of current generation in the renewable energy resources. Therefore, the integrated energy planning has become a requirement of any country for sustainable energy sector development since economic load dispatch models helps to enhance the stochastic availability of wind power. As the production of fossil fuels reached its peak in 2015 [2], it is now declining because, when consumed, fossil fuels cant be extracted and it takes a long geographical time to recover, keeping in view of our future demand of electricity and to
meet them alternative of fossil fuels are mandatory as wind. Wind energy is the fastgrowing renewable energy source because of its abundance in nature and maturity of its technology. The main driver of this transformation is the advent of wind as the strong price winner against fossil fuels and nuclear energy as a clean source of energy. Therefore ,wind
energy impacts are generally seen as more environmentally friendly than thermal energy source impacts.
The first important issue of wind energy is that the nature of wind, which is uncertain [3]. Because, it is highly dependent on weather and geographic situations. Because of this discontinuous nature, the power generated with wind farm is uncertain and difficult to dispatch. For this reason, renewable energy source, such as wind, is overestimated and underestimated making a mathematical model, which also quantifies the direct cost, Penalty cost and Reserve cost.
As, India rank 4thin global wind power installed capacity index [4].In present scenario, our dependency is growing towards power, and India will require 34 times as much energy than the total energy consumed today, in order to meet the energy requirement of such a fastgrowing economy renewable energy is one and only option because approximately 63% of total power is generated through fossil fuels, which pollutes the environment, but renewable energy accounts for approximately 33% of Indias primary consumption[5] and as of now reducing global 2 emissions is one of todays major challenges [6]. However, it seems more than uncertain in many parts of the world whether the objectives of the Paris Climate Agreement (PCA) will be achieved [7].
There are many options for power generation as Solar Energy, Super capacitors, Superconducting magnetic energy storage, Wind energy, Wind energy storage system, Hydro pump stations, Compressed air system and many more [8]. As, available technology is allowing the renewable industry to achieve important milestones. In general, the future of renewable energy, and especially wind energy, is very promising since electricity generated by wind turbines does not pollute the water we drink or the air we breathe, so wind power means less smog, less acid rain, and less emission of greenhouse gases. Because of this reason, Wind energy is truly well to fulfill the objectives of Paris Climate Agreement (PCA) and providing a motivation to move towards renewable energy sources.
The collective cost function of conventional thermal generators, with the consideration of renewable energy source, is considered as an economic factor. To achieve this economic factor, a dayahead scheduling of wind is performed, for 24 hours. Moreover, an economic optimization is obtained by
newly developed algorithm named as Particle Swarm Optimization (PSO) introduced by Kennedy and Eberhart [9]. The aim of this paper is to combine wind generators with traditional thermal generators to solve the problem of economic dispatch using a mathematical model.
Here, the objective of economic load dispatch (ELD) is to
,, penalty cost function;
,, reserve cost function; D total demand.
=
=
Taking a clear view of objective function which specifies each term as the 1st the term 1 ( ) is the traditional sum of fuel
costs of the thermal generator units. 2ndterm 1 ,() is
minimize the operating cost of generation from integration of classic thermal generators and wind farm including the
the direct cost of wind generator.
=
The 3rdterm
underestimation and overestimation, having certain physical constraints. The mathematical equation [10] which specifies the
1 ,,(, ) accounts for penalty cost of available
=
=
=
=
wind power (underestimation) and the last term
model considered and taken as objective function is as follows:
1 ,,( , )
signifies reserve cost because of
=1 () + 1 , + 1 ,,(, ) +
overestimation of available wind power.
For classic thermal generators, a quadratic function is taken into
=
=
=
1 ,,( ,)
=
(1)
account [11] and this function is applied for most of cases and defined as:
() = ()2 + + ; (6)
Wind power Output forecasting
Dynamic

PROPOSED MODEL
Unit commitment of thermal Generators in the system
Where is the generated (input) power from conventional generators, and , , are the operating cost coefficient of
conventional thermal units, which are obtained from inputoutput of generators and relies on type of fuel used as per objective function equation.
For wind generators, a linear cost function is taken into account in which the wind power generation cost , () may not exist if the operator owns the wind Powered generators, because the output of the wind powered generators is restricted by an upper and lower limit, decided by the system operator based on the
constraints( e.g.
the power otput constraints
Of thermal units)
Power dispatch
annexure for optimal operation of system, but it could be assumed as a payback cost or maintenance cost and defined as:
, () = (7)
Where,
scheduled wind power from the ith wind generator;
Making power generation schedules according to the dispatch model
direct cost coefficient for the ith wind generator.
Load Center
Overestimation represents the deficit of wind power because, if a certain amount of wind power is a considered and that power is not found at the considered time, power must be purchased from an alternate source and it is defined as:
,,( ,) = ,( , );
Fig.1. Proposed model
=
( )
() ; (8)
The entire model is based on the summation of generated power
Where,
, 0
through conventional thermal generators and wind farm, should be equal to losses as follows:
, reserve cost( overestimation)coefficient for the ith wind powered generator;
=1
+
+
=1
=
(2)
(w) wind power pdf.
And underestimation represents the surplus of wind power
because, if available wind power is actually more than what was
And = + (3)
Subjected to , , (4)
0 , (5)
Where,
considered, that power will be wasted, and it is reasonable for the system operator to pay money to the wind power producer for the waste of available wind power and it is defined as:
,,(, ) = ,(, );

number of conventional thermal generators;
=
,
( )
( )
(); (9)

number of wind generators;
power form the conventional generators;
Where,
,
scheduled wind power from wind generator;
, penalty cost (underestimation) coefficient for
,
available wind power from the ith wind generator;
the wind generator.
, rated wind power from the ith wind generator;
cost function for the ith thermal generator;
, cost function for the ith wind generator;


SCENARIO ANALYSIS METHOD FOR
WIND FORECASTING
The generation of wind power varies from traditional thermal generation due to the stochastic nature of wind [12]. Wind power forecasting therefore plays a key role in resolving the problems of supply balancing and demand. Precise wind power projections minimize the need for additional energy balance and reserve wind power integration capability. Tools for wind power prediction allow better transmission, planning and thermal generators device involvement [13]. There are three steps in the forecasting of wind power:

Wind speed of a model

Wind power measurement

Prediction of outputs.

Shortterm Forecast: From 1hr to several hrs ahead.

Medium term Forecast: From several hrs to 1 week ahead.

Longterm Forecast: From 1 week to 1 year ahead.
In this paper wind power forecasting done by using scenario analysis method.

Shortterm Forecast [15], which focuses on wind hourly data. So, the wind power data of one year is picked from NREL.

Since, nature of wind is uncertain, Due to which the curve of forecasted wind would be nonlinear.

Normal Distribution
The normal distribution is an extraordinary symmetrical distribution approach used to make comparisons among
Power
Cut in speed
Rated speed
Cut out Speed
different kinds of statistical decisions having bell curve [16]. This curve says that mainly results lie close to centre as we move away from centre its frequency decreases. Sometimes is also known as natural distribution.
= 222 ; (12)
Where,
e exponential constant;
mean or location parameter;
standard deviation ;
Fig.2. Variation of output wind power with wind speed
Power equation of wind turbine is described as below:
3 3
3 3
() ()
2 variance or scale parameter.
Bell curve relies on two parameters and which controls the xaxis and yaxis respectively. In this, the curve is completely
() = { () ; < () < } ; (10)
33
Where,
< () <
symmetrical about . This expresses that is equal to expectation. It can also be used to determine the percentage of
data present above or below of assumed data.
wind speed;
rated wind speed;
cutin wind speed;
rated wind power.
Rated wind power can be defined as:
TABLE 1. Probability of wind power.
= 1
3 ; (11)
S. No. 
Forecast error 
Rate of wind power 
Probability 
1. 
3 
70% 
0.0063 
2. 
2 
80% 
0.0607 
3. 
90% 
0.2416 

4. 
0 
0 
0.3834 
5. 
+ 
110% 
0.2416 
6. 
+2 
120% 
0.0607 
7. 
+3 
130% 
0.0063 
S. No. 
Forecast error 
Rate of wind power 
Probability 
1. 
3 
70% 
0.0063 
2. 
2 
80% 
0.0607 
3. 
90% 
0.2416 

4. 
0 
0 
0.3834 
5. 
+ 
110% 
0.2416 
6. 
+2 
120% 
0.0607 
7. 
+3 
130% 
0.0063 
2
Where,
rated wind power;
power coefficient;
a air density;
efficiency of turbine generator;
area swept by turbine blade;
rated speed of wind.
One of the ultimate objectives of a model of prediction of wind
power is to estimate the output of the wind as early and as accurately as possible. Because overall detailed prediction of wind power reduces wind insecurity's financial and technological danger power generation.
There are various methods for forecasting the production of wind power. These approaches are categorized as centered on time scales and according to numerous methodologies available in the literature [14]. They are exactly divided in four categories:

Ultrashortterm forecast: Form a few minutes to 1hr ahead.
Fig.3. Wind power curves under different scenarios.

METHODOLOGY
In this objective function, an optimization tool, which determines location and velocity of particles in multi dimensional search space, and whose each particle adjust itself making best use of its own experience and that of its neighbor named as PSO, and is being used [17]. In PSO conventional terminology parameters are defined as:
Let and denotes particle location and velocity in multi dimensional search space respectively. Therefore, particle location is defined as:
= (1, 2, 3, . , ); (13)
= 1,2,3,4, . , ; (14)
Where,
m no of units in multi dimensional space.
Previous portion of particle is taken into account and represented as:
= ( 1, 2, 3, . , ); (15) Here the best position of particle is taken as ( ) and in
Computational parameters used in PSO are given below

No. of particles = 25
No. of iteration / ( ) = 1000

Acceleration constant, c1 =2 & c2 = 2
Weighting factor, = 0.4 & = 0.9

Change in velocity for each individual,
= 0.5 ; (21)
Vpdmin = 0.5Pdmin; (22)
Start
Read data from text file; demand, available and predicted values of wind power
Calculate the difference between renewable power and total demand
Initialize the parameters of PSO (c1, c2, w, x, v, and itermax)
particular dimension () the index of best particle among all
the particles is denoted by (), and it is also termed as global best. The velocity of particle is denoted as:
= (1, 2, 3, . , ); (16) Now, modified velocity and location of each particle can be demonstrated by using following formulas:
+1 = + 1 1 ( ) + 2 2 ( ); (17)
+1 = + +1; (18)
= 1,2,3, . , ; = 1,2,3 . , ; (19)
Where,
n population size;
m number of units;
t pointer of iterations;
w inertia weight factor; c1, c2 acceleration constant;
U1, U2 uniform random values in the range [0, 1];
location of particle i at iteration t;
velocity of particle i at iteration t.
The acceleration constants 1 and 2 on basis of past experience were often set to be 2.0 which () denotes the weighting of the stochastic acceleration terms that pull each particle toward and positions.
Proper selection of inertia weight factor () is helpful to
increase the rate of convergence of the standard PSO algorithm and suitable to find the optimal solution. As originally developed inertia weight factor(), very often decrease from
0.9 to 0.4 linearly during a run. In general, the inertia weight () proposed in velocity equation is represented by following equation:
Calculate transmission losses
Start Iteration Iter = 0
Calculate the objective function value
Update the gbest and pbest values
Update particle velocity and location Using (14) and (15)
Iteration, Iter = Iter+1
Is this stopping Criterion satisfied NO
YES
gbest value is the solution of ELD with renewable energy problem
Print Result
Stop
=
Where,
() ; (20)
Fig.4. Flow chart of PSO
maximum number of iterations; iter current number of iterations.
Cost Comparison
1500
1000
500
0
1 2 3 4 5 6 7 8 9 10 11 12
Time internal (hour)
Cost Comparison
1500
1000
500
0
1 2 3 4 5 6 7 8 9 10 11 12
Time internal (hour)
TABLE 2. Generating unit power and coefficients
0.000140 0.000017 0.000015 0.000019 0.000026 0.000022
0.000017 0.000060 0.000013 0.000016 0.000015 0.000020
0.000015 0.000013 0.000065 0.000017 0.000024 0.000019
=
0.000019 0.000016 0.000017 0.000071 0.000030 0.000025
0.000026 0.000015 0.000024 0.000030 0.000069 0.000032
[0.000022 0.000020 0.000019 0.000025 0.000032 0.000085]Unit
(M
W)
(MW
)
($/)
($/MW)
($)
1
10
125
0.15240
38.53973
756.7988
6
2
10
150
0.10587
46.15916
451.3251
3
3
35
225
0.02803
40.39655
1049.997
7
4
35
210
0.03546
38.30553
1243.531
1
5
130
235
0.02111
36.32782
1658.559
6
6
125
315
0.01799
38.27041
1356.659
2
TABLE 3: Best power output for 6 thermal generators and single integrated wind farm
Unit Output
Conventional
PSO
P1(MW)
23.90
23.84
P2(MW)
10.00
10.00
P3(MW)
95.63
95.57
P4(MW)
100.70
100.52
P5(MW)
202.82
202.78
P6(MW)
182.02
181.52
1()
46
45.7
Total power output(MW)
34013.07
34011.1
Total generation cost(MW)
32096.58
32094.69
Power losses(MW)
15.07
14.23


RESULT
In this proposed work, an economic load dispatch model is developed with the help of integration of wind farm. To investigate system operations, one of the modern optimization technique named as Particle Swarm Optimization (PSO) has been used. The results are demonstrated on the basis of calculation of overestimation and underestimation cost factor separately. The optimization problem is numerically solved for a scenario involving six thermal generators units and wind farm.
Cost Estimation Curve
Fig.5. Curve between Cost ($) and Time (hr)
Cost ($/hr)
Cost ($/hr)
Fig.6. Cost comparison without penalty factor and with penalty factor respectively

CONCLUSION

In this paper, an economic load dispatch model is developed with the help of integration of wind farm. Forecasting of wind power has been done by Scenario Analysis method and uncertain nature of wind is responsible for the calculation of overestimation and underestimation factors and model results are also demonstrated on the basis of calculation of overestimation and underestimation cost factor separately. The problem is numerically solved by using one of modern optimization algorithm (PSO), for a scenario involving six thermal generators units and wind farm. The results are also found to be most reliable and environmental friendlier.
REFERENCES

Talaq, J. H., Ferial ElHawary, and M. E. ElHawary. "A summary of environmental/economic dispatch algorithms." IEEE Transactions on Power Systems 9.3 (1994): 15081516.

National Renewable Energy Laboratory(NREL). http://www.nrel.gov/

Liu, Xian, and Wilsun Xu. "Economic load dispatch constrained by wind power availability: A hereandnow approach." IEEE Transactions on sustainable energy 1.1 (2010): 29.

GWEC (Global Wind Energy Council). "Global Wind Statistics 2017." Global Wind Energy Council Rep. (2018): 4.

Kumar, Ashwani, et al. "Renewable energy in India: current status and future potentials." Renewable and sustainable energy reviews 14.8 (2010): 24342442.

Poulter, Benjamin, et al. "Contribution of semiarid ecosystems to interannual variability of the global carbon cycle." Nature 509.7502 (2014): 600603.

Salawitch,Ross J., and Timothy P. Canty. Paris climate agreement: Beacon of hope. Springer International Publishing, 2017.

Luo, Fengji, et al. "Coordinated operational planning for wind farm with battery energy storage system." IEEE Transactions on Sustainable Energy 6.1 (2015): 253262.

Park, JongBae, et al. "An improved particle swarm optimization for nonconvex economic dispatch problems." IEEE Transactions on Power Systems 25.1 (2009): 156166.

Hetzer, John, C. Yu David, and Kalu Bhattarai. "An economic dispatch model incorporating wind power." IEEE Transactions on energy conversion 23.2 (2008): 603611.

J.J Grainger and W.D. Stevension, Jr., Power System Analysis, New York: McGrawHill, 1994.

Hetzer, John, C. Yu David, and Kalu Bhattarai. "An economic dispatch model incorporating wind power." IEEE Transactions on energy conversion 23.2 (2008): 603611.

METHAPRAYOON, Kittipong, et al. An integration of ANN wind power estimation into unit commitment considering the forecasting uncertainty. IEEE Transactions on Industry Applications, 2007, 43.6: 14411448.

MA, XiYuan; SUN, YuanZhang; FANG, HuaLiang. Scenario generation of wind power based on statistical uncertainty and variability. IEEE Transactions on Sustainable Energy, 2013, 4.4: 894 904.

ZHOU, Yun; YAN, Zheng; LI, Naihu. A novel state of charge feedback strategy in wind power smoothing based on shortterm forecast and scenario analysis. IEEE Transactions on Sustainable Energy, 2016, 8.2: 870879.

XU, Jian, et al. Stochastic optimal scheduling based on scenario analysis for wind farms. IEEE Transactions on Sustainable Energy, 2017, 8.4: 15481559.

Gaing, ZweLee. "Particle swarm optimization to solving the economic dispatch considering the generator constraints." IEEE Transactions on power systems 18.3 (2003): 11871195