CN105259409A - 一种基于过零时域特性的电力信号频率计算方法 - Google Patents
一种基于过零时域特性的电力信号频率计算方法 Download PDFInfo
- Publication number
- CN105259409A CN105259409A CN201510664622.9A CN201510664622A CN105259409A CN 105259409 A CN105259409 A CN 105259409A CN 201510664622 A CN201510664622 A CN 201510664622A CN 105259409 A CN105259409 A CN 105259409A
- Authority
- CN
- China
- Prior art keywords
- zero crossing
- power signal
- zero
- actual
- electric power
- 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
Landscapes
- Measuring Frequencies, Analyzing Spectra (AREA)
- Remote Monitoring And Control Of Power-Distribution Networks (AREA)
Abstract
本发明提供一种基于过零时域特性的电力信号频率计算方法,包括:数据预处理的步骤;寻找实际过零点的步骤;寻找理论过零点的步骤;计算电力信号频率的步骤。本发明基于电力信号频率的稳态特性和过零时域特性,能够提高含高次谐波的电力信号频率的计算精度,满足了供配电系统的故障诊断和用户端的电能平衡测试及节电方案设计等对频率精度的要求。
Description
技术领域
本发明涉及电力信号频率计算技术领域,具体是一种基于过零时域特性的电力信号频率计算方法。
背景技术
电力信号的频率稳态特性指的是电力信号稳态频率变化持续时间需要4~10s(参见文献《现代电网频率控制应用技术》,高翔等编写,2010年,中国电力出版社,P30与P141),也即可以认为在4~10s这段时间内电力信号频率是恒定的。
电力信号的过零时域特性指的是(1)不含谐波的电力信号的实际过零点只有一个(从波峰到相邻的波谷或者从波谷到相邻的波峰这个过程),或者由正到负,或者由负到正,理论过零点与实际过零点个数相同,如图1所示;(2)不含谐波的电力信号在两个理论过零点之间,不会再出现实际过零点,如图2所示;(3)含高次谐波的电力信号的实际过零点不止一个(从波峰到相邻的波谷或者从波谷到相邻的波峰这个过程),但不管有多少个,其个数一定是奇数,理论过零点只有一个,如图3所示;含高次谐波的电力信号在两个理论过零点之间,可能会因谐波含量过大而出现实际过零点,其个数一定是偶数,如图4所示。
现有的电力信号综合分析软件、电能质量测试分析系统在分析含高次谐波的电力信号时,频率计算误差较大,不能满足用户的需求。
发明内容
本发明的目的在于提供一种基于过零时域特性的电力信号频率计算方法,提高含高次谐波的电力信号频率的计算精度。
本发明的技术方案为:
一种基于过零时域特性的电力信号频率计算方法,包括以下步骤:
(1)在采样时间q*T内对待测电力信号数据进行均匀采样,获得q*N个采样点xi,其中,q表示小于200的自然数,T=0.02s,N表示每个T内的采样点数,xi表示位置编号为i的采样点的电力信号值,i=0、1、2、…、q*N-1;
(2)依次对两相邻采样点的电力信号值进行比较,若两者的数值符号相反,则将两者中的绝对值较小者作为实际过零点,并依次提取所有实际过零点的位置编号,构成实际过零点集合;
(3)基于电力信号的过零时域特性,从实际过零点中寻找理论过零点,并依次提取所有理论过零点的位置编号,构成理论过零点集合;
(4)采用以下公式计算待测电力信号的频率:
其中,f表示待测电力信号的频率,M≥2,表示理论过零点的个数,n1表示理论过零点集合中的第一个元素,即第一个理论过零点的位置编号,nM表示理论过零点集合中的最后一个元素,即最后一个理论过零点的位置编号,△ts=T/N,表示采样间隔。
所述的基于过零时域特性的电力信号频率计算方法,所述步骤(1)中,q=55,N=256。
所述的基于过零时域特性的电力信号频率计算方法,所述步骤(3)中,所述基于电力信号的过零时域特性,从实际过零点中寻找理论过零点,具体包括:
a、在实际过零点集合中,依次判断两相邻实际过零点的位置编号之间的绝对差值是否大于间隔限值,若是,则说明相应的两相邻实际过零点属于不同的理论过零点,若否,则说明相应的两相邻实际过零点属于同一理论过零点;
b、属于同一理论过零点的所有实际过零点构成一个组,依次判断每组实际过零点的个数是偶数还是奇数,若是偶数,则说明该组实际过零点处不存在理论过零点,若是奇数,则说明该组实际过零点处存在理论过零点,将该组实际过零点的位置编号的平均值取整数,作为相应理论过零点的位置编号。
所述的基于过零时域特性的电力信号频率计算方法,所述步骤a中,所述间隔限值由以下公式给出:
其中,L表示两相邻实际过零点之间的间隔限值,N表示每个T内的采样点数,{}表示取整数。
由上述技术方案可知,本发明基于电力信号频率的稳态特性和过零时域特性,能够提高含高次谐波的电力信号频率的计算精度,满足了供配电系统的故障诊断和用户端的电能平衡测试及节电方案设计等对频率精度的要求。
附图说明
图1是标准正弦电力信号实际过零点与理论过零点,实线圆圈中存在一个理论过零点;从波峰到相邻的波谷或者从波谷到相邻的波峰,标准正弦电力信号只存在一个实际过零点,并且与理论过零点重合;
图2是标准正弦电力信号两个理论过零点之间无实际过零点,实线圆圈中存在一个理论过零点;虚线圆圈中,即两个理论过零点之间,无实际过零点;
图3是含高次谐波的电力信号的实际过零点与理论过零点,实线圆圈中存在一个理论过零点;从波峰到相邻的波谷或者从波谷到相邻的波峰,存在奇数个实际过零点;
图4是含高次谐波的电力信号两个理论过零点之间有偶数个实际过零点,实线圆圈中存在一个理论过零点;虚线圆圈中,即两个理论过零点之间,存在偶数个实际过零点;
图5是本发明的方法流程图;
图6是本发明具体实施例的程序实现流程图。
具体实施方式
如图5所示,一种基于过零时域特性的电力信号频率计算方法,包括以下步骤:
S1、数据预处理:
标准正弦电力信号的周期T=0.02s,根据电力信号的频率稳态特性,通常取1s左右的数据进行分析,本实施例中,将采样时间取为1.1s,即55*T。设每个T内的采样点数为N,则在采样时间55*T内对待测电力信号数据进行均匀采样,可以获得55*N个采样点,记为xi,i=0、1、2、…、55*N-1,表示位置编号为i的采样点的电力信号值,N通常取值为256。
S2、寻找实际过零点:
依次对两相邻采样点的电力信号值进行比较,若某两相邻采样点的数值符号相反,即一个为正,一个为负,则说明该两相邻采样点处存在实际过零点,将该两相邻采样点中的绝对值较小者作为实际过零点,依次将所有实际过零点的位置编号存放在实际过零点集合中。
注:实际过零点集合中存放的是位置编号,也即步骤S1中的i的取值,如若经过判断发现x64、x67、x69为实际过零点,则依次将64、67、69存放在实际过零点集合中。
S3、寻找理论过零点:
在实际过零点集合中,依次判断两相邻实际过零点的位置编号之间的绝对差值l是否大于间隔限值L,若是,则说明相应的两相邻实际过零点属于不同的理论过零点,若否,则说明相应的两相邻实际过零点属于同一理论过零点。
两相邻实际过零点之间的间隔限值L由以下公式给出:
其中,N表示每个T内的采样点数,{}表示取整数。
属于同一理论过零点的所有实际过零点构成一个组,依次判断每组实际过零点的个数是偶数还是奇数,若是偶数,则说明该组实际过零点处不存在理论过零点,若是奇数,则说明该组实际过零点处存在理论过零点,将该组实际过零点的位置编号的平均值取整数,作为相应理论过零点的位置编号。然后,依次将所有理论过零点的位置编号存放在理论过零点集合中。
S4、计算电力信号频率:
通过以上步骤S1~S3,在55*N个采样点中,若找到的理论过零点个数为M(M≥2),则可以采用以下公式计算待测电力信号的频率f:
其中,n1表示理论过零点集合中的第一个元素,即第一个理论过零点的位置编号,nM表示理论过零点集合中的最后一个元素,即最后一个理论过零点的位置编号,△ts=0.02/N,表示采样间隔。
本发明具体实施例的程序实现,如图6所示。
以上所述实施方式仅仅是对本发明的优选实施方式进行描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通技术人员对本发明的技术方案作出的各种变形和改进,均应落入本发明的权利要求书确定的保护范围内。
Claims (4)
1.一种基于过零时域特性的电力信号频率计算方法,其特征在于,包括以下步骤:
(1)在采样时间q*T内对待测电力信号数据进行均匀采样,获得q*N个采样点xi,其中,q表示小于200的自然数,T=0.02s,N表示每个T内的采样点数,xi表示位置编号为i的采样点的电力信号值,i=0、1、2、…、q*N-1;
(2)依次对两相邻采样点的电力信号值进行比较,若两者的数值符号相反,则将两者中的绝对值较小者作为实际过零点,并依次提取所有实际过零点的位置编号,构成实际过零点集合;
(3)基于电力信号的过零时域特性,从实际过零点中寻找理论过零点,并依次提取所有理论过零点的位置编号,构成理论过零点集合;
(4)采用以下公式计算待测电力信号的频率:
其中,f表示待测电力信号的频率,M≥2,表示理论过零点的个数,n1表示理论过零点集合中的第一个元素,即第一个理论过零点的位置编号,nM表示理论过零点集合中的最后一个元素,即最后一个理论过零点的位置编号,△ts=T/N,表示采样间隔。
2.根据权利要求1所述的基于过零时域特性的电力信号频率计算方法,其特征在于,所述步骤(1)中,q=55,N=256。
3.根据权利要求1所述的基于过零时域特性的电力信号频率计算方法,其特征在于,所述步骤(3)中,所述基于电力信号的过零时域特性,从实际过零点中寻找理论过零点,具体包括:
a、在实际过零点集合中,依次判断两相邻实际过零点的位置编号之间的绝对差值是否大于间隔限值,若是,则说明相应的两相邻实际过零点属于不同的理论过零点,若否,则说明相应的两相邻实际过零点属于同一理论过零点;
b、属于同一理论过零点的所有实际过零点构成一个组,依次判断每组实际过零点的个数是偶数还是奇数,若是偶数,则说明该组实际过零点处不存在理论过零点,若是奇数,则说明该组实际过零点处存在理论过零点,将该组实际过零点的位置编号的平均值取整数,作为相应理论过零点的位置编号。
4.根据权利要求3所述的基于过零时域特性的电力信号频率计算方法,其特征在于,所述步骤a中,所述间隔限值由以下公式给出:
其中,L表示两相邻实际过零点之间的间隔限值,N表示每个T内的采样点数,{}表示取整数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510664622.9A CN105259409B (zh) | 2015-10-13 | 2015-10-13 | 一种基于过零时域特性的电力信号频率计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510664622.9A CN105259409B (zh) | 2015-10-13 | 2015-10-13 | 一种基于过零时域特性的电力信号频率计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105259409A true CN105259409A (zh) | 2016-01-20 |
CN105259409B CN105259409B (zh) | 2017-11-28 |
Family
ID=55099175
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510664622.9A Active CN105259409B (zh) | 2015-10-13 | 2015-10-13 | 一种基于过零时域特性的电力信号频率计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105259409B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107817421A (zh) * | 2017-10-24 | 2018-03-20 | 南京航天银山电气有限公司 | 一种波形频率的快速稳定判断模块和判断方法 |
CN108195021A (zh) * | 2017-12-29 | 2018-06-22 | 奥克斯空调股份有限公司 | 一种电机电源频率判断方法、装置及空调 |
CN108957133A (zh) * | 2018-08-13 | 2018-12-07 | 国网浙江省电力有限公司绍兴供电公司 | 变电站电容器采样频率动态评估和选取方法 |
CN109342814A (zh) * | 2018-09-05 | 2019-02-15 | 苏州光格设备有限公司 | 用于确定电缆局放工频相位频率的方法、装置及存储介质 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003098203A (ja) * | 2001-09-20 | 2003-04-03 | Asahi Kasei Microsystems Kk | 群遅延差測定方法およびその測定装置 |
CN1815248A (zh) * | 2005-02-02 | 2006-08-09 | 艾默生网络能源系统有限公司 | 交流电频率监测方法 |
CN101871965A (zh) * | 2010-06-11 | 2010-10-27 | 威胜集团有限公司 | 电力正弦信号过零时间、频率、相位差的检测方法 |
CN102095929A (zh) * | 2010-12-17 | 2011-06-15 | 南京邮电大学 | 一种快速测量交流电信号频率的方法 |
CN102116798A (zh) * | 2011-03-07 | 2011-07-06 | 深圳市锐能微科技有限公司 | 一种电网频率测量方法及装置 |
CN102879639A (zh) * | 2012-09-13 | 2013-01-16 | 华中科技大学 | 一种电力系统中频率的实时测量方法 |
CN103197129A (zh) * | 2013-04-10 | 2013-07-10 | 安徽节源节能科技有限公司 | 电力信号过零点检测方法 |
CN103575979A (zh) * | 2012-07-26 | 2014-02-12 | 南京邮电大学 | 一种数字化测量交流电频率的方法 |
CN104360144A (zh) * | 2014-12-08 | 2015-02-18 | 广东美的环境电器制造有限公司 | 确定交流信号过零点的方法及系统 |
-
2015
- 2015-10-13 CN CN201510664622.9A patent/CN105259409B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003098203A (ja) * | 2001-09-20 | 2003-04-03 | Asahi Kasei Microsystems Kk | 群遅延差測定方法およびその測定装置 |
CN1815248A (zh) * | 2005-02-02 | 2006-08-09 | 艾默生网络能源系统有限公司 | 交流电频率监测方法 |
CN101871965A (zh) * | 2010-06-11 | 2010-10-27 | 威胜集团有限公司 | 电力正弦信号过零时间、频率、相位差的检测方法 |
CN102095929A (zh) * | 2010-12-17 | 2011-06-15 | 南京邮电大学 | 一种快速测量交流电信号频率的方法 |
CN102116798A (zh) * | 2011-03-07 | 2011-07-06 | 深圳市锐能微科技有限公司 | 一种电网频率测量方法及装置 |
CN103575979A (zh) * | 2012-07-26 | 2014-02-12 | 南京邮电大学 | 一种数字化测量交流电频率的方法 |
CN102879639A (zh) * | 2012-09-13 | 2013-01-16 | 华中科技大学 | 一种电力系统中频率的实时测量方法 |
CN103197129A (zh) * | 2013-04-10 | 2013-07-10 | 安徽节源节能科技有限公司 | 电力信号过零点检测方法 |
CN104360144A (zh) * | 2014-12-08 | 2015-02-18 | 广东美的环境电器制造有限公司 | 确定交流信号过零点的方法及系统 |
Non-Patent Citations (1)
Title |
---|
孙晖 等: "基于过零点一极点估计的瞬时频率幅度算法", 《电子与信息学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107817421A (zh) * | 2017-10-24 | 2018-03-20 | 南京航天银山电气有限公司 | 一种波形频率的快速稳定判断模块和判断方法 |
CN108195021A (zh) * | 2017-12-29 | 2018-06-22 | 奥克斯空调股份有限公司 | 一种电机电源频率判断方法、装置及空调 |
CN108957133A (zh) * | 2018-08-13 | 2018-12-07 | 国网浙江省电力有限公司绍兴供电公司 | 变电站电容器采样频率动态评估和选取方法 |
CN109342814A (zh) * | 2018-09-05 | 2019-02-15 | 苏州光格设备有限公司 | 用于确定电缆局放工频相位频率的方法、装置及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN105259409B (zh) | 2017-11-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wichakool et al. | Smart metering of variable power loads | |
CN107025365B (zh) | 一种用于用户侧的非侵入负荷辨识方法 | |
CN105259409A (zh) | 一种基于过零时域特性的电力信号频率计算方法 | |
CN101539596A (zh) | 电网频率监测的方法 | |
CN106154040B (zh) | 一种用于计算电网单点等值时等效谐波阻抗方法 | |
CN104091056B (zh) | 一种实验室设备使用寿命预测系统和方法 | |
CN103983849B (zh) | 一种实时高精度的电力谐波分析方法 | |
CN106093565A (zh) | 一种基于稳态特征波形匹配的电能分项计量方法及装置 | |
Lee et al. | New power quality index in a distribution power system by using RMP model | |
CN104407198A (zh) | 一种用于dvr装置中检测sag信号的方法和系统 | |
CN108333426A (zh) | 基于傅氏算法的电力系统频率测量方法 | |
Vanfretti et al. | A phasor measurement unit based fast real-time oscillation detection application for monitoring wind-farm-to-grid sub–synchronous dynamics | |
Najafi et al. | An improved sag detection approach based on modified Goertzel algorithm | |
CN103543331B (zh) | 一种计算电信号谐波和间谐波的方法 | |
CN108982954B (zh) | 适用于馈线终端的计算相电压幅值与相位的方法及系统 | |
CN110389251A (zh) | 一种用于电网电压跌落检测的瞬时电压dq分解方法 | |
Tan et al. | High‐accurate estimation method of sub‐synchronous and super‐synchronous harmonic parameters on PMU | |
CN101399444A (zh) | 电容器投切对动态电压波动影响的特征评价方法 | |
CN103267896B (zh) | 周期信号初相角的测量方法 | |
Lin | Development of interharmonics identification using enhanced‐fft algorithm | |
Xian-chun et al. | Design of three-phase multi-purpose standard electric energy meter | |
CN108229095A (zh) | 油中溶解气体体积分数的预测方法及终端设备 | |
Danilov et al. | Operational Identification of Resistances of Wires of 380 V Distribution Networks by Automated Accounting Systems | |
CN110389252A (zh) | 一种用于电网电压跌落的αβ检测方法 | |
Plienis et al. | Design of IOT-Based Framework for Evaluation of Energy Efficiency in Power Transformers |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |