張杰愷,李輝,鄧吉利,賀勝,曹林,何蓓
(1.國網(wǎng)重慶江北供電公司,重慶401123;2.重慶大學(xué)輸配電裝備及系統(tǒng)安全與新技術(shù)國家重點(diǎn)實(shí)驗(yàn)室,重慶400044;3.國網(wǎng)重慶市電力公司,重慶401123)
Rogowski線圈電子式電流互感器不含鐵心,消除了磁飽和、鐵磁諧振問題,具有頻率響應(yīng)范圍寬,暫態(tài)響應(yīng)范圍大,測量精度高等優(yōu)點(diǎn)[1]。電子式電流互感器作為系統(tǒng)信息采集的主要設(shè)備,將逐漸取代傳統(tǒng)的電磁式互感器,在電力系統(tǒng)的電能計量、測量和繼電保護(hù)中發(fā)揮重要作用[2]。積分環(huán)節(jié)作為Rogowski線圈電子式電流互感器的關(guān)鍵環(huán)節(jié)之一,不同廠家制作上存在差異,使得相差比差不同,影響電能計量、測量和繼電保護(hù)的準(zhǔn)確性。此外,積分環(huán)節(jié)往往導(dǎo)致“拖尾效應(yīng)”,當(dāng)電網(wǎng)中電流突變,富含直流分量或諧波分量時,不同積分算法的動態(tài)響應(yīng)能力可能不同。因此全面比較研究不同積分算法下電子式互感器的傳變特性,對電子式互感器的優(yōu)化設(shè)計和運(yùn)行監(jiān)測以及適應(yīng)性保護(hù)的研究都具有重要現(xiàn)實(shí)意義。
目前,現(xiàn)有學(xué)者對電子式電流互感器的積分環(huán)節(jié)提出了較多不同算法,如模擬積分[3]、數(shù)字積分矩形算法[4]、梯形算法[5]、Simpson算法[6]、雙線性變換算法[7-8]和 Al-Alaoui[9]算法等。由于模擬器件受時漂和溫漂的影響較大,在實(shí)際中很少采用。文獻(xiàn)[5]設(shè)計了基于STM32F103梯形算法的數(shù)字積分器,但并沒有與其他積分算法進(jìn)行比較與分析。文獻(xiàn)[6]研究了3種常用的矩形算法、梯形算法、Simpson算法的電子式互感器特性,但是文獻(xiàn)僅涉及了不同積分算法幅頻特性的比較,沒有深入分析含不同積分算法的互感器頻帶特性,對其動態(tài)響應(yīng)能力的比較沒有涉及。文獻(xiàn)[7]研究了雙線性變換算法的幅頻特性和諧波特性,但是僅從精度上與其它積分算法進(jìn)行了比較。文獻(xiàn)[10]設(shè)計了對Rogowski線圈電流互感器的精確度、飽和特性、頻率特性和暫態(tài)特性的測試方案,但是文中僅涉及互感器等效電路本身,沒有考慮數(shù)字積分環(huán)節(jié)對其傳變特性的影響。綜上所述,對現(xiàn)有不同積分算法的電子式電流互感器整體的頻帶特性還需要完善,對其動態(tài)響應(yīng)能力的比較還有待深入研究。
本文以Rogowski線圈電子式電流互感器為例,對常用的矩形、梯形、Simpson和雙線性變換積分算法進(jìn)行比較分析,全面評估不同積分算法可能對其傳變特性的影響。首先,分別獲取Rogowski線圈等效電路和積分算法的傳遞函數(shù);其次,對各個環(huán)節(jié)幅頻和相頻特性進(jìn)行分析,獲取并比較了含不同積分算法的電子式互感器整體的頻帶特性;最后,以某實(shí)際110 kV變電站用電子式電流互感器參數(shù)為研究對象進(jìn)行分析,針對穩(wěn)態(tài)和信號突變的動態(tài)情況,對含不同積分算法的電子式互感器傳變特性進(jìn)行仿真。
所采用的Rogowski線圈電子式電流互感器的結(jié)構(gòu)如圖1所示。被測電流經(jīng)Rogowski線圈輸出微分后的信號,再通過A/D采樣和數(shù)字積分器輸出與被測電流成一定比例的信號,最后通過合并單元作進(jìn)一步處理。
圖1 Rogowski線圈電子式電流互感器結(jié)構(gòu)Fig.1 Structure of Rogowski coil electronic current transformers
Rogowski線圈的等效電路如圖2所示,其中L0、R0、C0分別為Rogowski線圈的自感、內(nèi)阻和雜散電容,Rs為取樣電阻,M為Rogowski線圈的互感系數(shù),e(t)為 Rogowski線圈的感應(yīng)電勢,u0(t)為 Rogowski線圈的輸出電壓,I為被測電流。
圖2 Rogowski線圈的等效電路圖Fig.2 Equivalent circuit diagram of Rogowski coil
通過基本電路分析和拉氏變換可得到Rogowski線圈的傳遞函數(shù):
用于電力系統(tǒng)的Rogowski線圈電子式電流互感器常工作于外積分方式,此時Rogowski線圈的輸出電壓與被測電流滿足微分關(guān)系,即:
由上式可知,為了還原被測電流,就需要對輸出電壓進(jìn)行積分。
數(shù)字積分通過軟件算法來實(shí)現(xiàn),且基本不受溫度和其他因素的影響,因而得到廣泛的應(yīng)用。常用的數(shù)字積分算法主要有矩形算法、梯形算法和Simpson算法。根據(jù)它們的z傳遞函數(shù),依次得到其遞推公式如下。
矩形算法:
梯形算法:
Simpson算法:
式中 T為采樣間隔,n表示采樣時刻,X(n)與Y(n)分別表示n時刻積分環(huán)節(jié)的輸入和輸出,下標(biāo)r、t和s分別表示矩形算法、梯形算法和Simpson算法。
為了抑制低頻信號的干擾,最近文獻(xiàn)提出了雙線性變換積分算法[7-8]。其主要思想是:使s平面的虛軸映射為z平面的單位圓,從而實(shí)現(xiàn)模擬域到數(shù)字域的映射。s與z的關(guān)系公式為:
利用雙線性變換,將公式(6)代入某模擬積分公式中:
并得到對應(yīng)的遞推公式為:
式中 α=-TRf/(TR+2RRfC),β=-(TR-2RRfC)/(TR+2RRfC),R、C和Rf分別為對應(yīng)模擬積分的電阻、電容和反饋電阻值。
為了全面分析比較不同積分算法的電子式電流互感器的傳變特性,本文采用式(3)~式(5)和式(8)的數(shù)字積分傳遞函數(shù)對Rogowski線圈電子式電流互感器進(jìn)行頻帶特性以及不同信號源下傳變特性的研究。
按照前述Rogowski線圈的等效電路和第3節(jié)的仿真參數(shù)可得Rogowski線圈環(huán)節(jié)頻帶特性,其Bode圖如圖3所示。
從相頻特性可以看出,Rogowski線圈的輸出電壓與被測電流約有90°的相位差,存在微分關(guān)系。
按照4種數(shù)字積分的傳遞函數(shù)和第3節(jié)的仿真參數(shù),可得4種積分算法的頻帶特性,其Bode圖如圖4所示。
圖3 Rogowski線圈環(huán)節(jié)的Bode圖Fig.3 Bode diagram of Rogowski coil
圖4 4種數(shù)字積分算法的Bode圖Fig.4 Bode diagram of four kinds of digital integrators
從上述幅頻特性可以看出,這4種不同積分算法的對數(shù)幅頻曲線斜率均為-20 dB/dec,但是它們在低頻響應(yīng)上存在差異,特別是雙線性變換積分算法由于是根據(jù)模擬積分變換而得,因此它具有與模擬積分相類似的頻帶特性,即在頻率f<1 Hz區(qū)域?qū)?shù)幅頻曲線斜率為0 dB/dec,而常用的數(shù)字積分算法對數(shù)幅頻曲線斜率為-20 dB/dec。
另外,從上述相頻特性可以看出,這4種不同積分算法都達(dá)到了-90°的相位延遲,但是它們在低頻響應(yīng)和高頻區(qū)域上仍然存在差異。在低頻區(qū)域,采用的雙線性變換積分算法的相位延遲幾乎為0°,表示在低頻區(qū)域,該算法喪失積分功能;在高頻區(qū)域,矩形算法的相位則隨著頻率的增大而增大。
因此,從上述Bode圖可以看出,不同數(shù)字積分在低頻和高頻區(qū)域存在一定的差異,其中雙線性變換積分算法反映在低頻方面的差異;而矩形算法反映在高頻區(qū)域。
基于上述4種不同積分算法的傳遞函數(shù),聯(lián)合Rogowski線圈環(huán)節(jié)傳遞函數(shù),以及第3節(jié)的仿真參數(shù)可得含積分環(huán)節(jié)的電子式電流互感器整體傳遞函數(shù)的頻帶特性,其Bode圖如圖5所示。
圖5 含不同積分算法的電子式電流互感器整體Bode圖Fig.5 Bode diagram of electronic current transformer with different integrators
從圖5中可以看出,在頻率10 Hz~10 kHz的范圍內(nèi),4種積分算法的互感器幅頻特性和相頻特性表現(xiàn)為水平直線,可見在上述頻帶范圍內(nèi),不同積分算法時電子式電流互感器輸出信號的幅值和相位不受頻率變化的影響,即能夠準(zhǔn)確地還原被測電流中的各種頻率分量,頻帶特性較好。然而,當(dāng)頻率小于10 Hz時,相對其他積分算法,采用雙線性變換積分算法時的還原特性會更容易受到原始信號低頻分量的影響;而當(dāng)頻率大于10 kHz時在高頻區(qū)域,上述4種積分算法的互感器輸出的相位則會隨著頻率的增大而增大,從而很大程度上影響其還原特性。
為了較全面的分析不同積分算法Rogowski線圈電子式電流互感器傳變特性及其影響,本文以實(shí)際110 kV數(shù)字化變電站中的Rogowski線圈電子式電流互感器為例,其仿真用的主要參數(shù)如下[11]:一次額定電流為I=600 A,M=3.369μH,L0=21.56 mH,R0=270Ω,C0=0.12 nF,Rs=10 kΩ,采樣時間 T=0.1 ms,其中當(dāng)采用雙線性變換積分算法時取R=30 kΩ,C=0.1μF,Rf=1 MΩ?;赑LECS仿真平臺搭建了含上述4種不同數(shù)字積分算法的Rogowski線圈電子式電流互感器的模型,下面針對穩(wěn)態(tài)信號、積分初值不確定,以及信號突變或富含直流、諧波信號不同情況,對不同積分算法的互感器傳變特性進(jìn)行仿真。
現(xiàn)選擇初始相位角為0°,頻率為50 Hz,標(biāo)幺值為1的額定電流 i(t)=sin(100πt)作為輸入信號,4種不同積分算法的互感器輸入-輸出波形如圖6所示。
圖6 穩(wěn)態(tài)特性下互感器傳變特性比較Fig.6 Comparisons of the transducer transfer characteristics under the steady-state input condition
從結(jié)果可以看出,4種積分算法的互感器穩(wěn)態(tài)響應(yīng)良好,均能較好的還原被測波形。其比差和相差如表1所示,通過表中計算結(jié)果可知,4種積分算法的電流互感器在額定電流下都能滿足較高的準(zhǔn)確度要求。
表1 額定電流下互感器的準(zhǔn)確度Tab.1 Accuracy of transformers at rated current
由于電子式電流互感器還原被測電流需要對Rogowski線圈的輸出電壓進(jìn)行積分處理,即:
因此需要確定數(shù)字積分的初值。為了驗(yàn)證不同積分初值對上述積分算法的魯棒性,現(xiàn)假定輸入電流存在 90°初始相位,即 i(t)=sin(100πt+90°),而積分算法仍假設(shè)初值為0,圖7得到了4種積分算法下Rogowski線圈電子式電流互感器傳變特性。
圖7 積分初值不確定下互感器傳變特性比較Fig.7 Comparisons of the transducer transfer aracteristicswith the unknown initial integration points
從圖7可以看出,采用雙線性變換積分算法時互感器可以在較短時間內(nèi)的還原輸入電流,具有積分初值不確定下魯棒性好的優(yōu)點(diǎn)。而采用矩形算法、梯形算法和Simpson算法的互感器輸出則受積分初值不確定下影響大,輸出信號含有50%幅值的直流分量,導(dǎo)致其傳變特性輸出信號幅值下移50%,還原信號能力明顯降低。
由于電網(wǎng)中的功率控制或負(fù)荷變化等各種因素下可使互感器的輸入電流突然增加或者降低,為了仿真分析電流幅值突變時對不同積分算法下互感器傳變特性的動態(tài)跟蹤能力比較,現(xiàn)假定輸入電流為i(t)=sin(100πt),在 t=0.044 s時刻電流幅值從 1 p.u.突變?yōu)?.5 p.u.,在t=0.104 s時刻電流幅值又從1.5 p.u.突變回1 p.u.,仿真比較4種不同積分算法時互感器動態(tài)響應(yīng)情況,如圖8所示。
圖8 輸入信號幅值突變下互感器傳變特性比較Fig.8 Comparisons of the transducer transfer characteristicswhen the input amplitude is changed
從圖8中的仿真波形可以看出,4種不同積分算法的互感器在電流幅值突變情況下輸出波形與輸入波形吻合,能較好地還原輸入信號。但是從其動態(tài)響應(yīng)能力可以看出(觀察放大波形),4種不同積分算法的互感器動態(tài)響應(yīng)和幅值跟蹤方面仍存在一定的差異,與其余3中積分算法比較,Simpson算法的幅值突變的跟蹤響應(yīng)能力較差。
當(dāng)線路中發(fā)生過電壓或故障時,由于電感或電容的存在,往往使得電流不能突變而導(dǎo)致線路中存在直流分量的情況。此外,由于Rogowski線圈輸出信號弱,在外界干擾下也可能導(dǎo)致積分環(huán)節(jié)的輸入信號中含有直流分量。為了比較分析富含直流分量時互感器傳變特性,下面對Rogowski線圈輸入信號和積分環(huán)節(jié)輸入信號富含直流分量時的情況進(jìn)行仿真,如圖9所示。其中,圖9(a)為輸入電流i(t)=sin(100πt),在 t=0.1 s開始含有 1 p.u.衰減時間常數(shù)T=0.05 s的直流分量時傳感器傳變特性的仿真結(jié)果;圖9(b)為假定在t=0.1 s時積分環(huán)節(jié)輸入信號含10%幅值的穩(wěn)定直流量干擾下傳感器傳變特性的仿真結(jié)果。
圖9 含不同直流分量時互感器傳變特性比較Fig.9 Comparisons of the transducer transfer characteristicswhen the input has different DC signals
從圖9(a)可以看出,當(dāng)Rogowski線圈輸入含衰減直流分量時,不同積分算法基本上都能還原其輸入信號的變化,但相比之下,采用雙線性變換積分算法時還原能力相對較差。
從圖9(b)可以看出,當(dāng)積分環(huán)節(jié)輸入含穩(wěn)定直流分量干擾情況下,雙線性變換積分算法受直流量干擾的影響相對較小,還原特性好,而其余3種算法受直流量干擾的影響較大,導(dǎo)致其互感器輸出信號隨著時間推移不斷增大,呈現(xiàn)明顯發(fā)散的傳變特性。
為了模擬線路中存在諧波分量時互感器的動態(tài)性能,現(xiàn)假定線路中存在2~13次的諧波分量,即輸入電流為 i(t)=sin(100πt),并在 t=0.1 s開始輸入幅值分別為 1/2、1/3、1/4、1/5、1/6、1/7、1/8、1/9、1/10、1/11、1/12和1/13的2~13次諧波分量,得到圖10所示仿真波形。從圖中可以看到,采用4種積分算法的互感器均能很好的還原被測電流中的低次諧波分量。
圖10 含2-13次諧波分量時互感器傳變特性比較Fig.10 Comparisons of the transducer transfer characteristicswhen the input has 2-13 order harmonic components
另外,當(dāng)線路中電感、電容存在導(dǎo)致諧振可能會存在高次諧波情況,為了分析不同積分算法下互感器對線路中高次諧波分量的影響,現(xiàn)假定輸入電流為 i(t)=sin(100πt),在 t=0.1 s開始輸入頻率為106Hz、幅值為0.1的高次諧波,采用不同積分算法時其互感器傳變特性如圖11所示。
從圖中可以看出,采用上述4種積分算法的互感器還原高次諧波信號的能力都較差,這與圖5中含不同積分算法時互感器整體輸出的頻帶特性分析結(jié)果一致,即當(dāng)Rogowski線圈輸入信號含高次諧波分量時,上述不同積分算法并不能很好反映其輸出高次諧波特性。
圖11 含高次諧波分量時互感器傳變特性比較Fig.11 Comparisons of the transducer transfer characteristicswhen the input has a high-order harmonic component
以Rogowski線圈電子式電流互感器為例,通過建立傳遞函數(shù)和仿真模型,對常用的矩形、梯形、Simpson和雙線性變換積分算法下的互感器頻帶特性和傳變特性時域仿真進(jìn)行了比較分析,得到以下結(jié)論:
(1)從頻帶特性上和輸入信號富含諧波分量時的仿真比較可以看出,在頻率10 Hz~10 kHz的范圍內(nèi),上述4種積分算法的電子式電流互感器幅值和相位不受頻率變化的影響,即能夠較準(zhǔn)確地還原被測電流中的各種頻率分量,頻帶特性較好,但在低頻以及高頻率信號下還原能力較差;
(2)從穩(wěn)態(tài)特性的比較分析可以看出,4種積分算法的互感器還原輸入電流的能力都較好,在額定電流下都能滿足較高的準(zhǔn)確度要求。但在積分初值不確定的情況下,采用雙線性變換積分算法的互感器傳變特性較好,魯棒性較高;
(3)從不同輸入信號突變的動態(tài)跟蹤能力的比較分析可以看出,當(dāng)輸入信號幅值突變時,不同積分算法的互感器能較好地還原特性,相比Simpson算法存在較大的誤差;當(dāng)輸入信號含直流分量時,采用雙線性變換積分算法的互感器傳變特性較好,尤其是在積分環(huán)節(jié)輸入有直流分量情況。
從上述結(jié)論可以看出,采用不同積分算法時電子式電流互感器的頻帶特性和穩(wěn)態(tài)傳變特性較好,并沒有顯著的差別;但在動態(tài)特性方面,尤其在積分初值不確定、輸入信號富含直流分量和高頻諧波分量時,能較好地適應(yīng)其抗干擾能力和動態(tài)跟蹤能力方面的數(shù)字積分算法還有待進(jìn)一步研究。