張詩琪，何立新，陳子敬，鐘文強，王庭剛

(貴州電網(wǎng)有限責(zé)任公司 貴陽供電局，貴州 貴陽 550000))

為了滿足電力能源可持續(xù)發(fā)展要求，光伏、水電等新能源開始接入配電網(wǎng)，達到了節(jié)能環(huán)保的目的[1]。新能源發(fā)電具有間歇性，造成配電網(wǎng)運行不穩(wěn)定的因素越來越多，各種配電風(fēng)險問題也接踵而來[2]。提升配電網(wǎng)運行的穩(wěn)定性，必須開展配電網(wǎng)風(fēng)險態(tài)勢預(yù)測。提前了解即將發(fā)生的風(fēng)險事故，制定最優(yōu)應(yīng)對方案，最大程度降低配電網(wǎng)風(fēng)險造成的負面影響。近年來，隨著人們對配電網(wǎng)運行風(fēng)險預(yù)測越來越重視，各種預(yù)測方法也得以推廣應(yīng)用。

文獻[3]通過極限學(xué)習(xí)機對歷史數(shù)據(jù)進行學(xué)習(xí)，進行風(fēng)險態(tài)勢預(yù)測。在自適應(yīng)學(xué)習(xí)過程中，應(yīng)用改進粒子群優(yōu)化算法，對不同迭代次數(shù)下的慣性權(quán)重和學(xué)習(xí)因子進行調(diào)整，確保輸出的預(yù)測結(jié)果更加準(zhǔn)確。但該方法計算復(fù)雜度較高。文獻[4]從配電網(wǎng)云翔風(fēng)險影響因素入手，建立安全態(tài)勢評價體系。在層次分析法和熵權(quán)法的共同作用下，為評價體系中每項指標(biāo)進行賦權(quán)。通過加權(quán)計算，得出配電網(wǎng)安全態(tài)勢評價值，并將其輸入深度神經(jīng)網(wǎng)絡(luò)模型進行計算，得到風(fēng)險態(tài)勢預(yù)測結(jié)果。但該預(yù)測方法穩(wěn)定性較差。文獻[5]依托信息熵為核心的大數(shù)據(jù)挖掘技術(shù)，檢測歷史數(shù)據(jù)中存在的異常數(shù)據(jù)。再通過蟻群優(yōu)化支持向量機算法進行自適應(yīng)學(xué)習(xí)，輸出風(fēng)險態(tài)勢預(yù)測結(jié)果。但是，該方法預(yù)測結(jié)果誤差較大。

考慮到常規(guī)的預(yù)測方法存在較多不足之處，文中針對配電網(wǎng)風(fēng)險態(tài)勢預(yù)測問題，提出一種以關(guān)聯(lián)規(guī)則為核心的預(yù)測方法。通過建立評估指標(biāo)體系、態(tài)勢數(shù)據(jù)處理、挖掘關(guān)聯(lián)規(guī)則、預(yù)測風(fēng)險態(tài)勢四個環(huán)節(jié)，快速得到貼合實際情況的風(fēng)險態(tài)勢預(yù)測結(jié)果。

1 以關(guān)聯(lián)規(guī)則為基礎(chǔ)設(shè)計配電網(wǎng)風(fēng)險態(tài)勢預(yù)測方法

1.1 建立配電網(wǎng)風(fēng)險評估指標(biāo)體系

深入分析每種配電故障發(fā)生的概率，以及故障造成的經(jīng)濟損失，定義配電網(wǎng)風(fēng)險評估指標(biāo)體系[6]，作為配電網(wǎng)風(fēng)險態(tài)勢預(yù)測的基礎(chǔ)環(huán)節(jié)。根據(jù)配電網(wǎng)運行特點可知，過電壓、低電壓、線路過載、失負荷，是配電網(wǎng)出現(xiàn)的主要風(fēng)險形式，文章針對這四種情況，定義四項主要風(fēng)險評估指標(biāo)。其中，配電網(wǎng)過電壓運行風(fēng)險計算公式為：

(1)

