郭永鑫，肖洪光，楊 璐

(國(guó)網(wǎng)吉林省電力有限公司，吉林 長(zhǎng)春 130062)

0 引言

微電網(wǎng)作為充分發(fā)揮分布式能源效用的主要形式，能夠滿(mǎn)足當(dāng)前能源需求的快速提升。其中，以變壓器為核心的微電網(wǎng)結(jié)構(gòu)成為當(dāng)前研究領(lǐng)域的關(guān)注熱點(diǎn)[1]。微電網(wǎng)變壓器受干擾后易產(chǎn)生電壓波動(dòng)問(wèn)題，嚴(yán)重影響微電網(wǎng)供電質(zhì)量，直接或間接導(dǎo)致的經(jīng)濟(jì)損失越來(lái)越嚴(yán)重[2]。抑制變壓器電壓波動(dòng)需以準(zhǔn)確的電壓波動(dòng)估計(jì)為基礎(chǔ)[3]。

作為機(jī)器學(xué)習(xí)的主要組成結(jié)構(gòu)，深度學(xué)習(xí)網(wǎng)絡(luò)通過(guò)學(xué)習(xí)樣本內(nèi)在規(guī)律完成信息識(shí)別[4]。深度學(xué)習(xí)網(wǎng)絡(luò)中的置信網(wǎng)絡(luò)具有特征提取與分類(lèi)功能[5]，可應(yīng)用于變壓器電壓波動(dòng)估計(jì)。為此，提出了基于深度學(xué)習(xí)網(wǎng)絡(luò)的微電網(wǎng)變壓器電壓波動(dòng)估計(jì)與抑制方法，采用深度學(xué)習(xí)網(wǎng)絡(luò)完成特征提取功能，以自動(dòng)生成特征替代人工設(shè)置特征，準(zhǔn)確估計(jì)電壓波動(dòng)情況；在電壓出現(xiàn)波動(dòng)條件下，通過(guò)組合抑制方法抑制其波動(dòng)，保障微電網(wǎng)的穩(wěn)定運(yùn)行。

1 微電網(wǎng)變壓器電壓波動(dòng)抑制方法

1.1 微電網(wǎng)系統(tǒng)結(jié)構(gòu)

微電網(wǎng)系統(tǒng)結(jié)構(gòu)包括主電網(wǎng)、變壓器、光伏陣列、風(fēng)機(jī)、蓄能系統(tǒng)等[6]。變壓器包含三個(gè)層級(jí)[7]，分別是H橋型交流電源(alternating current，AC)與直流電源(direct current，DC)變流器級(jí)聯(lián)構(gòu)建的輸入級(jí)、雙向全橋DC-DC變換器變流器并聯(lián)構(gòu)建的隔離級(jí)、完成電壓源型逆變的輸出級(jí)。三個(gè)層級(jí)分別具有連接主電網(wǎng)不同電壓等級(jí)、變換電壓、提供低壓交流接口引出微電網(wǎng)交流母線的功能。微電網(wǎng)交流母線利用DC-AC、AC-DC-AC變換器連接風(fēng)機(jī)與光伏等分布式能源。具有輔助本地供電功能的蓄能系統(tǒng)利用DC-AC變換器與微電網(wǎng)連接。微電網(wǎng)系統(tǒng)結(jié)構(gòu)框架如圖1所示。

圖1 微電網(wǎng)系統(tǒng)結(jié)構(gòu)框架

變壓器受干擾后易產(chǎn)生電壓波動(dòng)問(wèn)題。引起微電網(wǎng)變壓器電壓波動(dòng)的主要因素如表1所示?；谖㈦娋W(wǎng)變壓器電壓波動(dòng)主要因素造成的電壓波動(dòng)波形，采用深度學(xué)習(xí)網(wǎng)絡(luò)中的置信網(wǎng)絡(luò)構(gòu)建模型，提取波形特征，根據(jù)特征估計(jì)電壓電壓波動(dòng)情況。

表1 微電網(wǎng)變壓器電壓波動(dòng)的主要因素

1.2 變壓器電壓波動(dòng)估計(jì)

1.2.1 深度學(xué)習(xí)網(wǎng)絡(luò)結(jié)構(gòu)

