CN106802368A - 一种基于频域插值的广域电网相量测量方法 - Google Patents
一种基于频域插值的广域电网相量测量方法 Download PDFInfo
- Publication number
- CN106802368A CN106802368A CN201710043642.3A CN201710043642A CN106802368A CN 106802368 A CN106802368 A CN 106802368A CN 201710043642 A CN201710043642 A CN 201710043642A CN 106802368 A CN106802368 A CN 106802368A
- Authority
- CN
- China
- Prior art keywords
- frequency
- power network
- network signal
- signal
- delta
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R25/00—Arrangements for measuring phase angle between a voltage and a current or between voltages or currents
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
本发明公开了一种基于频域插值的广域电网相量测量方法,首先,对电网信号进行采样并利用三阶最大旁瓣衰减窗对采样信号进行加权,然后对加权后的信号做离散傅里叶变换得到频谱函数,通过三谱线插值算法实现频率偏移值的计算,最后得到电网信号的频率,幅值和相角,实现电网信号相量的测量。该方法可以确定采样频率、窗函数等对测量结果的影响,能有效提高电网信号检测精度,而且计算量小。
Description
技术领域
本发明涉及信号测量领域,具体是一种基于频域插值的广域电网相量测量方法。
背景技术
随着电力系统的快速发展,电力系统的规模越来越大,运行的环境也越来越复杂,维持电力网络的安全、健康、可持续发展成为一个亟待解决的问题。为解决复杂、庞大的电力网络日益突出的问题,电力系统需要运用创新的方法和技术手段来确保电网的运行和控制能力,保证电网能在一个安全稳定的环境中运行。电网相量测量实际上是对电网信号的频率、幅值和相角这三个特征量的测量。电网向量作为反映电力系统运行特性的重要参数,及时、准确地测量系统的相量参数,可以预测系统是否将失去稳定,并通过切机、切负荷控制等保证系统的安全运行,防止事故进一步的扩大,从而导致系统的崩溃。因此实现电网相量的快速准确测量是非常必需且意义深远的。
早期的相量测量方法是将交流信号的波形送到控制中心,进行表面的分析和比较,这种方法不严谨且存在时延,测量精度低。随着国内网学者对相量测量的深入研究,已提出了很多电压和电流的相量测量方法,如:过零检测法、基于泰勒展开模型的算法、基于神经网络的算法、最小二乘多项式曲线拟合法、解析法、FFT等。其中,FFT算法因其易于微机实现而被广泛运用。而本专利就是选择的FFT算法。然而当被测信号频率发生变化导致非同步采样或对信号进行非整数周期截断时,FFT存在栅栏效应和频谱泄露现象,使得相量测量的误差较大。为此,可以选择性能优良的窗函数及对计算结果进行插值修正可减小非同步采样或数据非整周期截断引起的误差。
发明内容
本发明的目的是为了针对现有技术存在的不足,提供了一种简单易懂、精度更高、计算速度更快的相量测量方法。
本发明提出的解决方案为:首先,对电网信号进行采样并用三阶最大旁瓣衰减窗对采样信号进行加权,然后对加权后的信号做离散傅里叶变换得到频谱函数,通过三谱线插值算法实现频率偏移值的计算,最后得到电网信号的频率,幅值和相角,实现电网信号相量的测量。
本发明提出的基于频域插值的广域电网相量测量方法包括以下步骤:用采样频率fs对电网信号进行采样,得到长度为M的离散采样序列x(m);利用三阶最大旁瓣衰减窗w(m)对离散化的电网信号进行加权得到加权序列xw(m)=x(m)w(m),然后对加权序列做离散傅里叶变换得到频谱函数Xw(m),通过三谱线插值算法计算信号频率偏移值δ,由频率偏移值计算电网信号的频率f0,幅值A0和相角实现电网相量的测量。
所述的方法,关于频率分量偏移值δ的求解,由频谱函数其中函数W(m)的计算表达式为通过搜索|Xw(m)|,其中||表示取模值,得到幅值最大的谱线的位置l及其左右两根谱线位置,即l-1和l+1,由插值表达式η:继而得到频率偏移量δ的计算式:
进而可计算电网信号的频率为
所述的方法,电网信号的幅值A0的计算式为:
所述的方法,电网信号的相角的计算式为:
其中Phase{}是求相角运算。
本发明基于频域插值的广域电网相量测量方法,具有简单易懂、计算速度更快、精度更高等特点,采用加窗插值方法计算信号的频率,从而降低栅栏效应和频谱泄露对频率计算的影响,选用FFT算法使得计算更加快速,易于微机实现。
附图说明
图1是基于频域插值的广域电网相量测量方法的流程图。
具体实施方式
本发明提出了一种基于频域插值的广域电网相量测量方法。以下结合附图作详细说明:
本发明实现基于频域插值的广域电网相量测量方法的流程图如图1所示,其具体步骤包括:
第一步,设置固定采样频率fs=4000Hz,对输入信号进行离散采样,得到采样序列x(m);
第二步,用长度M=1024的三阶最大旁瓣衰减窗对离散序列进行加权,得到xw(m)=x(m)w(m);
第三步,对加权信号xw(m)进行快速傅里叶变换,得到频谱函数|Xw(m)|;
第四步,查找频谱函数|Xw(m)|,得到幅值最大的谱线的位置l=14及其左右两根谱线位置,即13和15,相应的频谱幅值|Xw(l)|=81.78,|Xw(l-1)|=61.6和|Xw(l+1)|=48.08,由插值表达式计算得到η=-0.1657;
第五步,由频率分量偏移值计算得到δ=-0.1488;
第六步,由频率,幅值,相角的计算公式
得到输入信号的频率f0=50.2Hz,幅值A0=220,相角
Claims (4)
1.一种基于频域插值的广域电网相量测量方法,其特征在于:用采样频率fs对电网信号进行采样,得到长度为M的离散采样序列x(m);利用三阶最大旁瓣衰减窗w(m)对离散化的电网信号进行加权得到加权序列xw(m)=x(m)w(m),然后对加权序列做离散傅里叶变换得到频谱函数Xw(m),通过三谱线插值算法计算信号频率偏移值δ,由频率偏移值计算电网信号的频率f0,幅值A0和相角实现电网相量的测量。
2.根据权利要求1所述的方法,其特征在于,所述频率分量偏移值δ的求解,由频谱函数其中函数W(m)的计算表达式为通过搜索|Xw(m)|,其中||表示取模值,得到幅值最大的谱线的位置l及其左右两根谱线位置,即l-1和l+1,由插值表达式:继而获得到频率偏移量δ的计算式:
进而可计算电网信号的频率为
3.根据权利要求1所述的方法,其特征在于,所述电网信号的幅值A0的计算式为:
4.根据权利要求1所述的方法,其特征在于,所述电网信号的相角的的计算式为:
其中Phase{}是求相角运算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710043642.3A CN106802368B (zh) | 2017-01-19 | 2017-01-19 | 一种基于频域插值的广域电网相量测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710043642.3A CN106802368B (zh) | 2017-01-19 | 2017-01-19 | 一种基于频域插值的广域电网相量测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106802368A true CN106802368A (zh) | 2017-06-06 |
CN106802368B CN106802368B (zh) | 2019-10-01 |
Family
ID=58987038
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710043642.3A Active CN106802368B (zh) | 2017-01-19 | 2017-01-19 | 一种基于频域插值的广域电网相量测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106802368B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108390698A (zh) * | 2018-03-16 | 2018-08-10 | 贵州电网有限责任公司 | 一种基于插值fft算法的电力线载波参数测量方法 |
CN109324225A (zh) * | 2018-11-15 | 2019-02-12 | 国网山东省电力公司电力科学研究院 | 一种用于cvt在线监测的高精度电流采样系统 |
CN109884388A (zh) * | 2019-04-02 | 2019-06-14 | 华中科技大学 | 一种基于半周期移相法的电网频率计量、测量装置及方法 |
CN109946512A (zh) * | 2019-04-17 | 2019-06-28 | 贵州电网有限责任公司 | 一种改进频域插值的动态功率分析方法 |
CN110133564A (zh) * | 2019-06-02 | 2019-08-16 | 国网江西省电力有限公司电力科学研究院 | 一种电流传感器相位特性的同步测试方法及装置 |
CN111624400A (zh) * | 2020-04-29 | 2020-09-04 | 中国人民解放军军事科学院国防科技创新研究院 | 正弦信号频率测量方法 |
CN112051444A (zh) * | 2020-08-18 | 2020-12-08 | 许继集团有限公司 | 一种基于抛物线插值迭代的抗谐波频率校准方法及装置 |
CN115420946A (zh) * | 2022-08-18 | 2022-12-02 | 南瑞集团有限公司 | 一种考虑谐波电压的电网电压频率测量方法和装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102680948A (zh) * | 2012-05-15 | 2012-09-19 | 东南大学 | 一种线性调频信号调频率和起始频率估计方法 |
CN103675758A (zh) * | 2013-12-05 | 2014-03-26 | 东南大学 | 一种双曲调频信号周期斜率和起始频率估计方法 |
CN104897961A (zh) * | 2015-06-17 | 2015-09-09 | 中南民族大学 | 基于互乘法窗函数的三谱线插值fft谐波分析方法及系统 |
CN105137175A (zh) * | 2015-08-06 | 2015-12-09 | 湖南大学 | 一种加窗短时傅里叶变换三点插值动态频率测量方法 |
-
2017
- 2017-01-19 CN CN201710043642.3A patent/CN106802368B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102680948A (zh) * | 2012-05-15 | 2012-09-19 | 东南大学 | 一种线性调频信号调频率和起始频率估计方法 |
CN103675758A (zh) * | 2013-12-05 | 2014-03-26 | 东南大学 | 一种双曲调频信号周期斜率和起始频率估计方法 |
CN104897961A (zh) * | 2015-06-17 | 2015-09-09 | 中南民族大学 | 基于互乘法窗函数的三谱线插值fft谐波分析方法及系统 |
CN105137175A (zh) * | 2015-08-06 | 2015-12-09 | 湖南大学 | 一种加窗短时傅里叶变换三点插值动态频率测量方法 |
Non-Patent Citations (3)
Title |
---|
D BELEGA 等: ""Estimation of the Multifrequency Signal Parameters by Interpolated DFT Method with Maximum Sidelobe Decay"", 《 2007 4TH IEEE WORKSHOP ON INTELLIGENT DATA ACQUISITION AND ADVANCED COMPUTING SYSTEMS: TECHNOLOGY AND APPLICATIONS》 * |
温和 等: ""改进加窗插值FFT动态谐波分析算法及应用"", 《电工技术学报》 * |
牛胜锁 等: ""基于三谱线插值FFT的电力谐波分析算法"", 《中国电机工程学报》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108390698A (zh) * | 2018-03-16 | 2018-08-10 | 贵州电网有限责任公司 | 一种基于插值fft算法的电力线载波参数测量方法 |
CN108390698B (zh) * | 2018-03-16 | 2021-08-10 | 贵州电网有限责任公司 | 一种基于插值fft算法的电力线载波参数测量方法 |
CN109324225A (zh) * | 2018-11-15 | 2019-02-12 | 国网山东省电力公司电力科学研究院 | 一种用于cvt在线监测的高精度电流采样系统 |
CN109884388B (zh) * | 2019-04-02 | 2021-05-18 | 华中科技大学 | 一种基于半周期移相法的电网频率计量、测量装置及方法 |
CN109884388A (zh) * | 2019-04-02 | 2019-06-14 | 华中科技大学 | 一种基于半周期移相法的电网频率计量、测量装置及方法 |
CN109946512A (zh) * | 2019-04-17 | 2019-06-28 | 贵州电网有限责任公司 | 一种改进频域插值的动态功率分析方法 |
CN109946512B (zh) * | 2019-04-17 | 2019-12-03 | 贵州电网有限责任公司 | 一种改进频域插值的动态功率分析方法 |
CN110133564B (zh) * | 2019-06-02 | 2021-04-09 | 国网江西省电力有限公司电力科学研究院 | 一种电流传感器相位特性的同步测试方法及装置 |
CN110133564A (zh) * | 2019-06-02 | 2019-08-16 | 国网江西省电力有限公司电力科学研究院 | 一种电流传感器相位特性的同步测试方法及装置 |
CN111624400A (zh) * | 2020-04-29 | 2020-09-04 | 中国人民解放军军事科学院国防科技创新研究院 | 正弦信号频率测量方法 |
CN111624400B (zh) * | 2020-04-29 | 2021-10-19 | 中国人民解放军军事科学院国防科技创新研究院 | 正弦信号频率测量方法 |
CN112051444A (zh) * | 2020-08-18 | 2020-12-08 | 许继集团有限公司 | 一种基于抛物线插值迭代的抗谐波频率校准方法及装置 |
CN112051444B (zh) * | 2020-08-18 | 2023-02-10 | 许继集团有限公司 | 一种基于抛物线插值迭代的抗谐波频率校准方法及装置 |
CN115420946A (zh) * | 2022-08-18 | 2022-12-02 | 南瑞集团有限公司 | 一种考虑谐波电压的电网电压频率测量方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN106802368B (zh) | 2019-10-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106802368A (zh) | 一种基于频域插值的广域电网相量测量方法 | |
CN103308766A (zh) | 一种基于凯撒自卷积窗双谱线插值fft谐波分析方法及其装置 | |
CN103576002B (zh) | 一种容性绝缘设备介质损耗角的计算方法 | |
CN103308804A (zh) | 基于快速k-s变换电能质量扰动信号时频参数提取方法 | |
CN101701984B (zh) | 基于三项系数Nuttall窗插值FFT的基波与谐波检测方法 | |
CN105137180B (zh) | 基于六项余弦窗四谱线插值的高精度谐波分析方法 | |
CN101216512A (zh) | 一种非正弦周期信号实时高精度检测方法 | |
CN101806832A (zh) | 一种低频率信号的频率测量方法 | |
CN105044456A (zh) | 一种基于正交子带的电网瞬时频率测量与跟踪方法 | |
CN103728663B (zh) | 一种时频分析方法 | |
CN106324340B (zh) | 一种同步相量和频率测量动态性能的方法 | |
CN111046327A (zh) | 适用于低频振荡与次同步振荡辨识的Prony分析方法 | |
Stojanović et al. | The algorithm for directional element without dead tripping zone based on digital phase comparator | |
CN107490397A (zh) | 高精度自适应滤波fbg光谱快速寻峰方法 | |
CN103091545A (zh) | 一种频率无关的正弦信号相量半波计算方法 | |
CN104849545A (zh) | 一种电力系统频率测量方法及测量装置 | |
CN102967779A (zh) | 一种输电线路分布参数的辨识方法 | |
CN105445614B (zh) | 一种基于小波分析的双端行波故障定位的方法和系统 | |
CN102519597A (zh) | 一种傅里叶变换光谱仪相位校正切趾方法 | |
CN102508022A (zh) | 采用最优乘子牛顿算法检测电网频率的方法 | |
CN106446405B (zh) | 一种集成电路器件神经网络建模样本选择方法及装置 | |
CN103543331B (zh) | 一种计算电信号谐波和间谐波的方法 | |
CN101718816B (zh) | 基于四项系数Nuttall窗插值FFT的基波与谐波检测方法 | |
CN106053937A (zh) | 一种基于fft+ft的基波测频方法 | |
Sahoo et al. | An enhanced travelling wave-based fault detection and location estimation technique for series compensated transmission network |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |