任建文， 何培成

(華北電力大學 新能源電力系統(tǒng)國家重點實驗室，河北 保定 071003)

基于多支路切除和背離路徑的潮流轉移識別算法

任建文， 何培成

(華北電力大學 新能源電力系統(tǒng)國家重點實驗室，河北 保定 071003)

給出了一種基于多支路切除和路徑搜索的潮流轉移快速識別算法。該方法根據(jù)故障支路切除后的潮流轉移特征，首先深入分析了轉移到非故障線路上的潮流組成，其次在一定的節(jié)點深度內搜索切除支路兩端節(jié)點間的前K條最短路徑集，最終識別出受潮流轉移影響較大的線路構成輸電斷面。文中在切除支路的潮流轉移估算、K值選取和斷面判據(jù)方面作了重要改進，以適應于發(fā)生多支路切除情況下潮流轉移的識別；重新建立了一種K的取值方案，并采用背離路徑算法確定最短路徑集，在減少分析計算量的同時有效避免了輸電斷面漏選問題。IEEE30節(jié)點系統(tǒng)的仿真驗證了本文算法的有效性。

電力系統(tǒng)；潮流轉移；輸電斷面；路徑搜索；背離路徑

0 引 言

電網中由故障支路切除引起的潮流轉移是近年來世界各國發(fā)生輸電線路連鎖切除事故的重要原因[1-3]，嚴重危害著系統(tǒng)的安全穩(wěn)定運行。在此背景下，目前電網調度中心嚴格按照《電力系統(tǒng)安全穩(wěn)定導則》要求安排電網運行方式，運行中嚴格執(zhí)行“N-1”準則[4]，如設備檢修，機組調控等操作，在系統(tǒng)運行狀態(tài)發(fā)生變化時及時消除安全隱患，這樣可以使由單一元件開斷引發(fā)支路連鎖跳閘的幾率越來越??；但是，國內大部分電網還沒有達到能夠抵御同一斷面多條支路切除嚴重故障的強度[5]，這主要是因為一方面校驗計算復雜，另一方面若要滿足多支路切除安全校驗需提高線路的熱穩(wěn)定極限功率。由于電網缺乏對多支路切除情況的應對措施，一旦發(fā)生此類故障時任由后備保護切除線路，大大增加了輸電線路連鎖跳閘的幾率。所以，對多條支路切除后的輸電斷面的行為分析，即快速找出潮流轉移較大的支路并判斷其上潮流是否越限，對預防線路連鎖跳閘具有重要意義。

文獻[6]提出潮流轉移因子的概念，將多支路切除潮流轉移等效為切除支路間潮流的來回折返過程，推導得到潮流轉移的計算式，有效對轉移到非故障支路的潮流進行估算。文獻[7]給出了一種基于直流公式并結合故障后網絡拓撲參數(shù)來求解潮流轉移的方法，然而此方法需要修改節(jié)點阻抗矩陣，雖然采用矩陣求逆定理簡化求解過程，但計算量偏大，不滿足快速性要求。文獻[8]提出有功轉移系數(shù)的概念，并結合廣度優(yōu)先搜索算法對切除支路的潮流轉移區(qū)域進行界定，將區(qū)域中有功潮流轉移系數(shù)大于某一特定值的線路加入輸電斷面，但搜索限值由人為經驗而定，其適應性較差。文獻[9]指出受潮流轉移影響嚴重的支路均包含在切除支路兩端節(jié)點間的若干條最短路徑中，并采用路徑回溯法搜索輸電斷面，但多支路切除潮流轉移的規(guī)律與單一支路切除情況差異較大，無法直接套用。

本文在文獻[6]基礎上，對轉移到非故障支路上的潮流組成進行分析；依據(jù)圖搜索中背離路徑算法，在一定的節(jié)點深度內搜索切除支路兩端節(jié)點間的前K條最短路徑，確定潮流轉移較大的區(qū)域范圍。文中結合了多支路切除的潮流轉移特征，重新建立一種K的取值方案，并采取邊搜索邊修正K值的策略，避免搜索范圍設置不當而漏選支路的問題；考慮電網的實時運行狀態(tài)，給出了輸電斷面的選取原則，最終識別出實際潮流增加較大的支路集，且能成功找出潮流反向增大的線路。

