楊洋，肖湘寧，陶順，陳鵬偉

(新能源電力系統(tǒng)國家重點實驗室(華北電力大學)，北京市 102206)

考慮衰減直流分量的dq-120改進算法及其在混合仿真中的應(yīng)用

楊洋，肖湘寧，陶順，陳鵬偉

(新能源電力系統(tǒng)國家重點實驗室(華北電力大學)，北京市 102206)

該文分析了故障電流含有衰減直流分量時傳統(tǒng)dq-120算法的誤差形成機理及表現(xiàn)形式。在此基礎(chǔ)上，結(jié)合Prony算法提出了一種dq-120改進算法，該算法基于故障電流中含有直流分量的假設(shè)條件，對故障電流的dq-120計算結(jié)果進行模型降階，簡化了傳統(tǒng)Prony算法的計算量并提升了dq-120算法對直流分量的免疫水平。借助PSCAD/EMTDC仿真軟件，搭建了含有一回直流線路的IEEE 39節(jié)點系統(tǒng)機電-電磁暫態(tài)混合仿真模型。仿真結(jié)果表明，與傳統(tǒng)的dq-120算法相比，改進算法明顯提高了故障期間相量提取精度，從而有效提升了機電-電磁暫態(tài)混合仿真故障期間的仿真精度以及整個仿真過程的仿真精度。

衰減直流分量；dq-120；Prony算法；混合仿真；PSCAD

0 引 言

基于單相變換平均值法的dq-120算法在電力電子控制中得到了廣泛應(yīng)用[1-3]。經(jīng)過坐標變換，可以將正弦量變換到旋轉(zhuǎn)坐標系，從而得到直流形式表示的電壓或電流相量的實部和虛部。由于計算量小，該方法也應(yīng)用在機電-電磁暫態(tài)混合實時仿真領(lǐng)域，并取得一定效果[4-5]。當接口處發(fā)生接地故障時，故障電流中往往伴隨直流分量，經(jīng)坐標變換后疊加到原來的直流量中，影響故障期間相量的提取精度，進而影響機電側(cè)的求解。為了濾除故障電流中的衰減直流分量，提高故障期間電流相量的提取精度，專家和學者進行了一些研究，主要有兩種思路。一種思路是在對電流進行處理前，通過數(shù)字濾波技術(shù)對獲得的電流波形進行處理[6]，濾除電流中所含有的非周期分量成分。另一種思路是針對計算后得到的結(jié)果進行處理，通過另外增加1至2個采樣點估計出衰減直流分量的參數(shù)從而得到精確估計結(jié)果[7]。前一種方法往往需要設(shè)計濾波器參數(shù)，濾波的效果取決于所假定的模型，因此適應(yīng)性不足。后一種方法需要對得到的結(jié)果進行適當?shù)耐蒲?，公式復雜程度與采用的提取基波分量算法相關(guān)。

近年來Prony算法在相量提取方面得到一些應(yīng)用[8-11]，通過設(shè)定合適的模型階數(shù)以及數(shù)據(jù)窗長，可以較為精確地保證相量提取精度。但其主要弊端是窗長選取困難，計算量大，因而難以適用于實時計算場合。

本文將分析含有衰減直流分量情況下dq-120算法的誤差產(chǎn)生機理及表現(xiàn)形式，借鑒Prony算法的思路，通過適當簡化，提出一種考慮衰減直流分量的dq-120改進算法。該算法僅需要故障發(fā)生一個工頻周期后5個采樣點即可準確提取出故障電流的基波幅值和相位，在減少計算量的同時保證了電流含有衰減直流分量時的相量提取精度。借助PSCAD/EMTDC仿真軟件搭建含有一回直流線路的IEEE 39節(jié)點系統(tǒng)機電-電磁暫態(tài)混合仿真模型，對本文提出的改進的dq-120算法進行驗證。

1 dq-120算法誤差產(chǎn)生機理及表現(xiàn)形式

首先，以電壓為例論述dq-120算法的基本原理。采用單相dq變換法先將每一相的電氣量的幅值和相角求出。設(shè)單相電壓信號為

(1)

其中X=Umcosφ，Y=Umsinφ，ω是基波角頻率。

令

uq(t)=2u(t)×cosωt

(2)

ud(t)=2u(t)×sinωt

(3)

將上述兩式分別化簡，得：

uq(t)=2(Xsinωtcosωt+Ycos2ωt)=Xsin2ωt+Y(1+cos2ωt)=Y+Ycos2ωt+Xsin2ωt

