粒子群優(yōu)化模擬電荷法在HGIS工頻電場計(jì)算中的應(yīng)用

李永明，汪小莞，楊勃，于金剛，徐祿文

(1. 重慶大學(xué) 輸配電裝備及系統(tǒng)安全與新技術(shù)國家重點(diǎn)實(shí)驗(yàn)室，重慶 400044；2. 國網(wǎng)重慶市電力公司電力科學(xué)研究院，重慶 401123)

0 引 言

我國疆域廣大能源中心和負(fù)荷中心相距遙遠(yuǎn)，為了提高輸電能力，實(shí)現(xiàn)大功率的中、遠(yuǎn)距離輸電，超高壓變電站大量增加[1]。為了緩解環(huán)境和占地的壓力，超高壓500 kV變電站的主設(shè)備選型尤為重要。氣體絕緣半封閉高壓開關(guān)(Half Gas Insulated Switch，HGIS)兼顧了敞開式和封閉式開關(guān)設(shè)備的優(yōu)點(diǎn)，具有極高的運(yùn)行穩(wěn)定性和抗干擾能力，且占地面積小，維護(hù)方便，完全滿足目前電網(wǎng)的建設(shè)需求[2]。

由于HGIS采用常規(guī)的空氣絕緣懸吊方式架設(shè)，同時(shí)為方便進(jìn)出線和懸吊母線，架設(shè)了大量金屬構(gòu)架。這些因素導(dǎo)致HGIS設(shè)備區(qū)出現(xiàn)較高的工頻電場，甚至大于一般的空氣絕緣敞開式開關(guān)設(shè)備變電站[3-4]。因此研究分析500 kV HGIS工頻電場的分布規(guī)律，對HGIS變電站的電磁環(huán)境評估和優(yōu)化設(shè)計(jì)等方面具有重要的現(xiàn)實(shí)意義[5]。目前國內(nèi)外多采取現(xiàn)場測量、仿真計(jì)算的方法進(jìn)行研究分析?，F(xiàn)場測量易受外界環(huán)境影響，需要消耗大量的人力物力，且只能檢測有限的點(diǎn)，難以全面地獲取現(xiàn)場工頻電場的整體分布規(guī)律。因此數(shù)值計(jì)算的方法顯得極具優(yōu)越性，且隨著計(jì)算機(jī)技術(shù)的發(fā)展，大規(guī)模的數(shù)值運(yùn)算不再難以實(shí)現(xiàn)，因而受到更多學(xué)者的青睞[6-8]。模擬電荷法作為常用的工頻電場數(shù)值計(jì)算方法，原理簡單，便于編程，具有較高的計(jì)算精度[9-10]。但在合理地確定模擬電荷的類型、位置、數(shù)量的關(guān)鍵技術(shù)上，不可避免地與操作者的直接經(jīng)驗(yàn)相關(guān)，因人而異的直接經(jīng)驗(yàn)使模擬電荷法難以實(shí)現(xiàn)高精度的數(shù)值分析[11-12]。針對這些弊端文章將粒子群優(yōu)化算法與模擬電荷法結(jié)合形成粒子群優(yōu)化模擬電荷法(PSO-CSM)，算法將自動尋求模擬電荷的最佳位置和量值，提高工頻電場數(shù)值計(jì)算的精度和效率。

1 粒子群優(yōu)化模擬電荷法

1.1 模擬電荷法

模擬電荷法基于電磁場的唯一性定理，將電極表面連續(xù)分布的自由電荷或介質(zhì)分界面上連續(xù)分布的束縛電荷，用位于計(jì)算場域邊界之外的虛設(shè)的一組數(shù)量有限、布置在一定幾何位置上、離散化的模擬電荷予以等值替代，并將原來具有邊界的非均勻媒質(zhì)空間變換成無限大的均勻媒質(zhì)空間[13]?；讷@取的模擬電荷量值，再根據(jù)電場疊加原理，可近似求得原連續(xù)分布電荷在空間任意一點(diǎn)處所產(chǎn)生的電場[14-15]。

1.2 粒子群算法

美國的Kennedy和Eberhar等學(xué)者受鳥群覓食行為的啟發(fā)，于1995年提出了粒子群優(yōu)化算法(Particle Swarm Optimization, PSO)[16]。粒子在可行域內(nèi)移動尋找目標(biāo)，粒子的移動方向和距離取決于自身的速度。粒子每完成一次移動，便根據(jù)自身搜索到的“最優(yōu)點(diǎn)”(Pbest)以及所有相鄰粒子到目前為止找到的“最好點(diǎn)”(Gbest)修正當(dāng)前的速度，然后向著距離目標(biāo)最近的方向進(jìn)行下一次移動[17-19]。速度和位置更新公式分別為：

