賴勝杰 夏成軍 紀煥聰 王澤青
(1.華南理工大學電力學院 廣州 510640 2.廣東省新能源電力系統(tǒng)智能運行與控制企業(yè)重點實驗室 廣州 510663)
我國配電網(wǎng)遵循“閉環(huán)設(shè)計,開環(huán)運行”的原則[1]。當配電網(wǎng)線路長時間高負荷運行或設(shè)備計劃檢修時,需對負荷進行轉(zhuǎn)供,目前廣泛采用“先通后斷”的方式,即合環(huán)轉(zhuǎn)供電[2-3]。合環(huán)轉(zhuǎn)供電可在供電不間斷的情況下,實現(xiàn)負荷轉(zhuǎn)移,滿足電力用戶對供電可靠性的要求[4-5]。但合環(huán)過程中可能會產(chǎn)生較大的穩(wěn)態(tài)環(huán)流及沖擊電流,導致繼電保護誤動、電力設(shè)備過載,造成合環(huán)失敗[6]。因此,合環(huán)穩(wěn)態(tài)電流與沖擊電流計算是進行合環(huán)分析的基礎(chǔ),而配電網(wǎng)合環(huán)模型則是進行合環(huán)電流計算的關(guān)鍵。
現(xiàn)有的合環(huán)穩(wěn)態(tài)電流分析模型基本圍繞穩(wěn)態(tài)環(huán)流計算展開。文獻[7]在計算合環(huán)前聯(lián)絡(luò)開關(guān)兩側(cè)的電壓和穩(wěn)態(tài)環(huán)流時直接忽略負荷的作用,僅保留合環(huán)環(huán)路中的阻抗,形成簡化的合環(huán)模型,算出的合環(huán)穩(wěn)態(tài)電流可能會有較大誤差。文獻[8]將饋線上的負荷處理為一個集中負荷,接在饋線的首端或末端,形成四種合環(huán)分析模型,然后考慮合環(huán)電流最嚴重的情形,進行保守的合環(huán)風險判斷。文獻[9]考慮了三種特殊的負荷分布情形,即均勻分布、遞增分布、遞減分布,分析了不同負荷分布情形下的合環(huán)穩(wěn)態(tài)電流。文獻[8-9]分析了負荷分布對合環(huán)前聯(lián)絡(luò)開關(guān)兩側(cè)電壓的影響,而在穩(wěn)態(tài)環(huán)流的計算模型中,沒有考慮負荷的分流作用。
對于合環(huán)沖擊電流,文獻[10]建立復頻域下的戴維南等效電路,得到合環(huán)支路的電流表達式,再通過拉普拉斯反變換,得到時域下的全電流表達式。文獻[11]建立時域下的戴維南等效電路,利用三要素公式得到合環(huán)支路以及非合環(huán)支路全電流表達式,進而求得各支路的沖擊電流。文獻[12]建立時域下戴維南等效電路的微分方程,得到最嚴重情況下的沖擊電流,即用穩(wěn)態(tài)電流幅值乘以沖擊系數(shù)。文獻[10-12]求解沖擊電流時所用的分析模型都是忽略負荷阻抗的簡化戴維南等效電路,而事實上沖擊電流的衰減時間常數(shù),不僅與合環(huán)環(huán)路的阻抗有關(guān),還受負荷阻抗的影響。文獻[13]提出了一種基于最佳頻率法求解合環(huán)過程中衰減時間常數(shù)的方法,但所述求解方法比較繁瑣,需對合環(huán)前后最佳頻率下的網(wǎng)絡(luò)進行潮流計算。
饋線負荷具有隨機特性,隨時間有較大的波動性[14-17],一般難以獲取饋線每個負荷的實時數(shù)據(jù),處理起來相對困難,故現(xiàn)有的合環(huán)模型基本都忽略負荷。隨著實時同步測量技術(shù)的發(fā)展,通過微型同步相量測量裝置(micro-Phasor Measurement Unit,μPMU)可實時獲取配電網(wǎng)節(jié)點相量數(shù)據(jù),μPMU 以其優(yōu)異性能在配電網(wǎng)具有廣闊的應(yīng)用前景[18-21]。
本文考慮負荷對合環(huán)穩(wěn)態(tài)與暫態(tài)過程的影響,基于饋線首末端少量μPMU 的實時量測信息對負荷進行等效處理,提出計及負荷等值阻抗的合環(huán)轉(zhuǎn)供電分析模型。同時,基于最佳頻率法[13],結(jié)合本文所提的合環(huán)分析模型,給出了求解合環(huán)沖擊電流的一種有效方法,以提高沖擊電流計算的準確性,避免繁瑣的潮流計算過程,最后通過算例驗證了本文所提模型及方法的有效性。
圖1 是典型的合環(huán)示意圖,圖中S1、S2分別為饋線a、b 的出線開關(guān),S 為聯(lián)絡(luò)開關(guān),Sai、Sbj(i,j=1,2,···,n)分別為饋線a、b 所帶負荷。當母線1發(fā)生故障或檢修時,就需要母線2 轉(zhuǎn)帶母線1 的負荷。具體操作流程是:先閉合聯(lián)絡(luò)開關(guān),再斷開S1,將饋線a 的負荷轉(zhuǎn)移到饋線b,以此來保證不間斷供電,這種閉合聯(lián)絡(luò)開關(guān)的方式就是合環(huán)轉(zhuǎn)供電操作[22]。
圖1 典型合環(huán)示意圖 Fig.1 Schematic diagram of a typical loop closing
配電網(wǎng)饋線負荷眾多,往往成梳狀結(jié)構(gòu)[23]。隨著經(jīng)濟的發(fā)展,城市配電網(wǎng)負荷快速增長[24]。當配電網(wǎng)饋線所帶負荷較重時,由于負荷越大所對應(yīng)的負荷阻抗越小,合環(huán)時其分流作用不可忽視,故建立合環(huán)電流分析模型時需將負荷阻抗考慮在內(nèi)。但如果建立每個饋線負荷的負荷模型,則工作量很大,不利于合環(huán)電流的分析計算。采用等效負荷模型[25]可以將一段饋線之間的所有負荷等效為一個負荷,能在簡化計算的同時反映原有的負荷特性,提高合環(huán)模型的準確度。
圖2 是一段簡單梳狀饋線示意圖,UA、UB分別為饋線首末端節(jié)點A、B 的電壓,θA、θB為兩個電壓的相角,PA+jQA、PB+jQB分別為饋線首末端的功率,Pi+jQi為節(jié)點i處的負荷功率,Lij為節(jié)點i和節(jié)點j之間的線路長度,饋線AB 的總長為L,饋線單位長度(km)的阻抗z=r+jx。
圖2 簡單梳狀饋線示意圖 Fig.2 Schematic diagram of simple comb feeder
圖 3 是經(jīng)過負荷等效后的示意圖,其中PK+jQK為等效負荷節(jié)點K 的功率,IK∠θIK為流過等效負荷的電流相量,IA∠θIA為節(jié)點A 流向節(jié)點K 的電流相量,IB∠θIB為節(jié)點K 流向節(jié)點B 的電流相量,L1、L2分別為等效負荷節(jié)點距離節(jié)點A、B 的長度。
圖3 饋線AB 的等效負荷 Fig.3 Equivalent load of feeder AB
負荷等效過程實際上是根據(jù)節(jié)點A、B 的電壓相量和功率推導出等效負荷的大小和位置。隨著μPMU 在配電網(wǎng)的發(fā)展,實時獲取配電網(wǎng)節(jié)點電壓相量成為可能,而要實現(xiàn)饋線AB 的負荷等效,僅需在節(jié)點A、B 兩處配置μPMU 即可。
負荷等效的原則是:由節(jié)點A 電壓降落公式得到的等效負荷節(jié)點電壓UK和由節(jié)點B 電壓降落公式得到的等效負荷節(jié)點電壓相等。其中
這里可以通過迭代法求出L1和L2的值,給定一個誤差函數(shù)和迭代步長ΔL,通過不斷修正L1的長度可以使得誤差函數(shù)滿足給定精度ε,從而確定等效負荷的位置K,計算流程如圖4 所示。
圖4 求解L1、L2 的流程 Fig.4 Flow chart for solving L1 and L2
由圖3 可得
可以求出等效負荷功率SK,即
以上是梳狀饋線的等效負荷求解過程。如果饋線AB 間帶有分支饋線,相當于在分支節(jié)點處接入一個大小與分支饋線首端功率相等的負荷,同樣可以根據(jù)梳狀饋線的方法進行負荷等效。
負荷等效過程是在保證節(jié)點A、B 的電壓和功率不變的前提下進行的,所以等效負荷在節(jié)點A、B 產(chǎn)生的效果與實際負荷分布產(chǎn)生的效果是等效的,不影響合環(huán)電流的計算。等效負荷的大小和位置可以反映饋線中負荷的實際分布情況,如果等效負荷較大且靠近源端,則該段饋線負荷較重且分布偏向源端;反之,饋線負荷較輕且分布偏向末端。當所分析的配電網(wǎng)是負荷密集的大型城市配電網(wǎng)時,負荷對合環(huán)的影響不能忽視,等效負荷就比較適用。
在求得等效負荷參數(shù)后,由功率公式可得到等效負荷的阻抗ZK=RK+jXK,即
也可直接根據(jù)歐姆定律求得,即
該負荷等值阻抗是基于實時量測信息等效,能反映饋線實時負荷分布情況,由此得到饋線AB 的等效電路,如圖5 所示。圖5 中,R1+jX1為饋線L1段的阻抗,R2+jX2為饋線L2段的阻抗。
圖5 饋線AB 的等效電路 Fig.5 Equivalent circuit of feeder AB
基于上述方法,將饋線上的負荷處理為一個等效負荷,并用其等值阻抗代替,將10kV 以上的電源和阻抗折算到10kV 側(cè),得到計及負荷等值阻抗的合環(huán)模型如圖6 所示。
圖6 計及負荷等值阻抗的合環(huán)模型 Fig.6 Loop closing model considering load equivalent impedance
根據(jù)圖6 中的合環(huán)模型,可以求出其ab 端口的戴維南等效電路,如圖7 所示。
圖7 戴維南等效電路 Fig.7 Thevenin equivalent circuit
等值阻抗Zeq的計算需要將負荷等值阻抗計算在內(nèi),具體計算公式為
式中,Zeqa為聯(lián)絡(luò)開關(guān)a 側(cè)等值阻抗;Zeqb為聯(lián)絡(luò)開關(guān)b 側(cè)等值阻抗。當聯(lián)絡(luò)開關(guān)S 閉合且達到穩(wěn)態(tài)時,可求出穩(wěn)態(tài)環(huán)流,即
圖8 合環(huán)后穩(wěn)態(tài)等效電路 Fig.8 Steady-state equivalent circuit after closing loop
圖9 利用疊加原理求解合環(huán)后饋線電流 Fig.9 Using superposition principle to solve the feeder current after closing loop
求解合環(huán)過程中的沖擊電流,需要對圖7 中戴維南等效電路進行暫態(tài)分析,將寫成Emsin(ω t+α)的形式,將Zeq寫成Req+jωLeq的形式,列寫微分方程為
式中,Em、ω、α分別為戴維南等效電壓源的幅值、角頻率和初相角;Req、Leq分別為戴維南等效電阻和電感;ic為合環(huán)支路ab 的全電流。假設(shè)在t=0 時刻合上聯(lián)絡(luò)開關(guān),求解式(13)中的周期分量與非周期分量,相加可得ic為
式中,ip、iap分別為ic的周期分量和非周期分量;Ic為周期分量有效值;φ為戴維南等效電路的阻抗角;C為非周期分量初始值;Ta為衰減時間常數(shù)。(ω t+α-φ)對應(yīng)的相量就是式(10)中求出的,由換路定則可知。
對于衰減時間常數(shù)Ta的計算,現(xiàn)有文獻的方法大都是在工頻下進行求解[10-12],即先求出合環(huán)點處的戴維南等值阻抗Zeq=Req+jωLeq(ω為工頻下的角頻率),則Ta=Leq/Req,該方法多適用于各支路R/X相差不大的網(wǎng)絡(luò),且對于網(wǎng)狀電力網(wǎng)絡(luò)的計算精度不是很高。
最佳頻率法能在各支路R/X相差較大的情況下仍保持Ta有較高精確度,在網(wǎng)狀電力網(wǎng)絡(luò)幾種求解Ta的方法中作為推薦方法[26],當要處理的網(wǎng)絡(luò)是大型城市電網(wǎng)時,使用最佳頻率法能有更高的精確性。最佳頻率法求解Ta的基本思路是:對于合環(huán)后某個時刻,取一個特定的頻率,即最佳頻率fc,將合環(huán)網(wǎng)絡(luò)中所有額定頻率fn下的阻抗轉(zhuǎn)換為fc下的阻抗,求出最佳頻率下的戴維南等值阻抗Zeq(fc)=Req(fc)+jXeq(fc),則衰減時間常數(shù)Ta為
求解沖擊電流時,對于50Hz 的電網(wǎng),最佳頻率fc通常取20Hz[26]。本文求解Ta的方法是先對負荷進行等效處理,得到等效負荷參數(shù)后換算成相應(yīng)的等值阻抗,然后求解Ta,具體求解步驟如下:
(1)按前述方法進行負荷處理,得到計及負荷等值阻抗的合環(huán)模型,如圖6 所示。
(2)將圖6 中所有額定頻率下的阻抗轉(zhuǎn)換成最佳頻率下的阻抗。
(3)利用式(9)計算最佳頻率下的戴維南等值阻抗Zeq(fc)。
(4)根據(jù)式(15)計算衰減時間常數(shù)Ta。
求解合環(huán)電流最大瞬時值icmax,只需調(diào)用Matlab 中的max 函數(shù),即可快速求出式(14)中的最大值。對應(yīng)于50Hz 的頻率,周期為0.02s,實際合環(huán)操作中,控制具體的合閘時刻較為困難,而不同的合閘時刻對應(yīng)不同的非周期分量初值C,所以關(guān)注最嚴重情況下的合環(huán)電流最大值,即合環(huán)沖擊電流。通過對合環(huán)過程的暫態(tài)分析[27],如果假定t=0時刻是最嚴重的合閘時刻,則有α-φ=-π/2,由此得出最嚴重情況下合環(huán)支路ab 的全電流表達式為
利用max 函數(shù)求得式(16)的最大值與最大值出現(xiàn)的精確時刻,該最大值即為支路ab 的合環(huán)沖擊電流icim。
假定上述求得的最大值時刻為tc,則支路ab 合環(huán)電流最大有效值計算公式為
對于流過饋線開關(guān)的沖擊電流,應(yīng)用疊加原理求解。將ic1的全電流與合環(huán)前饋線a 的全電流相加,即可得到合環(huán)后饋線a 的全電流ia表達式;將合環(huán)前饋線b 的全電流減去ic2的全電流,即可得到合環(huán)后饋線b 的全電流ib表達式。假設(shè)在t=0 時刻合上聯(lián)絡(luò)開關(guān),有
式中,Ia0、θa分別為合環(huán)前饋線a 電流的有效值和初相角;Ib0、θb分別為合環(huán)前饋線b 電流的有效值和初相角;Ic1、θc1分別為分布到饋線a 首端電流的有效值和初相角;Ic2、θc2分別為分布到饋線b 首端的電流I˙c2的有效值和初相角;C1、C2分別為ia、ib非周期分量初始值,,。式(18)中前兩項之和對應(yīng)的相量實際上就是合環(huán)后饋線a 的電流相量。式(19)中前兩項之和對應(yīng)的相量為合環(huán)后饋線b 的電流相量。
假定t=0 時刻是最嚴重的合閘時刻,考慮最嚴重情況下的饋線合環(huán)沖擊電流。顯然,當合環(huán)支路中流過最嚴重情況下的環(huán)流ic時,即式(16)中的表達式,有最嚴重的饋線合環(huán)沖擊電流。此時,式(18)中相量最嚴重情況下初相角應(yīng)取-π/2,。合環(huán)前饋線a 電流相量最嚴重情況下的初相角,則可得饋線a 最嚴重情況下的全電流表達式為
同理,可求出饋線b 最嚴重情況下的全電流為
調(diào)用Matlab 里的max 函數(shù),可以求出式(20)和式(21)t=0 時刻后的最大絕對值(饋線b 絕對值最大時ib為負值),即最嚴重情況下饋線a、b 的合環(huán)沖擊電流。這里如果直接用穩(wěn)態(tài)電流幅值乘沖擊系數(shù)來求解最嚴重情況下饋線a 的合環(huán)沖擊電流可能會產(chǎn)生較大誤差,因為式(20)中非周期分量初始值可能與饋線a 合環(huán)后穩(wěn)態(tài)電流幅值相差較大,饋線b 也是如此。
將式(20)中的前兩項相加得到
式中,Ia為饋線a 合環(huán)后穩(wěn)態(tài)電流有效值;θk為式(20)中前兩項相加后電流的初相角,當θa-θc1較大時,θk會偏離-π/2 較大角度,如圖10 所示,此時不滿足合環(huán)電流最大瞬時值發(fā)生在t=0.01s 左右的條件。
圖10 式(20)與式(22)中的電流相量圖 Fig.10 The current phasor diagram in equations(20) and (22)
通過max 函數(shù)找出式(22)最大值出現(xiàn)的精確時刻,設(shè)最大值時刻為ta,則饋線a 合環(huán)電流最大有效值計算公式為
同理,饋線b 合環(huán)電流最大有效值為
式中,tb為饋線b 電流出現(xiàn)最大負值的時刻;Ib為合環(huán)后饋線b 穩(wěn)態(tài)電流有效值。通常,合環(huán)后一側(cè)饋線電流會增大,而另一側(cè)減小,但當環(huán)流很大時,可能會造成一側(cè)潮流反向,導致兩側(cè)電流都增大。若饋線b 合環(huán)后穩(wěn)態(tài)電流小于合環(huán)前穩(wěn)態(tài)電流,可不計算饋線b 的沖擊電流與合環(huán)電流最大有效值。
以圖11 中某具體合環(huán)網(wǎng)絡(luò)為例,進行合環(huán)電流分析,并與PSCAD 的仿真結(jié)果進行對比。網(wǎng)絡(luò)參數(shù)如下:T1、L11、T3折算到 10kV 側(cè)的阻抗為(0.023 2+j0.453 7)Ω,T2、L12、T4折算到10kV 側(cè)的阻抗為(0.023 2+j0.376 7)Ω,圖11 中的饋線均采用YJV22 型電纜,阻抗為(0.097 8+j0.091 1)Ω/km,每段饋線的長度已標注在圖11 中,負荷的分布情況如圖11 所示,功率單位為MV·A,母線上所接的功率為其他饋線的功率之和。合環(huán)前母線1 的線電壓相量為10.60∠-30.34° kV,母線2 的線電壓相量為10.46∠-35.87° kV。
圖11 某具體的合環(huán)網(wǎng)絡(luò) Fig.11 A specific loop closing network
在進行PSCAD 建模時,負荷用相應(yīng)的負荷阻抗表示。在本算例中,仿真時長設(shè)置為2s,通過對合環(huán)暫態(tài)過程的分析,當合環(huán)支路ab 的周期電流分量取幅值時,為最嚴重合閘時刻,其中一個最嚴重合閘時刻為0.998 2s,控制聯(lián)絡(luò)開關(guān)在0.998 2s 閉合。進行合環(huán)電流計算,首先,對兩側(cè)饋線的負荷進行等效處理,按照圖4 所示的流程分別計算出兩側(cè)饋線L1和L2的值,其中,a 側(cè)饋線L1初值取2km,b 側(cè)饋線L1初值取3.5km,并計算出兩側(cè)饋線的等效負荷功率。根據(jù)程序運行結(jié)果和計算結(jié)果,將兩側(cè)饋線的等效負荷參數(shù)記錄在表1 中。
表1 兩側(cè)饋線的等效負荷參數(shù) Tab.1 Equivalent load parameters of feeders on both sides
然后利用本文所提的合環(huán)模型進行合環(huán)電流的計算,并和其他方法對比。合環(huán)穩(wěn)態(tài)電流計算值與仿真實測值的對比見表2,其中均勻分布模型、遞增分布模型、遞減分布模型均為考慮負荷阻抗時的計算結(jié)果。
表2 合環(huán)穩(wěn)態(tài)電流的計算結(jié)果及對比 Tab.2 Calculation results and comparison of loop closing steady-state current
為了比較使用最佳頻率法求解的沖擊電流值與現(xiàn)有文獻中采用工頻求解衰減時間常數(shù)所得沖擊電流值的精確性,在各種不同負荷處理形式下,分別利用現(xiàn)有文獻中工頻求解衰減時間常數(shù)的方法和最佳頻率法求解沖擊電流,并將計算結(jié)果分別記錄在表3 和表4 中。這里計算的合環(huán)沖擊電流與合環(huán)電流最大有效值均是最嚴重情況下的值。在本算例中,由于饋線b 合環(huán)后的穩(wěn)態(tài)電流減小,故無需對饋線b 進行沖擊電流計算。
表3 工頻下合環(huán)沖擊電流計算結(jié)果及對比 Tab.3 Calculation results and comparison of loop closing surge current at power frequency
表4 最佳頻率下合環(huán)沖擊電流計算結(jié)果及對比 Tab.4 Calculation results and comparison of loop closing surge current at the optimal frequency
(續(xù))
由表2~表4 的數(shù)據(jù)分析可得出如下結(jié)論:
(1)各方法計算出來的穩(wěn)態(tài)電流以及沖擊電流值,基本上為負誤差,這在實際運用中需要加以注意,為保證安全合環(huán),可以將計算結(jié)果乘以一個稍大于1 的安全系數(shù)來避免出現(xiàn)負誤差。對于利用本文模型與方法的計算結(jié)果而言,穩(wěn)態(tài)電流以及沖擊電流icim、iaim可以乘以1.005 來避免出現(xiàn)負誤差,而合環(huán)電流最大有效值Icim、Iaim的計算結(jié)果,可以乘以1.05 來避免出現(xiàn)負誤差。
(2)由表3 和表4 的誤差對比可知,對于同一負荷處理形式而言,與現(xiàn)有文獻中工頻求解衰減時間常數(shù)得出的沖擊電流值相比,最佳頻率法求解的沖擊電流誤差更小(忽略負荷阻抗時,兩者計算結(jié)果相同),采用基于本文合環(huán)模型的最佳頻率法求解沖擊電流,與仿真值最接近,驗證了本文最佳頻率法求解沖擊電流方法的有效性。
(3)由表2 和表4 的誤差數(shù)據(jù)可知,采用本文的合環(huán)模型及計算方法,穩(wěn)態(tài)電流及沖擊電流icim、iaim的誤差絕對值在0.3%以內(nèi),合環(huán)電流最大有效值Icim、Iaim的誤差絕對值在3%以內(nèi);忽略負荷阻抗的穩(wěn)態(tài)電流及沖擊電流icim、iaim的誤差絕對值多在2%以上,最高超過4%,合環(huán)電流最大有效值Icim、Iaim的誤差絕對值最大接近5%;其他考慮負荷分布模型的穩(wěn)態(tài)電流及沖擊電流icim、iaim的誤差絕對值多在1%~2%之間,合環(huán)電流最大有效值Icim、Iaim的誤差絕對值多在3%~4%之間。驗證了本文所提合環(huán)模型的有效性。
針對現(xiàn)有配電網(wǎng)合環(huán)模型在負荷處理上存在的不足,本文考慮了饋線負荷對合環(huán)穩(wěn)態(tài)與暫態(tài)過程的影響,基于饋線首末端少量μPMU 的實時量測信息,利用等效負荷模型將饋線上的負荷等效為一個集中負荷,再用等效負荷相應(yīng)的等值阻抗代替,提出了計及負荷等值阻抗的配電網(wǎng)合環(huán)分析模型。在該模型的基礎(chǔ)上,進行合環(huán)穩(wěn)態(tài)電流及沖擊電流的計算。在穩(wěn)態(tài)電流的求解過程中,對疊加定理的應(yīng)用原理進行了詳細分析。對于合環(huán)沖擊電流,重點研究了最嚴重合閘時刻下的計算方法,同時,根據(jù)本文的合環(huán)模型,提出了一種基于最佳頻率法求解沖擊電流的有效方法,可以提高沖擊電流計算的準確性,且不需要進行復雜的潮流計算過程。算例分析結(jié)果表明,采用本文合環(huán)模型及方法計算合環(huán)電流準確性較高,驗證了本文所提模型及計算方法的有效性。
新型電力系統(tǒng)背景下,隨著分布式電源、電動汽車等大量接入配電網(wǎng),需要進一步研究考慮分布式電源出力不確定性以及電動汽車充電隨機性的配電網(wǎng)合環(huán)分析模型。