CN103884910A - 一种适用于频率偏移的电力***相量计算方法 - Google Patents
一种适用于频率偏移的电力***相量计算方法 Download PDFInfo
- Publication number
- CN103884910A CN103884910A CN201410142911.8A CN201410142911A CN103884910A CN 103884910 A CN103884910 A CN 103884910A CN 201410142911 A CN201410142911 A CN 201410142911A CN 103884910 A CN103884910 A CN 103884910A
- Authority
- CN
- China
- Prior art keywords
- phasor
- frequency
- sequence
- phase
- long
- 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
- Measurement Of Resistance Or Impedance (AREA)
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
本发明公开了一种适用于频率偏移的电力***相量计算方法,包括以下步骤:对三相信号进行同步定时间间隔离散化采样,对采样值序列进行低通滤波,并对***频率进行实时精确测量,得到***频率f;根据不同的***状态以及***频率范围,分别选择定窗长傅里叶变换算法、变窗长傅里叶变换算法、相量修正正序补偿算法或者变窗长正序补偿算法来计算相量的幅值和相角。本发明有益效果:只取一个周期的采样数据,内存占用小;三相数据可以并行运算,计算速度快;精度非常高,相角误差在0.1°以内,幅值误差在0.1%以内。适合于能够得到准确频率值、要求相量测量精度高、计算速度快的场合。该方法适用于当前智能电网对相量的测量要求。
Description
技术领域
本发明属于电力***中的相量测量技术领域,具体涉及一种电力***中适用于频率偏移情况下的高精度相量计算方法。
背景技术
电力***中的相量是一个重要的观测量,它包括幅值和相位。电力***安全自动装置要实时测量其值,并根据电压、电流的相量变化做出相应的控制措施。当前智能化变电站中普遍采用PMU对电力***的电压、电流相量进行测量,计算方法几乎均是基于傅里叶变换来计算幅值和相位。由于电力***实际的频率是波动的,不能持续稳定地保持在工频50Hz,传统的锁相处理技术已经不再适用于一块A/D转换器对频率不同的多路相量进行A/D采样。同时《电力***实时动态监测***技术规范》要求,PMU采样频率是常量,当频率偏移50Hz时,采样频率不是信号频率的整数倍,在这种情况下使用传统傅里叶算法会产生很大的误差。因此必须对傅里叶变换所得结果进行修正,减小测量误差。
随着电力***不断地发展完善,对相量测量的精度要求也大大地提高。而当前的相量算法要么为了满足精度,使傅里叶变换后的修正非常麻烦,计算量很大;要么为了满足计算速度,牺牲测量精度。目前仍没有合适的算法能够在顾及计算成本的基础上,满足实际的相量测量精度,并且适用于严重频率偏移以及各种***运行方式情况。
发明内容
本发明的目的就是为了解决上述问题,提出了一种适用于频率偏移的电力***相量计算方法,本方法利用电力***中的三相采样数据,根据不同的***状态以及频率范围,选择相应的修正公式以及采样窗长,实现高精度的相量快速实时测量。
为了实现上述目的,本发明采用如下技术方案:
一种适用于频率偏移的电力***相量计算方法,包括以下步骤:
(1)对三相信号进行同步定时间间隔离散化采样,对采样值序列进行低通滤波,并对***频率进行实时精确测量,得到***频率f。
(2)根据不同的***状态以及***频率范围,分别选择定窗长傅里叶变换算法、变窗长傅里叶变换算法、相量修正正序补偿算法或者变窗长正序补偿算法来计算相量的幅值和相角。
所述步骤(2)的具体方法为:
***频率值在50.00Hz±0.01Hz范围内时,直接采用离散傅里叶变换算法:将单相采样序列进行离散傅里叶变换,选取采样窗长为工频下每周期固定采样点数N,得到的相量X即为测量结果。
当电力***三相平衡,即***只包含正序分量,且***频率在50.00Hz±0.01Hz所述范围以外时,采用相量修正正序补偿算法:将三相采样序列分别进行离散傅里叶变换,选取采样窗长为工频下每周期固定采样点数N,得到对应同一时刻的三个相量XA’、XB’、XC’,根据三相相量求取正序分量X1;对向量X1进行修正,得到最终的测量结果X。
当电力***三相平衡,即***只包含正序分量,且***频率在50.00Hz±0.01Hz所述范围以外时,还可以采用变窗长离散傅里叶变换算法:将单相采样值序列进行离散傅里叶变换,选取采样窗长为M,得到的相量X即为测量结果,其中M为正整数。
所述步骤(2)的具体方法为:
当电力***三相不平衡,即***包含正序分量、负序分量和/或零序分量,且***频率满足时,采用变窗长离散傅里叶变换算法:将单相采样值序列进行离散傅里叶变换,选取采样窗长为M,得到的相量X即为测量结果,其中M为正整数。
所述步骤(2)的具体方法为:
根据***频率f确定采样窗长M,将三相采样序列分别进行采样窗长为M的离散傅里叶变换,得到对应同一时刻的三个相量XA’、XB’、XC’。
根据三个相量XA’、XB’、XC’求取正序分量X1,对向量X1进行修正,修正后的相量X即为测量结果。
所述将采样序列进行离散傅里叶变换方法为:
其中,N为采样窗长,i为采样点序号,xi+k为第i+k个采样点的值,k为采样窗内的采样点序号,取值从1-N到0。
所述求取正序分量的方法为:假设以A相为基准相,则
其中,XA’、XB’、XC’为三相采样序列分别经离散傅里叶变换得到的向量。
所述对正序向量X1进行修正的方法为:
其中,Δf为***当前频率与额定频率f0的差,即Δf=f-f0;N为工频下每周期固定采样点数。
所述采样窗长M的确定方法为:
其中,f为***频率,f0为额定频率,N为工频下每周期固定采样点数。
本发明的有益效果是:
本发明方法简单,只取一个周期的采样数据,内存占用小;三相数据可以并行运算,计算速度快;精度非常高,相角误差在0.1°以内,幅值误差在0.1%以内。适合应用于能够得到准确频率值、要求相量测量精度高以及计算速度快的场合;适用于当前智能电网对相量的测量要求。
附图说明
图1为本发明的方法流程图。
具体实施方式:
下面结合附图与实施例对本发明做进一步说明:
定义π为圆周率常数;定义运算ea为常数2.71828183的a次方;定义运算cos(a)为对a作余弦运算;定义运算Round(a)为对a先进行四舍五入,再取整运算。
定义整点频率为在定时间间隔离散采样下一个周期正好包含整数个采样点所对应的频率值,也就是一个周期正好包含M个采样点,则采样窗长为M(M为正整数),M对应的整点频率为其中N为工频下每周期采样点数,f0=50.00Hz;例如采样频率为4000Hz,一个周期恰好包含78个采样点,则其对应的整点频率就是51.282Hz。
定义变窗长为在***频率偏移额定频率时,为了使采样窗中的数据更好地反映一个周波的信号,令采样窗长为最接近实际信号周期的整数。例如频率f=45Hz时,一个周包含88.89个采样点,则采样窗长M应取89。
步骤一:对三相信号进行同步定时间间隔(Ts)离散化采样,所得采样值序列先经过低通滤波,以消除原始输入信号中高频分量的干扰;并对***频率进行实时精确测量,得到***频率。
步骤二:当***频率值在50.00Hz±0.01Hz范围内时,直接采用上述离散傅里叶变换法。将单相采样值序列进行离散傅里叶变换,采样窗长为工频下每周期采样点数N,得到的相量X即为测量结果。
步骤三:当***三相平衡,即只包含正序分量;同时***频率在步骤二中所述范围以外时,采用相量修正正序补偿算法,简称正序补偿算法。将三相采样序列各自进行离散傅里叶变换,采样窗长为工频下每周期采样点数N。得到对应同一时刻的三个相量XA’、XB’、XC’,根据三相相量求取正序分量,得到相量XA1。
若以A相作为基准相,则
步骤四:当三相***不平衡,即包含正序分量、负序分量和/或零序分量;同时***频率在整点频率附近,满足时,其中M为正整数,采用变窗长离散傅里叶变换算法。将单相采样值序列进行离散傅里叶变换,采样窗长为M,得到的相量X即为测量结果。
步骤五:当三相***不平衡,即包含正序分量、负序分量和/或零序分量;同时***频率在步骤四中所述范围以外时,采用变窗长与正序补偿算法相结合。首先,根据***频率f确定采样窗长,窗长然后,将三相采样序列各自进行采样窗长为M的离散傅里叶变换,得到对应同一时刻的三个相量XA’、XB’、XC’。最后,同步骤三一样,求取正序分量XA1,再进行修正,修正后的相量XA即为A相的测量结果。
步骤四中的方法在步骤三中***三相平衡的时候也可以使用。
下面实施例中以三相正弦波电压信号为例,具体说明进行相量测量的过程,但本实施例中的待测信号并不限定为电压信号,电力***中的其他三相信号如电流等均可适用。
实施例1:
假设经过低通滤波后的单相采样值序列表示为X(i)=1.34*cos(2*f*π*i*0.00025),i=0,1,2,3,…。N取80,即每周波80点采样,采样间隔TS=0.00025秒。如无特殊说明,计算变量取double型,相位的单位为角度,频率的单位为Hz,小数截断做四舍五入处理。
假设当前***处于三相平衡状态,***频率为47Hz,当前采样时刻是(100×TS),则当f=47Hz时,与此频率最近的整点频率范围是:M=85,47Hz不在整点频率范围内。根据算法规则,选择步骤三。
X(i)A=1.34*cos(2*f*π*i*0.00025)
对X(i)A,X(i)B,X(i)C分别进行离散傅里叶变换,N取80。得到XA’、XB’、XC’三个相量值,利用下式求取以A相为基准相的正序分量:
将XA1再进行修正,得到 修正后的相量值与实际相量值幅值误差为0,相角误差为0,满足测量精度要求。
实施例2:
假设经过低通滤波后的单相采样值序列表示为X(i)=1.34*cos(2*f*π*i*0.00025),i=0,1,2,3,…。N取80,即每周波80点采样,采样间隔TS=0.00025秒。如无特殊说明,计算变量取double型,相位的单位为角度,频率的单位为Hz,小数截断做四舍五入处理。
假设当前***处于三相不平衡状态,存在5%的负序分量,负序超前正序30°。***频率为47Hz,当前采样时刻是(100×TS)。
将XA1再进行修正,得到 修正后的相量值与实际相量值幅值误差小于0.05%,相角误差为0.0009°,满足测量精度要求。
通过对本算法大量仿真实验,在负序分量为10%的情况下,本算法依然能够满足幅值测量最大误差0.1%,相角测量最大误差0.1°的要求。符合当前智能电网对相量的测量要求。
实施例3:
假设经过低通滤波后的单相采样值序列表示为X(i)=1.34*cos(2*f*π*i*0.00025),i=0,1,2,3,…。N取80,即每周波80点采样,采样间隔TS=0.00025秒。如无特殊说明,计算变量取double型,相位的单位为角度,频率的单位为Hz,小数截断做四舍五入处理。
假设当前***频率为50.01Hz,采样时刻是(100×TS),则当f=50.01Hz时,根据算法规则,选择步骤二。选择采样窗长为默认工频每周期采样点数N=80。
实施例4:
假设经过低通滤波后的单相采样值序列表示为X(i)=1.34*cos(2*f*π*i*0.00025),i=0,1,2,3,…。N取80,即每周波80点采样,采样间隔TS=0.00025秒。如无特殊说明,计算变量取double型,相位的单位为角度,频率的单位为Hz,小数截断做四舍五入处理。
假设当前***频率为48.20Hz,当前采样时刻是(100×TS),则当f=48.20Hz时,与此频率最近的整点频率范围是:M=83,48.20Hz在整点频率范围内。根据算法规则,选择步骤四。选择采样窗长为83的离散傅里叶变换。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
Claims (10)
1.一种适用于频率偏移的电力***相量计算方法,其特征是,包括以下步骤:
(1)对三相信号进行同步定时间间隔离散化采样,对采样值序列进行低通滤波,并对***频率进行实时精确测量,得到***频率f;
(2)根据电力***三相是否平衡以及电力***频率范围的不同,分别选择定窗长傅里叶变换算法、变窗长傅里叶变换算法、相量修正正序补偿算法或者变窗长正序补偿算法来计算相量的幅值和相角。
2.如权利要求1所述的一种适用于频率偏移的电力***相量计算方法,其特征是,所述步骤(2)的具体方法为:
***频率值在50.00Hz±0.01Hz范围内时,直接采用离散傅里叶变换算法:将单相采样序列进行离散傅里叶变换,选取采样窗长为工频下每周期固定采样点数N,得到的相量X即为测量结果。
3.如权利要求1所述的一种适用于频率偏移的电力***相量计算方法,其特征是,所述步骤(2)的具体方法为:
当电力***三相平衡,即***只包含正序分量,且***频率在50.00Hz±0.01Hz所述范围以外时,采用相量修正正序补偿算法:将三相采样序列分别进行离散傅里叶变换,选取采样窗长为工频下每周期固定采样点数N,得到对应同一时刻的三个相量XA’、XB’、XC’,根据三相相量求取正序分量X1;对向量X1进行修正,得到最终的测量结果X。
4.如权利要求1或3所述的一种适用于频率偏移的电力***相量计算方法,其特征是,当电力***三相平衡,即***只包含正序分量,且***频率在50.00Hz±0.01Hz所述范围以外时,还可以采用变窗长离散傅里叶变换算法:将单相采样值序列进行离散傅里叶变换,选取采样窗长为M,得到的相量X即为测量结果,其中M为正整数。
7.如权利要求2-6所述的任一种适用于频率偏移的电力***相量计算方法,其特征是,所述将采样序列进行离散傅里叶变换方法为:
其中,N为采样窗长,i为采样点序号,xi+k为第i+k个采样点的值,k为采样窗内的采样点序号,取值从1-N到0。
8.如权利要求3或6所述的一种适用于频率偏移的电力***相量计算方法,其特征是,所述求取正序分量的方法为:
假设以A相为基准相,则
其中,XA’、XB’、XC’为三相采样序列分别经离散傅里叶变换得到的向量。
9.如权利要求3或6所述的一种适用于频率偏移的电力***相量计算方法,其特征是,所述对正序向量X1进行修正的方法为:
其中,Δf为***当前频率与额定频率f0的差,即Δf=f-f0;N为工频下每周期固定采样点数。
10.如权利要求4-6所述的任一种适用于频率偏移的电力***相量计算方法,其特征是,所述采样窗长M的确定方法为:
其中,f为***频率,f0为额定频率,N为工频下每周期固定采样点数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410142911.8A CN103884910B (zh) | 2014-04-10 | 2014-04-10 | 一种适用于频率偏移的电力***相量计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410142911.8A CN103884910B (zh) | 2014-04-10 | 2014-04-10 | 一种适用于频率偏移的电力***相量计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103884910A true CN103884910A (zh) | 2014-06-25 |
CN103884910B CN103884910B (zh) | 2016-06-01 |
Family
ID=50953913
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410142911.8A Expired - Fee Related CN103884910B (zh) | 2014-04-10 | 2014-04-10 | 一种适用于频率偏移的电力***相量计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103884910B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104793053A (zh) * | 2015-04-22 | 2015-07-22 | 福州大学 | 一种基于dft的同步相量相角测量方法 |
CN108614155A (zh) * | 2018-05-31 | 2018-10-02 | 许继集团有限公司 | 一种加入汉明窗的同步相量测量方法及*** |
CN114167134A (zh) * | 2022-02-14 | 2022-03-11 | 青岛鼎信通讯股份有限公司 | 一种低压配电网数据缺失情况下的特征信号强度检测方法 |
CN114461973A (zh) * | 2021-12-23 | 2022-05-10 | 北京四方继保工程技术有限公司 | 一种频率偏移时dft计算相角的补偿方法与*** |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5406495A (en) * | 1993-02-01 | 1995-04-11 | Systems Analysis And Integration, Inc. | Substation load distribution monitor system |
CN1477401A (zh) * | 2003-07-18 | 2004-02-25 | 清华大学 | 高精度的同步相量测量方法 |
CN101587147A (zh) * | 2009-06-25 | 2009-11-25 | 中国电力科学研究院 | 一种同步相量测量装置进行相量修正的方法 |
CN101655519A (zh) * | 2009-09-14 | 2010-02-24 | 国电南京自动化股份有限公司 | 数字化变电站测控装置交流采样数据处理方法 |
CN102128975A (zh) * | 2010-12-22 | 2011-07-20 | 四川省电力公司 | 电压稳定在线监测相量数据测量装置及相量测量方法 |
CN102680785A (zh) * | 2012-04-26 | 2012-09-19 | 中国人民解放军海军工程大学 | 基于自适应变窗的同步相量测量方法 |
-
2014
- 2014-04-10 CN CN201410142911.8A patent/CN103884910B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5406495A (en) * | 1993-02-01 | 1995-04-11 | Systems Analysis And Integration, Inc. | Substation load distribution monitor system |
CN1477401A (zh) * | 2003-07-18 | 2004-02-25 | 清华大学 | 高精度的同步相量测量方法 |
CN101587147A (zh) * | 2009-06-25 | 2009-11-25 | 中国电力科学研究院 | 一种同步相量测量装置进行相量修正的方法 |
CN101655519A (zh) * | 2009-09-14 | 2010-02-24 | 国电南京自动化股份有限公司 | 数字化变电站测控装置交流采样数据处理方法 |
CN102128975A (zh) * | 2010-12-22 | 2011-07-20 | 四川省电力公司 | 电压稳定在线监测相量数据测量装置及相量测量方法 |
CN102680785A (zh) * | 2012-04-26 | 2012-09-19 | 中国人民解放军海军工程大学 | 基于自适应变窗的同步相量测量方法 |
Non-Patent Citations (3)
Title |
---|
刘炳旭等: "基于固定采样频率的相量精确测量算法", 《山东大学学报(工学版)》, vol. 35, no. 6, 31 December 2005 (2005-12-31), pages 46 - 49 * |
姚文轩等: "非同步采样下电力***相量测量修正算法", 《仪器仪表学报》, vol. 34, no. 7, 31 July 2013 (2013-07-31), pages 1518 - 1523 * |
蒋超: "电力***同步相量测量算法的研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》, no. 4, 15 April 2012 (2012-04-15) * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104793053A (zh) * | 2015-04-22 | 2015-07-22 | 福州大学 | 一种基于dft的同步相量相角测量方法 |
CN104793053B (zh) * | 2015-04-22 | 2017-10-20 | 福州大学 | 一种基于dft的同步相量相角测量方法 |
CN108614155A (zh) * | 2018-05-31 | 2018-10-02 | 许继集团有限公司 | 一种加入汉明窗的同步相量测量方法及*** |
CN114461973A (zh) * | 2021-12-23 | 2022-05-10 | 北京四方继保工程技术有限公司 | 一种频率偏移时dft计算相角的补偿方法与*** |
CN114167134A (zh) * | 2022-02-14 | 2022-03-11 | 青岛鼎信通讯股份有限公司 | 一种低压配电网数据缺失情况下的特征信号强度检测方法 |
CN114167134B (zh) * | 2022-02-14 | 2022-04-26 | 青岛鼎信通讯股份有限公司 | 一种低压配电网数据缺失情况下的特征信号强度检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103884910B (zh) | 2016-06-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8198904B2 (en) | Synchrophasor measuring device and inter-bus-line phase angle difference measurement unit using the same | |
CN102435844B (zh) | 一种频率无关的正弦信号相量计算方法 | |
CN102095934B (zh) | 一种交流电信号的相位差测量方法 | |
CN103884910B (zh) | 一种适用于频率偏移的电力***相量计算方法 | |
CN103941088A (zh) | 一种基于三相信号的电力***频率快速测量方法 | |
CN102902879A (zh) | 一种基于fpga硬件dft递推的同步相量计算方法 | |
Mokeev | Filter synthesis for PMU | |
CN104597320A (zh) | 一种适用于多个频率交流信号计算的方法 | |
CN102236048A (zh) | 一种电力***相量频率测量方法 | |
CN103063913A (zh) | 用于傅里叶变换的频率跟踪方法 | |
CN104502707A (zh) | 一种基于三次样条插值的电力***同步相量测量方法 | |
CN103091545A (zh) | 一种频率无关的正弦信号相量半波计算方法 | |
CN108333426A (zh) | 基于傅氏算法的电力***频率测量方法 | |
CN103969508B (zh) | 一种实时高精密的电力谐波分析方法及装置 | |
CN102095935B (zh) | 一种交流电信号的瞬时相位测量方法 | |
CN102170143A (zh) | 一种微机自动准同期并列实现方法 | |
CN103293379B (zh) | 基于有效值的apf谐波检测方法及其直流侧电压的控制方法 | |
Bliznyuk et al. | The Technique for Accelerated Evaluation of the Instantaneous Values of Power System Performance Parameters | |
CN110927452B (zh) | 一种基于瞬时无功功率的相位差测量方法及装置 | |
CN107359872B (zh) | 一种基于动态调频技术的多路软件锁相装置及锁相方法 | |
Kunac et al. | Grid Voltage Amplitude and Frequency Real-Time Estimation using Linear Kalman Filter | |
Zhang et al. | An innovative timestamp-based convolution integral method in synchrophasor estimation within digital relays | |
Rens et al. | The effect of real network phase disturbances on the calculation of IEC 61000-4-30 parameters | |
WO2021082036A1 (zh) | 一种电力***频率测量方法、母线电压校正方法及装置 | |
Kitzig et al. | Proof of Concept for Three-Phase Extension of a Mains Frequency Estimation Algorithm |
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: 20160601 Termination date: 20210410 |
|
CF01 | Termination of patent right due to non-payment of annual fee |