Power Systems Laboratory: Submitted by
Power Systems Laboratory: Submitted by
POWER
SYSTEMS
LABORATORY
REPORT
Submitted by:
RAJESH ROY
17-13-018
6th Semester
EE Dept
NIT SILCHAR
2|Page POWER SYSTEM LABORATORY
CONTENTS
1. (a) Solution of Economic Load Dispatch problem neglecting both
transmission losses and generating limits. ------------------------ (3)
Where F t=Total Fuel cost of all the generating units of the system.
F n= Fuel cost of the ‘nth’ generating units.
P D= Total demand of the System load.
Pn=Real power generating of the ‘nth’ generating unit.
Therefore the condition for optimum operating is,
(IC)1= (IC)2=…………………………….=(IC)n=Incremental Fuel Cost
PROCEDURE:
Step 1: Read the generating limits, Real power demand i.e. PD cost coefficients of
all the generating units present in the power system.
Step 2: The value of the incremental fuel cost (λ) & the change in the value of the
incremental fuel cost i.e. Δλ is assumed in the range of 0.01.
Step 3: Read the Convergence criterion i.e Ɛ
Step 4: Calculate the real power generation of the generators for the system.
4|Page POWER SYSTEM LABORATORY
Step 7: Check ΔP<Ɛ (Convergence Criterion) and if satisfies then ‘stop else’, then
change the value of λ i.e. increase or decrease the incremental fuel cost by Δλ.
Step 8: Go to Step 4.
SYSTEM INFORMATION:
Generator Data:
PD = 800 MW
Generator No α Β γ
1 510 7.2 0.00142
2 310 7.8 0.00194
5
3 18 7.9 0.00482
7
FLOWCHART: START
START
SET n=1
No SET n = n+1
Yes
CALCULATE ΔP = |∑Pn – PD|
PRINT
GENERATION AND
ITS COST
5|Page POWER SYSTEM LABORATORY
IS ΔP> Ɛ Yes
No
END
ADJUST THE VALUE OF λ AND
CODE: EVALUATE λ = λ + Δλ
clc;
clear all;
close all;
Pd=800;
alpha=[510 310 18];
beta=[7.2 7.85 7.97];
gamma=[0.00142 0.00194 0.00482];
lambda=12;
deltaP=1;
iteration=0;
while(abs(deltaP)>0.001)
iteration=iteration+1
for i=1:3
P(i)=(lambda-beta(i))/(2*gamma(i))
end
Pg=sum(P)
deltaP=(Pd-Pg);
j=0;
for i=1:3
j=j+(1/(2*gamma(i)));
end
deltaLambda=deltaP/j;
lambda=lambda+(deltaLambda/2)
for i=1:3
FC(i)=alpha(i)+(beta(i)*P(i))+(gamma(i)*P(i)*P(i));
end
TFC(iteration)=sum(FC);
end
fprintf('The Generated Power are: ');
display(P);
fprintf('The value of Lambda is: ');
display(lambda);
fprintf('The number of Iterations is: ');
display(iteration);
fprintf('The Total Cost is: ');
display(TFC);
plot(TFC);
title('ECONOMIC LOAD DISPATCH NEGLECTING BOTH TRANSMISSION LOSSES
AND GENERATING LIMITS');
6|Page POWER SYSTEM LABORATORY
xlabel('ITERATIONS');
ylabel('TOTAL FUEL COST');
grid on
RESULTS:
The values of Lambda, P1, P2, P3 and total fuel cost throughout the iterations are:
GRAPH:
DISCUSSIONS:
The iterative method is used for solving economic load dispatch problems. In this
case, we considered no generating limits or transmission losses. The iteration
terminated after 23 counts for the initial assumption. The lambda after the iteration
came out to be Rs. 8.6678/MWh. The power generation for unit 1 is 516.8362
MW, for unit 2 is 210.7770 MW and for unit 3 is 72.3874 MW. The minimum fuel
cost came out to be Rs. 7281.50 per hour.
8|Page POWER SYSTEM LABORATORY
Min F t=∑ Fn
n =1
N
PROCEDURE:
Step 1: Read the generating limits, Real power demand i.ePD, cost coefficients of
all the generating units present in the power system.
9|Page POWER SYSTEM LABORATORY
Step 2: The value of the incremental fuel cost (λ) & the change in the value of the
incremental fuel cost i.e. Δλ is assumed in the range of 0.01.
Step 3: Read the Convergence criterion i.e Ɛ.
Step 4: Calculate the real power generation of the generators of the system.
Step 5: Check inequalities and take corrective actions.
N N
SYSTEM INFORMATION:
Generator Data:
PD = 800 MW
Generator No Pmax(MW) Pmin(MW) α β γ
1 500 150 510 7.2 0.00142
2 400 100 310 7.8 0.00194
5
3 200 50 18 7.9 0.00482
7
10 | P a g e POWER SYSTEM LABORATORY
START
FLOWCHART:
SET n=1
No Yes
CHECK IF 𝑃n > 𝑃n, max SET 𝑃n=𝑃n, max
No
CHECK IF 𝑃n < 𝑃n, min Yes SET 𝑃n=𝑃n, min
No
SET n = n+1
Yes
CODE:
clc;
clear all;
close all;
Pd=800;
alpha=[510 310 18];
beta=[7.2 7.85 7.97];
gamma=[0.00142 0.00194 0.00482];
Pmax=[500 400 200];
Pmin=[150 100 50];
lambda=10;
deltaP=1;
iteration=0;
while abs(deltaP)>0.01
iteration=iteration+1;
for i=1:3
P(i)=(lambda-beta(i))/(2*gamma(i))
if(P(i)>Pmax(i))
P(i)=Pmax(i);
end
if(P(i)<Pmin(i))
P(i)=Pmin(i);
end
end
Pg=sum(P);
deltaP=Pd-Pg;
k=0;
for j=1:3
k=k+1/(2*gamma(j));
end
deltaLambda=deltaP/k;
lambda=lambda+(deltaLambda/2)
for i=1:3
FC(i)=alpha(i)+beta(i)*P(i)+gamma(i)*P(i)*P(i);
end
TFC(iteration)=sum(FC);
end
fprintf('The Generated Power are: ');
display(P);
fprintf('The value of Lambda is: ');
display(lambda);
fprintf('The number of Iterations is: ');
display(iteration);
fprintf('The Total Cost is: ');
display(TFC);
plot(TFC);
title('ECONOMIC LOAD DISPATCH CONSIDERING GENERATING LIMITS AND NEGLECTING
BOTH TRANSMISSION LOSSES ');
xlabel('ITERATIONS');
ylabel('TOTAL FUEL COST');
grid on
12 | P a g e POWER SYSTEM LABORATORY
RESULTS:
The values of Lambda, P1, P2, P3 and total fuel cost throughout the iterations are:
GRAPH:
DISCUSSIONS:
The iterative method is used for solving economic load dispatch problems. In this
case, we considered generating limits but no transmission losses. The iteration
terminated after 39 counts for the initial assumption. The lambda after the iteration
came out to be Rs. 8.7144/MWh. The power generation for unit 1 is 500 MW, for
unit 2 is 222.7875 MW and for unit 3 is 77.2215 MW. The minimum fuel cost
came out to be Rs. 7282.37 per hour.
14 | P a g e POWER SYSTEM LABORATORY
EXPERIMENT NO.2
TITLE: Solution of Economic Load Dispatch problem considering both
transmission losses and generating limits.
Min F t=∑ Fn
n =1
N
PROCEDURE:
Step 1: Read the generating limits, Cost coefficients of all generating units and
Real power demand i.e PD and system loss coefficient of all the generating units
present in the power system.
Step 2: The value of the incremental fuel cost (λ) & the change in the value of the
incremental fuel cost i.eΔλ is assumed in the range of 0.01.
Step 3: Read the Convergence criterion i.e Ɛ.
Step 4: Calculate the real power generation of the generators of the system.
Step 5: Calculate the real power loss (PL).
Step 6: Check inequalities and take corrective actions.
N N
Step 7: Calculate ∑ Pn=PD-PL and ∑ Pn - PD- PL = ΔP.
n =1 n =1
SYSTEM INFORMATION:
Generator Data:
PD = 800 MW
Generator No Pmax(MW) Pmin(MW) α β γ
1 600 150 510 7.2 0.00142
2 400 100 310 7.8 0.00194
5
3 200 50 18 7.9 0.00482
7
[
B= 0.0093 0.0228 0.0017
0.0028 0.0031 0.0015 ]
B0=[ 0.0003 0.0031 0.0015 ]
B00 = 0.00030523
16 | P a g e POWER SYSTEM LABORATORY
START
FLOWCHART:
DETERMINE 𝑃n CORRESPONDING TO λ
SET K=0
SET n=1
Solve for 𝑃n
𝑃n =
fn
1− −2 Σ m ≠n B mn Pm
λ
Fnn
+2B
CHECKλ IF 𝑃n >nn 𝑃n, max
Yes
SET 𝑃n=𝑃n, max
No
CHECK IF 𝑃n < 𝑃n, min Yes SET 𝑃n=𝑃n, min
No
n = n+1
No Yes
IS | 𝑃ik - 𝑃ik-1 | < Ɛ
17 | P a g e POWER SYSTEM LABORATORY
K=K+1
CODE:
clc;
clear all;
close all;
alpha=[510 310 18];
beta=[7.2 7.85 7.97];
gamma=[0.00142 0.00194 0.00482];
Pmax=[500 400 200];
Pmin=[150 100 50];
B = [0.0218 0.0093 0.0028;
0.0093 0.0228 0.0017;
0.0028 0.0031 0.0015] * 0.01;
B0 = [0.0003 0.0031 0.0015];
B00 = 0.00030523;
lambda=14;
Pd=800;
P = [516.8359 210.7768 72.3873];
deltaP=1;
iteration=0;
while abs(deltaP)>0.00001
iteration=iteration+1;
for i = 1:3
BijPj = 0;
18 | P a g e POWER SYSTEM LABORATORY
for j = 1:3
if j ~= i
BijPj = BijPj + B(i,j).*P(j);
end
end
P(i) = (lambda - 2.*lambda.*BijPj - lambda.*B0(i) - beta(i)) /
(2.*(gamma(i) + lambda.*B(i, i)))
end
for i = 1:3
if (P(i)>Pmax(i))
P(i) = Pmax(i);
end
if (P(i)<Pmin(i))
P(i) = Pmin(i);
end
end
for i = 1:3
FC(i) = alpha(i) + beta(i)*P(i) + gamma(i)*P(i)*P(i);
end
P_Tran = transpose(P);
Pl = P*B*P_Tran + B0*P_Tran + B00;
Pg = sum(P);
deltaP = abs(Pg - Pd - Pl);
j = 0;
for i = 1:3
j = j + (gamma(i) + B(i,i)*beta(i)) / (2*(gamma(i) +
lambda*B(i,i))^2);
end
deltaLambda = deltaP/j;
if(Pg>Pd + Pl)
lambda = lambda - deltaLambda
end
if(Pg<Pd + Pl)
lambda = lambda + deltaLambda
end
for i=1:3
FC(i)=alpha(i)+beta(i)*P(i)+gamma(i)*P(i)*P(i);
end
TFC(iteration) = sum(FC);
end
fprintf('The value of Lambda is: ');
disp(lambda);
fprintf('The Generated Power are: ');
disp(P);
fprintf('The number of Iterations is: ');
disp(iteration);
fprintf('The Power Lost is: ');
disp(Pl);
fprintf('The Total Generated power is: ');
disp(Pg);
fprintf('The Total costs is: ')
disp(TFC);
plot(TFC);
title('ECONOMIC LOAD DISPATCH CONSIDERING BOTH GENERATING LIMITS AND
TRANSMISSION LOSSES ');
xlabel('ITERATIONS');
ylabel('TOTAL FUEL COST');
19 | P a g e POWER SYSTEM LABORATORY
grid on
RESULTS:
The values of Lambda, P1, P2, P3 and total fuel cost throughout the iterations are:
GRAPH:
DISCUSSIONS:
The iterative method is used for solving economic load dispatch problems. In this
case, we considered both generating limits and transmission losses. The iteration
terminated after 40 counts for the initial assumption. The lambda after the iteration
came out to be Rs. 11.1522/MWh. The power generation for unit 1 is 427.2084
MW, for unit 2 is 257.1718 MW and for unit 3 is 200 MW. The power lost in
transmission is 84.3802 MW. The minimum fuel cost came out to be Rs. 8106.96
per hour.