陳 劍 周宇植

?

基于PMU的分布式次同步振蕩在線辨識方法

陳 劍1周宇植2

（1. 國家電網(wǎng)江蘇省電力公司生產(chǎn)技能培訓中心，江蘇蘇州 215004；2. 國家電網(wǎng)江蘇省電力公司蘇州供電公司，江蘇蘇州 215004）

為減輕WAMS主站端的軟硬件開銷，同時也為了避免WAMS主站端對PMU數(shù)據(jù)進行采樣時產(chǎn)生的頻率混疊現(xiàn)象，本文在分析了WAMS/PMU通信架構(gòu)和同步相量測量的算法原理的基礎(chǔ)上，提出了在PMU裝置上實現(xiàn)分布式次同步振蕩在線參數(shù)的辨識方法。本文深入分析和比較了兩種不同的參數(shù)辨識方法，基于相量數(shù)據(jù)的次同步振蕩參數(shù)辨識方法和基于原始采樣值的次同步振蕩參數(shù)辨識方法，討論了基于相量數(shù)據(jù)的次同步振蕩參數(shù)辨識方法的局限性，指出了基于原始采樣值的次同步振蕩參數(shù)辨識方法，既可避免對次同步振蕩分量的誤判，提高次同步振蕩分析監(jiān)測的準確性和適應(yīng)性，同時又更易于實現(xiàn)超同步振蕩等擴展性應(yīng)用需求。

次同步振蕩；分布式測量；FFT；同步相量測量單元

我國經(jīng)濟的不斷發(fā)展促進了電網(wǎng)規(guī)模的不斷擴大，大容量、遠距離輸電是我國電力系統(tǒng)發(fā)展的客觀需求。為提高線路輸送能力、控制并行線路之間的功率分配、增強電力系統(tǒng)暫態(tài)穩(wěn)定性，電力系統(tǒng)中安裝了串聯(lián)補償裝置、高壓直流（HVDC）控制系統(tǒng)、動態(tài)無功補償裝置等設(shè)備。然而，這些設(shè)備和技術(shù)的發(fā)展也帶來了次同步振蕩問題[1]。電網(wǎng)次同步振蕩對發(fā)電機組和系統(tǒng)穩(wěn)定危害很大，當次同步振蕩頻率與發(fā)電機組扭振的固有頻率滿足特定互補關(guān)系時可引發(fā)發(fā)電機組軸系扭振迅速增大甚至機組軸系損壞事故。及時發(fā)現(xiàn)并采取措施抑制系統(tǒng)中的次同步振蕩[2-4]，是保障電力系統(tǒng)安全、穩(wěn)定運行亟待解決的關(guān)鍵問題之一。

基于同步相量測量技術(shù)的廣域測量系統(tǒng)（wide area measurement system, WAMS）在電網(wǎng)安全穩(wěn)定運行中扮演著越來越重要的角色，它在動態(tài)監(jiān)測、動態(tài)安全分析和控制、狀態(tài)估計等方面取得了廣泛的應(yīng)用[5]，同步相量測量單元（phasor measure unit，PMU）作為WAMS子站系統(tǒng)的核心擔負著原始數(shù)據(jù)采集、測量和傳輸?shù)娜蝿?wù)。尤其是近年來，隨著低頻振蕩、次同步振蕩等現(xiàn)象的增加，PMU為各網(wǎng)省級調(diào)度中心在基于WAMS系統(tǒng)基礎(chǔ)上分析各振蕩現(xiàn)象時提供了高精度、高密度原始數(shù)據(jù)源，因此日益取得了廣泛重視。但隨著電網(wǎng)規(guī)模的不斷擴大，PMU布點數(shù)目的不斷增多，使得網(wǎng)絡(luò)結(jié)構(gòu)日益復雜，一方面在WAMS主站端進行低頻振蕩、次同步振蕩分析將會大大增加WAMS主站端的軟硬件開銷，另一方面由于WAMS主站對PMU數(shù)據(jù)采樣存在的頻率混疊現(xiàn)象，使得在WAMS主站分析次同步現(xiàn)象存在一定的弊端[6]。

