張仁武，帥定新

（1.攀枝花學(xué)院交通與汽車工程學(xué)院，攀枝花 617000；2.攀枝花學(xué)院電氣信息工程學(xué)院，攀枝花 617000）

諧波諧振的s域模態(tài)分析法

張仁武1，帥定新2

（1.攀枝花學(xué)院交通與汽車工程學(xué)院，攀枝花 617000；2.攀枝花學(xué)院電氣信息工程學(xué)院，攀枝花 617000）

為了減少頻率掃描的計算時間，提高計算精度，提出了利用s域模態(tài)分析確定諧波諧振頻率的方法。通過牛頓-拉夫遜迭代法對電力系統(tǒng)在s域中傳遞函數(shù)極點的求解，可得到電力系統(tǒng)中諧振頻率及大小程度等信息。在每次迭代求解過程中，只需更新復(fù)頻率s值即能確定諧振的幅度等信息，不需要頻率掃描，從而節(jié)省計算時間并提高計算精度。測試系統(tǒng)分析表明，s域模態(tài)分析法對于諧波諧振的分析簡單有效，快速準確，是分析諧波諧振及諧振治理的有效工具。

諧波諧振；牛頓-拉夫遜迭代法；復(fù)頻率；測試；s域模態(tài)分析

隨著電力電子技術(shù)的快速發(fā)展，大量具有非線性特性的電力設(shè)備如電弧爐、電力機車、變頻設(shè)備等投入運行，電力系統(tǒng)中的諧波污染、諧振問題也隨之越來越嚴重。諧波諧振是電力系統(tǒng)諧波引起的危害之一，易造成因諧波過電壓引發(fā)的設(shè)備故障問題[1-2]。諧波諧振本質(zhì)是由系統(tǒng)中的容性元件和感性元件之間頻繁的無功功率交換造成。

目前在工業(yè)電力系統(tǒng)，如風電系統(tǒng)、電氣化鐵道、煤礦電力系統(tǒng)等均發(fā)生諧波諧振現(xiàn)象[3-5]，但鮮有文獻能提出對該現(xiàn)象有效分析并治理的工具[6-9]。目前，研究諧波諧振的方法有3種：①頻譜分析，頻譜分析對驅(qū)動點阻抗進行頻率掃描可識別諧振存在及確定諧振頻率，但頻譜分析不能提供更多有效的諧振信息[10]；②模態(tài)分析，模態(tài)分析可識別諧振存在并確定諧振頻率，進而準確定位出諧波諧振區(qū)域、元件敏感度等[11-12]，在電力系統(tǒng)中得到廣泛應(yīng)用；③s域模態(tài)分析，s域模態(tài)分析不僅能得到全面的諧波諧振信息，而且計算精度與計算時間都優(yōu)于模態(tài)分析。

然而，頻譜分析和模態(tài)分析均需對頻率進行掃描，理論上來說當頻率精度提高時，將成倍增加計算量。對于模態(tài)分析而言，每次需對節(jié)點導(dǎo)納矩陣進行更新及特征值求解，嚴重限制了模態(tài)分析在分析大型電力網(wǎng)絡(luò)及頻率精度高等場景的應(yīng)用。節(jié)點導(dǎo)納矩陣在s域中計算時，不需要對每個頻率點進行掃描，只需對趨于奇異性的節(jié)點導(dǎo)納矩陣進行迭代計算，求解諧振頻率精度高、計算速度快。

因此，本文將s域分析法應(yīng)用于電力系統(tǒng)的諧波諧振分析，通過建立s域下的系統(tǒng)節(jié)點導(dǎo)納矩陣，利用牛頓-拉夫遜N-R（Newton-Raphson）迭代法進行極點/零點的求取，從而確定系統(tǒng)的并聯(lián)/串聯(lián)諧振頻率及諧振幅度。

1 電力系統(tǒng)s域模型

1.1 s域的節(jié)點導(dǎo)納矩陣

對某個電力網(wǎng)絡(luò)，可建立關(guān)系式為

