UGC - 7 - Leela Kumari
UGC - 7 - Leela Kumari
org (ISSN-2349-5162)
Abstract: Economic Load dispatch is one of the challenging problems in power system planning, operation and
control. Single area economic load dispatch problem of electric power system is confined to electric utilities of
similar area. In the proposed research, the authors has presented the mathematical formulation of single area
economic load dispatch problem incorporating renewable energy sources, plug-in electric vehicles (PEVs) and
battery electric vehicles (BEVs). This mathematical formulation of single area dynamic dispatch will be helpful to
the researchers, who are working in single area economic load dispatch problems with electric vehicles (EVs) and
Renewable energy sources (RES) as one of the their research objectives.
Keywords: BEVs, Economic Load Dispatch, PEVs, Renewable Energy Sources.
1. INTRODUCTION
In the expansion of different sectors of economy electric power performs major role in the survival. Actually, every
economy completely relays on power as a fundamental feeding. This is directed to guide in order to multiply more
number of power stations and their capabilities and improved its number in the power lines which conjugates the
generating stations to the load centers. To enhance reliability and economy, interconnections between systems are
also been increased. While depending on the power transmitted, the transmission voltage should provide flexibility
in system operation by long-term system requirement. The main aim of the electric power companies is to the
provide the high-quality, reliable electric power supply to the load center /distribution station at most economical
way, subject to the fulfillment of all the power generating units constraints.
The Economic Load Dispatch (ELD) problem may be stated very concisely as determining the loadings of the
various thermal generators available in the power system for which the total cost of generation is minimized,
maintaining the power balance equality constraint and the generation capacity inequality constraints. The
complete mathematical formulation of Dynamic Economic dispatch has been divided into the sub-sections: (i)
Classical Economic Load Dispatch (ii) Economic Load Dispatch with Renewable Energy Sources (iii)
Economic Load Dispatch with BEVs and PEVs (iv) Economic Load dispatch with BEVs/PEVs and Renewable
Energy Sources.
JETIRDW06220 Journal of Emerging Technologies and Innovative Research (JETIR) www.jetir.org 1367
© 2019 JETIR January 2019, Volume 6, Issue 1 www.jetir.org (ISSN-2349-5162)
ABBREVIATIONS
3. MATHEMATICAL FORMULATION
The foremost objective of the single area dynamic dispatch is to minimize the total fuel cost of the power generating
units subject to the fulfillments of different constraints. The overall objective function of the single area dynamic
dispatch problem can be categorized to the following sub-sections [1-6]:
The mathematical formulation of conventional single area dynamic dispatch for one hour can be represented as:
NG
F ( PG ) [an ( PnG ) 2 bn PnG cn )]
n 1 (1)
The dispatch of power generating units for ‘H’ Hours may be represented as:
H
NG
F ( P ) [an ( PnG ) 2 bn PnG cn )]
G
(2)
h 1 n 1
This eqn.(2) represented the exact mathematical formulation for Dynamic Dispatch. The hour ‘h’ may be varied for 1
to H-th Hour for time varying load demand.
JETIRDW06220 Journal of Emerging Technologies and Innovative Research (JETIR) www.jetir.org 1368
© 2019 JETIR January 2019, Volume 6, Issue 1 www.jetir.org (ISSN-2349-5162)
To establish the output of the power of an online generating units the economic load dispatch problem aims to
congregate the system load at least price whilst fulfilling the system constraints. So as to attain correct dispatch
outcomes, a cubical function is used for modeling the unit cost.
NG
F ( PG ) [an ( PnG )3 bn ( PnG ) 2 cn PnG d n )] (3a)
n 1
The Cubical ELD with Valve Point effect can be represented as:
NG
F ( PG ) [an ( PnG )3 bn ( PnG ) 2 cn PnG d n )] n sin( n ( PnG(min) PnG ) ) (3b)
n 1
The heat and power economic load dispatch problem of a system is to resolve the unit heat and power production. The
mathematical formulation for heat and power ELD may be described as:
NG
FPower ( PnG ) cn bn PnG an ( PnG )
2
(4a)
n 1
NG
FHeat ( PnG ) g n hn PnG qn ( PnG )
2
(4b)
n 1
NG
FOverall ( PnG ) cn bn PnG an ( PnG ) g n PnG hn ( PnG ) qn ( PnG )
2 2 2
(4c)
n 1
and the objective function for heat and power ELD considering valve point loading effects can be reframed as:
NG
FPower ( PnG ) cn bn PnG an ( PnG ) n sin( n ( PnG(min) PnG ) )
2
(5a)
n 1
NG
FHeat ( PnG ) g n hn PnG qn ( PnG ) n sin( n ( PnG(min) PnG ) )
2
(5b)
n 1
NG
FOverall ( PnG ) cn bn PnG an ( PnG ) g n PnG hn ( PnG ) qn ( PnG ) n sin( n ( PnG(min) PnG ) )
2 2 2
n 1
(6)
All the above mentioned objective functions are subjected to the following equality and inequality constraints:
NG
P
n 1
n
G
P Demand P Loss (7)
JETIRDW06220 Journal of Emerging Technologies and Innovative Research (JETIR) www.jetir.org 1369
© 2019 JETIR January 2019, Volume 6, Issue 1 www.jetir.org (ISSN-2349-5162)
NG NG
P Loss P
n 1 m 1
G G
n Bnm Pm (8)
if Bi 0 and B00 matrices for loss coefficients are given, then the above equation can be modified as:
NG
P Loss PnG Bnm PmG P
n 1
n
G
Bi 0 B00 (9)
The real power generation of each generator is to be controlled inside its particular upper and lower operating limits.
max[ PnG(max) ,(URn PnG )] PnG min[ PnG(max) ,( PnGo DRn )] n 1, 2,3,..., NG (14)
3.2 Single Area Dynamic Load Dispatch Considering Renewable Energy Sources
The overall objective function for the Dynamic dispatch considering renewable energy sources may be described
as:
NG
F ( PG ) [an ( PnG ) 2 bn PnG cn )] (16a)
n 1
The dispatch of power generating units for ‘H’ Hours may be represented as:
H
NG
F ( P G ) [an ( PnG ) 2 bn PnG cn )] (16b)
h 1 n 1
JETIRDW06220 Journal of Emerging Technologies and Innovative Research (JETIR) www.jetir.org 1370
© 2019 JETIR January 2019, Volume 6, Issue 1 www.jetir.org (ISSN-2349-5162)
This eqn.(16b) represented the exact mathematical formulation for Dynamic Dispatch. The hour ‘h’ may be varied for
1 to H-th Hour for time varying load demand.
To establish the output of the power of an online generating units the economic load dispatch problem aims to
congregate the system load at least price whilst fulfilling the system constraints. So as to attain correct dispatch
outcomes, a cubical function is used for modeling the unit cost [2-4].
NG
F ( PG ) [an ( PnG )3 bn ( PnG ) 2 cn PnG d n )] (17a)
n 1
The Cubical ELD with Valve Point effect can be represented as:
NG
F ( PG ) [an ( PnG )3 bn ( PnG ) 2 cn PnG d n )] n sin( n ( PnG(min) PnG ) ) (17b)
n 1
NG
FPower ( PnG ) cn bn PnG an ( PnG )
2
(18a)
n 1
NG
FHeat ( PnG ) g n hn PnG qn ( PnG )
2
(18b)
n 1
NG
FOverall ( PnG ) cn bn PnG an ( PnG ) g n PnG hn ( PnG ) qn ( PnG )
2 2 2
(18c)
n 1
and the objective function for heat and power ELD considering valve point loading effects can be reframed as:
NG
FPower ( PnG ) cn bn PnG an ( PnG ) n sin( n ( PnG(min) PnG ) )
2
(19a)
n 1
NG
FHeat ( PnG ) g n hn PnG qn ( PnG ) n sin( n ( PnG(min) PnG ) )
2
(19b)
n 1
NG
FOverall ( PnG ) cn bn PnG an ( PnG ) g n PnG hn ( PnG ) qn ( PnG ) n sin( n ( PnG(min) PnG ) )
2 2 2
n 1
(20)
All the above mentioned objective functions are subjected to the following equality and inequality constraints:
NG
P
n 1
n
G
P Demand P Loss (21)
NG
P
n 1
G
n P
Re newable
P Demand P Loss (22)
JETIRDW06220 Journal of Emerging Technologies and Innovative Research (JETIR) www.jetir.org 1371
© 2019 JETIR January 2019, Volume 6, Issue 1 www.jetir.org (ISSN-2349-5162)
Where,
NG NG
P Loss P
n 1 m 1
G G
n Bnm Pm (23)
if Bi 0 and B00 matrices for loss coefficients are given, then the above equation can be modified as:
NG
P Loss PnG Bnm PmG P
n 1
n
G
Bi 0 B00 (24)
max[ PnG(max) ,(URn PnG )] PnG min[ PnG(max) ,( PnGo DRn )] n 1, 2,3,..., NG (29)
The overall objective function for the Dynamic dispatch considering PEVs and BEVs may be represented as:
NG
F ( PG ) [an ( PnG ) 2 bn PnG cn )] (31a)
n 1
The dispatch of power generating units for ‘H’ Hours may be represented as:
H
NG
F ( P G ) [an ( PnG ) 2 bn PnG cn )] (31b)
h 1 n 1
JETIRDW06220 Journal of Emerging Technologies and Innovative Research (JETIR) www.jetir.org 1372
© 2019 JETIR January 2019, Volume 6, Issue 1 www.jetir.org (ISSN-2349-5162)
This eqn.(31b) represented the exact mathematical formulation for Dynamic Dispatch. The hour ‘h’ may be varied for
1 to H-th Hour for time varying load demand.
To establish the output of the power of an online generating units the economic load dispatch problem aims to
congregate the system load at least price whilst fulfilling the system constraints. So as to attain correct dispatch
outcomes, a cubical function is used for modeling the unit cost.
NG
F ( PG ) [an ( PnG )3 bn ( PnG ) 2 cn PnG d n )] (32a)
n 1
The Cubical ELD with Valve Point effect can be represented as:
NG
F ( PG ) [an ( PnG )3 bn ( PnG ) 2 cn PnG d n )] n sin( n ( PnG(min) PnG ) ) (32b)
n 1
NG
FPower ( PnG ) cn bn PnG an ( PnG )
2
(33a)
n 1
NG
FHeat ( PnG ) g n hn PnG qn ( PnG )
2
(33b)
n 1
NG
FOverall ( PnG ) cn bn PnG an ( PnG ) g n PnG hn ( PnG ) qn ( PnG )
2 2 2
(33c)
n 1
and the objective function for heat and power ELD considering valve point loading effects can be reframed as:
NG
FPower ( PnG ) cn bn PnG an ( PnG ) n sin( n ( PnG(min) PnG ) )
2
(34a)
n 1
NG
FHeat ( PnG ) g n hn PnG qn ( PnG ) n sin( n ( PnG(min) PnG ) )
2
(34b)
n 1
NG
FOverall ( PnG ) cn bn PnG an ( PnG ) g n PnG hn ( PnG ) qn ( PnG ) n sin( n ( PnG(min) PnG ) )
2 2 2
n 1
(35)
The overall fuel cost function considering BEVs and PEVs has to be reduced within the following constraints:
NG
P
n 1
n
G
P Demand P Loss (36)
JETIRDW06220 Journal of Emerging Technologies and Innovative Research (JETIR) www.jetir.org 1373
© 2019 JETIR January 2019, Volume 6, Issue 1 www.jetir.org (ISSN-2349-5162)
NG NPEVs NBEVs
n 1
PnG P Demand P Loss n 1
P PEVs P
n 1
BEVs
(37)
NG NPEVs NBEVs
n 1
PnG P Demand P Loss n 1
P PEVs P
n 1
BEVs
(38)
Where,
NG NG
P Loss P
n 1 m 1
G G
n Bnm Pm (39)
if Bi 0 and B00 matrices for loss coefficients are given, then the above equation can be modified as:
NG
P Loss PnG Bnm PmG P
n 1
n
G
Bi 0 B00 (40)
The real power generation of each generator is to be controlled inside its particular upper and lower operating limits.
max[ PnG(max) ,(URn PnG )] PnG min[ PnG(max) ,( PnGo DRn )] n 1, 2,3,..., NG (45)
JETIRDW06220 Journal of Emerging Technologies and Innovative Research (JETIR) www.jetir.org 1374
© 2019 JETIR January 2019, Volume 6, Issue 1 www.jetir.org (ISSN-2349-5162)
3.4 Single Area Dynamic Dispatch Considering PEVs and BEVs and Renewable Energy Sources
The overall objective function for the Dynamic dispatch considering PEVs, BEVs and renewable energy sources
may be represented as:
NG
F ( P ) [an ( PnG ) 2 bn PnG cn )]
G
(47a)
n 1
The dispatch of power generating units for ‘H’ Hours may be represented as:
H
NG
F ( P G ) [an ( PnG ) 2 bn PnG cn )] (47b)
h 1 n 1
To establish the output of the power of an online generating units the economic load dispatch problem aims to
congregate the system load at least price whilst fulfilling the system constraints. So as to attain correct dispatch
outcomes, a cubical function is used for modeling the unit cost.
NG
F ( PG ) [an ( PnG )3 bn ( PnG ) 2 cn PnG d n )] (48a)
n 1
The Cubical ELD with Valve Point effect can be represented as:
NG
F ( P ) [an ( PnG )3 bn ( PnG ) 2 cn PnG d n )] n sin( n ( PnG(min) PnG ) )
G
(48b)
n 1
3.4.2 Heat and Power ELD with PEVs, BEVs and RES
NG
FPower ( PnG ) cn bn PnG an ( PnG )
2
(49a)
n 1
NG
FHeat ( PnG ) g n hn PnG qn ( PnG )
2
(49b)
n 1
NG
FOverall ( P ) cn bn PnG an ( PnG ) g n PnG hn ( PnG ) qn ( PnG )
2 2 2
G
n (49c)
n 1
and the objective function for heat and power ELD considering valve point loading effects can be reframed as [8-11]:
NG
FPower ( PnG ) cn bn PnG an ( PnG ) n sin( n ( PnG(min) PnG ) )
2
(50a)
n 1
NG
FHeat ( PnG ) g n hn PnG qn ( PnG ) n sin( n ( PnG(min) PnG ) )
2
(50b)
n 1
NG
FOverall ( PnG ) cn bn PnG an ( PnG ) g n PnG hn ( PnG ) qn ( PnG ) n sin( n ( PnG(min) PnG ) )
2 2 2
n 1
(51)
JETIRDW06220 Journal of Emerging Technologies and Innovative Research (JETIR) www.jetir.org 1375
© 2019 JETIR January 2019, Volume 6, Issue 1 www.jetir.org (ISSN-2349-5162)
The overall fuel cost function has to be minimized within the following constraints:
NG
P
n 1
n
G
P Demand P Loss (52)
NG
P
n 1
G
n P
Re newable
P Demand P Loss (53)
The power balance constraints during charging phase for PEVs and BEVs can be mathematical described as per the
following eqn.
NG NPEVs NBEVs
n 1
PnG P Re newable P Demand P Loss n 1
P PEVs P
n 1
BEVs
(54)
The power balance constraints during discharging phase for PEVs and BEVs can be mathematical described as per
the following eqn.(2).
NG NPEVs NBEVs
n 1
PnG P Re newable P Demand P Loss
n 1
P PEVs P
n 1
BEVs
(55)
Where,
NG NG
P Loss P
n 1 m 1
G G
n Bnm Pm (56)
if Bi 0 and B00 matrices for loss coefficients are given, then the above equation can be modified as:
NG
P Loss PnG Bnm PmG P
n 1
n
G
Bi 0 B00 (57)
JETIRDW06220 Journal of Emerging Technologies and Innovative Research (JETIR) www.jetir.org 1376
© 2019 JETIR January 2019, Volume 6, Issue 1 www.jetir.org (ISSN-2349-5162)
max[ PnG(max) ,(URn PnG )] PnG min[ PnG(max) ,( PnGo DRn )] n 1, 2,3,..., NG (62)
4. CONCLUSION
In the proposed research, the authors has successfully presented the mathematical formulation of single area dynamic
dispatch problem considering plug-in electric vehicles, battery electric vehicles and renewable energy sources (solar
and wind power), which is one of the challenging problems in power system operation control and planning. Single
area economic load dispatch problem of electric power system is limited to electric power generating units of same
area. The proposed mathematical formulation of single area dynamic dispatch will be helpful to the researchers, who
are working in single area economic load dispatch problems with electric vehicles (EVs) i.e. battery electric vehicles,
plug-in electric vehicles and Renewable energy sources (RES) as one of the their research objectives.
REFERENCES
[1] J. S. KOTHARI, D. P., DHILLON, power system optimization, Second. PHI Learning, 2014.
[2] M. Singh and J. S. Dhillon, “Electrical Power and Energy Systems Multiobjective thermal power dispatch using
opposition-based greedy heuristic search,” vol. 82, pp. 339–353, 2016.
[3] N. Jap, J. S. Dhillon, and D. P. Kothari, “Synergic predator-prey optimization for economic thermal power dispatch
problem,” vol. 43, pp. 298–311, 2016.
[4] V. K. Kamboj, S. K. Bath, and J. S. Dhillon, “Solution of non-convex economic load dispatch problem using Grey
Wolf Optimizer,” Neural Comput. Appl., 1934.
[5] D. Singh and J. S. Dhillon, “Ameliorated Grey Wolf Optimization for Economic Load Dispatch,” Energy, 2018.
[6] B. Qu, B. Qiao, Y. Zhu, J. Liang, and L. Wang, “Dynamic power dispatch considering electric vehicles and wind
power using decomposition based multi-objective evolutionary algorithm,” Energies, vol. 10, no. 12, 2017.
[7] P. You, S. Member, Z. Yang, and M. Chow, “Optimal Cooperative Charging Strategy for a Smart Charging Station
of Electric Vehicles,” pp. 1–11, 2015.
[8] S. Pal, S. Member, R. Kumar, and S. Member, “Electric Vehicle Scheduling Strategy in Residential Demand
Response Programs with Neighbor,” vol. 3203, no. c, 2017.
[9] W. Gu, H. Yu, W. Liu, J. Zhu, and X. Xu, “Demand response and economic dispatch of power systems considering
large-scale plug-in hybrid electric vehicles/electric vehicles (phevs/evs): A review,” Energies, vol. 6, no. 9, pp.
4394–4417, 2013.
[10] P. Goyal, A. Sharma, S. Vyas, and R. Kumar, “Customer and Aggregator Balanced Dynamic Electric Vehicle
Charge Scheduling in a Smart Grid Framework,” pp. 276–283, 2016.
[11] M. R. Andervazh and S. Javadi, “Emission-economic dispatch of thermal power generation units in the presence
of hybrid electric vehicles and correlated wind power plants,” IET Gener. Transm. Distrib., vol. 11, no. 9, pp.
2232–2243, 2017.
JETIRDW06220 Journal of Emerging Technologies and Innovative Research (JETIR) www.jetir.org 1377