一種電-氣耦合多能流網(wǎng)絡狀態(tài)估計方法
【專利摘要】本發(fā)明涉及一種電?氣耦合多能流網(wǎng)絡估計方法，屬于綜合能源系統(tǒng)的運行和控制技術領域。本方法彌補了氣網(wǎng)在狀態(tài)估計方面的空白，形成了一種考慮電網(wǎng)和氣網(wǎng)間耦合關系的電?氣耦合多能流狀態(tài)估計。與單獨狀態(tài)估計相比，本方法有利于提高網(wǎng)絡的估計效果，尤其在電?氣耦合多能流網(wǎng)絡端口量測不準確的情況下，本方法具有明顯的優(yōu)越性。目前氣網(wǎng)的自動化水平尚且不能滿足本方法的需求，本方法的應用需要在氣網(wǎng)安裝更多的量測設施。成功在電?氣耦合多能流網(wǎng)絡應用本方法后，可以得到精度更高的全局一致基態(tài)潮流解，為后續(xù)的電?氣耦合多能流網(wǎng)絡在線安全評估、優(yōu)化調度打下基礎。
【專利說明】
一種電-氣耦合多能流網(wǎng)絡狀態(tài)估計方法
技術領域
[0001] 本發(fā)明涉及一種電-氣耦合多能流網(wǎng)絡估計方法，屬于綜合能源系統(tǒng)的運行和控 制技術領域。
【背景技術】
[0002] 由于燃氣發(fā)電機的成本較低，對環(huán)境的影響較小，燃氣站的建造周期較短等優(yōu)勢， 在世界范圍內，天然氣應用中用于發(fā)電的比例逐漸上升，在電力系統(tǒng)能源供給中也占據(jù)了 越來越重要的地位，因此，電力網(wǎng)絡與天然氣網(wǎng)絡之間的耦合程度逐漸加深，相互之間的影 響與依賴性也愈來愈強。鑒于能源資源條件的限制，我國的能源消費結構在未來相當長的 一段時間內仍然以煤炭為主，但是環(huán)境保護及可持續(xù)發(fā)展戰(zhàn)略都要求我國必須加快改善能 源結構的步伐。在一次能源消費結構中，提高天然氣等可再生的清潔能源的比重，而降低煤 炭的比重。
[0003] 目前，對于電-氣耦合多能流網(wǎng)絡的建模和潮流計算已經(jīng)有了一系列的研究成果。 然而，由于氣網(wǎng)自動化水平低，實現(xiàn)像電網(wǎng)一樣自動化運行還有很長一段距離。為實現(xiàn)電-氣耦合多能流網(wǎng)絡的自動化，狀態(tài)估計是一項基礎技術，而氣網(wǎng)在狀態(tài)估計方面的研究還 非常空白。
[0004] 狀態(tài)估計可以充分利用系統(tǒng)的量測冗余，更準確地了解網(wǎng)絡的運行狀態(tài)，是系統(tǒng) 運行中的基礎部分。在電力系統(tǒng)中，很早就已經(jīng)引入了狀態(tài)估計，目前有關電力系統(tǒng)狀態(tài)估 計的研究已趨于成熟，不論是計算速度還是壞數(shù)據(jù)辨識等方面都有相應的解決措施。對電-氣耦合多能流網(wǎng)絡引入電-氣聯(lián)合狀態(tài)估計，可以得到精度更高的全局一致基態(tài)潮流解，為 后續(xù)的電-氣耦合多能流網(wǎng)絡在線安全評估、優(yōu)化調度打下基礎。

【發(fā)明內容】

