GAMS系列分享13——综合能源系统——包含储能的单能源枢纽模型
? ? ?本節(jié)介紹一個包含儲能的 能源樞紐模型的日前優(yōu)化調(diào)度。
? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?電儲能
? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 24個時(shí)刻
? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?CHP
原來也寫過綜合能源例子的分享，但是編程格式什么的很亂，可讀性差。
。。。。。。。（l）
現(xiàn)在比原來強(qiáng)了呦?,
? ? ? ? ? 過了很久讀寫的例子，
? ? ? ? ? ? ? ? ? ? ?可讀性還是很強(qiáng)的，
? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 前提，對著模型=讀程序
? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ding? ding? lai? le
目錄
1，物理模型
2，數(shù)學(xué)模型
3，代碼
3，求解結(jié)果
4，難點(diǎn)分析
5，思路點(diǎn)撥
?
1，物理模型
? ? ? ? ? ? ? ? ? ??
2，數(shù)學(xué)模型
? ? ? ??
? ? ? ? ? ? ? ? ? ? ? ? ???
? ?式10.12a ，目標(biāo)函數(shù)保證系統(tǒng)運(yùn)行成本最低
? ?式10.12b，? 滿足電負(fù)荷需求
? ?式10.12e-i，電儲能約束
? ?式10.12k,? ? ?滿足熱負(fù)荷
? ? 式10.12m,? ?滿足冷負(fù)荷需求
? ?其他約束較簡單，不再介紹。
??
3，代碼
? ? ?根據(jù)建立的數(shù)學(xué)模型，編寫gams程序，如下。
? ? ? 絕對可以運(yùn)行的！！！！
$ontextEC 能源樞紐24小時(shí)經(jīng)濟(jì)調(diào)度。 $offtext *集合 set t /t1*t24/ *數(shù)據(jù)----------------------------------------------------------------------------------------------- ** 負(fù)荷大小 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 ; **各轉(zhuǎn)換設(shè)備的效率 Scalar eta_ee 變壓器效率 /0.98/ eta_ghf 燃?xì)鉄釥t的效率/0.90/ eta_hc 熱轉(zhuǎn)冷的效率 /0.95/ eta_ge chp機(jī)組的效率 /0.45 / eta_gh chp機(jī)組的效率 /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; **目標(biāo)函數(shù) eq1 .. cost =e=sum(t,data(t,'lamda')*E(t)+12*G(t)); **電負(fù)荷，熱，冷負(fù)荷需求滿足 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;*模型創(chuàng)建和求解—————————————————————————————————————————— 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;*輸出結(jié)果 display cost.l,g.l,e.l;?
3，求解結(jié)果
? ? ? ? ? ? ? ? ? ? ? ? ? ?? ? ? ?
4，難點(diǎn)分析
eq8(t) ? .. ?SOC(t) =e=SOC0$(ord(t)=1)+SOC(t-1)$(ord(t)>1)+(E1(t)*eta_c-Edch(t)/eta_d);?
主要是儲能這一部分要寫成這種形式。。
多看兩遍就行了。
多看兩遍
多看
多
。
。
。
動手寫寫
多動手寫寫
初學(xué)者動手寫寫嘍
5，思路點(diǎn)撥
? ? ? ?如何編程呢？初學(xué)者（我自己也算，在這里分享也算督促自己學(xué)習(xí)）應(yīng)該這樣，我自己在分享的過程中，也整理思路。到底一個程序怎么快速完成
? ? 1，首先要有詳盡的數(shù)學(xué)模型。
? ? ? ? ? ?目標(biāo)函數(shù)，
? ? ? ? ? ?約束條件
? ? ? ? ? ?能有多詳細(xì)就多詳細(xì)
? ? ? ? ? ?最好再配個圖
? ? 2，按照GAMS的格式進(jìn)行編程
? ? ? ? ? 集合：t
? ? ? ? ? ?數(shù)據(jù)：負(fù)荷數(shù)據(jù)，耦合元件效率等參數(shù)
? ? ? ? ? ?變量
? ? ? ? ? ?方程
? ? ? ? ? ? 模型
? ? ? ?整體框架寫出來，再慢慢填充
? ? ? ?3，寫不出來就多看例子，一般寫不出來，是頭腦里面沒有清晰的物理模型，或編程沒思路。
? ? ? ? ? ?還需要多看別人的例子。
?
?第一個單能源樞紐模型的例子介紹到這里了。
有不懂的多多交流，一般晚上我可能會回復(fù)的。
?
搜索“GAMS系列分享”，查看GAMS在電力系統(tǒng)（綜合能源）中的應(yīng)用！！！
搜索“GAMS系列分享”，查看GAMS在電力系統(tǒng)（綜合能源）中的應(yīng)用！！！
搜索“GAMS系列分享”，查看GAMS在電力系統(tǒng)（綜合能源）中的應(yīng)用！！！
?
?
總結(jié)
以上是生活随笔為你收集整理的GAMS系列分享13——综合能源系统——包含储能的单能源枢纽模型的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。
- 上一篇： mm_struct
- 下一篇： windows 20003 扩展安装后不