(4)

ud(t)=2(Xsin2ωt+Ycosωtsinωt)=X(1-cos2ωt)+Ysin2ωt=X-Xcos2ωt+Ysin2ωt

(5)

對兩個信號取基波半個周期(或半個周期的整數(shù)倍)的平均值，則可以由此求出X和Y的值，從而可以得到電壓幅值為

(6)

相位為

(7)

又由正負零序電壓和三相電壓的關(guān)系

(8)

可求得基波正序電壓量為

(9)

(10)

Ur=Ucosφ，是實軸分量，Ui=Ucosφ是虛軸分量，U是由dq變換求得的相電壓幅值。

(11)

(12)

其中

(13)

(14)

式中：U1r是基波正序相量的實軸分量；U1i是基波正序相量的虛軸分量。

設(shè)電氣量中含有衰減直流分量，分析其對dq-120方法提取的影響。假定電壓信號為

u(t)=Um1sin(ωt+φ1)+Udce-τt

(15)

按照式(4)和(5)對上式分別進行處理可得

uq(t)=Y+Ycos2ωt+Xsin2ωt+2Udce-τtcosωt

(16)

ud(t)=X-Xcos2ωt+Ysin2ωt+2Udce-τtsinωt

(17)

在[t,t+T]范圍內(nèi)進行積分可得

Y′=∫tt+Tuq(t)dt=

∫tt+T(Y+Ycos2ωt+Xsin2ωt+2Udce-τtcosωt)dt=

Y-U0e-τtsin(ωt+φ1)

(18)

X′=∫tt+Tud(t)dt=

∫tt+T(X-Xcos2ωt+Ysin2ωt+2Udce-τtsinωt)dt=

X+U0e-τtcos(ωt+φ1)

(19)

2 改進的dq-120算法

由式(18)和式(19)可知，故障發(fā)生后，從得到的Y′和X′中準確地提取出Y和X是計算故障后電流相量的關(guān)鍵所在。應(yīng)用歐拉公式，式(18)和式(19)可以表示為

(20)

式中X、Y、C1、C2、C3,、C4均為常數(shù)，均為e的含t的指數(shù)冪。以X為例論述其直流分量的求取過程。

設(shè)N個采樣點為X′(0),X′(1),X′(2)，…，X′(N-1)則有

(21)

令z1、z2、z3為以下代數(shù)方程的根

z3+a1z2+a2z+a3=0

(22)

為了求出式(22)中的系數(shù)，令式(21)的第1個方程乘以系數(shù)a3，第2個方程乘以a2，第3個方程乘以a1可得：

(23)

對式(23)求和并利用式(22)可得

X′(0)a3+X′(1)a2+X′(2)a1+X′(3)=0

(24)

同理，令式(21)的第2個方程乘以系數(shù)a3，第3個方程乘以a2，第4個方程乘以a1可得：

(25)

同樣可得

X′(1)a3+X′(2)a2+X′(3)a1+X′(4)=0

(26)

由于X′含有直流分量，因此也是式(22)的根

a1+a2+a3+1=0

(27)

聯(lián)立式(24)，(26)，(27)并寫成矩陣形式，有

(28)

其中，X′(0)—X′(4)為5個采樣值，原則上最終計算結(jié)果與采樣周期無關(guān)。式(28)有3個方程，3個未知數(shù)，因而可解。求解可得：

(29)

對于Q3(z)=z3+a1z2+a2z+a3=0，可以看成某個實矩陣特征多項式，即

Q3(λ)=λ3+a1λ2+a2λ+a3=0

(30)

因此，求方程的全部根的問題就變成了求矩陣全部特征值的問題，可以驗證，上述特征多項式所對應(yīng)的矩陣為

(31)

矩陣B是上H矩陣，可以直接用QR方法求出全部特征值，得到所有的取值，依次代入式(21)的前3個方程，即可以得到包括X在內(nèi)的所有系數(shù)，而X即為對應(yīng)的直流分量。同理，可以求出Y。

3 適應(yīng)改進dq-120算法的交互時序

