I think you must many thing change . kind if bus I means kind 1 = slack , 2 =PV and 3= PQ if you wnat put generator on a bus , you must change change type too, to add generator , and add cost function
On Thu, Sep 13, 2018 at 7:40 PM Michal Polecki <[email protected]> wrote: > Dear All, > > I've wrote a script which is checking longest possible number of terms in > most. In my case it was 152 terms (for 118 bus case with 13 wind farms). As > I previously written it is not clear to me why the Gurobi is chosing the > solution in which there is no Startup and shutdown. In case Pmin of > convetional generation greater than Pload it is curtailing wind generation > (which gencost is zero). I've attached plut_uc of the result file. > > This is options I have. > > mpopt = mpoption('verbose', verbose); > mpopt = mpoption(mpopt, 'out.gen', 1); > mpopt = mpoption(mpopt, 'model', 'AC'); > mpopt = mpoption(mpopt, 'most.skip_prices',1); > mpopt = mpoption(mpopt, 'most.solver', 'gurobi'); > mpopt = mpoption(mpopt, 'most.uc.run', '1'); > mpopt = mpoption(mpopt, 'gurobi.method', -1); %% automatic > % mpopt = mpoption(mpopt, 'gurobi.method', 0); %% primal simplex > % mpopt = mpoption(mpopt, 'gurobi.method', 1); %% dual simplex > % mpopt = mpoption(mpopt, 'gurobi.method', 2); %% barrier > mpopt = mpoption(mpopt, 'gurobi.threads', 2); > mpopt = mpoption(mpopt, 'gurobi.opts.MIPGap', 0); > mpopt = mpoption(mpopt, 'gurobi.opts.MIPGapAbs', 0); > > Best Regards, > Michał Połecki > > On 09/10/18 21:15, *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 > > -- > Michał Połecki
