陳晶，劉明群，覃日升，李長松，胡清泉，李周龍，王家華

（1.云南電力試驗研究院＜集團＞有限公司，昆明650217；2.云南電力科學(xué)研究院，昆明650217；3.四川大學(xué)電氣信息學(xué)院，成都610065；4.云南臨淪供電局，云南 臨淪677000）

0 引 言

連鎖故障是近年來的研究熱點［1］。國內(nèi)外學(xué)者提出了很多連鎖故障分析理論和模型，從不同角度揭示了連鎖故障的演變。Dobson等學(xué)者利用復(fù)雜系統(tǒng)理論［2］，從宏觀演變趨勢上研究連鎖故障的變化，提出的電網(wǎng)自組織臨界態(tài)可以為調(diào)度人員提供一定的預(yù)警參考信息。文獻［3］認為保護系統(tǒng)在電網(wǎng)不正常狀態(tài)下的隱性故障是造成線路連鎖開斷的原因，利用隱性故障模型刻畫連鎖故障的蔓延。也有學(xué)者利用可信性理論［4］建立連鎖故障風(fēng)險評估框架，該方法可以一定程度上克服隨機性刻畫連鎖故障時存在的抽樣次數(shù)與精度的矛盾。

上述研究大多通過標準節(jié)點系統(tǒng)或?qū)嶋H系統(tǒng)離線仿真，未反映實時運行狀態(tài)，不能為調(diào)度人員提供在線預(yù)警操作的依據(jù)。而從離線分析到在線分析的關(guān)鍵是能從EMS系統(tǒng)中快速解析出電網(wǎng)模型數(shù)據(jù)。之前電網(wǎng)模型交換規(guī)范采用的CIM/XML模型具有冗余信息量大、效率較低的缺點［5］，不能很好滿足在線要求。為此，國家電力調(diào)度通信中心提出了E語言［6］，將CIM/XML模型優(yōu)化為CIM/E模型，能夠高效地描述模型數(shù)據(jù)，減小了文件大小，提升了文件解析和模型效驗的效率，能夠滿足較大規(guī)模電網(wǎng)對模型在線交換的要求［7-8］。然而結(jié)合CIM/E模型文件和連鎖故障現(xiàn)有理論方法的研究還比較少。

因此，本文著重研究了CIM/E文件和連鎖故障現(xiàn)有理論方法的融合，開發(fā)了基于CIM/E的連鎖故障在線辨識模塊。首先對CIM/E文件進行了正確性效驗，確保能夠讀取到電網(wǎng)模型數(shù)據(jù)。進而，基于E語言語法開發(fā)了讀取規(guī)則，將CIM/E模型數(shù)據(jù)解析到程序模塊中，便于后續(xù)潮流計算和拓撲分析。以上述為基礎(chǔ)，結(jié)合連鎖故障分析的可信性測度方法，識別出當(dāng)前運行狀態(tài)下風(fēng)險較大的連鎖故障序列，并將序列識別結(jié)果寫成E語言格式文件存于指定位置，為調(diào)度人員提供預(yù)警參考。本文開發(fā)的程序模塊已成功運用到西南某地區(qū)電網(wǎng)。

1 CIM/E模型

CIM/E文件中電網(wǎng)模型為物理連接模型，包含了區(qū)域、基準電壓、廠站、電壓等級、間隔、斷路器、刀閘、母線段、同步發(fā)電機、交流線段、負荷、變壓器、變壓器繞組、變壓器分接頭類型、并聯(lián)補償器、串聯(lián)補償器、遙測、遙信等22類對象，每一類對象的模型包含其標識、路徑、拓撲連接關(guān)系、狀態(tài)等信息。圖1即為斷路器類，“mRID”是原系統(tǒng)中的標識，在整個模型文件中唯一。“name”是原系統(tǒng)中的設(shè)備中文名稱。“I_node”、“J_node”、“Substation”是斷路器的拓撲連接信息，“status”是斷路器的量測狀態(tài)。

圖1 斷路器模型Fig.1 Model of breaker

2 CIM/E文件的讀取使用

2.1 CIM/E文件正確性效驗

為了確保CIM/E文件的正確性，提取數(shù)據(jù)前對CIM/E文件進行正確性效驗。檢查指定路徑位置是否存在對應(yīng)的 CIM/E文件，若找不到，則會提示“CIM/E文件不存在”的錯誤。找到之后，則再對該文件進行格式效驗，若格式錯誤，則會將錯誤的地方一次列出來，以便修改。格式正確將不會報出錯誤提示。

2.2 CIM/E文件數(shù)據(jù)提取與處理

完成上一步正確性效驗后，便將電網(wǎng)模型數(shù)據(jù)從CIM/E文件中提取到C/C++（本文開發(fā)的程序模塊基于C/C++語言）中，22類對象數(shù)據(jù)分別存入對應(yīng)名字的結(jié)構(gòu)體中。

為確保正確讀取，基于E語言語法［7］，開發(fā)了讀取規(guī)則：遇到“＜類名”開始讀取這一類數(shù)據(jù)，直到遇到“＜/”結(jié)束讀取，并開始讀取下一類數(shù)據(jù)；“@”、“//”字符分別為屬性名起始和注釋起始，遇到后直接跳過，直到遇到“#”開始的行，將該行中空格劃分的內(nèi)容分別存入到相應(yīng)結(jié)構(gòu)體中對應(yīng)名字的變量中。

此時存入結(jié)構(gòu)體中的數(shù)據(jù)還比較分散，還要通過數(shù)據(jù)處理生成便于潮流計算的BUS、LINE數(shù)組和便于拓撲分析的LSTA數(shù)組。BUS數(shù)組中包含節(jié)點信息，如發(fā)電機出力、負荷、節(jié)點電壓等，LINE數(shù)組包含支路信息，如線路電阻、線路電抗、電納、變壓器支路變比等。BUS數(shù)組中數(shù)據(jù)主要提取自BusbarSection、SynchronousMachine、Load等結(jié)構(gòu)體，LINE數(shù)組中數(shù)據(jù)主要提取自 ACLineSegment、ACLineDot、PowerTransformer等結(jié)構(gòu)體。LSTA數(shù)組存儲支路上斷路器的狀態(tài)，亦表征相關(guān)支路的運行狀態(tài)，提取自Breaker結(jié)構(gòu)體，并將在程序拓撲分析中起到重要作用。

2.3 拓撲分析

利用上述步驟從CIM/E文件獲取線路運行狀態(tài)OSTATE數(shù)組，進行電網(wǎng)孤島識別［9］。建立BISLAND與LISLAND兩個數(shù)組，分別存儲節(jié)點和線路的孤島號。兩個數(shù)組行數(shù)為1，列數(shù)分別表示節(jié)點數(shù)和線路數(shù)，均初始化為0，拓撲分析完成后，BISLAND數(shù)組中不存在零值，因為每一個節(jié)點都會被劃分到一個孤島中，LISLAND數(shù)組中可能存在零值，因為有線路斷開的情況。建立指標INUM、BUSN，INUM表示系統(tǒng)中當(dāng)前孤島號（初始化為1），BUSN表示下一搜索節(jié)點（初始化為1）。孤島識別的步驟如圖2所示。

（1）連續(xù)性搜索

從母線號1開始搜索，搜索母線上是否存在尚未被搜索過（LISLAND中對應(yīng)元素為0）并且在運（OSTATE中對應(yīng)元素為1，1為在運）的線路，若存在，則將該線路孤島號置為當(dāng)前孤島號（INUM），并檢查另一端母線是否被搜索過（BISLAND是否為0），若否，將其置為下一搜索母線（BUSN置為該母線號），此母線孤島號置為當(dāng)前孤島號（INUM）。

圖2 孤島辨識流程圖Fig.2 Flow chart of island identification

（2）廣度搜索

