基于社交網(wǎng)絡(luò)影響力的連鎖故障關(guān)鍵線路辨識(shí)

郭 琦， 郝乾鵬， 劉 軍， 孟凡成， 胡 博， 薛艷軍

(1.內(nèi)蒙古電力(集團(tuán))有限責(zé)任公司電力調(diào)度控制分公司， 內(nèi)蒙古 呼和浩特 010020； 2.輸配電裝備及系統(tǒng)安全與新技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室(重慶大學(xué))， 重慶 400044； 3.北京清大科越股份有限公司， 北京 100084)

1 引言

近二十年來各國(guó)電網(wǎng)的多數(shù)大停電事故均涉及輸電系統(tǒng)的連鎖故障，例如，印度大停電事故和南澳大停電事故中均出現(xiàn)了線路的連鎖跳閘，造成大規(guī)模電力負(fù)荷的長(zhǎng)時(shí)間停運(yùn)[1,2]。為了有效降低大停電事故風(fēng)險(xiǎn)，準(zhǔn)確辨識(shí)電力系統(tǒng)的薄弱環(huán)節(jié)是非常有必要的，通過運(yùn)行或規(guī)劃手段來降低可能誘發(fā)連鎖故障的輸電元件故障率。

為有效識(shí)別電力系統(tǒng)的關(guān)鍵線路，現(xiàn)有研究分別從多個(gè)角度提出了輸電元件的排序方法。一方面，部分文獻(xiàn)采用電網(wǎng)介數(shù)等結(jié)構(gòu)指標(biāo)衡量元件的重要性。文獻(xiàn)[3]通過分析電源與負(fù)荷的潮流路徑構(gòu)建加權(quán)電氣介數(shù)，用于量化線路對(duì)系統(tǒng)最大輸電能力的影響。文獻(xiàn)[4]建立了潮流轉(zhuǎn)移度指標(biāo)，可以量化分析輸電容量的安全裕度。為了進(jìn)一步考慮連鎖故障中的頻率穩(wěn)定和電壓穩(wěn)定，文獻(xiàn)[5]采用動(dòng)態(tài)潮流推導(dǎo)了線路元件的安穩(wěn)影響指標(biāo)。另一方面，為更全面地考慮連鎖故障的傳播機(jī)理，部分文獻(xiàn)采用線路故障傳播信息構(gòu)建相關(guān)性網(wǎng)絡(luò)，然后基于復(fù)雜網(wǎng)絡(luò)理論辨識(shí)關(guān)鍵線路。文獻(xiàn)[6-8]在采用潮流轉(zhuǎn)移模擬連鎖故障后，將故障鏈組合成故障網(wǎng)絡(luò)，通過統(tǒng)計(jì)故障網(wǎng)絡(luò)中各元件的出入度等網(wǎng)絡(luò)參數(shù)辨識(shí)薄弱環(huán)節(jié)。

以上研究均為辨識(shí)連鎖故障中的關(guān)鍵線路提供了有效思路，但在分析中沒有考慮輸電元件對(duì)連鎖故障的重疊影響。由于調(diào)控資源有限，電網(wǎng)運(yùn)行人員通常希望所考慮的關(guān)鍵線路集合能夠最大化體現(xiàn)連鎖故障的總影響。如果所辨識(shí)的關(guān)鍵線路對(duì)故障傳播過程的影響重疊度過高，關(guān)鍵元件的代表性會(huì)大幅降低。為此，本文基于社交網(wǎng)絡(luò)的最大影響力分析理論，建立了連鎖故障的故障傳播網(wǎng)絡(luò)，從中選取最具影響力的關(guān)鍵線路集合，可更全面地代表連鎖故障的風(fēng)險(xiǎn)源，為大停電風(fēng)險(xiǎn)管控提供決策參考。

2 輸電線路的連鎖故障模型

