夏 芃，張 倩，王群京，王 璨

(1.安徽大學(xué)電氣工程及自動化學(xué)院，合肥 230601；2.安徽大學(xué)工業(yè)節(jié)電與電能質(zhì)量控制協(xié)同創(chuàng)新中心，合肥 230601;3.教育部電能質(zhì)量工程研究中心(安徽大學(xué))，合肥 230601；4.工業(yè)節(jié)電與用電安全安徽省重點實驗室(安徽大學(xué))，合肥 230601；5.國網(wǎng)安徽省電力有限公司，合肥 230601)

分布式電源并網(wǎng)導(dǎo)致配電網(wǎng)結(jié)構(gòu)趨于復(fù)雜，配電網(wǎng)電壓波動過大[1]、網(wǎng)損較大[2]等問題日趨嚴(yán)重。因此，如何保障電力系統(tǒng)穩(wěn)定且安全運行，并減小配電網(wǎng)電壓波動和網(wǎng)損是國內(nèi)外學(xué)者研究的熱點。配電網(wǎng)無功優(yōu)化是指在電力系統(tǒng)能穩(wěn)定且安全運行的前提情況下通過針對配電網(wǎng)發(fā)電機端電壓的調(diào)節(jié)[3]、變壓器分接頭次數(shù)[4]以及增加無功補償裝置[5-6]等方式來減少電壓波動并降低網(wǎng)損。

為解決大規(guī)模分布式電源并網(wǎng)所引發(fā)的相關(guān)問題，目前中外學(xué)者主要從改變無功優(yōu)化模型以及無功優(yōu)化算法等方面進行研究。文獻[7]通過建立系統(tǒng)運行成本和用戶滿意度為目標(biāo)的配電網(wǎng)無功優(yōu)化模型，提高了計算精度和速度。文獻[8]建立了混合整數(shù)半定規(guī)劃模型，先將傳統(tǒng)無功優(yōu)化模型轉(zhuǎn)化成半定規(guī)劃模型，在此基礎(chǔ)上增加離散變量使其變成混合整數(shù)半定規(guī)劃模型，使配電網(wǎng)無功優(yōu)化模型的求解更加方便。上述改進雖然和配電網(wǎng)原始模型相比計算速度得到提高，但改進的方法過于依賴配電網(wǎng)的模型，當(dāng)配電網(wǎng)的模型過于復(fù)雜時，該方法存在局限性。文獻[9]將禁忌搜索算法和粒子群算法相融合，對傳統(tǒng)的粒子群算法進行改進，并運用到配電網(wǎng)無功優(yōu)化策略當(dāng)中。文獻[10]通過構(gòu)造初始信息矩陣和蟻群算法相結(jié)合，提高了蟻群算法的計算效率，上述算法改進雖然和傳統(tǒng)算法相比尋優(yōu)精度和計算精度有所提高，但仍未擺脫算法計算時間長、步驟較為復(fù)雜的特點。

近些年來，大數(shù)據(jù)技術(shù)在電力系統(tǒng)中的應(yīng)用越來越廣泛，從而給配電網(wǎng)無功優(yōu)化的方法帶來了一種新的途徑。目前大數(shù)據(jù)技術(shù)大多用于電力系統(tǒng)負(fù)荷預(yù)測以及故障診斷領(lǐng)域[11-13]，但是部分學(xué)者也將大數(shù)據(jù)技術(shù)和配電網(wǎng)無功優(yōu)化應(yīng)用領(lǐng)域相結(jié)合。文獻[14]運用大數(shù)據(jù)技術(shù)對風(fēng)力發(fā)電廠并入光伏后的電壓進行無功控制。文獻[15]將支持向量機的方法運用到遺傳算法當(dāng)中并加以改進，得出了基于數(shù)據(jù)挖掘技術(shù)的無功優(yōu)化策略。文獻[16]將負(fù)荷分布匹配與熵權(quán)理論相結(jié)合，運用熵權(quán)法從歷史優(yōu)化方案中選擇最優(yōu)無功優(yōu)化控制方案。文獻[17]將粒子群算法對極限學(xué)習(xí)機進行優(yōu)化，并運用到配電網(wǎng)的調(diào)壓策略中。上述基于數(shù)據(jù)挖掘技術(shù)的無功優(yōu)化控制策略對配電網(wǎng)具體模型依賴程度低，均取得了較好的無功優(yōu)化控制效果。

