CN113447713A - 一种基于傅式快速高精度的电力***频率测量方法及装置 - Google Patents
一种基于傅式快速高精度的电力***频率测量方法及装置 Download PDFInfo
- Publication number
- CN113447713A CN113447713A CN202110713847.4A CN202110713847A CN113447713A CN 113447713 A CN113447713 A CN 113447713A CN 202110713847 A CN202110713847 A CN 202110713847A CN 113447713 A CN113447713 A CN 113447713A
- Authority
- CN
- China
- Prior art keywords
- frequency
- sampling
- input signal
- fourier
- calculating
- 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
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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
本发明公开了一种基于傅式快速高精度的电力***频率测量方法及装置,所述方法包括:设定输入信号u(t),所述输入信号u(t)是频率为fx的正弦波电压信号;以输入信号的中心频率为中心设置M个采样率,周波采样N次;使用M个采样率同时对输入信号u(t)进行采样,采样数据进行缓存;对输入信号u(t)进行采样过程中,每间隔一个周波计算一次频率;使用中心频率下的采样数据,计算得出中间频率f1;从M个采样率中选取与中间频率f1最接近的采样率对应的采样数据组;利用傅式测频算法计算最终的频率f;本发明的优点在于:实现了计算精度和计算速度之间较好的统一,满足一次调频的要求。
Description
技术领域
本发明涉及电力***频率测量技术领域,更具体涉及一种基于傅式快速高精度的电力***频率测量方法及装置。
背景技术
电力***频率是重要的电能质量指标之一。频率是反映电力***运行特性的重要参数,通过及时准确地测量频率,可以预测电力***是否将失去稳定,进行一调频、切机、切负荷控制来保证***的安全运行。电力***测频需要从扰动信号中提取频率,测频精度需要达到0.003hz,测频时间需要控制在100ms以内,以便实现一次调频功能。但是现有技术电力***频率测量方法计算精度和计算速度之间不能较好的统一,不能同时保证频率测量的计算精度和计算速度,导致难以满足一次调频的要求。
2005年1月中国电机工程学报第26卷第2期,公开了一种基于傅氏算法的高精度测频方法,该文献仔细研究了正弦信号经傅氏算法变换后的结果,发现随着数据窗的推移,傅氏算法得到的相量实部和虚部满足一个恒等式,由此得到一种新的测频方法,为了提高谐波情况下测频的精度,以前一次频率测量值为基础进行迭代,在对采样数据插值的基础上,使得迭代很快的收敛,结果表明,算法在各种情况下都具有很高的计算精度。但是其虽然保证了计算精度,对于计算速度没有进行深入研究,不能同时保证频率测量的计算精度和计算速度,计算精度和计算速度之间不能较好的统一,难以满足一次调频的要求。
发明内容
本发明所要解决的技术问题在于现有技术测频方法存在计算精度和计算速度之间不能较好的统一,难以满足一次调频的要求的问题。
本发明通过以下技术手段实现解决上述技术问题的:一种基于傅式快速高精度的电力***频率测量方法,所述方法包括:
步骤A:设定输入信号u(t),所述输入信号u(t)是频率为fx的正弦波电压信号,输入信号u(t)包含电力***中的谐波分量;
步骤B:以输入信号的中心频率为中心设置M个采样率,周波采样N次,其中M和N均为大于1的正整数;
步骤C:使用M个采样率同时对输入信号u(t)进行采样,采样数据进行缓存;
步骤D:对输入信号u(t)进行采样过程中,每间隔一个周波计算一次频率;
步骤E:当频率变化时,使用中心频率下的采样数据,按照傅式测频算法计算得出中间频率f1;
步骤F:从M个采样率中选取与中间频率f1最接近的采样率对应的采样数据组;
步骤G:根据M个采样率中与中间频率f1最接近的采样率对应的采样数据组,利用傅式测频算法计算最终的频率f。
本发明使用多个采样率,对同一个输入信号进行同步采样,选取中心频率下采样数据进行傅式测频计算,得出中间频率f1,再根据中间频率f1,选取合适采样率数据进行傅式计算得到最终结果频率f,通过先计算中间频率再根据中间频率计算最终的频率的方式提高整体计算精度,计算过程简单,算法不需要迭代寻优,计算速度大大提高,因此本发明实现了计算精度和计算速度之间较好的统一,满足一次调频的要求。
进一步地,所述步骤A中谐波含量包括3×f0和5×f0,f0=50hz,f0为输入信号u(t)的中心频率。
更进一步地,所述步骤B包括:
按照额定输入频率fs={44,45,46,47,48,49,50,51,52,53,54,55}设置12个采样率,周波采样256次。
更进一步地,所述步骤E包括:
使用50hz采用率对输入信号u(t)进行采样过程中,当频率变化时,根据两个相邻数据窗间的相角差利用公式fb=ΔQ/2πΔt求出变化后的频率,其中,ΔQ表示两个相邻数据窗间的相角差,Δt表示两个相邻数据窗间的时间间隔;
每个采样点计算一次角度差,依次利用上述求取变化后的频率的公式计算得到一组临时频率fbuf,对该组临时频率fbuf,进行周期抽取后平均,得到中间频率f1。
更进一步地,所述数据窗的长度为一个采样周期。
本发明还提供一种基于傅式快速高精度的电力***频率测量装置,所述装置包括:
输入信号设定模块,用于设定输入信号u(t),所述输入信号u(t)是频率为fx的正弦波电压信号,输入信号u(t)包含电力***中的谐波分量;
采样率设置模块,用于以输入信号的中心频率为中心设置M个采样率,周波采样N次,其中M和N均为大于1的正整数;
采样模块,用于使用M个采样率同时对输入信号u(t)进行采样,采样数据进行缓存;
频率计算模块,用于对输入信号u(t)进行采样过程中,每间隔一个周波计算一次频率;
中间频率获取模块,用于当频率变化时,使用中心频率下的采样数据,按照傅式测频算法计算得出中间频率f1;
采样数据选取模块,用于从M个采样率中选取与中间频率f1最接近的采样率对应的采样数据组;
测频模块,用于根据M个采样率中与中间频率f1最接近的采样率对应的采样数据组,利用傅式测频算法计算最终的频率f。
进一步地,所述输入信号设定模块中谐波含量包括3×f0和5×f0,f0=50hz,f0为输入信号u(t)的中心频率。
更进一步地,所述采样率设置模块还用于:
按照额定输入频率fs={44,45,46,47,48,49,50,51,52,53,54,55}设置12个采样率,周波采样256次。
更进一步地,所述中间频率获取模块还用于:
使用50hz采用率对输入信号u(t)进行采样过程中,当频率变化时,根据两个相邻数据窗间的相角差利用公式fb=ΔQ/2πΔt求出变化后的频率,其中,ΔQ表示两个相邻数据窗间的相角差,Δt表示两个相邻数据窗间的时间间隔;
每个采样点计算一次角度差,依次利用上述求取变化后的频率的公式计算得到一组临时频率fbuf,对该组临时频率fbuf,进行周期抽取后平均,得到中间频率f1。
更进一步地,所述数据窗的长度为一个采样周期。
本发明的优点在于:本发明使用多个采样率,对同一个输入信号进行同步采样,选取中心频率下采样数据进行傅式测频计算,得出中间频率f1,再根据中间频率f1,选取合适采样率数据进行傅式计算得到最终结果频率f,通过先计算中间频率再根据中间频率计算最终的频率的方式提高整体计算精度,计算过程简单,算法不需要迭代寻优,计算速度大大提高,因此本发明实现了计算精度和计算速度之间较好的统一,满足一次调频的要求。
附图说明
图1为本发明实施例所公开的一种基于傅式快速高精度的电力***频率测量方法的算法流程图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例1
如图1所示,一种基于傅式快速高精度的电力***频率测量方法,所述方法包括:
步骤A:设定输入信号u(t),所述输入信号u(t)是频率为fx的正弦波电压信号,输入信号u(t)包含电力***中常见的谐波分量,谐波含量为3×f0,5×f0等,f0=50hz,f0为输入信号u(t)的中心频率。
步骤B:以输入信号的中心频率为中心设置M个采样率,周波采样N次,其中M和N均为大于1的正整数;具体过程为:按照额定输入频率fs={44,45,46,47,48,49,50,51,52,53,54,55}设置12个采样率,周波采样256次。
步骤C:使用上述12个采样率同时对输入信号u(t)进行采样,采样数据进行缓存;
步骤D:对输入信号u(t)进行采样过程中,每间隔一个周波计算一次频率;具体过程为:每隔20ms计算一次频率,每次使用3个周波数据,整体频率响应时间在60-80ms之间。计算频率的方法可以采用现有技术的多种频率计算方法,其中一种方法为利用相角差进行计算,具体公式参考步骤E。
步骤E:当频率变化时,使用中心频率下的采样数据,按照傅式测频算法计算得出中间频率f1;具体过程为:使用50hz采用率对输入信号u(t)进行采样过程中,当频率变化(步骤D已经计算采样过程的频率,所以这里能够得到频率变化)时,根据两个相邻数据窗间的相角差利用公式fb=ΔQ/2πΔt求出变化后的频率,其中,ΔQ表示两个相邻数据窗间的相角差,Δt表示两个相邻数据窗间的时间间隔;因此只要求出一段时间Δt内的角度变化量,就可以得到***的频率值。
每个采样点(本发明周波采样256次,因此每个周期内采样点256个)计算一次角度差,依次利用上述求取变化后的频率的公式计算得到一组临时频率fbuf,对该组临时频率fbuf,进行周期抽取后平均,得到中间频率f1。
步骤F:从12个采样率中选取与中间频率f1最接近的采样率对应的采样数据组;
步骤G:根据12个采样率中与中间频率f1最接近的采样率对应的采样数据组,利用傅式测频算法也即上述求取变化后的频率的公式计算最终的频率f。所述数据窗的长度为一个采样周期。
通过以上技术方案,本发明使用多个采样率,对同一个输入信号进行同步采样,选取中心频率下采样数据进行傅式测频计算,得出中间频率f1,再根据中间频率f1,选取合适采样率数据进行傅式计算得到最终结果频率f,通过先计算中间频率再根据中间频率计算最终的频率的方式提高整体计算精度,计算过程简单,算法不需要迭代寻优,计算速度大大提高,因此本发明实现了计算精度和计算速度之间较好的统一,满足一次调频的要求。
实施例2
基于本发明实施例1,本发明实施例2还提供一种基于傅式快速高精度的电力***频率测量装置,所述装置包括:
输入信号设定模块,用于设定输入信号u(t),所述输入信号u(t)是频率为fx的正弦波电压信号,输入信号u(t)包含电力***中的谐波分量;
采样率设置模块,用于以输入信号的中心频率为中心设置M个采样率,周波采样N次,其中M和N均为大于1的正整数;
采样模块,用于使用M个采样率同时对输入信号u(t)进行采样,采样数据进行缓存;
频率计算模块,用于对输入信号u(t)进行采样过程中,每间隔一个周波计算一次频率;
中间频率获取模块,用于当频率变化时,使用中心频率下的采样数据,按照傅式测频算法计算得出中间频率f1;
采样数据选取模块,用于从M个采样率中选取与中间频率f1最接近的采样率对应的采样数据组;
测频模块,用于根据M个采样率中与中间频率f1最接近的采样率对应的采样数据组,利用傅式测频算法计算最终的频率f。
具体的,所述输入信号设定模块中谐波含量包括3×f0和5×f0,f0=50hz,f0为输入信号u(t)的中心频率。
更具体的,所述采样率设置模块还用于:
按照额定输入频率fs={44,45,46,47,48,49,50,51,52,53,54,55}设置12个采样率,周波采样256次。
更具体的,所述中间频率获取模块还用于:
使用50hz采用率对输入信号u(t)进行采样过程中,当频率变化时,根据两个相邻数据窗间的相角差利用公式fb=ΔQ/2πΔt求出变化后的频率,其中,ΔQ表示两个相邻数据窗间的相角差,Δt表示两个相邻数据窗间的时间间隔;
每个采样点计算一次角度差,依次利用上述求取变化后的频率的公式计算得到一组临时频率fbuf,对该组临时频率fbuf,进行周期抽取后平均,得到中间频率f1。
更具体的,所述数据窗的长度为一个采样周期。
以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (10)
1.一种基于傅式快速高精度的电力***频率测量方法,其特征在于,所述方法包括:
步骤A:设定输入信号u(t),所述输入信号u(t)是频率为fx的正弦波电压信号,输入信号u(t)包含电力***中的谐波分量;
步骤B:以输入信号的中心频率为中心设置M个采样率,周波采样N次,其中M和N均为大于1的正整数;
步骤C:使用M个采样率同时对输入信号u(t)进行采样,采样数据进行缓存;
步骤D:对输入信号u(t)进行采样过程中,每间隔一个周波计算一次频率;
步骤E:当频率变化时,使用中心频率下的采样数据,按照傅式测频算法计算得出中间频率f1;
步骤F:从M个采样率中选取与中间频率f1最接近的采样率对应的采样数据组;
步骤G:根据M个采样率中与中间频率f1最接近的采样率对应的采样数据组,利用傅式测频算法计算最终的频率f。
2.根据权利要求1所述的一种基于傅式快速高精度的电力***频率测量方法,其特征在于,所述步骤A中谐波含量包括3×f0和5×f0,f0=50hz,f0为输入信号u(t)的中心频率。
3.根据权利要求2所述的一种基于傅式快速高精度的电力***频率测量方法,其特征在于,所述步骤B包括:
按照额定输入频率fs={44,45,46,47,48,49,50,51,52,53,54,55}设置12个采样率,周波采样256次。
4.根据权利要求3所述的一种基于傅式快速高精度的电力***频率测量方法,其特征在于,所述步骤E包括:
使用50hz采用率对输入信号u(t)进行采样过程中,当频率变化时,根据两个相邻数据窗间的相角差利用公式fb=ΔQ/2πΔt求出变化后的频率,其中,ΔQ表示两个相邻数据窗间的相角差,Δt表示两个相邻数据窗间的时间间隔;
每个采样点计算一次角度差,依次利用上述求取变化后的频率的公式计算得到一组临时频率fbuf,对该组临时频率fbuf,进行周期抽取后平均,得到中间频率f1。
5.根据权利要求4所述的一种基于傅式快速高精度的电力***频率测量方法,其特征在于,所述数据窗的长度为一个采样周期。
6.一种基于傅式快速高精度的电力***频率测量装置,其特征在于,所述装置包括:
输入信号设定模块,用于设定输入信号u(t),所述输入信号u(t)是频率为fx的正弦波电压信号,输入信号u(t)包含电力***中的谐波分量;
采样率设置模块,用于以输入信号的中心频率为中心设置M个采样率,周波采样N次,其中M和N均为大于1的正整数;
采样模块,用于使用M个采样率同时对输入信号u(t)进行采样,采样数据进行缓存;
频率计算模块,用于对输入信号u(t)进行采样过程中,每间隔一个周波计算一次频率;
中间频率获取模块,用于当频率变化时,使用中心频率下的采样数据,按照傅式测频算法计算得出中间频率f1;
采样数据选取模块,用于从M个采样率中选取与中间频率f1最接近的采样率对应的采样数据组;
测频模块,用于根据M个采样率中与中间频率f1最接近的采样率对应的采样数据组,利用傅式测频算法计算最终的频率f。
7.根据权利要求6所述的一种基于傅式快速高精度的电力***频率测量装置,其特征在于,所述输入信号设定模块中谐波含量包括3×f0和5×f0,f0=50hz,f0为输入信号u(t)的中心频率。
8.根据权利要求7所述的一种基于傅式快速高精度的电力***频率测量装置,其特征在于,所述采样率设置模块还用于:
按照额定输入频率fs={44,45,46,47,48,49,50,51,52,53,54,55}设置12个采样率,周波采样256次。
9.根据权利要求8所述的一种基于傅式快速高精度的电力***频率测量装置,其特征在于,所述中间频率获取模块还用于:
使用50hz采用率对输入信号u(t)进行采样过程中,当频率变化时,根据两个相邻数据窗间的相角差利用公式fb=ΔQ/2πΔt求出变化后的频率,其中,ΔQ表示两个相邻数据窗间的相角差,Δt表示两个相邻数据窗间的时间间隔;
每个采样点计算一次角度差,依次利用上述求取变化后的频率的公式计算得到一组临时频率fbuf,对该组临时频率fbuf,进行周期抽取后平均,得到中间频率f1。
10.根据权利要求9所述的一种基于傅式快速高精度的电力***频率测量装置,其特征在于,所述数据窗的长度为一个采样周期。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110713847.4A CN113447713B (zh) | 2021-06-25 | 2021-06-25 | 一种基于傅式快速高精度的电力***频率测量方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110713847.4A CN113447713B (zh) | 2021-06-25 | 2021-06-25 | 一种基于傅式快速高精度的电力***频率测量方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113447713A true CN113447713A (zh) | 2021-09-28 |
CN113447713B CN113447713B (zh) | 2023-03-07 |
Family
ID=77812944
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110713847.4A Active CN113447713B (zh) | 2021-06-25 | 2021-06-25 | 一种基于傅式快速高精度的电力***频率测量方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113447713B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106918741A (zh) * | 2017-03-02 | 2017-07-04 | 浙江大学 | 应用于频率宽范围波动电网的自适应采样的相位差校正法 |
EP3199956A1 (en) * | 2016-01-28 | 2017-08-02 | General Electric Technology GmbH | Apparatus for determination of the frequency of an electrical signal and associated method |
CN107782966A (zh) * | 2016-08-24 | 2018-03-09 | 施耐德电器工业公司 | 确定交流信号的频率 |
CN108333426A (zh) * | 2017-12-25 | 2018-07-27 | 南京丰道电力科技有限公司 | 基于傅氏算法的电力***频率测量方法 |
US20190204370A1 (en) * | 2013-08-16 | 2019-07-04 | Arcteq Relays Oy | Method and system and computer program for measuring alternating-current system quantities |
CN111354365A (zh) * | 2020-03-10 | 2020-06-30 | 苏宁云计算有限公司 | 一种纯语音数据采样率识别方法、装置、*** |
WO2020220804A1 (zh) * | 2019-04-28 | 2020-11-05 | 哈尔滨工业大学 | 一种可变采样率的傅里叶分析方法 |
-
2021
- 2021-06-25 CN CN202110713847.4A patent/CN113447713B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20190204370A1 (en) * | 2013-08-16 | 2019-07-04 | Arcteq Relays Oy | Method and system and computer program for measuring alternating-current system quantities |
EP3199956A1 (en) * | 2016-01-28 | 2017-08-02 | General Electric Technology GmbH | Apparatus for determination of the frequency of an electrical signal and associated method |
CN107782966A (zh) * | 2016-08-24 | 2018-03-09 | 施耐德电器工业公司 | 确定交流信号的频率 |
CN106918741A (zh) * | 2017-03-02 | 2017-07-04 | 浙江大学 | 应用于频率宽范围波动电网的自适应采样的相位差校正法 |
CN108333426A (zh) * | 2017-12-25 | 2018-07-27 | 南京丰道电力科技有限公司 | 基于傅氏算法的电力***频率测量方法 |
WO2020220804A1 (zh) * | 2019-04-28 | 2020-11-05 | 哈尔滨工业大学 | 一种可变采样率的傅里叶分析方法 |
CN111354365A (zh) * | 2020-03-10 | 2020-06-30 | 苏宁云计算有限公司 | 一种纯语音数据采样率识别方法、装置、*** |
Non-Patent Citations (2)
Title |
---|
张义辉: "一种符合电力***动态频率变化的软件测量算法", 《重庆科技学院学报(自然科学版)》 * |
邹正华等: "电能质量监测装置中数据采集模块的设计与实现", 《低压电器》 * |
Also Published As
Publication number | Publication date |
---|---|
CN113447713B (zh) | 2023-03-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107271768B (zh) | 一种最小二乘拟合动态频率测量方法 | |
CN103575980B (zh) | ***频率测量方法、同步相量测量方法及设备 | |
EP2133702A1 (en) | Synchronous phasor measuring device and inter-bus phase angle difference measuring device employing it | |
CN103941088A (zh) | 一种基于三相信号的电力***频率快速测量方法 | |
CN102253282B (zh) | 纳托尔窗函数连续频谱内插电力谐波参数获取方法 | |
CN103018555B (zh) | 一种高精度的电力参数软件同步采样方法 | |
CN108333426A (zh) | 基于傅氏算法的电力***频率测量方法 | |
CN101261293A (zh) | 基于自适应滤波器的电力稳态信号跟踪测量法 | |
Mokeev | Filter synthesis for PMU | |
CN110222309A (zh) | 一种基于鲁棒容积卡尔曼滤波的发电机动态估计方法 | |
CN101572441A (zh) | 一种微机自动准同期功能实现方法 | |
Abdolkhalig et al. | Phasor measurement based on IEC 61850-9-2 and Kalman–Filtering | |
CN102346219B (zh) | 电压源逆变器接入点电压三相软件锁相环相位检测方法 | |
CN104007408A (zh) | 一种pmu动态性能的在线检测方法及装置 | |
CN103884910B (zh) | 一种适用于频率偏移的电力***相量计算方法 | |
CN113447713B (zh) | 一种基于傅式快速高精度的电力***频率测量方法及装置 | |
CN103293379B (zh) | 基于有效值的apf谐波检测方法及其直流侧电压的控制方法 | |
CN115656625A (zh) | 一种基于动态自适应阶梯波的动态信号测量方法及*** | |
EP3199956B1 (en) | Apparatus for determination of the frequency of an electrical signal and associated method | |
CN104198809B (zh) | 一种电力***多频率振荡的频率测量方法 | |
Kunac et al. | Grid Voltage Amplitude and Frequency Real-Time Estimation using Linear Kalman Filter | |
CN110927452B (zh) | 一种基于瞬时无功功率的相位差测量方法及装置 | |
Avalos-Almazan et al. | Real-Time Phasor Estimation via the Taylor-Fourier's Subspace | |
Singh et al. | Implementation of Goertzel algorithm for low cost synchrophasor design | |
CN104833878B (zh) | 基于同步相量测量装置的电力***功率振荡判断方法 |
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 |