電力系統(tǒng)的連鎖故障涉及到潮流轉(zhuǎn)移和穩(wěn)定分析等多種復(fù)雜機(jī)理[9,10]。目前，大多數(shù)的連鎖故障模型在故障傳播過程中主要刻畫潮流轉(zhuǎn)移的風(fēng)險(xiǎn)，其中由ORNL、PSerc和Alaska三個(gè)機(jī)構(gòu)提出的OPA(ORNL-PSerc-Alaska)模型最為經(jīng)典[11]。OPA模型的核心是從重載線路中隨機(jī)抽取后續(xù)故障，模擬輸電網(wǎng)的連鎖故障。本文采用OPA模型生成連鎖故障的樣本集合，用于構(gòu)建連鎖故障的社交網(wǎng)絡(luò)圖及辨識(shí)關(guān)鍵線路。生成連鎖故障樣本的過程主要如下：

步驟1：隨機(jī)模擬負(fù)荷狀態(tài)，其中每個(gè)節(jié)點(diǎn)的負(fù)荷Di隨機(jī)生成如下：

(1)

步驟2：根據(jù)隨機(jī)抽樣后的負(fù)荷水平，采用直流最優(yōu)潮流確定系統(tǒng)運(yùn)行狀態(tài)。

步驟3：采用初始故障概率P0，隨機(jī)抽取初始故障元件。

步驟4：跳開故障線路后，采用直流最優(yōu)潮流更新系統(tǒng)運(yùn)行狀態(tài)。

步驟5：若出現(xiàn)重載線路(負(fù)載率大于閾值β)，采用故障概率P1跳開重載線路。

步驟6：若步驟5生成新的線路故障，返回步驟4；否則，至步驟7。

步驟7：輸出故障事件記錄，統(tǒng)計(jì)負(fù)荷損失。

需要說明的是，除了以O(shè)PA模型為代表的穩(wěn)態(tài)機(jī)理模型外，現(xiàn)有研究還提出了考慮電網(wǎng)安全穩(wěn)定的動(dòng)態(tài)機(jī)理模型。由于在不同的機(jī)理模型中，輸電元件的停運(yùn)均是多數(shù)大停電事件的主導(dǎo)因素，在將停運(yùn)元件劃分為不同的故障階段后，本文所提出的方法同樣適用于考慮更加復(fù)雜機(jī)理的連鎖故障模型。例如，進(jìn)一步考慮無功特性和信息系統(tǒng)的影響機(jī)理后[12,13]，所提方法仍可根據(jù)樣本分析故障傳播影響，辨識(shí)關(guān)鍵線路。

3 基于社交網(wǎng)絡(luò)影響力的關(guān)鍵線路辨識(shí)

3.1 連鎖故障的社交網(wǎng)絡(luò)圖

為分析媒體廣告、病毒等在社會(huì)群體內(nèi)的傳播模式，國(guó)內(nèi)外學(xué)者們提出了多種社交網(wǎng)絡(luò)的信息傳播模型，如獨(dú)立級(jí)聯(lián)模型和線性閾值模型等[14]。獨(dú)立級(jí)聯(lián)模型憑借通用性較好、可解釋性高等優(yōu)點(diǎn)而被廣泛采用。因此，本文采用獨(dú)立級(jí)聯(lián)模型構(gòu)建連鎖故障的社交網(wǎng)絡(luò)圖。

獨(dú)立級(jí)聯(lián)模型主要包括三個(gè)要素：節(jié)點(diǎn)、有向邊和激活概率。其中，節(jié)點(diǎn)之間通過有向邊連接，而每條邊的首節(jié)點(diǎn)被激活后，末節(jié)點(diǎn)有一定概率被進(jìn)一步激活。以邊(u,v)為例，當(dāng)節(jié)點(diǎn)u被激活后，節(jié)點(diǎn)v被已激活節(jié)點(diǎn)u激活的概率為p(u,v)。圖1給出了一個(gè)含4個(gè)節(jié)點(diǎn)的獨(dú)立級(jí)聯(lián)模型，圖中每條邊上標(biāo)注的數(shù)值即為激活概率。

圖1 獨(dú)立級(jí)聯(lián)模型示意Fig.1 Illustrative diagram of independent cascade model

