劉青,查虹麗,馬龍雄,周寧馨

(西安科技大學(xué)電氣與控制工程學(xué)院,710054,西安)

地磁暴誘發(fā)的地磁擾動會在中性點接地的交流電網(wǎng)中引起地磁感應(yīng)電流(GIC)。準(zhǔn)直流性質(zhì)的GIC可能引起變壓器直流偏磁,對設(shè)備及系統(tǒng)的安全運行產(chǎn)生影響,甚至引發(fā)大面積停電事故[1-2]。中國電網(wǎng)發(fā)展迅速,具有電壓等級高、網(wǎng)絡(luò)規(guī)模大、交直流混聯(lián)等特點,GIC可能會對電網(wǎng)導(dǎo)致嚴(yán)重的致災(zāi)風(fēng)險,因此研究GIC的治理方法尤為重要。

關(guān)于電網(wǎng)GIC的抑制,多位學(xué)者做了大量的研究工作。文獻(xiàn)[3-4]研究了電容隔直裝置對GIC的抑制效果。針對電容隔直裝置風(fēng)險轉(zhuǎn)移引起相鄰變電站GIC增大的問題,文獻(xiàn)[5-6]提出了小電阻抑制方法,對加裝電阻后的GIC計算方法的效果進(jìn)行了分析。文獻(xiàn)[7]指出電阻限流器裝置簡單,對系統(tǒng)其他裝置不會帶來明顯影響,并且系統(tǒng)對電阻裝置的要求只有兩個:電阻阻值一般在幾歐姆左右;電阻對并聯(lián)放電間隙、避雷器進(jìn)行保護(hù)。

GIC抑制措施的安裝位置選擇是全局優(yōu)化問題,最早多采用窮舉法進(jìn)行試驗,比如在高風(fēng)險節(jié)點、系統(tǒng)拐點處重點抑制。但是,對于節(jié)點較多的大電網(wǎng),很難列舉完所有方案,且計算效率較低。文獻(xiàn)[8-9]以電容裝置數(shù)量最少為單一目標(biāo)進(jìn)行智能優(yōu)化,但忽略了其他條件的相互制約。文獻(xiàn)[10-11]用多目標(biāo)粒子群算法分別對GIC影響下無功補償裝置的安裝位置以及直流接地極入地電流作用下電容隔直裝置的安裝位置進(jìn)行了優(yōu)化,能夠一定程度上反映實際系統(tǒng)的優(yōu)化目標(biāo)之間的相互制約性,但二者分別只是與未治理的節(jié)點電壓和未治理的直流量進(jìn)行了比較,并未提出評價指標(biāo)對所提方法進(jìn)行驗證。

另外,以往多以均勻電場作用下的GIC衡量治理方案的效果,如文獻(xiàn)[12]用關(guān)于GIC的4個指標(biāo)評估治理方案的效果,但這些均是基于穩(wěn)態(tài)GIC的指標(biāo),且與目標(biāo)函數(shù)有重復(fù)。為更全面地評價優(yōu)化效果,本文研究時變電場下動態(tài)模型暫態(tài)量的變化規(guī)律,并提出3個指標(biāo)對優(yōu)化方案的效果進(jìn)行評價。

綜上所述,目前關(guān)于GIC的優(yōu)化措施存在電容方案的風(fēng)險轉(zhuǎn)移和高成本、評價指標(biāo)不夠全面的問題。因此,本文針對這些問題,基于多目標(biāo)離散粒子群算法,研究采用電阻限流裝置的GIC抑制措施優(yōu)化方案,并采用時變電場作用下能夠反映暫態(tài)過程的多個指標(biāo)對抑制效果進(jìn)行評價。

1 傳統(tǒng)隔直裝置存在的問題

圖1是地磁感應(yīng)電流Benchmark標(biāo)準(zhǔn)算例的單線圖。Benchmark標(biāo)準(zhǔn)算例由美國電力科學(xué)院牽頭、多國科學(xué)家合作提出,用來測試GIC算法,線路參數(shù)和節(jié)點參數(shù)詳見文獻(xiàn)[13]。

