1) My question was whether you checked the MOST results for the cases that run
but do not include any startup or shutdown of units to confirm that the total
generation is tracking the load profile you are attempting to specify. If there
is an error somewhere in the inputs that prevents the load profile being
applied as expected, it would result in a fixed commitment schedule like this.
2) That means “infeasible or unbounded”.
Ray
> On Sep 10, 2018, at 2:13 PM, Michal Polecki <[email protected]>
> wrote:
>
> Thank you for your answer.
>
> 1) I have already checked if the load can be covered by generation and if the
> ramps of generators are suffcient. I'll do it one more time, thank you.
>
> 2) Gurobi gives exitflag -4.
>
> Best Regards,
> Michał Połecki
>
> On 09/10/18 16:20, Ray Zimmerman <[email protected]> wrote:
>>
>> Did you check that the load profile is working correctly? Is the total
>> generation dispatch following the load profile you expected?
>>
>> Also, what do you mean “it is impossible for MATPOWER to calculate” for the
>> 24 hour horizons? Does Gurobi run forever, find that the problem is
>> infeasible, or what?
>>
>> Ray
>>
>>
>>
>> > On Sep 10, 2018, at 8:24 AM, Michal Polecki <[email protected]>
>> > wrote:
>> >
>> > Dear Creators and Users of Matpower,
>> >
>> > I have a problem with Unit Commiment problem. I've studied all examples
>> > attached to Matpower and also I've created some more cases before going to
>> > my large problem. It seemed to give me correct results. I'm using Gurobi.
>> > Now I'm trying to run 118 bus case with 13 additional wind farms. What I
>> > have observed is, that MATPOWER is sticking to the first plan of work (no
>> > startups or shutdowns though cost of it is 0). What is more is some cases
>> > it is impossible for MATPOWER to calculate long period e.i nt = 24. I have
>> > written simple script to cut this long period into 4 6-hour periods. Then
>> > Matpower was able to calculate each of them. Plot_UC of each 6-hour period
>> > is attached.
>> >
>> > Do you have an idea what could be problem?
>> >
>> >
>> > nt = 6;
>> > mpc = loadcase('MMcase118v5');
>> > Pg = mpc.gen(:,PG);
>> >
>> > for i = 1:4
>> > result = MM_118bus_UC(nt,Xwind,Xload,cw,cl,Pg);
>> > if result.results.f>0
>> > res(j) = result;
>> > Pg = result.results.Pc([1:19],nt);
>> > cw = cw+nt;
>> > cl = cl+nt;
>> > j= j+1;
>> > else
>> > break
>> > end
>> > end
>> >
>> > Michał Połecki
>> >
>> > --
>> > Michał Połecki <P1.jpg><P2.jpg><P3.jpg><P4.jpg>
>>
>>
>>
> --
> Michał Połecki
