SlideShare a Scribd company logo
IJRET: International Journal of Research in Engineering and Technology eISSN: 2319-1163 | pISSN: 2321-7308
__________________________________________________________________________________________
Volume: 03 Special Issue: 07 | May-2014, Available @ http://www.ijret.org 320
LOCATIONAL MARGINAL PRICING FRAMEWORK IN SECURED
DISPATCH SCHEDULING UNDER CONTINGENCY CONDITION
R.Manikamdan1
, M.Bhoopathi2
, R.Saravanakumar3
1
PG Student, Department of EEE, Jayaram College of Engineering and Technology, TamilNadu, India
2
Assistant Professor, Department of EEE, Jayaram College of Engineering and Technology, TamilNadu, India
3
Assistant Professor, Department of EEE, Sudharsan Engineering college, TamilNadu, India
Abstract
This paper is to design the locational marginal pricing [LMP] under security constrained condition. The pricing framework of the
regulated power market has same for normal and the contingency condition. So whenever the maximum power used from the
customer, the massive blackouts occurs in the power system network due to exceeds the transmission limits. In deregulated power
market having the LMP pricing method. The LMP is mainly used for contingency condition pricing of load at each location in secured
manner. And also reduce the cost of the minimum load usage customer. Line outage and generation outage is consider in contingency
analysis in security constrained optimal power flow and to calculate LMP in each location (bus) in IEEE-14 bus system by using
power world simulator.
Keywords: Blackouts, Contingency, Line Outage, Generation Outage, Locational Marginal Pricing.
----------------------------------------------------------------------***----------------------------------------------------------------------
1. INTRODUCTION
The recent vertically integrated power market consists of
bundled generation, transmission, distribution [1]. And only
one seller and many buyers in the regulated power system.
The pricing of power transmission is constant at all conditions
and there is no incentive for minimum load user of each
location. So we introduce the deregulation of power system
network to lowering the utility rates, customer specific
services, and encourage the renewable sources. In deregulated
power market is mainly for reduction of cost and to design the
pricing rates of all the load entities and give some incentives
to the minimum load usage customer for demand response
improvement [2]. Due to open access and lower cost
transmission in deregulation, the competition is occurring in
the transmission network. So the transmission of power
exceeds the transfer limits, that time contingency create in the
power system network [1, 2].
The contingency can be classified as line outage, generator
outage, transformer outage etc. under contingency condition
the security of the power system is collapsed, that time
massive blackouts is occurred. So security analysis of power
system network is important task in deregulation [3]. Before
security –constrained optimal power flow we have to find the
contingency analysis form to predict the outages in power
system. And to include the contingency analysis form in the
security-constrained optimal power flow solution [2]
The power system is highly non-linear system which operates
in a constantly changing environment such as load, generator
output, topology and key operating parameters which changes
continuously. Due to secured limits of the power system
network, the cost will be varied by adding the congestion cost
in contingency condition [4]. So the LMP pricing framework
is important in both normal and contingency condition. Most
of the security analysis is based on state estimation of power
system, but we introduce the new approach for security
analysis and LMP pricing calculation [5].
In this paper proposed to new method solution of marginal
cost of the each location (bus) in both normal and contingency
condition and we have to calculate the congestion cost in
contingency condition of each location(bus) is calculated
through security constrained optimal power floe using power
world simulator. The power world simulator is highly essential
tool for marginal cost evaluation in the easy manner. There is
no complicate equation design and coding. The output
response time is maximum than compared with others. Full
Newton’s method is used for the power flow solution and
binding the constraints also included. The mathematical
Problem formulation is in section II. The locational marginal
pricing framework algorithm in power world simulator is in
section III. The description of test system in section IV In
section V includes simulation results and description.
Conclusion from the results in section VI
2. PROBLEM FORMULATION
2.1 Power Flow Equation (N-R Method)
The power world simulator can be set to use a full Newton
solution or use a DC load flow method to analyze each
IJRET: International Journal of Research in Engineering and Technology eISSN: 2319-1163 | pISSN: 2321-7308
__________________________________________________________________________________________
Volume: 03 Special Issue: 07 | May-2014, Available @ http://www.ijret.org 321
contingency. The full Newton approach is not as fast as a DC
load flow, but the results tend to be significantly more
accurate and allow for gauging voltage/VAR effects. The
Newton solution method (also called Newton-Rapson method)
is more efficient for large power systems. The number of
iteration required to obtain a solution is independent of a
system size but more functional evaluation are required at
each iteration
Equation for bus Admittance matrix
Ii = 𝑌𝑛
𝑗 =1 i j Vj (A1)
In above equation j includes bus i expressing this equation in
polar form, we have
Ii = |𝑌𝑛
𝑗 =1 i j||Vj| θij + δj (A2)
The complex power at bus
Pi – Qi = Vi
*
Ii (A3)
Substituting from (2) for Ii in (3)
Pi – Qi = |Vi| ∠ -δi |𝑌𝑛
𝑗 =1 i j||Vj|∠ θij + δj (A4)
Separating the real and imaginary parats
𝑃i= 𝑉𝑖 𝑉𝑗 𝑌𝑖𝑗
𝑛
𝑗 =1 cos 𝜃𝑖𝑗 – 𝛿𝑖+ 𝛿𝑗 (A5)
𝑄𝑖
= − 𝑉𝑖 𝑉𝑗 𝑌𝑖𝑗
𝑛
𝑗=1 sin⁡(𝜃𝑖𝑗 – 𝛿𝑖+ 𝛿𝑗 ) (A6)
Equation (5) and (6) constitute of nonlinear algebraic equation
in terms of the independent variables, voltage magnitude in
per unit and phase angle in radians.
∆𝑃2
⋮
𝑘
∆𝑃 𝑛
𝑘
∆𝑄2
𝑘
⋮
∆𝑄 𝑛
𝑘
=
𝜕𝑃2
𝜕𝛿2
𝑘
⋯
𝜕𝑃2
𝜕𝛿 𝑛
𝑘
⋮ ⋱ ⋮
𝜕𝑃 𝑛
𝜕𝛿2
𝑘
⋯
𝜕𝑃 𝑛
𝜕𝛿 𝑛
𝑘
𝜕𝑃2
𝜕 𝑉2
𝑘
⋯
𝜕𝑃2
𝜕 𝑉𝑛
𝑘
⋮ ⋱ ⋮
𝜕𝑃2
𝜕 𝑉2
𝑘
⋯
𝜕𝑃 𝑛
𝜕 𝑉𝑛
𝑘
𝜕𝑄2
𝜕𝛿2
𝑘
⋯
𝜕𝑄2
𝜕𝛿 𝑛
𝑘
⋮ ⋱ ⋮
𝜕𝑄 𝑛
𝜕𝛿2
𝑘
⋯
𝜕𝑄 𝑛
𝜕𝛿 𝑛
𝑘
𝜕𝑄2
𝜕 𝑉2
𝑘
⋯
𝜕𝑄2
𝜕 𝑉2
𝑘
⋮ ⋱ ⋮
𝜕𝑄 𝑛
𝜕 𝑉2
𝑘
⋯
𝜕𝑄2
𝜕 𝑉𝑛
𝑘
∆𝛿2
⋮
𝑘
∆𝛿 𝑛
𝑘
∆|𝑉2
𝑘
|
⋮
∆|𝑉𝑛
𝑘
|
(A7)
In above equation, bus 1 is assumed to be slack bus. The
jacobian matrix gives the linearized relationship between
small changes in voltage angle ∆𝛿𝑖
(𝑘)
and voltage magnitude
Δ|𝑉𝑖
(𝑘)
| with small changes in real and reactive power Δ𝑃𝑖
(𝑘)
and Q i
(k)
elements of jacobian matrix are the partial
derivatives of (5) and (6) evaluated at ∆𝛿𝑖
(𝑘)
and Δ|𝑉𝑖
(𝑘)
|.
𝛥𝑃
𝛥𝑄
=
𝐽1 𝐽2
𝐽3 𝐽4
𝛥𝛿
𝛥|𝑉|
(A8)
Accordingly there are (n-1) real power constraints and (n-1-m)
reactive power constraints and the jacobian matrix is the order
of (2n-2-m) (2n-2-m).
J1 is the order of (n-1) x (n-1)
𝜕𝑃 𝑖
𝜕𝛿 𝑖
= 𝑉𝑖 𝑉𝑗 𝑌𝑖𝑗
𝑛
𝑗 ≠1 sin⁡(𝜃𝑖𝑗 – 𝛿𝑖+ 𝛿𝑗 ) (A9)
𝜕𝑃 𝑖
𝜕𝛿 𝑗
= − 𝑉𝑖 Vj Yij sin 𝜃𝑖𝑗 – 𝛿𝑖+ 𝛿𝑗 j≠ 1 (A10)
J2 is the order of (n-1) x (n-1-m)
𝜕𝑃 𝑖
𝜕|𝑉 𝑖|
= 2 𝑉𝑖 𝑌𝑖𝑖 𝑐𝑜𝑠𝜃𝑖𝑖 + 𝑉𝑖 𝑉𝑗 𝑌𝑖𝑗
𝑛
𝑗≠1 cos⁡(𝜃𝑖𝑗 – 𝛿𝑖+ 𝛿𝑗 )
(A11)
𝜕𝑃 𝑖
𝜕|𝑉 𝑗 |
= 𝑉𝑖 𝑌𝑖𝑗 𝑐 os 𝜃𝑖𝑗 – 𝛿𝑖+ 𝛿𝑗 j≠ 𝑖 (A12)
J3 is the order of (n-1-m) x (n-1)
𝜕𝑄 𝑖
𝜕𝛿 𝑖
= 𝑉𝑖 𝑉𝑗 𝑌𝑖𝑗
𝑛
𝑗 ≠1 cos 𝜃𝑖𝑗 – 𝛿𝑖+ 𝛿𝑗 (A13)
J4 is the order of (n-1-m) x (n-1-m)
𝜕𝑄 𝑖
𝜕𝛿 𝑗
= − 𝑉𝑖 Vj Yij cos⁡(𝜃𝑖𝑗 – 𝛿𝑖+ 𝛿𝑗 )j≠ 𝑖 (A14)
𝜕𝑄 𝑖
𝜕|𝑉 𝑖|
= 𝑉𝑖 𝑌𝑖𝑖 𝑠𝑖𝑛𝜃𝑖𝑖 − 𝑉𝑗 𝑌𝑖𝑗
𝑛
𝑗≠1 sin 𝜃𝑖𝑗 – 𝛿𝑖+ 𝛿𝑗 (A15)
𝜕𝑄 𝑖
𝜕|𝑉 𝑗 |
= − 𝑉𝑖 𝑌𝑖𝑗 𝑠𝑖𝑛⁡(𝜃𝑖𝑗 – 𝛿𝑖+ 𝛿𝑗 )j≠ 𝑖 (A16)
The terms Δ𝑃𝑖
(𝑘)
and Q i
(k)
are difference between the
schedule and calculated values, known as the power residuals,
given by
Δ𝑃𝑖
𝑘
= Pi
sch
− Pi
k
(A17)
Δ𝑄𝑖
𝑘
= Qi
sch
− Qi
k
(A18)
The new estimated for bus voltage i
𝛿𝑖
(𝑘+1)
= 𝛿𝑖
(𝑘)
+ ∆𝛿𝑖
(𝑘)
(A19)
IJRET: International Journal of Research in Engineering and Technology eISSN: 2319-1163 | pISSN: 2321-7308
__________________________________________________________________________________________
Volume: 03 Special Issue: 07 | May-2014, Available @ http://www.ijret.org 322
|𝑉𝑖
(𝑘+1)
| = |𝑉𝑖
𝑘
| + ∆|𝑉𝑖
(𝑘)
| (A20)
2.2 Security Constrained Optimal Power Flow
Objective function
Min f ( p ) (A21)
Subject to g (p) = 0
hmin ≤ h(p) ≤ hmax in normal condition constraints
h’min ≤ h’ (p) ≤ h’max in contingency condition constraints.
Security constrained optimal power flow solution (SCOPF)
will always a cost ≥ optimal power flow cost. If we ignore
losses, then we can say that an OPF solution differs from an
EDC solution only when a normal transmission constraint.
When normal flow moves from just < 100% to ≥ 100% of
continuous rating
SCOPF differs from an OPF solution only when contingency
constraint becomes binding occurs when post- contingency
flow moves from just < 100% to ≥ 100% of emergency rating.
Now let’s consider the SCOPF. Its problem statement is given
as problem Pp :
Min f ( x0, u0 ) (A22)
gk ( xk , uk ) =0 k=0,1,2….c
hk ( xk, uk ) = hk
max
k=0,1,2,….c
Notice that there are C contingencies to be addressed in the
SCOPF and that there are a complete new set of constraints for
each of these C contingencies observe. Each of contingency
related equality constraints is exactly like the original set of
equality constraints except if corresponds the system with an
element removed.
Each set of contingency related inequality-constraints is
exactly like the original set of inequality constraints except its
corresponds to the system with an element removed and
branch flow constraints and for voltage magnitudes, the limits
will be different.
Also notice that the constraints are a function of xk , the
voltage magnitudes and angles under the pre-contingency
(k=0) and contingency condition ( k > 1,2,….c) and u0.
2.3 Locational Marginal Pricing Calculation
Locational marginal pricing (LMPs) are determined from the
result of a security-constrained least-cost dispatch. It is a taxi
ride for MW. It may differs in the various location(bus). We
need two factors to deside the locational marginal pricing.
(i)Transmission congestion
(ii)Losses
The Locational Marginal Pricing (nodel price) at bus i can be
calculated using the following equation
LMP = marginalcost + congestioncost + lossescost
λi = λRef + λCongest + λLossi (A23)
λi = λRef - Li x λRef - 𝜇𝑗 𝑋 𝑆𝐹𝑗𝑖𝑗 (A24)
λLossi = (- Li x λRef) - losses from the reference bus to bus i
= (+ Li x λRef ) – losses from bus i to reference bus
λCongest = (- 𝜇𝑗 𝑋 𝑆𝐹𝑗𝑖𝑗 - congestion from reference bus
to bus i
= (+ 𝜇𝑗 𝑋 𝑆𝐹𝑗𝑖𝑗 - congestion from bus I to
reference bus.
3. LMP ALGORITHM IN POWER WORLD
SIMULATOR
3.1. Locational Marginal Pricing Algorithm
Step 1: Draw the simulink one line diagram in new case
window of power world simulator for the given power system
in edit mode.
Step 2: Set the cubic cost of each generation and to convert
piece wise linear cost.
Step 3: Save the case with apt name.
Step 4: Select tools in run mode and to solve the power flow
by using full N-R method.
Step 5: Open Add-ons Ribbon Tab
Step 6: To select the OPF case information of the
Dialog box and to select the all OPF area records.
a.) If OPF records “YES” that record is included in the
marginal cost calculation.
b.) If OPF records “NO” that records is not include in the
marginal cost calculation.
Step 7: To set all the OPF constraints and also include
common constraints.
Step 8: Open the SCOPF dialog box in the add-ons ribbon
tab.
Step 9: Run Full security constrained OPF under normal
condition (zero contingency in contingency analysis form)
Step10: To calculate marginal cost of each bus (location)
before contingency.
Step11: To view contingency analysis form in the SCOPF
dialog box.
Step12: Right click on label and select auto insert
contingencies through insert special option.
Step 13: Verify that single transmission line or transformer is
selected.
Step 14: If can limit the contingencies inserted to only those
meeting define filter.
Step 15: we want to insert contingencies for all branches and
generators so no filtering is desired.
Step 16: To check the following conditions
a.) Remove the checkmarks in use area/zone filters.
IJRET: International Journal of Research in Engineering and Technology eISSN: 2319-1163 | pISSN: 2321-7308
__________________________________________________________________________________________
Volume: 03 Special Issue: 07 | May-2014, Available @ http://www.ijret.org 323
b.) Verify no other options are selected.
Step 17: Click Do insert contingencies button to accept the all
contingencies.
Step 18: Click “ YES” to get the contingencies.
Step 19: Now the contingency analysis dialog shows
contingencies
a.) Right click on the list display on the contingency
tap and select insert special and click auto insert to
the local menu.
b.) Select single generating unit then click the do insert
contingency button. Click “YES” to complete.
Step 20: Click “Start Run” on the contingency tab and close
the tab.
Step 21: The contingency elements are include in the SCOPF
dialog box.
Step 21: To Run the Full Security Constrained optimal power
Flow.
Step 22: We get the marginal cost under contingency
condition. In addition congestion, transmission loss cost in
each location (bus)
4. TEST SYSTEM
The Locational marginal pricing of 14-bus test system is
shown below when the power flow is running on the power
world simulator. The percentage of power flow is mentioned
in power flow diagram. It consists of five generators for
dispatch of power.
Fig. 1 IEEE-14 bus system
5. SIMULATION RESULTS
Locational marginal pricing (LMP) is important in deregulated
power market under normal and contingency condition. So we
calculate the LMP by considering without contingency, line
outage, generation outage and also include all contingency.
First we find the contingency elements in IEEE-14 bus
system.
Table.1 Single contingency of line & generator
Label Violations
Max
Branch %
Min Voltage
L_0000011-
0000022C1
2 276.5 -
L_0000011-
0000055C1
1 127.5 -
L_0000022-
0000033C1
2 102.0 -
L_0000022-
0000044C1
1 102.0 -
L_0000022-
0000055C1
1 103.0 -
L_0000066-
00001313C1
2 - 0.898
L_0000099-
00001414C1
1 - 0.848
G_0000022U
1
1 103.6 -
Table.2 Multiple contingency of both line and Generator
Label Violations
Max
Branch
%
Min Voltage
G_0000022u1&
L_0000022-
0000033c1
3 115.7 -
G_0000022u1&
L_0000066-
00001313c1
3 104.6 0.894
G_0000033u1&
L_0000011-
0000022c1
4 325.3 0.891
G_0000033u1&
L_0000011-
0000055c1
3 140.2 -
G_0000033u1&
L_0000022-
0000033c1
6 118.0 0.739
G_0000066u1&
L_0000011-
0000022c1
3 283.4 0.883
G_0000066u1&
L_0000077-
5 - 0.842
IJRET: International Journal of Research in Engineering and Technology eISSN: 2319-1163 | pISSN: 2321-7308
__________________________________________________________________________________________
Volume: 03 Special Issue: 07 | May-2014, Available @ http://www.ijret.org 324
0000099c1
G_0000066u1&
L_0000099-
000001414c1
3 - 0.733
Table 3 LMP in Normal condition
Locational marginal pricing of each bus system is in the above
table. Ther is no congestion cost due to normal power flow
condition.
Table 4 LMP in Contingency Condition (only includes single
contingency).
Bus
No.
Area
name
Energy
cost
($/mwh)
Conge
cost
($/mw
h)
Loss
cost
($/mwh
)
LMP
($/mwh)
1 Top 5.86 0.00 0.00 5.86
2 Top 5.86 5.88 0.26 12.00
3 Top 5.86 593.05 0.78 599.69
4 Top 5.86 530.35 0.74 536.95
5 Top 5.86 389.45 0.64 395.94
6 Top 5.86 432.90 0.64 439.40
7 Top 5.86 510.31 0.75 516.92
8 Top 5.86 510.31 0.75 516.92
9 Top 5.86 499.13 0.76 505.74
10 Top 5.86 490.63 0.78 497.27
11 Top 5.86 463.86 0.73 470.46
12 Top 5.86 446.39 0.78 453.03
13 Top 5.86 456.67 0.86 463.40
14 Top 5.86 509.62 1.21 516.68
The congestion cost is occurred in each bus due to single
contingency of line outage and generation outage of each
elements
Table 5 LMP price (includes all contingencies)
Bus
No.
Area
name
Energy
cost
($/mwh)
Conge
Cost
($/mwh)
Losses
($/mwh)
LMP
($/mwh)
1 Top -1161.26 1167.12 0.00 5.86
2 Top -1161.26 1223.32 -50.06 12.00
3 Top -1161.26 1610.35 -152.80 296.29
4 Top -1161.26 3318.47 -144.14 2013.07
5 Top -1161.26 1123.77 -124.52 -162.01
6 Top -1161.26 14055.13 -124.74 12769.1
3
7 Top -1161.26 3572.12 -146.98 2263.83
8 Top -1161.26 3572.12 -146.98 2263.88
9 Top -1161.26 5166.58 -148.48 3856.84
10 Top -1161.26 13749.70 -152.99 12435.4
5
11 Top -1161.26 13944.48 -143.86 12639.3
6
12 Top -1161.26 14223.56 -152.48 12909.8
2
13 Top -1161.26 14038.94 -169.60 12708.0
7
14 Top -1161.26 12030.02 -237.59 10631.1
7
In security constrained condition all single and multiple
contingency elements are include in the power system. The
LMP value is maximum due to congestion cost.
Fig. 2 LMP Cost Curve In Each Bus (before contingencies)
Bus
No.
Area
name
Energy
cost
($/mwh)
Cong
cost
($/m
wh)
Loss
cost
($/mwh)
LMP
($/mwh)
1 Top 11.28 -0.00 0.00 11.28
2 Top 11.28 -0.00 0.72 12.00
3 Top 11.28 -0.00 1.72 12.99
4 Top 11.28 -0.00 1.60 12.87
5 Top 11.28 -0.00 1.39 12.66
6 Top 11.28 -0.00 1.39 12.67
7 Top 11.28 -0.00 1.63 12.90
8 Top 11.28 -0.00 1.63 12.90
9 Top 11.28 -0.00 1.64 12.92
10 Top 11.28 -0.00 1.68 12.96
11 Top 11.28 -0.00 1.59 12.86
12 Top 11.28 -0.00 1.67 12.94
13 Top 11.28 -0.00 1.84 13.11
14 Top 11.28 -0.00 2.51 13.79
IJRET: International Journal of Research in Engineering and Technology eISSN: 2319-1163 | pISSN: 2321-7308
__________________________________________________________________________________________
Volume: 03 Special Issue: 07 | May-2014, Available @ http://www.ijret.org 325
Fig.3 Lmp Cost Curve in Each Bus (only includes line outage
and generator outage)
Fig .4 Lmp Cost Curve In Each Bus (includes all
contingencies)
6. CONCLUSIONS
The state estimation based contingency analysis technique is
more complex. But Newton’s method based contingency
analysis and security constrained optimal power flow in
contingency condition is simple and more accuracy in power
world simulator. The LMP calculation of test system is easy in
power world simulation. In addition LMP is depends on
congestion cost of each bus system. Whenever the
contingency is occur in this system that time only congestion
cost is added in LMP. Otherwise there is no congestion cost in
the system. So the locational marginal pricing is reduced in the
normal condition. The LMP calculation is helpful for demand
response improvement under security condition in deregulated
power market. In future we can easily reduce the cost of the
power transmission and to improve the demand response in
secured manner by using the reserve option and to connect
from renewable generation to the transmission network. And
we get minimum congestion cost in each location. So we can
easily improve the demand response in the secured dispatch
scheduling in both normal and the contingency condition in
the power world simulator tool. It is user friendly software in
the pricing calculation of deregulated power market in the
power system network.
REFERENCES
[1]. Amit Kumar Roy “contingency analysis in power system”
Thapur University Patiala.2009.
[2]. R.Manikandan, M.Bhoopathi “ Contingency analysis in
deregulated power market” Jayaram College of Engineering
and Technology, Tamilnadu 2013.
[3]. Richard D.Christie, AnjanBose “load frequency control
issues in power system operations after deregulation”
University of Washington, 1995.
[4]. Saavedra, O.R., "Solving the security constrained optimal
power flow problem in a distributed computing environment,"
Generation, transmission and distribution, IEEE proceedings,
vol.143, No.6 pp 593-598, 1996.
[5]. Fangxing Li, Rui Bo, “ Congestion and Price Prediction
Under Load Variation,”IEEE Transaction on power system ,
Vol 24,No.2 May 2009.
BIOGRAPHIES
Manikandan R obtained his Bachelor
degree in Electrical & Electronics
Engineering from P.R.Engineering
College, Thanjavur in the year 2012 and
Master degree in Power Systems
Engineering doing from Jayaram College of
Engineering and technology, Thuraiyur,
India. His research area includes deregulation of power market
and Smart grid technologies.
IJRET: International Journal of Research in Engineering and Technology eISSN: 2319-1163 | pISSN: 2321-7308
__________________________________________________________________________________________
Volume: 03 Special Issue: 07 | May-2014, Available @ http://www.ijret.org 326
Bhoopathi M received his B.E. degree in
Electrical and Electronics Engineering from
Kumaraguru College of technology,
Bharathiyar University, Coimbatore, India
and M.E. degree in Power systems
engineering from Annamalai University,
Chidambaram, India. He is currently
working as an Assistant Professor in Department of Electrical
and Electronics Engineering, Jayaram college of engineering
and technology, Pagalavadi, Thuraiyur , India. His research
interest includes Restructured power system and smart grid
technologies.
Saravanakumar R received his B.E.
degree in Electrical and Electronics
Engineering from PSNA College of
engineering & technology, Dindigul, Anna
University, Chennai, India and M.E.
degree in Power systems engineering from
University College of engineering, BIT
campus, Trichirappalli, India. He is
currently working as an Assistant Professor in Department of
Electrical and Electronics Engineering,Sudharsan engineering
college, Sathiyamangalam, India. His research interest
includes power systems and Distributed Generation.