SUB1～SUB6,SUB8—變電站;G1～G7—發(fā)電機;T1～T15—變壓器;Sta7—開關(guān)站;1～8,11～19—母線編號。圖1 GIC-Benchmark標(biāo)準(zhǔn)算例單線圖

計算電網(wǎng)中的GIC主要有兩個步驟[14-15]:根據(jù)變電站的地理坐標(biāo)計算出各變電站之間的等效距離,將感應(yīng)地電場沿線路積分得到驅(qū)動GIC的等效電壓源;建立考慮多電壓等級的電網(wǎng)GIC全節(jié)點計算電路,采用基于電網(wǎng)節(jié)點導(dǎo)納矩陣的電路求解方法求解GIC。

假設(shè)某系統(tǒng)有N個節(jié)點,Ie(N×1階)為流入或流出接地節(jié)點的電流,則

(1)

式中:Ze為接地阻抗矩陣;Ytotal為網(wǎng)絡(luò)導(dǎo)納矩陣之和;V0為感應(yīng)電壓矩陣;Yn為導(dǎo)納矩陣,n為節(jié)點。節(jié)點i和j之間的電流為

(2)

定義Um和Yn都為N階方陣,將式(2)改寫為矩陣形式

In=-(V0+Um)Yn

(3)

In即為所求節(jié)點GIC。

文獻(xiàn)[16]規(guī)定三相三柱式的GIC限值應(yīng)不超過變壓器額定電流的0.7%。文獻(xiàn)[17]指出,對于500 kV變壓器,如果直流作用時間不超過30 min,20 A以下的直流電流將不會對變壓器造成危害。本文計算選取Benchmark標(biāo)準(zhǔn)算例,500 kV變壓器的中性點單相允許直流量為10 A,345 kV變壓器的中性點單相允許直流量為8 A,整個電網(wǎng)所允許的直流限值總和為384 A。GIC-Benchmark各變電站內(nèi)變壓器信息和中性點單相允許直流電流如表1所示。

表1 GIC-Benchmark各變電站內(nèi)變壓器信息和中性點單相允許直流電流

為了說明電容隔直裝置存在的問題,此處采用感應(yīng)地電場強度E=0.5 V/km的均勻電場分析兩種抑制方案。方案1為高風(fēng)險節(jié)點優(yōu)先安裝,根據(jù)GIC計算結(jié)果在GIC較大的SUB6和SUB8安裝電容隔直裝置。方案2為超標(biāo)就安裝,隔直裝置安裝位置為SUB5、SUB6、SUB8。

圖2對比了未治理、方案1和方案2這3種情況下各變電站對地GIC的絕對值,即|IGIC|。圖中圓圈大小表示該站點GIC的大小,虛線框表示GIC被動升高的站點。

(a)未治理

(b)方案1

(c)方案2圖2 3種情況下變電站對地GIC分布

從圖2可以看出,雖然兩種方案都阻斷了安裝隔直裝置變電站的GIC,但卻導(dǎo)致周邊的SUB3的GIC比未采取抑制措施時大幅增加。本文采用中性點串聯(lián)電阻的方案抑制GIC,提出基于多目標(biāo)離散粒子群算法的電阻優(yōu)化方案。

2 電阻措施的優(yōu)化模型及結(jié)果分析

2.1 優(yōu)化模型

本文對變壓器中性點串聯(lián)電阻的方案進(jìn)行優(yōu)化,約束條件為各變電站中性點GIC小于直流限值,優(yōu)化目標(biāo)為限流裝置數(shù)量最少以及全網(wǎng)變壓器中性點GIC總和最小,優(yōu)化模型為

(4)

式中:f1和f2為待優(yōu)化的兩個目標(biāo);xi為電阻裝置的安裝位置,是一串7位的0-1序列,0表示不安裝,1表示安裝,i=1,2,…,7為節(jié)點編號;sum(xi)為xi這7位0-1序列中1的數(shù)量之和,即需要安裝的裝置數(shù)量;Ia和Ib分別為345 kV和500 kV電壓等級下的直流限值。兩個目標(biāo)相互制約,其中一個目標(biāo)達(dá)到最優(yōu)的同時會導(dǎo)致另一個目標(biāo)成為不好的解甚至是最壞的解。因此,需要在滿足兩個目標(biāo)的所有解集中選取不壞的解或者Pareto最優(yōu)解[18]。

