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 .
