It looks to me like there may be either some network congestion or some ramping restrictions that are causing the wind output to be curtailed. It’s hard to tell without having a full working script with all input files here to replicate what you are seeing. Can you try eliminating the branch flow limits (set RATE_A to 0 for all branches) and see if that resolves it? And double-check the values of load-following ramp reserves in the wind xGenData inputs.
If you are not able to resolve it with these hints, feel free to submit this as a question in the MOST issue tracker <https://github.com/MATPOWER/most/issues> and attach the files needed to reproduce the problem. Ray > On Oct 19, 2017, at 12:40 PM, Vlad <[email protected]> wrote: > > Dear Ray, > > I'm doing several simulations, based on the t_most_uc.m > > I have the following initial data: > > %% bus data > > % bus_i type Pd Qd Gs Bs area Vm Va > baseKV zone Vmax Vmin > > mpc.bus = [ > > 1 3 0 0 0 0 1 1 0 > 220 1 1.05 0.95; > > 2 2 0 0 0 0 1 1 0 > 110 1 1.05 0.95; > > 3 2 0 0 0 0 1 1 0 > 220 1 1.05 0.95; > > ]; > > > > %% 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 > > mpc.gen = [ > > 1 465 0 50 -50 1 100 1 465 > 186 0 0 0 0 0 0 3 30 90 > 0 0; > > 1 450 0 50 -50 1 100 1 450 > 180 0 0 0 0 0 0 3 30 90 > 0 0; > > 2 72 0 10 -10 1 100 1 72 > 28.8 0 0 0 0 0 0 3 30 90 > 0 0; > > 3 -681 0 0 0 1 100 1 0 > -681 0 0 0 0 0 0 1 10 30 > 0 0; > > ]; > > > > %% branch data > > % fbus tbus r x b rateA rateB rateC ratio > angle status angmin angmax > > mpc.branch = [ > > 1 2 0.005 0.01 0 270 270 270 0 > 0 1 -360 360; > > 1 3 0.005 0.01 0 540 540 540 0 > 0 1 -360 360; > > 2 3 0.005 0.01 0 270 270 270 0 > 0 1 -360 360; > > ]; > > > > %% generator cost data > > % 1 startup shutdown n x1 y1 ... xn > yn > > % 2 startup shutdown n c(n-1) ... c0 > > mpc.gencost = [ > > 2 500 500 2 40 0; > > 2 400 400 2 30 0; > > 2 100 100 2 20 0; > > 2 0 0 2 100 0; > > ]; > > > > Wind data: > > %% 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 > > wind.gen = [ > > 3 0 0 10 -10 1 100 1 100 > 0 0 0 0 0 0 0 0 200 200 > 0 0; > > ]; > > > > }; > > > > wind.xgd_table.data = [ > > 2 1 0 0 1e-8 0 2e-8 0 1e-9 > 1e-9 1e-6 0 1e-6 0; > > ]; > > windprofile.values(:, :, 1) = [ > > 0.378 0.445 0.512; > > 0.427 0.502 0.577; > > 0.463 0.545 0.626; > > 0.493 0.580 0.667; > > 0.561 0.660 0.759; > > 0.537 0.632 0.727; > > 0.561 0.660 0.759; > > 0.583 0.686 0.789; > > 0.610 0.718 0.825; > > 0.641 0.754 0.867; > > 0.660 0.777 0.893; > > 0.689 0.811 0.932; > > 0.724 0.852 0.979; > > 0.739 0.870 1.000; > > 0.728 0.856 0.984; > > 0.695 0.818 0.940; > > 0.662 0.779 0.895; > > 0.606 0.713 0.820; > > 0.554 0.652 0.750; > > 0.511 0.602 0.692; > > 0.343 0.404 0.465; > > 0.188 0.221 0.254; > > 0.251 0.295 0.340; > > 0.452 0.531 0.611; > > ]; > > The problem is that the wind generation should be always commited on and the > power from the wind should go into the network without any regulation. Why > it's not like this? > > The target is to connect my wind power to the network and the other > generators shoud be commited as it should be to follow the load profile but > the power from the wind should go to the network without any bariers. > > ========== PG ========== > > Gen t = 1 t = 2 t = 3 t = 4 t = 5 t = 6 t = 7 t = 8 > t = 9 t =10 t =11 t =12 t =13 t =14 t =15 t =16 t > =17 t =18 t =19 t =20 t =21 t =22 t =23 t =24 > > ---- ------- ------- ------- ------- ------- ------- ------- ------- > ------- ------- ------- ------- ------- ------- ------- ------- > ------- ------- ------- ------- ------- ------- ------- ------- > > 1 186.00 186.00 186.00 186.00 186.00 186.00 186.00 186.00 > 186.00 186.00 186.00 186.00 186.00 186.00 186.00 186.00 > 186.00 186.00 186.00 186.00 186.00 186.00 186.00 186.00 > > 2 400.90 390.90 385.90 400.90 400.90 395.90 366.90 306.90 > 291.90 295.90 299.90 293.90 295.90 291.90 297.90 288.90 > 290.90 281.90 290.90 350.90 382.90 382.90 370.90 360.90 > > 3 72.00 72.00 72.00 72.00 72.00 72.00 72.00 72.00 > 72.00 72.00 72.00 72.00 72.00 72.00 72.00 72.00 > 72.00 72.00 72.00 72.00 72.00 72.00 72.00 72.00 > > 4 -681.00 -671.00 -666.00 -681.00 -681.00 -676.00 -647.00 -587.00 > -572.00 -576.00 -580.00 -574.00 -576.00 -572.00 -578.00 -569.00 > -571.00 -562.00 -571.00 -631.00 -663.00 -663.00 -651.00 -641.00 > > 5 22.10 22.10 22.10 22.10 22.10 22.10 22.10 22.10 > 22.10 22.10 22.10 22.10 22.10 22.10 22.10 22.10 > 22.10 22.10 22.10 22.10 22.10 22.10 22.10 22.10 > > Another problem is that when I change the power of wind from 100 to 500 for > ex. and I have my load profile with the maximum value of -681 MW but in the > output result I have the following: > > ========== PG ========== > > Gen t = 1 t = 2 t = 3 t = 4 t = 5 t = 6 t = 7 t = 8 > t = 9 t =10 t =11 t =12 t =13 t =14 t =15 t =16 t > =17 t =18 t =19 t =20 t =21 t =22 t =23 t =24 > > ---- ------- ------- ------- ------- ------- ------- ------- ------- > ------- ------- ------- ------- ------- ------- ------- ------- > ------- ------- ------- ------- ------- ------- ------- ------- > > 1 - - - - - - - - > - - - - - - - - - > - - - - - - - > > 2 450.00 450.00 450.00 450.00 450.00 450.00 450.00 404.50 > 389.50 393.50 397.50 391.50 393.50 389.50 395.50 386.50 > 388.50 379.50 388.50 448.50 450.00 450.00 450.00 450.00 > > 3 72.00 72.00 72.00 72.00 72.00 72.00 72.00 72.00 > 72.00 72.00 72.00 72.00 72.00 72.00 72.00 72.00 > 72.00 72.00 72.00 72.00 72.00 72.00 72.00 72.00 > > 4 -632.50 -632.50 -632.50 -632.50 -632.50 -632.50 -632.50 -587.00 > -572.00 -576.00 -580.00 -574.00 -576.00 -572.00 -578.00 -569.00 > -571.00 -562.00 -571.00 -631.00 -632.50 -632.50 -632.50 -632.50 > > 5 110.50 110.50 110.50 110.50 110.50 110.50 110.50 110.50 > 110.50 110.50 110.50 110.50 110.50 110.50 110.50 110.50 > 110.50 110.50 110.50 110.50 110.50 110.50 110.50 110.50 > > It seems that the load is also under regulation. In other cases I have my > load with the value of -323. And the wind output also not as expected. Can't > find the reason for that. > > > > With another values of branches and 100 MW wind I have the following: > > > > %% branch data > > % fbus tbus r x b rateA rateB rateC ratio > angle status angmin angmax > > mpc.branch = [ > > 1 2 0.005 0.01 0 30 30 30 0 > 0 1 -360 360; > > 1 3 0.005 0.01 0 540 540 540 0 > 0 1 -360 360; > > 2 3 0.005 0.01 0 60 60 60 0 > 0 1 -360 360; > > ]; > > > > ========== PG ========== > > Gen t = 1 t = 2 t = 3 t = 4 t = 5 t = 6 t = 7 t = 8 > t = 9 t =10 t =11 t =12 t =13 t =14 t =15 t =16 t > =17 t =18 t =19 t =20 t =21 t =22 t =23 t =24 > > ---- ------- ------- ------- ------- ------- ------- ------- ------- > ------- ------- ------- ------- ------- ------- ------- ------- > ------- ------- ------- ------- ------- ------- ------- ------- > > 1 - - - - - - - - > - - - - - - - - - > - - - - - - - > > 2 - - - - - - - - > - - - - - - - - - > - - - - - - - > > 3 72.00 72.00 72.00 72.00 72.00 72.00 72.00 72.00 > 72.00 72.00 72.00 72.00 72.00 72.00 72.00 72.00 > 72.00 72.00 72.00 72.00 72.00 72.00 72.00 72.00 > > 4 -94.10 -94.10 -94.10 -94.10 -94.10 -94.10 -94.10 -94.10 > -94.10 -94.10 -94.10 -94.10 -94.10 -94.10 -94.10 -94.10 > -94.10 -94.10 -94.10 -94.10 -94.10 -94.10 -94.10 -94.10 > > 5 22.10 22.10 22.10 22.10 22.10 22.10 22.10 22.10 > 22.10 22.10 22.10 22.10 22.10 22.10 22.10 22.10 > 22.10 22.10 22.10 22.10 22.10 22.10 22.10 22.10 > > > > Please, help with these cases. > > Vlad . >
