[發(fā)明專利]一種基于供熱相量模型的電-熱多能流系統(tǒng)優(yōu)化調(diào)度方法

有效

• 申請?zhí)枺?/td>: 202010090172.8
• 申請日: 2020-02-13
• 公開（公告）號: CN111340271B
• 公開（公告）日: 2022-04-08
• 發(fā)明（設(shè)計）人: 孫宏斌;郭慶來;王彬;陳瑜瑋;潘昭光
• 申請（專利權(quán)）人: 清華大學(xué)
• 主分類號: G06Q10/04
• 分類號: G06Q10/04;G06Q50/06
• 代理公司: 北京清亦華知識產(chǎn)權(quán)代理事務(wù)所(普通合伙) 11201
• 代理人: 羅文群

1.一種基于供熱相量模型的電-熱多能流系統(tǒng)優(yōu)化調(diào)度方法，其特征在于包括以下步驟：

(1)將電-熱多能流系統(tǒng)中供熱系統(tǒng)的負荷功率轉(zhuǎn)換成相量形式，其表達式如下：

其中，上標HL為熱負荷標識，為當供熱系統(tǒng)中第i個熱負荷的頻率為kΩ時的用熱功率相量，為供熱系統(tǒng)中第i個熱負荷在調(diào)度時刻τ_w的負荷功率，Ω為相量的基波頻率，K為相量基波頻率的階數(shù)，K的取值等于調(diào)度時刻數(shù)，并滿足KΩ＝24小時，Δτ為調(diào)度時間間隔；

(2)設(shè)定相量形式的電-熱多能流系統(tǒng)中供熱系統(tǒng)約束條件，包括：

(2-1)相量形式的供熱系統(tǒng)熱網(wǎng)管道熱量損失約束方程如下：

其中，B^H為供熱系統(tǒng)管道編號集合，為供熱系統(tǒng)中第b條管道的頻率為kΩ的首端溫度相量，為供熱系統(tǒng)中第b條管道的頻率為kΩ的末端溫度相量，t_am為供熱系統(tǒng)環(huán)境溫度，m_b為第b條管道的流量，L_b為第b條管道的長度，C_w為水的比熱容，比熱容的取值為4182焦耳/(千克·攝氏度)，λ_b為第b條管道單位長度的傳熱系數(shù)，λ_b從電-熱耦合多能流系統(tǒng)的能量管理系統(tǒng)中獲取；

(2-2)相量形式的供熱系統(tǒng)的多管道匯合點的溫度約束方程如下：

其中，上標HS為供熱系統(tǒng)中熱源的標識，為供熱系統(tǒng)中第i個熱源的頻率為kΩ的供熱功率相量，為供熱系統(tǒng)中第b條管道的頻率為kΩ的末端溫度相量，C_w為水的比熱容，比熱容的取值為4182焦耳/(千克·攝氏度)，m_b為第b條管道的流量，為供熱系統(tǒng)中第i個熱負荷的頻率為kΩ的用熱功率相量，為供熱系統(tǒng)中第n個節(jié)點的頻率為kΩ的溫度相量，為供熱系統(tǒng)熱源節(jié)點編號的集合，為供熱系統(tǒng)熱負荷節(jié)點編號的集合，為供熱系統(tǒng)中末端節(jié)點是第n個節(jié)點的管道編號的集合，為供熱系統(tǒng)中首端節(jié)點是第n個節(jié)點的管道編號的集合；

(2-3)相量形式的供熱系統(tǒng)的管道首端溫度約束方程如下：

其中，為供熱系統(tǒng)中第b條管道的頻率為kΩ的首端溫度相量，為供熱系統(tǒng)中第n個節(jié)點的頻率為kΩ的溫度相量，為供熱系統(tǒng)中首端節(jié)點是第n個節(jié)點的管道編號的集合，N^H為供熱系統(tǒng)節(jié)點編號集合；

(2-4)供熱系統(tǒng)熱源相量傅里葉反變換約束方程如下：

其中，上標HS為供熱系統(tǒng)中熱源標識，為供熱系統(tǒng)中第i個熱源在調(diào)度時刻τ_w的供熱功率，函數(shù)Re(·)表示對復(fù)數(shù)取實部，為供熱系統(tǒng)中第i個熱源的頻率為kΩ的供熱功率相量，I^HS為供熱系統(tǒng)熱源編號的集；

(2-5)供熱系統(tǒng)節(jié)點溫度歷史數(shù)據(jù)約束方程如下：

其中，為供熱系統(tǒng)中第n個節(jié)點在歷史調(diào)度時刻τ_w,his的節(jié)點溫度，函數(shù)Re(·)表示對復(fù)數(shù)取實部，為供熱系統(tǒng)中第n個節(jié)點的頻率為kΩ的溫度相量，N^H為供熱系統(tǒng)節(jié)點編號集合；

