張杰愷，李輝，鄧吉利，賀勝，曹林，何蓓

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

0 引 言

Rogowski線圈電子式電流互感器不含鐵心，消除了磁飽和、鐵磁諧振問(wèn)題，具有頻率響應(yīng)范圍寬，暫態(tài)響應(yīng)范圍大，測(cè)量精度高等優(yōu)點(diǎn)［1］。電子式電流互感器作為系統(tǒng)信息采集的主要設(shè)備，將逐漸取代傳統(tǒng)的電磁式互感器，在電力系統(tǒng)的電能計(jì)量、測(cè)量和繼電保護(hù)中發(fā)揮重要作用［2］。積分環(huán)節(jié)作為Rogowski線圈電子式電流互感器的關(guān)鍵環(huán)節(jié)之一，不同廠家制作上存在差異，使得相差比差不同，影響電能計(jì)量、測(cè)量和繼電保護(hù)的準(zhǔn)確性。此外，積分環(huán)節(jié)往往導(dǎo)致“拖尾效應(yīng)”，當(dāng)電網(wǎng)中電流突變，富含直流分量或諧波分量時(shí)，不同積分算法的動(dòng)態(tài)響應(yīng)能力可能不同。因此全面比較研究不同積分算法下電子式互感器的傳變特性，對(duì)電子式互感器的優(yōu)化設(shè)計(jì)和運(yùn)行監(jiān)測(cè)以及適應(yīng)性保護(hù)的研究都具有重要現(xiàn)實(shí)意義。

目前，現(xiàn)有學(xué)者對(duì)電子式電流互感器的積分環(huán)節(jié)提出了較多不同算法，如模擬積分［3］、數(shù)字積分矩形算法［4］、梯形算法［5］、Simpson算法［6］、雙線性變換算法［7-8］和 Al-Alaoui［9］算法等。由于模擬器件受時(shí)漂和溫漂的影響較大，在實(shí)際中很少采用。文獻(xiàn)［5］設(shè)計(jì)了基于STM32F103梯形算法的數(shù)字積分器，但并沒(méi)有與其他積分算法進(jìn)行比較與分析。文獻(xiàn)［6］研究了3種常用的矩形算法、梯形算法、Simpson算法的電子式互感器特性，但是文獻(xiàn)僅涉及了不同積分算法幅頻特性的比較，沒(méi)有深入分析含不同積分算法的互感器頻帶特性，對(duì)其動(dòng)態(tài)響應(yīng)能力的比較沒(méi)有涉及。文獻(xiàn)［7］研究了雙線性變換算法的幅頻特性和諧波特性，但是僅從精度上與其它積分算法進(jìn)行了比較。文獻(xiàn)［10］設(shè)計(jì)了對(duì)Rogowski線圈電流互感器的精確度、飽和特性、頻率特性和暫態(tài)特性的測(cè)試方案，但是文中僅涉及互感器等效電路本身，沒(méi)有考慮數(shù)字積分環(huán)節(jié)對(duì)其傳變特性的影響。綜上所述，對(duì)現(xiàn)有不同積分算法的電子式電流互感器整體的頻帶特性還需要完善，對(duì)其動(dòng)態(tài)響應(yīng)能力的比較還有待深入研究。

本文以Rogowski線圈電子式電流互感器為例，對(duì)常用的矩形、梯形、Simpson和雙線性變換積分算法進(jìn)行比較分析，全面評(píng)估不同積分算法可能對(duì)其傳變特性的影響。首先，分別獲取Rogowski線圈等效電路和積分算法的傳遞函數(shù)；其次，對(duì)各個(gè)環(huán)節(jié)幅頻和相頻特性進(jìn)行分析，獲取并比較了含不同積分算法的電子式互感器整體的頻帶特性；最后，以某實(shí)際110 kV變電站用電子式電流互感器參數(shù)為研究對(duì)象進(jìn)行分析，針對(duì)穩(wěn)態(tài)和信號(hào)突變的動(dòng)態(tài)情況，對(duì)含不同積分算法的電子式互感器傳變特性進(jìn)行仿真。

1 電子式電流互感器的建模

所采用的Rogowski線圈電子式電流互感器的結(jié)構(gòu)如圖1所示。被測(cè)電流經(jīng)Rogowski線圈輸出微分后的信號(hào)，再通過(guò)A/D采樣和數(shù)字積分器輸出與被測(cè)電流成一定比例的信號(hào)，最后通過(guò)合并單元作進(jìn)一步處理。

圖1 Rogowski線圈電子式電流互感器結(jié)構(gòu)Fig.1 Structure of Rogowski coil electronic current transformers

1.1 Rogowski線圈的等效電路