式中：Y（s）為電力系統(tǒng)的s域節(jié)點導(dǎo)納矩陣；x(s)為廣義狀態(tài)空間；u（s）為系統(tǒng)的輸入；y（s）為系統(tǒng)的輸出；B和C分別為與u（s）和x(s)維度匹配的系數(shù)矩陣。

為便于分析網(wǎng)絡(luò)對于某個節(jié)點注入電流量的響應(yīng)，可通過式（1）和式（2），得到特定的單輸入單輸出系統(tǒng)的關(guān)系式

式中：b為行向量，其第j個元素為1，其余元素為0；c為行向量，其第i個元素為1，其余元素為0；uj表示系統(tǒng)在第j個節(jié)點的注入電流。當在系統(tǒng)的節(jié)點j注入電流時，可得節(jié)點i的輸出電壓，則輸入阻抗為

輸入阻抗Zij（s）代表節(jié)點i對節(jié)點j的互阻抗，其物理意義是節(jié)點i對節(jié)點j的影響。

1.2 s域模型以及關(guān)于s的導(dǎo)數(shù)

節(jié)點導(dǎo)納矩陣Y（s）是對電力系統(tǒng)的結(jié)構(gòu)及參數(shù)特性的描述，電力系統(tǒng)中的元件可分為RLC元件、傳輸線路、變壓器等元件。

系統(tǒng)的大部分元件可等效為如圖1所示的RLC串聯(lián)或并聯(lián)的形式。通過在時域中的拉普拉斯變換，可得出電阻、電感、電容在s域中的模型。

圖1 RLC支路模型Fig.1 Model of RLC branch