當(dāng)采用社交網(wǎng)絡(luò)分析連鎖故障的傳播影響時(shí)，獨(dú)立級(jí)聯(lián)模型的節(jié)點(diǎn)對(duì)應(yīng)于線路元件，而有向邊則對(duì)應(yīng)于連鎖故障的傳播方向。為了估計(jì)不同線路元件之間的故障傳播概率，本文采用文獻(xiàn)[15]的參數(shù)估計(jì)方法，從OPA模型生成的連鎖故障鏈中，估計(jì)線路u的故障激活線路v發(fā)生故障的概率，表示為p(u,v)。在社交網(wǎng)絡(luò)中，根據(jù)不同個(gè)體被激活的時(shí)間順序，可將廣告、病毒等媒介的傳播過程劃分為不同的代，并在分析傳播影響時(shí)，一般假設(shè)每一代的個(gè)體僅被其上一代的個(gè)體所影響[14]。因此，本文中也根據(jù)連鎖故障的傳播過程將輸電元件劃分為不同階段，估計(jì)前后階段中不同元件之間的故障傳播概率，具體步驟可參見文獻(xiàn)[15]，此處不再贅述。

根據(jù)所估計(jì)的連鎖故障傳播概率p(u,v)，即可構(gòu)建連鎖故障的社交網(wǎng)絡(luò)圖，用于線路故障的影響分析。不同于文獻(xiàn)[6]中的連鎖故障圖，本文的連鎖故障網(wǎng)絡(luò)并不假設(shè)故障鏈的每個(gè)環(huán)節(jié)僅包括1個(gè)線路故障，具有更好的通用性。

3.2 基于最大影響力的關(guān)鍵線路辨識(shí)

基于社交網(wǎng)絡(luò)的影響力分析可識(shí)別信息傳播的關(guān)鍵個(gè)體，實(shí)現(xiàn)最優(yōu)的影響力控制，如最大化媒體廣告的推廣效果，或最小化病毒的傳播范圍[14]。類似地，為了更有效地控制連鎖故障的傳播，可選擇在連鎖故障傳播網(wǎng)絡(luò)中影響力最大的線路元件集合。因此，本節(jié)提出了一種基于最大影響力分析的關(guān)鍵線路辨識(shí)方法。

首先，在連鎖故障傳播網(wǎng)絡(luò)中，計(jì)算線路u對(duì)線路i的影響力I(u,i)。在社交網(wǎng)絡(luò)中，廣告信息、病毒等媒質(zhì)的傳播概率一般均較低，所以在社交影響力的分析中多考慮二跳的情形。類似地，考慮到輸電系統(tǒng)的連鎖故障為小概率事件，故障鏈的發(fā)生概率隨著傳播長(zhǎng)度迅速降低，因此本文綜合考慮故障鏈的傳播長(zhǎng)度和發(fā)生概率，計(jì)算影響力I(u,i)為：

(2)

式中，ωu和ωi分別為線路u和線路i的權(quán)重系數(shù)；N為輸電線路的總數(shù)量；式(2)右側(cè)的第一項(xiàng)為線路i直接被線路u激活的概率；第二項(xiàng)為線路j受線路u影響下經(jīng)二跳傳播后被激活的概率。

為描述線路u對(duì)傳播連鎖故障的作用，可進(jìn)一步定義線路u的總影響力R(u)為：

(3)

記連鎖故障傳播中的關(guān)鍵線路集合為S。在選定集合S時(shí)，若不考慮關(guān)鍵線路的影響力重疊效應(yīng)，可采用式(4)計(jì)算集合S在連鎖故障社交網(wǎng)絡(luò)上的總影響力f(S)：

(4)

由于社交網(wǎng)絡(luò)的最大影響力分析屬于非決定性多項(xiàng)式(Non-deterministic Polynomial，NP)問題[14]，在選擇令f(S)最大化的關(guān)鍵線路時(shí)，常采用貪心算法求解。具體地，首先初始化S為空集，然后采用貪心策略逐個(gè)添加線路至S，在每次選擇線路時(shí)，均保證使當(dāng)前S的總影響力f(S)增長(zhǎng)量最大，即表示為：

(5)

由于式(2)～式(5)的計(jì)算過程不考慮影響力的重疊，關(guān)鍵線路辨識(shí)結(jié)果易出現(xiàn)“富人俱樂部效應(yīng)”，即影響力較大的線路容易存在緊密連接，導(dǎo)致集合S中的關(guān)鍵線路可能發(fā)生聚集，使得S無法更有效地表征連鎖故障的傳播，進(jìn)而影響故障阻斷的決策效果。

為考慮連鎖故障傳播中線路影響的重疊性，本文提出一種新的影響力計(jì)算方法，具體如下：

首先，根據(jù)當(dāng)前的關(guān)鍵線路集S更新線路k(k不屬于S)的權(quán)重為：

(6)

式(6)表示根據(jù)當(dāng)前集合S內(nèi)線路的影響范圍，降低線路k的權(quán)重值ωk。其含義為：當(dāng)S內(nèi)的線路對(duì)線路k存在較高影響時(shí)，可以通過集合S等效控制連鎖故障中線路k的故障。以一種特殊情況為例，若在連鎖故障樣本集中線路k的所有故障均由集合S內(nèi)的線路傳播導(dǎo)致，那么根據(jù)式(6)計(jì)算的ω′k近似為0。這使得在選擇新的線路時(shí)，即使線路k的故障頻率很高，也不會(huì)優(yōu)先考慮線路k加入S。

同時(shí)，在計(jì)算線路k對(duì)其他線路的影響時(shí)，對(duì)集合S內(nèi)的線路u，更新相應(yīng)的傳播概率為：

p′(k,u)=(1-α)p(k,u) ?u∈S

(7)

式中，α為取值在[0,1]的影響折扣因子。

式(7)表示當(dāng)選擇線路u加入集合S后，在社交網(wǎng)絡(luò)上降低附近線路對(duì)線路u的影響能力，以保證最終選擇的關(guān)鍵線路具有更全面的代表性。

每次選擇新的線路后，分別根據(jù)式(6)和式(7)更新線路的權(quán)重系數(shù)和故障傳播概率。采用貪心策略[16]重復(fù)k次以上線路選擇過程，可得到對(duì)連鎖故障傳播影響最大的Top-k線路集合，同時(shí)大幅降低Top-k線路的影響范圍重疊。

3.3 關(guān)鍵線路辨識(shí)算法流程

基于3.1節(jié)的社交網(wǎng)絡(luò)和3.2節(jié)的影響力分析，即可選擇電力系統(tǒng)連鎖故障的關(guān)鍵線路集合。圖2給出了基于社交網(wǎng)絡(luò)影響力分析的連鎖故障關(guān)鍵線路辨識(shí)算法流程。

圖2 基于社交網(wǎng)絡(luò)影響力的連鎖故障關(guān)鍵線路辨識(shí)流程Fig.2 Flow chart for critical branch identification of cascading failure based on influence analysis of social network

4 算例分析

本節(jié)采用某省級(jí)電網(wǎng)驗(yàn)證所提出算法的有效性。該省級(jí)電網(wǎng)共有110 kV以上線路904條。采用OPA模型生成200 000個(gè)連鎖故障樣本。其中，參照文獻(xiàn)[13]，第2節(jié)中OPA模型的參數(shù)γ、β、P0和P1取值分別為2.0、0.90、0.001和0.30。在計(jì)算影響力時(shí)，所有線路的權(quán)重ω均設(shè)為1。

4.1 連鎖故障的社交網(wǎng)絡(luò)圖及關(guān)鍵線路辨識(shí)

首先，考慮該電網(wǎng)的某典型運(yùn)行場(chǎng)景A，其中負(fù)荷水平為28 862 MW，新能源出力為3 416 MW。采用連鎖故障的樣本集合生成連鎖故障的社交網(wǎng)絡(luò)，如圖3所示。圖3中，每個(gè)頂點(diǎn)均對(duì)應(yīng)于電力系統(tǒng)的一個(gè)輸電線路，而每一條有向邊均對(duì)應(yīng)于一種故障傳播鏈接，即在該有向邊的首端頂點(diǎn)故障后，末端頂點(diǎn)發(fā)生連鎖故障。

圖3 某實(shí)際電網(wǎng)連鎖故障的社交網(wǎng)絡(luò)Fig.3 Social network for cascading failure of actual grid

從圖3中可見，盡管連鎖故障涉及到大量的線路，只有少量線路處于社交網(wǎng)絡(luò)的中心，對(duì)連鎖故障的傳播起關(guān)鍵作用。分別采用以下4種方法辨識(shí)連鎖故障的關(guān)鍵線路集合：

方法1：考慮不同線路的影響重疊性，采用所提方法辨識(shí)影響力排序前20的關(guān)鍵線路，其中折扣因子設(shè)為0。

方法2：不考慮不同線路的影響重疊性，采用所提方法辨識(shí)影響力排序前20的關(guān)鍵線路。

方法3：參考文獻(xiàn)[6]，以激活概率p(u,v)為權(quán)重，采用加權(quán)出度選取排序前20的關(guān)鍵線路，即易影響周邊元件的關(guān)鍵線路。

方法4：參考文獻(xiàn)[6]，同樣以激活概率p(u,v)為權(quán)重，采用加權(quán)入度選取排序前20的關(guān)鍵線路，即易受周邊元件影響的關(guān)鍵線路。表1對(duì)比了不同方法辨識(shí)的關(guān)鍵影響線路。

表1 關(guān)鍵線路辨識(shí)結(jié)果對(duì)比Tab.1 Comparison of identified critical branches

從表1中可見，在所辨識(shí)的關(guān)鍵線路中，方法2和方法3的重合度較高，且結(jié)果與方法4差異較大。這是因?yàn)榉椒?和方法3均主要量化線路對(duì)其他輸電元件的故障影響，而方法4主要辨識(shí)易受其他元件影響的關(guān)鍵線路。另一方面，雖然方法1也量化線路對(duì)其他元件的故障影響，但由于計(jì)及了線路影響的重疊，方法1的影響力排序相較于方法2發(fā)生明顯變化。

為了對(duì)比方法1和方法2的結(jié)果差異，圖4和圖5分別給出了所辨識(shí)關(guān)鍵線路在圖3所示的社交網(wǎng)絡(luò)上的連接關(guān)系。

圖5 方法2所辨識(shí)關(guān)鍵線路在社交網(wǎng)絡(luò)上的連接關(guān)系Fig.5 Connection of critical branched identified by method 2 on social network

對(duì)比圖4和圖5的網(wǎng)絡(luò)連接關(guān)系可以發(fā)現(xiàn)，方法1所辨識(shí)關(guān)鍵線路之間的耦合關(guān)系相對(duì)稀疏。其中，連接到圖4的每個(gè)節(jié)點(diǎn)的平均邊數(shù)為3.10，而圖5則為9.95。結(jié)果表示：方法2的辨識(shí)結(jié)果中存在明顯的“富人俱樂部效應(yīng)”，即由于緊密連接關(guān)系的存在，激活圖5中所有節(jié)點(diǎn)所要求的故障數(shù)量遠(yuǎn)低于圖4。同時(shí)，圖4中線路L417、L535和L629與其他17條關(guān)鍵線路之間不存在連接關(guān)系，這表明在故障傳播中，其余17條線路的故障對(duì)這3條線路產(chǎn)生的影響較小。

以方法1中排序前20的關(guān)鍵線路為例，圖6對(duì)比了考慮重疊性前后的關(guān)鍵線路影響力。從圖6中可見，計(jì)及影響重疊效應(yīng)時(shí)，排序靠后的線路影響力有所降低，且降低幅度差異較大。例如，若不計(jì)及重疊性，線路L150的影響力排序?yàn)榈?，但由于L150與排序第1的L774的影響范圍重疊度較高，在計(jì)及L774后L150的重要程度會(huì)有較明顯的降幅，線路L198的關(guān)鍵性會(huì)超過線路L150。

