I’m not sure why, but it appears that the OPF is probably not converging
successfully. Did you check results.success to see that it converged? I am not
able to reproduce your result. When I execute the 3 lines of code you have
below, I get:
================================================================================
| Reserves |
================================================================================
Gen Bus Status Reserves Price
# # (MW) ($/MW) Included in Zones ...
---- ----- ------ -------- -------- ------------------------
1 1 1 25.00 2.00 1
2 2 1 15.00 2.00 1
3 22 1 0.00 2.00 1
4 27 1 0.00 2.00 1
5 23 1 19.39 5.50 1, 2
6 13 1 0.61 5.50 1, 2
--------
Total: 60.00 Total Cost: $177.81
Zone Reserves Price
# (MW) ($/MW)
---- -------- --------
1 60.00 2.00
2 20.00 3.50
================================================================================
| Reserve Limits |
================================================================================
Gen Bus Status Rmin mu Rmin Reserves Rmax Rmax mu Pmax mu
# # ($/MW) (MW) (MW) (MW) ($/MW) ($/MW)
---- ----- ------ -------- -------- -------- -------- -------- --------
1 1 1 - 0.00 25.00 25.00 0.10 0.00
2 2 1 - 0.00 15.00 25.00 - 0.00
3 22 1 1.00 0.00 0.00 25.00 - -
4 27 1 2.00 0.00 0.00 25.00 - 0.00
5 23 1 - 0.00 19.39 25.00 0.00 0.50
6 13 1 0.00 0.00 0.61 25.00 - -
--------
Total: 60.00
… which looks like a correct result.
--
Ray Zimmerman
Senior Research Associate
B30 Warren Hall, Cornell University, Ithaca, NY 14853 USA
phone: (607) 255-9645
On Jul 29, 2014, at 10:05 AM, Sohail Khan <[email protected]> wrote:
> Dear all,
>
> I have two questions:
>
> Question 1
> While running the OPF extension of the fixed zonal reserves extension of OPF
> in the default settings that is simulating following code
>
> mpc = loadcase('t_case30_userfcns');
> mpc = toggle_reserves(mpc, 'on');
> results = runopf(mpc);
>
> leading to the result:
> Zone Reserves Price
> # (MW) ($/MW)
> ---- -------- --------
> 1 49.7353975328242405556000.00
> 2 16.58 36.71
>
> But the zonal requirement is:
> %% reserve requirements for each zone in MW
> mpc.reserves.req = [60; 20];
>
> Then according to equation 6.5 of the manual the additional constraint 'Rreq'
> should be satisfied but it is not as shown above. Even while decrease the
> zone requirements to
> mpc.reserves.req = [10; 10]; the result does not satisfy the constraint
>
> Zone Reserves Price
> # (MW) ($/MW)
> ---- -------- --------
> 1 19.12 13.25
> 2 7.842960557094756162200000.00
>
> Why is it so?
>
> Question 2
> My second question is that while the reserve cost is different for different
> generators
> mpc.reserves.cost = [ 1.9; 2; 3; 4; 5; 5.5 ];
>
> the reserve distribution among generators doesn't reflect the difference in
> the amount of reserve power assigned such that reserves allocated per
> generator 8.29 MW.
>
> Why is it the same?
>
> Please guide me.
>
> Thank you!
> Sohail Khan
