基于PSCAD的閃變源定位數字仿真

吳瑾櫻陳錦植蘇文博

（1. 福州大學電氣工程與自動化學院，福州　350116；2. 國網福建省電力有限公司寧德供電公司，福建 寧德　352000）

?

基于PSCAD的閃變源定位數字仿真

吳瑾櫻1陳錦植2蘇文博1

（1. 福州大學電氣工程與自動化學院，福州350116；2. 國網福建省電力有限公司寧德供電公司，福建 寧德352000）

摘要閃變污染會降低電網的電能質量，因此定位電力系統中的閃變源，劃分電力系統中系統側與用戶側的閃變水平十分必要。本文介紹間諧波有功功率流向、危害評估的基本原理，以及其它識別閃變的工程方法，以PSCAD為平臺搭建仿真模型，結合Matlab編程實現電網中的閃變源定位。對各種定位方法的結果及各自的特點進行對比分析，總結出這些方法在不同工況下各自的優勢和劣勢。

關鍵詞：電能質量；閃變；間諧波；閃變功率；定位

電能是一種能源形態的特殊商品，在安全經濟的基礎上要保證質量[1]。配網中的非線性負荷是導致電能質量污染的主要因素，電能質量惡化不僅影響系統運行，也會影響用戶生產，所以，一旦閃變發生，首先要搜索污染負荷，確定它在網絡中的位置。

閃變是人眼對照度的視覺反映，長時間閃變影響人的正常工作和生活，而且會進一步引起電網更大范圍的電壓質量問題[2]。閃變儀檢測到的是公共連接點的電壓閃變強度，這是背景干擾與多個負荷共同作用的結果，單純檢測電壓閃變不能直接判斷閃變干擾的來源[3]。因此通過公共連接點的測量，尋找出準確、實用的對配電網閃變源定位的方法，并確定公共連接點閃變污染責任的指標具有重要意義。

本文重點分析了有功功率流向法和注入水平評估法的基本原理，并介紹了其他閃變源識別的方法。利用PSCAD仿真軟件建立波動負荷仿真模型，通過Matlab計算波動功率和注入水平指標，比較各結果的符號或大小進行全網閃變源定位。最后以IEEE-33節點系統模型中接入波動負荷為例，結合兩種方法進行定位，驗證了各方法的有效性和存在的局限。

1　基本原理

1.1有功功率識別閃變源

閃變的數學模型通常是用幅值調制的形式給出

式中，mu(t)、mi(t)分別為電壓、電流調幅波信號；A、ωc為工頻幅值和角頻率。

若電網中出現關于50Hz對稱的成對的間諧波，可以通過計算間諧波功率的符號在整個電網中搜尋干擾源的位置。

間諧波有功功率定義[4]：

式中，θuh為間諧波電壓相量的相角，θih為間諧波電流相量的相角。

同時，文獻[5]中給出了閃變功率的定義式：

間諧波有功或者閃變功率流向的參考方向如圖1所示，功率流向的判據見表1。

圖1　基波功率參考方向

表1　試樣說明

1.2電壓電流比率法

普通負荷在運行過程中電壓和電流均保持正弦，閃變干擾使得波形畸變，通過變化的電壓與電流之間的尋找數學關系，也可以定性分析是否存在閃變，判斷來源。

圖2　等效電路

由圖2電路知識可知，Z1不變時，PCC處的電壓隨著Z2的增大而增大，流過PCC的電流也減??；Z2不變時，PCC處的電壓隨著Z1的增大而增大，而此時流過PCC 的電流減小。