作為具有多隱含層的神經(jīng)網(wǎng)絡(luò)，深度學(xué)習(xí)網(wǎng)絡(luò)中反向傳播(back propagation，BP)神經(jīng)網(wǎng)絡(luò)和受限玻爾茲曼機(jī)(restricted boltzmann machine，RBM)的數(shù)量分別為一層和多層。通過(guò)逐層貪婪學(xué)習(xí)算法調(diào)整不同層級(jí)之間的連接權(quán)重。無(wú)監(jiān)督學(xué)習(xí)性能極強(qiáng)，能夠最大限度地防止類(lèi)別劃分性能較差以及過(guò)度擬合缺陷的出現(xiàn)[8]。圖2所示為深度學(xué)習(xí)網(wǎng)絡(luò)結(jié)構(gòu)。

圖2 深度學(xué)習(xí)網(wǎng)絡(luò)結(jié)構(gòu)圖

圖2中：x1、x2、xn0、f1、f2、fN和l1、l2、lnl分別為可見(jiàn)層輸入情況、相鄰單元層間連接權(quán)重向量和標(biāo)簽層不同單元情況；yx,z和ftp分別為第x個(gè)隱含層第z個(gè)單元情況和BP神經(jīng)網(wǎng)絡(luò)權(quán)重向量。

作為基于能量函數(shù)的概率建模方法，受限玻爾茲曼機(jī)由可視層和隱含層共同組成[9]。其主要功能為利用棧式逐層堆疊學(xué)習(xí)初始數(shù)據(jù)內(nèi)復(fù)雜化準(zhǔn)則。深度學(xué)習(xí)網(wǎng)絡(luò)訓(xùn)練由三個(gè)環(huán)節(jié)組成，分別是模型初始化環(huán)節(jié)、正向無(wú)監(jiān)督學(xué)習(xí)環(huán)節(jié)和反向有監(jiān)督學(xué)習(xí)環(huán)節(jié)。

1.2.2 深度學(xué)習(xí)網(wǎng)絡(luò)計(jì)算原理

設(shè)(v,y)為給定狀態(tài)，通過(guò)式(1)能夠描述深度學(xué)習(xí)網(wǎng)絡(luò)可見(jiàn)層和隱含層單元的統(tǒng)一設(shè)定能量函數(shù)：

(1)

式中：λ=(f,ε,φ)為深度學(xué)習(xí)網(wǎng)絡(luò)參數(shù)；fij為第i個(gè)可見(jiàn)層單元與第j個(gè)隱含層單元之間的權(quán)重，只有鄰近兩個(gè)層間具有連權(quán)重接；εi為第i個(gè)可見(jiàn)層單元；φj為第j個(gè)隱含層單元的偏置；vi、yi均為狀態(tài)參數(shù)；n為可見(jiàn)層單元數(shù)量；m為隱含層單元數(shù)量。

深度學(xué)習(xí)網(wǎng)絡(luò)狀態(tài)概率具有正則分布特性[10]，因此可利用式(2)描述隨機(jī)狀態(tài)的統(tǒng)一概率分布：

(2)

(3)

(4)

式中：s為訓(xùn)練樣數(shù)量。

利用深度學(xué)習(xí)網(wǎng)絡(luò)對(duì)數(shù)似然函數(shù)L(λ)最大化處理，獲取深度學(xué)習(xí)網(wǎng)絡(luò)參數(shù)λ，同時(shí)擬合訓(xùn)練樣本。此時(shí)，可視層輸入數(shù)據(jù)的特征為隱藏層：

(5)

將比散度算法應(yīng)用于對(duì)數(shù)似然函數(shù)中[11]，求解λ。以λ內(nèi)包含的fij為例：

(6)

式中：Dvi,yj,1為訓(xùn)練樣本參數(shù)；Dvi,yj,2為深度學(xué)習(xí)網(wǎng)絡(luò)定義的期望參數(shù)。兩者均可通過(guò)吉布斯采樣獲取。

在此基礎(chǔ)上，可優(yōu)化fij，得到：

fij=βfij+α(Dvi,yj,1-Dvi,yj,2)

(7)

式中：β為動(dòng)量參數(shù)；α為學(xué)習(xí)率參數(shù)。

εi與φj同樣可采用上述過(guò)程實(shí)施優(yōu)化。

1.2.3 基于深度學(xué)習(xí)網(wǎng)絡(luò)的電壓波動(dòng)估計(jì)

