GAMS系列分享13——综合能源系统——包含储能的单能源枢纽模型
本节介绍一个包含储能的 能源枢纽模型的日前优化调度。
电储能
24个时刻
CHP
原来也写过综合能源例子的分享，但是编程格式什么的很乱，可读性差。
。。。。。。。（l）
现在比原来强了呦 ,
过了很久读写的例子，
可读性还是很强的，
前提，对着模型=读程序
ding ding lai le
目录
1，物理模型
2，数学模型
3，代码
3，求解结果
4，难点分析
5，思路点拨
1，物理模型
2，数学模型
式10.12a ，目标函数保证系统运行成本最低
式10.12b， 满足电负荷需求
式10.12e-i，电储能约束
式10.12k, 满足热负荷
式10.12m, 满足冷负荷需求
其他约束较简单，不再介绍。
3，代码
根据建立的数学模型，编写gams程序，如下。
绝对可以运行的！！！！
$ontext
EC 能源枢纽24小时经济调度。
$offtext
*集合
set t /t1*t24/
*数据-----------------------------------------------------------------------------------------------
** 负荷大小
Table data(t,*)
Dh De Dc lamda
t1 21.4 52.1 11.5 36.7
t2 23.2 66.7 13.7 40.4
t3 26.1 72.2 16 38.5
t4 26.7 78.4 21.4 38
t5 25.6 120.2 22 40.2
t6 26.4 83.5 30.8 38.6
t7 39.5 110.4 38.9 52.3
t8 47.3 124.3 46.8 67.3
t9 52.1 143.6 51 70.5
t10 49.1 149.3 48.9 66.2
t11 69.3 154.2 34.8 73.3
t12 62 147.3 32.7 60.8
t13 68 200.7 27.8 63.2
t14 68.6 174.4 32 70.8
t15 56.4 176.5 33.2 63.1
t16 41.3 136.1 34.1 52.5
t17 37.4 108.7 40.8 57
t18 25.4 96.9 43.6 49.2
t19 25.7 89.1 51.5 47.5
t20 21.9 82.5 43.1 49.5
t21 22.4 76.9 36.5 53.1
t22 24.6 66.8 27.7 51.6
t23 22.7 47.2 19.1 50.5
t24 22.6 64.7 11 36.4 ;
**各转换设备的效率
Scalar
eta_ee 变压器效率 /0.98/
eta_ghf 燃气热炉的效率/0.90/
eta_hc 热转冷的效率 /0.95/
eta_ge chp机组的效率 /0.45 /
eta_gh chp机组的效率 /0.35/
eta_c 电池充放电效率 /0.9/
eta_d /0.9/
** 耦合元件的特性，最大最小值，容量，充放电速率等
SOC0 /120/
SOCmax /600/
SOCmin /20 /
EchMin /0/
EdchMin /0 /
Chpmax /250/
CBmax /500/
Fmax /600/ ;
*变量———————————————————————————————————————————————
Variables cost ;
positive variables E(t),E1(t),E2(t),G(t),G1(t),G2(t),SOC(t),H1(t),H2(t),Edch(t);
binary variables idch(t),ich(t);
*方程
equations eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9,eq10,eq11,eq12,eq13;
**目标函数
eq1 .. cost =e=sum(t,data(t,'lamda')*E(t)+12*G(t));
**电负荷，热，冷负荷需求满足
eq2(t) .. eta_ee*E2(t)+Edch(t)+eta_ge*G1(t)=e=data(t,'De');
eq3(t) .. eta_gh*G1(t)+H1(t) =e=data(t,'Dh');
eq4(t) .. eta_hc*H2(t) =e= data(t,'Dc');
**输入端元件
eq5(t) .. E1(t) +E2(t) =e= E(t);
eq6(t) .. G1(t) +G2(t) =e= G(t);
**耦合元件
eq7(t) .. H1(t) +H2(t) =e= eta_ghf*G2(t);
eq8(t) .. SOC(t) =e=SOC0$(ord(t)=1)+SOC(t-1)$(ord(t)>1)+(E1(t)*eta_c-Edch(t)/eta_d);
eq9(t) .. E1(t) =l=ich(t)*0.2*SOCmax;
eq10(t) .. E1(t) =g=ich(t)*EchMin;
eq11(t) .. Edch(t)=l=idch(t)*0.2*SOCmax;
eq12(t) .. Edch(t)=g=idch(t)*EchMin;
eq13(t) .. idch(t)+ich(t)=l= 1;
*模型创建和求解——————————————————————————————————————————
model hub2 /all/
solve hub2 us mip min cost ;
*变量约束
G1.up(t) = Chpmax;
G2.up(t) = Fmax;
H2.up(t) = CBmax;
SOC.fx('t24') = SOC0;
*输出结果
display cost.l,g.l,e.l;
3，求解结果
4，难点分析
eq8(t) .. SOC(t) =e=SOC0$(ord(t)=1)+SOC(t-1)$(ord(t)>1)+(E1(t)*eta_c-Edch(t)/eta_d);
主要是储能这一部分要写成这种形式。。
多看两遍就行了。
多看两遍
多看
多
。
。
。
动手写写
多动手写写
初学者动手写写喽
5，思路点拨
如何编程呢？初学者（我自己也算，在这里分享也算督促自己学习）应该这样，我自己在分享的过程中，也整理思路。到底一个程序怎么快速完成
1，首先要有详尽的数学模型。
目标函数，
约束条件
能有多详细就多详细
最好再配个图
2，按照GAMS的格式进行编程
集合：t
数据：负荷数据，耦合元件效率等参数
变量
方程
模型
整体框架写出来，再慢慢填充
3，写不出来就多看例子，一般写不出来，是头脑里面没有清晰的物理模型，或编程没思路。
还需要多看别人的例子。
第一个单能源枢纽模型的例子介绍到这里了。
有不懂的多多交流，一般晚上我可能会回复的。
搜索“GAMS系列分享”，查看GAMS在电力系统（综合能源）中的应用！！！
搜索“GAMS系列分享”，查看GAMS在电力系统（综合能源）中的应用！！！
搜索“GAMS系列分享”，查看GAMS在电力系统（综合能源）中的应用！！！