代怡重,王曉晶,葛寬寬,肖 群,梁郭江,張紫薇
(1.新疆大學(xué)電氣工程學(xué)院,新疆 烏魯木齊 830047;2.國(guó)網(wǎng)石嘴山供電公司,寧夏 石嘴山 753000)
近年來(lái),風(fēng)能憑借其清潔、成本低、取用不盡等特點(diǎn),使得風(fēng)力發(fā)電日益受到重視。在我國(guó)“三北”地區(qū)有著大量的風(fēng)電資源,然而,在冬季供暖期為滿足供熱需求,熱電機(jī)組需要“以熱定電”的方式運(yùn)行,導(dǎo)致系統(tǒng)的調(diào)峰能力下降從而造成大量的棄風(fēng)[1-3];同時(shí),風(fēng)電的反調(diào)峰特性也是風(fēng)電難以消納的重要原因[4-6]。根據(jù)國(guó)家能源局發(fā)布的2017年風(fēng)電數(shù)據(jù)顯示我國(guó)年平均棄風(fēng)率達(dá)到了12%,其中甘肅、新疆、吉林、內(nèi)蒙古分別棄風(fēng)33%、29%、21%、15%[7]?!叭薄钡貐^(qū)的嚴(yán)重棄風(fēng)已經(jīng)影響到了我國(guó)風(fēng)電的發(fā)展速度。
國(guó)內(nèi)外學(xué)者當(dāng)前主要從熱電解耦、調(diào)節(jié)風(fēng)電出力特性、用戶需求響應(yīng)等方面研究風(fēng)電消納問(wèn)題。文獻(xiàn)[8]通過(guò)儲(chǔ)能裝置的時(shí)移特性調(diào)節(jié)風(fēng)電出力特性,增加夜間的風(fēng)電消納水平,同時(shí)利用儲(chǔ)能系統(tǒng)建立跨區(qū)域風(fēng)電消納模型,提升風(fēng)電跨區(qū)域消納水平。文獻(xiàn)[9]為解耦電熱耦合關(guān)系,提出一種電熱聯(lián)合儲(chǔ)能方法,利用儲(chǔ)能裝置的時(shí)移特性和電、熱系統(tǒng)耦合關(guān)系,促進(jìn)風(fēng)電消納。并引入棄風(fēng)懲罰,建立以系統(tǒng)運(yùn)行成本最小為目標(biāo)的綜合調(diào)度模型。文獻(xiàn)[10]建立了多目標(biāo)電熱綜合能源系統(tǒng)低碳經(jīng)濟(jì)調(diào)度模型,考慮需求側(cè)響應(yīng),通過(guò)分時(shí)電價(jià)優(yōu)化負(fù)荷需求曲線,有效削減負(fù)荷峰谷差,增大風(fēng)電上網(wǎng)空間。文獻(xiàn)[11]引入用戶舒適度,把熱負(fù)荷需求轉(zhuǎn)變?yōu)榭勺儏^(qū)間,作為系統(tǒng)運(yùn)行約束,在滿足用戶舒適度的同時(shí)調(diào)整熱電聯(lián)產(chǎn)機(jī)組出力,促進(jìn)風(fēng)電消納。
綜上可知,熱電聯(lián)合儲(chǔ)能系統(tǒng)、考慮用戶需求響應(yīng)優(yōu)化電負(fù)荷需求曲線、考慮用戶熱舒適度優(yōu)化熱負(fù)荷需求曲線的方法均能有效提高風(fēng)電消納。本文在以上方法的基礎(chǔ)上,首先建立了綜合考慮用戶用電舒適度和用電經(jīng)濟(jì)度的用電滿意度模型,然后根據(jù)用戶的熱舒適度指標(biāo)建立了熱用戶的熱滿意度模型,最后以用戶的滿意度最大化和系統(tǒng)運(yùn)行成本最小為目標(biāo)函數(shù),建立考慮用戶滿意度的熱電儲(chǔ)系統(tǒng)日前優(yōu)化調(diào)度模型。通過(guò)算例驗(yàn)證了該模型能夠進(jìn)一步提高系統(tǒng)的風(fēng)電接納空間,促進(jìn)風(fēng)電消納,降低系統(tǒng)的煤耗成本,同時(shí)提高用戶用電的經(jīng)濟(jì)性。
熱電儲(chǔ)系統(tǒng)由熱電廠、火電廠、風(fēng)電場(chǎng)、儲(chǔ)電裝置、電鍋爐、儲(chǔ)熱裝置以及電熱傳輸網(wǎng)絡(luò)構(gòu)成,熱電儲(chǔ)系統(tǒng)結(jié)構(gòu)如圖1所示。熱電儲(chǔ)系統(tǒng)通過(guò)儲(chǔ)電裝置的時(shí)移特性改變風(fēng)電出力,把負(fù)荷低谷期的棄風(fēng)電量?jī)?chǔ)存起來(lái),在負(fù)荷高峰期供給用戶,從而促進(jìn)風(fēng)電消納;通過(guò)儲(chǔ)熱裝置和電鍋爐改變熱負(fù)荷等效需求,在風(fēng)電出力高峰、熱負(fù)荷需求低谷期,電鍋爐消納棄風(fēng)把電功率轉(zhuǎn)換成熱功率供給用戶,在配合熱電聯(lián)產(chǎn)機(jī)組滿足供熱的情況下,把剩余熱量用儲(chǔ)熱裝置儲(chǔ)存起來(lái),在風(fēng)電低谷期供給用戶,降低熱電聯(lián)產(chǎn)機(jī)組的熱出力,解耦機(jī)組的電熱耦合關(guān)系,從而提高風(fēng)電的接納空間。
圖1 熱電儲(chǔ)系統(tǒng)結(jié)構(gòu)示意
電用戶為降低用電的費(fèi)用響應(yīng)分時(shí)電價(jià),從而改變用電習(xí)慣,降低用電的舒適度,因此定義用戶的用電滿意度為
θ1=δε1+(1-δ)ε2
(1)
式中,θ1為用戶用電滿意度;δ為用戶用電經(jīng)濟(jì)度權(quán)重;ε1為用戶用電經(jīng)濟(jì)度;ε2為用戶用電舒適度。
2.1.1 用電經(jīng)濟(jì)度
用戶響應(yīng)分時(shí)電價(jià),改變用電習(xí)慣,使得用電費(fèi)用發(fā)生改變。定義用戶需求響應(yīng)前后,用戶用電費(fèi)用的改變比為用戶的用電經(jīng)濟(jì)度。
(2)
其中,用戶響應(yīng)分時(shí)電價(jià)后的電負(fù)荷需求PL,t與分時(shí)電價(jià)的關(guān)系為
(3)
2.1.2 用電舒適度
電用戶在未響應(yīng)分時(shí)電價(jià)前,可以隨意根據(jù)自己的需求使用電,此時(shí)用戶的用電舒適度最高,響應(yīng)分時(shí)電價(jià)之后,用戶改變用電習(xí)慣,把如洗衣機(jī)、電動(dòng)車(chē)、電熱水器等柔性負(fù)荷轉(zhuǎn)移到谷時(shí)電價(jià)時(shí)段,從而降低用電舒適度。定義用戶用電舒適度為
(4)
通過(guò)調(diào)節(jié)用戶的熱滿意度,改變熱負(fù)荷的需求,從而降低系統(tǒng)的供熱成本。系統(tǒng)對(duì)用戶的熱滿意度的調(diào)節(jié),首先要保證用戶的熱舒適度。因此要把用戶的熱舒適度作為約束條件。
熱舒適度指標(biāo)PMV是衡量環(huán)境溫度、相對(duì)濕度、空氣流速、穿衣情況等因素對(duì)人體熱感覺(jué)影響程度的指標(biāo)。文獻(xiàn)[13]對(duì)熱舒適度模型進(jìn)行簡(jiǎn)化,僅將室內(nèi)溫度和相對(duì)濕度作為輸入?yún)?shù)
IPMV,t=0.5Tin,t-5.38ρt-8.772
(5)
式中,IPMV,t為t時(shí)刻人體熱舒適度指標(biāo)值;Tin,t為t時(shí)刻室內(nèi)溫度;ρt為t時(shí)刻室內(nèi)空氣相對(duì)濕度,取60%。
根據(jù)式(5)得出人體舒適度指標(biāo)值對(duì)應(yīng)的溫度如表1所示。
表1 PMV與體感、溫度對(duì)應(yīng)關(guān)系
為保證用戶的熱舒適度,我國(guó)規(guī)定IPMV取值在-0.9~0.9之間[14]。因此室溫在22.2 ℃~25.8 ℃之間人體熱舒適度較佳。根據(jù)人體熱舒適度對(duì)應(yīng)的溫度區(qū)間,求解用戶t時(shí)刻的熱滿意度
(6)
式中,ε3,t為用戶t時(shí)刻的熱滿意度;Tin,min為保證用戶熱舒適度的最低室內(nèi)溫度,取22.2 ℃;Tin,fit為最適室內(nèi)溫度,取24 ℃;Tin,max為保證用戶熱舒適度的最高室內(nèi)溫度,取25.8 ℃。
對(duì)系統(tǒng)調(diào)度周期內(nèi)熱用戶總的滿意度進(jìn)行歸一化,得出綜合熱滿意度為
(7)
式中,ε3min,t為t時(shí)刻用戶最小熱滿意度,取0;ε3max,t為t時(shí)刻用戶最大熱滿意度,取1。
室內(nèi)溫度和熱負(fù)荷需求關(guān)系為[15]
(8)
式中,HH,t為t時(shí)刻的熱負(fù)荷需求;S為供熱面積;μ為單位熱量損失,取1.037×104J/(m2℃);C為單位供熱面積下的熱容,取1.63×105J/(m2℃);Tout,t為t時(shí)刻的室外溫度。
考慮用戶滿意度的熱電儲(chǔ)系統(tǒng)日前優(yōu)化調(diào)度模型以用戶的電、熱滿意度最大和熱電儲(chǔ)系統(tǒng)的運(yùn)行成本最小構(gòu)建目標(biāo)函數(shù)
maxf1=θ1=δε1+(1-δ)ε2
(9)
(10)
式中,f1為用戶用電滿意度;f2為用戶熱滿意度;熱電儲(chǔ)系統(tǒng)的運(yùn)行成本包括火電機(jī)組發(fā)電成本、熱電聯(lián)產(chǎn)機(jī)組發(fā)電成本、儲(chǔ)電裝置運(yùn)維成本、電鍋爐運(yùn)維成本、儲(chǔ)熱裝置運(yùn)維成本以及棄風(fēng)懲罰成本。
(11)
(12)
綜上,系統(tǒng)的目標(biāo)函數(shù)為
maxF=ω1f1+ω2f2+ω3f3
(13)
式中,ω1、ω2、ω3為權(quán)重,且ω1+ω2+ω3=1。
本文研究系統(tǒng)的日前優(yōu)化調(diào)度,而機(jī)組的啟停時(shí)間比較長(zhǎng),因此不考慮機(jī)組的啟停成本。火電機(jī)組的發(fā)電成本僅考慮火電機(jī)組的煤耗成本,發(fā)電成本公式為
(14)
式中,N為火電機(jī)組臺(tái)數(shù);αi、βi、χi為第i臺(tái)火電機(jī)組的煤耗系數(shù);Ph,i,t為第i臺(tái)火電機(jī)組t時(shí)刻出力。
抽氣式熱電聯(lián)產(chǎn)機(jī)組純凝工況下的電功率為
Pe,j,t=Pchp,j,t+cvHchp,j,t
(15)
式中,Pchp,j,t為第j臺(tái)熱電聯(lián)產(chǎn)機(jī)組t時(shí)刻的發(fā)電功率;Hchp,j,t為第j臺(tái)熱電聯(lián)產(chǎn)機(jī)組t時(shí)刻的產(chǎn)熱功率;cv為固定進(jìn)氣量下每多發(fā)單位供熱量時(shí)電出力的減小值,取0.15。
熱電聯(lián)產(chǎn)機(jī)組的發(fā)電成本為
(16)
式中,M為熱電聯(lián)產(chǎn)機(jī)組的臺(tái)數(shù);aj、bj、cj為第j臺(tái)熱電機(jī)組的煤耗系數(shù)。
儲(chǔ)電裝置運(yùn)維成本為
Ccd=keQdcf
(17)
式中,ke為單位能量的運(yùn)維成本;Qdcf為運(yùn)行周期內(nèi)儲(chǔ)電裝置放電量。
電鍋爐運(yùn)行成本為
(18)
式中,kdgl為電鍋爐單位用電量的維護(hù)成本;Pdgl,t為電鍋爐t時(shí)刻的用電功率。
儲(chǔ)熱裝置運(yùn)行成本為
Ccr=kcrQcr
(19)
式中,kcr為儲(chǔ)熱裝置單位維護(hù)成本;Qcr為儲(chǔ)熱裝置運(yùn)行周期的儲(chǔ)熱量。
為促進(jìn)風(fēng)電的消納,引入棄風(fēng)懲罰機(jī)制,未消納的風(fēng)電要被懲罰。棄風(fēng)懲罰成本為
(20)
3.3.1 火電機(jī)組運(yùn)行約束
純凝工況下的機(jī)組出力約束為
Ph,i,min≤Ph,i,t≤Ph,i,max
(21)
式中,Ph,i,max、Ph,i,min分別為火電機(jī)組i的最大、最小發(fā)電功率。
機(jī)組爬坡約束
-ΔPh,i,down≤Ph,i,t-Ph,i,t-1≤ΔPh,i,up
(22)
式中,ΔPh,i,down為火電機(jī)組i向下爬坡的最大爬坡率;ΔPh,i,up為火電機(jī)組i向上爬坡的最大爬坡率。
機(jī)組旋轉(zhuǎn)備用約束
(23)
式中,Uu、Ud分別為常規(guī)火電機(jī)組的正、負(fù)旋轉(zhuǎn)備用;λu、λd為火電機(jī)組的正、負(fù)備用系數(shù);σu、σd分別為風(fēng)電并網(wǎng)的旋轉(zhuǎn)備用系數(shù)。
3.3.2 熱電聯(lián)產(chǎn)機(jī)組運(yùn)行約束
抽氣式熱電聯(lián)產(chǎn)機(jī)組電出力約束
(24)
式中,κj為熱電機(jī)組i電功率與熱功率的彈性比值,取0.75;cv為熱電比,即機(jī)組蒸汽輸入量恒定時(shí),每多發(fā)單位供熱量時(shí)電出力的減小值,取0.15。
機(jī)組熱出力約束
Hchp,j,min≤Hchp,j,t≤Hchp,j,max
(25)
式中,Hchp,j,min為熱電聯(lián)產(chǎn)機(jī)組最小熱出力;Hchp,j,max為熱電聯(lián)產(chǎn)機(jī)組最大熱出力。
機(jī)組爬坡約束
-ΔHchp,j,down≤Hchp,j,t-Hchp,j,t-1≤ΔHchp,j,up
(26)
式中,ΔHchp,j,down為熱電聯(lián)產(chǎn)機(jī)組j的熱功率向下爬坡的最大爬坡率;ΔHchp,j,up為熱電聯(lián)產(chǎn)機(jī)組j的熱功率向上爬坡的最大爬坡率。
-ΔPchp,j,down≤Pchp,j,t-Pchp,j,t-1≤ΔPchp,j,up
(27)
式中,ΔPchp,j,down為熱電聯(lián)產(chǎn)機(jī)組j的電功率向下爬坡的最大爬坡率;ΔPchp,j,up為熱電聯(lián)產(chǎn)機(jī)組j的電功率向上爬坡的最大爬坡率。
3.3.3 風(fēng)電出力約束
風(fēng)電出力約束為
(28)
3.3.4 電鍋爐運(yùn)行約束
電鍋爐運(yùn)行約束為
(29)
式中,Pdgl,t為電鍋爐t時(shí)刻的功率,為促進(jìn)風(fēng)電消納,電鍋爐只使用風(fēng)電和儲(chǔ)電裝置儲(chǔ)存的風(fēng)電,因此,Pdgl,t要小于Pw,t-Pdc,t;Pdgl,max為電鍋爐最大功率;Pdc,t>0儲(chǔ)電裝置處于充電狀態(tài),Pdc,t=0儲(chǔ)電裝置不工作,Pdc,t<0儲(chǔ)電裝置處于放電狀態(tài);Hdgl,t為電鍋爐t時(shí)刻的產(chǎn)熱功率;ψ1為電熱轉(zhuǎn)換系數(shù),取98%。
3.3.5 儲(chǔ)電裝置運(yùn)行約束
儲(chǔ)電裝置運(yùn)行約束為
(30)
式中,Pfd,max為儲(chǔ)電裝置最大放電速率;Pcd,max為儲(chǔ)電裝置最大充電速率;E0為初始容量;Et為儲(chǔ)電裝置t時(shí)刻的儲(chǔ)電量,同一個(gè)調(diào)度周期內(nèi)充電量要等于放電量,因此初始容量E0要等于T時(shí)刻的容量Et;Et-1為儲(chǔ)電裝置t-1時(shí)刻的儲(chǔ)電量;V1為儲(chǔ)電裝置充、放電效率,本文取92%。
3.3.6 儲(chǔ)熱裝置運(yùn)行約束
儲(chǔ)熱裝置運(yùn)行約束為
(31)
式中,Hfr,max為儲(chǔ)熱裝置最大放熱速率;Hxr,max為儲(chǔ)熱裝置最大儲(chǔ)熱速率;Hcr,t>0,儲(chǔ)熱裝置處于儲(chǔ)熱狀態(tài),Hcr,t=0儲(chǔ)熱裝置不工作,Hcr,t<0儲(chǔ)熱裝置處于放熱狀態(tài);S0為初始容量;St為儲(chǔ)熱裝置t時(shí)刻的蓄熱量,同一個(gè)調(diào)度周期內(nèi)儲(chǔ)熱裝置的儲(chǔ)熱量要等于放熱量,因此S0要等于St;St-1為儲(chǔ)熱裝置t-1時(shí)刻的蓄熱量;V2為儲(chǔ)熱裝置的儲(chǔ)、放熱效率,取92%。
3.3.7 功率平衡約束
電功率平衡公式為
(32)
熱功率平衡公式為
(33)
以北方某地區(qū)的運(yùn)行數(shù)據(jù)為例,冬季典型室外溫度曲線如圖2所示,風(fēng)電預(yù)測(cè)出力曲線如圖3所示。系統(tǒng)含有常規(guī)火電機(jī)組4臺(tái)(機(jī)組具體參數(shù)見(jiàn)表2);熱電聯(lián)產(chǎn)機(jī)組1臺(tái),其最大發(fā)電功率300 MW,最小發(fā)電功率90 MW,最大熱功率400 MW,最小熱功率100 MW,向上爬坡率150 MW·h,向下爬坡率150 MW·h,煤耗系數(shù)aj、bj、cj分別為0.008 13、27.79、515.5;電鍋爐最大功率為300 MW;儲(chǔ)熱裝置容量300 MW,初始容量50 MW,儲(chǔ)放熱功率50 MW·h;儲(chǔ)電裝置容量為300 MW,初始容量50 MW,儲(chǔ)放電功率50 MW·h;棄風(fēng)懲罰為490 元/(MW·h);系統(tǒng)調(diào)度周期為24 h(22∶00~次日22∶00);峰谷平電價(jià)如表3所示,采用分時(shí)電價(jià)前的電價(jià)為450元/(MW·h)。
表2 常規(guī)火電機(jī)組參數(shù)
表3 分時(shí)電價(jià)
圖2 冬季典型室外溫度曲線
圖3 風(fēng)電預(yù)測(cè)曲線
為驗(yàn)證考慮用戶滿意度的電熱儲(chǔ)系統(tǒng)日前優(yōu)化調(diào)度的合理性、有效性,設(shè)置以下3種運(yùn)行情景:①用戶不響應(yīng)分時(shí)電價(jià),不對(duì)用戶熱滿意度進(jìn)行調(diào)節(jié)(即價(jià)格彈性系數(shù)為0,室內(nèi)溫度恒溫24 ℃);②不考慮用戶滿意度,僅以系統(tǒng)運(yùn)行成本最低為目標(biāo)函數(shù);③用戶響應(yīng)分時(shí)電價(jià),以用戶電熱滿意度和系統(tǒng)運(yùn)行成本為優(yōu)化目標(biāo)。
3種情景的電負(fù)荷需求曲線如圖4所示。情景①用戶不響應(yīng)分時(shí)電價(jià),此時(shí)的電負(fù)荷需求曲線為原始負(fù)荷需求,電負(fù)荷峰谷差較大,22∶00~07∶00電負(fù)荷需求較低,而此時(shí)段是風(fēng)電高發(fā)時(shí)期,系統(tǒng)風(fēng)電接納空間很小,造成大量棄風(fēng)。情景②用戶響應(yīng)分時(shí)電價(jià),且不考慮用戶的用電滿意度,此時(shí)用戶最大響應(yīng)分時(shí)電價(jià),調(diào)整用電的時(shí)間段,把負(fù)荷高峰期的部分用電轉(zhuǎn)移到負(fù)荷低谷期的夜間使用,減少負(fù)荷峰谷差,為風(fēng)電高發(fā)期提供更多的風(fēng)電接納空間。情景③用戶響應(yīng)分時(shí)電價(jià),且調(diào)度時(shí)考慮用戶的用電滿意度,在兼顧用戶滿意度時(shí)提高風(fēng)電接納空間。
圖4 3種情景的電負(fù)荷需求曲線
3種情景的熱負(fù)荷需求曲線如圖5所示。情景①的熱負(fù)荷需求為原始熱負(fù)荷需求,此時(shí)用戶的熱滿意度最高為1。情景②、③在保證用戶舒適度較佳的溫度范圍內(nèi),調(diào)節(jié)用戶室內(nèi)溫度,從而改變熱負(fù)荷需求。在夜間風(fēng)電高發(fā)時(shí)期讓室內(nèi)溫度略高于24 ℃,使夜間熱負(fù)荷需求增多,從而增大風(fēng)電的消納空間,白天讓室內(nèi)溫度略低于24 ℃,減少風(fēng)電低谷期的熱負(fù)荷需求,從而降低儲(chǔ)熱裝置的運(yùn)行成本。情景②不考慮用戶熱滿意度,此時(shí)用戶的熱滿意度最低,犧牲用戶的熱滿意度獲得更多的風(fēng)電消納空間。3種情景的風(fēng)電出力曲線如圖6所示,情景①棄風(fēng)最多;情景②不考慮用戶滿意度風(fēng)電消納最多;情景③兼顧用戶滿意度較情景①也提高了風(fēng)電消納量。
圖5 3種情景的熱負(fù)荷需求曲線
圖6 3種情景的風(fēng)電出力曲線
3種情景的系統(tǒng)運(yùn)行結(jié)果如表5所示,情景②的煤耗成本、系統(tǒng)運(yùn)行總成本、用戶用電成本和棄風(fēng)率均是最高的;情景②的機(jī)組煤耗成本、系統(tǒng)運(yùn)行成本、用戶用電成本和棄風(fēng)率均是最低的,但是用戶的熱滿意度最低;情景③用戶的電、熱滿意度均較佳,比情景①發(fā)電側(cè)節(jié)約發(fā)電成本220 390元,為用戶側(cè)減少用電成本238 630元。
表5 3種情景仿真結(jié)果
針對(duì)“三北”地區(qū)大規(guī)模棄風(fēng)問(wèn)題,本文提出一種考慮用戶滿意度的熱電儲(chǔ)系統(tǒng)日前優(yōu)化調(diào)度方法。用戶滿意度包含用戶用電滿意度和用熱滿意度,其中用電滿意度又包含了用電舒適度和用電經(jīng)濟(jì)度,用戶通過(guò)響應(yīng)分時(shí)電價(jià),犧牲部分用電舒適度,減少用電成本提高了用電經(jīng)濟(jì)度,同時(shí)降低了用電負(fù)荷的峰谷差,增大了風(fēng)電上網(wǎng)空間;用熱滿意度根據(jù)熱舒適度指標(biāo),在保證用戶熱舒適的前提下,把室內(nèi)溫度的舒適區(qū)間轉(zhuǎn)換成用戶滿意度區(qū)間,通過(guò)適當(dāng)?shù)奶岣呋蛘呓档褪覂?nèi)溫度,從而提高風(fēng)電的接納空間,同時(shí)可以降低系統(tǒng)的煤耗成本。算例仿真結(jié)果表明:考慮用戶滿意度時(shí)熱電儲(chǔ)系統(tǒng)的棄風(fēng)率和運(yùn)行成本明顯降低,此時(shí)用戶的滿意度也處于較高的水平,并為電力系統(tǒng)日前調(diào)度提供數(shù)據(jù)參考。