CN109030938B - 一种基于正弦滤波的抗谐波测频方法和装置 - Google Patents
一种基于正弦滤波的抗谐波测频方法和装置 Download PDFInfo
- Publication number
- CN109030938B CN109030938B CN201710429103.3A CN201710429103A CN109030938B CN 109030938 B CN109030938 B CN 109030938B CN 201710429103 A CN201710429103 A CN 201710429103A CN 109030938 B CN109030938 B CN 109030938B
- Authority
- CN
- China
- Prior art keywords
- power system
- frequency
- frequency measurement
- signal
- low
- 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.)
- Active
Links
- 238000001914 filtration Methods 0.000 title claims abstract description 54
- 238000000691 measurement method Methods 0.000 title claims abstract description 20
- 238000005259 measurement Methods 0.000 claims abstract description 81
- 238000005070 sampling Methods 0.000 claims description 20
- 238000000034 method Methods 0.000 abstract description 20
- 238000004364 calculation method Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/02—Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage
- G01R23/04—Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage adapted for measuring in circuits having distributed constants
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
- G01R23/165—Spectrum analysis; Fourier analysis using filters
- G01R23/167—Spectrum analysis; Fourier analysis using filters with digital filters
Abstract
本发明提供一种基于正弦滤波的抗谐波测频方法和装置,其中方法包括如下步骤:对电力系统的频率信号进行低通滤波,得到电力系统的低频信号;对电力系统的低频信号进行正弦滤波,得到电力系统的测频信号;对电力系统的测频信号进行测频处理,得到电力系统的实时频率。本发明提供的技术方案,在对电力系统进行测频处理时,首先消除其中的高频信号,然后消除其中的低频谐波的干扰,所以测频结果更加准确。
Description
技术领域
本发明属于电力系统检测技术领域,具体涉及一种基于正弦滤波的抗谐波测频方法和装置。
背景技术
在电力系统中,频率是电能质量的重要指标之一,是反映电力系统运行状态的重要参数。一般情况下,电力系统的频率反映了电力系统中有功功率供需平衡的基本状态,它将随负荷波动在小范围内缓慢的变化。当电力系统处于稳定运行状态下时,发电机的输出功率与电力系统负荷和损耗的功率维持平衡,电力系统的频率为标称值。如果电力系统中有大容量的负荷或大功率的发电机投切,或者控制设备不完善,就会导致电力系统频率的偏移,对电力系统的稳定运行和用户设备的正常工作带来影响,所以对电力系统频率的检测非常重要。
目前测量电力系统频率的方法很多,如过零点测频、相量测频、最小二乘法测频等,当电力系统的频率偏移时,如果电力系统的波形比较纯净,这些方法都均能获得较好的测频效果。但是检测到的电力系统的真实波形,其中都会含有谐波,如此以来这些方法的测频结果都会出现一定的偏差,甚至出现严重的错误。
2005年4月的《电力系统及其自动化学报》记载了一篇名称为“电力系统频率测量的一种新方法”的论文,该论文公开了一种在被测信号进入AD采样之前,采用性能良好的带通或低通滤波器消除被测信号中的高频噪声和谐波。使用带通滤波器虽然能够得到较好的测频信号,但是高精度的带通滤波器设计比较困难,需要的成本比较高;采用低通滤波器虽然能够很好的滤除掉被测信号中的高频噪声,但是不能滤除其中的低频谐波,尤其是频率在30Hz-45Hz之间的谐波。所以经过低通滤波器滤波处理后的测频信号中依然存在干扰信号,导致测频结果不准确。
发明内容
本发明的目的在于提供一种基于正弦滤波的抗谐波测频方法和装置,用于解决电力系统存在谐波时测频结果不准确的问题。
为实现上述目的,本发明提供的技术方案是:
一种基于正弦滤波的抗谐波测频方法,包括如下步骤:
(1)对电力系统的频率信号进行低通滤波,得到电力系统的低频信号;
(2)对电力系统的低频信号进行正弦滤波,得到电力系统的测频信号;
(3)对电力系统的测频信号进行测频处理,得到电力系统的实时频率。
进一步的,在对电力系统的低频信号进行正弦滤波时,滤波公式为
其中x(n)为电力系统的低频信号,n为采样序号,N为装置采样率下的一周波点数,y(n)为电力系统的测频信号。
进一步的,所述步骤(3)中采用过零测频法对电力系统的测频信号进行测频处理,得到电力系统的实时频率。
进一步的,在采用过零侧频法对电力系统的测频信号进行测频处理时,通过线性插值算法得出电力系统测频信号的过零点。
一种基于正弦滤波的抗谐波测频装置,包括:
对电力系统的频率信号进行低通滤波,得到电力系统的低频信号的模块;
对电力系统的低频信号进行正弦滤波,得到电力系统的测频信号的模块;
对电力系统的测频信号进行测频处理,得到电力系统的实时频率的模块。
进一步的,在对电力系统的低频信号进行正弦滤波时,滤波公式为
其中x(n)为电力系统的低频信号,n为采样序号,N为装置采样率下的一周波点数,y(n)为电力系统的测频信号。
进一步的,采用过零测频法对电力系统的测频信号进行测频处理,得到电力系统的实时频率。
进一步的,在采用过零侧频法对电力系统的测频信号进行测频处理时,通过线性插值算法得出电力系统测频信号的过零点。
本发明的有益效果是:本发明提供的技术方案,在对电力系统进行测频处理时,首先对电力系统的频率信号进行低通滤波,消除其中的高频信号,然后再进行正弦滤波,消除其中的低频谐波的干扰,所以测频结果更加准确。
附图说明
图1为方法实施例中全周傅氏、傅氏余弦部分和傅氏正弦部分的幅频特性图;
图2为方法实施例中过零测频法的示意图。
具体实施方式
本发明的目的在于提供一种基于正弦滤波的抗谐波测频方法和装置,用于解决电力系统存在谐波时测频结果不准确的问题。
为实现上述目的,本发明提供的技术方案是:
一种基于正弦滤波的抗谐波测频方法,包括如下步骤:
(1)对电力系统的频率信号进行低通滤波,得到电力系统的低频信号;
(2)对电力系统的低频信号进行正弦滤波,得到电力系统的测频信号;
(3)对电力系统的测频信号进行测频处理,得到电力系统的实时频率。
下面结合附图对本发明的实施方式作进一步说明。
方法实施例:
本实施例提供一种基于正弦滤波的抗谐波测频方法,用于在对电力系统进行测频时,消除电力系统频率信号中的谐波,具体步骤如下:
(1)采用低通滤波器对电力系统的频率信号进行低通滤波,消除电力系统频率信号中的高频干扰信号,得到电力系统的低频信号;
我国电力系统的额定频率为50Hz,为了适应频率为50Hz的电力系统,在本实施例中按照滤除90Hz以上谐波和阶数尽量低的原则,运用切比雪夫滤波器设计方法,调节截止频率和衰减幅度,得到滤波系数,对电力系统的频率信号进行低通滤波,得到电力系统的低频信号;例如采样率为1200Hz,低通滤波器的通带截至频率为3Hz,阻带截至频率为88Hz,通带最大波纹衰减为1.8967dB,阻带最大衰减55dB,阶数为24阶,具体系数见表1。
表1
(2)对电力系统的低频信号进行正弦滤波,消除其中的低频谐波,得到电力系统的测频信号;
在本实施例中,利用傅氏正弦部分具有比全周傅氏更好的低通滤波特性,在电力系统频率偏移过低时,削弱谐波对测频精度的影响;图1中为全周傅氏、傅氏余弦部分和傅氏正弦部分的幅频特性,对他们的幅频特性对比分析,可知傅氏正弦部分在50-100Hz有更好的低通滤波特性;
在低通滤波的基础上,对信号使用傅氏正弦滤波,滤波的公式如下:
其中x(n)为电力系统的低频信号,n为采样序号,N为装置采样率下的一周波点数,即装置采样率下,一个周期的波形内的点数,y(n)为电力系统的测频信号;
(3)对电力系统的测频信号进行测频处理,得到电力系统的实时频率;
根据经过傅氏正弦滤波的电力系统测频信号,得到电力系统实时频率的方法很多,在本实施例中采用计算量较小的过零侧频法进行测频;采样装置在对电力系统的频率信号进行采样时,不一定能采集到电力系统频率信号的过零点,所以需要对电力系统的测频信号进行插值运算,寻找过零点;在本实施例中,采用线性插值寻找电力系统测频信号中的过零点,如图2所示,计算公式为
式中,x(tn)和x(tn+1)满足关系式x(tn)·x(tn+1)<0,tn和tn+1为相邻两个采样点的采样时刻,t0为过零点时刻;
在采用过零侧频法进行测频时,根据3个连续的过零点即可确定当前测频信号的频率,计算公式为:
式中,t0-1为一周波的第一个过零点,t0-3为一周波的第三个过零点。
作为其他实施方式,可将低通滤波和正弦滤波合并,得到一组新的系数,合并方式如下:
为了验证上述基于正弦滤波的抗谐波测频方法的效果,本实施例中还利用上述方法得到表2中数据,记录不同频率下的频率测量误差,信号模型如下:
xsig=A·sin(w·t+phase1)+0.2·A·sin(h·w·t+phase2)
式中,A为信号幅值,w为信号频率,h谐波次数,phase1为信号相位,phase2为谐波相位。
表2中的数据,phase1取0度,phase2取60度,A取单位值。从表中可以看出本发明的测频精度在30~60Hz的宽范围内,可以达到0.01Hz以下的绝对精度。
表2
装置实施例:
本实施例提供一种基于正弦滤波的抗谐波测频装置,包括:
对电力系统的频率信号进行低通滤波,得到电力系统的低频信号的模块;
对电力系统的低频信号进行正弦滤波,得到电力系统的测频信号的模块;
对电力系统的测频信号进行测频处理,得到电力系统的实时频率的模块。
本实施例所提供的一种基于正弦滤波的抗谐波测频装置,其中各模块并不是硬件模块,而是按照上述方法实施例所提供的方法进行编程得到的软件模块,运行在相应的处理器或者控制器中,能够存储在可移动或者固定的存储装置中。
Claims (6)
2.根据权利要求1所述的一种基于正弦滤波的抗谐波测频方法,其特征在于,所述步骤(3)中采用过零测频法对电力系统的测频信号进行测频处理,得到电力系统的实时频率。
3.根据权利要求2所述的一种基于正弦滤波的抗谐波测频方法,其特征在于,在采用过零测 频法对电力系统的测频信号进行测频处理时,通过线性插值算法得出电力系统测频信号的过零点。
5.根据权利要求4所述的一种基于正弦滤波的抗谐波测频装置,其特征在于,采用过零测频法对电力系统的测频信号进行测频处理,得到电力系统的实时频率。
6.根据权利要求5所述的一种基于正弦滤波的抗谐波测频装置,其特征在于,在采用过零测 频法对电力系统的测频信号进行测频处理时,通过线性插值算法得出电力系统测频信号的过零点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710429103.3A CN109030938B (zh) | 2017-06-08 | 2017-06-08 | 一种基于正弦滤波的抗谐波测频方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710429103.3A CN109030938B (zh) | 2017-06-08 | 2017-06-08 | 一种基于正弦滤波的抗谐波测频方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109030938A CN109030938A (zh) | 2018-12-18 |
CN109030938B true CN109030938B (zh) | 2021-05-11 |
Family
ID=64628480
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710429103.3A Active CN109030938B (zh) | 2017-06-08 | 2017-06-08 | 一种基于正弦滤波的抗谐波测频方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109030938B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH03205572A (ja) * | 1990-01-08 | 1991-09-09 | Matsushita Electron Corp | 半導体デバイスの測定方法および測定装置 |
US7424065B2 (en) * | 2002-10-22 | 2008-09-09 | Walker Harold R | Apparatus and method for ultra narrow band wireless communications |
CN102288820A (zh) * | 2011-08-05 | 2011-12-21 | 上海理工大学 | 基于锁相环和神经网络相结合的谐波检测方法 |
CN102457250A (zh) * | 2010-10-20 | 2012-05-16 | Tcl集团股份有限公司 | 一种采集数据滤波处理方法及装置 |
CN103023292A (zh) * | 2012-11-27 | 2013-04-03 | 深圳市晶福源电子技术有限公司 | 一种特定谐波消除系统及其方法 |
CN104991104A (zh) * | 2015-07-09 | 2015-10-21 | 广东电网有限责任公司电力科学研究院 | 电力信号的幅值检测方法和系统 |
CN105510719A (zh) * | 2015-11-30 | 2016-04-20 | 云南电网有限责任公司玉溪供电局 | 三相电网谐波阻抗测量方法 |
Family Cites Families (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008047944A (ja) * | 2004-12-01 | 2008-02-28 | Advantest Corp | Da変換器の試験方法、da変換器の試験装置およびda変換器 |
US8150643B1 (en) * | 2004-12-21 | 2012-04-03 | Battelle Energy Alliance, Llc | Method of detecting system function by measuring frequency response |
CN100371724C (zh) * | 2005-02-26 | 2008-02-27 | 重庆大学 | 消弧线圈接地系统接地故障线路选择方法 |
CN101115037B (zh) * | 2007-08-13 | 2011-12-14 | 北京天碁科技有限公司 | 一种数字接收机及其本振频率初始误差校准的方法 |
CN101441233A (zh) * | 2008-12-02 | 2009-05-27 | 湖南海兴电器有限责任公司 | 基于Kaiser窗双谱线插值FFT的基波与谐波检测方法 |
CN101458279B (zh) * | 2008-12-09 | 2010-11-03 | 中国科学院国家授时中心 | 一种能提高精度的差拍法频率测量方法 |
CN101505060B (zh) * | 2009-03-18 | 2010-11-10 | 天津大学 | 电网中利用窄带滤波实现多采样率信号的频率转换方法 |
CN101915874B (zh) * | 2010-07-20 | 2013-01-09 | 北海市深蓝科技发展有限责任公司 | 一种基于傅立叶变换的谐波检测方法 |
CN102116796B (zh) * | 2010-12-10 | 2013-03-06 | 天津凯发电气股份有限公司 | 基于eTPU的交流信号同步采样和频率测量装置及其工作方法 |
CN102253281A (zh) * | 2011-04-26 | 2011-11-23 | 江苏斯菲尔电气股份有限公司 | 一种基于软件抗谐波干扰并自动换相测量电网频率的方法 |
CN102495282A (zh) * | 2011-12-27 | 2012-06-13 | 国电南瑞科技股份有限公司 | 一种交流电压频率的测量方法 |
CN103941088A (zh) * | 2014-04-10 | 2014-07-23 | 山东大学 | 一种基于三相信号的电力系统频率快速测量方法 |
CN103913631A (zh) * | 2014-04-30 | 2014-07-09 | 乐山一拉得电网自动化有限公司 | 用于数显表的电网频率测量方法和测量系统 |
CN104198811B (zh) * | 2014-08-18 | 2017-02-01 | 广东电网公司电力科学研究院 | 低频信号频率测量方法和装置 |
CN105004920B (zh) * | 2015-07-10 | 2017-11-17 | 国网天津市电力公司 | 傅里叶修正系数频率测量方法 |
CN105044456B (zh) * | 2015-07-21 | 2018-04-13 | 电子科技大学 | 一种基于正交子带的电网瞬时频率测量与跟踪方法 |
CN105629060B (zh) * | 2015-12-24 | 2018-05-29 | 电子科技大学 | 基于最优基带滤波的电网频率测量方法和装置 |
CN105629064B (zh) * | 2015-12-30 | 2018-10-02 | 浙江埃菲生能源科技有限公司 | 一种正弦波低通滤波提取方法 |
-
2017
- 2017-06-08 CN CN201710429103.3A patent/CN109030938B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH03205572A (ja) * | 1990-01-08 | 1991-09-09 | Matsushita Electron Corp | 半導体デバイスの測定方法および測定装置 |
US7424065B2 (en) * | 2002-10-22 | 2008-09-09 | Walker Harold R | Apparatus and method for ultra narrow band wireless communications |
CN102457250A (zh) * | 2010-10-20 | 2012-05-16 | Tcl集团股份有限公司 | 一种采集数据滤波处理方法及装置 |
CN102288820A (zh) * | 2011-08-05 | 2011-12-21 | 上海理工大学 | 基于锁相环和神经网络相结合的谐波检测方法 |
CN103023292A (zh) * | 2012-11-27 | 2013-04-03 | 深圳市晶福源电子技术有限公司 | 一种特定谐波消除系统及其方法 |
CN104991104A (zh) * | 2015-07-09 | 2015-10-21 | 广东电网有限责任公司电力科学研究院 | 电力信号的幅值检测方法和系统 |
CN105510719A (zh) * | 2015-11-30 | 2016-04-20 | 云南电网有限责任公司玉溪供电局 | 三相电网谐波阻抗测量方法 |
Non-Patent Citations (5)
Title |
---|
Comparison of Nonuniform Transmission Lines With Gaussian and Exponential Impedance Profiles for Z-Pinch;Rui Zhang 等;《IEEE TRANSACTIONS ON PLASMA SCIENCE》;20121231;第40卷(第12期);第3395-3398页 * |
不同傅氏算法对衰减分量的抑制性能分析;张涛 等;《继电器》;20040630;第32卷(第11期);第25-28页 * |
傅氏算法的滤波特性分析;金德福 等;《继电器》;20051130;第33卷(第21期);第45-50页 * |
基于调整滤波器初相位的新型半波傅氏算法;安林 等;《电力自动化设备》;20041031;第24卷(第10期);第19-23页 * |
电力系统数字保护算法的仿真分析;冯巧玲 等;《郑州轻工业学院学报(自然科学版)》;20070831;第22卷(第4期);第49-51页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109030938A (zh) | 2018-12-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109142863B (zh) | 一种电力系统测频方法及系统 | |
CN106679741B (zh) | 基于涡街流量计抗干扰信号的处理方法及系统 | |
CN101915874B (zh) | 一种基于傅立叶变换的谐波检测方法 | |
CN107179440B (zh) | 一种在线自适应频率变化的次同步振荡间谐波提取方法 | |
CN103257271A (zh) | 一种基于stm32f107vct6的微电网谐波与间谐波检测装置及检测方法 | |
CN203287435U (zh) | 一种基于stm32f107vct6的微电网谐波与间谐波检测装置 | |
de Carvalho et al. | A PLL-based multirate structure for time-varying power systems harmonic/interharmonic estimation | |
CN107576851B (zh) | 基于旋转电流相量的系统谐波阻抗测量方法 | |
CN109669072B (zh) | 一种配电网的自适应同步相量量测方法 | |
CN104316768A (zh) | 一种三相不平衡扰动源定位的负序阻抗参数估算方法 | |
CN109061345B (zh) | 适用于电力系统的有效值测量方法与系统 | |
CN104101785A (zh) | 一种四端法高值电容阻抗测量装置及其测量方法 | |
CN107085134A (zh) | 一种电压闪变的检测方法 | |
Kaiser et al. | Estimation of power systems amplitudes, frequencies, and phase characteristics using energy operators | |
CN111505375A (zh) | 一种频率测量方法及装置 | |
CN111948454B (zh) | 一种抗衰减直流分量干扰的同步相量测量方法及系统 | |
CN109030938B (zh) | 一种基于正弦滤波的抗谐波测频方法和装置 | |
CN103592514A (zh) | 一种新型谐波高精度检测方法 | |
CN108334822B (zh) | 基于电动汽车充电非线性负荷特征的卡尔曼和修正小波变换滤波方法 | |
CN113358922B (zh) | 非稳态电网信号的基波、谐波、间谐波电能功率测量方法 | |
Serov et al. | Method to reduce the measurement error of the spectrum by the demodulation technique | |
CN111624408B (zh) | 基于共生多元泛函计算的电力线路对地电容实时测量方法 | |
Zarraga | On-line extraction of modal characteristics from power system measurements based on Hilbert-Huang analysis | |
CN112557750B (zh) | 一种计量装置互感器二次回路信号频率测量方法及系统 | |
Zhou et al. | Study on Self-Monitoring Scheme of the Accuracy Deviation for the Merging Unit |
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 |