GAMS系列分享8--电力系统最优潮流—24母线最优潮流

      本节重点：   24母线电力系统

                                      基于直流潮流

                                                       最优潮流计算

 

   编程参考：GAMS系列分享6--电力系统最优潮流—3母线最优潮流.

                    物理模型，数学模型都一模一样，编程只是数据量大一点，一模一样

   重点：一模一样

   

*DC-opf five bus

*1，集合--------------------------------------------------------------------------------------------
Sets bus /1*24/,slack(bus) /13/,gen /g1*g12/;
scalars sbase /100/;
alias (bus,node);

*2，数据输入----------------------------------------------------------------------------------------
table gendata(gen,*)
      pmax    pmin    b     costst  costsd  RU    RD    SU   SD    UT     DT   uini   U0   S0
g1    400     100     5.47  0        0      47    47    105  108   1      1    1      5    0
g2    400     100     5.47  0        0      47    47    106  112   1      1    1      6    0
g3    152     30.4    13.32 1430.4   1430.4 14    14    43   45    8      4    1      2    0
g4    152     30.4    13.32 1430.4   1430.4 14    14    44   57    8      8    1      2    0
g5    155     54.25   16    0        0      21    21    65   77    8      8    0      0    2
g6    155     54.25   10.52 312      312    21    21    66   73    8      8    1      10   0
g7    310     108.5   10.52 624      624    21    21    112  125   8      8    1      10   0
g8    350     140     10.89 2298     2298   28    28    154  162   8      8    1      5    0
g9    350     75      20.7  1725     1725   49    49    77   80    8      8    0      0    2
*g10   591     206.85  20.93 3056.7   3056.7 21    21    213  228   12     10   0      0    8
g11   60      12      26.11 437      437    7     7     19   31    4      2    0      0    1
g12   300     0       0     0        0      35    35    315  326   0      0    1      2    0;

set GBconect(bus,gen)
/18 . g1
 21 . g2
 1  . g3
 2  . g4
15 . g5
 16 . g6
 23 . g7
 23 . g8
 7  . g9
*13 . g10
 15 . g11
 22 . g12 /;
table busData(bus,*)
     pd      qd
1    108     22
2    97      20
3    180     37
4    74      15
5    71      14
6    136     28
7    125     25
8    171     35
9    175     36
10   195     40
13   265     54
14   194     39
15   317     64
16   100     20
18   333     68
19   181     37
20   128     26 ;
set conex
/
1 . 2
1 . 3
1 . 5
2 . 4
2 . 6
3 . 9
3 . 24
4 . 9
5 . 10
6 . 10
7 . 8
8 . 9
8 . 10
9 . 11
9 . 12
10. 11
10. 12
11. 13
11. 14
12. 13
12. 23
13. 23
14. 16
15. 16
15. 21
15. 24
16. 17
16. 19
17. 18
17. 22
18. 21
19. 20
20. 23
21. 22
/;
conex(bus,node)$conex(node,bus)=1;
table branch(bus,node,*)
        r          x          b        limit
1 . 2   0.0026     0.0139     0.4611   175
1 . 3   0.0546     0.2112     0.0572   175
1 . 5   0.0218     0.0845     0.0229   175
2 . 4   0.0328     0.1267     0.0343   175
2 . 6   0.0497     0.1920     0.0520   175
3 . 9   0.0308     0.1190     0.0322   175
3 . 24  0.0023     0.0839     0.0000   400
4 . 9   0.0268     0.1037     0.0281   175
5 . 10  0.0228     0.0883     0.0239   175
6 . 10  0.0139     0.0605     2.4590   175
7 . 8   0.0159     0.0614     0.0166   175
8 . 9   0.0427     0.1651     0.0447   175
8 . 10  0.0427     0.1651     0.0447   175
9 . 11  0.0023     0.0839     0.0000   400
9 . 12  0.0023     0.0839     0.0000   400
10. 11  0.0023     0.0839     0.0000   400
10. 12  0.0023     0.0839     0.0000   400
11. 13  0.0061     0.0476     0.0999   500
11. 14  0.0054     0.0418     0.0879   500
12. 13  0.0061     0.0476     0.0999   500
12. 23  0.0124     0.0966     0.2030   500
13. 23  0.0111     0.0865     0.1818   500
14. 16  0.0050     0.0389     0.0818   500
15. 16  0.0022     0.0173     0.0364   500
15. 21  0.0032     0.0245     0.2060   1000
15. 24  0.0067     0.0519     0.1091   500
16. 17  0.0033     0.0259     0.0545   500
16. 19  0.0030     0.0231     0.0485   500
17. 18  0.0018     0.0144     0.0303   500
17. 22  0.0135     0.1053     0.2212   500
18. 21  0.0017     0.0130     0.1090   1000
19. 20  0.0026     0.0198     0.1666   1000
20. 23  0.0014     0.0108     0.0910   1000
21. 22  0.0087     0.0678     0.1424   500  ;
* 研究线路传输容量对最优潮流的影响
*branch(bus,node,'limit') = 0.7*branch(bus,node,'limit');
* 研究线路开断对最优潮流的影响
*branch('19','20','limit') = 0;
*branch('12','23','limit') = 0;


branch(bus,node,'x')$(branch(bus,node,'x')=0)=branch(node,bus,'x');
branch(bus,node,'limit')$(branch(bus,node,'limit')=0)=branch(node,bus,'limit');
branch(bus,node,'bij')$conex(bus,node) = 1/branch(bus,node,'x');

*3,定义变量-----------------------------------------------------------------------------------------
Variables Pij(bus,node),pg(gen),delta(bus),of;

*4,方程声明和定义-----------------------------------------------------------------------------------
Equations const1,const2,const3;
const1(bus,node)$conex(bus,node)..pij(bus,node)=e=branch(bus,node,'bij')*(delta(bus)-delta(node));
const2(bus)..+sum(gen$gbconect(bus,gen),pg(gen))-busdata(bus,'pd')/sbase
                  =e=sum(node$conex(node,bus),pij(bus,node));
const3..OF=g=sum(gen,pg(gen)*gendata(gen,'b')*sbase);
model loadflow /all/;

*5,变量约束，变量初值——---------------------------------------------------------------------------
pg.lo(gen)=gendata(gen,'Pmin')/sbase;
pg.up(gen)=gendata(gen,'pmax')/sbase;
delta.up(bus)=pi/2;
delta.lo(bus)=-pi/2;
delta.fx(slack) =0;
pij.up(bus,node)$conex(bus,node)=branch(bus,node,'limit')/sbase;
pij.lo(bus,node)$conex(bus,node)=-branch(bus,node,'limit')/sbase;
*模型求解和，结果展示-------------------------------------------------------------------------------
solve loadflow min of us lp;


parameter report(bus,*),congestioncost;
report(bus,'gen(MW)') = sum(gen$gbconect(bus,gen),pg.l(gen))*sbase;
report(bus,'angle')   = delta.l(bus);
report(bus,'LMP($/MW)') = const2.m(bus)/sbase;
congestioncost = sum((bus,node),pij.l(bus,node)*(-const2.m(bus)+const2.m(node)))/2;
display report,pij.l,congestioncost;

 

欢迎观看系列文章，目前更新到电力系统的最优潮流

欢迎观看系列文章，目前更新到电力系统的最优潮流

欢迎观看系列文章，目前更新到电力系统的最优潮流（全部手码）

 

程序可以执行，已经验证过！！！！！

 

最好自己编程实现一下。。