More Related Content

PDF
Explicit model predictive control of fast dynamic system
PDF
Comparison of cascade P-PI controller tuning methods for PMDC motor based on ...
PDF
GWO-based estimation of input-output parameters of thermal power plants
PDF
Security constrained optimal load dispatch using hpso technique for thermal s...
PDF
40220140503006
PDF
Firefly Algorithm to Opmimal Distribution of Reactive Power Compensation Units
PDF
IRJET- A Genetic based Stochastic Approach for Solving Thermal Unit Commitmen...
PDF
Security Constrained UCP with Operational and Power Flow Constraints
Explicit model predictive control of fast dynamic system
Comparison of cascade P-PI controller tuning methods for PMDC motor based on ...
GWO-based estimation of input-output parameters of thermal power plants
Security constrained optimal load dispatch using hpso technique for thermal s...
40220140503006
Firefly Algorithm to Opmimal Distribution of Reactive Power Compensation Units
IRJET- A Genetic based Stochastic Approach for Solving Thermal Unit Commitmen...
Security Constrained UCP with Operational and Power Flow Constraints

What's hot (17)

PDF
Design of Quadratic Optimal Regulator for DC Motor
PDF
Multi objective economic load dispatch using hybrid fuzzy, bacterial
PDF
Centralized Optimal Control for a Multimachine Power System Stability Improve...
PDF
Prakash narendra
PDF
Comparison of backstepping, sliding mode and PID regulators for a voltage inv...
PDF
HYBRID PARTICLE SWARM OPTIMIZATION FOR SOLVING MULTI-AREA ECONOMIC DISPATCH P...
PDF
IJERD (www.ijerd.com) International Journal of Engineering Research and Devel...
PDF
Dynamic Economic Dispatch Assessment Using Particle Swarm Optimization Technique
PDF
NOVEL PSO STRATEGY FOR TRANSMISSION CONGESTION MANAGEMENT
PDF
Fuzzy logic based direct torque control of induction motor with space vector ...
PDF
Iaetsd estimation of damping torque for small-signal
PDF
Fuzzy logic based direct torque control of induction motor with space vector ...
PPTX
Divyesh_14Me63r02(new)
PDF
Reduction of Active Power Loss byUsing Adaptive Cat Swarm Optimization
PDF
Gy3312241229
PDF
The Neural Network-Combined Optimal Control System of Induction Motor
Design of Quadratic Optimal Regulator for DC Motor
Multi objective economic load dispatch using hybrid fuzzy, bacterial
Centralized Optimal Control for a Multimachine Power System Stability Improve...
Prakash narendra
Comparison of backstepping, sliding mode and PID regulators for a voltage inv...
HYBRID PARTICLE SWARM OPTIMIZATION FOR SOLVING MULTI-AREA ECONOMIC DISPATCH P...
IJERD (www.ijerd.com) International Journal of Engineering Research and Devel...
Dynamic Economic Dispatch Assessment Using Particle Swarm Optimization Technique
NOVEL PSO STRATEGY FOR TRANSMISSION CONGESTION MANAGEMENT
Fuzzy logic based direct torque control of induction motor with space vector ...
Iaetsd estimation of damping torque for small-signal
Fuzzy logic based direct torque control of induction motor with space vector ...
Divyesh_14Me63r02(new)
Reduction of Active Power Loss byUsing Adaptive Cat Swarm Optimization
Gy3312241229
The Neural Network-Combined Optimal Control System of Induction Motor
Ad