深度學(xué)習(xí)網(wǎng)絡(luò)中，不同隱含層特征向量維度與其單元數(shù)量相同。隱含層數(shù)量越多，模型提取的特征精度越高，電壓波動(dòng)情況估計(jì)越準(zhǔn)確。

電壓波動(dòng)估計(jì)的實(shí)現(xiàn)流程是通過(guò)采集微電網(wǎng)變壓器電壓波動(dòng)波形時(shí)域數(shù)據(jù)，經(jīng)過(guò)數(shù)據(jù)預(yù)處理后建立批量輸入深度置信網(wǎng)絡(luò)模型。結(jié)合RBM正向無(wú)監(jiān)督學(xué)習(xí)數(shù)據(jù)各層的特征，采用逐層貪婪算法修正各層單元的偏置值和隱層權(quán)重，并對(duì)順序輸入的不同電壓波動(dòng)波形進(jìn)行準(zhǔn)確分類(lèi)和標(biāo)志。利用BP神經(jīng)網(wǎng)絡(luò)對(duì)訓(xùn)練樣本進(jìn)行反向微調(diào)，完成模型訓(xùn)練。根據(jù)輸入電網(wǎng)的數(shù)據(jù)，可以準(zhǔn)確地提取出微電網(wǎng)變壓器的電壓波動(dòng)特征。

1.3 電壓波動(dòng)抑制

采用深度學(xué)習(xí)網(wǎng)絡(luò)提取電壓波動(dòng)波形特征完成電壓波動(dòng)估計(jì)。若電壓出現(xiàn)波動(dòng)，則需對(duì)其進(jìn)行抑制，保障微電網(wǎng)穩(wěn)定運(yùn)行。在采用比例諧振控制器作為抑制變壓器電壓波動(dòng)的基礎(chǔ)上[12]，考慮到變壓器的可用補(bǔ)償容量，可得比例諧振控制器控制下變壓器輸出電壓為：

(8)

式中：EU(γ)為電壓控制器傳遞函數(shù)；EI(γ)為電流內(nèi)環(huán)比例增益函數(shù)；K為調(diào)制增益參數(shù)；U2為輸出級(jí)參考電壓參數(shù)；i1為輸出電流參數(shù)。

根據(jù)式(8)可知，利用電壓控制器抑制電壓波動(dòng)過(guò)程中，輸出電流i1正向增益EIR(γ)的變化將導(dǎo)致變壓器等效阻抗發(fā)生變化。

(9)

U1(γ)=EU(γ)U2-x(γ)i2-(1-O)x(γ)i2

(10)

式中：O為補(bǔ)償系數(shù)參數(shù)；i2為諧波分量參數(shù)；輸出層基波等效輸出阻抗參數(shù)。

O值能夠?qū)崟r(shí)抑制變壓器電壓波動(dòng)，其取值范圍在0～1之間。其值為0和1時(shí)，分別為變壓器當(dāng)前為未補(bǔ)償和全補(bǔ)償狀態(tài)。O值與諧波阻抗呈反比關(guān)系。O值的提升可使電壓器電壓波動(dòng)得到抑制，其值可根據(jù)變壓器可補(bǔ)償容量和實(shí)際諧波輸出容量得到。

2 實(shí)例驗(yàn)證

以某微電網(wǎng)為研究對(duì)象，在監(jiān)測(cè)系統(tǒng)中采集1 000組電壓數(shù)據(jù)。以800組實(shí)測(cè)變壓器電壓波動(dòng)數(shù)據(jù)作為訓(xùn)練樣本，微電網(wǎng)功率缺額和變壓器投切數(shù)據(jù)分別為180組和170組，短路故障數(shù)據(jù)為450組(單相接地、兩相短路和三項(xiàng)短路分別為150組)。剩余200組作為測(cè)試樣本，不同類(lèi)別數(shù)據(jù)數(shù)量一致。

2.1 微電網(wǎng)變壓器電壓波動(dòng)估計(jì)測(cè)試