普遍認為機電-電磁暫態(tài)混合仿真的誤差主要來源為模型誤差，交互誤差和數(shù)據(jù)傳輸誤差，對相量提取方法引入的誤差認識不足[12]。當接口處發(fā)生故障后，傳統(tǒng)的相量提取算法不能解決直流分量問題，因此，文獻[13]在發(fā)生故障后通過增大交互周期，借助直流分量的自然衰減，得到穩(wěn)態(tài)的故障電流相量，來減小其對相量提取算法的影響；文獻[14]提出了并行轉(zhuǎn)串行的時序，故障發(fā)生后，機電側(cè)停止運行，電磁側(cè)運行若干交互步長后得到穩(wěn)定相量值，再傳遞給機電側(cè)進行計算。無論哪種方式，都犧牲了故障后的分析效率。借助本文所提的相量提取算法，可以一定程度解決這個問題。適應(yīng)于本文相量提取算法的交互時序如圖1所示。

圖1 改進的交互時序Fig.1 Improved interface protocol

假設(shè)T1時刻在電磁側(cè)的接口處發(fā)生三相接地故障，此時，電磁暫態(tài)過程或者機電暫態(tài)過程均沒有開始進行故障后的計算，機電暫態(tài)計算過程照常將接口信息送入電磁暫態(tài)計算過程(如圖1步驟(1)),電磁暫態(tài)計算過程獲取信息后，進行半個周波及Δtdelay的計算(如圖1中步驟(2))，其中Δtdelay=4Δt，再將計算得到的電流相量值傳遞給機電側(cè)(如圖1中步驟(3))，由于Δtdelay=0.2 ms<<10 ms，也就是因為相量計算引起的時間延時在機電側(cè)看來可以忽略不計，因此機電側(cè)接收到電磁側(cè)傳遞的等值信息之后可以正常進行機電計算，而無需進行修正。借助改進的交互時序，可以在保證電磁側(cè)相量提取精度的前提下盡量保證混合仿真的分析效率。

4 仿真驗證

在PSCAD/EMTDC中建立如圖2所示基于PSCAD+C架構(gòu)的機電-電磁暫態(tài)仿真模型[15]，將新英格蘭39節(jié)點系統(tǒng)中的第38節(jié)點的發(fā)電機和變壓器替換為直流輸電線路，直流線路的逆變側(cè)接入29節(jié)點，并選擇29母線為接口母線。直流模型參數(shù)及控制方式參見Cigre標準直流模型[16]。發(fā)電機采用六階模型，計及勵磁和調(diào)速。負荷采用恒阻抗模型。直流系統(tǒng)在電磁暫態(tài)側(cè)仿真，交流系統(tǒng)在機電側(cè)系統(tǒng)仿真。在每個交互周期，機電側(cè)向電磁側(cè)提供戴維南等值電勢，電磁側(cè)向機電側(cè)提供正序基波電流相量。機電側(cè)仿真步長設(shè)定為10 ms，電磁側(cè)仿真步長設(shè)定為50 μs。交互周期設(shè)定為10 ms。

圖2 修改后的IEEE 39節(jié)點系統(tǒng)Fig.2 Improved IEEE39-node system

分別采用傳統(tǒng)dq-120方法和改進dq-120算法對該算例進行仿真。在直流線路的逆變側(cè)設(shè)置經(jīng) 0.01 Ω電阻三相接地短路故障，故障持續(xù)時間為 0.1 s。選取34號發(fā)電機的轉(zhuǎn)速和機端電壓為機電側(cè)主要觀察對象，選取接口有功功率為電磁側(cè)的主要觀察對象。仿真結(jié)果對比如圖3所示。

由圖3(a)圖可知，由于接口電流在故障期間含有衰減直流分量，采用dq-120算法，得到的電流相量幅值中存在衰減的正弦分量。采用改進的dq-120算法后，故障期間電流相量中的衰減正弦分量得到抑制，故障期間得到相對穩(wěn)定的電流相量。這種改進反映在故障期間機電側(cè)求解中，從圖3(b)可以看出，采用改進算法得到的發(fā)電機轉(zhuǎn)速擺動的最大幅值更接近全電磁；從圖3(c)可以看出，采用改進算法前，發(fā)電機機端電壓在故障期間存在波動，而采用改進算法后，故障期間電壓的波動得到抑制。由于故障期間機電側(cè)的求解精度得到提升，從圖3(d)圖可以看出，故障清除之后，電磁側(cè)接口有功功率的精度也有所提升。因此，混合仿真整體的仿真精度均得到提升。

5 結(jié) 論

本文分析了含有衰減直流分量情況下dq-120算法的誤差產(chǎn)生機理及表現(xiàn)形式，結(jié)合Prony算法提出了一種dq-120改進算法并應(yīng)用到機電-電磁暫態(tài)仿真中，得到以下結(jié)論。