Viewers also liked (20)

PDF
Economical placement of shear walls in a moment resisting frame for earthquak...
PDF
Application of ibearugbulem’s model for optimizing granite concrete mix
PDF
Jit dynamic cryptosystem
PDF
Design of dual master i2 c bus controller
PDF
A novel way of verifiable redistribution of the secret in a multiuser environ...
PDF
Effect of tuck loop in bursting strength of single jersey knitted fabrics
PDF
Seismic response of reinforced concrete structure by using different bracing ...
PDF
Design of digital signature verification algorithm using relative slope method
PDF
Pile response due to earthquake induced lateral spreading
PDF
Adaptive approach to retrieve image affected by impulse noise
PDF
Detection and identification of chemical agent using atomic absorption spectr...
PDF
Performance improvement of energy aware and adaptive routing protocols for ma...
PDF
Comparison of friction stirs welding technique with conventional welding methods
PDF
Issues in optimizing the performance of wireless sensor networks
PDF
Prevention against new cell counting attack against tor network
PDF
Measuring effort for modifying software package as
PDF
An innovative way for computerized smith chart generation and transmission li...
PDF
Density and compaction characteristics of wma using additives
PDF
Real time object tracking and learning using template matching
PDF
System of quasilinear equations of reaction diffusion
Economical placement of shear walls in a moment resisting frame for earthquak...
Application of ibearugbulem’s model for optimizing granite concrete mix
Jit dynamic cryptosystem
Design of dual master i2 c bus controller
A novel way of verifiable redistribution of the secret in a multiuser environ...
Effect of tuck loop in bursting strength of single jersey knitted fabrics
Seismic response of reinforced concrete structure by using different bracing ...
Design of digital signature verification algorithm using relative slope method
Pile response due to earthquake induced lateral spreading
Adaptive approach to retrieve image affected by impulse noise
Detection and identification of chemical agent using atomic absorption spectr...
Performance improvement of energy aware and adaptive routing protocols for ma...
Comparison of friction stirs welding technique with conventional welding methods
Issues in optimizing the performance of wireless sensor networks
Prevention against new cell counting attack against tor network
Measuring effort for modifying software package as
An innovative way for computerized smith chart generation and transmission li...
Density and compaction characteristics of wma using additives
Real time object tracking and learning using template matching
System of quasilinear equations of reaction diffusion
Ad