對于a、b兩個個體,如果a的所有目標(biāo)都不比b的差,并且a至少有一個目標(biāo)比b的好,那么a支配b,或者a非劣于b。Pareto最優(yōu)解定義為:對于最小化多目標(biāo)問題,n個目標(biāo)分量fi(i=1,…,n)組成的向量為F(x)=[f1(x),f2(x),…,fn(x)],當(dāng)且僅當(dāng)a非劣于b,且不存在另一個決策變量b使F(b)支配F(a)時,決策變量a為Pareto最優(yōu)解。

2.2 多目標(biāo)離散粒子群算法的計算流程

根據(jù)優(yōu)化目標(biāo)和多目標(biāo)之間的Pareto支配關(guān)系,利用多目標(biāo)離散粒子群算法進(jìn)行優(yōu)化計算。

為了讓適應(yīng)度越大的粒子越容易尋優(yōu),將關(guān)于電阻限流裝置數(shù)量和變壓器中性點GIC總和的適應(yīng)度函數(shù)用目標(biāo)函數(shù)的倒數(shù)形式進(jìn)行計算

F1(i)=1/sum(xi)

(5)

(6)

初始篩選非劣解集是指初始化后,某粒子不受其他粒子支配時,將其放入初始非劣解集中,并在更新前從非劣解集中選一個粒子作為群體最優(yōu)粒子;更新非劣解集是指新粒子不受其他粒子或當(dāng)前非劣解集中粒子的支配時,將新粒子放入非劣解集中,每次粒子更新前都從非劣解集中隨機選擇一個粒子作為最優(yōu)粒子。

群體中的每一個粒子更新和調(diào)整自己的速度和位置,公式為

vij(t+1)=wvij(t)+c1r1(t)[pij(t)-xij(t)]+

c2r2(t)[pgj(t)-xij(t)]

(7)

Xij(t+1)=Xij(t)+vij(t+1)

(8)

式中:w為慣性權(quán)重,防止算法早熟收斂;c1、c2為學(xué)習(xí)因子,表示粒子速度轉(zhuǎn)向的偏好程度,均取為1.5;r1、r2為[0,1]區(qū)間上的隨機數(shù),增加了粒子飛行的隨機性;vij為第i個粒子在第j維搜索空間的飛行速度,vij∈[-10,10];Xi為第i個粒子的當(dāng)前位置,i=1,2,…,N;pij為第i個粒子在第j維搜索空間經(jīng)歷的最優(yōu)位置;pgj為整個群體在第j維搜索空間經(jīng)歷的最優(yōu)位置。算法的終止條件為迭代過程中取得最大的適應(yīng)度,即獲得粒子群最佳位置。

裝置的安裝數(shù)量和安裝位置是離散問題,因此采用離散化思想,對變電站節(jié)點進(jìn)行二進(jìn)制(0-1)編碼,0表示不安裝,1表示安裝,粒子按照適應(yīng)度函數(shù)進(jìn)化的方向?qū)?yōu),最終1的數(shù)量就是需要安裝限流裝置的數(shù)量,1的位置就是限流裝置安裝的位置。位置更新等式為

s(vij)=1/(1+exp(-vij))

(9)

(10)

式中r是從[0,1]的均勻分布中產(chǎn)生的隨機數(shù)。

為了使算法不陷入局部最優(yōu)解,本文在算法中利用了帶慣性策略的計算公式[19]

(11)

式中:wmax=0.9為最大權(quán)重;wmin=0.4為最小權(quán)重;Tmax為最大進(jìn)化代數(shù);t為當(dāng)前迭代次數(shù)。