采用本文方法進(jìn)行研究對(duì)象變壓器電壓波動(dòng)估計(jì)，構(gòu)建深度學(xué)習(xí)網(wǎng)絡(luò)。模型輸入層單元數(shù)量同變壓器電壓波動(dòng)監(jiān)測(cè)節(jié)點(diǎn)數(shù)量一致。隱含層共6層，各層內(nèi)均包含8個(gè)單元，第六層隱含層單元數(shù)量與提取特征數(shù)量一致。輸出層單元數(shù)量同電壓波動(dòng)數(shù)據(jù)類(lèi)型一致。由此確定，試驗(yàn)過(guò)程中本文方法所用深度學(xué)習(xí)網(wǎng)絡(luò)共包含8層，總節(jié)點(diǎn)數(shù)量為1個(gè)輸入層×180+6個(gè)隱含層×8+1個(gè)輸出層×5=233個(gè)。設(shè)定模型參數(shù)：受限玻爾茲曼機(jī)迭代次數(shù)上限、學(xué)習(xí)率、動(dòng)量參數(shù)分別為60、0.2和1。在深度學(xué)習(xí)網(wǎng)絡(luò)內(nèi)輸入預(yù)處理后的訓(xùn)練樣本數(shù)據(jù)。由于模型第六層隱含層單元數(shù)量為8，因此模型提取特征數(shù)量為8個(gè)。特征提取結(jié)果如表2所示。由表2可知，相同類(lèi)型特征數(shù)據(jù)一致度高，而不同類(lèi)型間特征數(shù)據(jù)差異顯著。這說(shuō)明本文方法能夠有效提取電壓波形特征，且特征顯性突出，有利于提升電壓波動(dòng)估計(jì)精度。

表2 特征提取結(jié)果

根據(jù)提取的特征對(duì)測(cè)試樣本進(jìn)行估計(jì)，結(jié)果如表3所示。由表3可知，本文方法估計(jì)精度達(dá)到98.5%，充分說(shuō)明本文方法能夠準(zhǔn)確進(jìn)行電壓波動(dòng)估計(jì)。

表3 測(cè)試樣本估計(jì)結(jié)果

2.2 微電網(wǎng)變壓器電壓波動(dòng)抑制測(cè)試

在估計(jì)正確的測(cè)試樣本中隨機(jī)選取一組樣本，采用本文方法進(jìn)行波動(dòng)抑制。設(shè)定補(bǔ)償系數(shù)值為0.8，壓波動(dòng)抑制結(jié)果如圖3所示。

分析圖3(a)可知，初始電壓波形中由于電壓波動(dòng)導(dǎo)致波形畸變顯著，電壓波形整體畸變率達(dá)到18.8%。

分析圖3(b)可知，采用本文方法抑制電壓波動(dòng)后，電壓波形與初始電壓波形相比更為平滑整體，波形畸變率降至2.9%，下降幅度為15.9%。由此說(shuō)明本文方法能夠顯著抑制電壓波動(dòng)，保障微電網(wǎng)穩(wěn)定運(yùn)行。

圖3 電壓波動(dòng)抑制結(jié)果

2.3 補(bǔ)償系數(shù)設(shè)定分析

本文方法在不同補(bǔ)償系數(shù)設(shè)定值條件下，初始電壓波形畸變率下降幅度和電壓波動(dòng)抑制過(guò)程中的能耗，補(bǔ)償系數(shù)設(shè)定如表4所示。由表4可知，隨著補(bǔ)償系數(shù)由0.2上升至0.9，畸變率下降幅度和能耗均隨之表現(xiàn)出明顯的上升趨勢(shì)。當(dāng)補(bǔ)償系數(shù)為0.8時(shí)，畸變率下降幅度達(dá)到15.9%，能耗為4.76 J；當(dāng)補(bǔ)償系數(shù)升至0.9時(shí)，畸變率下降幅度雖然上升0.7%，但能耗卻提升2.68 J，畸變率下降幅度與能耗提升不匹配。因此，補(bǔ)償系數(shù)確定為0.8。

表4 補(bǔ)償系數(shù)設(shè)定

3 結(jié)論

本文研究深度學(xué)習(xí)網(wǎng)絡(luò)的微電網(wǎng)變壓器電壓波動(dòng)估計(jì)與抑制方法，利用深度學(xué)習(xí)網(wǎng)絡(luò)中的置信網(wǎng)絡(luò)提取電壓波動(dòng)特征估計(jì)電壓波動(dòng)情況，并提出一種組合抑制方法抑制電壓波動(dòng)。試驗(yàn)結(jié)果顯示，本文方法能夠準(zhǔn)確估計(jì)電壓波動(dòng)情況，并顯著改善電壓波形畸變率，具有較高應(yīng)用價(jià)值。