Please see the base case file in attachment
Natakorn
2015-02-13 0:57 GMT+07:00 Abhyankar, Shrirang G. <[email protected]>:
> Can you please send your base case file.
>
> Shri
>
> From: Natakorn Thasnas <[email protected]>
> Reply-To: MATPOWER discussion forum <[email protected]>
> Date: Fri, 13 Feb 2015 00:21:15 +0700
> To: MATPOWER discussion forum <[email protected]>
> Subject: Re: PV Curve
>
> Hi Shri
>
> I tried using a smaller continuation step size as your suggestion (step
> size=0.05). And I tried decreasing step size to 0.002,
> I still got the same result, the PV curve is not full nose curve.
>
> Best Regards,
>
> Natakorn Thasnas
> Ph.D Student
> Khonkaen University, Thailand
>
> 2015-02-12 20:41 GMT+07:00 Abhyankar, Shrirang G. <[email protected]>:
>
>> Try using a smaller continuation step size.
>> mpopt = mpoption(mpopt, 'cpf.stop_at', 'full', 'cpf.step', 0.05);
>>
>> Shri
>>
>> > On Feb 12, 2015, at 4:03 AM, "Natakorn Thasnas" <
>> [email protected]> wrote:
>> >
>> > Hi All
>> >
>> > I am trying to run a continuation power flow with the following code
>> >
>> > define_constants;
>> > mpopt = mpoption('out.all', 0, 'verbose', 2);
>> > mpopt = mpoption(mpopt, 'cpf.stop_at', 'full', 'cpf.step', 0.2);
>> > mpopt = mpoption(mpopt, 'cpf.plot.level', 2 ,'cpf.plot.bus', 14);
>> > mpcb = loadcase('case14result'); % load base case
>> > mpct = mpcb; % set up target case with
>> > mpct.gen(:, [PG QG]) = mpcb.gen(:, [PG QG])*2.5;
>> > mpct.bus(:,[PD QG]) = mpcb.bus(:,[PD QG])*2.5;
>> > results = runcpf(mpcb, mpct, mpopt);
>> >
>> > After running it, I got the warning message "corrector did not converge
>> in 1000 iterations" and the PV curve is not full nose curve as showing in
>> attached file (c14.jpg)
>> >
>> > <c14.jpg>
>> >
>> > Then I increased "max_it" ​from 1000 to 10000. I got the warning
>> messages as following
>> >
>> > Warning: Matrix is close to singular or badly scaled.
>> > Results may be inaccurate. RCOND = 1.944067e-018.
>> > > In cpf_corrector at 143
>> > In runcpf at 268
>> > In cpfrun at 55
>> > step 48 : lambda = -7.770, corrector did not converge in 10000
>> iterations
>> >
>> >
>> > Could you please suggest me how to fix it to get "the full nose curve"?
>> >
>> > Thank you in advance.
>> >
>> > Best Regards,
>> >
>> > Natakorn Thasnas
>> > Ph.D Student
>> > Khonkaen University, Thailand
>>
>>
>
function mpc = case14test
%CASE14TEST
%% 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 lam_P lam_Q mu_Vmax mu_Vmin
mpc.bus = [
1 3 0 0 0 0 1 1.05999209
0 0 1 1.06 0.94 36.7229 -0.0945 584.4041 0.0136;
2 2 21.7 12.7 0 0 1 1.04072114
-4.02176403 0 1 1.06 0.94 38.3586 0.0003 0.0849 0.0298;
3 2 94.2 19 0 0 1 1.01557791
-9.92501388 0 1 1.06 0.94 40.5740 0.0008 0.0368 0.0662;
4 1 47.8 -3.9 0 0 1 1.01441589
-8.66430744 0 1 1.06 0.94 40.1894 0.1204 0.0359 0.0444;
5 1 7.6 1.6 0 0 1 1.01631851
-7.42772074 0 1 1.06 0.94 39.6599 0.2081 0.0375 0.0398;
6 2 11.2 7.5 0 0 1 1.05989599
-12.6881908 0 1 1.06 0.94 39.7327 0.0007 55.5000 0.0136;
7 1 0 0 0 0 1 1.04629698
-11.1878723 0 1 1.06 0.94 40.1709 0.1202 0.1194 0.0207;
8 2 0 0 0 0 1 1.05997098
-10.4151518 0 1 1.06 0.94 40.1692 0.0003 71.2138 0.0136;
9 1 29.5 16.6 0 19 1 1.0436341
-12.9972796 0 1 1.06 0.94 40.1656 0.1965 0.1000 0.0216;
10 1 9 5.8 0 0 1 1.03906422
-13.2328806 0 1 1.06 0.94 40.3172 0.3094 0.0781 0.0219;
11 1 3.5 1.8 0 0 1 1.04592126
-13.0904801 0 1 1.06 0.94 40.1547 0.2287 0.1162 0.0179;
12 1 6.1 1.6 0 0 1 1.04471806
-13.5318927 0 1 1.06 0.94 40.3782 0.2129 0.0910 0.0160;
13 1 13.5 5.8 0 0 1 1.03984873
-13.5818675 0 1 1.06 0.94 40.5747 0.3541 0.0576 0.0173;
14 1 14.9 5 0 0 1 1.02380687
-14.2740305 0 1 1.06 0.94 41.1970 0.5715 0.0144 0.0253;
];
%% 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 mu_Pmax mu_Pmin mu_Qmax mu_Qmin
mpc.gen = [
1 194.316746 0.0618531273 10 0 1.05999209
100 1 332.4 0 0 0 0 0 0 0
0 0 0 0 0 0.0001 0.0000 0.0002 0.0946;
2 36.7182776 23.6635636 50 -40 1.04072114
100 1 140 0 0 0 0 0 0 0
0 0 0 0 0 0.0000 0.0004 0.0002 0.0000;
3 28.7515402 24.1151879 40 0 1.01557791
100 1 100 0 0 0 0 0 0 0
0 0 0 0 0 0.0000 0.0009 0.0006 0.0001;
6 0.00901839562 11.5078717 24 -6 1.05989599
100 1 100 0 0 0 0 0 0 0
0 0 0 0 0 0.0000 0.2675 0.0005 0.0001;
8 8.49088543 8.28549621 24 -6 1.05997098
100 1 100 0 0 0 0 0 0 0
0 0 0 0 0 0.0000 0.0006 0.0002 0.0001;
];
%% branch data
% fbus tbus r x b rateA rateB rateC ratio
angle status angmin angmax Pf Qf Pt Qt mu_Sf mu_St
mu_angmin mu_angmax
mpc.branch = [
1 2 0.01938 0.05917 0.0528 0 0 0 0
0 1 -360 360 129.6607 -6.3204 -126.7590 9.3541
0.0000 0.0000 0.0000 0.0000;
1 5 0.05403 0.22304 0.0492 0 0 0 0
0 1 -360 360 64.6560 6.3823 -62.6055 -3.2227 0.0000
0.0000 0.0000 0.0000;
2 3 0.04699 0.19797 0.0438 0 0 0 0
0 1 -360 360 55.5872 0.4828 -54.2431 0.5492 0.0000
0.0000 0.0000 0.0000;
2 4 0.05811 0.17632 0.034 0 0 0 0
0 1 -360 360 48.9142 -0.4710 -47.6296 0.7783 0.0000
0.0000 0.0000 0.0000;
2 5 0.05695 0.17388 0.0346 0 0 0 0
0 1 -360 360 37.2759 1.5976 -36.5389 -3.0083 0.0000
0.0000 0.0000 0.0000;
3 4 0.06701 0.17103 0.0128 0 0 0 0
0 1 -360 360 -11.2054 4.5660 11.3047 -5.6312 0.0000
0.0000 0.0000 0.0000;
4 5 0.01335 0.04211 0 0 0 0 0
0 1 -360 360 -49.1666 11.5740 49.4976 -10.5299
0.0000 0.0000 0.0000 0.0000;
4 7 0 0.20912 0 0 0 0 0.978
0 1 -360 360 22.8501 -3.9914 -22.8501 5.0373 0.0000
0.0000 0.0000 0.0000;
4 9 0 0.55618 0 0 0 0 0.969
0 1 -360 360 14.8414 1.1703 -14.8414 -0.0455 0.0000
0.0000 0.0000 0.0000;
5 6 0 0.25202 0 0 0 0 0.932
0 1 -360 360 42.0469 15.1609 -42.0469 -10.9269
0.0000 0.0000 0.0000 0.0000;
6 11 0.09498 0.1989 0 0 0 0 0
0 1 -360 360 6.0877 4.5535 -6.0388 -4.4512 0.0000 0.0000
0.0000 0.0000;
6 12 0.12291 0.25581 0 0 0 0 0
0 1 -360 360 7.6515 2.6593 -7.5797 -2.5098 0.0000 0.0000
0.0000 0.0000;
6 13 0.06615 0.13027 0 0 0 0 0
0 1 -360 360 17.1167 7.7219 -16.9091 -7.3130 0.0000
0.0000 0.0000 0.0000;
7 8 0 0.17615 0 0 0 0 0
0 1 -360 360 -8.4909 -8.0648 8.4909 8.2855 0.0000 0.0000
0.0000 0.0000;
7 9 0 0.11001 0 0 0 0 0
0 1 -360 360 31.3410 3.0276 -31.3410 -2.0313 0.0000
0.0000 0.0000 0.0000;
9 10 0.03181 0.0845 0 0 0 0 0
0 1 -360 360 6.4866 3.2131 -6.4713 -3.1725 0.0000 0.0000
0.0000 0.0000;
9 14 0.12711 0.27038 0 0 0 0 0
0 1 -360 360 10.1958 2.9580 -10.0643 -2.6782 0.0000
0.0000 0.0000 0.0000;
10 11 0.08205 0.19207 0 0 0 0 0
0 1 -360 360 -2.5287 -2.6275 2.5388 2.6512 0.0000 0.0000
0.0000 0.0000;
12 13 0.22092 0.19988 0 0 0 0 0
0 1 -360 360 1.4797 0.9098 -1.4736 -0.9043 0.0000 0.0000
0.0000 0.0000;
13 14 0.17093 0.34802 0 0 0 0 0
0 1 -360 360 4.8827 2.4173 -4.8357 -2.3218 0.0000 0.0000
0.0000 0.0000;
];
%%----- OPF Data -----%%
%% generator cost data
% 1 startup shutdown n x1 y1 ... xn
yn
% 2 startup shutdown n c(n-1) ... c0
mpc.gencost = [
2 0 0 3 0.0430292599 20 0;
2 0 0 3 0.25 20 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
2 0 0 3 0.01 40 0;
];
%% bus names
mpc.bus_name = {
'Bus 1 HV';
'Bus 2 HV';
'Bus 3 HV';
'Bus 4 HV';
'Bus 5 HV';
'Bus 6 LV';
'Bus 7 ZV';
'Bus 8 TV';
'Bus 9 LV';
'Bus 10 LV';
'Bus 11 LV';
'Bus 12 LV';
'Bus 13 LV';
'Bus 14 LV';
};