計算開始時使用較大的慣性權(quán)重保證算法在全局范圍內(nèi)探測到較好的區(qū)域,后期使用較小的慣性權(quán)重保證算法在極值點周圍做精細(xì)的搜索,從而增大收斂于全局最優(yōu)解的概率。多目標(biāo)離散粒子群算法的流程如圖3所示。

圖3 多目標(biāo)離散粒子群算法流程

粒子數(shù)為50,最大迭代數(shù)為100。感應(yīng)電場E在0°～360°范圍內(nèi)旋轉(zhuǎn),強度為0.5 V/km。若電場與正東方向的夾角為α,則電場在正東和正北方向上的分量為

EE=Ecosα;EN=Esinα

(12)

文獻(xiàn)[20]規(guī)定,電阻型限流裝置阻值一般不超過3.5 Ω。本文計算中電阻為3 Ω,在目標(biāo)函數(shù)的計算程序中將接地電阻公式加上相應(yīng)的限流電阻值,關(guān)聯(lián)矩陣保持不變。在優(yōu)化程序中調(diào)用目標(biāo)函數(shù),計算各變電站節(jié)點每一角度下的GIC,并判斷每個變電站的360個GIC是否都小于各變電站節(jié)點的GIC限值。如果小于限值,則不需要安裝電阻限流裝置,反之則安裝,直至所有節(jié)點每一角度下的GIC都在直流限值以內(nèi)。

2.3 優(yōu)化結(jié)果及分析

多次運行程序,電阻限流裝置數(shù)和變壓器中性點GIC總和隨迭代數(shù)的變化如圖4所示?？梢钥闯?電阻限流裝置數(shù)和中性點GIC總和呈反比,最終裝置數(shù)收斂于3,GIC總和收斂于223.11 A。

圖4 電阻限流裝置數(shù)和變壓器中性點GIC總和隨迭代次數(shù)的變化

計算得到兩個目標(biāo)函數(shù)的Pareto最優(yōu)解,如圖5所示。實線表示Pareto最優(yōu)解的變化趨勢,陰影區(qū)域內(nèi)為可行解,可行解的兩個目標(biāo)函數(shù)無法同時滿足非劣。從變化趨勢可以看出,兩個目標(biāo)函數(shù)滿足雙目標(biāo)函數(shù)最小化的規(guī)律,一個目標(biāo)函數(shù)增大勢必會造成另一個目標(biāo)函數(shù)的減小。點A、B、C、D均為Pareto最優(yōu)解,這4個解的兩個目標(biāo)之間互不支配、互不占優(yōu),并且可以支配陰影面內(nèi)的任何一點。當(dāng)電阻限流裝置數(shù)為1、2時,雖然GIC總量分別為377.06、271.67 A,都在限值以下,但是未滿足算法的約束條件,因此這兩個解也僅屬于可行解。算法約定首先滿足電阻限流裝置的數(shù)量最少,因此A點為最終解,即在SUB5、SUB6、SUB8安裝電阻限流裝置,IGIC分別為0、42.5、48、55.6、13.6、20.9、42.51 A,IGIC總和為223.11 A。

圖5 兩個目標(biāo)函數(shù)的Pareto最優(yōu)解

圖4和圖5表明,GIC-Benchmark電網(wǎng)在0.5 V/km的地電場作用下,最少需要3臺電阻限流裝置,才能將各個變電站在每一角度下的GIC限制在可承受的直流量內(nèi),并且GIC總和最小。

圖6 E=0.5 V/km時不同治理方案下的GIC總和

旋轉(zhuǎn)電場在變電站和輸電線路中產(chǎn)生的GIC具有一定的對稱性,電場方向與線路走向夾角越小則GIC越大[21]。由于本算例和正東向、正北向的夾角為0°時產(chǎn)生的GIC最大,因此這里只分析正東和正北電場產(chǎn)生的GIC總和的治理效果。GIC-Benchmark在0.5 V/km的電場作用下,未治理、傳統(tǒng)電容方案和電阻優(yōu)化方案的中性點GIC總和如圖6所示,圖中GIC均為三相絕對值?？梢钥闯?兩種電容治理方案都會使SUB3的GIC超標(biāo),這是風(fēng)險轉(zhuǎn)移的結(jié)果;采用電阻優(yōu)化方案治理后,電網(wǎng)中所有變電站的GIC均在可承受范圍內(nèi),可達(dá)到全局最優(yōu)。