圖6 考慮影響重疊性前后的關(guān)鍵線路影響力對(duì)比Fig.6 Comparison of influence of critical branches before and after consideration of influence overlapping

最后，為驗(yàn)證所辨識(shí)關(guān)鍵線路的有效性，依次選取前k(k=1，2,…，20)個(gè)關(guān)鍵線路作為連鎖故障的初始故障，模擬連鎖故障的傳播情況。對(duì)每個(gè)初始故障集，均模擬1 000次連鎖故障，并統(tǒng)計(jì)平均線路停運(yùn)數(shù)量作為影響力衡量指標(biāo)。圖7給出了對(duì)應(yīng)的平均線路停運(yùn)數(shù)量結(jié)果。

圖7 關(guān)鍵線路作為初始故障引起的平均線路停運(yùn)數(shù)量Fig.7 Average number of branch outages with critical branches as initial outages

從圖7中可見，當(dāng)初始故障選擇的關(guān)鍵線路數(shù)量較少時(shí)，方法1、方法2和方法3所誘發(fā)連鎖故障的數(shù)量比較接近，均高于方法4誘發(fā)的線路停運(yùn)數(shù)。這是因?yàn)榉椒?識(shí)別的線路更易受其他故障影響，而非誘發(fā)新的線路停運(yùn)。但是，當(dāng)初始故障數(shù)量大于10以后，所提方法誘發(fā)的線路停運(yùn)數(shù)量明顯高于其余方法。特別地，當(dāng)初始故障數(shù)量為20時(shí)，方法2和方法3所引起線路停運(yùn)總數(shù)也接近20，而方法1則為30。這表明由于“富人俱樂部效應(yīng)”的存在，方法2和方法3的關(guān)鍵線路構(gòu)成緊密聯(lián)系，這些線路內(nèi)部存在較高的故障傳播概率，僅部分線路故障即可激活這些關(guān)鍵線路，但是對(duì)系統(tǒng)內(nèi)其余線路的影響則并非最大。方法1考慮了影響重疊性對(duì)關(guān)鍵線路辨識(shí)的遮蔽效應(yīng)，可更全面地辨識(shí)連鎖故障傳播網(wǎng)絡(luò)的關(guān)鍵線路集合。

4.2 影響折扣因子的靈敏度分析

為分析影響折扣因子α對(duì)關(guān)鍵線路辨識(shí)結(jié)果的影響，分別將α設(shè)為0、0.1、0.2、0.3、0.4、0.5、0.6、0.7、0.8、0.9和1.0。在不同的折扣因子取值下，依次辨識(shí)關(guān)鍵線路，并按照4.1節(jié)的方法驗(yàn)證所辨識(shí)關(guān)鍵線路的有效性。將以上11個(gè)場(chǎng)景依次記錄為場(chǎng)景1～場(chǎng)景11。另外，為了說明式(6)更新線路影響權(quán)重對(duì)關(guān)鍵線路辨識(shí)的影響，在折扣因子α設(shè)為1.0的前提下，去除式(6)的計(jì)算過程，并將其記為場(chǎng)景12。

圖8中對(duì)比了結(jié)果中對(duì)應(yīng)的平均線路停運(yùn)數(shù)量。從圖8可見，當(dāng)考慮α取不同值的時(shí)候，由于場(chǎng)景1～場(chǎng)景11均采用式(6)更新了線路權(quán)重，各個(gè)場(chǎng)景下誘發(fā)的連鎖故障數(shù)量較為接近，而且α取較大的折扣系數(shù)時(shí)連鎖故障數(shù)量更加穩(wěn)定。這是因?yàn)樵诟戮€路權(quán)重后，場(chǎng)景1～場(chǎng)景11所辨識(shí)的線路集合重合度很大，主要體現(xiàn)為線路之間的相對(duì)排序發(fā)生小幅換位。

圖8 不同折扣因子下平均線路停運(yùn)數(shù)量Fig.8 Average number of branch outages under different discount factors