1 多支路切除潮流轉移分析

1.1 支路開斷分布因子

若電網中發(fā)生支路切除，由切除支路l轉移到支路k上的潮流量(以初始潮流方向為正)可通過l對k的支路開斷分布因子Dk-l表示：

(1)

其中：Pl為切除支路l的初始有功潮流；ΔPk為轉移到支路k上的潮流。

由式(1)可知，若已知分布因子即可計算切除支路對任意支路的潮流轉移量，而Dk-l可由系統(tǒng)的節(jié)點阻抗矩陣計算得到，其計算式如下：

(2)

其中：Xk-l、Xl-l、xk、xl只與電網的拓撲參數(shù)有關[10]。對于一個確定的網絡，支路l對任意支路k的分布因子Dk-l隨之確定。

1.2 多支路切除潮流轉移模型

假設初始潮流分別Pa與Pb的支路a、b發(fā)生故障切除，支路k為網絡中除支路a、b外的任意一條支路。根據(jù)文獻[6]思想：多支路切除潮流轉移可以分解為切除支路間無窮次潮流的來回折返過程。首先，支路a、b的初始潮流分別依據(jù)分布因子Dk-a、Dk-b轉移到支路k上，此外Pa按照分布因子Db-a轉移到支路b，Pb按照分布因子Da-b轉移到支路a，完成一次支路間潮流交換過程；其次，上述支路間交換的潮流替代了切除支路的初始潮流，繼續(xù)對網絡進行轉移，依此類推，直至潮流全部轉移到網絡中其他支路上，如圖1所示。

圖1 多支路切除的潮流轉移模型Fig.1 Power flow transferring modle of multiple line tripping

1.3 轉移到支路k上的潮流成分分析

(3)

(1)Pb=0時

此時有如下關系式：

(4)

(5)

由式(4)與式(5)推得：

(6)

(2)Pa=0時

此時有如下關系式：

(7)

(8)

由式(7)與式(8)推得：

(9)

(10)

(11)

將式(10)與(11)代入式(3)可得：

(12)

由式(12)可得：發(fā)生多支路切除時，支路k上的轉移潮流可看成多條以等效潮流為初始潮流的支路單獨開斷后轉移到支路k上潮流的疊加，將多支路切除簡化為多條單一支切除情況，避免了直流法中對節(jié)點阻抗矩陣的大量修正計算，并實現(xiàn)了對潮流轉移的快速估算。

2 潮流轉移的路徑搜索方法

2.1 前K最短路徑法

由電路基本方程可知：當發(fā)生支路切除時，可將切除支路用附加電流源(與支路切除前的電流等量、反向)替代后的等值無源網絡中的潮流分布來近似切除支路的潮流轉移。若對附加電源搜索其前K條最短閉合路徑回路，則這些回路包含該等值網絡中絕大多數(shù)電流較大的輸電線路，也即包含了受潮流轉移影響較大的支路，理論上說明了前K最短路徑法搜索輸電斷面的可行性。

2.2 改進的K值選取方案

(13)

本文采用的取值方案：首先由ΔPk1確定支路a的基本路徑數(shù)Ka，其次把ΔPk2作為輔助判據(jù)，修正Ka取值，將經修正后的Ka作為支路a的所求路徑數(shù)。其步驟如下：

步驟1：確定基本路徑數(shù)Ka

ΔPk1是支路a直接對支路k進行轉移的潮流，與支路a發(fā)生單支路切除的轉移潮流分布規(guī)律是一致的，可類比于單支路切除時的K值求取方法：

(14)

其中：β為潮流轉移系數(shù)的閾值，考慮算法的效率，一般取0.3[9]。

(15)

結合式(14)與式(15)可得：

(16)

步驟2：修正Ka

由最短路徑可知，除去前Ka條最短路徑外，轉移到支路a兩端節(jié)點間的第Ka+1條最短路徑上的潮流是其余最短路徑中的最大值。因此，只需判斷第Ka+1條最短路徑中加入分量ΔPk2后，其上轉移潮流是否越限即可，修正過程如下所述：

①搜索第Ka+1條最短路徑。

②找出該路徑中與前Ka條路徑不關聯(lián)的線路構成新增線路集，計算新增線路集中任意線路t對支路a的潮流轉移系數(shù)ΔPt/Pa，若滿足下式，則令Ka=Ka+1。

(17)

式中：ΔPt為轉移到支路t上的潮流量；ε為潮流轉移系數(shù)的上限系數(shù)。當ΔPt/Pa大于ε時，說明ΔPk2占ΔPt的比例較大，即支路t有更大的幾率包含在支路b兩端節(jié)點間的前若干條最短路徑中，為了避免重復搜索，故不修正Ka。經過大量仿真驗證，當ε=0.5時，能夠有效抑制Ka的增大且搜索結果包含所有受潮流轉移影響較大的支路。

③重復執(zhí)行①、②直至第Ka+1條最短路徑上沒有滿足式(17)的新增線路。

同理可得支路b的路徑數(shù)Kb的取值。

綜合上述改進方案：考慮了多支路切除潮流轉移特征，并合理的在單支路切除時K值選取原則的基礎上進行改進，易于實現(xiàn)；所有參數(shù)均能通過初始條件求得，說明本方案可行；一般A值不會太大，滿足快速性要求。

3 快速識別流程

目前已積累許多特點各異的前K最短路徑算法，其中背離路徑算法通過求取第K-1條路徑的背離路徑來確定其第K條最短路徑[11]，與本文K值的修正過程相符，且每派生一條背離路徑只需執(zhí)行一次Dijstra算法，算法的時間復雜度較小，更易滿足搜索的實時性要求。綜合上述優(yōu)點，本文采用背離路徑算法對發(fā)生多支路切除時受潮流轉移影響嚴重的線路進行識別，背離路徑的相關概念詳見文獻[11]，本文在此不作過多闡述。識別算法流程按照下述模塊化進行設計。

(1)潮流轉移分析模塊

該模塊的功能是提供網絡拓撲結構與切除支路的等效潮流，為路徑搜索提供初始數(shù)據(jù)，其流程如圖2所示。

圖2 潮流轉移分析流程圖Fig.2 Flow chart of analysis on power flow transferring

(2)路徑搜索模塊

該模塊的功能是搜索切除支路兩端節(jié)點間的若干條最短路徑，對受多支路切除潮流轉移影響較大的區(qū)域進行界定。由于電網中受潮流轉移影響較大的元件不超過全網的20%[12]，所以沒必要進行全網分析，只需在一定的節(jié)點深度T內進行搜索即可，進而提高算法效率，搜索流程如圖3所示。

圖3 路徑搜索流程圖Fig.3 Flow chart of path search

(3)潮流轉移識別模塊

該模塊功能是結合電網的實時運行狀態(tài)，識別出路徑集中潮流實際增加較大的線路；將潮流增量大于某一特定值的線路加入輸電斷面。

(18)

其中：ΔPk.0為支路k的潮流增量；Pi，min為最小的切除支路潮流；γ為有功增加系數(shù)閾值，本文取0.3。

然而，電網中潮流實際增加較大輸電線路存在兩種情況：①其上轉移潮流與初始潮流方向一致；②其上轉移潮流與初始潮流方向相反。對于不同情況，線路上的潮流增量是不同的。

對于情況1：線路上的潮流增量即為轉移到其上的潮流量，即ΔPk.0=ΔPk。代入式(18)可得輸電斷面的判據(jù)：

(19)

對于情況2：轉移到線路的潮流ΔPk<0，僅僅反應潮流減小，為了反應潮流增量，文獻[13]給出了該情況下潮流增加的判據(jù)：

(20)

式中：Pk為支路k的初始潮流。

(21)

因此，該模塊的流程設計如圖4所示。

圖4 潮流轉移識別流程圖Fig.4 Flow chart of recognition of power flow transferring

4 仿真算例

對IEEE30節(jié)點系統(tǒng)進行仿真驗證，將不能與線路構成閉合回路的發(fā)電機與負荷節(jié)點從網絡中剔除，得到的系統(tǒng)簡化拓撲如圖5所示。K等于4是相對保守的估計，一般第K條最短路徑中的新增支路的潮流轉移系數(shù)達不到1/K[14]。因此，基準值K取3不會影響潮流轉移的識別結果。