Similar to Locational marginal pricing framework in secured dispatch scheduling under contingency condition (20)

PDF
Optimal placement of distributed power flow controller for loss reduction usi...
PDF
40220140504003
PDF
F43022431
PDF
IRJET- Optimal Placement and Size of DG and DER for Minimizing Power Loss and...
PPT
P1111145969
PDF
Model Order Reduction of an ISLANDED MICROGRID using Single Perturbation, Dir...
PDF
Congestion Management in Power System by Optimal Location And Sizing of UPFC
PDF
Explicit model predictive control of fast dynamic system
PDF
Intelligent fault diagnosis for power distribution systemcomparative studies
PDF
IRJET- Location Identification for FACTs Device
PDF
IRJET- Optimal Placement and Size of DG and DER for Minimizing Power Loss and...
PDF
Power loss reduction in radial distribution system by using plant growth simu...
PDF
11.power loss reduction in radial distribution system by using plant growth s...
DOCX
Final report Review
PDF
Economic Load Dispatch Optimization of Six Interconnected Generating Units Us...
PDF
Passerine swarm optimization algorithm for solving optimal reactive power dis...
PDF
Hysteresis Current Control Based Shunt Active Power Filter for Six Pulse Ac/D...
PDF
Power quality improvement in a weak bus system using
PDF
Transmission Loss Minimization Using Optimization Technique Based On Pso
Optimal placement of distributed power flow controller for loss reduction usi...
40220140504003
F43022431
IRJET- Optimal Placement and Size of DG and DER for Minimizing Power Loss and...
P1111145969
Model Order Reduction of an ISLANDED MICROGRID using Single Perturbation, Dir...
Congestion Management in Power System by Optimal Location And Sizing of UPFC
Explicit model predictive control of fast dynamic system
Intelligent fault diagnosis for power distribution systemcomparative studies
IRJET- Location Identification for FACTs Device
IRJET- Optimal Placement and Size of DG and DER for Minimizing Power Loss and...
Power loss reduction in radial distribution system by using plant growth simu...
11.power loss reduction in radial distribution system by using plant growth s...
Final report Review
Economic Load Dispatch Optimization of Six Interconnected Generating Units Us...
Passerine swarm optimization algorithm for solving optimal reactive power dis...
Hysteresis Current Control Based Shunt Active Power Filter for Six Pulse Ac/D...
Power quality improvement in a weak bus system using
Transmission Loss Minimization Using Optimization Technique Based On Pso