配電網(wǎng)無功優(yōu)化是典型的高維非線性函數(shù)問題，而深度學(xué)習(xí)理論可直接確定輸入和輸出之間的隱含關(guān)系。其中，深度置信網(wǎng)絡(luò)(deep belief networks，DBN)的網(wǎng)絡(luò)模型結(jié)構(gòu)較為簡單且訓(xùn)練難度較低，常用于求解高維非線性函數(shù)問題?，F(xiàn)利用配電網(wǎng)運行的光伏、負(fù)荷、溫度、光照強度以及風(fēng)速5個數(shù)據(jù)構(gòu)造隨機矩陣構(gòu)成樣本集，運用粒子群優(yōu)化 (particle swarm optimization,PSO)對DBN的初始權(quán)重進行優(yōu)化，通過建立PSO-DBN模型得出配電網(wǎng)運行特征與無功優(yōu)化策略的函數(shù)關(guān)系，將本文方法用于改進后的IEEE33節(jié)點，并和粒子群優(yōu)化方法、改進的粒子群優(yōu)化方法以及DBN對比，以驗證本文方法對電壓波動和網(wǎng)損的調(diào)節(jié)效果。

1 隨機矩陣?yán)碚?/h2>

隨機矩陣?yán)碚揫18]是一種常用的大數(shù)據(jù)分析方法，其主要應(yīng)用于故障診斷以及變電站二次評估。當(dāng)隨機矩陣的行數(shù)和列數(shù)的比值保持恒定時，其經(jīng)驗分布函數(shù)滿足單環(huán)定律、半圓律等。主要以單環(huán)定律進行研究。

1.1 單環(huán)定理

i=1,2,…,N,j=1,2,…,T

(1)

(2)

式(2)中：U為酉矩陣。

(3)

根據(jù)式(4)將矩陣積標(biāo)準(zhǔn)化得

(4)

(5)

式(5)中：λ為矩陣的特征值。

1.2 線性特征值統(tǒng)計量

線性特征值統(tǒng)計量是對隨機矩陣特征值分布的特點具體形式的表現(xiàn),而線性特征值統(tǒng)計量定義為

(6)

式(6)中：s為統(tǒng)計函數(shù)對應(yīng)的映射值；λi(i=1,2,…,n)為隨機矩陣的特征值；φ(·)為一個線性特征統(tǒng)計函數(shù)，不同的線性統(tǒng)計函數(shù)的選擇即可得到不同的線性特征統(tǒng)計量。

1.3 特征指標(biāo)的提取

采用光伏、負(fù)荷、溫度、環(huán)境以及風(fēng)速共5種數(shù)據(jù)，每種數(shù)據(jù)用單環(huán)定理算出隨機矩陣的特征值，然后通過特征值計算每種隨機矩陣的平均譜半徑、最大/最小譜半徑、圓環(huán)外/上/內(nèi)特征根的分布概況。隨機矩陣的模以及協(xié)方差共8個統(tǒng)計特征，因而5種數(shù)據(jù)共構(gòu)成40種統(tǒng)計特征。每個統(tǒng)計特征的具體公式為

(7)

(8)

(9)

(10)

(11)

(12)

式中：N為矩陣特征值的個數(shù)；rMSR、Rmax、Rmin分別為平均譜半徑、最大以及最小譜半徑；P1、P2、P3分別為圓環(huán)外/上/內(nèi)的特征值占總特征值比例；S1、S2、S3分別為分布在圓環(huán)外、上、內(nèi)的特征值個數(shù)。

(13)

(14)

2 DBN理論

DBN由多層受限玻爾茲曼機(restricted Boltzmann machine,RBM)和最外層的BP網(wǎng)絡(luò)構(gòu)造而成[19]，具體結(jié)構(gòu)如圖1所示。

圖1 DBN結(jié)構(gòu)Fig.1 DBN structure

DBN的訓(xùn)練步驟主要有RBM的預(yù)訓(xùn)練以及DBN參數(shù)反向微調(diào)。預(yù)訓(xùn)練即通過輸入自下而上逐層訓(xùn)練RBM，如此進行反復(fù)訓(xùn)練，不斷優(yōu)化網(wǎng)絡(luò)中的模型參數(shù)，使網(wǎng)絡(luò)達到局部最優(yōu)。而參數(shù)的反向微調(diào)即以DBN模型的最后一層為輸入，運用BP神經(jīng)網(wǎng)絡(luò)自下而上地微調(diào)整個DBN的網(wǎng)絡(luò)模型參數(shù)，從而實現(xiàn)全局參數(shù)最優(yōu)。

2.1 RBM預(yù)訓(xùn)練

RBM作為概率神經(jīng)網(wǎng)絡(luò)的分支，主要由可視層以及隱藏層所構(gòu)成的?？梢晫觱=(v1,v2,…,vm)以及隱含層h=(h1,h2,…,hn)的聯(lián)合概率能量函數(shù)為

(15)

式(15)中：θ={ω,a,b}為RBM網(wǎng)絡(luò)模型相關(guān)參數(shù)；vi和ai分別為可視層神經(jīng)元的當(dāng)前狀態(tài)以及該層網(wǎng)絡(luò)的偏置；hj和bj分別為隱含層神經(jīng)元的當(dāng)前狀態(tài)以及偏置；ωij為每層之間所需連接的權(quán)值。

在RBM中定義任意一組關(guān)于v和h的聯(lián)合概率分布為

(16)

而由聯(lián)合概率分布可得v和h的條件分布為

(17)

(18)

RBM的訓(xùn)練的最終目的即找到合適的參數(shù)θ={ω,a,b}，進而使RBM的能量誤差函數(shù)達到最小值。

2.2 DBN的反向微調(diào)

當(dāng)每層的RBM參數(shù)調(diào)整完成之后，卻很難保證整體DBN的參數(shù)達到最佳。因此，需要運用BP算法針對整個DBN網(wǎng)絡(luò)參數(shù)進行微調(diào)，具體微調(diào)過程是利用反向誤差將權(quán)值進行更新，因而獲得了整體DBN的最佳網(wǎng)絡(luò)參數(shù)。

在參數(shù)調(diào)整過程當(dāng)中，將RBM訓(xùn)練的最后一層作為輸入，利用BP算法對RBM預(yù)訓(xùn)練過程中產(chǎn)生的訓(xùn)練誤差自下而上傳到RBM層中，進而調(diào)節(jié)整個DBN網(wǎng)絡(luò)，大大減少了訓(xùn)練時間。反向微調(diào)過程中參數(shù)更新公式為

(19)

式(19)中：ωn為每層網(wǎng)絡(luò)之間的權(quán)值；n為RBM更新迭代次數(shù)；η為學(xué)習(xí)率；m為動量系數(shù)。

3 基于PSO-DBN無功優(yōu)化策略

3.1 粒子群優(yōu)化算法及其改進

在建立DBN模型之初，DBN網(wǎng)絡(luò)初始權(quán)重是通過隨機賦值所得出，若賦值不恰當(dāng)會使得DBN模型在訓(xùn)練中會出現(xiàn)局部收斂現(xiàn)象,因而采用粒子群算法對DBN網(wǎng)絡(luò)初始權(quán)重進行優(yōu)化。

粒子群優(yōu)化算法作為一種啟發(fā)式算法，其算法本身是通過對每個粒子的速度Vi和位置Xi進行更新，當(dāng)滿足結(jié)束條件時，即可獲得最優(yōu)解，粒子的速度以及位置更新公式為

(20)

(21)

式中：m=1,2,…,M；ω為粒子的權(quán)值；c1和c2為學(xué)習(xí)因子；r1和r2為在0～1所取的隨機數(shù)；t為迭代次數(shù)；Pim和Gim為個體以及群體極值。

因為傳統(tǒng)的粒子群算法易陷入局部最優(yōu)解的狀況。因此，對其改進的方法也有很多。目前大多數(shù)學(xué)者主要采用調(diào)整權(quán)重、限制學(xué)習(xí)因子等方法對粒子群算法進行改進。

主要使用對慣性權(quán)重的改進，通過線性權(quán)重遞減法，具體公式為

(22)

式(22)中：ωmax、ωmin分別為粒子群算法初始權(quán)重的最大值以及最小值；T為最高迭代次數(shù)；k為當(dāng)前迭代次數(shù)。