式中：t——目標(biāo)時刻；S——配電網(wǎng)運行總時間；?1——過電壓風(fēng)險；?1,t——目標(biāo)時刻過電壓風(fēng)險。

與過電壓相對應(yīng)的，配電網(wǎng)運行過程中還存在低電壓運行風(fēng)險，其計算公式可以表示為：

(2)

式中：?2——低電壓風(fēng)險；?2,t——目標(biāo)時刻t的低電壓運行風(fēng)險。

除此之外，線路過載情況和失負荷情況，是評估配電網(wǎng)風(fēng)險狀態(tài)的另外兩項關(guān)鍵指標(biāo)，二者計算公式分別為：

(3)

(4)

式中：?3——線路過載風(fēng)險；?4——失負荷風(fēng)險；?3,t、?4,t——目標(biāo)時刻的線路過載風(fēng)險和失負荷風(fēng)險。

綜合分析上述提出的風(fēng)險評估指標(biāo)，可以將綜合風(fēng)險指標(biāo)表示為：

(5)

式中：w1、w2、w3、w4表示權(quán)重系數(shù)，可以根據(jù)配電網(wǎng)實際運行要求來取值，但需要保證四項權(quán)重系數(shù)之和為1。通過上述運算，可以結(jié)合配電網(wǎng)歷史運行數(shù)據(jù)，得出風(fēng)險態(tài)勢數(shù)據(jù)。

1.2 配電網(wǎng)態(tài)勢數(shù)據(jù)預(yù)處理

采用層次化安全威脅模型[7]，可以描述配電網(wǎng)風(fēng)險的層次關(guān)系，見圖1。

圖1 配電網(wǎng)層次化安全威脅模型Fig.1 Hierarchical security threat model for distribution network

根據(jù)圖1可知，原始態(tài)勢數(shù)據(jù)中除了風(fēng)險評估數(shù)據(jù)外，還存在網(wǎng)絡(luò)攻擊數(shù)據(jù)、流量監(jiān)測數(shù)據(jù)等。考慮配電網(wǎng)面對一次嚴重的網(wǎng)絡(luò)攻擊，可能會產(chǎn)生大量重復(fù)報警，導(dǎo)致配電網(wǎng)態(tài)勢數(shù)據(jù)中存在大量冗余信息[8]，直接進行定量計算會降低預(yù)測效率，同時影響風(fēng)險態(tài)勢預(yù)測結(jié)果的準(zhǔn)確性。為此，文中提出在配電網(wǎng)風(fēng)險態(tài)勢預(yù)測過程中，設(shè)置一個態(tài)勢預(yù)處理模塊，建立一個數(shù)據(jù)庫保存聚類處理后的數(shù)據(jù)，作為關(guān)聯(lián)規(guī)則挖掘的基礎(chǔ)。為了便于分析，將每個風(fēng)險報警數(shù)據(jù)表示為以下數(shù)學(xué)公式：

E={Q，B，F(xiàn)，G，T，H}

(6)

式中：E——配電網(wǎng)運行風(fēng)險報警數(shù)據(jù)，Q——源地址，B——源端口，F(xiàn)——目的地址，G——目的端口，T——報警時間，H——報警類型。

設(shè)置一個相異度矩陣，用來計算兩個報警數(shù)據(jù)之間的相似程度，得到量化結(jié)果。通常情況下，相異度計算需要先計算報警數(shù)據(jù)各個屬性的相異度。

(7)

式中：a,b——兩個風(fēng)險報警數(shù)據(jù)；d——相異度；λ——數(shù)據(jù)中屬性數(shù)量，l——屬性。對比公式(7)計算結(jié)果和給定的相異度水平值，當(dāng)其小于相異度水平，此時兩個報警數(shù)據(jù)可以歸納為一類數(shù)據(jù)，按這種計算方式，完成對所有風(fēng)險態(tài)勢數(shù)據(jù)的聚類處理。

1.3 設(shè)計關(guān)聯(lián)規(guī)則挖掘算法

