復(fù)雜電網(wǎng)下基于CSDFT-MAF-PLL的電壓同步方法

陳 偉 ，何 山 ，2

（1.新疆大學(xué)電氣工程學(xué)院，烏魯木齊 830049；2.可再生能源發(fā)電與并網(wǎng)控制教育部工程研究中心，烏魯木齊 830049）

以風(fēng)電、光伏為代表的新能源發(fā)電并網(wǎng)運行，需準確獲取電網(wǎng)電壓的相角、頻率和幅值等信息，以提升系統(tǒng)故障恢復(fù)時間及并網(wǎng)能力[1-2]。鎖相環(huán)PLL（phase-locked loop）作為并網(wǎng)的關(guān)鍵環(huán)節(jié)，直接影響電力轉(zhuǎn)換系統(tǒng)的穩(wěn)定性與可靠性[3]，其快速精準的鎖相能力對電力電子裝置并網(wǎng)具有十分重要的作用。

同步旋轉(zhuǎn)坐標系鎖相環(huán)SRF-PLL（synchronous reference frame PLL）在電力電子和電力系統(tǒng)中應(yīng)用廣泛[4-5]，理想電網(wǎng)下，能快速準確地提取相位信息；復(fù)雜電網(wǎng)下，需降低鎖相環(huán)帶寬，以提高諧波抑制能力，但響應(yīng)速度下降[6]。文獻[7]提出SRFPLL存在相位與頻率檢測緊密耦合問題，諧波或任何類型的故障失衡引起的誤差嚴重影響鎖相性能。文獻[8]提出附加前置解耦單元的準一階鎖相環(huán)，結(jié)合級聯(lián)延遲消去法，增強系統(tǒng)抗諧波能力，該前置濾波環(huán)節(jié)只適合頻率波動小的系統(tǒng)，復(fù)雜電網(wǎng)條件下效果不佳。

為克服SRF-PLL在抑制直流分量及諧波方面的不足，文獻[9]利用解耦雙同步參考坐標系鎖相環(huán)DDSRF-PLL（decoupled double SRF-PLL）提取三相電壓正、負序分量，利用解耦模塊消除擾動，但該方法結(jié)構(gòu)復(fù)雜，給系統(tǒng)帶來一定延遲。

為克服DDSRF-PLL結(jié)構(gòu)復(fù)雜和延遲問題，文獻[10]提出雙二階廣義積分器鎖相環(huán)DSOGI-PLL（double second-order generalized integrator PLL），電網(wǎng)畸變程度低時，能快速有效鎖定相位，但在含直流偏置及多次諧波等復(fù)雜電網(wǎng)下，提取的同步信號會出現(xiàn)較大波動。

為克服DSOGI正交信號發(fā)生器DSOGI-QSG（DSOGI-quadrature signals generator）在復(fù)雜電網(wǎng)下提取同步信號困難問題，文獻[11]提出結(jié)構(gòu)簡便、濾波程度高的滑動平均濾波鎖相環(huán)MAF-PLL（moving average filter PLL），滑窗寬度的選取基準直接影響濾波效果及動態(tài)特性。為克服該弊端，文獻[12]提出雙滑動平均濾波鎖相環(huán)，并附加角頻率重構(gòu)模塊以消除鎖相誤差，該滑窗寬度設(shè)置較大，提升了動態(tài)特性，但該鎖相環(huán)設(shè)計只針對單相系統(tǒng)。

滑動離散傅立葉算法SDFT（sliding discrete Fourier transform）是諧波測量中使用最廣泛的算法之一。電力系統(tǒng)的諧波成分通常是時變的，特別是在較長傳輸距離的弱配電網(wǎng)中，直接進行頻譜分析會由于泄漏和柵欄效應(yīng)導(dǎo)致結(jié)果不準確。文獻[13]針對SDFT頻譜泄漏問題，提出建立泄漏狀態(tài)下的相角或幅值關(guān)系，間接計算電網(wǎng)頻率，但復(fù)雜系統(tǒng)對關(guān)系式本身會造成影響。

針對上述問題，文中設(shè)計思路如下：首先，為解決常規(guī)SDFT容易出現(xiàn)頻譜泄漏問題，提出固定滑窗-變滑窗結(jié)構(gòu)的級聯(lián)型滑動離散傅立葉變換方法，后級SDFT滑窗寬度依據(jù)實時檢測的電網(wǎng)頻率調(diào)整，利用DFT相角特性，生成每一相虛擬正交信號，提取基波信號；其次，設(shè)計附加相位補償MAF正、負序分離模塊，既能準確濾除負序分量，又能避免因MAF階數(shù)過高造成的相位延遲；最后，仿真驗證方法有效性。