串聯(lián)RLC支路以及并聯(lián)RLC支路的復(fù)頻域?qū)Ъ{模型為

式（5）和式（6）關(guān)于s的一階導(dǎo)數(shù)分別為

對于分布參數(shù)傳輸線路模型而言，可等效為圖2的π模型線路，π模型的阻抗和導(dǎo)納分別為

式中：l為單位傳輸線路長度；γ為傳輸線路傳播常數(shù)；Yc為傳輸線路特征導(dǎo)納。Ys應(yīng)被加入到Y(jié)（s）對角線元素中；而Yp的負值應(yīng)被加入到Y(jié)（s）的非對角線元素中。傳輸線路的導(dǎo)納函數(shù)是關(guān)于單位傳輸線路長度l、傳輸線路傳播常數(shù)γ、傳輸線路特征導(dǎo)納Yc的函數(shù)。在單位長度下的傳輸線路傳播常數(shù)γ以及傳輸線路特征導(dǎo)納Yc分別為為單位長度的傳輸線路阻抗，Yu為單位長度的傳輸線路導(dǎo)納。

圖2 傳輸線路等效π模型Fig.2 Equivalent π model of transmission line

Ys、Yp中線路導(dǎo)納Yc關(guān)于s的一階導(dǎo)數(shù)為

2 電力系統(tǒng)傳遞函數(shù)極點計算

由定義可知系統(tǒng)的極點與系統(tǒng)輸入輸出無關(guān)。故令系統(tǒng)輸出yi=1，系統(tǒng)輸入uj=0，由式（3）得系統(tǒng)傳遞函數(shù)為

式中，G(s)為系統(tǒng)的二階傳遞函數(shù)。

將式（13）運用N-R迭代法進行迭代，可得

式中：λ為系統(tǒng)傳遞函數(shù)G（s）的極點；k為迭代次數(shù)；Δλ為λ的增量，即

將系統(tǒng)傳遞函數(shù)G（s）及其導(dǎo)數(shù)代入式（14）中，可得迭代增量Δλ為

通過式（1）、式（2）表示的簡單線性系統(tǒng)，可求出式（16）中的cTY-1與Y-1b分別為

將式（17）、式（18）代入式（16）中經(jīng)運算，可得

式中：λj為已找到的極點；Rj為對應(yīng)的留數(shù)（j為已找到極點的指數(shù)）。

若所求解的極點為重極點或者多重極點，則式（21）可表達為

需要注意到的是在N-R迭代法過程中，當Δλ小于某個給定誤差（如10-6）時，則停止迭代，否則更新λ，繼續(xù)迭代，直到滿足迭代精度要求為止。迭代初值的選取對迭代的快慢以及迭代的收斂有著重要的關(guān)系，根據(jù)文獻[13]提出的方法，可消除初值對迭代的影響。為消除已求出極點對迭代的影響，在每一次迭代過程中，可先將已求出的極點濾除，再繼續(xù)迭代求解，即

式中：m為多重極點個數(shù)，m＞1；Rm為傳遞函數(shù)的留數(shù)。

再通過式（16）～式（19）的計算，可以得出修正后的迭代增量Δλ為

傳遞函數(shù)的留數(shù)，根據(jù)定義，可得

將式（23）迭代求出的值代入式（24）中，有

化簡為

由式（23）可求出傳遞函數(shù)的留數(shù)，即諧振幅度為

3 算例驗證

3.1 三母線測試系統(tǒng)

圖3 三母線測試系統(tǒng)Fig.3 Three-bus test system

圖4 測試系統(tǒng)頻譜分析結(jié)果Fig.4 Spectrum analysis results of the test system

圖5 三母線測試系統(tǒng)模態(tài)阻抗Fig.5 Modal impedances of the three-bus test system

表1 測試系統(tǒng)模態(tài)分析結(jié)果Tab.1 Modal analysis results of the test system

對圖3的三母線測試系統(tǒng)分別使用文中所提到的3種方法進行諧波諧振分析。系統(tǒng)中的3條母線可被激勵產(chǎn)生或觀察到諧波諧振。通過對母線2注入1 p.u.的電流，得到如圖4的頻譜分析。雖然頻譜分析可得出測試系統(tǒng)的諧振頻率，但難以確定每條母線所對應(yīng)的諧振頻率，也即無法定位具體諧振位置。對測試系統(tǒng)進行模態(tài)分析可得出圖5及表1，通過模態(tài)分析可知母線1的諧振頻率為51.8 p.u.，母線3的諧振頻率為9.8 p.u.。最后對測試系統(tǒng)進行s域模態(tài)分析，通過N-R迭代計算可得測試系統(tǒng)傳遞函數(shù)極點值，對極點值進行處理（虛數(shù)部分除以2π）可得測試系統(tǒng)諧振頻率及諧振位置，如表2所示。在N-R迭代求解極點計算過程中，極點所對應(yīng)的留數(shù)也可一同求出，留數(shù)可以表征諧振幅度，這說明s域模態(tài)分析可確定更多的諧振信息。雖然模態(tài)分析可得出測試系統(tǒng)的諧振頻率和諧振位置，但由于模態(tài)分析需對頻率一一掃描，耗時較長，共耗時0.14 s，而s域模態(tài)分析共用時0.001 6 s，耗時明顯少于模態(tài)分析所用時間。測試系統(tǒng)只含有三母線，當分析多母線系統(tǒng)時，s域模態(tài)分析在計算時間上的優(yōu)勢會更加明顯。

表2 三母線測試系統(tǒng)極點計算Tab.2 Calculation of poles in the three-bus test system

3.2 測試系統(tǒng)

為了驗證s域模態(tài)分析方法的性能，本文用如下的測試系統(tǒng)進行s域模態(tài)分析。圖6為一典型的工業(yè)電力系統(tǒng)，基本的網(wǎng)路結(jié)構(gòu)和線路參數(shù)均摘自IEEE Standard 519—1992[14]。電力系統(tǒng)中每一個輸出負載都并聯(lián)了功率補償電容，而電容會引起諧波諧振。圖6所示的電力系統(tǒng)，可由文中前述的建模方法建立等效電路。

圖6 五母線測試系統(tǒng)Fig.6 Five-bus test system

系統(tǒng)采用工頻50 Hz，系統(tǒng)中各元件的參數(shù)在表3中給出，分別為變壓器T1等效電感和電阻，L1以及R1分別為輸出負載端的等效電感和電阻。

表3 測試系統(tǒng)參數(shù)Tab.3 Parameters of the test system

由模態(tài)分析法，可得出測試系統(tǒng)在頻率范圍1～250 Hz內(nèi)有4個諧振峰值，如圖7所示，母線1的諧振頻率為143 p.u.，母線2的諧振頻率為16 p.u.，母線3的諧振頻率為65 p.u.，母線4的諧振頻率為172 p.u.。由s域模態(tài)分析，得出表4及圖8。圖8是由測試系統(tǒng)傳遞函數(shù)極點值畫出，通過圖8與圖7的對比，可得出s域模態(tài)分析得出的諧振頻率與模態(tài)分析得出的結(jié)果比例幾乎一致，圖8還可揭示出諧振幅度。模態(tài)分析所用時間為0.094 s，而s域共用時間為0.015 s，耗時明顯減少。

圖7 測試系統(tǒng)模態(tài)阻抗Fig.7 Modal impedances of the test system

表4 測試系統(tǒng)極點計算Tab.4 Calculation of poles in the test system

圖8 測試系統(tǒng)傳遞函數(shù)極點Fig.8 Poles of the test system’s transfer function

4 結(jié) 語

通過在s域?qū)﹄娏ο到y(tǒng)建立節(jié)點導(dǎo)納矩陣，計算傳遞函數(shù)極點求出電力系統(tǒng)的諧振頻率并與模態(tài)分析法進行對比。結(jié)果表明s域的求解方法不僅能夠揭示出電力系統(tǒng)諧振信息，并且計算時間明顯少于現(xiàn)有的其他方法，這是由于s域模態(tài)分析不需要對頻率掃描，而是直接對電力系統(tǒng)傳遞函數(shù)的極點進行迭代計算，從而得出電力系統(tǒng)的諧振頻率及幅度信息。s域模態(tài)分析可用于任何大小、任何拓撲的電力系統(tǒng)諧振分析，在分析多節(jié)點的電力系統(tǒng)時，計算時間會大幅度減少，s域模態(tài)分析是分析諧波諧振及諧振治理的有效工具。