當(dāng)兩端母線均被搜索過時，連續(xù)性搜索失效，利用廣度搜索找到連續(xù)性搜索下一節(jié)點。從母線號1開始掃描完所有母線，檢查母線孤島號是否為當(dāng)前孤島號（BISLAND=INUM）并且其上有未搜索過的在運線路（LISLAND對應(yīng)元素為0并且OSTATE對應(yīng)元素為1），若滿足，置上述線路孤島號為當(dāng)前孤島號（INUM），并將最后一條線路上未被搜索到的母線號設(shè)為下一搜索節(jié)點（BUSN置為該母線號），返回步驟（1）連續(xù)性搜索。

（3）完整性搜索

當(dāng)連續(xù)性搜索和廣度搜索都找不到下一搜索節(jié)點時，從母線號1開始依次掃描，檢查是否存在未搜索過的母線（BISLAND為0）的母線，一旦發(fā)現(xiàn)存在，以該母線為新搜索母線（BUSN置為該母線號），當(dāng)前孤島號（INUM）加1，新搜索母線的孤島號置為改變后的孤島號，跳到步驟（1），繼續(xù)搜索，重復(fù)步驟（1）、（2）、（3），直至完整性搜索找不出新的母線。

3 連鎖故障序列辨識

從CIM/E文件中提取數(shù)據(jù)并處理完成之后，結(jié)合可信性測度方法，對連鎖故障序列進行辨識。

3.1 連鎖故障的可信性測度

文獻［4］將可信性測度［10］運用到電網(wǎng)連鎖故障風(fēng)險評估中。在精度偏差允許范圍內(nèi)，能更快速地獲取連鎖故障序列。

電網(wǎng)連鎖故障a發(fā)生的可信性測度Cr（a）可表示為：

其中：

式中ac為a的對立事件，可能性測度Pos（a）描述a發(fā)生的容易程度，必要性測度Nec（a）描述ac發(fā)生的不可能程度。

由可信性測度的四條公理［12］可知，連鎖故障a發(fā)生的可能性測度Pos（a）和必要性測度Nec（a）分別為：

式中aj為連鎖故障第j個階段，算子∧為取小運算。連鎖故障第j個階段aj發(fā)生的可能性測度和其補集發(fā)生的可能性測度通過聯(lián)合可信性分布函數(shù)定義求得，具體求解方法參加文獻［11］。

3.2 連鎖故障的嚴重性指標

嚴重性指標用于刻畫連鎖故障的后果。但連鎖故障的后果是廣泛的，電網(wǎng)直接的負荷損失、間接的工業(yè)生產(chǎn)損失等都可謂其后果。在線辨識連鎖故障序列，使用全面反映連鎖故障的指標是不現(xiàn)實的，因此，本文使用負荷損失作為嚴重性指標。負荷是電網(wǎng)供需平衡中重要一環(huán)，負荷損失量易于統(tǒng)計，并且電網(wǎng)采取的措施多與負荷有關(guān)。

3.3 序列辨識

遍歷所有線路為初始線路，以隱性故障為傳播機理搜尋下一層開斷線路（因多重隱性故障發(fā)生的幾率遠小于單重隱性故障，忽略其影響），當(dāng)電網(wǎng)損失負荷超過20%或者潮流不收斂時，判定為發(fā)生災(zāi)難性事故，連鎖故障終止。本文利用風(fēng)險理論量化連鎖故障的災(zāi)難后果，可表示為：

4 序列辨識結(jié)果輸出

將辨識出風(fēng)險值較大的前n個連鎖故障序列（n可由調(diào)度人員設(shè)定）和辨識程序運行情況寫成E語言文件，存放在指定位置，供調(diào)度員使用。

圖3為E語言輸出文件中辨識程序運行情況。第一行是電網(wǎng)名稱和文件生成時間，便于不同版本管理。統(tǒng)計信息以單列式結(jié)構(gòu)呈現(xiàn)，包含四類信息：辨識成功與否、辨識失敗原因、辨識出的序列個數(shù)以及辨識耗時。辨識成功以“Yes”顯示，同時失敗原因為“0”。辨識失敗以“No”提示，并通過“1-2”告知調(diào)度員辨識失敗原因以便調(diào)整。辨識序列個數(shù)可根據(jù)調(diào)度人員需求設(shè)置，辨識時間為模塊運行時間。

