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程序,如下。
? ? ? 絕對可以運行的!!!!
$ontextEC 能源樞紐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在電力系統(綜合能源)中的應用!!!
?
?
總結
以上是生活随笔為你收集整理的GAMS系列分享13——综合能源系统——包含储能的单能源枢纽模型的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: mm_struct
- 下一篇: GDOI2017小结