1 復(fù)雜電網(wǎng)多故障誤差分析

以SRF-PLL為例，分析系統(tǒng)在出現(xiàn)三相電壓不平衡、含直流偏移及多次諧波分量下的動態(tài)性能。假設(shè)：

（1）輸入信號為理想三相正弦波；

（2）在變量發(fā)生變化之前，系統(tǒng)已達到穩(wěn)態(tài)；

（3）推導(dǎo)過程在時域下進行，系統(tǒng)視為零狀態(tài)系統(tǒng)。

式中：vd、vq為 d、q坐標下電壓分量；vd_dc、vq_dc為直流偏移分量在d、q坐標下產(chǎn)生的電壓分量；vd_U、vq_U為電壓相位不平衡在d、q坐標下產(chǎn)生的電壓分量；vd_h、vq_h為諧波分量在d、q坐標下產(chǎn)生的電壓分量；Vdc_α、Vdc_β為α、β坐標下直流偏移分量；θ為檢測的三相電壓相位角；ε為相位誤差。

由式（2）可知：由直流偏移引起的誤差具有與電網(wǎng)電壓相同的頻率分量；三相電壓不平衡誤差分量為2倍頻擾動分量；電網(wǎng)非3倍奇次諧波分量會產(chǎn)生6倍頻擾動誤差。理想電網(wǎng)下，三相電壓信號通過坐標變換轉(zhuǎn)換為dq坐標下的電壓分量，SRFPLL對信號進行處理，可快速精確地跟蹤相位，但其結(jié)構(gòu)無法克服在三相電壓不平衡、含直流偏移及多次諧波污染下的擾動誤差。

2 CSDFT-QSG頻率自適應(yīng)設(shè)計

2.1 CSDFT-QSG設(shè)計

DFT算法可提取電網(wǎng)輸入信號中的基波信號，消除鎖相外環(huán)中的諧波分量及直流分量，但需要計算完整一個周期才更新一次數(shù)據(jù)?；跇颖緮?shù)據(jù)的相似性，由前一階段窗口函數(shù)的頻譜特性所生成的遞推關(guān)系，可推算出當(dāng)前時刻的頻譜。SDFT算法依據(jù)固定窗口隨時間滑動提取樣本數(shù)據(jù)[14]。

由式（4）、（5）可知，輸入電網(wǎng)電壓信號，生成其虛擬正交信號，提取基波信號，具體結(jié)構(gòu)如圖1所示。

圖1 SDFT-QSG結(jié)構(gòu)Fig.1 Structure of SDFT-QSG

圖1中，vα(k) 、vqα(k)為前級結(jié)構(gòu)輸出的虛擬正交信號。

由式（3）～（5）得SDFT-QSG的z域傳遞函數(shù)為[16]

傳遞函數(shù)伯德圖如圖2所示。

圖2 SDFT-QSG伯德圖Fig.2 Bode diagram of SDFT-QSG

由圖2可知：SDFT-QSG在指定頻率處，幅值增益衰減大，理論上可完全濾除畸變電網(wǎng)下多次諧波及直流偏移分量。

由文獻[15]可知，若SDFT采樣頻率非電網(wǎng)電壓基波頻率的N倍，系統(tǒng)將會發(fā)生頻譜泄漏，影響采樣精度，且SDFT-QSG提取的基波信號也會產(chǎn)生偏差，諧波與直流分量無法完全消除，系統(tǒng)會存在穩(wěn)態(tài)誤差。

針對此問題，設(shè)計一種固定滑窗寬度-變滑窗寬度相結(jié)合機制的級聯(lián)型滑動窗離散傅立葉變換結(jié)構(gòu)。

由圖2可知，vα(k)在基波頻率處與v(k)同相位，故以 vα(k)作為后級設(shè)計輸入端，其結(jié)構(gòu)如圖3所示。

圖3 CSDFT-QSG結(jié)構(gòu)Fig.3 Structure of CSDFT-QSG

2.2 后級頻率自適應(yīng)設(shè)計

由式（3）可知，通過改變樣本周期點數(shù)N和系統(tǒng)的采樣周期Ts都可改變?yōu)V波器的窗口長度Tω，以達到后級SDFT-QSG頻率自適應(yīng)。改變系統(tǒng)采樣周期，需每次周期迭代時，都依據(jù)電網(wǎng)電壓反饋的實時頻率更新采樣頻率 fs，數(shù)字處理器需要較大存儲空間，實現(xiàn)復(fù)雜。故此，采用更新樣本周期點數(shù)N，即