圖5 IEEE30節(jié)點系統(tǒng)簡化拓撲圖Fig.5 The simplified topological graph of IEEE30 bus-system

表1中基本路徑數(shù)為5條，第6條路徑中新出現(xiàn)支路L6-10的潮流轉移系數(shù)為0.48，滿足式(17)，說明L16-17的等效潮流轉移到L6-10的分量使得其潮流轉移系數(shù)超過閾值，則對Ka進行修正；繼續(xù)搜索第7條路徑，然而第7條路徑中新增支路L4-2與L2-6的潮流轉移系數(shù)分別為0.13與0.14，不滿足式(18)，不需要修正Ka。同理表2中Kb的修正過程與表1類似，最終經修正后Ka=6、Kb=5。

將路徑集中的線路依據(jù)斷面選取原則進行識別，其結果如表3所示：

表1 對L18-19的路徑搜索數(shù)據(jù)

表2 對L16-17的路徑搜索數(shù)據(jù)

表3 L18-19和L16-17切除事故的識別結果

表4 對L28-27的路徑搜索數(shù)據(jù)

表5 對L12-15的路徑搜索數(shù)據(jù)

表6 L28-27和L12-15切除事故的識別結果

由結果可知，算例1、算例2中分別搜索到10條、15條潮流增量較大的線路，相較于精確潮流計算所得的有功增加量，說明本文對轉移到非故障支路上的潮流計算方法能有效地近似線路的實際有功增加值，為進一步判斷支路切除后輸電斷面內線路是否過載提供參考依據(jù)。并且有效找出潮流反向增大的支路，如表6中線路L21-22與L25-27。

對比分析算例1和算例2：從等效潮流與初始潮流的比值A可以看出，算例1中潮流間的耦合強度相較于算例2而言更為嚴重，導致算例1經修正后的路徑數(shù)大于算例2的數(shù)值。在極端情況下，如開斷線路為L15-23與L27-29時，切除線路間的開斷分布因子為0，比值A=1，則開斷潮流間沒有耦合關系，此時轉移潮流在網絡中的分布和單支路切除時是一致的。因此，系數(shù)A可以反映該支路與其他斷線支路間的潮流耦合強度，并能夠進一步預判此切除事故的嚴重程度。A值越接近1，切除支路的開斷潮流在網絡中的分布越接近單支路切除情況，隨著A值增大，其開斷潮流的分布與單支路切除情況時的差異也越大，事故也越嚴重。

本文與文獻[8]中所述的潮流轉移搜索方法對比可得：(1)從理論基礎上，文獻[8]指出發(fā)生多支路切除后，有功潮流沿著其他沒有斷開的線路由切除線路的始端流向各斷線支路的末端，則需要搜索所有斷線支路始端節(jié)點與末端節(jié)點間的若干最短路徑，增大了分析計算量，而本文經過潮流分析計算，將多支路切除轉化為多條相互獨立的單支路切除故障，只需搜索斷線支路兩端節(jié)點間若干最短路徑即可；(2)從搜索方法上，文獻[8]未對流向其他斷線支路的潮流進行定量分析，從而采用經驗系數(shù)作為搜索的終止判據(jù)，而本文給出了合理的K值修正方案，邊搜索邊修正終止判據(jù)，從而在減少冗余計算的同時，提高搜索效率。

5 結 論

本文將路徑搜索引入多支路切除的潮流轉移識別算法中，是一種新的嘗試。文中通過分析計算，在獲取潮流估算式的同時得到經修正后的路徑數(shù)，進而結合路徑搜索與實時運行狀態(tài)識別出潮流增加較大的支路，所提算法具有如下特點：

(1)潮流轉移分析模塊采用了極限的思想，保持了初始網絡的拓撲結構，避免對節(jié)點阻抗矩陣進行修改，大大減少矩陣運算量，提高算法效率。

(2)考慮了輸電線路的初始潮流，從而能夠識別得到實時的潮流增加較大支路集，所得結果更加符合實際電網運行情況。

(3)采用邊搜索邊修正K值的方法，避免搜索范圍設置不當而漏選支路的問題。