[0005] 本發(fā)明的目的是提出一種電_氣耦合多能流網(wǎng)絡估計方法，彌補現(xiàn)有領域研究的 空白，實現(xiàn)對電-氣耦合多能流網(wǎng)絡運行狀態(tài)更準確的監(jiān)控，得到的全局一致基態(tài)潮流解， 為后續(xù)的電-氣耦合多能流網(wǎng)絡在線安全評估、優(yōu)化調度打下基礎。
[0006] 本發(fā)明提出的電-氣耦合多能流網(wǎng)絡估計方法，包括以下步驟：
[0007] (1)建立一個電-氣耦合多能流狀態(tài)估計的目標函數(shù)，如下： -? 「弋⑷、、
[0008] mm(( - ,) R { -,,、）） Lzd LW」 LZJ LM~)」
[0009] 其中：Z(5為電網(wǎng)量測值，包括電-氣耦合多能流網(wǎng)絡的電網(wǎng)中第i個節(jié)點的電壓幅 值Vi、有功功率Pi和無功功率Qi以及電網(wǎng)第i個節(jié)點和第j個節(jié)點間支路ij的有功功率Pij和 無功功率QU，上述電網(wǎng)量測值從電-氣耦合多能流網(wǎng)絡的電網(wǎng)的數(shù)據(jù)采集與監(jiān)視控制系統(tǒng) 中采集 ;Zg為氣網(wǎng)量測值，包括電-氣耦合多能流網(wǎng)絡的氣網(wǎng)中第k個節(jié)點的壓強pk和注出流 量L k以及氣網(wǎng)中第k個節(jié)點和第1個節(jié)點間管道kl的流量fkl，上述氣網(wǎng)量測值從電-氣耦合 多能流網(wǎng)絡的氣網(wǎng)的數(shù)據(jù)采集與監(jiān)視控制系統(tǒng)中采集;^為電網(wǎng)狀態(tài)量，ejT，％、 I分別為電網(wǎng)第i個節(jié)點的電壓幅值和相角，Xg為氣網(wǎng)狀態(tài)量，xg=[pk]T，pk為氣網(wǎng)第k個節(jié) 點的壓強，R為量測值的協(xié)方差矩陣，T為矩陣轉置，為電網(wǎng)潮流計算方程： 'd 唞- % + 5".如說-
[0010] Q^V^V^niO^O^-B^om-OS) 4 cos(4 ―沒,.）十 5" sin(6^..-沒,..).） Qli=yr(Bli-h) + sin(6? - - 5 cos(6> -Oj))
[0011]上述電網(wǎng)潮流計算方程中，Gl偽電網(wǎng)節(jié)點導納矩陣Y第i行、第j列元素的實部，Bu 為電網(wǎng)節(jié)點導納矩陣Y第i行、第j列元素的虛部，電網(wǎng)節(jié)點導納矩陣Y從電網(wǎng)調度中心獲取； b。為電網(wǎng)支路的等效對地電納，從電網(wǎng)調度中心獲取，
[0012] hg為氣網(wǎng)潮流計算方程：
[0014] 上述氣網(wǎng)潮流計算方程中
，妒1為第k個節(jié)點 和第1個節(jié)點之間管道kl的管道直徑，Lkl為第k個節(jié)點和第1個節(jié)點之間管道kl的管道長 度，y C為天然氣比重，Ta為天然氣平均溫度，Tn為天然氣標準狀態(tài)下的溫度，取值為288開爾 文，Pn為天然氣標準狀態(tài)下的壓力，Pn的取值為〇. IMpa，Zg為天然氣平均可壓縮系數(shù)，Zg的取 值范圍為〇~1，F(xiàn)為管道內壁的摩擦系數(shù)
》Re為雷諾數(shù)，Re的取值范 圍在3500以上，Ef為管道的效率系數(shù)，Ef的取值為0.92,當（/<-/，）20時，氣網(wǎng)潮流計算方 程中的sgnP (pk，pi) = 1，當(ft2 -內2) < 〇時，氣網(wǎng)潮流計算方程中的sgnP (pk，pi) = -1;
[0015] (2)建立電-氣耦合多能流網(wǎng)絡狀態(tài)估計的約束條件，包括：
[0016] (2-1)電-氣耦合多能流網(wǎng)絡中氣網(wǎng)的壓縮機耗能約束：
[0018] 其中：Pf、Pt分別為氣網(wǎng)中壓縮機的首端節(jié)點和末端節(jié)點的壓強，BHP為壓縮機的能 耗，壓縮機耗能由電動機或燃氣輪機提供，為壓縮機的入口流量，n。為壓縮機的總效率， 取值范圍為〇~l，c k為該壓縮機的多變系數(shù)，與壓縮機中的氣體壓縮過程有關，取值范圍為 1.2 ~1.5;
[0019] 若由電動機提供氣網(wǎng)中的壓縮機耗能，則電動機消耗電-氣耦合多能流網(wǎng)絡中電 網(wǎng)節(jié)點的有功功率滿足以下關系：
[0020] Pcomp = Kc〇mPBHP
[0021] 其中：？。。_為電動機消耗電網(wǎng)節(jié)點的有功功率，1(。。_為壓縮機電能消耗系數(shù)，由壓 縮機的出廠說明書獲得；
[0022]若由燃氣輪機提供氣網(wǎng)中的壓縮機耗能，則將該約束放入步驟（1)中氣網(wǎng)潮流計 算方程中，燃氣輪機消耗電-氣耦合多能流網(wǎng)絡中氣網(wǎng)節(jié)點的燃氣滿足以下關系：
[0023] fc = ac+PcBHP+y c(BHP)2
[0024] 其中：fc為燃氣輪機消耗氣網(wǎng)節(jié)點的燃氣，ac、&和y c分別為氣網(wǎng)中壓縮機的氣體 消耗系數(shù)，由壓縮機的出廠說明書獲得；
[0025] (2-2)電-氣耦合多能流網(wǎng)絡中的燃氣輪機約束：
[0026] PTl，Chfl-CJ;-u，CjTur
[0027] 其中，fTur為燃氣輪機的燃氣流量，PTur為燃氣輪機的有功功率輸出，&、(： 2和(：3分別 為燃氣輪機的燃氣系數(shù)，由燃氣輪機的出廠說明書獲取；
[0028] (3)利用拉格朗日乘子法將步驟(2)中的約束條件和步驟(1)中的目標函數(shù)構成一 個拉格朗日函數(shù)，拉格朗日函數(shù)的構造方法已在數(shù)學界的最優(yōu)化理論中有了詳細說明，利 用牛頓-拉夫遜法，求解該拉格朗日函數(shù)，得到電-氣耦合多能流網(wǎng)絡狀態(tài)估計的結果。
[0029] 本發(fā)明提出的電-氣耦合多能流網(wǎng)絡狀態(tài)估計方法，其特點和效果是:本方法彌補 了氣網(wǎng)在狀態(tài)估計方面的空白，形成了一種考慮電網(wǎng)和氣網(wǎng)間耦合關系的電-氣耦合多能 流狀態(tài)估計。與單獨狀態(tài)估計相比，本方法有利于提高網(wǎng)絡的估計效果，尤其在電-氣耦合 多能流網(wǎng)絡端口量測不準確的情況下，本方法具有明顯的優(yōu)越性。目前氣網(wǎng)的自動化水平 尚且不能滿足本方法的需求，本方法的應用需要在氣網(wǎng)安裝更多的量測設施。成功在電-氣 耦合多能流網(wǎng)絡應用本方法后，可以得到精度更高的全局一致基態(tài)潮流解，為后續(xù)的電-氣 耦合多能流網(wǎng)絡在線安全評估、優(yōu)化調度打下基礎。
【具體實施方式】
[0030] 本發(fā)明提出的電-氣耦合多能流網(wǎng)絡估計方法，包括以下步驟：
[0031 ] (1)建立一個電-氣耦合多能流狀態(tài)估計的目標函數(shù)，如下：
[0032] mm((-,,、）及（一"、）） W L~」lAh)」
[0033] 其中：Z(5為電網(wǎng)量測值，包括電-氣耦合多能流網(wǎng)絡的電網(wǎng)中第i個節(jié)點的電壓幅 值Vi、有功功率Pi和無功功率Qi以及電網(wǎng)第i個節(jié)點和第j個節(jié)點間支路ij的有功功率Pij和 無功功率QU，上述電網(wǎng)量測值從電-氣耦合多能流網(wǎng)絡的電網(wǎng)的數(shù)據(jù)采集與監(jiān)視控制系統(tǒng) 中采集 ;Zg為氣網(wǎng)量測值，包括電-氣耦合多能流網(wǎng)絡的氣網(wǎng)中第k個節(jié)點的壓強pk和注出流 量L k以及氣網(wǎng)中第k個節(jié)點和第1個節(jié)點間管道kl的流量fkl，上述氣網(wǎng)量測值從電-氣耦合 多能流網(wǎng)絡的氣網(wǎng)的數(shù)據(jù)采集與監(jiān)視控制系統(tǒng)中采集;為電網(wǎng)狀態(tài)量，ejT，％、 I分別為電網(wǎng)第i個節(jié)點的電壓幅值和相角，xg為氣網(wǎng)狀態(tài)量，xg=[Pk]T，p k為氣網(wǎng)第k個節(jié) 點的壓強，R為量測值的協(xié)方差矩陣，T為矩陣轉置，為電網(wǎng)潮流計算方程： ^ ^Zcos^-O^ + B^mW-e^)
[0034] cos(^ - ^,)) p" = f (i?" - A )十!//,.(G〃. sin⑷一) - /i". c〇S(6 -弋））
[0035]上述電網(wǎng)潮流計算方程中，Gi偽電網(wǎng)節(jié)點導納矩陣Y第i行、第j列元素的實部， 為電網(wǎng)節(jié)點導納矩陣Y第i行、第j列元素的虛部，電網(wǎng)節(jié)點導納矩陣Y從電網(wǎng)調度中心獲取； b。為電網(wǎng)支路的等效對地電納，從電網(wǎng)調度中心獲取，
[0036] hg為氣網(wǎng)潮流計算方程：
[0038] 上述氣網(wǎng)潮流計算方程中：
，妒1為第k個節(jié)點 和第1個節(jié)點之間管道kl的管道直徑，Lkl為第k個節(jié)點和第1個節(jié)點之間管道kl的管道長度， y C為天然氣比重，Ta為天然氣平均溫度，Tn為天然氣標準狀態(tài)下的溫度，取值為288開爾文， Pn為天然氣標準狀態(tài)下的壓力，Pn的取值為〇. IMpa，Zg為天然氣平均可壓縮系數(shù)，Zg的取值 范圍為〇~1，F(xiàn)為管道內壁的摩擦系數(shù)，
，Re為雷諾數(shù)，Re的取值范圍 在3500以上，Ef為管道的效率系數(shù)，Ef的取值為0.92，當（凡:2 〇時，氣網(wǎng)潮流計算方程 中的sgnP(pk，pi) = 1，當(八-仍2) <_ 〇時，氣網(wǎng)潮流計算方程中的sgnP(pk，pi) = -1;
[0039] (2)建立電-氣耦合多能流網(wǎng)絡狀態(tài)估計的約束條件，包括：
[0040] (2-1)電-氣耦合多能流網(wǎng)絡中氣網(wǎng)的壓縮機耗能約束：
[0042]其中：pf、pt分別為氣網(wǎng)中壓縮機的首端節(jié)點和末端節(jié)點的壓強，BHP為壓縮機的 能耗，壓縮機耗能由電動機或燃氣輪機提供，為壓縮機的入口流量，n。為壓縮機的總效 率，取值范圍為〇~i，ck為該壓縮機的多變系數(shù)，與壓縮機中的氣體壓縮過程有關，取值范 圍為1.2~1.5;
[0043] 若由電動機提供氣網(wǎng)中的壓縮機耗能，則電動機消耗電-氣耦合多能流網(wǎng)絡中電 網(wǎng)節(jié)點的有功功率滿足以下關系：
[0044] Pcomp = Kc〇mPBHP
[0045] 其中：P_PS電動機消耗電網(wǎng)節(jié)點的有功功率，1(。。_為壓縮機電能消耗系數(shù)，由壓 縮機的出廠說明書獲得；
[0046] 若由燃氣輪機提供氣網(wǎng)中的壓縮機耗能，則將該約束放入步驟（1)中氣網(wǎng)潮流計 算方程中，燃氣輪機消耗電-氣耦合多能流網(wǎng)絡中氣網(wǎng)節(jié)點的燃氣滿足以下關系：
[0047] fc = ac+PcBHP+y c(BHP)2
[0048]其中：fc為燃氣輪機消耗氣網(wǎng)節(jié)點的燃氣，ac、&和Y c分別為氣網(wǎng)中壓縮機的氣體 消耗系數(shù)，由壓縮機的出廠說明書獲得；
[0049] (2-2)電-氣耦合多能流網(wǎng)絡中的燃氣輪機約束：
[0050]
[0051] 其中，fTur為燃氣輪機的燃氣流量，PTur為燃氣輪機的有功功率輸出，&、(： 2和(：3分別 為燃氣輪機的燃氣系數(shù)，由燃氣輪機的出廠說明書獲取；
[0052] (3)利用拉格朗日乘子法將步驟(2)中的約束條件和步驟(1)中的目標函數(shù)構成一 個拉格朗日函數(shù)，拉格朗日函數(shù)的構造方法已在數(shù)學界的最優(yōu)化理論中有了詳細說明，利 用牛頓-拉夫遜法，求解該拉格朗日函數(shù)，得到電-氣耦合多能流網(wǎng)絡狀態(tài)估計的結果。
【主權項】
1. 一種電-氣耦合多能流網(wǎng)絡狀態(tài)估計方法，其特征在于該方法包括以下步驟： (1) 建立一個電-氣耦合多能流狀態(tài)估計的目標函數(shù)，如下：其中：Ze為電網(wǎng)量測值，包括電-氣耦合多能流網(wǎng)絡的電網(wǎng)中第i個節(jié)點的電壓幅值Vi、 有功功率Pi和無功功率Qi以及電網(wǎng)第i個節(jié)點和第j個節(jié)點間支路ij的有功功率Pij和無功 功率QU，上述電網(wǎng)量測值從電-氣耦合多能流網(wǎng)絡的電網(wǎng)的數(shù)據(jù)采集與監(jiān)視控制系統(tǒng)中采 集;z g為氣網(wǎng)量測值，包括電-氣耦合多能流網(wǎng)絡的氣網(wǎng)中第k個節(jié)點的壓強pk和注出流量Lk 以及氣網(wǎng)中第k個節(jié)點和第1個節(jié)點間管道kl的流量fkl，上述氣網(wǎng)量測值從電-氣耦合多能 流網(wǎng)絡的氣網(wǎng)的數(shù)據(jù)采集與監(jiān)視控制系統(tǒng)中采集; Xe3為電網(wǎng)狀態(tài)量，Xe3=[Vl 別為電網(wǎng)第i個節(jié)點的電壓幅值和相角，xg為氣網(wǎng)狀態(tài)量，xg=[Pk]T，pk為氣網(wǎng)第k個節(jié)點的 壓強，R為量測值的協(xié)方差矩陣，T為矩陣轉置，為電網(wǎng)潮流計算方程：上述電網(wǎng)潮流計算方程中，Gi偽電網(wǎng)節(jié)點導納矩陣Y第i行、第j列元素的實部幾偽電 網(wǎng)節(jié)點導納矩陣Y第i行、第j列元素的虛部，電網(wǎng)節(jié)點導納矩陣Y從電網(wǎng)調度中心獲取;b。為 電網(wǎng)支路的等效對地電納，從電網(wǎng)調度中心獲取， hg為氣網(wǎng)潮流計算方程：上述氣網(wǎng)潮流計算方程中：Dkl為第k個節(jié)點和第1 個節(jié)點之間管道kl的管道直徑，Lkl為第k個節(jié)點和第1個節(jié)點之間管道kl的管道長度，丫^為 天然氣比重，Ta為天然氣平均溫度，Tn為天然氣標準狀態(tài)下的溫度，取值為288開爾文，？"為 天然氣標準狀態(tài)下的壓力，p n的取值為〇. IMpa，Zg為天然氣平均可壓縮系數(shù)，Zg的取值范圍 為0~1，F(xiàn)為管道內壁的摩擦系勠Re為雷諾數(shù)，Re的取值范圍在 3500以上，Ef為管道的效率系數(shù)，Ef的取值為0.92，當(么2 - 0時，氣網(wǎng)潮流計算方程中 的sgnp (pk，pi) = 1，當) _< _〇時，氣網(wǎng)潮流計算方程中的sgnp (pk，pi) = -1; (2) 建立電-氣耦合多能流網(wǎng)絡狀態(tài)估計的約束條件，包括： (2-1)電-氣耦合多能流網(wǎng)絡中氣網(wǎng)的壓縮機耗能約束：其中：pf、pt分別為氣網(wǎng)中壓縮機的首端節(jié)點和末端節(jié)點的壓強，BHP為壓縮機的能耗， 壓縮機耗能由電動機或燃氣輪機提供，於為壓縮機的入口流量，η。為壓縮機的總效率，取值 范圍為0~l，Ck為該壓縮機的多變系數(shù)，與壓縮機中的氣體壓縮過程有關，取值范圍為1.2 ~1 · 5; 若由電動機提供氣網(wǎng)中的壓縮機耗能，則電動機消耗電-氣耦合多能流網(wǎng)絡中電網(wǎng)節(jié) 點的有功功率滿足以下關系： Pcomp 一 KcompBHP 其中：p_PS電動機消耗電網(wǎng)節(jié)點的有功功率，1(。。_為壓縮機電能消耗系數(shù)，由壓縮機 的出廠說明書獲得； 若由燃氣輪機提供氣網(wǎng)中的壓縮機耗能，則將該約束放入步驟（1)中氣網(wǎng)潮流計算方 程中，燃氣輪機消耗電-氣耦合多能流網(wǎng)絡中氣網(wǎng)節(jié)點的燃氣滿足以下關系： fc = ac+0cBHP+yc(BHP)2 其中：fc為燃氣輪機消耗氣網(wǎng)節(jié)點的燃氣，ac、說和γ c分別為氣網(wǎng)中壓縮機的氣體消耗 系數(shù)，由壓縮機的出廠說明書獲得； (2 - 2)電-氣耦合多能流網(wǎng)絡中的燃氣輪機約束：其中，fTUr為燃氣輪機的燃氣流量，PTur為燃氣輪機的有功功率輸出，&、(： 2和(：3分別為燃 氣輪機的燃氣系數(shù)，由燃氣輪機的出廠說明書獲取； (3)利用拉格朗日乘子法將步驟(2)中的約束條件和步驟（1)中的目標函數(shù)構成一個拉 格朗日函數(shù)，拉格朗日函數(shù)的構造方法已在數(shù)學界的最優(yōu)化理論中有了詳細說明，利用牛 頓-拉夫遜法，求解該拉格朗日函數(shù)，得到電-氣耦合多能流網(wǎng)絡狀態(tài)估計的結果。
【文檔編號】H02J3/06GK105958531SQ201610364751
【公開日】2016年9月21日
【申請日】2016年5月27日
【發(fā)明人】孫宏斌, 郭慶來, 王彬, 董今妮, 潘昭光
【申請人】清華大學