3.2 PSO優(yōu)化DBN網(wǎng)絡(luò)的無功優(yōu)化策略

基于PSO-DBN網(wǎng)絡(luò)的無功優(yōu)化流程圖如圖2所示，其具體優(yōu)化步驟如下。

圖2 PSO-DBN無功優(yōu)化流程圖Fig.2 PSO-DBN reactive power optimization flowchart

步驟1對配電網(wǎng)歷史時刻的負(fù)荷、光伏、環(huán)境、溫度以及風(fēng)速分別構(gòu)建各時刻的5種隨機矩陣。

步驟2根據(jù)前文所述的特征指標(biāo)提取方法，提取出40個歷史特征指標(biāo)集，進行歸一化處理并得到樣本訓(xùn)練集和測試集。

步驟3確定DBN網(wǎng)絡(luò)中隱含層層數(shù)及隱含層中神經(jīng)元的個數(shù)，以便對粒子維度進行確定。

步驟4對粒子群各參數(shù)進行設(shè)置，即粒子的群體規(guī)模、學(xué)習(xí)因子、慣性權(quán)重以及最大迭代次數(shù)。并將DBN網(wǎng)絡(luò)中各層之間的連接權(quán)值作為PSO的向量。

步驟5計算粒子的適應(yīng)度函數(shù)值f。其函數(shù)為

(23)

式(23)中：N為樣本數(shù)目；m為粒子維數(shù)；pij和tij分別為第i個樣本中第j維粒子的重構(gòu)值和實際值。f和個體極值Pbest進行比較，若f>Pbest，則將當(dāng)前適應(yīng)度值替換掉,否則將個體極值保留。

步驟6將每個粒子的適應(yīng)度值大小f和群體極值gbest相比較，若f>gbest，則將當(dāng)前適應(yīng)度值替換掉,否則將群體極值保留。

步驟7根據(jù)式(20)和式(21)將粒子的速度以及位置進行更新。

步驟8利用PSO優(yōu)化后獲得的群體極值的各維數(shù)值作為DBN網(wǎng)絡(luò)的初始權(quán)重。再進行DBN模型的訓(xùn)練和參數(shù)微調(diào)，直至DBN訓(xùn)練結(jié)束，DBN的網(wǎng)絡(luò)模型建立，從而得到無功優(yōu)化策略。

4 算例分析

4.1 算例介紹及樣本數(shù)據(jù)集的構(gòu)造

以IEEE33節(jié)點系統(tǒng)作為研究對象，對系統(tǒng)的22節(jié)點和33節(jié)點增加無功補償，在16-17支路、19-20支路、24-25支路、26-27支路進行變壓器分接頭變比的調(diào)節(jié)，TA1～TA4分別為4個載調(diào)壓器，SVC代表無功補償，其改進后的節(jié)點拓?fù)鋱D如圖3所示。

圖3 改進IEEE33節(jié)點拓?fù)浣Y(jié)構(gòu)圖Fig.3 Improved IEEE33 node topology diagram

配電網(wǎng)運行狀態(tài)數(shù)據(jù)源于安徽省安慶市電網(wǎng)公司數(shù)據(jù)集。利用配電網(wǎng)兩個月的負(fù)荷、光伏、環(huán)境、溫度以及風(fēng)速構(gòu)造5種隨機矩陣。研究每小時配電網(wǎng)的無功優(yōu)化策略，構(gòu)造5種不同的高維隨機矩陣，提取特征統(tǒng)計量作為歷史輸入，從當(dāng)?shù)仉娋W(wǎng)數(shù)據(jù)庫中得到兩個月中每小時的無功優(yōu)化控制策略作為歷史輸出，構(gòu)成1 440個歷史樣本的輸入和輸出標(biāo)簽,放入PSO-DBN模型中進行訓(xùn)練，訓(xùn)練過程中設(shè)置DBN網(wǎng)絡(luò)共有4層結(jié)構(gòu)，兩層隱含層神經(jīng)元的數(shù)量可分別設(shè)為20和10。粒子群優(yōu)化算法參數(shù)設(shè)置如下：ωmax=0.9，ωmin=0.5，學(xué)習(xí)因子c1=c2=2,T=1 000。接著，選取前文所述兩個月之外的某日作為待優(yōu)化時刻，利用隨機矩陣?yán)碚摣@取該天每小時共24個特征統(tǒng)計量作為樣本輸入，再將輸入運用到訓(xùn)練好的DBN模型中得到當(dāng)前時刻的無功優(yōu)化策略。圖4和圖5表明了該日通過本文方法得出改進后IEEE33節(jié)點系統(tǒng)的變壓器調(diào)節(jié)擋位和無功補償量。再將所得策略進行仿真計算。