考慮到配電網(wǎng)內(nèi)包含較多電力節(jié)點，為了進行更加全面的風(fēng)險預(yù)測，文中將整個配電網(wǎng)劃分為多個子區(qū)域，分別針對每個區(qū)域的子網(wǎng)采集風(fēng)險態(tài)勢數(shù)據(jù)，并利用圖2所示的并行挖掘原理，對多個子網(wǎng)同步挖掘關(guān)聯(lián)規(guī)則。

圖2 關(guān)聯(lián)規(guī)則并行挖掘的工作模式Fig.2 Working mode of parallel mining association rules

關(guān)聯(lián)規(guī)則的挖掘本質(zhì)上是從海量風(fēng)險態(tài)勢數(shù)據(jù)中提取頻繁項集[9]，挖掘存在潛在關(guān)聯(lián)的信息，即可得到關(guān)聯(lián)規(guī)則。將風(fēng)險態(tài)勢數(shù)據(jù)描述為包含多個項的項集，每個項存在對應(yīng)的事務(wù)，在此基礎(chǔ)上可以將項集的支持度計數(shù)表示為：

ω(X)=|{f1|X≤fi,fi?G}|

(8)

式中：ω——支持度計數(shù)；X——項集；fi——項i對應(yīng)的事務(wù)；G——事務(wù)二元組，也是特定項集所有事務(wù)的集合。

根據(jù)頻繁項集提取結(jié)果，得到關(guān)聯(lián)規(guī)則。之后，計算支持度和置信度，用來描述關(guān)聯(lián)規(guī)則的強度，兩種度量形式計算公式為：

(9)

(10)

式中：A、B——兩個項集；χ——支持度；C——置信度；p——條件概率；I——包含多個事務(wù)的配電網(wǎng)風(fēng)險態(tài)勢數(shù)據(jù)集。

在事務(wù)集內(nèi)尋找滿足最小支持度、最小置信度要求的規(guī)則，作為關(guān)聯(lián)規(guī)則挖掘結(jié)果。通常情況下，挖掘關(guān)聯(lián)規(guī)則需要經(jīng)歷兩個環(huán)節(jié)，其一是總結(jié)所有滿足最小支持度閾值的項集，得到頻繁項集。再根據(jù)置信度要求，得到強關(guān)聯(lián)規(guī)則挖掘結(jié)果。

1.4 基于關(guān)聯(lián)規(guī)則獲取態(tài)勢預(yù)測結(jié)果

根據(jù)挖掘的關(guān)聯(lián)規(guī)則，可以描述不同電力節(jié)點之間的依賴關(guān)系，以此為基礎(chǔ)，在配電網(wǎng)風(fēng)險態(tài)勢預(yù)測時，只要確定某一節(jié)點的運行狀態(tài)，就可以根據(jù)關(guān)聯(lián)規(guī)則判斷下一時刻其他節(jié)點的運行狀態(tài)，從而明確配電網(wǎng)風(fēng)險態(tài)勢?；陉P(guān)聯(lián)規(guī)則進行風(fēng)險態(tài)勢預(yù)測，模糊推理是預(yù)測過程中的關(guān)鍵環(huán)節(jié)，文中采用可加性標(biāo)準(zhǔn)模型[10]，對輸入數(shù)據(jù)和關(guān)聯(lián)規(guī)則庫進行運算，預(yù)測未來配電網(wǎng)風(fēng)險態(tài)勢，見圖3。

圖3 配電網(wǎng)風(fēng)險態(tài)勢預(yù)測流程圖Fig.3 Flow chart of distribution network risk situation prediction

圖3中輸入原始數(shù)據(jù)與挖掘出的所有關(guān)聯(lián)規(guī)則進行并聯(lián)匹配，根據(jù)輸入數(shù)據(jù)與關(guān)聯(lián)規(guī)則的匹配程度，輸出對應(yīng)的項集，作為關(guān)聯(lián)規(guī)則后件。匯總所有輸出項集進行加權(quán)融合運算，得到最終的配電網(wǎng)風(fēng)險態(tài)勢預(yù)測結(jié)果。其中，通過模糊推理得到的輸出項集可以表示為：