實際采集到的電網(wǎng)頻率通常不是整數(shù)，即使為整數(shù)，經(jīng)式（7）運算，N也可能不為整數(shù)，需要進行取整運算，實現(xiàn)過程為

2.3 SDFT-QSG動態(tài)性能分析

經(jīng)前文分析，前級SDFT采樣點數(shù)N取128，SDFT-QSG對諧波及直流偏移分量的抑制能力強。由文獻[15]可知，系統(tǒng)動態(tài)響應(yīng)能力與樣本點數(shù)密切相關(guān)，樣本點數(shù)越大，系統(tǒng)動態(tài)響應(yīng)速度越慢，這與鎖相環(huán)快速精準的鎖相能力相悖，故需要驗證SDFT結(jié)構(gòu)動態(tài)性能。SOGI-QSG為常規(guī)正交信號發(fā)生器，具有較好的諧波抑制能力與動態(tài)響應(yīng)能力，但對直流偏移分量抑制不佳。對SDFT與SOGI進行階躍響應(yīng)對比，驗證SDFT的動態(tài)響應(yīng)能力、抑制直流偏移能力。仿真結(jié)果如圖4所示。結(jié)果表明，SOGI-QSG的α與qα信號動態(tài)響應(yīng)時間均優(yōu)于SDFT，但SDFT也在一個基波周期20 ms達到穩(wěn)態(tài)。SDFT-QSG的α與qα信號超調(diào)量較低，且穩(wěn)態(tài)輸出為零，能夠較好地濾除直流偏移分量；而SOGIQSG的qα信號穩(wěn)態(tài)輸出不為零，無法濾除直流偏移分量，其α與qα信號都存在穩(wěn)態(tài)誤差。綜合考量，SDFT-QSG動態(tài)性能更優(yōu)。

圖4 動態(tài)響應(yīng)對比Fig.4 Comparison of dynamic response

3 MAF正、負序分量分離設(shè)計

不平衡電網(wǎng)電壓存在負序分量，將SDFT-QSG提取到的正交信號代入正、負序分量計算模塊PNSC（positive negative sequence calculator）原則上可消除負序分量，但電網(wǎng)畸變程度高將導(dǎo)致SDFT頻譜泄漏嚴重，PNSC交叉分解矩陣會導(dǎo)致相位檢測出現(xiàn)嚴重波動，影響頻率精度，故設(shè)計滑動平均濾波方法進行正、負序分離。

3.1 MAF原理分析

滑動平均濾波器MAF（moving average filter）的離散域傳遞函數(shù)為[11]

式中：Nm為MAF采樣點數(shù)（也稱濾波器階數(shù)），Nm=Tωmfsm，Tωm為MAF窗口寬度，其大小既影響濾波性能，又影響系統(tǒng)動態(tài)響應(yīng)速度，Tωm越大，MAF濾波效果越好，但動態(tài)響應(yīng)速度越慢；fsm為滑動平均濾波器采樣頻率。

取 Nm=100，fsm=10 kHz，則Tωm=0.01 s，相應(yīng)的伯德圖如圖5所示。

圖5 MAF伯德圖Fig.5 Bode diagram of MAF

式中，Tsm為MAF采樣周期。

由式（10）及圖5可知，MAF可抑制特定頻次的諧波分量，但對直流分量無抑制作用。

3.2 電壓負序分量擾動分析

設(shè)αβ兩相靜止坐標系下的電壓表達式為

式中：U+、U-分別為正、負序基波電壓幅值；φ-為負序分量相對正序分量的初始相位角。

對式（11）進行Park變換，得dq兩相旋轉(zhuǎn)坐標系下的電壓表達式為

當(dāng)鎖相環(huán)鎖定電網(wǎng)電壓正序分量旋轉(zhuǎn)位置時，即θ=ωct，式（12）改寫為

由式（13）可知：經(jīng)坐標變換，正序分量轉(zhuǎn)化為直流分量，負序分量轉(zhuǎn)化為2倍頻分量?；冸娋W(wǎng)下，負序分量對鎖相環(huán)產(chǎn)生較大擾動，使其無法準確檢測出正序分量相位信息。

3.3 負序分量消除設(shè)計