3 電阻優(yōu)化方案的效果評價

3.1 變化電場下暫態(tài)響應(yīng)的規(guī)律研究

為了進(jìn)一步對比傳統(tǒng)電容和電阻優(yōu)化兩種抑制方案在變化電場作用下對GIC的抑制效果,選擇傳統(tǒng)隔直方案中的方案2(在SUB5、SUB6、SUB8安裝電容隔直裝置)和電阻限流優(yōu)化方案,與未治理的GIC進(jìn)行比較驗證。在Powerworld Simulator軟件中建立Benchmark標(biāo)準(zhǔn)算例模型,如圖7所示。圖中:G1為平衡機,單位為MV的量為有功功率,單位為MV·A的量為無功功率,各發(fā)電機功率范圍如表2所示。部分暫態(tài)參數(shù)如下[22]:交軸同步電抗Xq=0.5 Ω;直軸同步電抗Xd=2.1 Ω;交軸暫態(tài)電抗Xqp=0.5 Ω;直軸暫態(tài)電抗Xdp=0.2 Ω;直軸次暫態(tài)電抗Xdpp=0.18 Ω;暫態(tài)開路時間常數(shù)Tdop=7.0 s,Tqop=0.75 s;暫態(tài)短路時間常數(shù)Tdopp=0.035 s,Tqopp=0.05 s。

圖7 Powerworld Simulator軟件中建立的Benchmark仿真電路

表2 發(fā)電機功率范圍

為了與第2節(jié)的分析條件保持一致,基于線路和電場之間的夾角為0°時產(chǎn)生最大GIC的事實,設(shè)置地電場方向為正東方向(北向同理),在第0 s開始作用,考慮3種地磁擾動情況:①E在1 s內(nèi)從0增加到1 V/km;②E在2 s內(nèi)從0增加到1 V/km;③E在1 s內(nèi)從0增加到0.5 V/km。

圖8 情況①和③下的G1無功出力

治理前,為了得出不同變化率的E對暫態(tài)響應(yīng)的影響,在情況①②③下計算發(fā)電機無功出力和母線電壓變化,選擇平衡機G1和母線2的暫態(tài)響應(yīng)規(guī)律研究。對比情況①③得到圖8和圖9,對比情況①②得到圖10和圖11。從圖8～11可知:在感應(yīng)地電場影響下,發(fā)電機母線電壓迅速下降,最低可至0.63(標(biāo)幺值),低于0.95(標(biāo)幺值)的持續(xù)時間約為1.2 s;平衡機G1在受到擾動時為調(diào)節(jié)整個電網(wǎng)的無功出力,會大幅增加自身的無功輸出,電場強度越大、變化率越大,則發(fā)電機無功出力越多,超出限值時刻越早。

情況③下3種方案的發(fā)電機無功出力、母線電壓響應(yīng)如圖12和圖13所示?？梢钥闯?在未治理和電容治理時,發(fā)電機G1都超出了無功出力限值,采用電容治理超出限值的時刻更早,驗證了個別節(jié)點經(jīng)電容治理后會導(dǎo)致風(fēng)險轉(zhuǎn)移;采用電阻優(yōu)化治理后,所有發(fā)電機的無功出力均在限值以內(nèi);采用電阻治理方案,母線2電壓降幅最小,且跌落至低于0.95(標(biāo)幺值)的時間最短。

圖9 情況①和③下母線2的電壓響應(yīng)

圖10 情況①和②下G1的無功出力

圖11 情況①和②下母線2的電壓響應(yīng)

圖12 3種方案下G1的無功出力

圖13 3種方案下母線2的電壓響應(yīng)

3.2 優(yōu)化方案的效果評價指標(biāo)

(1)發(fā)電機無功剩余量Qrest。發(fā)電機的無功出力會在磁暴來臨時瞬時增大,無功儲備剩余量減小。因此,用發(fā)電機無功剩余量指標(biāo)Qrest判斷系統(tǒng)無功儲備的充裕性。Qrest公式為