(11)

式中：r——關(guān)聯(lián)規(guī)則；R——關(guān)聯(lián)規(guī)則集；u——輸入數(shù)據(jù)；υ——權(quán)重;α——模糊項集，αr(u)——輸入數(shù)據(jù)在模糊項集上的隸屬程度;β——關(guān)聯(lián)規(guī)則的規(guī)則后件。通過上述運算，得出配電網(wǎng)風(fēng)險態(tài)勢預(yù)測結(jié)果，作為提升配電網(wǎng)運行穩(wěn)定性的依據(jù)。

2 實驗

2.1 配電網(wǎng)模型

為了提升實驗的真實性，實驗參考IEEE33節(jié)點電力系統(tǒng)，建立一個配電網(wǎng)模型，在實驗室內(nèi)展開風(fēng)險態(tài)勢預(yù)測。配電網(wǎng)模型包括四條饋線，具體結(jié)構(gòu)見圖4。實際運行時，設(shè)置功率基準(zhǔn)值為12 000 kVA，電壓為15 kV。

從圖4可以明顯看出，四條饋線的首節(jié)點分別為1、19、23、26，尾節(jié)點分別為18、22、25和33。正常運行過程中，四條饋線的部分運行參數(shù)見表1。

圖4 IEEE33節(jié)點配電網(wǎng)模型Fig.4 IEE33 node distribution network model

表1 配電網(wǎng)模型主要參數(shù)統(tǒng)計表Tab.1 Statistical table main parameters of distribution network model

按照上述設(shè)置運行配電網(wǎng)模型，再應(yīng)用文中所提方法，針對配電網(wǎng)模型進行風(fēng)險態(tài)勢預(yù)測。

2.2 參數(shù)設(shè)置

依托于關(guān)聯(lián)規(guī)則的態(tài)勢預(yù)測方法應(yīng)用時，最小支持度、最小置信度的取值，直接影響了最終預(yù)測結(jié)果的準(zhǔn)確率。為了保證預(yù)測結(jié)果誤差更小，在實驗開始之初，運行配電網(wǎng)模型獲取風(fēng)險態(tài)勢數(shù)據(jù)，得到三個實驗數(shù)據(jù)集。探索不同參數(shù)條件下，每一個數(shù)據(jù)集的風(fēng)險態(tài)勢預(yù)測準(zhǔn)確率變化情況，得到圖5。

根據(jù)圖5(a)可知，配電網(wǎng)風(fēng)險態(tài)勢預(yù)測準(zhǔn)確率會隨著最小支持度的增加，呈現(xiàn)出先上升后下降的變化趨勢。實際變化過程中，當(dāng)最小支持度參數(shù)設(shè)置為0.6時，在三個數(shù)據(jù)集的預(yù)測準(zhǔn)確率均處于最高值，隨后準(zhǔn)確率開始逐步下降。因此實驗過程中設(shè)置最小支持度取值為0.6。

從圖5(b)中可以看出，前期最小支持度的取值變化，對預(yù)測準(zhǔn)確率的影響不大，因此直接從最小置信度取值為0.5開始研究。仔細觀察可以看出，最小置信度取值為0.85之前，預(yù)測準(zhǔn)確率雖然變化幅度較小，但總體表現(xiàn)為增長趨勢。而在0.85以后，預(yù)測準(zhǔn)確率開始出現(xiàn)明顯下滑。故實驗中設(shè)置最小置信度為0.85。

圖5 不同參數(shù)條件下風(fēng)險態(tài)勢預(yù)測準(zhǔn)確率變化Fig.5 Change of risk situation prediction accuracy under different parameters

雖然支持度和置信度取值較大時，關(guān)聯(lián)規(guī)則才能擁有較高的可信度，但可信度提升的同時規(guī)則數(shù)目也會急劇減少，導(dǎo)致基于關(guān)聯(lián)規(guī)則進行風(fēng)險態(tài)勢預(yù)測時，會遺漏部分信息，導(dǎo)致預(yù)測出現(xiàn)誤差。因此，支持度和置信度的取值，不能直接選擇最大值，而需要根據(jù)實際預(yù)測情況，選取最合適的數(shù)值。

