Hi,
I'm trying to run a simple PF with the newton method, i attached my case
files with this message.
And i also attached 2 pictures of the network cases.
I'm working with the case30.m that i modify, i just wanted to use the
line specifications.
The network now have just 16 bus, 17 branch and 1 Generator(also slack
bus), and the network is radial.
1st problem: The powerflow for this network ("case30_1dep.m") doesn't
converge and i can't understand why, how can it be explain?
But it run normally if i add a generator ("case30_1dep_ok.m") that i
take from the original "case30.m".
So i thought that maybe the slack bus genreation limit is reached but
even with modyfing the slack bus specification in the "case30_1dep.m",
it's still doesn't convergent.
2nd question: What kind of physical specification (V, P and Q limits
for example) can affect the convergence.
Thank you in advance for the time you will take to read my questions
Best regards
Abdelkrim Ali Zazou
.
The problem is that this loadflow don't converge and i can't understand
whey,
can you give me some explaination about this issue.
The problem
function mpc = case30
%CASE30 Power flow data for 30 bus, 6 generator case.
% Please see CASEFORMAT for details on the case file format.
%
% Based on data from ...
% Alsac, O. & Stott, B., "Optimal Load Flow with Steady State Security",
% IEEE Transactions on Power Apparatus and Systems, Vol. PAS 93, No. 3,
% 1974, pp. 745-751.
% ... with branch parameters rounded to nearest 0.01, shunt values divided
% by 100 and shunt on bus 10 moved to bus 5, load at bus 5 zeroed out.
% Generator locations, costs and limits and bus areas were taken from ...
% Ferrero, R.W., Shahidehpour, S.M., Ramesh, V.C., "Transaction analysis
% in deregulated power systems using game theory", IEEE Transactions on
% Power Systems, Vol. 12, No. 3, Aug 1997, pp. 1340-1347.
% Generator Q limits were derived from Alsac & Stott, using their Pmax
% capacities. V limits and line |S| limits taken from Alsac & Stott.
% MATPOWER
% $Id: case30.m,v 1.12 2010/03/10 18:08:13 ray Exp $
%% MATPOWER Case Format : Version 2
mpc.version = '2';
%%----- Power Flow Data -----%%
%% system MVA base
mpc.baseMVA = 100;
%% bus data
% bus_i type Pd Qd Gs Bs area Vm Va
baseKV zone Vmax Vmin
mpc.bus = [
3 1 2.4 1.2 0 0 1 1 0
135 1 1.05 0.95;
4 1 7.6 1.6 0 0 1 1 0
135 1 1.05 0.95;
6 1 0 0 0 0 1 1 0
135 1 1.05 0.95;
7 1 22.8 10.9 0 0 1 1 0
135 1 1.05 0.95;
8 1 30 30 0 0 1 1 0
135 1 1.05 0.95;
9 1 0 0 0 0 1 1 0
135 1 1.05 0.95;
10 1 5.8 2 0 0 3 1 0
135 1 1.05 0.95;
11 1 0 0 0 0 1 1 0
135 1 1.05 0.95;
12 1 11.2 7.5 0 0 2 1 0
135 1 1.05 0.95;
13 3 0 0 0 0 2 1 0
135 1 1.5 0.5;
14 1 6.2 1.6 0 0 2 1 0
135 1 1.05 0.95;
15 1 8.2 2.5 0 0 2 1 0
135 1 1.05 0.95;
16 1 3.5 1.8 0 0 2 1 0
135 1 1.05 0.95;
17 1 9 5.8 0 0 2 1 0
135 1 1.05 0.95;
18 1 3.2 0.9 0 0 2 1 0
135 1 1.05 0.95;
19 1 9.5 3.4 0 0 2 1 0
135 1 1.05 0.95;
20 1 2.2 0.7 0 0 2 1 0
135 1 1.05 0.95;
28 1 0 0 0 0 1 1 0
135 1 1.05 0.95;
];
%% generator data
% bus Pg Qg Qmax Qmin Vg mBase status Pmax
Pmin Pc1 Pc2 Qc1min Qc1max Qc2min Qc2max ramp_agc ramp_10
ramp_30 ramp_q apf
mpc.gen = [
13 70 0 300 300 2 300 1 300
0 0 0 0 0 0 0 0 0 0
0 0;
];
%% branch data
% fbus tbus r x b rateA rateB rateC ratio
angle status angmin angmax
mpc.branch = [
3 4 0.01 0.04 0 130 130 130 0
0 1 -360 360;
4 6 0.01 0.04 0 90 90 90 0
0 1 -360 360;
6 7 0.03 0.08 0.01 130 130 130 0
0 1 -360 360;
6 8 0.01 0.04 0 32 32 32 0
0 1 -360 360;
6 10 0 0.56 0 32 32 32 0
0 1 -360 360;
9 11 0 0.21 0 65 65 65 0
0 1 -360 360;
9 10 0 0.11 0 65 65 65 0
0 1 -360 360;
4 12 0 0.26 0 65 65 65 0
0 1 -360 360;
12 13 0 0.14 0 65 65 65 0
0 1 -360 360;
12 14 0.12 0.26 0 32 32 32 0
0 1 -360 360;
12 15 0.07 0.13 0 32 32 32 0
0 1 -360 360;
12 16 0.09 0.2 0 32 32 32 0
0 1 -360 360;
15 18 0.11 0.22 0 16 16 16 0
0 1 -360 360;
19 20 0.03 0.07 0 32 32 32 0
0 1 -360 360;
10 20 0.09 0.21 0 32 32 32 0
0 1 -360 360;
10 17 0.03 0.08 0 32 32 32 0
0 1 -360 360;
6 28 0.02 0.06 0.01 32 32 32 0
0 1 -360 360;
];
function mpc = case30
%CASE30 Power flow data for 30 bus, 6 generator case.
% Please see CASEFORMAT for details on the case file format.
%
% Based on data from ...
% Alsac, O. & Stott, B., "Optimal Load Flow with Steady State Security",
% IEEE Transactions on Power Apparatus and Systems, Vol. PAS 93, No. 3,
% 1974, pp. 745-751.
% ... with branch parameters rounded to nearest 0.01, shunt values divided
% by 100 and shunt on bus 10 moved to bus 5, load at bus 5 zeroed out.
% Generator locations, costs and limits and bus areas were taken from ...
% Ferrero, R.W., Shahidehpour, S.M., Ramesh, V.C., "Transaction analysis
% in deregulated power systems using game theory", IEEE Transactions on
% Power Systems, Vol. 12, No. 3, Aug 1997, pp. 1340-1347.
% Generator Q limits were derived from Alsac & Stott, using their Pmax
% capacities. V limits and line |S| limits taken from Alsac & Stott.
% MATPOWER
% $Id: case30.m,v 1.12 2010/03/10 18:08:13 ray Exp $
%% MATPOWER Case Format : Version 2
mpc.version = '2';
%%----- Power Flow Data -----%%
%% system MVA base
mpc.baseMVA = 100;
%% bus data
% bus_i type Pd Qd Gs Bs area Vm Va
baseKV zone Vmax Vmin
mpc.bus = [
3 1 2.4 1.2 0 0 1 1 0
135 1 1.05 0.95;
4 1 7.6 1.6 0 0 1 1 0
135 1 1.05 0.95;
6 1 0 0 0 0 1 1 0
135 1 1.05 0.95;
7 1 22.8 10.9 0 0 1 1 0
135 1 1.05 0.95;
8 1 30 30 0 0 1 1 0
135 1 1.05 0.95;
9 1 0 0 0 0 1 1 0
135 1 1.05 0.95;
10 1 5.8 2 0 0 3 1 0
135 1 1.05 0.95;
11 1 0 0 0 0 1 1 0
135 1 1.05 0.95;
12 1 11.2 7.5 0 0 2 1 0
135 1 1.05 0.95;
13 3 0 0 0 0 2 1 0
135 1 1.5 0.5;
14 1 6.2 1.6 0 0 2 1 0
135 1 1.05 0.95;
15 1 8.2 2.5 0 0 2 1 0
135 1 1.05 0.95;
16 1 3.5 1.8 0 0 2 1 0
135 1 1.05 0.95;
17 1 9 5.8 0 0 2 1 0
135 1 1.05 0.95;
18 1 3.2 0.9 0 0 2 1 0
135 1 1.05 0.95;
19 1 9.5 3.4 0 0 2 1 0
135 1 1.05 0.95;
20 1 2.2 0.7 0 0 2 1 0
135 1 1.05 0.95;
28 1 0 0 0 0 1 1 0
135 1 1.05 0.95;
22 3 0 0 0 0 3 1 0 135
1 1.1 0.95;
];
%% generator data
% bus Pg Qg Qmax Qmin Vg mBase status Pmax
Pmin Pc1 Pc2 Qc1min Qc1max Qc2min Qc2max ramp_agc ramp_10
ramp_30 ramp_q apf
mpc.gen = [
13 70 0 300 300 2 300 1 300
0 0 0 0 0 0 0 0 0 0
0 0;
22 21.59 0 62.5 -15 1 100 1 50 0
0 0 0 0 0 0 0 0 0 0
0;
];
%% branch data
% fbus tbus r x b rateA rateB rateC ratio
angle status angmin angmax
mpc.branch = [
3 4 0.01 0.04 0 130 130 130 0
0 1 -360 360;
4 6 0.01 0.04 0 90 90 90 0
0 1 -360 360;
6 7 0.03 0.08 0.01 130 130 130 0
0 1 -360 360;
6 8 0.01 0.04 0 32 32 32 0
0 1 -360 360;
6 10 0 0.56 0 32 32 32 0
0 1 -360 360;
9 11 0 0.21 0 65 65 65 0
0 1 -360 360;
9 10 0 0.11 0 65 65 65 0
0 1 -360 360;
4 12 0 0.26 0 65 65 65 0
0 1 -360 360;
12 13 0 0.14 0 65 65 65 0
0 1 -360 360;
12 14 0.12 0.26 0 32 32 32 0
0 1 -360 360;
12 15 0.07 0.13 0 32 32 32 0
0 1 -360 360;
12 16 0.09 0.2 0 32 32 32 0
0 1 -360 360;
15 18 0.11 0.22 0 16 16 16 0
0 1 -360 360;
19 20 0.03 0.07 0 32 32 32 0
0 1 -360 360;
10 20 0.09 0.21 0 32 32 32 0
0 1 -360 360;
10 17 0.03 0.08 0 32 32 32 0
0 1 -360 360;
6 28 0.02 0.06 0.01 32 32 32 0
0 1 -360 360;
10 22 0.07 0.15 0 32 32 32 0 0
1 -360 360;
];