(1)

(2)

式中ω為慣性權(quán)重；d=1, 2, …,D；i=1, 2, …,n；k為當(dāng)前迭代次數(shù)；Vid為粒子速度；c1與c2為加速度因子，為非負(fù)常數(shù)；r1與r2是0～1之間的隨機(jī)數(shù)。PSO屬于群體智能優(yōu)化算法，因此對于多變量、非線性優(yōu)化問題的尋優(yōu)有著較大的優(yōu)勢，極具工程應(yīng)用意義[20-22]。

1.3 粒子群優(yōu)化模擬電荷法的基本思路

傳統(tǒng)的模擬電荷法主要依據(jù)電極外形特點(diǎn)選擇模擬電荷的類型、位置和數(shù)量,使得所設(shè)置的模擬電荷的整體在給定電極表面的總電位有可能滿足電極表面等電位要求。電位方程式為：

[Pm×n][Qn×1]=[φm×1]

(3)

對于形狀不規(guī)則的電極，在設(shè)置模擬電荷的位置和量值時(shí)通常需反復(fù)多次試驗(yàn)且不能保證最優(yōu)。因此文中嘗試將粒子群優(yōu)化算法與模擬電荷法結(jié)合，形成基于粒子群算法的優(yōu)化模擬電荷法。利用粒子群算法快速自動尋優(yōu)的特點(diǎn)在保證精度的前提下使用較少的模擬電荷自動獲取合理的布點(diǎn)方案和模擬電荷量值。

粒子群優(yōu)化模擬電荷法通過粒子群算法產(chǎn)生初始化粒子群，包括位置參量粒子和電荷量參數(shù)粒子。將每個粒子內(nèi)數(shù)值分別傳遞給模擬電荷法中模擬電荷的位置和電荷量。對于每個粒子，由初始模擬電荷位置參數(shù)計(jì)算電位系數(shù)矩陣，并與初始電荷量相乘，獲取該粒子本次對應(yīng)的φm×1值，如式(4)所示：

(4)

對于每個粒子，由所得匹配點(diǎn)電位，計(jì)算平均電位誤差，求得結(jié)果即為適應(yīng)度值并記憶。并執(zhí)行設(shè)定的誤差指標(biāo)計(jì)算，將誤差指標(biāo)值回傳給粒子群算法進(jìn)行計(jì)算。然后對比個體極值與所有走過的路徑的適應(yīng)度值。若個體極值最小則保留，若走過的路徑里存在小于目前個體極值的，則將該路徑賦值給個體極值。

對于全部粒子，對比所有適應(yīng)度值與群體極值。若群體極值最小則保留，若某個粒子存在更小的適應(yīng)度值，則將該值賦值給群體極值。通過式(1)和式(2)更新該粒子的速度和位置，再次計(jì)算各粒子適應(yīng)度值，并檢查是否達(dá)到終止條件。若未到達(dá)終止條件則更新粒子并繼續(xù)上述循環(huán)，直至達(dá)到終止條件。

粒子群優(yōu)化模擬電荷法的算法流程圖如圖1所示。

圖1 粒子群優(yōu)化模擬電荷法流程圖

2 粒子群優(yōu)化模擬電荷法的驗(yàn)證

分別采用粒子群優(yōu)化模擬電荷法與傳統(tǒng)模擬電荷法計(jì)算球-板電極模型的電場分布，并將計(jì)算結(jié)果進(jìn)行對比，從而驗(yàn)證粒子群優(yōu)化模擬電荷法的有效性。

2.1 球-板模型

如圖2所示，球狀電極位于地面之上，球心距離地面高度h0=5 cm，球狀電極半徑為r0=1 cm。大地為水平面，電位為0 V。球體表面電位為等勢面，φ=1 V。

圖2 球-板模型及其優(yōu)化模擬電荷配置示意圖

2.2 粒子群優(yōu)化模擬電荷法的對比驗(yàn)證

采用粒子群優(yōu)化模擬電荷法計(jì)算球型電極周圍的電場分布。優(yōu)化算法以匹配點(diǎn)平均電位誤差最小為優(yōu)化目標(biāo)，在球體表面設(shè)置m=100個匹配點(diǎn)，約束條件為計(jì)算場域。由此構(gòu)造的優(yōu)化問題為：

(5)

式中φvi為點(diǎn)電荷在第i個匹配點(diǎn)的電位疊加之和。優(yōu)化變量為{x}T=[q1,q2,z1,z2]。

在球體電極表面選擇20個校驗(yàn)點(diǎn)。傳統(tǒng)模擬電荷法與優(yōu)化模擬電荷法校驗(yàn)點(diǎn)誤差對比如圖3所示。