圖3 E語言輸出文件Fig.3 E-language output file

結(jié)果輸出文件給調(diào)度員提供風(fēng)險值由大到小的前n個序列，圖4為單個序列結(jié)果展示。首先列出了序列結(jié)果以及其風(fēng)險值，并在下面依次給出了每一階段斷開線路的相關(guān)信息以及斷開后的負荷損失情況，以便調(diào)度員直觀了解連鎖故障的發(fā)展過程。

圖4 風(fēng)險值較大序列Fig.4 Sequence with great risk

5 連鎖故障序列在線預(yù)警框架

連鎖故障序列在線預(yù)警框架如圖5所示。

（1）CIM/E文件正確性效驗。若文件不存在或格式不正確，則會給出相應(yīng)的錯誤提示；

（2）CIM/E文件讀取與處理。將CIM/E文件的數(shù)據(jù)讀入存儲到C/C++結(jié)構(gòu)體中。進而從不同的結(jié)構(gòu)體取出數(shù)據(jù)，形成所需要的拓撲數(shù)據(jù)和運行數(shù)據(jù)的數(shù)組，以便后續(xù)拓撲分析和潮流計算；

（3）遍歷所有線路為初始故障，以繼電保護隱性故障為傳播機制，若潮流不收斂或丟失負荷超過20%則判定連鎖故障結(jié)束?？尚判詼y度便于快速計算連鎖故障的風(fēng)險后果；

（4）將程序運行情況和辨識序列結(jié)果寫成E語言文件，保存在指定位置。E語言文件中的時間信息便于版本管理。

圖5 基于CIM/E的連鎖故障在線預(yù)警流程圖Fig.5 Online warning flow chart of cascading failure based on CIM/E

6 實際地區(qū)電網(wǎng)應(yīng)用

該地區(qū)電網(wǎng)位于我國西南部，為典型送端電網(wǎng)，圖6為其系統(tǒng)結(jié)構(gòu)圖。該系統(tǒng)有57條線路、52條母線，其中電源母線21條，負荷母線18條，母線14、15、22上負荷為外部電網(wǎng)等效而得。

圖6 地區(qū)實際電網(wǎng)Fig.6 A regional practical power grid system

辨識輸出結(jié)果中最大風(fēng)險值序列展示如圖7所示。L31由于故障等原因斷開以后，系統(tǒng)一個外送通道關(guān)閉，母線22上的負荷失去電源支撐，負荷損失嚴重。進而隱性故障造成L28斷開，系統(tǒng)解列為兩個孤島分別運行，母線14、15上的外部負荷失掉來自東北部發(fā)電機群的出力，系統(tǒng)運行狀態(tài)繼續(xù)惡化。此后系統(tǒng)重要的并行外送通道L26、L27斷開，系統(tǒng)解列為三個孤島，母線14、15上的外部負荷進一步失掉東南部的發(fā)電機群出力，僅能靠西部發(fā)電機群勉強支撐，母線14、15上的外送負荷被迫大幅減少，并且損失負荷已超過20%，災(zāi)難性事故發(fā)生。由此可見，本應(yīng)用開發(fā)識別出的連鎖故障序列結(jié)果符合地區(qū)電網(wǎng)運行實際。

圖7 風(fēng)險值最大序列Fig.7 Sequence with greatest risk

7 結(jié)束語

本文提出的基于CIM/E的連鎖故障在線辨識模塊，成功應(yīng)用于西南某地區(qū)電網(wǎng)。仿真結(jié)果表明，該應(yīng)用得到的風(fēng)險較大的連鎖故障序列符合地區(qū)電網(wǎng)運行實際情況，能夠為調(diào)度員提供采取預(yù)防措施的依據(jù)。