進(jìn)一步結(jié)合場(chǎng)景12的結(jié)果可以發(fā)現(xiàn)，如果不考慮式(6)對(duì)線路影響權(quán)重的更新，最后的連鎖故障數(shù)量明顯低于場(chǎng)景11。同時(shí)，場(chǎng)景1～場(chǎng)景12中所有的連鎖故障數(shù)量均高于圖7中不考慮影響重疊性的情形。這表明所提出的線路影響量化策略可以有效地提升薄弱環(huán)節(jié)的辨識(shí)效果。

4.3 多運(yùn)行場(chǎng)景下的電網(wǎng)關(guān)鍵線路辨識(shí)分析

為分析新能源滲透率、負(fù)荷水平等典型運(yùn)行場(chǎng)景差異的影響，除了4.1節(jié)中的運(yùn)行場(chǎng)景A外，考慮另外兩個(gè)典型場(chǎng)景B和場(chǎng)景C。其中，場(chǎng)景B和場(chǎng)景C的負(fù)荷水平分別為24 828 MW和29 506 MW，新能源出力水平則分別為2 273 MW和10 560 MW。

對(duì)場(chǎng)景A、場(chǎng)景B、場(chǎng)景C分別生成連鎖故障樣本后，在α設(shè)為0的前提下，依次辨識(shí)關(guān)鍵線路集合。圖9中對(duì)比了不同運(yùn)行場(chǎng)景下關(guān)鍵線路誘發(fā)的平均線路停運(yùn)數(shù)量。

圖9 不同運(yùn)行場(chǎng)景下平均線路停運(yùn)數(shù)量Fig.9 Average number of branch outages under different operation scenarios

由圖9可見，三個(gè)運(yùn)行場(chǎng)景下發(fā)生連鎖故障的風(fēng)險(xiǎn)明顯不同。這是由于關(guān)鍵線路的辨識(shí)不僅由電網(wǎng)的電氣結(jié)構(gòu)決定，還與電網(wǎng)的實(shí)時(shí)運(yùn)行狀態(tài)相關(guān)。特別地，場(chǎng)景C中平均線路停運(yùn)數(shù)量分布的斜率近似為1，即對(duì)應(yīng)的初始故障未誘發(fā)新的連鎖故障。這是因?yàn)楸M管場(chǎng)景C的負(fù)荷水平最高，但新能源處于大發(fā)狀態(tài)，系統(tǒng)可調(diào)度電源充裕度明顯高于場(chǎng)景A和場(chǎng)景B，而在OPA模型中，每次故障后調(diào)度全網(wǎng)機(jī)組實(shí)現(xiàn)功率再平衡，因此場(chǎng)景C的輸電網(wǎng)潮流分布更均勻，降低了連鎖故障風(fēng)險(xiǎn)。此外，場(chǎng)景A和場(chǎng)景B的新能源出力接近，但場(chǎng)景A下的負(fù)荷水平高于場(chǎng)景B，因此多數(shù)情況下場(chǎng)景A的連鎖故障風(fēng)險(xiǎn)高于場(chǎng)景B。

5 結(jié)論

為有效分析輸電系統(tǒng)連鎖故障的傳播特性，本文基于社交網(wǎng)絡(luò)理論建立了連鎖故障的傳播網(wǎng)絡(luò)，通過計(jì)及重疊效應(yīng)的社交網(wǎng)絡(luò)影響力分析，提出了一種關(guān)鍵線路辨識(shí)方法。以某實(shí)際電網(wǎng)為例，基于OPA模型的連鎖故障樣本建立連鎖故障傳播網(wǎng)絡(luò)，關(guān)鍵線路的辨識(shí)結(jié)果表明：在連鎖故障的傳播過程中，僅少量的線路對(duì)故障傳播存在支配性影響；在考慮不同關(guān)鍵線路影響的重疊性后，所提方法可以更全面地表征故障的傳播模式，有效降低“富人俱樂部效應(yīng)”導(dǎo)致的遮蔽效應(yīng)。

下一步研究將考慮建立多樣化的折扣指標(biāo)和更復(fù)雜的連鎖故障機(jī)理，在連鎖故障關(guān)鍵線路辨識(shí)中更全面地刻畫輸電元件的影響力。