因此，本文通過對PMU采樣的原始電流數(shù)據(jù)進行頻譜分析，從而準確計算出次同步振蕩的頻率成分，實現(xiàn)了次同步振蕩參數(shù)的辨識與告警。并與文獻[6]提出的“對PMU測量的相量數(shù)據(jù)進行頻譜分析的方法”進行了對比，尤其深入分析了在幅值調(diào)制的動態(tài)條件下兩種辨識方法的區(qū)別，指出了通過對原始采樣數(shù)據(jù)進行分析，既可避免對次同步振蕩分量的誤判，提高次同步振蕩分析監(jiān)測的準確性和適應(yīng)性，同時又更易于實現(xiàn)超同步等擴展性應(yīng)用需求。

1 WAMS/PMU系統(tǒng)構(gòu)架

文獻[7]闡述了國內(nèi)的WMAS/PMU系統(tǒng)采用通信網(wǎng)絡(luò)結(jié)構(gòu)，其中，PMU子站采用TCP/IP協(xié)議通過電力數(shù)據(jù)調(diào)度網(wǎng)按照WAMS主站下發(fā)的傳輸頻率向WAMS主站以恒定頻率實時發(fā)送動態(tài)數(shù)據(jù)，其網(wǎng)絡(luò)結(jié)構(gòu)如圖1所示。

圖1 WAMS/PMU系統(tǒng)構(gòu)架

在該架構(gòu)下，WAMS主站對PMU的數(shù)據(jù)也是一個數(shù)據(jù)采樣頻譜分析的過程，其采樣率即為PMU上送同步相量的速率。IEEE Std C37.118—2011規(guī) 定[8]，其傳輸速率最高為100Hz。但實際工程應(yīng)用中，國內(nèi)大多采用了50Hz甚至25Hz的傳輸速率。以25Hz為例，根據(jù)香農(nóng)采樣定律，WAMS主站在取得PMU數(shù)據(jù)后，只能分析出不超過12.5Hz的次同步振蕩分量。因此，這極大限制了在WAMS主站分析次同步振蕩現(xiàn)象的應(yīng)用前景。

2 次同步振蕩對PMU的測量影響

電力系統(tǒng)次同步振蕩是電力系統(tǒng)中的電氣系統(tǒng)與發(fā)電機組以低于同步頻率的某個或多個自然振蕩頻率交換能量的一種不正常運行狀態(tài)。電力系統(tǒng)中安裝的串聯(lián)補償裝置、高壓直流（HVDC）控制系統(tǒng)、動態(tài)無功補償裝置等設(shè)備均有可能引起次同步振蕩問題[9]。文獻[6]指出，無論是何種產(chǎn)生機理，對于PMU裝置而言，產(chǎn)生次同步振蕩時均表現(xiàn)為所測量的電流信號中疊加了次同步振蕩分量。其數(shù)學模型為

式中，0為基波頻率；為各次同步振蕩分量的頻率；1為基波信號初相角；為各次同步振蕩分量的初相角；為基波的幅值；DA為各次同步振蕩分量幅值。頻率范圍通常在低頻振蕩頻率（0.2～2.5Hz）和同步頻率（50Hz）之間。

文獻[10]中，作者對使用DFT方法計算相量的方法做了詳細的推導，文中指出對于含有諧波輸入的信號為

其相量計算結(jié)果中會存在0+和0的頻譜分量，其中0頻譜分量的幅值最大。由此，可以推導出，對于式（1）所示的含有次同步分量的輸入信號而言，其PMU相量結(jié)果中將存在0+和0的頻譜分量，這為使用PMU相量結(jié)果進行FFT分析，從而得到次同步分量的頻率和幅值奠定了理論基礎(chǔ)。

3 基于相量數(shù)據(jù)的次同步振蕩參數(shù)辨識

3.1 辨識流程

文獻[6]基于相量數(shù)據(jù)實現(xiàn)了次同步振蕩參數(shù)的辨識和告警，其辨識流程如圖2所示。

圖2 次同步振蕩分析辨識流程

本文對文獻[6]提出的“基于相量數(shù)據(jù)分析次同步振蕩參數(shù)”的方法進行了如下Matlab仿真。其中根據(jù)FFT基本理論，F(xiàn)FT的頻率分辨率為