Rogowski線圈的等效電路如圖2所示，其中L0、R0、C0分別為Rogowski線圈的自感、內(nèi)阻和雜散電容，Rs為取樣電阻，M為Rogowski線圈的互感系數(shù)，e（t）為 Rogowski線圈的感應(yīng)電勢(shì)，u0（t）為 Rogowski線圈的輸出電壓，I為被測(cè)電流。

圖2 Rogowski線圈的等效電路圖Fig.2 Equivalent circuit diagram of Rogowski coil

通過(guò)基本電路分析和拉氏變換可得到Rogowski線圈的傳遞函數(shù)：

用于電力系統(tǒng)的Rogowski線圈電子式電流互感器常工作于外積分方式，此時(shí)Rogowski線圈的輸出電壓與被測(cè)電流滿足微分關(guān)系，即：

由上式可知，為了還原被測(cè)電流，就需要對(duì)輸出電壓進(jìn)行積分。

1.2 不同積分算法的傳遞函數(shù)

數(shù)字積分通過(guò)軟件算法來(lái)實(shí)現(xiàn)，且基本不受溫度和其他因素的影響，因而得到廣泛的應(yīng)用。常用的數(shù)字積分算法主要有矩形算法、梯形算法和Simpson算法。根據(jù)它們的z傳遞函數(shù)，依次得到其遞推公式如下。

矩形算法：

梯形算法：

Simpson算法：

式中 T為采樣間隔，n表示采樣時(shí)刻，X（n）與Y（n）分別表示n時(shí)刻積分環(huán)節(jié)的輸入和輸出，下標(biāo)r、t和s分別表示矩形算法、梯形算法和Simpson算法。

為了抑制低頻信號(hào)的干擾，最近文獻(xiàn)提出了雙線性變換積分算法［7-8］。其主要思想是：使s平面的虛軸映射為z平面的單位圓，從而實(shí)現(xiàn)模擬域到數(shù)字域的映射。s與z的關(guān)系公式為：

利用雙線性變換，將公式（6）代入某模擬積分公式中：

并得到對(duì)應(yīng)的遞推公式為：

式中 α=-TRf/（TR+2RRfC），β=-（TR-2RRfC）/（TR+2RRfC），R、C和Rf分別為對(duì)應(yīng)模擬積分的電阻、電容和反饋電阻值。

為了全面分析比較不同積分算法的電子式電流互感器的傳變特性，本文采用式（3）～式（5）和式（8）的數(shù)字積分傳遞函數(shù)對(duì)Rogowski線圈電子式電流互感器進(jìn)行頻帶特性以及不同信號(hào)源下傳變特性的研究。

2 含積分環(huán)節(jié)的Rogowski線圈互感器的頻帶特性

2.1 Rogowski線圈環(huán)節(jié)的頻帶特性

按照前述Rogowski線圈的等效電路和第3節(jié)的仿真參數(shù)可得Rogowski線圈環(huán)節(jié)頻帶特性，其Bode圖如圖3所示。

從相頻特性可以看出，Rogowski線圈的輸出電壓與被測(cè)電流約有90°的相位差，存在微分關(guān)系。

2.2 不同積分算法的頻帶特性

按照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種不同積分算法的對(duì)數(shù)幅頻曲線斜率均為-20 dB/dec，但是它們?cè)诘皖l響應(yīng)上存在差異，特別是雙線性變換積分算法由于是根據(jù)模擬積分變換而得，因此它具有與模擬積分相類似的頻帶特性，即在頻率f＜1 Hz區(qū)域?qū)?shù)幅頻曲線斜率為0 dB/dec，而常用的數(shù)字積分算法對(duì)數(shù)幅頻曲線斜率為-20 dB/dec。

另外，從上述相頻特性可以看出，這4種不同積分算法都達(dá)到了-90°的相位延遲，但是它們?cè)诘皖l響應(yīng)和高頻區(qū)域上仍然存在差異。在低頻區(qū)域，采用的雙線性變換積分算法的相位延遲幾乎為0°，表示在低頻區(qū)域，該算法喪失積分功能；在高頻區(qū)域，矩形算法的相位則隨著頻率的增大而增大。

因此，從上述Bode圖可以看出，不同數(shù)字積分在低頻和高頻區(qū)域存在一定的差異，其中雙線性變換積分算法反映在低頻方面的差異；而矩形算法反映在高頻區(qū)域。

2.3 含積分環(huán)節(jié)的電子式互感器整體頻帶特性

基于上述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)為水平直線，可見(jiàn)在上述頻帶范圍內(nèi)，不同積分算法時(shí)電子式電流互感器輸出信號(hào)的幅值和相位不受頻率變化的影響，即能夠準(zhǔn)確地還原被測(cè)電流中的各種頻率分量，頻帶特性較好。然而，當(dāng)頻率小于10 Hz時(shí)，相對(duì)其他積分算法，采用雙線性變換積分算法時(shí)的還原特性會(huì)更容易受到原始信號(hào)低頻分量的影響；而當(dāng)頻率大于10 kHz時(shí)在高頻區(qū)域，上述4種積分算法的互感器輸出的相位則會(huì)隨著頻率的增大而增大，從而很大程度上影響其還原特性。

3 仿真分析

為了較全面的分析不同積分算法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Ω，采樣時(shí)間 T=0.1 ms，其中當(dāng)采用雙線性變換積分算法時(shí)取R=30 kΩ，C=0.1μF，Rf=1 MΩ。基于PLECS仿真平臺(tái)搭建了含上述4種不同數(shù)字積分算法的Rogowski線圈電子式電流互感器的模型，下面針對(duì)穩(wěn)態(tài)信號(hào)、積分初值不確定，以及信號(hào)突變或富含直流、諧波信號(hào)不同情況，對(duì)不同積分算法的互感器傳變特性進(jìn)行仿真。

3.1 工頻穩(wěn)態(tài)下的傳變特性比較

現(xiàn)選擇初始相位角為0°，頻率為50 Hz，標(biāo)幺值為1的額定電流 i（t）=sin（100πt）作為輸入信號(hào)，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)良好，均能較好的還原被測(cè)波形。其比差和相差如表1所示，通過(guò)表中計(jì)算結(jié)果可知，4種積分算法的電流互感器在額定電流下都能滿足較高的準(zhǔn)確度要求。

表1 額定電流下互感器的準(zhǔn)確度Tab.1 Accuracy of transformers at rated current

3.2 積分初值不確定情況下的動(dòng)態(tài)響應(yīng)特性比較

由于電子式電流互感器還原被測(cè)電流需要對(duì)Rogowski線圈的輸出電壓進(jìn)行積分處理，即：

因此需要確定數(shù)字積分的初值。為了驗(yàn)證不同積分初值對(duì)上述積分算法的魯棒性，現(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可以看出，采用雙線性變換積分算法時(shí)互感器可以在較短時(shí)間內(nèi)的還原輸入電流，具有積分初值不確定下魯棒性好的優(yōu)點(diǎn)。而采用矩形算法、梯形算法和Simpson算法的互感器輸出則受積分初值不確定下影響大，輸出信號(hào)含有50%幅值的直流分量，導(dǎo)致其傳變特性輸出信號(hào)幅值下移50%，還原信號(hào)能力明顯降低。

3.3 輸入信號(hào)幅值突變下的動(dòng)態(tài)響應(yīng)特性比較

由于電網(wǎng)中的功率控制或負(fù)荷變化等各種因素下可使互感器的輸入電流突然增加或者降低，為了仿真分析電流幅值突變時(shí)對(duì)不同積分算法下互感器傳變特性的動(dòng)態(tài)跟蹤能力比較，現(xiàn)假定輸入電流為i（t）=sin（100πt），在 t=0.044 s時(shí)刻電流幅值從 1 p.u.突變?yōu)?.5 p.u.，在t=0.104 s時(shí)刻電流幅值又從1.5 p.u.突變回1 p.u.，仿真比較4種不同積分算法時(shí)互感器動(dòng)態(tài)響應(yīng)情況，如圖8所示。

圖8 輸入信號(hào)幅值突變下互感器傳變特性比較Fig.8 Comparisons of the transducer transfer characteristicswhen the input amplitude is changed

從圖8中的仿真波形可以看出，4種不同積分算法的互感器在電流幅值突變情況下輸出波形與輸入波形吻合，能較好地還原輸入信號(hào)。但是從其動(dòng)態(tài)響應(yīng)能力可以看出（觀察放大波形），4種不同積分算法的互感器動(dòng)態(tài)響應(yīng)和幅值跟蹤方面仍存在一定的差異，與其余3中積分算法比較，Simpson算法的幅值突變的跟蹤響應(yīng)能力較差。

3.4 輸入信號(hào)富含直流分量時(shí)動(dòng)態(tài)響應(yīng)特性比較

當(dāng)線路中發(fā)生過(guò)電壓或故障時(shí)，由于電感或電容的存在，往往使得電流不能突變而導(dǎo)致線路中存在直流分量的情況。此外，由于Rogowski線圈輸出信號(hào)弱，在外界干擾下也可能導(dǎo)致積分環(huán)節(jié)的輸入信號(hào)中含有直流分量。為了比較分析富含直流分量時(shí)互感器傳變特性，下面對(duì)Rogowski線圈輸入信號(hào)和積分環(huán)節(jié)輸入信號(hào)富含直流分量時(shí)的情況進(jìn)行仿真，如圖9所示。其中，圖9（a）為輸入電流i（t）=sin（100πt），在 t=0.1 s開(kāi)始含有 1 p.u.衰減時(shí)間常數(shù)T=0.05 s的直流分量時(shí)傳感器傳變特性的仿真結(jié)果；圖9（b）為假定在t=0.1 s時(shí)積分環(huán)節(jié)輸入信號(hào)含10%幅值的穩(wěn)定直流量干擾下傳感器傳變特性的仿真結(jié)果。

