齊四清,任春雷,張 悅,李鵬飛
(國網內蒙古東部電力有限公司信息通信分公司,內蒙古 呼和浩特 010010)
隨著能源互聯網、泛在電力物聯網建設進程的不斷加快,電力系統的規模變得越來越大,結構也變得越來越復雜,監管地市、省級電網甚至全國互聯大電網系統,已經成為我國電力系統發展的必然趨勢[1-2]。電力行業關系到整個國民經濟的發展,電力系統的安全穩定運行也與國民生活密切相關。隨著電網建設密集度加大,電網運行以及管控方式也逐漸變得更加復雜。然而,電網的安全以及穩定性影響因素也逐漸增加,這就需要提高技術措施來保證電力系統能夠安全經濟穩定地運行。
在現代電網中,各區域相互聯系、緊密相關,尤其當太陽能、風能、地熱能、自備發電站等一系列分布式微電網系統接入電力系統之后,一個局部的小擾動或者異常運行都可能引發整個電力系統發生連鎖反應,嚴重的話還會導致大面積停電、電力系統崩潰解列等惡性事故的發生。近年來,國際國內就頻繁發生了多起電網大規模連鎖故障,對整個國民經濟、社會生活都造成了嚴重的破壞[3]。
例如,在2003 年8 月14 日起因于俄亥俄州一條345 kV 線路跳閘發生的北美五大湖停電事故,紐約州80%供電中斷[4];2005 年8 月18 日起因于一條線路故障所導致的關聯主變以及過載設備跳閘,造成印尼雅加達城市出現“黑城”現場[5];在2005年9 月26 日,海南受到“達維”臺風的影響,海南電網受到嚴重的破壞,導致了海南部分市區出現大面積停電[6];除此之外,還有2018 年巴西的“3.21”大范圍停電,2019 年“3.7”委內瑞拉大規模停電,2019年“7.13”美國紐約曼哈頓大面積停電等事故。
截至目前,大停電事故仍然是電力系統研究的熱點,而過往研究表明,在大面積停電的初始階段往往是少量元件的相繼故障,但在事故擴大階段則與系統中的薄弱環節有著密切的聯系[7]。在故障發生后,電力系統中的關鍵輸電線路起到了推波助瀾的作用,加快了電力系統的崩潰速度。因此,尋找關鍵輸電線路并研究其運行規律,對預防連鎖故障蔓延、大規模電網調度控制以及提高電力系統穩定性方面有著重大意義。
目前國內外關于關鍵輸電線路辨識主要包含兩個方向,一是基于潮流計算及穩定分析的研究,結合復雜系統自組織臨界性和模式搜索策略,如OPA 模型[8]、CASCADE 模型[9];二是基于復雜電網拓撲理論的研究,從電網拓撲結構及參數角度分析故障連鎖的根本原因,如小世界網絡模型[10]。早先時候利用加權介數指標和介數值[11]來辨識關鍵輸電線路,但前提是假設功率是在最短傳輸路徑上流動,不符合電網的實際運行情況,無法體現電網中功率流動特點。后來為了彌補以往介數指標潮流只是沿最短電氣距離傳播的不足,定義最大流傳輸貢獻度[12]為電力系統關鍵輸電線路和節點辨識依據,解決了目前研究的有功流傳輸距離上面的研究缺陷,考慮了線路最大有功傳輸容量的約束,提出了有向加權網絡,但未考慮關鍵輸電線路受到擾動之后對其余線路潮流分布影響的不足。為了更加符合電力系統運行的特性,提出了基于綜合介數指標與電網傳輸效率[13]的關鍵線路識別方法,綜合考慮了線路功率傳輸裕度以及故障之后潮流轉移等電網的實際物理特性,科學量化了關鍵輸電線路在電網的結構與運行之中的作用。但目前關于關鍵輸電線路辨識方法的研究都沒有考慮關鍵輸電線路的運行規律。
隨著電網結構越來越復雜以及建設規模越來越大,電網分析的難度也變得越來越大。為了將復雜問題簡單化,眾多的專家學者進行研究,如文獻[14]研究了一類重要的不可微非線性電阻網絡的直流分析問題,包括分段線性和連續可微兩種情況。引入一個局部函數來推廣單調性的概念,并給出了它的一些有用性質,利用局部函數研究了這類電阻網絡中直流工作點的存在唯一性。文獻[15]給出了非線性電阻互易或非互易電路中電壓或電流源到電阻電壓或電流的傳輸特性的單調性判據。文獻[16]基于解耦的有功相位角潮流方程,計算給定相位角安全限值下的發電裕度,結合等效電路和單調關系直接識別給定負載或發電機功率變化的臨界傳輸線,但這些文獻對線路單調性的研究還是不足。文獻[17]結合電網載流不均性和網絡響應單調性,提出一定運行模式下,關鍵輸電線路具有主宰特性,利用關鍵元件的單調性來簡化電力網絡的分析問題,卻沒有研究不同運行模式下的關鍵輸電線路的單調性。
運行模式是研究電力系統規律的關鍵之一,它是具有標樣特征指標向量的運行電力系統的類別形式[18],由某時刻發電、負荷以及網絡結構總體構成,具有綜合性,反映了系統節點、局部和整體的不同層次的關系。建立運行模式的目的是提取出典型的電力系統運行狀態,從而分析關鍵輸電線路具有的單調性規律。定義P為運行模式向量,即發電、負荷節點有功注入功率,因運行模式與發電、負荷以及結構相關,因此可用式(1)表示:

式中:PG表示發電機節點注入有功功率向量;PD表示負荷節點注入有功功率的向量。
設電力系統中發電節點、負荷節點之間的比例系數為ω,一定運行模式下發電機、負荷之間傳輸的有功功率大小為P′,則式(1)可以表示為:

實際上,電網中的每一個負荷節點都具有規律性,即通過每天的負荷曲線就可以簡單觀測到負荷的基本規律。當然通過負荷模型,比如AR 預測法,可以表示負荷的時序規律性。
假設電網結構保持不變,在每隔Δt=5 min 較短的時間內,每個負荷節點僅僅發生小擾動,那么在t+Δt時刻負荷節點的有功功率可以表示為:

由式(5)同樣可以看出,在較短的時間范圍內,發電機節點j在t+Δt時刻的有功變化量可以近似線性化表示為:

那么,根據式(3)~式(6),在t+Δt時刻,電力系統的運行模式可以表示為式(1)和式(2)的形式:

由式(7)可見,當在t時刻,基態(即發電與負荷之間輸送有功功率為Pt,或者式(2)中的P′)下,辨識出關鍵輸電線路后,那么,在t+Δt時刻負荷僅僅發生小擾動情況下,只要辨識式(7)中的系數就可以進一步確定關鍵輸電線路,以及其安全性。
電力系統關鍵輸電線路就是對于電網輸電其至關重要的輸電通路,該輸電通路承擔比例較高的有功輸電。按此認識,可以給出關鍵輸電線路的辨識模型:

式中:Pm表示發電機節點m的有功出力;Pn表示負荷節點n的有功負荷;{min(Pm,Pn)}表示求取二者的最小值;Pij(m,n)中的ij表示輸電線路首節點為i,末節點為j,該式子含義是發電機節點m輸送有功功率經過線路ij,且輸送了Pij(m,n)有功功率至負荷節點n;P(m,n)表示發電機節點m向負荷節點n輸送的總有功功率。
由式(10)可以看出,發電機節點m向負荷節點n輸送有功功率過程中,要經過多條輸電線路。同理,設G個發電機節點,D個負荷節點,那么G個發電機節點向D個負荷節點輸送有功功率過程中,若按照潮流追蹤原理,可以發現要經過輸電線路ij多次,由此定義輸電線路ij的關鍵性辨識指標為:

通過式(11)可以看出,輸電線路ij被利用的次數越多、功率越大,顯示了其關鍵性,這也是實際電網中考量輸電線路關鍵性的核心要素。
要計算式(11),核心是計算式(10),即計算式(10)中的Pij(m,n)和P(m,n),按照傳統潮流追蹤算法可得P(m,n)計算模型為:

式中:PDn表示負荷節點n的有功負荷;Pn表示節點n的有功注入功率;PGm表示發電機節點m有功出力;Aunm表示潮流追蹤中的逆序矩陣中的nm元素,可以表示為:

Pij(m,n)的計算模型為:

式中:Adin表示潮流追蹤過程中的順序分配矩陣,其元素計算方式為:

將式(12)和式(14)代入式(11)可得:

按照第1.2 節的方式可以辨識出任意時刻,例如t時刻電力系統的關鍵輸電線路ij。然而,負荷在t+Δt時刻將發生小擾動變化,則發電機有功出力也將隨之變化,從而關鍵輸電線路ij上的潮流也將變化,若在t+Δt時刻,輸電線路ij仍然是關鍵輸電線路,則說明關鍵輸電線路ij在[t,t+Δt]時間段內,具有單調性。
當然,該單調性的核心是與式(7)中的負荷基值及其變化量相關,即:

通過數學知識可知,單調性必定是指某一區間,由此設關鍵輸電線路ij的單調性表示為MPij,則其定義域為:

隨著時間的推移,按照式(7)的方式可以進一步獲得定義域的表達式。
當然,若在某時刻t+kΔt(k為某一正整數),按照式(7)獲得關鍵定義域后,但是經過第1.2 節部分辨識,發現輸電線路ij不是關鍵輸電線路,此時的關鍵輸電線路是cl,則記錄此時的關鍵輸電線路定義域:

電力系統中包含若干個負荷節點,每個負荷節點都有其對應的日負荷曲線,負荷節點k的某月第1 周7 天的日負荷曲線如圖1 所示。

圖1 某周7 天的日負荷曲線
由圖1 可見,某周內的7 天負荷曲線波形具有相似性,即走向、趨勢都一致,即在一天的24 h 內,對應同一個橫坐標,7 天負荷的縱坐標幅值稍微有所變化。
由此可知,在電力系統配電網全時空量測環境下,通過每天記錄電力系統中每個負荷節點日負荷曲線,按照1.2 節方法辨識關鍵輸電線路,記錄關鍵輸電線路及其定義域,就可以以年為單位積累若干年的關鍵輸電線路及其單調性規律。
但是負荷的變化有其特點,如圖1 中在0:00—5:00 時間范圍內,負荷變化較小,此時屬于一種運行模式;在6:00—11:00,此時負荷變化較大,應該屬于另一種運行模式,以此類推。
經過全時空量測,以及按照1.2 節關鍵輸電線路辨識方法,設以天為單位,記錄配電網Td個運行模式向量:

式中:每個元素表示一種運行模式,該運行模式按照式(18)或者式(19)形式具有對應的關鍵輸電線路定義域。也就是說,式(20)中任何一種運行模式,對應一種關鍵輸電線路及其定義域。
經過若干年的積累,可以形成式(20)的大數據樣本,表示為:

獲取單調性大數據樣本進而預測的目的就是為了針對未來負荷的異常變化而導致的運行模式變化,進而導致關鍵輸電線路變化。
深度學習是目前機器學習領域的一個高級別研究方向,它通過學習數據樣本的內在規律和特征,進而能夠預測得未來規律。雙向長短時記憶網絡(Bi-directional Long Short-Term Memory,BLSTM)是在LSTM 基礎上發展而來的,它能夠根據輸入時序樣本數據,采用雙向準則判斷,挖掘多種輸入序列之間的潛在關系,進而獲得所需的預測規律。
傳統的LSTM 結構如圖2 所示,其由包含輸入門xt、遺忘門ft、輸出門Ct等構成。

圖2 LSTM 結構圖
其工作原理為:
首先遺忘門的輸出為當前輸入xt、前一時刻隱含信息ht-1作為輸入,經過作用函數σ的計算后得到:

式中:wfx、wfh為需要訓練的權值;bf為擾動量。
由圖2 可見,it對輸入的數據進行加工篩選,從而減少計算維度:

式中:wix、wih也為需要訓練的權值;bi為隨機量。
經過式(21)和式(22)的計算,可以得到更新的信息:

式中:wcx和wch為需要訓練的權值;bc為隨機量。
經過更新以及LSTM 原有的積累,可以得到t時刻LSTM 計算值:

根據圖2,進而可以得到當前輸出值:

式中:wy為需要訓練的權值;by為隨機量。
由上述可見,LSTM 能夠解決歷史數據學習功能,但是沒有充分利用未來的信息,比如未來負荷預測,因此采用BLSTM 能夠解決該問題。
BLSTM 是在LSTM 結構的基礎上,添加了前、后兩個時間序列相反的隱含層;然后將其連接到一個輸出,前向的隱含層為式(25),后向的隱含層為:

根據BLSTM 能夠將未來負荷預測的信息一起進行學習,從而有助于預測未來關鍵輸電的單調性規律。由此,設負荷未來預測列向量為,那么BLSTM 的輸入向量為:

式中:MP表示式(20)中的運行模式向量;表示運行模式向量相關的定義域向量;表示與MP運行模式相對應的歷史負荷值向量。
將式(28)作為輸入,進行式(21)~式(27)計算,可預測得到未來關鍵輸電線路單調性規律。
為了驗證上述算法的有效性,以如圖3 所示的IEEE5 節點系統進行算法驗證,該系統包含5 個節點,3 條輸電線路。

圖3 5 節點系統仿真圖
算法的驗證過程中,在Windows7 系統中使用MATLAB 軟件進行程序編寫和分析。
圖3 所示為IEEE5 節點系統結構圖,節點1 為PV 節點,節點2 為平衡節點,節點3~5 為PQ 節點。具體參數見表1。

表1 IEEE5 節點系統參數(標幺值)
在初始狀態時,發電機向負荷輸送的有功功率為:

保持初始狀態不變,即運行模式不變,改變發電機向負荷傳輸的有功功率ω,關鍵輸電線路,如表2所示。

表2 輸電線路有功潮流(標幺值)
由表中可知,在運行模式不變的情況下,隨著發電向負荷傳輸有功功率,輸電線路L3-4 始終是關鍵輸電線路,直至達到輸送的最大有功功率為10.38,由此可見,在[10,10.38]區間里輸電線路L3-4始終是關鍵輸電線路,即其具有單調性。
進一步采用內蒙古電網某季度的典型日負荷曲線探討在不同運行模式下,預測關鍵輸電線路具有單調性。
按照負荷曲線變化趨勢,把一天24 h 的運行模式以5 min 為間隔進行劃分,即可得到288 個點負荷,一般來說,電力系統5 min 內的運行模式保持不變,故而可給定電力系統288 個不同的運行模式,3個負荷共288 種組合方式,某實際電網某季度日負荷曲線如圖4 所示。

圖4 典型日負荷曲線
接著利用PowerWorld 進行分時段潮流仿真,插入288 個時間點,每個時間段間隔5 min,將每個時間點負荷數據輸入時段有功負荷進行仿真,查看實例信息中的電力系統潮流數據。
以某實際電網某季度共90 d 負荷數據利用MATLAB 按照第1.2 節進行關鍵輸電線路辨識,獲得大數據樣本。以某季度前89 d 的大數據樣本,按照式(21)~式(28)進行學習訓練,使用最后一天數據對學習預測的規律進行算法驗證。
隨著不同時間段各個負荷的變化,各個節點發電機有功出力也隨之變化,即運行模式不斷發生變化,從而導致關鍵輸電線路在電網中也跟著發生變化,且在研究的時段中,可以發現線路L3-4、L3-5在一定程度上制約了其他線路的潮流,這2 條線路隨著運行模式的變化,關鍵性輪流變化。
通過預測發現:在0:00—6:00 時間區間內,輸電線路L3-5 始終是關鍵輸電線路;而在6:00—23:55時間內,關鍵輸電線路變成了L3-4,這說明輸電線路L3-4 和L3-5 是從發電到負荷的重要輸電走廊,這與實際情況一致。
本文從考慮配電網輸電安全的角度出發,建立了配電網關鍵輸電線路單調性規律提取方法。通過在全時空配電網潮流量測的基礎上,建立了關鍵輸電線路的辨識模型,并以此為大數據樣本,結合負荷的超前預測,使用BLSTM 學習該樣本,進而獲得未來關鍵輸電線路及其單調性。通過實例驗證,結論如下:
(1)一種運行模式會對應著一條關鍵輸電線路,由于每條線路傳輸的有功功率限值不變,導致關鍵線路會發生變化,而關鍵輸電線路又與負荷變化相關,故而探尋關鍵輸電線路規律的關鍵主要就是負荷和發電的變化情況。
(2)關鍵輸電線路隨著負荷發生單調變化,在運行模式發生改變時,關鍵輸電線路在一定時間區間內保持不變,因負荷變化導致的潮流轉移會使得關鍵輸電線路發生改變,故而關鍵輸電線路具有局部單調性。