由上述分析，欲滿足精準鎖相的條件是：uq=0，而MAF對直流分量無抑制作用，故不會干擾系統(tǒng)判斷。負序分量通過坐標變換轉(zhuǎn)化為2倍頻分量，設(shè)計的核心是衰減100 Hz處頻率，由圖5可知，在100 Hz處，增益為-160 dB，滿足既定要求。

MAF負序分量消除原理如圖6所示。

圖6 MAF負序分量消除結(jié)構(gòu)Fig.6 Structure of MAF negative-sequence component elimination

圖6中，基于MAF負序分量消除具體步驟如下：

（1）通過坐標變換將abc三相靜止坐標轉(zhuǎn)換為dq兩相同步旋轉(zhuǎn)坐標；

（2）將輸出的q軸分量經(jīng)過滑動平均濾波，濾除二倍頻分量；

（3）根據(jù)輸出的q軸信號進行SRF-PLL相位檢測控制，得到相角信息，并通過反饋環(huán)節(jié)為Park變換提供相位參考，通過PI控制調(diào)節(jié)，使直流分量為0，完成鎖相。

3.4 PI參數(shù)設(shè)計

建立電壓含擾動分量的小信號模型，如圖7所示。

圖7 MAF-PLL小信號模型結(jié)構(gòu)Fig.7 Structure of MAF-PLL small signal model

故 ξ＞1，根據(jù)參考文獻[17]，對鎖相環(huán)暫態(tài)過程與穩(wěn)定性進行綜合考量，本文取ξ=2.4，kp=83.33，ki=2 893.5。

3.5 基于相位補償MAF動態(tài)性能設(shè)計

為提高MAF的動態(tài)性能，需對MAF引起的相位延遲進行補償。實現(xiàn)這一目標的簡單方法是在MAF控制回路中加入相位超前補償器[18]。

借鑒文獻[18]補償器的傳遞函數(shù)，設(shè)計MAF傳遞函數(shù)的逆函數(shù)，以有效補償MAF引起的相位延遲。傳遞函數(shù)為

式中：kpc為傳遞函數(shù)幅值增益系數(shù)；r∈[0 1）為衰減因子。

帶有相位補償器的MAF，其衰減因子r取0.90、0.95、0.99的頻率響應(yīng)如圖8所示。r為0時，系統(tǒng)傳遞函數(shù)不變，r為1時，相位補償控制器成為全通濾波器。結(jié)合MAF傳遞函數(shù)，權(quán)衡動態(tài)響應(yīng)速度與濾除由負序分量引起的二倍頻振蕩，本文r取0.99，可有效補償由MAF引起的相位延遲。

圖8 基于相位補償?shù)腗AF伯德圖Fig.8 Bode diagram of MAF based on phase compensation

由圖8可知，r=0.99系統(tǒng)動態(tài)特性更好。

綜上分析，設(shè)計的鎖相環(huán)結(jié)構(gòu)如圖9所示。

圖9 CSDFT-MAF-PLL結(jié)構(gòu)Fig.9 Structure of CSDFT-MAF-PLL

4 仿真驗證與分析

為驗證所提方法正確性與可行性，在Matlab/Simulink中搭建仿真模型。將DSOGI-PLL、SDFTPLL和CSDFT-MAF-PLL進行對比分析驗證。為均衡DSOGI-PLL濾波性能和動態(tài)響應(yīng)時間，阻尼系數(shù)設(shè)為0.707。CSDFT-QSG前級端N為128，采樣頻率fs為12.8 kHz。系統(tǒng)參數(shù)為：額定電壓有效值Vg為220 V，額定頻率為50 Hz，初始相位為0°。仿真工況設(shè)置為4種，如表1所示，具體故障類型設(shè)置為5種，如表2所示。

表1 仿真工況設(shè)計Tab.1 Design of simulation conditions

表2 故障類型設(shè)計Tab.2 Design of fault types

4.1 工況一

在1.4～1.6 s，進行工況一故障測試，仿真結(jié)果如圖10所示。

圖10 工況一仿真波形Fig.10 Simulation waveforms in Condition 1

由圖10可見：DSOGI-PLL頻率超調(diào)大，峰值為59.76 Hz，在既定時間內(nèi)不能趨于穩(wěn)定；SDFT-PLL超調(diào)了3.62 Hz；CSDFT-MAF-PLL超調(diào)量僅為1.5 Hz。SDFT-PLL相位誤差最大達22.5°，說明PNSC結(jié)構(gòu)對SDFT有干擾。

4.2 工況二