圖9 含不同直流分量時(shí)互感器傳變特性比較Fig.9 Comparisons of the transducer transfer characteristicswhen the input has different DC signals

從圖9（a）可以看出，當(dāng)Rogowski線圈輸入含衰減直流分量時(shí)，不同積分算法基本上都能還原其輸入信號(hào)的變化，但相比之下，采用雙線性變換積分算法時(shí)還原能力相對(duì)較差。

從圖9（b）可以看出，當(dāng)積分環(huán)節(jié)輸入含穩(wěn)定直流分量干擾情況下，雙線性變換積分算法受直流量干擾的影響相對(duì)較小，還原特性好，而其余3種算法受直流量干擾的影響較大，導(dǎo)致其互感器輸出信號(hào)隨著時(shí)間推移不斷增大，呈現(xiàn)明顯發(fā)散的傳變特性。

3.5 輸入信號(hào)富含諧波分量時(shí)動(dòng)態(tài)響應(yīng)特性比較

為了模擬線路中存在諧波分量時(shí)互感器的動(dòng)態(tài)性能，現(xiàn)假定線路中存在2～13次的諧波分量，即輸入電流為 i（t）=sin（100πt），并在 t=0.1 s開(kāi)始輸入幅值分別為 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種積分算法的互感器均能很好的還原被測(cè)電流中的低次諧波分量。

圖10 含2-13次諧波分量時(shí)互感器傳變特性比較Fig.10 Comparisons of the transducer transfer characteristicswhen the input has 2-13 order harmonic components

另外，當(dāng)線路中電感、電容存在導(dǎo)致諧振可能會(huì)存在高次諧波情況，為了分析不同積分算法下互感器對(duì)線路中高次諧波分量的影響，現(xiàn)假定輸入電流為 i（t）=sin（100πt），在 t=0.1 s開(kāi)始輸入頻率為106Hz、幅值為0.1的高次諧波，采用不同積分算法時(shí)其互感器傳變特性如圖11所示。

從圖中可以看出，采用上述4種積分算法的互感器還原高次諧波信號(hào)的能力都較差，這與圖5中含不同積分算法時(shí)互感器整體輸出的頻帶特性分析結(jié)果一致，即當(dāng)Rogowski線圈輸入信號(hào)含高次諧波分量時(shí)，上述不同積分算法并不能很好反映其輸出高次諧波特性。

圖11 含高次諧波分量時(shí)互感器傳變特性比較Fig.11 Comparisons of the transducer transfer characteristicswhen the input has a high-order harmonic component

4 結(jié)束語(yǔ)

以Rogowski線圈電子式電流互感器為例，通過(guò)建立傳遞函數(shù)和仿真模型，對(duì)常用的矩形、梯形、Simpson和雙線性變換積分算法下的互感器頻帶特性和傳變特性時(shí)域仿真進(jìn)行了比較分析，得到以下結(jié)論：

（1）從頻帶特性上和輸入信號(hào)富含諧波分量時(shí)的仿真比較可以看出，在頻率10 Hz～10 kHz的范圍內(nèi)，上述4種積分算法的電子式電流互感器幅值和相位不受頻率變化的影響，即能夠較準(zhǔn)確地還原被測(cè)電流中的各種頻率分量，頻帶特性較好，但在低頻以及高頻率信號(hào)下還原能力較差；

（2）從穩(wěn)態(tài)特性的比較分析可以看出，4種積分算法的互感器還原輸入電流的能力都較好，在額定電流下都能滿足較高的準(zhǔn)確度要求。但在積分初值不確定的情況下，采用雙線性變換積分算法的互感器傳變特性較好，魯棒性較高；

（3）從不同輸入信號(hào)突變的動(dòng)態(tài)跟蹤能力的比較分析可以看出，當(dāng)輸入信號(hào)幅值突變時(shí)，不同積分算法的互感器能較好地還原特性，相比Simpson算法存在較大的誤差；當(dāng)輸入信號(hào)含直流分量時(shí)，采用雙線性變換積分算法的互感器傳變特性較好，尤其是在積分環(huán)節(jié)輸入有直流分量情況。

從上述結(jié)論可以看出，采用不同積分算法時(shí)電子式電流互感器的頻帶特性和穩(wěn)態(tài)傳變特性較好，并沒(méi)有顯著的差別；但在動(dòng)態(tài)特性方面，尤其在積分初值不確定、輸入信號(hào)富含直流分量和高頻諧波分量時(shí)，能較好地適應(yīng)其抗干擾能力和動(dòng)態(tài)跟蹤能力方面的數(shù)字積分算法還有待進(jìn)一步研究。