CN102253282A - 纳托尔窗函数连续频谱内插电力谐波参数获取方法 - Google Patents
纳托尔窗函数连续频谱内插电力谐波参数获取方法 Download PDFInfo
- Publication number
- CN102253282A CN102253282A CN2011101549953A CN201110154995A CN102253282A CN 102253282 A CN102253282 A CN 102253282A CN 2011101549953 A CN2011101549953 A CN 2011101549953A CN 201110154995 A CN201110154995 A CN 201110154995A CN 102253282 A CN102253282 A CN 102253282A
- Authority
- CN
- China
- Prior art keywords
- frequency
- harmonic
- formula
- value
- window function
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Measuring Frequencies, Analyzing Spectra (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
纳托尔(Nuttall)窗函数连续频谱内插电力谐波参数获取方法适用于对电力网电压、电流的谐波分析与监测。首先使用线性调频Z变换(Chirp-ZTransform)或称CZT从含有谐波的电力信号中高精度提取基波信号参数(幅值、频率和相位)。然后从被分析电力信号中减去基波信号后加Nuttall窗函数截断电力信号,并用FFT计算出剩余信号的频谱。再依据基波频率精确计算出各谐波的频率值。最后依据各谐波频率对Nuttall窗函数在频域内插值,精确计算出各电力谐波的参数。本发明与加Nuttall窗FFT双谱线插值拟合分析电力谐波方法有基本同等的估值精度,而计算量约为1/2。
Description
技术领域
本发明涉及一种电网电压和电流波形畸变的分析和自动监测算法,可用于各种电网电压和电流波形畸变的分析仪器和自动监测装置。属于电力测量和自动化技术领域。
背景技术
随着电力电子技术和器件的发展,非线性负荷在电力系统中的应用越来越广泛,电力系统谐波污染日益严重,谐波已成为影响电能质量的主要问题。对谐波分量参数的高精度估计将有利于电能质量的评估和采取相应的必要治理措施。
快速傅立叶变换(FFT)是谐波分析最快捷的工具。但是, FFT精确分析频谱的前提是保证对信号的同步采样和整周期截断。实际电网频率通常在工频附近波动的,因此而造成非同步采样和非整数周期截断,这将产生频谱泄漏和谱间干扰,使谱分析产生误差。这一问题的解决通常有2条思路:一是通过锁相环技术(硬件或软件)来解决同步采样和整数周期截断问题。由于电网频率并非恒定值,而锁相环响应需要时间,因而不能保证完全同步采样。普遍采用的另一思路是通过选择谱能量主要集中在主瓣,旁瓣普能量小、且幅值衰减快的窗函数,以减小谱间干扰,即频谱的长范围泄漏;通过双谱线间插值修正,以减小栅栏效应,进而提高谐波估计精度。许多学者采用加窗插值法都有效地提高了谐波估计的精度[1~6],但随着插值修正曲线拟合函数的阶次增高及谐波含有次数的增多,谐波估计精度提高的同时计算量大量增加。本发明将提出谐波高精度估计的另一条思路,Nuttall窗函数连续频谱内插精确计算电力谐波参数方法。
参考文献:
[1] H. Xue and R. Yang,Optimal interpolating windowed discrete Fourier transform algorithms for harmonic analysis in power systems [J],IEE Proceedings of Generation, Transmission and Distribution, Vol.150, No.5, September 2003:583-587
[2] 庞浩,李东霞,俎云霄等,应用FFT进行电力统谐波分析的改进算法[J],中国电机工程学报,2003,23(6):50-54
[3] 曾博,滕召胜,高云鹏,王一,基于Rife-Vincent窗的高准确度电力谐波相量计算方法[J],电工技术学报,2009,24(8):154-158
[4] 曾博,滕召胜,温和,卿柏元,莱夫一文森特窗插值FFT谐波分析方法[J],中国电机工程学报,2009,29(10):115-120
[5] 卿柏元,滕召胜,高云鹏,温和,基于Nuttall窗双谱线插值FFT的电力谐波分析方法[J],中国电机工程学报,2008,28(25):153-157
[6] Reljin I, Reljin B, Papic V. Extremely flat-top windows for harmonic analysis [J]. IEEE Transactions on Instrumentation and Measurement, 2007, 56(3),1025-1041
发明内容
技术问题:本发明的目的是提供一种用于电网电压和电流波形畸变的分析和电
力谐波参数自动监测的纳托尔窗函数连续频谱内插电力谐波参数获取方法,可用于各种电网电压和电流波形畸变的分析仪器和自动监测装置。
技术方案:本发明的纳托尔(Nuttall)窗函数连续频谱内插电力谐波参数获取方法采用如下步骤:
这里,M为线性调频Z变换时在频域内的抽样点数;k ’ 为M个X(k)中取得
步骤b. 对Nuttall窗函数连续频谱在频域内抽样求得各次电力谐波的校正系数;
式4中:为数字角频率;为一预定值,,,为电网基波额定频率;Ts为采样周期,采样频率等于电网基波额定频率的2i倍,i取自然正整数,i=1、2、…;为以电网p次谐波信号的实际数字角频率在Nuttall窗函数连续频谱上在频域内抽样值,p为自然正整数;
有益效果:本发明提供了一种电力谐波分析的新方法。先用CZT高精度地估计出含有谐波的电力信号中的基波参数(幅值、频率和相位),进而,通过Nuttall窗函数频域内插值的方法精确计算出各次谐波的参数。通过同一电力信号的仿真试验证明,它与各种现有的加Nuttall窗插值估计电力谐波的分析方法有基本等同的估计精度,而本文提出的电力谐波估计方法不需要对各次谐波进行双谱线内插反拟合计算,计算量约为现有的这些算法的计算量的1/2,在计算速度上有明显的优势。因此,是一种很有实用价值的电力谐波高精度估计的方法。
附图说明
图1 快速CZT算法流程。
图2同步采样和整周期截断后的频谱。
图 3非同步采样和非整周期截断后的频谱。
具体实施方式
为了实现上述目的, 本发明的实施可直接用分压器或从电压互感器PT二次侧取得电网的母线电压信号、从电流互感器CT取得电流信号,经过适当的信号调理后送达信号采样入口。
步骤a.采样被分析电力信号电压或电流,并按图1流程计算其线性调频Z变换(Chirp-Z Transform)或称CZT值,取自然正整数,再由式1、式2和式3分别计算出电力信号的基波参数,幅值、频率和相位;
估计基波幅值: 式1
设被分析电力信号为:
采样被分析电力信号得:
Ts为采样周期。
根据标准,电网最大允许频率偏差为0.5Hz,因此,可设待分析频段的起始频率为θ=2×49.5,终止频率为θ+(M-1)φ=2×50.5,即分析带宽为50.5-49.5=1Hz。如果信号采样频率为f s =6400Hz,若取M=1280(10个周波),则频率采样间隔φ=2/639=2×7.82473×10-4Hz。若取模最大值所对应的频率作为基波频率估计值f 1 ,频率测量的误差为|△f|≤(φ/2)/2=3.90930×10-4Hz。按照以上方法确定了计算参数后,可按图1计算流程,由以下步骤计算得基波的频率f 1 、幅值A 1 和相位:
步骤b. 对Nuttall窗函数连续频谱在频域内抽样求得各次电力谐波的校正系数;
式4中:为数字角频率;为一预定值,,,为电网基波额定频率,如图2所示;Ts为采样周期,采样频率等于电网基波额定频率的2i倍,i取正整数,i=1、2、…; 为以电网p次谐波信号的实际数字角频率,p为自然正整数,也在Nuttall窗函数连续频谱上在频域内抽样值,如图3所示;
Nuttall窗函数为:
可预定以电网的基波额定频率的2i倍(i取自然数,i=1、2、…)采样,即采样频率。这样,电网基波和谐波的额定频率是与Nuttall窗函数同步采样并能整周期截断的。则等于p次电力谐波信号与Nuttall窗函数同步采样并整周期截断的数字角频率,则以角频率在Nuttall窗函数连续频谱上在频域内抽样求得:
Claims (1)
1. 一种纳托尔窗函数连续频谱内插电力谐波参数获取方法,其特征是该方法采用如下步骤:
估计基波相位值: 式3
这里,M为线性调频Z变换时在频域内的抽样点数;k ’ 为M个X(k)中取得
式4中:为数字角频率;为一预定值,,,为电网基波额定频率;Ts为采样周期,采样频率等于电网基波额定频率的2i倍,i取自然正整数,i=1、2、…;为以电网p次谐波信号的实际数字角频率在Nuttall窗函数连续频谱上在频域内抽样值,p为自然正整数;
P次谐波的相位为: 式6
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110154995.3A CN102253282B (zh) | 2011-06-10 | 2011-06-10 | 纳托尔窗函数连续频谱内插电力谐波参数获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110154995.3A CN102253282B (zh) | 2011-06-10 | 2011-06-10 | 纳托尔窗函数连续频谱内插电力谐波参数获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102253282A true CN102253282A (zh) | 2011-11-23 |
CN102253282B CN102253282B (zh) | 2014-07-16 |
Family
ID=44980654
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110154995.3A Expired - Fee Related CN102253282B (zh) | 2011-06-10 | 2011-06-10 | 纳托尔窗函数连续频谱内插电力谐波参数获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102253282B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102539915A (zh) * | 2012-01-06 | 2012-07-04 | 中国矿业大学 | 时延傅立叶变换测频法精确计算电力谐波参数方法 |
CN105137198A (zh) * | 2015-09-15 | 2015-12-09 | 福州大学 | 一种基于Nuttall窗-五点变换FFT的介损测量新方法 |
CN108107393A (zh) * | 2017-12-06 | 2018-06-01 | 中国矿业大学 | 一种谐波分析中谱峰可信度判断方法 |
CN108107269A (zh) * | 2017-12-07 | 2018-06-01 | 中国矿业大学 | 一种谐波分析中幅值求解方法 |
CN110988469A (zh) * | 2019-11-14 | 2020-04-10 | 中国矿业大学 | 一种快速谐波检测方法 |
CN111693774A (zh) * | 2020-05-06 | 2020-09-22 | 南方电网科学研究院有限责任公司 | 一种输电网的谐波测量方法和装置 |
CN112485522A (zh) * | 2020-12-09 | 2021-03-12 | 国网四川省电力公司电力科学研究院 | 基于电能数据感知的平顶窗函数同步相量测量方法及装置 |
CN114659791A (zh) * | 2022-02-28 | 2022-06-24 | 广东机电职业技术学院 | 汽轮机故障检测方法、系统、装置和存储介质 |
CN114705913A (zh) * | 2022-04-12 | 2022-07-05 | 上海赢双电机有限公司 | 一种旋转变压器的谐波分析方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101113995A (zh) * | 2007-08-29 | 2008-01-30 | 湖南大学 | 基于Nuttall窗双峰插值FFT的基波与谐波检测方法 |
CN101701984A (zh) * | 2009-11-23 | 2010-05-05 | 浙江大学 | 基于三项系数Nuttall窗插值FFT的基波与谐波检测方法 |
CN101718816A (zh) * | 2009-11-23 | 2010-06-02 | 浙江大学 | 基于四项系数Nuttall窗插值FFT的基波与谐波检测方法 |
-
2011
- 2011-06-10 CN CN201110154995.3A patent/CN102253282B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101113995A (zh) * | 2007-08-29 | 2008-01-30 | 湖南大学 | 基于Nuttall窗双峰插值FFT的基波与谐波检测方法 |
CN101701984A (zh) * | 2009-11-23 | 2010-05-05 | 浙江大学 | 基于三项系数Nuttall窗插值FFT的基波与谐波检测方法 |
CN101718816A (zh) * | 2009-11-23 | 2010-06-02 | 浙江大学 | 基于四项系数Nuttall窗插值FFT的基波与谐波检测方法 |
Non-Patent Citations (2)
Title |
---|
H. XUE等: "Optimal interpolating windowed discrete Fourier trabsform algorithms for harmonic analysis in power system", 《IEE PROC.-GENER. TRANSM. DISTRIB》, vol. 150, no. 5, 30 September 2003 (2003-09-30) * |
卿柏元等: "基于Nuttall窗双谱线插值FFT的电力谐波分析方法", 《中国电机工程学报》, vol. 28, no. 25, 5 September 2008 (2008-09-05), pages 153 - 158 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102539915A (zh) * | 2012-01-06 | 2012-07-04 | 中国矿业大学 | 时延傅立叶变换测频法精确计算电力谐波参数方法 |
CN105137198A (zh) * | 2015-09-15 | 2015-12-09 | 福州大学 | 一种基于Nuttall窗-五点变换FFT的介损测量新方法 |
CN108107393A (zh) * | 2017-12-06 | 2018-06-01 | 中国矿业大学 | 一种谐波分析中谱峰可信度判断方法 |
CN108107393B (zh) * | 2017-12-06 | 2019-11-05 | 中国矿业大学 | 一种谐波分析中谱峰可信度判断方法 |
CN108107269A (zh) * | 2017-12-07 | 2018-06-01 | 中国矿业大学 | 一种谐波分析中幅值求解方法 |
CN108107269B (zh) * | 2017-12-07 | 2019-11-05 | 中国矿业大学 | 一种谐波分析中幅值求解方法 |
CN110988469A (zh) * | 2019-11-14 | 2020-04-10 | 中国矿业大学 | 一种快速谐波检测方法 |
CN110988469B (zh) * | 2019-11-14 | 2021-03-16 | 中国矿业大学 | 一种快速谐波检测方法 |
CN111693774A (zh) * | 2020-05-06 | 2020-09-22 | 南方电网科学研究院有限责任公司 | 一种输电网的谐波测量方法和装置 |
CN112485522A (zh) * | 2020-12-09 | 2021-03-12 | 国网四川省电力公司电力科学研究院 | 基于电能数据感知的平顶窗函数同步相量测量方法及装置 |
CN114659791A (zh) * | 2022-02-28 | 2022-06-24 | 广东机电职业技术学院 | 汽轮机故障检测方法、系统、装置和存储介质 |
CN114659791B (zh) * | 2022-02-28 | 2023-07-04 | 广东机电职业技术学院 | 汽轮机故障检测方法、系统、装置和存储介质 |
CN114705913A (zh) * | 2022-04-12 | 2022-07-05 | 上海赢双电机有限公司 | 一种旋转变压器的谐波分析方法 |
CN114705913B (zh) * | 2022-04-12 | 2024-01-09 | 上海赢双电机科技股份有限公司 | 一种旋转变压器的谐波分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102253282B (zh) | 2014-07-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102338827B (zh) | 用于电网波形畸变分析和电力谐波参数自动监测方法 | |
CN102331526B (zh) | 汉宁窗函数连续频谱内插电力谐波参数获取方法 | |
CN102253282B (zh) | 纳托尔窗函数连续频谱内插电力谐波参数获取方法 | |
CN102539915B (zh) | 时延傅立叶变换测频法精确计算电力谐波参数方法 | |
Zhan et al. | A Clarke transformation-based DFT phasor and frequency algorithm for wide frequency range | |
CN105223418B (zh) | 次同步和超同步谐波相量的测量方法及测量装置 | |
Ren et al. | Real-time power system frequency and phasors estimation using recursive wavelet transform | |
CN103454497B (zh) | 基于改进加窗离散傅立叶变换的相位差测量方法 | |
CN103995181B (zh) | 数字化变电站电能质量谐波分析方法 | |
Zhang et al. | Analysis of white noise on power frequency estimation by DFT-based frequency shifting and filtering algorithm | |
CN107247182B (zh) | 一种基于量测相量数据的间谐波分量还原方法 | |
CN102818930B (zh) | 一种高精度快速计算电力谐波参数的方法 | |
WO2016138764A1 (zh) | 一种提升pmu同步相量测量精度的方法 | |
CN103353558A (zh) | 一种电能质量监测方法 | |
CN102520245A (zh) | 基于三次样条插值波形重构的微网谐波及间谐波分析方法 | |
CN103869162A (zh) | 一种基于时域准同步的动态信号相量测量方法 | |
CN103575984A (zh) | 基于凯塞窗双谱线插值fft的谐波分析方法 | |
CN108490257B (zh) | 一种基于频谱曲线拟合的短时窗间谐波测量方法 | |
CN104020352A (zh) | 一种适用于m类pmu单元的同步相量测量方法 | |
Rivas et al. | Adaptive current harmonic estimation under fault conditions for smart grid systems | |
Zhang et al. | Shifting window average method for phasor measurement at offnominal frequencies | |
CN109507480A (zh) | 一种邻近基波/谐波的间谐波检测方法和装置 | |
CN106645919A (zh) | 基于三相瞬时功率的电网全谱功率振荡相量同步测量方法 | |
CN102495285B (zh) | 对称窗函数功率重心估计电力谐波参数的方法 | |
Abdolkhalig et al. | Phasor measurement based on IEC 61850-9-2 and Kalman–Filtering |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140716 Termination date: 20170610 |
|
CF01 | Termination of patent right due to non-payment of annual fee |