More from eSAT Publishing House (20)

PDF
Likely impacts of hudhud on the environment of visakhapatnam
PDF
Impact of flood disaster in a drought prone area – case study of alampur vill...
PDF
Hudhud cyclone – a severe disaster in visakhapatnam
PDF
Groundwater investigation using geophysical methods a case study of pydibhim...
PDF
Flood related disasters concerned to urban flooding in bangalore, india
PDF
Enhancing post disaster recovery by optimal infrastructure capacity building
PDF
Effect of lintel and lintel band on the global performance of reinforced conc...
PDF
Wind damage to trees in the gitam university campus at visakhapatnam by cyclo...
PDF
Wind damage to buildings, infrastrucuture and landscape elements along the be...
PDF
Shear strength of rc deep beam panels – a review
PDF
Role of voluntary teams of professional engineers in dissater management – ex...
PDF
Risk analysis and environmental hazard management
PDF
Review study on performance of seismically tested repaired shear walls
PDF
Monitoring and assessment of air quality with reference to dust particles (pm...
PDF
Low cost wireless sensor networks and smartphone applications for disaster ma...
PDF
Coastal zones – seismic vulnerability an analysis from east coast of india
PDF
Can fracture mechanics predict damage due disaster of structures
PDF
Assessment of seismic susceptibility of rc buildings
PDF
A geophysical insight of earthquake occurred on 21 st may 2014 off paradip, b...
PDF
Effect of hudhud cyclone on the development of visakhapatnam as smart and gre...
Likely impacts of hudhud on the environment of visakhapatnam
Impact of flood disaster in a drought prone area – case study of alampur vill...
Hudhud cyclone – a severe disaster in visakhapatnam
Groundwater investigation using geophysical methods a case study of pydibhim...
Flood related disasters concerned to urban flooding in bangalore, india
Enhancing post disaster recovery by optimal infrastructure capacity building
Effect of lintel and lintel band on the global performance of reinforced conc...
Wind damage to trees in the gitam university campus at visakhapatnam by cyclo...
Wind damage to buildings, infrastrucuture and landscape elements along the be...
Shear strength of rc deep beam panels – a review
Role of voluntary teams of professional engineers in dissater management – ex...
Risk analysis and environmental hazard management
Review study on performance of seismically tested repaired shear walls
Monitoring and assessment of air quality with reference to dust particles (pm...
Low cost wireless sensor networks and smartphone applications for disaster ma...
Coastal zones – seismic vulnerability an analysis from east coast of india
Can fracture mechanics predict damage due disaster of structures
Assessment of seismic susceptibility of rc buildings
A geophysical insight of earthquake occurred on 21 st may 2014 off paradip, b...
Effect of hudhud cyclone on the development of visakhapatnam as smart and gre...

Recently uploaded (20)

PDF
Evaluating the Democratization of the Turkish Armed Forces from a Normative P...
PPTX
Lecture Notes Electrical Wiring System Components
PDF
R24 SURVEYING LAB MANUAL for civil enggi
PDF
BMEC211 - INTRODUCTION TO MECHATRONICS-1.pdf
PPTX
Sustainable Sites - Green Building Construction
PDF
Mohammad Mahdi Farshadian CV - Prospective PhD Student 2026
PPTX
CARTOGRAPHY AND GEOINFORMATION VISUALIZATION chapter1 NPTE (2).pptx
PDF
The CXO Playbook 2025 – Future-Ready Strategies for C-Suite Leaders Cerebrai...
PPT
CRASH COURSE IN ALTERNATIVE PLUMBING CLASS
PPTX
CH1 Production IntroductoryConcepts.pptx
PPTX
FINAL REVIEW FOR COPD DIANOSIS FOR PULMONARY DISEASE.pptx
PDF
Model Code of Practice - Construction Work - 21102022 .pdf
PDF
PRIZ Academy - 9 Windows Thinking Where to Invest Today to Win Tomorrow.pdf
PDF
Digital Logic Computer Design lecture notes
PPTX
MET 305 2019 SCHEME MODULE 2 COMPLETE.pptx
PDF
composite construction of structures.pdf
PDF
Embodied AI: Ushering in the Next Era of Intelligent Systems
PPTX
Recipes for Real Time Voice AI WebRTC, SLMs and Open Source Software.pptx
PPTX
UNIT 4 Total Quality Management .pptx
PDF
Mitigating Risks through Effective Management for Enhancing Organizational Pe...
Evaluating the Democratization of the Turkish Armed Forces from a Normative P...
Lecture Notes Electrical Wiring System Components
R24 SURVEYING LAB MANUAL for civil enggi
BMEC211 - INTRODUCTION TO MECHATRONICS-1.pdf
Sustainable Sites - Green Building Construction
Mohammad Mahdi Farshadian CV - Prospective PhD Student 2026
CARTOGRAPHY AND GEOINFORMATION VISUALIZATION chapter1 NPTE (2).pptx
The CXO Playbook 2025 – Future-Ready Strategies for C-Suite Leaders Cerebrai...
CRASH COURSE IN ALTERNATIVE PLUMBING CLASS
CH1 Production IntroductoryConcepts.pptx
FINAL REVIEW FOR COPD DIANOSIS FOR PULMONARY DISEASE.pptx
Model Code of Practice - Construction Work - 21102022 .pdf
PRIZ Academy - 9 Windows Thinking Where to Invest Today to Win Tomorrow.pdf
Digital Logic Computer Design lecture notes
MET 305 2019 SCHEME MODULE 2 COMPLETE.pptx
composite construction of structures.pdf
Embodied AI: Ushering in the Next Era of Intelligent Systems
Recipes for Real Time Voice AI WebRTC, SLMs and Open Source Software.pptx
UNIT 4 Total Quality Management .pptx
Mitigating Risks through Effective Management for Enhancing Organizational Pe...

Locational marginal pricing framework in secured dispatch scheduling under contingency condition

  • 1. IJRET: International Journal of Research in Engineering and Technology eISSN: 2319-1163 | pISSN: 2321-7308 __________________________________________________________________________________________ Volume: 03 Special Issue: 07 | May-2014, Available @ http://www.ijret.org 320 LOCATIONAL MARGINAL PRICING FRAMEWORK IN SECURED DISPATCH SCHEDULING UNDER CONTINGENCY CONDITION R.Manikamdan1 , M.Bhoopathi2 , R.Saravanakumar3 1 PG Student, Department of EEE, Jayaram College of Engineering and Technology, TamilNadu, India 2 Assistant Professor, Department of EEE, Jayaram College of Engineering and Technology, TamilNadu, India 3 Assistant Professor, Department of EEE, Sudharsan Engineering college, TamilNadu, India Abstract This paper is to design the locational marginal pricing [LMP] under security constrained condition. The pricing framework of the regulated power market has same for normal and the contingency condition. So whenever the maximum power used from the customer, the massive blackouts occurs in the power system network due to exceeds the transmission limits. In deregulated power market having the LMP pricing method. The LMP is mainly used for contingency condition pricing of load at each location in secured manner. And also reduce the cost of the minimum load usage customer. Line outage and generation outage is consider in contingency analysis in security constrained optimal power flow and to calculate LMP in each location (bus) in IEEE-14 bus system by using power world simulator. Keywords: Blackouts, Contingency, Line Outage, Generation Outage, Locational Marginal Pricing. ----------------------------------------------------------------------***---------------------------------------------------------------------- 1. INTRODUCTION The recent vertically integrated power market consists of bundled generation, transmission, distribution [1]. And only one seller and many buyers in the regulated power system. The pricing of power transmission is constant at all conditions and there is no incentive for minimum load user of each location. So we introduce the deregulation of power system network to lowering the utility rates, customer specific services, and encourage the renewable sources. In deregulated power market is mainly for reduction of cost and to design the pricing rates of all the load entities and give some incentives to the minimum load usage customer for demand response improvement [2]. Due to open access and lower cost transmission in deregulation, the competition is occurring in the transmission network. So the transmission of power exceeds the transfer limits, that time contingency create in the power system network [1, 2]. The contingency can be classified as line outage, generator outage, transformer outage etc. under contingency condition the security of the power system is collapsed, that time massive blackouts is occurred. So security analysis of power system network is important task in deregulation [3]. Before security –constrained optimal power flow we have to find the contingency analysis form to predict the outages in power system. And to include the contingency analysis form in the security-constrained optimal power flow solution [2] The power system is highly non-linear system which operates in a constantly changing environment such as load, generator output, topology and key operating parameters which changes continuously. Due to secured limits of the power system network, the cost will be varied by adding the congestion cost in contingency condition [4]. So the LMP pricing framework is important in both normal and contingency condition. Most of the security analysis is based on state estimation of power system, but we introduce the new approach for security analysis and LMP pricing calculation [5]. In this paper proposed to new method solution of marginal cost of the each location (bus) in both normal and contingency condition and we have to calculate the congestion cost in contingency condition of each location(bus) is calculated through security constrained optimal power floe using power world simulator. The power world simulator is highly essential tool for marginal cost evaluation in the easy manner. There is no complicate equation design and coding. The output response time is maximum than compared with others. Full Newton’s method is used for the power flow solution and binding the constraints also included. The mathematical Problem formulation is in section II. The locational marginal pricing framework algorithm in power world simulator is in section III. The description of test system in section IV In section V includes simulation results and description. Conclusion from the results in section VI 2. PROBLEM FORMULATION 2.1 Power Flow Equation (N-R Method) The power world simulator can be set to use a full Newton solution or use a DC load flow method to analyze each
  • 2. IJRET: International Journal of Research in Engineering and Technology eISSN: 2319-1163 | pISSN: 2321-7308 __________________________________________________________________________________________ Volume: 03 Special Issue: 07 | May-2014, Available @ http://www.ijret.org 321 contingency. The full Newton approach is not as fast as a DC load flow, but the results tend to be significantly more accurate and allow for gauging voltage/VAR effects. The Newton solution method (also called Newton-Rapson method) is more efficient for large power systems. The number of iteration required to obtain a solution is independent of a system size but more functional evaluation are required at each iteration Equation for bus Admittance matrix Ii = 𝑌𝑛 𝑗 =1 i j Vj (A1) In above equation j includes bus i expressing this equation in polar form, we have Ii = |𝑌𝑛 𝑗 =1 i j||Vj| θij + δj (A2) The complex power at bus Pi – Qi = Vi * Ii (A3) Substituting from (2) for Ii in (3) Pi – Qi = |Vi| ∠ -δi |𝑌𝑛 𝑗 =1 i j||Vj|∠ θij + δj (A4) Separating the real and imaginary parats 𝑃i= 𝑉𝑖 𝑉𝑗 𝑌𝑖𝑗 𝑛 𝑗 =1 cos 𝜃𝑖𝑗 – 𝛿𝑖+ 𝛿𝑗 (A5) 𝑄𝑖 = − 𝑉𝑖 𝑉𝑗 𝑌𝑖𝑗 𝑛 𝑗=1 sin⁡(𝜃𝑖𝑗 – 𝛿𝑖+ 𝛿𝑗 ) (A6) Equation (5) and (6) constitute of nonlinear algebraic equation in terms of the independent variables, voltage magnitude in per unit and phase angle in radians. ∆𝑃2 ⋮ 𝑘 ∆𝑃 𝑛 𝑘 ∆𝑄2 𝑘 ⋮ ∆𝑄 𝑛 𝑘 = 𝜕𝑃2 𝜕𝛿2 𝑘 ⋯ 𝜕𝑃2 𝜕𝛿 𝑛 𝑘 ⋮ ⋱ ⋮ 𝜕𝑃 𝑛 𝜕𝛿2 𝑘 ⋯ 𝜕𝑃 𝑛 𝜕𝛿 𝑛 𝑘 𝜕𝑃2 𝜕 𝑉2 𝑘 ⋯ 𝜕𝑃2 𝜕 𝑉𝑛 𝑘 ⋮ ⋱ ⋮ 𝜕𝑃2 𝜕 𝑉2 𝑘 ⋯ 𝜕𝑃 𝑛 𝜕 𝑉𝑛 𝑘 𝜕𝑄2 𝜕𝛿2 𝑘 ⋯ 𝜕𝑄2 𝜕𝛿 𝑛 𝑘 ⋮ ⋱ ⋮ 𝜕𝑄 𝑛 𝜕𝛿2 𝑘 ⋯ 𝜕𝑄 𝑛 𝜕𝛿 𝑛 𝑘 𝜕𝑄2 𝜕 𝑉2 𝑘 ⋯ 𝜕𝑄2 𝜕 𝑉2 𝑘 ⋮ ⋱ ⋮ 𝜕𝑄 𝑛 𝜕 𝑉2 𝑘 ⋯ 𝜕𝑄2 𝜕 𝑉𝑛 𝑘 ∆𝛿2 ⋮ 𝑘 ∆𝛿 𝑛 𝑘 ∆|𝑉2 𝑘 | ⋮ ∆|𝑉𝑛 𝑘 | (A7) In above equation, bus 1 is assumed to be slack bus. The jacobian matrix gives the linearized relationship between small changes in voltage angle ∆𝛿𝑖 (𝑘) and voltage magnitude Δ|𝑉𝑖 (𝑘) | with small changes in real and reactive power Δ𝑃𝑖 (𝑘) and Q i (k) elements of jacobian matrix are the partial derivatives of (5) and (6) evaluated at ∆𝛿𝑖 (𝑘) and Δ|𝑉𝑖 (𝑘) |. 𝛥𝑃 𝛥𝑄 = 𝐽1 𝐽2 𝐽3 𝐽4 𝛥𝛿 𝛥|𝑉| (A8) Accordingly there are (n-1) real power constraints and (n-1-m) reactive power constraints and the jacobian matrix is the order of (2n-2-m) (2n-2-m). J1 is the order of (n-1) x (n-1) 𝜕𝑃 𝑖 𝜕𝛿 𝑖 = 𝑉𝑖 𝑉𝑗 𝑌𝑖𝑗 𝑛 𝑗 ≠1 sin⁡(𝜃𝑖𝑗 – 𝛿𝑖+ 𝛿𝑗 ) (A9) 𝜕𝑃 𝑖 𝜕𝛿 𝑗 = − 𝑉𝑖 Vj Yij sin 𝜃𝑖𝑗 – 𝛿𝑖+ 𝛿𝑗 j≠ 1 (A10) J2 is the order of (n-1) x (n-1-m) 𝜕𝑃 𝑖 𝜕|𝑉 𝑖| = 2 𝑉𝑖 𝑌𝑖𝑖 𝑐𝑜𝑠𝜃𝑖𝑖 + 𝑉𝑖 𝑉𝑗 𝑌𝑖𝑗 𝑛 𝑗≠1 cos⁡(𝜃𝑖𝑗 – 𝛿𝑖+ 𝛿𝑗 ) (A11) 𝜕𝑃 𝑖 𝜕|𝑉 𝑗 | = 𝑉𝑖 𝑌𝑖𝑗 𝑐 os 𝜃𝑖𝑗 – 𝛿𝑖+ 𝛿𝑗 j≠ 𝑖 (A12) J3 is the order of (n-1-m) x (n-1) 𝜕𝑄 𝑖 𝜕𝛿 𝑖 = 𝑉𝑖 𝑉𝑗 𝑌𝑖𝑗 𝑛 𝑗 ≠1 cos 𝜃𝑖𝑗 – 𝛿𝑖+ 𝛿𝑗 (A13) J4 is the order of (n-1-m) x (n-1-m) 𝜕𝑄 𝑖 𝜕𝛿 𝑗 = − 𝑉𝑖 Vj Yij cos⁡(𝜃𝑖𝑗 – 𝛿𝑖+ 𝛿𝑗 )j≠ 𝑖 (A14) 𝜕𝑄 𝑖 𝜕|𝑉 𝑖| = 𝑉𝑖 𝑌𝑖𝑖 𝑠𝑖𝑛𝜃𝑖𝑖 − 𝑉𝑗 𝑌𝑖𝑗 𝑛 𝑗≠1 sin 𝜃𝑖𝑗 – 𝛿𝑖+ 𝛿𝑗 (A15) 𝜕𝑄 𝑖 𝜕|𝑉 𝑗 | = − 𝑉𝑖 𝑌𝑖𝑗 𝑠𝑖𝑛⁡(𝜃𝑖𝑗 – 𝛿𝑖+ 𝛿𝑗 )j≠ 𝑖 (A16) The terms Δ𝑃𝑖 (𝑘) and Q i (k) are difference between the schedule and calculated values, known as the power residuals, given by Δ𝑃𝑖 𝑘 = Pi sch − Pi k (A17) Δ𝑄𝑖 𝑘 = Qi sch − Qi k (A18) The new estimated for bus voltage i 𝛿𝑖 (𝑘+1) = 𝛿𝑖 (𝑘) + ∆𝛿𝑖 (𝑘) (A19)
  • 3. IJRET: International Journal of Research in Engineering and Technology eISSN: 2319-1163 | pISSN: 2321-7308 __________________________________________________________________________________________ Volume: 03 Special Issue: 07 | May-2014, Available @ http://www.ijret.org 322 |𝑉𝑖 (𝑘+1) | = |𝑉𝑖 𝑘 | + ∆|𝑉𝑖 (𝑘) | (A20) 2.2 Security Constrained Optimal Power Flow Objective function Min f ( p ) (A21) Subject to g (p) = 0 hmin ≤ h(p) ≤ hmax in normal condition constraints h’min ≤ h’ (p) ≤ h’max in contingency condition constraints. Security constrained optimal power flow solution (SCOPF) will always a cost ≥ optimal power flow cost. If we ignore losses, then we can say that an OPF solution differs from an EDC solution only when a normal transmission constraint. When normal flow moves from just < 100% to ≥ 100% of continuous rating SCOPF differs from an OPF solution only when contingency constraint becomes binding occurs when post- contingency flow moves from just < 100% to ≥ 100% of emergency rating. Now let’s consider the SCOPF. Its problem statement is given as problem Pp : Min f ( x0, u0 ) (A22) gk ( xk , uk ) =0 k=0,1,2….c hk ( xk, uk ) = hk max k=0,1,2,….c Notice that there are C contingencies to be addressed in the SCOPF and that there are a complete new set of constraints for each of these C contingencies observe. Each of contingency related equality constraints is exactly like the original set of equality constraints except if corresponds the system with an element removed. Each set of contingency related inequality-constraints is exactly like the original set of inequality constraints except its corresponds to the system with an element removed and branch flow constraints and for voltage magnitudes, the limits will be different. Also notice that the constraints are a function of xk , the voltage magnitudes and angles under the pre-contingency (k=0) and contingency condition ( k > 1,2,….c) and u0. 2.3 Locational Marginal Pricing Calculation Locational marginal pricing (LMPs) are determined from the result of a security-constrained least-cost dispatch. It is a taxi ride for MW. It may differs in the various location(bus). We need two factors to deside the locational marginal pricing. (i)Transmission congestion (ii)Losses The Locational Marginal Pricing (nodel price) at bus i can be calculated using the following equation LMP = marginalcost + congestioncost + lossescost λi = λRef + λCongest + λLossi (A23) λi = λRef - Li x λRef - 𝜇𝑗 𝑋 𝑆𝐹𝑗𝑖𝑗 (A24) λLossi = (- Li x λRef) - losses from the reference bus to bus i = (+ Li x λRef ) – losses from bus i to reference bus λCongest = (- 𝜇𝑗 𝑋 𝑆𝐹𝑗𝑖𝑗 - congestion from reference bus to bus i = (+ 𝜇𝑗 𝑋 𝑆𝐹𝑗𝑖𝑗 - congestion from bus I to reference bus. 3. LMP ALGORITHM IN POWER WORLD SIMULATOR 3.1. Locational Marginal Pricing Algorithm Step 1: Draw the simulink one line diagram in new case window of power world simulator for the given power system in edit mode. Step 2: Set the cubic cost of each generation and to convert piece wise linear cost. Step 3: Save the case with apt name. Step 4: Select tools in run mode and to solve the power flow by using full N-R method. Step 5: Open Add-ons Ribbon Tab Step 6: To select the OPF case information of the Dialog box and to select the all OPF area records. a.) If OPF records “YES” that record is included in the marginal cost calculation. b.) If OPF records “NO” that records is not include in the marginal cost calculation. Step 7: To set all the OPF constraints and also include common constraints. Step 8: Open the SCOPF dialog box in the add-ons ribbon tab. Step 9: Run Full security constrained OPF under normal condition (zero contingency in contingency analysis form) Step10: To calculate marginal cost of each bus (location) before contingency. Step11: To view contingency analysis form in the SCOPF dialog box. Step12: Right click on label and select auto insert contingencies through insert special option. Step 13: Verify that single transmission line or transformer is selected. Step 14: If can limit the contingencies inserted to only those meeting define filter. Step 15: we want to insert contingencies for all branches and generators so no filtering is desired. Step 16: To check the following conditions a.) Remove the checkmarks in use area/zone filters.
  • 4. IJRET: International Journal of Research in Engineering and Technology eISSN: 2319-1163 | pISSN: 2321-7308 __________________________________________________________________________________________ Volume: 03 Special Issue: 07 | May-2014, Available @ http://www.ijret.org 323 b.) Verify no other options are selected. Step 17: Click Do insert contingencies button to accept the all contingencies. Step 18: Click “ YES” to get the contingencies. Step 19: Now the contingency analysis dialog shows contingencies a.) Right click on the list display on the contingency tap and select insert special and click auto insert to the local menu. b.) Select single generating unit then click the do insert contingency button. Click “YES” to complete. Step 20: Click “Start Run” on the contingency tab and close the tab. Step 21: The contingency elements are include in the SCOPF dialog box. Step 21: To Run the Full Security Constrained optimal power Flow. Step 22: We get the marginal cost under contingency condition. In addition congestion, transmission loss cost in each location (bus) 4. TEST SYSTEM The Locational marginal pricing of 14-bus test system is shown below when the power flow is running on the power world simulator. The percentage of power flow is mentioned in power flow diagram. It consists of five generators for dispatch of power. Fig. 1 IEEE-14 bus system 5. SIMULATION RESULTS Locational marginal pricing (LMP) is important in deregulated power market under normal and contingency condition. So we calculate the LMP by considering without contingency, line outage, generation outage and also include all contingency. First we find the contingency elements in IEEE-14 bus system. Table.1 Single contingency of line & generator Label Violations Max Branch % Min Voltage L_0000011- 0000022C1 2 276.5 - L_0000011- 0000055C1 1 127.5 - L_0000022- 0000033C1 2 102.0 - L_0000022- 0000044C1 1 102.0 - L_0000022- 0000055C1 1 103.0 - L_0000066- 00001313C1 2 - 0.898 L_0000099- 00001414C1 1 - 0.848 G_0000022U 1 1 103.6 - Table.2 Multiple contingency of both line and Generator Label Violations Max Branch % Min Voltage G_0000022u1& L_0000022- 0000033c1 3 115.7 - G_0000022u1& L_0000066- 00001313c1 3 104.6 0.894 G_0000033u1& L_0000011- 0000022c1 4 325.3 0.891 G_0000033u1& L_0000011- 0000055c1 3 140.2 - G_0000033u1& L_0000022- 0000033c1 6 118.0 0.739 G_0000066u1& L_0000011- 0000022c1 3 283.4 0.883 G_0000066u1& L_0000077- 5 - 0.842
  • 5. IJRET: International Journal of Research in Engineering and Technology eISSN: 2319-1163 | pISSN: 2321-7308 __________________________________________________________________________________________ Volume: 03 Special Issue: 07 | May-2014, Available @ http://www.ijret.org 324 0000099c1 G_0000066u1& L_0000099- 000001414c1 3 - 0.733 Table 3 LMP in Normal condition Locational marginal pricing of each bus system is in the above table. Ther is no congestion cost due to normal power flow condition. Table 4 LMP in Contingency Condition (only includes single contingency). Bus No. Area name Energy cost ($/mwh) Conge cost ($/mw h) Loss cost ($/mwh ) LMP ($/mwh) 1 Top 5.86 0.00 0.00 5.86 2 Top 5.86 5.88 0.26 12.00 3 Top 5.86 593.05 0.78 599.69 4 Top 5.86 530.35 0.74 536.95 5 Top 5.86 389.45 0.64 395.94 6 Top 5.86 432.90 0.64 439.40 7 Top 5.86 510.31 0.75 516.92 8 Top 5.86 510.31 0.75 516.92 9 Top 5.86 499.13 0.76 505.74 10 Top 5.86 490.63 0.78 497.27 11 Top 5.86 463.86 0.73 470.46 12 Top 5.86 446.39 0.78 453.03 13 Top 5.86 456.67 0.86 463.40 14 Top 5.86 509.62 1.21 516.68 The congestion cost is occurred in each bus due to single contingency of line outage and generation outage of each elements Table 5 LMP price (includes all contingencies) Bus No. Area name Energy cost ($/mwh) Conge Cost ($/mwh) Losses ($/mwh) LMP ($/mwh) 1 Top -1161.26 1167.12 0.00 5.86 2 Top -1161.26 1223.32 -50.06 12.00 3 Top -1161.26 1610.35 -152.80 296.29 4 Top -1161.26 3318.47 -144.14 2013.07 5 Top -1161.26 1123.77 -124.52 -162.01 6 Top -1161.26 14055.13 -124.74 12769.1 3 7 Top -1161.26 3572.12 -146.98 2263.83 8 Top -1161.26 3572.12 -146.98 2263.88 9 Top -1161.26 5166.58 -148.48 3856.84 10 Top -1161.26 13749.70 -152.99 12435.4 5 11 Top -1161.26 13944.48 -143.86 12639.3 6 12 Top -1161.26 14223.56 -152.48 12909.8 2 13 Top -1161.26 14038.94 -169.60 12708.0 7 14 Top -1161.26 12030.02 -237.59 10631.1 7 In security constrained condition all single and multiple contingency elements are include in the power system. The LMP value is maximum due to congestion cost. Fig. 2 LMP Cost Curve In Each Bus (before contingencies) Bus No. Area name Energy cost ($/mwh) Cong cost ($/m wh) Loss cost ($/mwh) LMP ($/mwh) 1 Top 11.28 -0.00 0.00 11.28 2 Top 11.28 -0.00 0.72 12.00 3 Top 11.28 -0.00 1.72 12.99 4 Top 11.28 -0.00 1.60 12.87 5 Top 11.28 -0.00 1.39 12.66 6 Top 11.28 -0.00 1.39 12.67 7 Top 11.28 -0.00 1.63 12.90 8 Top 11.28 -0.00 1.63 12.90 9 Top 11.28 -0.00 1.64 12.92 10 Top 11.28 -0.00 1.68 12.96 11 Top 11.28 -0.00 1.59 12.86 12 Top 11.28 -0.00 1.67 12.94 13 Top 11.28 -0.00 1.84 13.11 14 Top 11.28 -0.00 2.51 13.79
  • 6. IJRET: International Journal of Research in Engineering and Technology eISSN: 2319-1163 | pISSN: 2321-7308 __________________________________________________________________________________________ Volume: 03 Special Issue: 07 | May-2014, Available @ http://www.ijret.org 325 Fig.3 Lmp Cost Curve in Each Bus (only includes line outage and generator outage) Fig .4 Lmp Cost Curve In Each Bus (includes all contingencies) 6. CONCLUSIONS The state estimation based contingency analysis technique is more complex. But Newton’s method based contingency analysis and security constrained optimal power flow in contingency condition is simple and more accuracy in power world simulator. The LMP calculation of test system is easy in power world simulation. In addition LMP is depends on congestion cost of each bus system. Whenever the contingency is occur in this system that time only congestion cost is added in LMP. Otherwise there is no congestion cost in the system. So the locational marginal pricing is reduced in the normal condition. The LMP calculation is helpful for demand response improvement under security condition in deregulated power market. In future we can easily reduce the cost of the power transmission and to improve the demand response in secured manner by using the reserve option and to connect from renewable generation to the transmission network. And we get minimum congestion cost in each location. So we can easily improve the demand response in the secured dispatch scheduling in both normal and the contingency condition in the power world simulator tool. It is user friendly software in the pricing calculation of deregulated power market in the power system network. REFERENCES [1]. Amit Kumar Roy “contingency analysis in power system” Thapur University Patiala.2009. [2]. R.Manikandan, M.Bhoopathi “ Contingency analysis in deregulated power market” Jayaram College of Engineering and Technology, Tamilnadu 2013. [3]. Richard D.Christie, AnjanBose “load frequency control issues in power system operations after deregulation” University of Washington, 1995. [4]. Saavedra, O.R., "Solving the security constrained optimal power flow problem in a distributed computing environment," Generation, transmission and distribution, IEEE proceedings, vol.143, No.6 pp 593-598, 1996. [5]. Fangxing Li, Rui Bo, “ Congestion and Price Prediction Under Load Variation,”IEEE Transaction on power system , Vol 24,No.2 May 2009. BIOGRAPHIES Manikandan R obtained his Bachelor degree in Electrical & Electronics Engineering from P.R.Engineering College, Thanjavur in the year 2012 and Master degree in Power Systems Engineering doing from Jayaram College of Engineering and technology, Thuraiyur, India. His research area includes deregulation of power market and Smart grid technologies.
  • 7. IJRET: International Journal of Research in Engineering and Technology eISSN: 2319-1163 | pISSN: 2321-7308 __________________________________________________________________________________________ Volume: 03 Special Issue: 07 | May-2014, Available @ http://www.ijret.org 326 Bhoopathi M received his B.E. degree in Electrical and Electronics Engineering from Kumaraguru College of technology, Bharathiyar University, Coimbatore, India and M.E. degree in Power systems engineering from Annamalai University, Chidambaram, India. He is currently working as an Assistant Professor in Department of Electrical and Electronics Engineering, Jayaram college of engineering and technology, Pagalavadi, Thuraiyur , India. His research interest includes Restructured power system and smart grid technologies. Saravanakumar R received his B.E. degree in Electrical and Electronics Engineering from PSNA College of engineering & technology, Dindigul, Anna University, Chennai, India and M.E. degree in Power systems engineering from University College of engineering, BIT campus, Trichirappalli, India. He is currently working as an Assistant Professor in Department of Electrical and Electronics Engineering,Sudharsan engineering college, Sathiyamangalam, India. His research interest includes power systems and Distributed Generation.