[1]谷毅，趙玉柱，張國威，等（Gu Yi，Zhao Yuzhu，Zhang Guowei，et al）.關(guān)于500 kV東明開關(guān)站啟動調(diào)試期間發(fā)生電壓諧振的分析（Discussion on voltage resonance occurred in Bao-Ren 5231 transmission line during startup of Dongming switching station）[J].電網(wǎng)技術(shù)（Power System Technology），2002，26（12）：71-74．

[2]趙書強，馬燕峰，賀春（Zhao Shuqiang，Ma Yanfeng，He Chun）.抑制諧波的配電網(wǎng)無功優(yōu)化規(guī)劃（Reactive power optimal planning of distribution network with harmonic mitigation）[J].電網(wǎng)技術(shù)（Power System Technology），2004，28（6）：71-75．

[3]王麗芬（Wang Lifen）.礦井供電系統(tǒng)中諧波與諧振過電壓現(xiàn)象分析（Harmonics and harmonic over-voltage analysis in coal mine supply system）[J].煤（Coal），2010，19（7）：102-103．

[4]何正友，胡海濤，方雷，等（He Zhengyou，Hu Haitao，F(xiàn)ang Lei，et al）.高速鐵路牽引供電系統(tǒng)諧波及其傳輸特性研究（Research on the harmonic in high-speed railway traction power supply system and its transmission characteristic）[J].中國電機工程學(xué)報（Proceedings of the CSEE），2011，31（16）：55-62．

[5]沈啟平，姚駿（Shen Qiping，Yao Jun）.電網(wǎng)諧波條件下雙饋感應(yīng)風電系統(tǒng)改進控制（Improved control of doubly fed induction generator wind turbine with series gridside converter under distorted grid voltage conditions）[J].電力系統(tǒng)及其自動化學(xué)報（Proceedings of the CSU-EPSA），2015，27（5）：7-13．