在1.8～2.0 s，進行工況二故障測試，仿真結(jié)果如圖11所示。

圖11 工況二仿真波形Fig.11 Simulation waveforms in Condition 2

由圖11可見：相較于工況一，DSOGI-PLL波動劇烈程度降低，但峰值增量提升，說明DSOGI結(jié)構(gòu)對頻率跳變與直流偏移量敏感程度高；CSDFTMAF-PLL相位誤差最小達到-1°，比工況一提升約3.5°，故其對諧波注入較為敏感。

4.3 工況三

在2.2～2.4 s，進行工況三故障測試，仿真結(jié)果如圖12所示。

由圖12可見：DSOGI-PLL頻率動態(tài)峰值達到65 Hz，說明頻率跳變和多次諧波同時并存，對DSOGI模塊影響劇烈，而所設(shè)計的鎖相環(huán)超調(diào)量為0.27 Hz，并無明顯變化，說明所提策略魯棒性強。

圖12 工況三仿真波形Fig.12 Simulation waveforms in Condition 3

4.4 工況四

在2.6～2.8 s，進行工況四故障測試，仿真結(jié)果如圖13所示。

由圖13可見，工況四最為復(fù)雜，相較于工況三，CSDFT-MAF-PLL穩(wěn)態(tài)相位誤差約為0.059°，比工況三提高約0.02°，頻率變化不明顯。DSOGI-PLL頻率誤差大，且波動明顯。

圖13 工況四仿真波形Fig.13 Simulation waveforms in Condition 4

綜上分析：數(shù)字濾波器采樣數(shù)據(jù)較大，SDFTPLL及CSDFT-MAF-PLL動態(tài)響應(yīng)時間均滯后DSOGI-PLL，但在規(guī)定時限內(nèi)；對比工況一、三，DSOGI-PLL相位穩(wěn)態(tài)波動基本一致，但工況三下的相位誤差明顯高于工況一，可見DSOGI結(jié)構(gòu)對頻率跳變很敏感，而設(shè)計的CSDFT結(jié)構(gòu)降低了頻率跳變對頻譜泄漏的影響程度，相位誤差穩(wěn)定在±0.04°；工況二、四與工況一、三對比，直流偏移使DSOGI-PLL頻率與相位誤差無法在額定百分比下達到穩(wěn)態(tài)，而CSDFT結(jié)構(gòu)波形光滑，穩(wěn)態(tài)無誤差，諧波注入對其有一定影響，相位誤差增加，但在預(yù)期之內(nèi)。綜合4種工況，都發(fā)生幅值與相位跳變，系統(tǒng)存在負序分量，而仿真結(jié)果表明所設(shè)計的MAF負序分量消除模塊效果顯著。CSDFT-MAF-PLL鎖相動態(tài)響應(yīng)時間較快，并不因MAF的引入對系統(tǒng)造成很大延遲，說明所提相位補償控制器的有效性。相較于文獻[12]，穩(wěn)態(tài)精度及動態(tài)響應(yīng)速度均有明顯改善。

5 結(jié)論

針對SDFT-PLL在復(fù)雜電網(wǎng)條件下鎖相優(yōu)勢與不足展開研究，結(jié)合MAF優(yōu)良濾波性能，設(shè)計CSDFT-MAF-PLL。理論分析與仿真驗證得出以下結(jié)論：

（1）提出CSDFT-QSG結(jié)構(gòu)，顯著提升系統(tǒng)諧波抑制與直流偏移消除能力，并對頻率、幅值、相位跳變不敏感，以電網(wǎng)實時頻率作為后級SDFT滑窗寬度選取依據(jù)，從而減小SDFT因頻譜泄漏導(dǎo)致的頻率穩(wěn)態(tài)誤差；

（2）提出MAF負序分量消除設(shè)計，將三相電網(wǎng)電壓的正序分量送入SRF-PLL，負序分量轉(zhuǎn)化為二倍頻分量被抑制；

（3）提出相位補償控制器，克服MAF相位延遲的弊端，使所設(shè)計的鎖相環(huán)在各種工況下，達到穩(wěn)態(tài)的時間均保持在1.5個基波周期左右。

本文建立單機無窮大電網(wǎng)系統(tǒng)，進行了4種工況下鎖相并網(wǎng)性能研究，模型結(jié)構(gòu)相對單一。下一步研究可以結(jié)合雙饋風(fēng)電系統(tǒng)，開展鎖相環(huán)動態(tài)特性對系統(tǒng)慣性響應(yīng)及穩(wěn)定性分析。