圖4 IEEE33節(jié)點各支路變壓器檔位調(diào)節(jié)圖Fig.4 IEEE33 node branch transformer tap adjustment diagram

圖5 22節(jié)點和33節(jié)點無功補償Fig.5 Reactive compensation of22 and 33 node

4.2 無功優(yōu)化結(jié)果分析

為更好地驗證本文方法在無功優(yōu)化方面的有效性，將本文方法與PSO算法、改進PSO算法、DBN這3種優(yōu)化方法調(diào)壓所得的電壓波動和網(wǎng)損進行對比，以系統(tǒng)降損率和系統(tǒng)電壓偏差波動率為指標(biāo)，可更直觀地看出本文方法調(diào)壓策略的有效性。其中定義某時刻系統(tǒng)降損率為

(24)

式(24)中：eL為系統(tǒng)降損率；f0為系統(tǒng)該時刻未使用任何調(diào)壓方法所得的網(wǎng)損；f為系統(tǒng)該時刻使用調(diào)壓方法所得的網(wǎng)損。

定義某時刻的系統(tǒng)電壓偏差波動率為

(25)

式(25)中：eu為系統(tǒng)電壓偏差波動率；ΔUR為系統(tǒng)該時刻運用調(diào)壓算法調(diào)壓所得的電壓值；ΔUc為系統(tǒng)該時刻未使用任何調(diào)壓方法所得的電壓值。

運用式(24)和式(25)對無功優(yōu)化結(jié)果分析，對于待優(yōu)化的24個時刻，4種優(yōu)化方法的系統(tǒng)降損率如圖6和表1所示，系統(tǒng)電壓偏差波動率如圖7和表2所示。

圖6 降損率比較Fig.6 Comparison of loss reduction rate curve

表1 降損率分析Table 1 Analysis of loss reduction rate

圖7 電壓偏差波動率比較Fig.7 Comparison of voltage deviation volatility rate curve

表2 電壓偏差波動率分析Table 2 Analysis of voltage deviation volatility

由圖6及表1可知，在減少網(wǎng)損方面，本文方法在當(dāng)天的降損率曲線幅度要大于其他3種方法的曲線，而降損率均值和方差分別為7.94%和1.67%，均高于其他2種方法。由圖7及表2可知，在減小電壓波動方面，該方法的整體電壓偏差波動率曲線幅度要小于其他3種方法的曲線，而電壓偏差波動率均值和方差分別為0.94%和0.29%,均小于其他2種方法。

5 結(jié)論

針對配電網(wǎng)的無功優(yōu)化問題,提出了基于PSO-DBN的配電網(wǎng)無功優(yōu)化方法，并在IEEE33節(jié)點系統(tǒng)進行仿真實驗，得出以下結(jié)論。

(1)利用隨機矩陣?yán)碚搶ε潆娋W(wǎng)運行狀態(tài)數(shù)據(jù)進行特征提取，大大提高了數(shù)據(jù)利用率，可有效地體現(xiàn)出配電網(wǎng)的運行狀態(tài)，為后續(xù)配電網(wǎng)無功優(yōu)化策略的分析提供了重要數(shù)據(jù)來源。

(2)所建立的改進PSO-DBN模型可搭建配電網(wǎng)運行狀態(tài)與配電網(wǎng)無功優(yōu)化策略之間的關(guān)系，使得無功優(yōu)化效果得以提高。

(3)本文方法對配電網(wǎng)模型參數(shù)依賴程度較低，并不需要配電網(wǎng)的網(wǎng)架結(jié)構(gòu)的具體參數(shù)，僅需將配電網(wǎng)運行狀態(tài)合無功優(yōu)化策略之間建立深度學(xué)習(xí)模型，相較于傳統(tǒng)優(yōu)化方法無功優(yōu)化效果有所提高。