GAMS系列分享15—电力系统优化运行—动态经济调度

   四台火电机组，给电负荷供能。

                                 24小时

                                            求最优的机组出力分配

                                                   不考虑网络约束

本期看着有点水（哈哈哈哈），主要是为后期目标函数的线性化做准备的。

目录

1，物理模型

2，数学模型

3，代码

4，难点

---

---

1，物理模型

           四台火电机组给电负荷供能，求最优的各机组的出力大小。

                                                                 

2，数学模型

                      

             模型比较简单，不做过多介绍，下节将对目标函数进行线性化处理。

3，代码

$Title Aircraft allocation under uncertain demand  (AIRCRAF,SEQ=8)

$Ontext

   动态经济分配
   仅考虑火电机组的运行成本

$Offtext


*集合-----------------------------------------------------------------------------------------------
Sets  t  hours          /t1*t24/
      i  terminal units /g1*g4/;
table gendata(i,*)
       a     b      c      d      e       f     Pmin    Pmax    RU0    RD0
g1     0.12  14.80  89     1.2    -5      3     28      200     40     40
g2     0.17  16.57  83     2.3    -4.24   6.09  20      290     30     30
g3     0.15  15.55  100    1.1    -2.15   5.69  30      190     30     30
g4     0.19  16.21  70     1.1    -3.99   6.2   20      260     50     50 ;
Parameter demand(t)
/t1  510,t2  530,t3  516,t4  510,t5  515,t6  544,
 t7  646,t8  686,t9  741,t10 734,t11 748,t12 760,
 t13 754,t14 700,t15 686,t16 720,t17 714,t18 761,
 t19 727,t20 714,t21 618,t22 584,t23 578,t24 544/;

*变量-----------------------------------------------------------------------------------------------
Variables    OBJ            objecet value
             costThermal    cost of thermal units
             p(i,t)         power generated by thermal power plant
             EM             Emission calculation ;
*方程-----------------------------------------------------------------------------------------------
Equations Genconst3,Genconst4,costThermalcalc,balance,EMcalc;
costThermalcalc ..costThermal=e=sum((i,t),gendata(i,'a')*p(i,t)*p(i,t)+gendata(i,'b')*p(i,t)+gendata(i,'c'));
Genconst3(i,t)..p(i,t+1)-p(i,t)=l=gendata(i,'RU0');
Genconst4(i,t)..p(i,t-1)-p(i,t)=l=gendata(i,'RD0');
balance(t)..sum(i,p(i,t))=g=demand(t);
EMcalc..Em=e=sum((t,i),gendata(i,'d')*power(p(i,t),2)+gendata(i,'e')*p(i,t)+gendata(i,'f'));

*模型及输出-----------------------------------------------------------------------------------------
Model DEDcostbased /all/;
Solve DEDcostbased us QCP min costThermal;

p.up(i,t) = gendata(i,'Pmax');
p.lo(i,t) = gendata(i,'Pmin');

execute_unload "DEDcostbased.gdx" P.l;
execute 'gdxxrw.exe DEDcostbased.gdx var=p rng=pthermal!a1';
execute 'gdxxrw.exe DEDcostbased.gdx var=p rng=ptherma2!a7';
display          costThermal.l , Em.l;

4，难点

    看看就行了，没什么难点。

 

---

搜索“GAMS系列分享”，查看GAMS在电力系统（综合能源）的应用！！！

搜索“GAMS系列分享”，查看GAMS在电力系统（综合能源）的应用！！！

搜索“GAMS系列分享”，查看GAMS在电力系统（综合能源）的应用！！！