劉家寧,葉 鵬,葉 臻,程緒可,孫 峰
(1.國網(wǎng)吉林供電公司,吉林 吉林 132000;2.沈陽工程學(xué)院電力學(xué)院,遼寧 沈陽 110136;3.國網(wǎng)遼寧省電力有限公司電力科學(xué)研究院,遼寧 沈陽 110006)
三相負(fù)荷不平衡是配電網(wǎng)中常見的現(xiàn)象,其主要負(fù)面影響之一是導(dǎo)致附加損耗,在計算線損過程中,這部分附加損耗是不容忽視的[1-4]。通過線損的高低,能夠了解供電企業(yè)的精益化管理水平,以及降低線損措施的成效。本文結(jié)合三相不平衡,兼顧最大負(fù)荷電流和線損增量之間的關(guān)系,經(jīng)過研究和分析,確定三相電流不平衡度與線損增量的量化關(guān)系。考慮到試驗數(shù)據(jù)樣本較少,利用對抗生成網(wǎng)絡(luò)予以解決,在這一過程中驗證該網(wǎng)絡(luò)是充分可靠的,從定量的層面揭示出“三相不平衡度-線損增量系數(shù)”二者是如何對彼此造成影響,具有一定的現(xiàn)實價值[5-7]。
在配電網(wǎng)三相負(fù)荷不平衡狀態(tài)下,考慮到各相的負(fù)荷電流存在差異,此時會形成不同相之間的不平衡電流。該電流的存在,直接導(dǎo)致相線和中性線產(chǎn)生了本可以避免的損耗,總線損因此而提高。
對于配電系統(tǒng)而言,如何界定三相電流不平衡度,學(xué)術(shù)界并未達(dá)成一致,可以通過不同的方法計算出來,大部分學(xué)者都以三相平均電流作為指標(biāo),并未有效地兼顧電流最大最小值和不平衡度之間的關(guān)系。由于線損增量對三相電流的最大偏移幅值呈指數(shù)型增長關(guān)系,本文考慮三相電流中最大電流值對不平衡度的影響,則定義三相電流不平衡度:
(1)
Icdh=max(IAh,IBh,ICh)
(2)
(3)
(4)
式中:IAh、IBh、ICh(h=1,2,…,24)為24 h整點三相電流;Icdh(h=1,2,…,24) 為24 h整點三相電流的最大值;Ipjh為24 h整點三相電流的平均值;δh為24 h整點負(fù)荷電流不平衡度,δh∈[0,1];δ為24 h負(fù)荷電流的平均不平衡度。
三相不平衡線損增量系數(shù)計算定義為
(5)
(6)
(7)
式中:Ipj為24個整點三相電流平均值的均值;ΔPph為三相負(fù)荷平衡時的線路功率損耗;ΔPn為第n種情況下求出的線路功率損耗;Kδn為第n種情況下(n=1,2,3,…)的線損增量系數(shù)。
三相不平衡導(dǎo)致的線損增加與每相上的負(fù)荷息息相關(guān),而負(fù)荷的大小又體現(xiàn)在三相電流大小上。所以,現(xiàn)實中出現(xiàn)三相負(fù)荷不平衡的情況有3種,每種情況下的三相負(fù)荷不平衡線損增量系數(shù)是不同的,可以推導(dǎo)出來。
a.三相負(fù)荷中“一相為大負(fù)荷,一相為小負(fù)荷,一相為平均負(fù)荷”。
b.三相負(fù)荷中“一相為大負(fù)荷,兩相為小負(fù)荷”。
c.三相負(fù)荷中“一相為小負(fù)荷,兩相為大負(fù)荷”。
假設(shè)三相電路中每相的電阻值都為R,中性線的橫截面積一般為輸電線路的一半,因此電阻值為2R,線損計算公式為ΔP=I2r(r為所在線路上的電阻值)。
(8)
當(dāng)三相負(fù)荷平衡時線路損耗見式(6)。
根據(jù)兩者的結(jié)果進(jìn)行比較得到:
(9)
式中:Kδ1為第一種情況下線損的增量系數(shù),其含義是該情況下的線損值為三相負(fù)荷平衡狀態(tài)下的Kδ1倍。按照規(guī)程的相關(guān)內(nèi)容可知:在低壓主干線和主要分支的首段,三相負(fù)荷電流不平衡最高為20%。如果有δ=0.2,Kδ1=1.11,這意味著該情況導(dǎo)致線損提高了11%。
b.在上述第二種情況下,用(1+δ)Ipj表示大負(fù)荷相電流,則兩相小負(fù)荷相電流是(1-0.5δ)Ipj,中性線電流是1.5δIpj,此時可以計算出功率損耗,具體如下:
(10)
(11)
此種不平衡類型下,三相負(fù)荷不平衡時其線損值是三相負(fù)荷平衡時線損值的Kδ2倍。δ=0.2時,Kδ2=1.11,即三相負(fù)荷不平衡所引起的線損增加8%。
c.在上述第三種情況下。用(1+δ)Ipj表示兩相大負(fù)荷相電流,則小負(fù)荷相電流是(1-2δ)Ipj,中性線電流是3δIpj,此時可以計算出功率損耗,具體如下:
(12)
根據(jù)比較得到:
(13)
此種不平衡類型下,三相負(fù)荷不平衡時其線損值是三相負(fù)荷平衡時線損值的Kδ3倍。當(dāng)δ=0.2時,Kδ3=1.32,即三相負(fù)荷不平衡所引起的線損增加32%。
由于存在一定的假設(shè)條件,因此計算過程中需要足夠的樣本數(shù)據(jù)保證計算精度。而實際操作過程中由于數(shù)據(jù)的獲取較為困難,因此運(yùn)用對抗生成網(wǎng)絡(luò)增加樣本數(shù)據(jù)保證計算精度。
GAN(對抗生成網(wǎng)絡(luò))具有深度學(xué)習(xí)的功能,是過去幾年間最具潛力、最受關(guān)注的無監(jiān)督學(xué)習(xí)模型,2014年由美國學(xué)者Goodfellow lan提出。GAN是生成模型的一種,GAN的最終目標(biāo)是對真實數(shù)據(jù)的分布或者密度進(jìn)行預(yù)估,并能夠根據(jù)學(xué)到的知識生成新的數(shù)據(jù)。
GAN由2個網(wǎng)絡(luò)構(gòu)成:生成網(wǎng)絡(luò)和判別網(wǎng)絡(luò),其具體結(jié)構(gòu)見圖1。
圖1 對抗生成網(wǎng)絡(luò)結(jié)構(gòu)
生成網(wǎng)絡(luò)G:基于訓(xùn)練集的數(shù)據(jù),對其概率分布持續(xù)進(jìn)行學(xué)習(xí),從而把提供給該網(wǎng)絡(luò)的隨機(jī)噪音數(shù)據(jù)轉(zhuǎn)變成足以以假亂真的數(shù)據(jù)。圖1表明,生成網(wǎng)絡(luò)接收到噪音數(shù)據(jù),然后輸出數(shù)據(jù)給判別網(wǎng)絡(luò)D。
判別網(wǎng)絡(luò)D:能夠識別數(shù)據(jù)的真實性,在整個對抗網(wǎng)絡(luò)中,其作用是把生成網(wǎng)絡(luò)G形成的“假”數(shù)據(jù)和來自訓(xùn)練集的“真”數(shù)據(jù)區(qū)分開來。如圖1中所示,判別網(wǎng)絡(luò)D接收生成網(wǎng)絡(luò)G送過來的假數(shù)據(jù)進(jìn)行一個二分類的判定。當(dāng)生成網(wǎng)絡(luò)G產(chǎn)生的數(shù)據(jù)足夠逼真時,判別網(wǎng)絡(luò)就判定為真,此時生成網(wǎng)絡(luò)就可以制造試驗數(shù)據(jù)。
生成網(wǎng)絡(luò)與判別網(wǎng)絡(luò)是完全獨立的2個模型,訓(xùn)練這2個模型的方法就是單獨交替迭代優(yōu)化訓(xùn)練,整個訓(xùn)練過程通過相互競爭讓這2個網(wǎng)絡(luò)同時得到增強(qiáng),G和D構(gòu)成一個動態(tài)的“博弈過程”。對于對抗生成網(wǎng)絡(luò)來說,其最終優(yōu)化目標(biāo)函數(shù):
Ez~Pnoise(z)[log(1-D(G(Z)))]
(14)
式中:E(*)為分布函數(shù)的期望值;Pdate(x)為真實樣本的分布;Pnoise(z)為定義在低維的噪聲分布。
此目標(biāo)函數(shù)是最大和最小優(yōu)化,因此是不可能一次完成的,先對D進(jìn)行優(yōu)化,再對G進(jìn)行優(yōu)化,這實際上是2個優(yōu)化問題。將目標(biāo)函數(shù)拆解可得:
判別網(wǎng)絡(luò)D的優(yōu)化:
+Ez~Pnoise(z)[log(1-D(G(Z)))]
(15)
生成網(wǎng)絡(luò)G的優(yōu)化:
(16)
在訓(xùn)練過程中固定一方,更新另一方網(wǎng)絡(luò)參數(shù),交替迭代,使另一方的誤差最大,最后,G可以估計樣本數(shù)據(jù)的分布。生成網(wǎng)絡(luò)G隱式定義了概率分布Pnoise,實際目的是為了Pnoise收斂到數(shù)據(jù)Pdate的真實分布。這個交替的參數(shù)優(yōu)化過程可以用馬鞍圖來表示更加直觀,如圖2所示。
(a)網(wǎng)絡(luò)D的訓(xùn)練過程
(b)網(wǎng)絡(luò)G的訓(xùn)練過程圖2 GAN網(wǎng)絡(luò)訓(xùn)練過程
先確定生成網(wǎng)絡(luò)G,最大化判別網(wǎng)絡(luò)D,讓點沿著判別網(wǎng)絡(luò)D變大的方向移動(藍(lán)色箭頭),然后確定判別網(wǎng)絡(luò)D,最小化生成網(wǎng)絡(luò)G,讓點沿著生成網(wǎng)絡(luò)G變小的方向移動(紅色箭頭)。循環(huán)上述若干步后,達(dá)到期望的鞍點(理想最優(yōu)解)。
基于對抗生成網(wǎng)絡(luò)的三相不平衡線損量化分析首先采集多個臺區(qū)多天的整點時刻的三相四線電流值,考慮三相最大電流對不平衡度的影響,運(yùn)用式(4)得到日平均三相電流不平衡度。由于不平衡度由三相電流計算得到,依據(jù)三相電流大小將三相不平衡情況分成3類。
基于對抗生成網(wǎng)絡(luò)的三相不平衡線損量化分析主要分為以下4個步驟:
a.依據(jù)三相電流對三相不平衡情況進(jìn)行分類;
b.運(yùn)用對抗生成網(wǎng)絡(luò)增加每一類的試驗數(shù)據(jù);
c.利用補(bǔ)充后的試驗數(shù)據(jù)計算三相不平衡線損增量系數(shù);
d.構(gòu)建“三相不平衡度-線損增量系數(shù)”的量化關(guān)系。
具體計算流程如圖3所示。
圖3 “三相不平衡度-線損增量系數(shù)”量化關(guān)系計算流程
由某市所屬各分區(qū)供電公司提供的共計1150個10/0.4 kV變壓器低壓側(cè)1天24個整點時刻內(nèi)采集的整點電流數(shù)據(jù)。
臺區(qū)屬性為A和B,A類臺區(qū)分布于經(jīng)濟(jì)發(fā)達(dá),人口集中地區(qū)的城網(wǎng),負(fù)荷密度每km2大于10 MW;B類臺區(qū)分布于城市的建成區(qū)或規(guī)劃區(qū),包括經(jīng)濟(jì)開發(fā)區(qū)和工業(yè)園區(qū)所在地,負(fù)荷密度在5~10 MW之間。由于城市存在著大量時空分布不平衡的單相負(fù)荷,而工業(yè)園區(qū)中又存在著較多的沖擊性負(fù)荷,因此這2類臺區(qū)存在著較為嚴(yán)重的三相不平衡想象。
由于有些樣本臺區(qū)位于老舊城區(qū),新建開發(fā)區(qū)以及工業(yè)園區(qū)的原因,這些地方的智能電能表存在著故障損壞沒有及時更換以及還未能及時配套安裝等問題。首先對所得到的臺區(qū)三相電流數(shù)據(jù)進(jìn)行預(yù)處理,要先刪除一些空白數(shù)據(jù),去除明顯不符合業(yè)務(wù)實際的數(shù)據(jù),并將所剩數(shù)據(jù)整理在一張表中,共計1125個臺區(qū)。另外,根據(jù)式(1)和式(5)計算日平均電流,將三相電流與日平均電流比較,將三相不平衡的情況分成3類。如表1所示。
表1 三相不平衡臺區(qū)分類結(jié)果
抽取每一個不平衡度下各50個臺區(qū)三相平衡時的真實線損率,計算各臺區(qū)三相不平衡時的線損率,并與其真實的三相不平衡線損率進(jìn)行比較,將相對誤差在0.1%以內(nèi)設(shè)為準(zhǔn)確。得到使用對抗生成網(wǎng)絡(luò)前后的相對誤差折現(xiàn)圖如圖4所示。
(a)使用GAN相對誤差折線圖
(b)未使用GAN相對誤差折線圖圖4 相對誤差折線圖表
對抗生成網(wǎng)絡(luò)GAN使用前后準(zhǔn)確率對比如表2所示。
表2 使用GAN前后計算的準(zhǔn)確率對比 單位:%
由圖4可以發(fā)現(xiàn),相對于使用對抗生成網(wǎng)絡(luò)以后得到的相對誤差,未使用對抗生成網(wǎng)絡(luò)的相對誤差更為分散,計算誤差更大。因此,依據(jù)表2可以得出結(jié)論,若經(jīng)過對抗生成網(wǎng)絡(luò)有效增加樣本數(shù)據(jù)以后,再進(jìn)行線損增量系數(shù)的計算,可以提高計算的精確性。
本文提出了一種基于對抗生成網(wǎng)絡(luò)的三相不平衡線損量化分析方法。該方法主要利用了當(dāng)下比較流行的一種深度學(xué)習(xí)算法對抗生成網(wǎng)絡(luò)GAN,并將其運(yùn)用到電網(wǎng)中的三相不平衡線損量化計算中去。在本文中,考慮了三相電流的最大值對不平衡度的影響。對某市提供的1150條包含完整參數(shù)的供電線路進(jìn)行算例分析。通過使用對抗生成網(wǎng)絡(luò)增加樣本數(shù)據(jù)前后的平行試驗對比發(fā)現(xiàn),本文利用基于對抗生成網(wǎng)絡(luò)的線損量化分析方法可以有效解決原始數(shù)據(jù)偏少的問題,大大提高三相不平衡線損增量系數(shù)的計算精度,構(gòu)建起“三相不平衡度-線損增量系數(shù)”的量化關(guān)系,具有現(xiàn)實的指導(dǎo)意義。這也為其他試驗計算過程中,面臨原始數(shù)據(jù)偏少的問題提供了一種解決思路。