李芊芊,張江濱,李亞琳
(1.西安理工大學(xué)水利水電學(xué)院,陜西西安 710048;2.國網(wǎng)陜西省電力公司,陜西西安 710048)
適用于雙饋風(fēng)機的并網(wǎng)電壓波動與閃變檢測系統(tǒng)的研究
李芊芊1,張江濱1,李亞琳2
(1.西安理工大學(xué)水利水電學(xué)院,陜西西安710048;2.國網(wǎng)陜西省電力公司,陜西西安710048)
風(fēng)電因電能質(zhì)量問題影響使其發(fā)展受到制約,電壓波動及閃變是其中之一。使用目前的閃變儀對風(fēng)電機組并網(wǎng)引起電壓閃變進行檢測,發(fā)現(xiàn)在檢測低頻段波動引起的閃變時誤差較大,而這也是風(fēng)電并網(wǎng)產(chǎn)生閃變的主要頻段。因此,如何使低頻段檢測誤差減小是該文研究的意義。根據(jù)IEC標準,在Matlab平臺建立數(shù)字閃變檢測系統(tǒng)仿真模型,并利用雙線性方法求得系統(tǒng)的參數(shù)。針對上述情況,分析其誤差產(chǎn)生的原因,通過對檢測的瞬時閃變視感度S(t)引入?yún)?shù)校正以減小低頻檢測誤差,使其滿足風(fēng)電并網(wǎng)閃變值檢測的要求。最后將修正系數(shù)后的檢測系統(tǒng)運用在一個雙饋風(fēng)力發(fā)電系統(tǒng)實例中,證明該風(fēng)電場滿足我國電壓閃變值要求。
SIMULINK;雙饋風(fēng)力發(fā)電機;電壓波動;閃變檢測
風(fēng)力發(fā)電在飛速發(fā)展的同時,由于輸出功率的波動給接入電網(wǎng)帶來一定影響,使其發(fā)展受到制約。在風(fēng)機運轉(zhuǎn)過程當中,風(fēng)速的變化、塔影效應(yīng)等各種情況下可能會造成電壓波動;葉輪旋轉(zhuǎn)過程當中的轉(zhuǎn)矩處于非穩(wěn)定狀態(tài)、風(fēng)機起停和發(fā)電機的切換,會導(dǎo)致輸出功率波動,造成電壓波動和電壓閃變[1]。對電壓閃變進行精確檢測能夠更好地了解電能質(zhì)量狀況,以便采取合理的措施來控制或抑制閃變。因此,選擇合適的檢測方法,設(shè)計出規(guī)范化和標準化的閃變檢測系統(tǒng)具有實際的研究價值和意義[2]。目前,國外對于電壓閃變儀都有一定程度的研究,如日本有△V10檢測儀、法國有EDF檢測儀、德國有FGH檢測儀等[3],但這些閃變檢測儀主要適用于檢測電弧爐等較高頻率波動造成的閃變(0~25 Hz),并不適用于風(fēng)電并網(wǎng)引起的低頻率波動造成的閃變。
文獻[4]提出研制一種高精度閃變儀,提出了提高閃變檢測精度的措施與統(tǒng)計環(huán)節(jié)、采樣頻率和采樣設(shè)備有關(guān),但要以高精度的采樣設(shè)備和大量的采樣數(shù)據(jù)為前提,計算量大,成本高;文獻[5]提出IEC閃變儀在進行風(fēng)電場閃變檢測中誤差較大,采用校正環(huán)節(jié)校正閃變測試系統(tǒng)的輸出能夠提高精度,但文中只介紹了模擬閃變測試系統(tǒng)的仿真過程與校正方法,沒有提到數(shù)字閃變測試系統(tǒng)。文獻[6]介紹了閃變測量系統(tǒng)的數(shù)字化實現(xiàn)方法,對不同頻率、幅值的調(diào)幅波引起的波動進行仿真,但對低頻段測量的較大誤差沒有提出具體修正方法。文獻[7]提出了離散化的方法用于閃變計算,其低頻段測量誤差小,適用于風(fēng)電閃變的檢測,但高頻段測量誤差較大。
在上述文獻研究成果的基礎(chǔ)上,本文依據(jù)IEC推薦的閃變檢測系統(tǒng)原理,在Matlab平臺下搭建檢測模型,利用Bilinear函數(shù)得到數(shù)字閃變檢測系統(tǒng)各環(huán)節(jié)參數(shù),并對特定幅值和頻率的調(diào)制波進行閃變的檢測。針對低頻段檢測誤差較大的情況,分析其誤差產(chǎn)生的原因,通過參數(shù)校正的方法減小低頻檢測誤差,使其滿足風(fēng)電并網(wǎng)閃變值檢測的要求。
1.1平方檢測法檢測原理
目前檢測電壓波動有許多種方法,由于平方法檢測易于實現(xiàn)數(shù)字化,也被IEC推薦的閃變檢測原理所采用,因此本文選取平方檢測法。其中平方檢測法流程圖如圖1所示。
圖1 平方檢測法流程圖Fig.1 Square test flow chart
圖1中框1為用于調(diào)幅波解調(diào)的平方器???為用于模擬人眼對照度波動的頻率最大感知范圍,在0.05~35 Hz之間的帶通濾波器;框3為模仿人眼對頻率選擇的加權(quán)濾波器???、框5為模仿視神經(jīng)反映與記憶效應(yīng)的平方器以及一階平滑濾波???為統(tǒng)計和分析,對前面得到的瞬時視感度值S(t)采樣后進行統(tǒng)計評定,最后計算得出短時閃變值Pst。Pst計算方法主要有分段線性差值法、多項式曲線擬合法以及統(tǒng)計排序法等,可通過Matlab編程實現(xiàn)。
高通濾波器:
視感度加權(quán)濾波器:
一階低通濾波器,其中τ=300 ms:
1.3數(shù)字閃變測試系統(tǒng)的設(shè)計
為了測量精度和穩(wěn)定性的提高,需要將模擬閃變檢測系統(tǒng)轉(zhuǎn)換成數(shù)字閃變檢測系統(tǒng),轉(zhuǎn)換的方法多采用雙線性變化法,它實際上是由s域到z域的映射,體現(xiàn)在傳遞函數(shù)上即為由H(s)到H(z)的變換,是復(fù)變函數(shù)的映射變換。
選擇采樣頻率為400 Hz,利用雙線性變換法將式由s域化為z域,公式如下:
式中:n,d為數(shù)字濾波器中的分子系數(shù)和分母系數(shù);N,D為模擬濾波器中的分子系數(shù)和分母系數(shù)。變換時的采樣頻率為fs。公式中各環(huán)節(jié)系數(shù)按照降冪的順序排列。
取采樣頻率fs=400 Hz,利用Bilinear函數(shù)對數(shù)字閃變檢測系統(tǒng)各環(huán)節(jié)濾波器參數(shù)進行計算。為了保證數(shù)字閃變檢測系統(tǒng)的精確度,利用Matlab中的format long語句,對變換后的參數(shù)設(shè)置15位有效數(shù)字,計算結(jié)果如下:
35 Hz低通濾波器參數(shù)Nd
0.05 Hz高通濾波器參數(shù):
視感度加權(quán)濾波器參數(shù):
300 ms延時平滑濾波器參數(shù):
依據(jù)計算得到的數(shù)字濾波器的系數(shù),搭建出閃變數(shù)字檢測系統(tǒng)仿真模型,如圖2所示。
系數(shù)的值K取為1/0.002 5,系統(tǒng)輸出窗口顯示的是S(t)瞬時視感度的波形。
圖2 數(shù)字閃變檢測系統(tǒng)模型Fig.2 Digital flicker test system
2.1模型檢驗
為了檢驗所搭建的數(shù)字閃變檢測模型的參數(shù)是否正確,首先需要模擬含有某一波動頻率的電壓,再對其進行閃變檢測。選擇幅值為0.002 5 V,頻率為8.8 Hz的正弦調(diào)幅波以及幅值為220 V的工頻電壓進行調(diào)制,將調(diào)制波作為閃變檢測系統(tǒng)的輸入波形,構(gòu)造函數(shù)如下:
上述模型得到的仿真電壓波形如圖3—8所示。
圖3 8.8 Hz正弦調(diào)幅波Fig.3 8.8 Hz sine modulated wave
圖4 調(diào)制波形Fig.4 Modulation waveform
圖5 通過帶通濾波器的波形Fig.5 Waveform through band-pass filter
圖6 通過視感度加權(quán)濾波器的波形Fig.6 Waveform through weighted filter
圖7 S(t)波形Fig.7 The S(t)waveform
圖8 S(t)微觀波形Fig.8 The S(t)micro waveform
圖3為幅值0.002 5 V,頻率8.8 Hz的正弦調(diào)幅波波形;圖4為調(diào)幅波與220 V工頻調(diào)制后的波形;圖5為通過帶通濾波器的波形;圖6為通過視感度加權(quán)濾波器的波形;圖7為瞬時閃變視感度S(t)波形;圖8為瞬時閃變視感度S(t)放大后微小幅值范圍內(nèi)的波形。
如圖6所示,輸出波形的頻率和幅值是8.8 Hz,121 V,由此證明,經(jīng)過加權(quán)濾波后,輸出波形的頻率等于調(diào)幅波頻率,輸出波形的幅值等于調(diào)幅波幅值和工頻電壓幅值平方的乘積;如圖7所示,數(shù)字閃變檢測系統(tǒng)最終檢測的瞬時閃變視感度S(t)為1左右;如圖8所示,S(t)基本處于1上下3%波動,滿足IEC對閃變檢測系統(tǒng)檢測精度要求的范圍。由仿真波形可知數(shù)字閃變測試系統(tǒng)各環(huán)節(jié)參數(shù)的選擇基本正確。
2.2誤差修正
IEC標準通過大量試驗,得到S(t)=1相應(yīng)的不同頻率的電壓波動d(%),利用這些數(shù)據(jù)進行仿真,對所建立的數(shù)字閃變測試系統(tǒng)模型進行校驗。結(jié)果表明,在2.5~20 Hz的頻率之間,該數(shù)字檢測系統(tǒng)輸出的S(t)值在1±0.05之間,其相對誤差小于±5%。說明S(t)所設(shè)計的閃變測試系統(tǒng)效果良好。但在低頻段,仿真輸出的S(t)值偏差較大,見表1。
表1 S(t)=1下的各頻率正弦電壓波動值及該波動下對應(yīng)的S(t)計算值Tab.1 Sine voltage fluctuations under different frequencies when S(t)=1 and the corresponding S(t)calculated value
由于在低頻段中,閃變數(shù)字檢測系統(tǒng)輸出的S(t)具有比較大的誤差,而這一頻段正好位于風(fēng)電引起的電壓波動的頻段(3P頻率范圍一般在1~3 Hz之間),超出檢測精度標準。通過大量試驗發(fā)現(xiàn)檢測誤差主要來自視感度加權(quán)濾波環(huán)節(jié)。因此,本文重點研究提高閃變檢測系統(tǒng)精度的方法。本文采用最小二乘法擬合修正系數(shù),根據(jù)計算,得到各頻率下對應(yīng)波動幅度的調(diào)幅波檢測的S(t)所對應(yīng)的修正系數(shù)k值,如表2所示。
按照該表數(shù)據(jù)求出的校正曲線,見圖9。
通過運用最小二乘法的方法對表2中的K值進行數(shù)據(jù)擬合,得到式(7)如下所示,所對應(yīng)的多項式各項系數(shù)見表3。
表2 S(t)=1的各低頻率范圍內(nèi)S(t)檢測值及校正系數(shù)KTab.2 The(t)values and correction coefficient K under the low frequency ranges when(t)=1
圖9 對應(yīng)0.5~4.5 Hz范圍的S(t)校正曲線Fig.9 The modulation curve for S(t)within frequency 0.5~4.5 Hz
表3 擬合多項式的系數(shù)Tab.3 The coefficients of the polynomial
校正后的S(t)誤差曲線見圖10。
圖10 校正前后S(t)誤差曲線對比Fig.10 S(t)error curve comparison before and after correction
由此圖10可知,通過添加校正參數(shù)環(huán)節(jié),能夠有效減小閃變檢測系統(tǒng)在調(diào)幅波處于低頻段時的測量誤差,適用于風(fēng)電并網(wǎng)電壓閃變的檢測。
以某個裝機總量為9 MW的風(fēng)電場作為模型,如圖11所示。風(fēng)電場由6臺1.5 MW雙饋風(fēng)力發(fā)電機組組成,發(fā)電機出口電壓為690 V。將出口電壓作為數(shù)字閃變檢測系統(tǒng)的輸入電壓。
經(jīng)過校正系數(shù)后的數(shù)字閃變檢測得到短時閃變視感度曲線,如圖12所示,得知S(t)符合標準低頻段誤差5%的要求,因此,理論上該風(fēng)電場符合風(fēng)電引起的電壓閃變限值,可以接入電力系統(tǒng)。
圖11 雙饋風(fēng)機并網(wǎng)模型Fig.11 Doubly-fed fan grid model
利用Matlab/Simulink軟件搭建IEC標準的數(shù)字化檢測系統(tǒng),通過仿真得到以下結(jié)論:
1)當調(diào)幅波在高頻率范圍,S(t)值的檢測誤差較??;當調(diào)幅波在低頻率范圍,S(t)值的檢測誤差則較大。
2)通過分析可知誤差來源于視感度加權(quán)濾波環(huán)節(jié),添加參數(shù)校正環(huán)節(jié)可以提高低頻段檢測精度。
圖12 短時閃變視感度S(t)曲線Fig.12 Short-term flicker visual sensitivity S(t)curve
3)風(fēng)力發(fā)電引起的閃變頻率正好位于低頻段,以一個雙饋機并網(wǎng)風(fēng)電場為例,利用所設(shè)計的數(shù)字閃變儀進行檢測,通過仿真曲線顯示,該風(fēng)電場滿足國標對風(fēng)機并網(wǎng)電壓閃變限值的要求。
[1]肖湘寧.電能質(zhì)量分析與控制[M].北京:中國電力出版社,2004.
[2]孫樹勤.電壓波動與閃變[M].北京:中國電力出版社,1996.
[3]凌昊.電壓閃變檢測儀的研究[D].哈爾濱工業(yè)大學(xué),2010.
[4]鄧志.一種高精度閃變分析儀[J].智能電網(wǎng),2014,2(8):36-39.DENG Zhi.A high precision flickermeter[J].Smart Grid,2014,2(8):36-39(in Chinese).
[5]趙海翔,陳默子,戴慧珠,等.風(fēng)電并網(wǎng)引起閃變的測試系統(tǒng)仿真[J].太陽能學(xué)報,2005,26(1):28-33.ZHAO Haixiang,CHEN Mozi,DAI Huizhu.The test system simulation of flicker caused by wind power grid[J].Acta Energiae Solaris Sinica,2005,26(1):28-33(in Chinese).
[6]趙琳,楊乃琪,張秀峰.基于MATLAB/SIMULINK的電壓波動及閃變測量研究[J].自動化與儀器儀表,2012(05):30-32.ZHAO Lin,YANG Naiqi,ZHANG Xiufeng.Based on the MATLAB/SIMULINK of voltage fluctuation and flicker measurement research[J].Automation and Instrumentation,2012(05):30-32(in Chinese).
[7]賈秀芳.IEC閃變儀誤差分析及改進設(shè)計[J].電工技術(shù)學(xué)報,2006,21(11):121-126.JIA Xiufang.Error analyzing of IEC flickermeter and its improving design[J].Transactions of China Electrotechnical Society,2006,21(11):121-126(in Chinese).
[8]趙海翔,陳默子,戴慧珠.閃變值離散化計算的誤差分析及其應(yīng)用[J].電網(wǎng)技術(shù),2004,28(13):84-87.ZHAO Haixiang,CHEN Mozi,DAI Huizhu.Error analysis of discrete callculation method of flicker severity and its application[J].Power System Technology,2004,28(13):84-87(in Chinese).
[9]馬玉龍,劉連光,張建華,等.IEC閃變測量原理的數(shù)字化實現(xiàn)方法[J].中國電機工程學(xué)報,2001,21(11):92-95.MA Yulong,LIU Lianguang,ZHANG Jianhua,et al.Research of degital flickermeter based on IEC standard[J].Proceedings of the CSEE,2001,21(11):92-95(in Chinese).
[10]丁然.基于Matlab/Simulink的電壓閃變儀仿真[J].中國科技信息,2013(8):130-131,144.DING Ran.A simulation of flickermeter based on matlab/simulink[J].China Science and Technology Information,2013(8):130-131,144(in Chinese).
[11]張全明,劉會金,蘭泉妮,等.基于頻譜分析的間諧波閃變效應(yīng)計算[J].電力系統(tǒng)自動化,2009,33(9):67-71.ZHANG Quanming,LIU Huijin,LAN Quanni,et al.Based on the flicker effect harmonics spectrum analysis of the calculation[J].Automation of Electric Power Systems,2009,33(9):67-71(in Chinese).
[12]趙海翔.風(fēng)電引起的電壓波動和閃變研究[D].北京:中國電力科學(xué)研究院,2004.
[13]孫濤.風(fēng)力發(fā)電引起的電壓波動和閃變[J].電網(wǎng)技術(shù),2003,27(12):62-66,70.SUN Tao.Voltage fluctuation and flicker caused by wind power generation[J].Power System Technology,2003,27(12):62-66,70(in Chinese).
[14]吳義純.風(fēng)電引起的電壓波動與閃變的仿真研究[J].電網(wǎng)技術(shù),2009,33(20):125-130.WU Yichun.Simulation study on voltage fluctuations and flicker caused by wind farms[J].Power System Technology,2009,33(20):125-130(in Chinese).
(編輯李沈)
Research on Grid-Connected Voltage Fluctuation and Flicker Detection System Suitable for Doubly-Fed Fan
LI Qianqian1,ZHANG Jiangbin1,LI Yalin2
(1.Institute of Water Resources and Hydro-Electric Engineering,Xi'an University of Technology,Xi'an 710048,Shaanxi,China;2.State Grid Shaanxi Electric Power Company,Xi'an 710048,Shaanxi,China)
The development of wind power is constrained by its power quality,and the influence of voltage fluctuation and flicker is one of factors which affect the power quality.At present,the flicker meter is often used to detect the voltage flicker which is caused when wind turbines are connected to the power grid.It is found the error is bigger when the flicker caused by the low frequency fluctuation is detected,and this is also the main frequency band where flickers occur when wind turbines are connected to the power grid.Therefore,this paper aims at how to reduce the low frequency band detection error.According to the IEC standard,the digital flicker detection system simulation model is built in the Matlab platform,using the bi-linear method to obtain the parameters of the system.According to the above situation,the paper analyzes the causes of the error,and based on the instantaneous flicker visual sensitivity in detection of S(t)correction parameters are introduced to reduce the low frequency error so that the requirements for the wind power grid flicker value detection are met.Finally,the detection system with coefficients corrected is applied in the actual case of the doubly-fed wind power generation system,and it has proved that the wind farm satisfies the requirement of the voltage flicker value in our country.
SIMULINK;doubly-fed wind generator;voltage fluctuation;flicker detection
1674-3814(2015)12-0126-06
TM933
A
2015-09-16。
李芊芊(1990—),女,西安理工大學(xué)碩士研究生,從事電力系統(tǒng)測量、保護與控制方向研究;
張江濱(1956—),男,教授,研究方向為電力系統(tǒng)運行與控制方向。