2.3 風(fēng)險態(tài)勢預(yù)測結(jié)果

參數(shù)設(shè)置結(jié)束后，應(yīng)用文中提出的預(yù)測方法，與配電網(wǎng)模型1 d內(nèi)的運行風(fēng)險態(tài)勢進行預(yù)測，得到圖6所示的預(yù)測結(jié)果。

圖6 配電網(wǎng)風(fēng)險態(tài)勢預(yù)測結(jié)果Fig.6 Prediction results of distribution network risk situation

根據(jù)圖6可知，配電網(wǎng)模型存在的主要風(fēng)險是電壓越線風(fēng)險和失負荷風(fēng)險。其中，電壓越線風(fēng)險主要發(fā)生在17 h～24 h內(nèi)，在同一時間段，還存在嚴重的失負荷風(fēng)險。而低電壓風(fēng)險、線路過載風(fēng)險的計算值較小，可直接忽略?？傮w來看，配電網(wǎng)風(fēng)險在17 h～24 h較高，該時段需要提出針對性運維措施，其他時間風(fēng)險值較為穩(wěn)定，風(fēng)險問題出現(xiàn)概率極低。

根據(jù)上述預(yù)測結(jié)果可以看出，文中設(shè)計方法應(yīng)用后，可以得到未來時刻配電網(wǎng)風(fēng)險態(tài)勢，并準(zhǔn)確描述風(fēng)險出現(xiàn)的大概時間區(qū)段，便于提出應(yīng)用措施。

2.4 方法性能分析

進一步分析文中設(shè)計態(tài)勢預(yù)測方法的計算復(fù)雜度可知，文中采用了關(guān)聯(lián)規(guī)則算法，直接對海量數(shù)據(jù)進行處理，不需要復(fù)雜的非線性運算，使得該方法在實際應(yīng)用中必然會有較好的時間性能。因此，運行配電網(wǎng)模型獲取不同規(guī)模的數(shù)據(jù)，分別應(yīng)用文中設(shè)計方法、基于深度神經(jīng)網(wǎng)絡(luò)方法、基于大數(shù)據(jù)挖掘的方法進行預(yù)測，獲取不同方法的面對不同數(shù)據(jù)規(guī)模時的預(yù)測時間，對比不同方法的時間性能，見圖7。

圖7 不同風(fēng)險態(tài)勢預(yù)測方法的時間性能對比Fig.7 Time performance comparison of different risk situation prediction methods

根據(jù)圖7可知，隨著數(shù)據(jù)規(guī)模的增長，三種方法的預(yù)測時間均表現(xiàn)出增長趨勢，但是所提方法的預(yù)測時間明顯更短。以樣本規(guī)模為8×106為例，此時文中設(shè)計方法的預(yù)測時間為24 s，而其他兩種方法的預(yù)測時間分別為54 s、77 s。綜上所述，應(yīng)用關(guān)聯(lián)規(guī)則后，配電網(wǎng)風(fēng)險態(tài)勢預(yù)測時間縮短了55.56%、68.83%。

3 結(jié)束語

隨著配電網(wǎng)運行數(shù)據(jù)的大量增長，現(xiàn)有的風(fēng)險態(tài)勢預(yù)測方法難以快速得出預(yù)測結(jié)果。為了充分利用大數(shù)據(jù)的優(yōu)勢，提出應(yīng)用關(guān)聯(lián)規(guī)則算法設(shè)計一種具有優(yōu)越時間性能的態(tài)勢預(yù)測方法。該方法不需要進行復(fù)雜的非線性計算，而是依靠基于數(shù)據(jù)之間的關(guān)聯(lián)，預(yù)測未來時刻風(fēng)險態(tài)勢。從實驗驗證結(jié)果可以看出，所提方法更符合大數(shù)據(jù)環(huán)境，面對大規(guī)模數(shù)據(jù)時風(fēng)險態(tài)勢預(yù)測時間依舊較短。