(4)給出了反映潮流耦合強度的比例系數(shù)，通過計算該系數(shù)的取值可粗略判斷多支路切除事故的嚴重程度。

[1] 毛安家，張戈力，呂躍春，等．2011年9月8日美墨大停電事故的分析及其對我國電力調度運行管理的啟示[J]．電網技術，2012，36(4)：74-78．[2] 楊文輝，畢天姝，馬強，等．基于廣域電壓相角信息的輸電斷面快速識別方法[J]．電力系統(tǒng)保護與控制，2013，41(24)：58-63．

[3] 方勇杰．用緊急控制降低由輸電斷面開斷引發(fā)系統(tǒng)崩潰的風險對印度大停電事故的思考[J]．電力系統(tǒng)自動化，2013，37(4)：1-6．

[4] 梁志峰，葛睿，董昱，等．印度“7.30”、“7.31”大停電事故分析及對我國電網調度運行工作的啟示[J]．電網技術，2013，37(7)：1841-1848．

[5] 高翔，莊侃沁，孫勇．西歐電網“11.4”大停電事故的啟示[J]．電網技術，2007，31(1)：25-31．

[6] 徐慧明，畢天姝，黃少鋒，等．計及暫態(tài)過程的多支路切除潮流轉移識別算法研究[J]．中國電機工程學報，2007，27(16)：24-30．

[7] 袁曉丹，張會強．多支路開斷潮流轉移識別及防連鎖過載策略研究[J]．現(xiàn)代電力，2014，31(5)：74-79．

[8] 劉獻超，任建文，渠衛(wèi)東．基于廣度優(yōu)先法的多支路連鎖切除潮流轉移快速搜索算法[J]．電力系統(tǒng)自動化，2015，39(13)：69-74．

[9] 王增平，李剛，任建文．基于前K最短路徑的輸電斷面搜索新算法[J]．電工技術學報，2012，27(4)：193-201．

[10] 張伯明，陳壽孫，嚴正．高等電力網絡分析[M]．北京：清華大學出版社，2007．

[11] 傅俊偉，李興明，陳捷．基于背離路徑的Kth最短路徑實用搜索算法[J]．計算機技術與發(fā)展，2009，19(2)：120-126．

[12] 周德才，張保會，姚峰，等．基于圖論的輸電斷面快速搜索[J]．中國電機工程學報，2006，26(12)：32-38．

[13] 李莎，任建文．基于有功增加因子的潮流轉移快速搜索[J]．電網技術，2012，36(12)：176-181．

[14] 苗世洪，馬帥，尚亞男，等．基于割點和路徑搜索的輸電斷面快速識別方法[J]．電力系統(tǒng)自動化，2014，38(2)：39-45．

Recognition Algorithm of Power Flow Transferring Based on Multi-branch Removal and Deviation Path

REN Jianwen, HE Peicheng
(State Key Laboratory of Alternate Electrical Power System with Renewable Energy Sources, North China Electric Power University, Baoding 071003, China)

A recognition algorithm of power flow transferring based on multi-branch removal and path research is proposed in this paper. Firstly, an in-depth analysis on the component of power flow transferring on the non-faulted lines is presented on the basis of power flow transferring’s characteristics. Secondly, the former shortest paths represented by the letterKbetween two nodes of removed lines would be found within a certain node depth. Finally, the transmission sections whose lines are greatly influenced by the power flow transferring within the above paths are recognized. Considerable improvements of estimation of power flow transferring, optimal selection of the value ofKand criterion for transmission section are put forward which are suitable for recognition of power flow transferring with multi-branch removal. Moreover, a method to ascertain the value ofKis reinstituted and the shortest paths are determined based on deviation path which reduces calculated quantities and effectively avoids miss-selection of transmission section. The simulation results of IEEE30 node system verified the effectiveness of the proposed algorithm.

power system; power flow transferring; transmission section; path search; deviation path

2016-03-03.

國家自然科學基金資助項目(50837002).

10.3969/j.ISSN.1007-2691.2016.06.01

TM711

A

1007-2691(2016)06-0001-06

任建文(1961-)，男，教授，研究方向為電力系統(tǒng)保護與安全控制、電網調度自動化等；何培成(1991-)，男，碩士研究生，研究方向為電力系統(tǒng)保護與安全控制。