本文中，選取s=100Hz、=256，則頻率分辨率可達0.390625Hz。圖3是對含有幅值為0.1，頻率為15Hz的次同步分量信號的PMU相量數(shù)據(jù)進行FFT分析結(jié)果的Matlab仿真圖，其原始信號為

由圖3可知，PMU相量中含有35.55Hz的頻率分量，因此次同步振蕩的頻率為14.45Hz。

圖3 次同步振蕩參數(shù)辨識

3.2 幅值調(diào)制的影響

IEEE Std C37.118—2011標準對PMU的動態(tài)性能提出了明確的規(guī)范要求[11]，幅值調(diào)制是其中的一種，其定義為

式中，a為幅值調(diào)制深度；a為調(diào)制頻率。本文選擇如下信號對其PMU測量相量進行FFT頻譜分析，其幅值調(diào)制深度為0.1，調(diào)制頻率為5Hz。

其FFT頻譜分析結(jié)果如圖4所示。

圖4 幅值調(diào)制的FFT頻譜分析

根據(jù)FFT的頻譜分析結(jié)果，在相量數(shù)據(jù)中存在5.469Hz的頻譜分量，根據(jù)基于相量數(shù)據(jù)的次同步振蕩辨識流程，PMU裝置將認為系統(tǒng)中存在44.531Hz的次同步振蕩分量，但實際系統(tǒng)中則是5Hz的幅值調(diào)制信號。因此在幅值調(diào)制的條件下，根據(jù)PMU相量數(shù)據(jù)進行FFT頻譜分析從而辨識次同步振蕩分量參數(shù)會出現(xiàn)異常。由此可見，在幅值調(diào)制的條件下，基于相量數(shù)據(jù)的次同步振蕩參數(shù)辨識將會引起裝置的誤判，從而實現(xiàn)錯誤的次同步振蕩振蕩告警。

4 基于原始采樣值的次同步振蕩參數(shù)辨識

4.1 辨識流程

如前所述，對于PMU裝置而言，產(chǎn)生次同步振蕩時均表現(xiàn)為所測量的電流信號中疊加了次同步振蕩分量，因此可對電流輸入信號做FFT分析，從而計算出次同步振蕩分量的頻率和幅值。因此本文修改了文獻[6]提出的PMU相量計算模型，如圖5所示。

圖5 基于原始采樣值的次同步振蕩參數(shù)辨識

4.2 FFT頻譜分析

基于原始采樣值的次同步振蕩參數(shù)辨識方法其優(yōu)點是在進行FFT頻譜分析后，可直接根據(jù)其頻譜的頻率分量判斷出次同步振蕩的頻率分量，這是與基于PMU相量數(shù)據(jù)進行次同步振蕩參數(shù)辨識方法不同之處。也正是因為如此，幅值調(diào)制對基于原始采樣值的次同步振蕩參數(shù)辨識沒有任何影響。因此避免了對次同步振蕩分量的誤判，提高了次同步振蕩分析監(jiān)測的準確性和適應(yīng)性。本文選取的采樣率為200Hz，采樣點數(shù)512點，則根據(jù)式（3），可得其頻率分辨率為0.390625Hz。再根據(jù)香農(nóng)采樣定理該方法可分析輸入頻率最高為100Hz的輸入信號，因此，基于原始采樣值的FFT分析方法還可用于分析50～100Hz之間的超同步振蕩現(xiàn)象，若進一步提高采樣頻率，則可分析更高頻率的超同步振蕩分量。

本文對式（6）所示的信號直接進行了FFT頻譜分析，其結(jié)果圖6所示?？芍?，基于原始采樣值的次同步振蕩參數(shù)辨識方法正確的識別出了振蕩分量。

圖6 FFT頻譜分析結(jié)果

5 結(jié)論