圖3 誤差對比

由圖3可知，采用粒子群優(yōu)化模擬電荷法檢驗(yàn)點(diǎn)誤差更小。傳統(tǒng)模擬電荷法依據(jù)既定的電位系數(shù)矩陣和匹配點(diǎn)電位求解模擬電荷量，損失了對其它點(diǎn)的考慮，難以達(dá)到最優(yōu)配置。而從優(yōu)化角度選取最佳的模擬電荷量值，可均勻地照顧到更多的點(diǎn)從而提高計(jì)算精度。

3 500 kV HGIS配電裝置區(qū)工頻電場計(jì)算分析

3.1 HGIS配電裝置區(qū)模型建立

HGIS配電裝置區(qū)設(shè)備外形復(fù)雜多樣，需對設(shè)備進(jìn)行一定的簡化再建模。HGIS懸吊式管母線可視為一定長度的金屬空心圓管。母線連接線、設(shè)備進(jìn)出線、設(shè)備間相連的短接導(dǎo)線等均可簡化為有限長直線段導(dǎo)體。管體表面和導(dǎo)線表面為等電位面。對于設(shè)備基座，“A”型架構(gòu)、避雷器等設(shè)備，主要根據(jù)其外形特點(diǎn)，在內(nèi)部配置一定數(shù)量的模擬點(diǎn)電荷進(jìn)行等效替代。同時(shí)忽略配電裝置區(qū)內(nèi)絕緣介質(zhì)對區(qū)域內(nèi)工頻電場的影響?；谏鲜龊喕瓌t，建立如圖4所示的HGIS工頻電場計(jì)算基本單元的簡化模型。

圖4中HGIS金屬支座長0.3 m，寬0.2 m，高2.3 m，其他參數(shù)如表1所示。

圖4 HGIS基本單元簡化模型

表1 HGIS基本單元外形參數(shù)

3.2 場域分析

圖5為4回出線HGIS配電裝置區(qū)界面圖。

圖5 500 kV HGIS配電裝置區(qū)截面圖

如圖5所示，每串各相間距為7.5 m，縱向跨距約為27 m。母線直徑200 mm，壁厚10 mm，跨距40 m，懸吊高度20 m，相間距7.5 m。首尾端裝有均壓球，對地電位為500 kV，雙母線之間跨距28 m，以相同相序(ABC-ABC)配置。“A”型架構(gòu)共6個，底部跨距5.6 m，柱體直徑20 cm，垂直高度24 m，架構(gòu)良好接地。設(shè)備進(jìn)出線全部采用2*LGJ500/35型軟導(dǎo)線，出線高度為26 m，相間距為7.5 m，為簡化計(jì)算，進(jìn)出線分裂導(dǎo)線采用等效半徑，走線垂直或水平。

HGIS設(shè)備總縱向跨度為25.8 m，斷路氣室半徑0.5 m，設(shè)備外殼接地。以中間相HGIS設(shè)備中心為z軸建立空間直角坐標(biāo)系，x軸沿母線方向，取值范圍為(-20, 20) m，y取值范圍為(-80, 80) m，計(jì)算場點(diǎn)為地面1.5 m高度處平面。

3.3 HGIS設(shè)備對場域工頻電場的影響

為研究HGIS設(shè)備對配電裝置區(qū)工頻電場產(chǎn)生的影響，建立了“3+0”模式排布的一個完整串HGIS設(shè)備模型，如圖6所示，母線上模擬電荷配置基于優(yōu)化模擬電荷法。

圖6 一個完整串HGIS設(shè)備簡化計(jì)算模型

計(jì)算場域內(nèi)地面上方1.5 m高度處的工頻電場強(qiáng)度，所得結(jié)果如圖7所示。

圖7 HGIS設(shè)備存在時(shí)工頻電場等值線圖

可以看出，整個區(qū)域工頻電場強(qiáng)度不超過4 000 V/m，場強(qiáng)最大處為母線邊相下方。同時(shí)，HGIS減弱了所在區(qū)域附近的工頻電場，說明了HGIS設(shè)備能夠屏蔽掉其下方一部分工頻電場。HGIS設(shè)備周圍有幾處電場強(qiáng)度較高，是因?yàn)樵贖GIS設(shè)備邊角處場強(qiáng)發(fā)生畸變的緣故。

分別計(jì)算HGIS設(shè)備對相間道路地面上方1.5 m高度處和母線下方環(huán)道上方1.5 m高度處的工頻電場強(qiáng)度的影響。計(jì)算結(jié)果如圖8、圖9所示，在相間道路上，HGIS對工頻電場有明顯屏蔽作用，最大衰減了580 V/m。而在環(huán)道上，HGIS的存在僅使環(huán)道處工頻電場略微下降。