[6]葛少云，董佳霖，劉福潮，等（Ge Shaoyun，Dong Jialin，Liu Fuchao，et al）.配電網(wǎng)諧波治理降損效果實證方法（Empirical method of harmonic control loss reduction effect in distribution network）[J].電力系統(tǒng)及其自動化學(xué)報（Proceedings of the CSU-EPSA），2015，27（3）：1-4．

[7]張杰成，喬鳴忠，朱鵬，等（Zhang Jiecheng，Qiao Mingzhong，Zhu Peng，et al）.滑窗迭代DFT檢測諧波和無功電流的新算法（New method of harmonic and reactive current detection based on sliding-window iteration DFT）[J].電力系統(tǒng)及其自動化學(xué)報（Proceedings of the CSUEPSA），2015，27（12）：18-22．

[8]穆子龍，張雨晴（Mu Zilong，Zhang Yuqing）.HVDC系統(tǒng)送端孤島下諧波不穩(wěn)定對策探討（Research on countermeasures to harmonic instability problem at sending end of HVDC transmission system in isolated island operation）[J].電力系統(tǒng)及其自動化學(xué)報（Proceedings of the CSUEPSA），2015，27（4）：33-38．

[9]劉軍，黃純，戴栩生（Liu Jun，Huang Chun，Dai Xusheng）.間諧波檢測的快速自適應(yīng)離散廣義S變換方法（Inter-harmonics detection based on fast adaptive discrete generalized S-transform）[J].電力系統(tǒng)及其自動化學(xué)報（Proceedings of the CSU-EPSA），2015，27（7）：13-17，29.

[10]Huang Zhenyu，Xu Wilsun，Dinavahi V R.A practical harmonic resonance guideline for shunt capacitor applications[J].IEEE Trans on Power Delivery，2003，18（4）：1382-1387．

[11]Huang Zhenyu，Cui Yu，Xu Wilsun.Application of modal sensitivity for power system harmonic resonance analysis[J].IEEE Trans on Power Systems，2007，22（1）：222-231.

[12]胡海濤，何正友，錢澄浩，等（Hu Haitao，He Zhengyou，Qian Chenghao，et al）.基于模態(tài)分析的全并聯(lián)AT網(wǎng)動態(tài)諧波諧振研究（Modal analysis based research on dynamic harmonic resonance of all-parallel autotransformer traction supply network）[J].電網(wǎng)技術(shù)（Power System Technology），2012，36（1）：163-169．

[13]Gomes Sergio Jr，Martins Nelson，Portela Carlos.Sequential computation of transfer function dominant poles of sdomain system models[J].IEEE Trans on Power Systems，2009，24（2）：776-784．

[14]IEEE STD 519—1992，IEEE recommended practice and requirements for harmonic control in electric power systems[S].

s-Domain Modal Analysis Method for Harmonic Resonance

ZHANG Renwu1，SHUAI Dingxin2
（1.School of Transportation and Automotive Engineering，Panzhihua University，Panzhihua 617000，China；2.School of Electrical and Information Engineering，Panzhihua University，Panzhihua 617000，China）

To reduce the calculation time and improve the calculation accuracy of frequency scanning，a method of determining the harmonic resonance frequency is proposed in this paper by usings-domain modal analysis.By using Newton-Raphson method to seek the poles of the power system’s transfer function ins-domain，the proposed method can obtain the resonance frequency and its magnitude of the power system.In each iteration of the solution process，only thesdata of complex frequency needs to be updated to determine the information of resonance（e.g.，amplitude）without scanning the frequency，which saves the calculation time and improves the calculation accuracy.The analysis of a test system shows thats-domain modal analysis is simple，effective，fast and accurate for harmonic resonance，indicating that it is an effective tool for the analysis of harmonic resonance and resonance control.

harmonic resonance；Newton-Raphson method；complex frequency；test；s-domain modal analysis

TM714

A

1003-8930（2017）10-0113-05

10.3969/j.issn.1003-8930.2017.10.019

2016-01-20；

2017-08-02

張仁武（1965—），男，本科，副教授，研究方向為電力系統(tǒng)分析。Email：971903609@qq.com

帥定新（1979—），男，博士，副教授，研究方向為電力電子與電力傳動。Email：31775897@qq.com