在PCC點處計算電壓電流有效值序列Upcc(k）、Ipcc(k)，并進一步按式（5）進行計算。

式中，k=0,1,2,…,n。若式（5）計算結果大于0，則系統側為主要干擾源；反之，則用戶側為干擾源。

1.3電壓電流相關系數法

對于電壓電流信號，可通過用戶電壓電流測量數據計算關聯系數判斷是否發生閃變——若發生閃變，則其相關系數很低[6]，相關系數計算式為

式中，uk、ik分別為信號u和i的采樣值（k=1,2,…,n）。和i分別為uk、ik的數學期望值。

1.4閃變干擾用戶危害評估

干擾用戶的閃變危害評估以公共連接點的電壓、電流測量參數為觀測數據，評估單個用戶可能對電網電能質量造成的影響[7]。

圖3　評估場景

利用動態相量法將時域信號轉化為相量的形式，可以避免微分方程的繁瑣計算。對式（7）先后以工頻、閃變頻率為基頻計算一階動態相量[8]，得到：

IEC 61000-4-15[9]標準提供了對應瞬時閃變視感度S(t)=1的正弦電壓波動值。因此從式（8）可計算各個用戶單獨接入系統造成的危害值，逐個比較危害值的大小，可以對配網的干擾源進行定位。

2　閃變源定位的數字仿真

PSCAD是電力系統中廣泛應用的電磁暫態仿真軟件，能夠求解大型電力系統暫態過程。本文利用PSCAD搭建含閃變干擾負荷的簡單電力系統，并監測PCC點的電壓、電流數據，實現閃變源的定位。對于較為復雜的數值計算，如閃變注入水平評估所用的動態相量法，則采用Matlab實現。

2.1檢測元件建模

在PSCAD中建立基于間諧波功率的閃變源識別方法的檢測模型，如圖4所示。該檢測模型分為兩部分實現：第一部分實現間諧波檢測，如圖4（a）所示，通過FFT元件實現，分析檢測結果，選擇含量最大的間諧波頻率作為計算間諧波功率的頻率；第二部分實現間諧波功率計算，如圖4（b）所示，從第一部分的檢測結果中取得間諧波電壓、電流的幅值和相角計算間諧波功率。該檢測模型需要對每一條支路計算其間諧波功率，比較其大小，大的為閃變源支路。

圖4　基于間諧波功率的閃變源識別方法

在PSCAD中建立電壓電流比率法的檢測模型，如圖5所示。該模型首先計算電壓、電流的有效值，此后利用延時元件計算電壓（或電流）信號在一定時間間隔內的變化量，結果為一序列數值點，然后將電壓的變化量除以電流的變化量，結果即為電壓關于電流的變化率。若變化率大于0，則電壓隨著電流的增大而增大，反之，則電壓隨著電流的增大而減小。對每一條支路的電壓、電流監測數據都進行計算，對照表1分析計算結果，即可實現閃變源的定位。

圖5　電壓電流比率法檢測模型

由于電壓電流相關系數法和閃變干擾用戶危害評估的方法需要進行較為復雜的數值計算，因為本文利用Matlab軟件實現其計算過程。

2.2簡單工況

在PSCAD軟件中搭建10kV簡單供電系統，如圖6所示。系統側等效為受控電壓源和系統阻抗的并聯，通過改變頻率和幅值可以模擬系統側產生的背景干擾。用戶側建模時，恒定負荷可以等效為串聯阻抗的形式，而對于具有波動性的負荷，可以用時變電阻和時變電感來模擬。系統參數如下。

系統短路容量：200MVA；

電壓源：Usrms=10×(1+0.002cos(2πfst))；

1號負荷功率：S1=3×(1+0.1cos(2πf1t))；

2號負荷功率：S2=2×(1+0.1cos(2πf2t))；

功率因數：0.85。

其中，系統電壓和1、2號負荷均含有波動分量，其波動頻率分別為fs、f1、f2。系統監測量為PCC點的電壓和各條支路的電流，通過電壓、電流監測量可以計算波動功率和比率，也可以檢測負荷模型的r(t)和l(t)，并結合系統額定電壓、短路容量計算用戶的干擾評估指標。

圖6　仿真模型

仿真時，通過改變fs、f1、f2的值，可以模擬閃變干擾源在不同位置上的運行工況，包括是否含有背景干擾以及改變閃變源的位置等。本文設置了以下4種運行工況：

算例1：干擾全部來自系統側，fs=5Hz，f1=0Hz，f2=0Hz。算例2：干擾來自負荷1，fs=f2=0Hz，f1=10Hz。算例3：干擾同時來自負荷1和負荷2，fs=0Hz，f1=f2=10Hz。

算例4：干擾來自系統和負荷2，fs=10Hz，f1= 0Hz，f2=10Hz。

2.3IEEE-33節點算例

在PSCAD中對圖7所示的IEEE-33節點配電網標準測試系統進行仿真。

圖7　IEEE-33節點配電網

設定l1、l2為節點7、19接入的兩個波動負荷，如圖7所示。其中，線路l1流過的功率大小為：Sl1=10×(1+0.1cos(2πfl1t))；線路l2流過的功率大?。篠l2=10×(1+0.1cos(2πfl2t))，根據監測數據進行閃變源定位。

3　仿真結果分析

3.1簡單工況的定位仿真分析

算例1：fs=5Hz，f1=0Hz，f2=0Hz，此時干擾全部來自系統側。計算結果如表2、圖8（a）、圖8（b）。

算例2：fs=f2=0Hz，f1=10Hz計算結果如表3和圖9（a）、圖9（b）。

算例3：fs=0Hz，f1=f2=10Hz，計算結果如表4和圖10（a）、圖10（b）所示。

表2　計算結果

圖8　比率法計算結果

表3　計算結果

圖9　比率法計算結果

表4　計算結果

圖10　比率法計算結果

算例4：fs=10Hz，f1=0Hz，f2=10Hz，計算結果如表5和圖11（a）、圖11（b）所示。

表5　計算結果

以上分別列出了采用4種方法進行閃變源定位的分析計算結果。以算例1為例說明：從表2可以看出，成對的間諧波功率的符號全為正，由此可以得出系統側存在閃變源。PSCAD計算得到的電壓電流比率計算結果如圖8（a）所示，可以看出大部分的比率值大于0，在部分點上存在比率值小于0情況，“電壓有效值-電流有效值”曲線如圖8（b）所示，觀察該曲線發現在部分點上存在電壓隨著電流增大而減小的情況，與PSCAD所得結果一致。但兩幅圖上，電壓和電流有效值的變化趨勢是電壓隨著電流的增大而增大，據此可以得出系統側存在閃變干擾的結論。表2中，危害指標系統側遠大于負荷側，即系統側存在閃變源。

圖11　比率法計算結果

3種方法顯示的干擾源和實際設置一致，從而驗證了各種方法用來定位的有效性。電壓電流的相關性可以作為補充判據。算例4中可以發現，若系統側和負荷側同時存在同樣頻率的干擾，功率流反映出二者“中和”后的結果。危害指標反映出各個用戶和系統的危害值，哪條出線側的危害值大，就沿著該出線繼續評估。但是從表中也能看見，由于電流之間的關聯性，危害評估值會把恒定負荷判定為輕微干擾。電壓電流的比率結果在單頻率干擾下可以判斷是來于系統側還是負荷側。對于多干擾源存在，只能得到干擾相互綜合后的主導干擾源位置。

3.2IEEE-33節點系統定位分析

當波動頻率fl1=10Hz，fl2=20Hz時，仿真結果見表6。

表6　計算結果

在IEEE-33節點模型中，從表6可以看出，20Hz調幅波所對應頻率成分的間諧波功率在線路2-19及l1均為負，因此可以判定l1為20Hz調幅波的閃變干擾源支路，同理，l2為10Hz調幅波的閃變干擾源支路。閃變注入水平的計算結果發現，2號節點的兩條出線側存在閃變干擾，比較注入水平沿出線繼續尋找干擾源，在19號節點上發現l1閃變干擾最為嚴重，但是線路19-20也被判定存在輕微干擾，從定位結果來看，兩種方法在結果上與PSCAD仿真實驗所提供的情況相吻合，由此也驗證了以上閃變干擾源定位方法的可行性。

4　結論

本文介紹了目前閃變源定位和識別的工程方法。根據間諧波功率與基波功率的流向關系來判斷閃變源的位置，通過注入水平評估來比較各個出線側的閃變情況，從而實現定位。通過數字仿真在PSCAD中建立含有閃變干擾情況的配電網，基于Matlab軟件編寫閃變源定位的分析程序，通過仿真實驗來驗證以上所提方法的正確性和存在的問題。

參考文獻

[1] 肖湘寧, 徐永海. 電能質量問題剖析[J]. 電網技術, 2001, 25(3): 66-69.

[2] 肖湘寧. 電能質量分析與控制[M]. 北京: 中國電力出版社, 2010.

[3] 邵振國, 吳丹岳, 薛禹勝. 閃變干擾源定位研究中的若干問題[J]. 高電壓技術, 2009, 35(7): 1595-1599.

[4] 王建勛, 劉會金. 間諧波閃變效應計算及閃變源識別方法[J]. 電力系統自動化, 2011, 35(12): 52-58, 75.

[5] Axelberg P, Bollen M, Gu I. Trace of flicker sources by using the quantity of flicker power[J]. IEEE Transactions on Power Delivery, 2008, 23(1): 465-471.

[6] 齊勇, 江亞群, 黃純, 等. 諧波源和閃變源的識別算法[J]. 電力自動化設備, 2011, 31(11): 60-63, 68.

[7] 邵振國, 吳丹岳, 林焱. 波動負荷的閃變危害評估[J]. 電力系統保護與控制, 2009, 37(14): 36-40.

[8] 蘇清梅, 張嫣, 邵振國. 基于波動模型及虛擬場景評價的閃變污染用戶危害評估[J]. 電力系統保護與控制, 2013, 41(23): 97-102.

[9] IEC 61000-4-15. Testing and measurement techniques: flickermeter—functional and design specifications[S].

吳瑾櫻（1991-），女，福建省三明市人，碩士研究生，研究方向為電能質量。

Digital Simulation of Flicker Source Location Using PSCAD

Wu Jinying1Chen Jinzhi2Su Wenbo1
(1. College of Electrical Engineering and Automation, Fuzhou University, Fuzhou350116；2. Electric Power Bureau of Ningde, Ningde, Fujian352000)

Abstract Flicker disturbance will reduce the power quality, therefore it is quite necessary to locate the flicker source and divide the flicker responsibility from customer and supplier. Firstly, this paper introduces the basic theory about the power-direction method, the injection level assessment method and other methods for identifying the flicker source. Then, simulation of flicker source localization with PSCAD was carried on combining with Matlab programming. Finally, the characteristics of different methods are analyzed with comparison of their corresponding simulation results and their advantages and disadvantages under different conditions are pointed out.

Keywords：power quality; flicker; interharmonics; flicker power; localization

作者簡介