(2-6)供熱系統(tǒng)節(jié)點溫度邊界約束方程如下：

其中，t _n 為供熱系統(tǒng)中第n個節(jié)點的溫度下限，為供熱系統(tǒng)中第n個節(jié)點的溫度上限，函數(shù)Re(·)表示對復(fù)數(shù)取實部，為供熱系統(tǒng)中第n個節(jié)點的頻率為kΩ的溫度相量，N^H為供熱系統(tǒng)節(jié)點編號集合；

(2-7)供熱系統(tǒng)中熱電聯(lián)產(chǎn)機組約束方程如下：

其中，上標ES為電源標識，為供熱系統(tǒng)中第i個熱電聯(lián)產(chǎn)機組調(diào)度時刻τ_w的發(fā)電功率，為供熱系統(tǒng)中第i個熱電聯(lián)產(chǎn)機組調(diào)度時刻τ_w的供熱功率，P_κ,i為第i個熱電聯(lián)供機組運行可行域近似多邊形的第κ個頂點的橫坐標，Q_κ,i為第i個熱電聯(lián)產(chǎn)機組運行可行域近似多邊形的第κ個頂點的縱坐標，為第i個熱電聯(lián)供機組在調(diào)度時刻τ_w的第κ個合系數(shù)，NK_i為第i個熱電聯(lián)供機組的運行可行域近似多邊形的頂點個數(shù)，上述構(gòu)成熱電聯(lián)供機組運行可行域近似多邊形從熱電聯(lián)供機組的出廠說明書中獲取，I^CHP為供熱系統(tǒng)中熱電聯(lián)產(chǎn)機組的編號集合；

(3)設(shè)定電-熱多能流系統(tǒng)中電力系統(tǒng)的約束條件，包括：

(3-1)電力系統(tǒng)直流潮流約束方程如下：

其中，上標ES為電源標識，為電力系統(tǒng)中第i個發(fā)電機組在調(diào)度時刻τ_w的發(fā)電功率，為電力系統(tǒng)中第n個節(jié)點在調(diào)度時刻τ_w的電負荷功率，I^ES為電力系統(tǒng)發(fā)電機組編號的集合，N^E為電力系統(tǒng)節(jié)點編號集合；F_b為電力系統(tǒng)中第b條線路的功率上限，Φ_b,n為電力系統(tǒng)中第n個節(jié)點和第b條線路間的轉(zhuǎn)移分布因子，為電力系統(tǒng)中第n個節(jié)點上發(fā)電機組的集合，B^E為電力系統(tǒng)的線路集合；

(3-2)電力系統(tǒng)中發(fā)電機組的約束方程如下：

其中，上標TU為電力系統(tǒng)中除熱電聯(lián)產(chǎn)機組以外的其他發(fā)電機組的標識，p_i為電力系統(tǒng)中第i個發(fā)電機組的功率下限，為電力系統(tǒng)中第i個發(fā)電機組的功率上限，為電力系統(tǒng)中第i個發(fā)電機組在調(diào)度時刻τ_w的發(fā)電功率，I^TU為電力系統(tǒng)中發(fā)電機組的編號集合；

(4)建立一個電-熱多能流系統(tǒng)優(yōu)化調(diào)度的目標函數(shù)如下：

其中，為供熱系統(tǒng)中第i個熱電聯(lián)產(chǎn)機組在調(diào)度時刻τ_w的運行成本，為電力系統(tǒng)第i個發(fā)電機組在調(diào)度時刻τ_w的運行成本，I^CHP為供熱系統(tǒng)中熱電聯(lián)產(chǎn)機組的編號集合，I^TU為電力系統(tǒng)中發(fā)電機組的編號集合，和的具體表達式如下：

其中，a_0,i、a_1,i、a_2,i、a_3,i、a_4,i和a_5,i為第i個熱電聯(lián)產(chǎn)機組/發(fā)電機組的成本系數(shù)，a_0,i、a_1,i、a_2,i、a_3,i、a_4,i和a_5,i從電-熱多能流系統(tǒng)的能量管理系統(tǒng)中獲取，為第i個熱電聯(lián)產(chǎn)機組或發(fā)電機組在調(diào)度時刻τ_w的發(fā)電功率，為第i個熱電聯(lián)產(chǎn)機組在調(diào)度時刻τ_w的供熱功率；

(5)采用內(nèi)點法，求解由步驟(4)目標函數(shù)和步驟(2)、步驟(3)的約束條件組成的優(yōu)化模型，得到電-熱多能流系統(tǒng)中發(fā)電機組的發(fā)電功率、熱電聯(lián)產(chǎn)機組的發(fā)電功率和供熱功率，作為電-熱多能流系統(tǒng)的優(yōu)化調(diào)度參數(shù)，實現(xiàn)基于供熱相量模型的電-熱多能流系統(tǒng)優(yōu)化調(diào)度。