(1)故障期間電流中含有的衰減直流分量影響dq-120算法的提取精度。

(2)經(jīng)過另外5個采樣點，改進的dq-120算法可以準確提取出故障發(fā)生后故障電流的相量。

圖3 傳統(tǒng)方法dq-120方法與本文方法的 仿真結(jié)果對比Fig.3 Simulation results comparison between improved and traditional dq-120 algorithms

仿真結(jié)果表明，經(jīng)過對故障期間相量提取算法的改良，一方面可以直接提高故障期間機電側(cè)的仿真精度，另一方面，也會顯著提升故障后電磁側(cè)的仿真精度，從而可以進一步促進機電-電磁暫態(tài)仿真整體精度的提升。

[1]陸晶晶, 劉正富, 張劍,等. MMC型統(tǒng)一電能質(zhì)量控制器故障分析及保護策略[J]. 電力建設(shè), 2015, 36(5):31-36. LU Jingjing, LIU Zhengfu, ZHANG Jian, et al. Failure analysis and protection strategy for MMC-UPQC[J].Electric Power Construction, 2015, 36(5):31-36.

[2]胡靜, 趙成勇, 翟曉萌. 適用于MMC多端高壓直流系統(tǒng)的精確電壓裕度控制[J]. 電力建設(shè), 2013, 34(4):1-7. HU Jing, ZHAO Chengyong, ZHAI Xiaomeng. Precise voltage margin control for modular multilevel Converter based Multi-terminal HVDC System[J].Electric Power Construction,2013, 34(4):1-7.

[3]李燕平，李開成，王海元，等. 一種基于改進dq變換的電壓驟降檢測方法[J]. 華北電力大學學報(自然科學版)， 2013，40(3): 25-29. LI Yanping, LI Kaicheng, WANG haiyuan, el al. Voltage sag detection method based on improveddqtransformation[J].Journal of North China Electric Power University，2013，40(3): 25-29.

[4]賈旭東，李庚銀，趙成勇，等. 基于RTDS/CBuilder的電磁-機電暫態(tài)混合實時仿真方法[J]. 電網(wǎng)技術(shù)，2009，33(11): 33-38. JIA Xudong, LI Gengyin, ZHAO Chengyong, et al. Electromagnetic transient and electromechanical transient hybrid real-time simulation method based on RTDS/CBuilder[J]. Power system technology， 2009, 33(11): 33-38.

[5]李秋碩，張劍，肖湘寧，等. 基于RTDS的機電電磁暫態(tài)混合實時仿真及其在FACTS中的應(yīng)用[J]. 電工技術(shù)學報， 2012，27(3): 219-226. LI Qiushuo, ZHANG Jian, XIAO Xiangning, el al. Electromechanical-electromagnetic transient real-time simulation based on RTDS and its application to FACTS[J].Transactions of China Electrotechnical Society，2012，27(3): 219-226.

[6] GU J C, YU S L. Removal of DC offset in current and voltage signals using a novel Fourier filter algorithm[J]. IEEE Transactions on Power Delivery， 2000, 15(1): 73-79.

[7]馬磊，王增平，徐巖. 微機繼電保護中濾除衰減直流分量的算法研究[J]. 繼電器， 2005，33(17): 11-13. MA Lei, WANG Zengping, XU Yan. An algorithm for optimal parameter estimation of the failure rate of electrical equipment[J].Relay，2005,33(17): 11-13.

[8]丁藍，薛安成，李津，等. 基于窗口滑動改進Prony算法的電力系統(tǒng)低頻振蕩識別[J]. 電力系統(tǒng)自動化， 2010，34(22): 24-28. DING Lan, XUE Ancheng, LI Jin, el al. A moving-window Prony algorithm for power system low frequency oscillation identification[J]. Automation of Electric Power Systems,2010,34(22):24-28.

[9]竺煒, 唐穎杰, 周有慶,等. 基于改進Prony算法的電力系統(tǒng)低頻振蕩模式識別[J]. 電網(wǎng)技術(shù), 2009，33(5):44-47. ZHU Wei,TANG Yingjie, ZHOU Youqing, et al. Identification of power system low frequency oscillation mode based on improved Prony algorithm[J]. Power System Technology，2009，33(5):44-47.

[10]李大虎，曹一家. 基于模糊濾波和Prony算法的低頻振蕩模式在線辨識方法[J]. 電力系統(tǒng)自動化， 2007，31(1): 14-19. LI Dahu, CAO Yijia. An online identification method for power system low-frequency oscillation based on fuzzy filtering and Prony algorithm[J]. Power System Automation，2007，31(1): 14-19.