隨著電網(wǎng)規(guī)模地不斷擴大，電力系統(tǒng)中安裝的串聯(lián)補償裝置、高壓直流（HVDC）控制系統(tǒng)、動態(tài)無功補償裝置等設(shè)備有可能會引起電力系統(tǒng)次同步振蕩等現(xiàn)象。而隨著PMU布點數(shù)目的不斷增多，使得WAMS主站端的軟硬件開銷越來越多，因此提出了在PMU裝置上實現(xiàn)分布式的次同步振蕩參數(shù)在線辨識方法，指出了與基于相量數(shù)據(jù)的次同步振蕩參數(shù)辨識方法相比，基于原始采樣值的次同步振蕩參數(shù)辨識方法更適合進行次同步振蕩參數(shù)辨識，而且其擴展性更好。但對數(shù)據(jù)直接進行FFT分析，其結(jié)果精度尚不是很高，還需做進一步研究，以提高FFT算法的頻譜分析精度。

[1] 徐政, 羅惠群. 電力系統(tǒng)次同步振蕩問題的分析方法概述[J]. 電網(wǎng)技術(shù), 1999, 23(6): 36-39.

[2] 張學群, 康積濤, 李林, 等. 基于SVC利用遠端信號緩解次同步諧振研究[J]. 電氣技術(shù), 2010, 11(12): 21-23.

[3] 蔡婷婷, 徐青山, 劉瑜俊, 等. TCR型SVC抑制電力系統(tǒng)次同步振蕩研究[J]. 電測與儀表, 2015(15): 26-33.

[4] 李偉, 肖湘寧, 趙洋. 無功發(fā)生源抑制次同步振蕩的機理分析[J]. 電工技術(shù)學報, 2011, 26(4): 168-174.

[5] 許樹楷, 謝小榮, 辛耀中. 基于同步相量測量技術(shù)的廣域測量系統(tǒng)應(yīng)用現(xiàn)狀及發(fā)展前景[J]. 電網(wǎng)技術(shù), 2005, 29(2): 44-49.

[6] 張敏, 沈鍵, 侯明國, 等. 相量測量單元實現(xiàn)次同步振蕩在線辨識和告警的探討[J]. 電力系統(tǒng)自動化, 2016, 40(16): 143-152.

[7] 侯明國, 沈健, 檀永, 等. 相量測量裝置實時數(shù)據(jù)的三段式變頻傳輸方法[J]. 電力系統(tǒng)自動化, 2015(21): 93-98.

[8] IEEE Std C37.118 Standard for Synchrophasor Measurements for Power Systems[S]. 2011.

[9] 王亮, 謝小榮, 姜齊榮, 等. 大規(guī)模雙饋風電場次同步諧振的分析與抑制[J]. 電力系統(tǒng)自動化, 2014, 38(22): 26-31.

[10] 沈健, 周斌, 汪昀, 等. 提升PMU動態(tài)測量性能的若干方法[J]. 電力系統(tǒng)保護與控制, 2015, 43(19): 89-94.

[11] 畢天姝, 劉灝, 楊奇遜. PMU算法動態(tài)性能及其測試系統(tǒng)[J]. 電力系統(tǒng)自動化, 2014, 38(1): 62-67.

On Line Identification Method of Sub Synchronous Oscillation based on Distributed PMU

Chen Jian1Zhou Yuzhi2

(1. State Grid Jiangsu Technican Training Center, Suzhou, Jiangsu 215004; 2. State Grid Suzhou Power Supply Company, Suzhou, Jiangsu 215004)

In order to reduce the overhead of software and hardware of WAMS main station, and to avoid the frequency aliasing, this article imposed an method of online parameter identification on sub synchronous oscillation, based on the analysis of the communication architecture of WAMS/PMU system and the principle of synchronized phasor measurement algorithm. Two different parameter identification methods are analyzed and compared in this paper, including parameter identification method based on phasor data and parameter identification method based on the original sample value. The limitations of the parameter identification method based on phasor data are discussed in this paper. It pointed out that the parameter identification method based on the original sample value avoids miscarriage of justice of the sub synchronous oscillation, and improves the accuracy and adaptability of the analysis and monitoring of the sub synchronous oscillation.

sub synchronous oscillation；distributed measurement；FFT；PMU

陳 劍（1972-），女，本科，工程師，主要研究方向：智能變電站二次運維及變電站自動化。