圖8 相間道路上工頻電場對比圖

圖9 巡視環(huán)道上工頻電場對比圖

穿過HGIS設(shè)備套管中心距離地面2.8 m高度處的工頻電場等值線圖如圖10所示?？梢钥吹诫S著位置的升高，工頻電場強(qiáng)度整體升高。HGIS設(shè)備對此處工頻電場已無明顯屏蔽作用。而在設(shè)備兩端存在一些場強(qiáng)較高的電場畸變點(diǎn)，最大處約2 400 V/m。對于這些地方應(yīng)給予重視，在設(shè)備表面安置二次設(shè)備時(shí)，應(yīng)避開這些場強(qiáng)容易畸變的部位。

圖10 考慮HGIS設(shè)備影響時(shí)地面2.8 m處工頻電場等值線圖

圖11為有無“A”型架構(gòu)對配電裝置區(qū)巡視環(huán)道及相間道路1.5 m高度處工頻電場的分布情況。

從圖11(a)可以看出“A”型架構(gòu)對環(huán)道中間的工頻電場有較明顯的屏蔽作用，而對環(huán)道兩端無明顯影響。從圖11(b)可以看出“A” 型架構(gòu)(包括其桁梁)的存在使原本雙母線邊相下較高的工頻電場得到明顯削弱，而對于原本場強(qiáng)較弱的HGIS設(shè)備區(qū)相間則無明顯影響。此結(jié)論可為變電站工頻電場屏蔽提供一定參考，例如可通過在激勵源上方加裝屏蔽線以降低下方工頻電場強(qiáng)度。

圖11 “A”型架構(gòu)對環(huán)道和相間道路工頻電場影響

3.4 HGIS配電裝置區(qū)整體工頻電場

HGIS配電裝置區(qū)設(shè)備眾多，對區(qū)域內(nèi)工頻電場的影響各異，因此在進(jìn)行配電裝置區(qū)整體工頻電場的計(jì)算時(shí)，需要將所建主要模型全部考慮在內(nèi)。將HGIS設(shè)備、母線、設(shè)備支座、“A”型構(gòu)架、設(shè)備進(jìn)出線等主要設(shè)備組合在一起，計(jì)算地面1.5 m高度平面工頻電場強(qiáng)度。整體計(jì)算模型如圖12所示。

圖12 HGIS配電裝置區(qū)工頻電場整體計(jì)算模型

采用粒子群優(yōu)化模擬電荷法計(jì)算整個配電裝置區(qū)內(nèi)的工頻電場，計(jì)算結(jié)果如圖13所示。

圖13 500 kV HGIS配電裝置區(qū)電場整體分布

由圖13可知， 500 kV HGIS配電裝置區(qū)工頻電場較高的區(qū)域主要集中在雙母線內(nèi)側(cè)。對于HGIS設(shè)備區(qū)，由于區(qū)域內(nèi)存在較多的設(shè)備進(jìn)出線，且設(shè)備外形較為復(fù)雜，導(dǎo)致個別設(shè)備表面電場發(fā)生畸變，電場強(qiáng)度最大值為9.434 kV/m。

圖14為整體區(qū)域的電場等值線圖，可以看出HGIS兩邊相的兩端，以及“A”型架構(gòu)附近的工頻電場較高。由于各相設(shè)備進(jìn)出線存在相位差，中相HGIS設(shè)備附近場強(qiáng)較低，因此在中相HGIS設(shè)備下方安裝智能組件可有效降低工頻電場對其的干擾。

圖14 HGIS配電裝置區(qū)整體電場等值線圖

4 結(jié)束語

提出將粒子群優(yōu)化算法與模擬電荷法相結(jié)合形成粒子群優(yōu)化模擬電荷法，利用粒子群算法全局尋優(yōu)的特點(diǎn)有效地解決了傳統(tǒng)模擬電荷法布點(diǎn)困難及校驗(yàn)繁瑣的問題。并通過球-板電極算例驗(yàn)證了粒子群優(yōu)化模擬電荷法的有效性。然后利用粒子群優(yōu)化模擬電荷法計(jì)算分析了HGIS配電裝置區(qū)工頻電場的分布情況。通過分析計(jì)算結(jié)果發(fā)現(xiàn)HGIS設(shè)備的金屬外殼對其所在區(qū)域的工頻電場有一定的屏蔽作用。利用這一點(diǎn)，可將必要的變電站二次智能組件安裝于此，以減少工頻電場對敏感設(shè)備的電磁干擾。