[11]肖晉宇，謝小榮，胡志祥，等. 電力系統(tǒng)低頻振蕩在線辨識的改進Prony算法[J]. 清華大學學報(自然科學版)， 2004, 44(7): 883-887. XIAO Jinning, XIE Xiaorong, HU Zhixiang, et al. Improved Prony method for online identification of low-frequency oscillations in power systems[J]. Journal of Tsinghua University (Science and Technology),2004,44(7):883-887.

[12]張樹卿，梁旭，童陸園，等. 電力系統(tǒng)電磁/機電暫態(tài)實時混合仿真的關(guān)鍵技術(shù)[J]. 電力系統(tǒng)自動化， 2008，32(15): 89-96. ZHANG Shuqing, LIANG Xu, TONG Luyuan, et al. Key technologies of the power system electromagnetic/electromechanical real-time hybrid simulation[J]. Automation of Electric Power Systems,2008,32(15):89-96.

[13]岳程燕. 電力系統(tǒng)電磁暫態(tài)與機電暫態(tài)混合實時仿真的研究[D]. 北京：中國電力科學研究院, 2005. YUE Chengyan. Study of power system electromagnetic transient and electromechanical transient real-time hybrid simulation[D].Beijing: China electric power research institute,2005.

[14] VAN DER MEER A A, GIBESCU M, VAN DER MEIJDEN M A M M, et al. Advanced hybrid transient stability and EMT simulation for VSC-HVDC Systems[J]. IEEE Transactions on Power Delivery， 2015, 30(3): 1057-1066.

[15]李偉，楊洋，陳鵬偉，等. 電磁-機電暫態(tài)混合仿真誤差傳遞機理分析[J]. 南方電網(wǎng)技術(shù)， 2015,9(9): 92-97. LI Wei, YANG Yang, CHEN Pengwei, et al. Analysis on error mechanism for electromagnetic-electromechanical transient hybrid simulation[J]. Southern power system technology， 2015, 9(9): 92-97.

[16] SZECHTMAN M, WESS T, THIO C V. A benchmark model for HVDC system studies[C]//AC and DC Power Transmission, 1991., International Conference on. London: IET, 1991: 374-378.

(編輯 張小飛)

An Improveddq-120 Algorithm Considering Decaying DC Component and Its Application in Hybrid Simulation

YANG Yang, XIAO Xiangning, TAO Shun, CHEN Pengwei

( State Key Laboratory for Alternate Electrical Power System with Renewable Energy Sources (North China Electric Power University), Beijing 102206, China)

This paper analyzes the error forming principle of traditionaldq-120 algorithm when the fault current containing decaying DC component and its display form. On this basis, an improveddq-120 algorithm is proposed combined with Prony algorithm. The improved algorithm can reduce the model order for thedq-120 calculation results of fault current based on the assumption that the fault current contains DC component, simplify the calculation of traditional Prony algorithm and improve the immune level ofdq-120 algorithm to the DC component. With using PSCAD/EMTDC simulation software, we build an electromechanical-electromagnetic hybrid transient simulation model of IEEE 39-node system with a DC link. The simulation result shows that, compared with traditionaldq-120 algorithm, the improved algorithm improves the phasor extraction accuracy of current during fault, so as to improve the accuracy of electromechanical-electromagnetic hybrid transient simulation during fault as well as the whole simulation process.

decaying DC component;dq-120; Prony algorithm; hybrid simulation; PSCAD

中央高校基本科研業(yè)務(wù)費專項資金資助項目(2015XS22)

TM 743

A

1000-7229(2016)06-0043-06

10.3969/j.issn.1000-7229.2016.06.007

2016-03-22

楊洋(1989)，男，通信作者，博士研究生，主要研究方向為電力系統(tǒng)仿真與分析;

肖湘寧(1953)，男，教授，博士生導師，主要研究方向為電力電子技術(shù)、現(xiàn)代電能質(zhì)量和高壓直流輸電技術(shù);

陶順(1972)，女，副教授，從事智能配電網(wǎng)和電能質(zhì)量等方面的教學與科研工作;

陳鵬偉(1992)，男，博士研究生，研究方向為電力系統(tǒng)仿真技術(shù)與微電網(wǎng)技術(shù)。

Project supported by the Fundamental Research Funds for the Central Universities(2015XS22)