Thank you very much for answering so quickely,
thank's for those advices, i corrected the test case, the last branch
was missing,as well as the two isolated nodes were useless, i should
have seen it i'm sorry. (the case is edited by a script and the one i
work on wasn't the one i sent).
Now that this it is corrected, i have no more "nan" values, but still no
convergence...
What parameters can affects convergence? because when i look at dx some
values are around 1e+2 when the algo stops, maybe to far from the
solution? So an initial guess could help but i have no initial values to
put.
thank you for your help
abdelkrim
Le 18/12/2014 15:51, Ray Zimmerman a écrit :
I found multiple problems with your case file ‘radial_case_no_null.m’.
First of all, buses 30 and 36 are connected to each other, but not to
the rest of the system, and buses 92 and 93 are completely isolated.
Second, check the columns of your gen matrix. It looks to me like you
missed the VG column or something. Currently you have both GEN_STATUS
and PMAX set to 999999. And, by the way, if you mean for it to be
unbounded, simply use Inf, instead of 999999.
Ray
On Dec 18, 2014, at 4:49 AM, ALI ZAZOU Abdelkrim
<[email protected]> wrote:
Hi,
I'm actually stuck with trying to solve the load flow for a radial
distribution, and i need your help because i really don't know what to
do.
The system is composed of 91 nodes and 1 generator, i attached to this
message the test case with null and non null values.
This case don't converge to a result, and matlab tells me:
"Warning: Matrix is singular to working precision.
In newtonpf at 186
In runpf at 226"
So i fisrtly thought that null values in bus Pd, Qd was the problem
but even when i put small value to it, it don't solve the problem.
I tried to see where the values become null or infinite and it comes
from the computation of "dx=-(J/F)"
J and F have non null value but are big (~1e+6)
I also think that BR_R and BR_X are to small, so to be sure here is my
per unit calculation:
S_ref=10Mva, V_ref=20Kv
Z_ref=(20000^2)/(10000000/3)
And for P an Q values are in MW and MVar
And i also tried to run a continuation PF but it can't handle "nan"
values.
So what can you advice to me to try to find a solution ?
Thank you very much
Best regards
Abdelkrim Ali zazou
<radial_case_no_null.m><radial_case_null.m><radial_case.png>
function mpc = case_01_1_dep
%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';
mpc.baseMVA = 10;
%% Bus Data
mpc.bus= [
1 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
2 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0 ;
3 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0 ;
4 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
5 1 0.03375 0.01634587 0 0 1 1
25.84193 20 1 1.05 0 0;
6 1 0.0918 0.04446078 0 0 1 1
25.84193 20 1 1.05 0 0;
7 1 0.0594 0.02876874 0 0 1 1
25.84193 20 1 1.05 0 0;
8 1 0.0828 0.04010188 0 0 1 1
25.84193 20 1 1.05 0 0;
9 1 0.1044 0.05056324 0 0 1 1
25.84193 20 1 1.05 0 0;
10 1 0.0351 0.01699971 0 0 1 1
25.84193 20 1 1.05 0 0;
11 1 0.0189 0.00915369 0 0 1 1
25.84193 20 1 1.05 0 0;
12 1 0.0567 0.02746107 0 0 1 1
25.84193 20 1 1.05 0 0;
13 1 0.0981 0.04751201 0 0 1 1
25.84193 20 1 1.05 0 0;
14 1 0.0243 0.01176903 0 0 1 1
25.84193 20 1 1.05 0 0;
15 1 0.1602 0.07758842 0 0 1 1
25.84193 20 1 1.05 0 0;
16 1 0.0162 0.00784602 0 0 1 1
25.84193 20 1 1.05 0 0;
17 1 0.0666 0.03225586 0 0 1 1
25.84193 20 1 1.05 0 0;
18 1 0.027 0.0130767 0 0 1 1
25.84193 20 1 1.05 0 0;
19 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
20 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
21 1 0.0828 0.04010188 0 0 1 1
25.84193 20 1 1.05 0 0;
22 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
23 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
24 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0 ;
25 1 0.0108 0.00523068 0 0 1 1
25.84193 20 1 1.05 0 0;
26 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
27 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
28 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
29 1 0.0513 0.02484573 0 0 1 1
25.84193 20 1 1.05 0 0;
30 1 0.0216 0.01046136 0 0 1 1
25.84193 20 1 1.05 0 0;
31 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
32 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0 ;
33 1 0.1053 0.05099913 0 0 1 1
25.84193 20 1 1.05 0 0;
34 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
35 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
36 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
37 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
38 1 0.0711 0.03443531 0 0 1 1
25.84193 20 1 1.05 0 0;
39 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
40 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
41 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
42 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
43 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
44 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
45 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
46 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
47 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
48 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
49 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
50 1 0.0162 0.00784602 0 0 1 1
25.84193 20 1 1.05 0 0;
51 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
52 3 0.00001 0.00001 0 0 1 1.03 0 20
1 1.05 0 0 ;
53 1 0.1008 0.04881968 0 0 1 1
25.84193 20 1 1.05 0 0;
54 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0 ;
55 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
56 1 0.0324 0.01569204 0 0 1 1
25.84193 20 1 1.05 0 0;
57 1 0.0999 0.04838379 0 0 1 1
25.84193 20 1 1.05 0 0;
58 1 0.1323 0.06407583 0 0 1 1
25.84193 20 1 1.05 0 0;
59 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
60 1 0.0027 0.00130767 0 0 1 1
25.84193 20 1 1.05 0 0;
61 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0 ;
62 1 0.16875 0.08172938 0 0 1 1
25.84193 20 1 1.05 0 0;
63 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0 ;
64 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
65 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
66 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
67 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
68 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
69 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
70 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
71 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
72 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
73 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
74 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
75 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
76 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
77 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
78 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
79 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
80 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
81 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
82 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
83 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
84 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
85 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
86 1 0.00001 0.00001 0 0 1 1 25.8 20
1 1.05 0 0 ;
87 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
88 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
89 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
90 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
91 1 0.00001 0.00001 0 0 1 1 0 20
1 1.05 0 0;
];
%% generator data
mpc.gen= [
52 3 1.5 Inf -Inf 1 10 1 Inf -Inf;
];
%% Branch data
mpc.branch= [
1 2 0.00079091481 0.0005140946265000001 0 8 8 8 0 0 1 -360 360 ;
4 1 2.1957460133333335e-5 9.033132333333336e-6 0 6 6 6 0 0 1 -360 360 ;
2 3 0.0005720009866666666 0.0003718006413333333 0 8 8 8 0 0 1 -360 360 ;
3 63 1.085185e-6 7.0537025e-7 0 8 8 8 0 0 1 -360 360 ;
65 19 0.0004964366608499999 0.00018061624875 0 2 2 2 0 0 1 -360 360 ;
18 53 0.0004413146533333334 0.0002868545246666667 0 8 8 8 0 0 1 -360 360 ;
51 18 5.060353333333334e-6 3.289229666666667e-6 0 8 8 8 0 0 1 -360 360 ;
26 68 0.0034558943249500002 0.0027615593920833335 0 4 4 4 0 0 1 -360 360 ;
35 23 0.0007813979533333334 0.00032146118333333337 0 6 6 6 0 0 1 -360 360 ;
69 70 0.0057047391906 0.004558581545 0 4 4 4 0 0 1 -360 360 ;
71 26 0.0035307791776500006 0.0028213988862499997 0 4 4 4 0 0 1 -360 360 ;
35 72 0.0018416064926000003 0.0014716033616666666 0 4 4 4 0 0 1 -360 360 ;
73 35 0.0023422466793499995 0.0018716583054166664 0 4 4 4 0 0 1 -360 360 ;
74 6 0.0017254742886333334 0.0006277713108333334 0 2 2 2 0 0 1 -360 360 ;
56 34 0.0016937272865500002 0.0006162209462499999 0 2 2 2 0 0 1 -360 360 ;
12 75 0.0031468124616499996 0.0011448901887499998 0 2 2 2 0 0 1 -360 360 ;
36 30 4.3935405000000013e-5 2.855801325e-5 0 8 8 8 0 0 1 -360 360 ;
53 31 0.0005246659250000001 0.00034103285125 0 8 8 8 0 0 1 -360 360 ;
48 27 5.1925578200000006e-5 2.13617885e-5 0 6 6 6 0 0 1 -360 360 ;
72 33 0.0044105512957000005 0.0035244131358333334 0 4 4 4 0 0 1 -360 360 ;
73 50 0.00036866939974999994 0.0002945988354166666 0 4 4 4 0 0 1 -360 360 ;
76 47 0.0010931357150500002 0.0008735102745833333 0 4 4 4 0 0 1 -360 360 ;
76 73 0.0009360721964 0.0007480028966666667 0 4 4 4 0 0 1 -360 360 ;
19 5 0.0008050901360666667 0.00029291221166666666 0 2 2 2 0 0 1 -360 360 ;
37 38 0.0000001 0.0000001 0 2 2 2 0 0 1 -360 360 ;
38 77 0.0000001 0.0000001 0 2 2 2 0 0 1 -360 360 ;
39 37 0.0000001 0.0000001 0 2 2 2 0 0 1 -360 360 ;
77 33 0.00032999706500000007 0.00021449809225000004 0 8 8 8 0 0 1 -360 360 ;
34 32 0.00022064490123333333 9.077163658333333e-5 0 6 6 6 0 0 1 -360 360 ;
34 40 0.0008979513868 0.00032669748999999997 0 2 2 2 0 0 1 -360 360 ;
40 72 0.004124282800516667 0.001500518690416667 0 2 2 2 0 0 1 -360 360 ;
78 7 0.0011785434983833333 0.0004287840170833333 0 2 2 2 0 0 1 -360 360 ;
78 41 0.0009989464838166668 0.00036344206791666666 0 2 2 2 0 0 1 -360 360 ;
64 9 0.00060720191805 0.00022091545874999998 0 2 2 2 0 0 1 -360 360 ;
64 42 0.00018438561006666666 6.708416166666666e-5 0 2 2 2 0 0 1 -360 360 ;
80 52 3.309269183333333e-5 1.3614082083333334e-5 0 6 6 6 0 0 1 -360 360 ;
81 13 0.0004699604888333333 0.0001709835458333333 0 2 2 2 0 0 1 -360 360 ;
81 44 0.0034228753730166663 0.0012453288779166664 0 2 2 2 0 0 1 -360 360 ;
82 21 0.0018645889964000002 0.00067838477 0 2 2 2 0 0 1 -360 360 ;
82 81 0.004525359672166667 0.0016464406291666663 0 2 2 2 0 0 1 -360 360 ;
14 82 0.0017651154466833334 0.0006421937695833334 0 2 2 2 0 0 1 -360 360 ;
83 46 0.00026861798535 9.773003625e-5 0 2 2 2 0 0 1 -360 360 ;
16 83 0.0033262158496666664 0.0012101616916666666 0 2 2 2 0 0 1 -360 360 ;
66 55 0.0026072950037333334 0.0009486000533333333 0 2 2 2 0 0 1 -360 360 ;
66 47 0.0012992917181833335 0.0004727152820833334 0 2 2 2 0 0 1 -360 360 ;
84 85 0.0044402999367999995 0.00161549374 0 2 2 2 0 0 1 -360 360 ;
50 84 0.0010149622526333332 0.0003692690108333333 0 2 2 2 0 0 1 -360 360 ;
86 28 3.911652333333333e-6 1.4231583333333332e-6 0 2 2 2 0 0 1 -360 360 ;
70 46 0.0004901945154166666 0.00017834519791666664 0 2 2 2 0 0 1 -360 360 ;
69 44 0.001162468045483333 0.00042293535958333326 0 2 2 2 0 0 1 -360 360 ;
69 45 0.0010746510497666667 0.00039098530916666666 0 2 2 2 0 0 1 -360 360 ;
68 42 0.0010267242739 0.00037354833249999994 0 2 2 2 0 0 1 -360 360 ;
87 71 0.0007138366438666667 0.00025971187666666665 0 2 2 2 0 0 1 -360 360 ;
86 87 0.0028728646676166665 0.0010452210329166664 0 2 2 2 0 0 1 -360 360 ;
29 86 0.0023903865145333335 0.0008696832433333333 0 2 2 2 0 0 1 -360 360 ;
88 10 1.6498781e-6 6.002674999999998e-7 0 2 2 2 0 0 1 -360 360 ;
23 24 0.0007645104867333334 0.00031451380783333336 0 6 6 6 0 0 1 -360 360 ;
67 88 0.0036057164840666667 0.0013118511116666667 0 2 2 2 0 0 1 -360 360 ;
58 87 2.6781052666666663e-6 1.1017521666666667e-6 0 6 6 6 0 0 1 -360 360 ;
59 58 2.3557649899999998e-5 9.69143825e-6 0 6 6 6 0 0 1 -360 360 ;
89 41 0.0006974911335166666 0.00025376496541666666 0 2 2 2 0 0 1 -360 360 ;
89 67 0.0029338324597833333 0.0010674026620833333 0 2 2 2 0 0 1 -360 360 ;
43 89 0.0010194760047333331 0.0003709112283333333 0 2 2 2 0 0 1 -360 360 ;
60 61 8.525732666666668e-6 3.507421666666667e-6 0 6 6 6 0 0 1 -360 360 ;
60 62 0.00017986328333333335 0.00011691113416666666 0 8 8 8 0 0 1 -360 360 ;
54 62 0.0008368968083333333 0.0005439829254166667 0 8 8 8 0 0 1 -360 360 ;
54 90 3.958368333333333e-6 2.5729394166666664e-6 0 8 8 8 0 0 1 -360 360 ;
17 90 0.004442779884616666 0.0016163960079166664 0 2 2 2 0 0 1 -360 360 ;
90 83 0.0005929933864833333 0.0002157460345833333 0 2 2 2 0 0 1 -360 360 ;
74 43 0.0031021882478166666 0.0011286547679166667 0 2 2 2 0 0 1 -360 360 ;
65 74 0.0017193581812333333 0.0006255461158333333 0 2 2 2 0 0 1 -360 360 ;
91 57 0.0011836329674000002 0.00043063569499999995 0 2 2 2 0 0 1 -360 360 ;
91 65 0.0005486792172333333 0.00019962341583333333 0 2 2 2 0 0 1 -360 360 ;
11 91 0.0016498279156666666 0.0006002492416666667 0 2 2 2 0 0 1 -360 360 ;
75 21 0.0007215457030166667 0.00026251662791666665 0 2 2 2 0 0 1 -360 360 ;
75 48 0.0029021376381833334 0.0010558712820833331 0 2 2 2 0 0 1 -360 360 ;
15 63 0.0009936733289166667 0.0003615235604166667 0 2 2 2 0 0 1 -360 360 ;
63 45 0.0010615525539666667 0.00038621974416666664 0 2 2 2 0 0 1 -360 360 ;
19 20 0.0007126708825 0.0004632360925 0 8.00000000 8.00000000 8.00000000 0 0 1
-360 360 ;
64 88 0.0086973985 0.0031643341 0 2.00000000 2.00000000 2.00000000 0 0 1 -360
360 ;
66 49 0.004808182516666666 0.0017493387083333332 0 2.00000000 2.00000000
2.00000000 0 0 1 -360 360 ;
67 8 0.003954792516666667 0.0014388538916666666 0 2.00000000 2.00000000
2.00000000 0 0 1 -360 360 ;
68 25 0.00583548645 0.004663059608333334 0 4.00000000 4.00000000 4.00000000 0 0
1 -360 360 ;
69 25 0.0017030360833333335 0.0013608735083333332 0 4.00000000 4.00000000
4.00000000 0 0 1 -360 360 ;
70 22 0.000738570715 0.0005901820958333334 0 4.00000000 4.00000000 4.00000000 0
0 1 -360 360 ;
71 80 0.005592975516666667 0.004469272991666666 0 4.00000000 4.00000000
4.00000000 0 0 1 -360 360 ;
76 22 0.0008360097833333333 0.0006680443883333334 0 4.00000000 4.00000000
4.00000000 0 0 1 -360 360 ;
78 79 0.001037387175 0.00037742778666666665 0 2.00000000 2.00000000 2.00000000
0 0 1 -360 360 ;
80 31 0.0011012909333333334 0.0008019279691666666 0 4.00000000 4.00000000
4.00000000 0 0 1 -360 360 ;
36 84 5.9327178916666666667e-4 2.15847309166666667e-4 0 4 4 4 0 0 1 -360 360;
];