Qrest=Qlim-Qmax

(13)

式中:Qlim為發(fā)電機無功儲備的限值;Qmax為暫態(tài)過程中無功出力的最大值。

3種方案下的發(fā)電機無功剩余量如表3所示,正值代表剩余量,負(fù)值代表超標(biāo)量?？梢钥闯?隨著電場強度的增大,發(fā)電機無功出力增大,剩余量減小;在1 V/km的電場作用下,平衡機在未治理時無功出力大幅增加,超出無功限值1 366.21 MV·A,幾乎為無功限值的3倍,電容治理比未治理的無功出力超標(biāo)更多;所有發(fā)電機經(jīng)電阻優(yōu)化治理后出力均在限值以內(nèi)。若電場強度進(jìn)一步增加,可能導(dǎo)致某些發(fā)電機因無功儲備不夠而造成機端電壓下降,因此暫態(tài)指標(biāo)值得進(jìn)一步關(guān)注。

表3 3種方案下的發(fā)電機無功剩余量

(2)節(jié)點電壓偏移指數(shù)Uofffset。地磁暴會造成節(jié)點電壓下降,節(jié)點電壓偏移指數(shù)定義為

(14)

式中:Umax(i)為節(jié)點電壓下降的最大值;Ue(i)為電壓額定值;i為節(jié)點編號。Uoffset(i)越大說明該節(jié)點電壓偏移越嚴(yán)重,系統(tǒng)越不穩(wěn)定。

3種方案下母線2的電壓偏移如表4所示,可以看出:電場強度越大,則節(jié)點電壓偏移越嚴(yán)重,電壓暫降幅值越大,系統(tǒng)穩(wěn)定性越易受到威脅;經(jīng)電阻優(yōu)化治理后,電壓偏移減小。

表4 3種方案下母線2的電壓偏移

(3)節(jié)點電壓暫降持續(xù)時間Tsag。Tsag為母線電壓低于0.95(標(biāo)幺值)的持續(xù)時間,持續(xù)時間越長說明暫降越嚴(yán)重。Tsag定義為

Tsag=t2-t1

(15)

式中t1和t2為母線電壓與0.95(標(biāo)幺值)的交點。

3種方案下母線2的電壓暫降持續(xù)時間如表5所示,可以看出:電場強度越大,節(jié)點電壓暫降持續(xù)時間越長,系統(tǒng)越不穩(wěn)定;兩種電場變化情況下,暫降時間從長到短依次為未治理、傳統(tǒng)電容治理、電阻優(yōu)化治理,說明電阻優(yōu)化治理方案有效,系統(tǒng)更易恢復(fù)穩(wěn)定狀態(tài)。

表5 3種方案下母線2的電壓暫降持續(xù)時間

4 結(jié) 論

(1)對于GIC-Benchmark采用傳統(tǒng)電容隔直方案抑制GIC,結(jié)果表明,無論是高風(fēng)險節(jié)點安裝電容或者越限就安裝電容,都不可避免地會導(dǎo)致GIC風(fēng)險轉(zhuǎn)移,因此有必要對抑制措施進(jìn)行優(yōu)化配置。

(2)計算得到在0°～360°旋轉(zhuǎn)、0.5 V/km的均勻電場作用下的Pareto最優(yōu)解,最少需要3臺電阻限流裝置才能將每個變電站超標(biāo)的GIC降至安全范圍內(nèi),且GIC總量最小,為223.11 A。

(3)對采取基于多目標(biāo)離散粒子群算法的電阻方案的Benchmark電網(wǎng)進(jìn)行了暫態(tài)分析,結(jié)果表明:在1 V/km電場作用下,未治理時,發(fā)電機無功出力從超出無功限值的3倍降低至剩余量為限值的20%;節(jié)點電壓偏移指數(shù)最少降低了40%;電壓暫降持續(xù)時間下降了30%以上;對0.5 V/km電場擾動的改善效果更為明顯。電阻優(yōu)化方案可行且有效。