CN107045081A - 一种保护装置中实时频率的测量方法 - Google Patents

一种保护装置中实时频率的测量方法 Download PDF

Info

Publication number
CN107045081A
CN107045081A CN201710043217.4A CN201710043217A CN107045081A CN 107045081 A CN107045081 A CN 107045081A CN 201710043217 A CN201710043217 A CN 201710043217A CN 107045081 A CN107045081 A CN 107045081A
Authority
CN
China
Prior art keywords
frequency
real
angular frequency
electric signals
time
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
Application number
CN201710043217.4A
Other languages
English (en)
Other versions
CN107045081B (zh
Inventor
羊阳
李绪勇
黄真真
汪战魁
张瑞杰
汪相坤
陈永胜
王启文
寇文强
王洲洲
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Henan Intelligent Science And Technology Co Ltd That Permitted To Continue
Original Assignee
Henan Intelligent Science And Technology Co Ltd That Permitted To Continue
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Henan Intelligent Science And Technology Co Ltd That Permitted To Continue filed Critical Henan Intelligent Science And Technology Co Ltd That Permitted To Continue
Priority to CN201710043217.4A priority Critical patent/CN107045081B/zh
Publication of CN107045081A publication Critical patent/CN107045081A/zh
Application granted granted Critical
Publication of CN107045081B publication Critical patent/CN107045081B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/02Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage

Abstract

本发明公开一种保护装置中实时频率的测量方法,包括以下步骤:步骤1)、对被测电信号的整周波进行固定间隔采样;步骤2)、得到计算公式;步骤3)、进行N/2点数累加上述公式后半部分的值;步骤4)、推算得到角频率偏差;步骤5)、由角频率偏差计算得到实时角频率,转换后得到被测电信号的实时频率。本发明的实时频率测量方法是利用工频信号正弦波形特征,通过固定间隔采样值,进行曲线拟合累加,在固定工频基础上推导出频差,从而得出实际频率值,本发明的测量方法用较少的采样点数就可算出实际频率值,且中间没有较为复杂的运算逻辑,实现起来简单高效。实际应用中测频精度完全满足电力系统继电保护和测量的需要。

Description

一种保护装置中实时频率的测量方法
技术领域
本发明涉及电力系统参数测量技术领域,具体涉及一种保护装置中实时频率的测量方法。
背景技术
频率是电力系统的主要性能参数,国际上曾多次发生频率崩溃导致的大面积停电事故。电力系统功率不足将导致频率下降,通过自动发电系统来启动旋转备,使电力系统频率恢复正常。频率严重不足时,低频减负荷装置切除部分次要负荷,保持系统频率和对重要用户的供电。实施这些措施的前提是相应设备和装置能准确测量到系统的实时频率。系统实时频率是同步相量测量装置的主要参数,是系统稳定性判别和控制决策的依据之一。另外,系统发生异常或故障时,继电保护装置根据实时频率可以准确计算相关的电气量参数,如序分量、故障分量和阻抗等。
频率的测量主要分为硬件测量和软件测量,硬件测量频率一般由过零比较电路、方波形成电路和计数器构成,简单、方便,但需要增加硬件电路,不但增加硬件成本而且易受干扰。软件测频包括基于插值的cross算法、最小二乘方法、递推富氏算法、全波富氏算法、卡尔曼滤波算法等,但都存在这样那样的不足与缺陷。
发明内容
本发明的目的是提供一种保护装置中实时频率的测量方法,较少的采样点数就可算出实际频率值,且中间没有较为复杂的运算逻辑,实现起来简单高效。
本发明实现上述技术目的采用的技术方案是:一种保护装置中实时频率的测量方法,包括以下步骤:
步骤1)、对被测电信号的整周波进行固定间隔采样,每周波采样点数N,采样间隔T,得到被测电信号的固定间隔采样值{U0,U1,……,UN-1,UN,……,U2N-1,U2N,U2N+1,……;
步骤2)、定义工频信号为以下表示方式:
其中,S为幅值,w为角频率,△w角频率偏差,w 0 固定波角频率,初始相角;
根据上述公式进行推导变化后,得到:
步骤3)、进行N/2点数累加上述公式后半部分的值,为避免分母累加值太小,计算相差90度的2组值,取分母较大的那一组进行计算,假设得出的两组分母的累积加值分别是M1=UN/4 + UN/4+1 + …… + U3N/4-1和M2=UN/2 + UN/2+1 + …… + UN-1,所累积对应的分子是N1=U0 + U1 + …… + UN-1和N2=UN/4 + UN/4+1 + …… + U 5N/4-1。
步骤4)、对所得两组分母进行比较, 取分母较大的那一组进行计算:
if(M1 > M2),则:M = N1/M1;
else,则:M = N2/M2;
即可得出:,从而推算得到角频率偏差=
步骤5)、由可得到实时角频率,通过可计算得到被测电信号的实时频率。
有益效果
本发明的实时频率测量方法是利用工频信号正弦波形特征,通过固定间隔采样值,进行曲线拟合累加,在固定工频基础上推导出频差,从而得出实际频率值,本发明的测量方法用较少的采样点数就可算出实际频率值,且中间没有较为复杂的运算逻辑,实现起来简单高效。实际应用中测频精度完全满足电力系统继电保护和测量的需要。
具体实施方式
下面通过实施例对本发明做进一步详细说明,这些实施例仅用来说明本发明,并不限制本发明的范围。
一种保护装置中实时频率的测量方法,包括以下步骤:
步骤1)、对被测电信号(电力系统50赫兹电压信号)的整周波进行固定间隔采样,每周波采样点数N,采样间隔T,得到被测电信号的固定间隔采样值{U0,U1,……,UN-1,UN,……,U2N-1,U2N,U2N+1,……;
步骤2)、定义工频信号为以下表示方式:
根据上述公式进行推导变化后,得到:
步骤3)、进行N/2点数累加上述公式后半部分的值,为避免分母累加值太小,计算相差90度的2组值,取分母较大的那一组进行计算。假设得出的两组分母的累积加值分别是M1=UN/4 + UN/4+1 + …… + U3N/4-1和M2=UN/2 + UN/2+1 + …… + UN-1,所累积对应的分子是N1=U0 + U1 + …… + UN-1和N2=UN/4 + UN/4+1 + …… + U 5N/4-1。
步骤4)、对所得两组分母进行比较, 取分母较大的那一组进行计算:
if(M1 > M2)
{
M = N1/M1;
}
else
{
M = N2/M2
}
即可得出:, 从而推算得到角频率偏差=
步骤5)、由可得到实时角频率,通过可计算得到被测电信号的实时频率。
实际应用中,电力系统中工频信号频率一般是50HZ,Δw实际角频率与50赫兹角频率的差,N为50赫兹周波采样点数,T和N的值根据实际进行设置(本例采样24点采样),通过下面公式:
进行N/2点数累加,这样,分子相当于一周的累加,分母是半周的累加,有一定的滤波效果。为避免分母累加值太小,计算相差90度的2组值,取分母较大的那一组进行计算。
该算法的精度依托于N/2累加点数的截取,所以T固定周波采样间隔精度,N固定周波采样总点数对测频精度有较大的影响,但由于嵌入式技术的近年来突飞猛进的发展,模拟采样精度再不断提高,这种要求在普通单片机中已很容易实现;同时考虑到谐波分量对测频精度的影响,为此添加了常规带通滤波器算法,可很好的滤掉直流分量和谐波分量,进一步提高了此算法的可靠性和准确性。
以上所述,仅是本发明的较佳实施例而已,并非对本发明作任何形式上的限制,虽然本发明已以较佳实施例揭露如上,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的技术内容做出出些许更动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。

Claims (1)

1.一种保护装置中实时频率的测量方法,其特征在于:包括以下步骤:
步骤1)、对被测电信号的整周波进行固定间隔采样,每周波采样点数N,采样间隔T,得到被测电信号的固定间隔采样值{U0,U1,……,UN-1,UN,……,U2N-1,U2N,U2N+1,……;
步骤2)、定义工频信号为以下表示方式:
其中,S为幅值,w为角频率,△w角频率偏差,w 0 固定波角频率,初始相角;
根据上述公式进行推导变化后,得到:
步骤3)、进行N/2点数累加上述公式后半部分的值,为避免分母累加值太小,计算相差90度的2组值,取分母较大的那一组进行计算,假设得出的两组分母的累积加值分别是M1=UN/4 + UN/4+1 + …… + U3N/4-1和M2=UN/2 + UN/2+1 + …… + UN-1,所累积对应的分子是N1=U0 + U1 + …… + UN-1和N2=UN/4 + UN/4+1 + …… + U 5N/4-1
步骤4)、对所得两组分母进行比较, 取分母较大的那一组进行计算:
if(M1 > M2),则:M = N1/M1;
else,则:M = N2/M2;
即可得出:,从而推算得到角频率偏差
步骤5)、由可得到实时角频率,通过可计算得到被测电信号的实时频率。
CN201710043217.4A 2017-01-21 2017-01-21 一种保护装置中实时频率的测量方法 Active CN107045081B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710043217.4A CN107045081B (zh) 2017-01-21 2017-01-21 一种保护装置中实时频率的测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710043217.4A CN107045081B (zh) 2017-01-21 2017-01-21 一种保护装置中实时频率的测量方法

Publications (2)

Publication Number Publication Date
CN107045081A true CN107045081A (zh) 2017-08-15
CN107045081B CN107045081B (zh) 2020-02-11

Family

ID=59543412

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710043217.4A Active CN107045081B (zh) 2017-01-21 2017-01-21 一种保护装置中实时频率的测量方法

Country Status (1)

Country Link
CN (1) CN107045081B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109725197A (zh) * 2018-12-25 2019-05-07 西电通用电气自动化有限公司 一种多处理单元继电保护系统频率跟踪方法及系统
CN110187174A (zh) * 2019-06-28 2019-08-30 北京智芯微电子科技有限公司 快速定点频偏测量方法
CN110297127A (zh) * 2019-05-28 2019-10-01 许昌许继软件技术有限公司 一种交流信号滤波方法及装置
CN111257634A (zh) * 2020-03-26 2020-06-09 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 一种电网频率实时高精度测量新方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0777547A (ja) * 1993-07-16 1995-03-20 Fuji Electric Co Ltd 正弦波交流信号の周波数検出方法
JPH08184618A (ja) * 1994-12-27 1996-07-16 Fuji Electric Co Ltd 正弦波交流信号の周波数検出方法
CN101241151A (zh) * 2008-02-28 2008-08-13 江苏省电力试验研究院有限公司 基于幅值线性变化模型的实时频率测量方法
CN103983847A (zh) * 2014-06-12 2014-08-13 福州大学 一种同步相量测量中基于rls的自适应频率跟踪测量方法
CN105044456A (zh) * 2015-07-21 2015-11-11 电子科技大学 一种基于正交子带的电网瞬时频率测量与跟踪方法
CN105629060A (zh) * 2015-12-24 2016-06-01 电子科技大学 基于最优基带滤波的电网频率测量方法和装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0777547A (ja) * 1993-07-16 1995-03-20 Fuji Electric Co Ltd 正弦波交流信号の周波数検出方法
JPH08184618A (ja) * 1994-12-27 1996-07-16 Fuji Electric Co Ltd 正弦波交流信号の周波数検出方法
CN101241151A (zh) * 2008-02-28 2008-08-13 江苏省电力试验研究院有限公司 基于幅值线性变化模型的实时频率测量方法
CN103983847A (zh) * 2014-06-12 2014-08-13 福州大学 一种同步相量测量中基于rls的自适应频率跟踪测量方法
CN105044456A (zh) * 2015-07-21 2015-11-11 电子科技大学 一种基于正交子带的电网瞬时频率测量与跟踪方法
CN105629060A (zh) * 2015-12-24 2016-06-01 电子科技大学 基于最优基带滤波的电网频率测量方法和装置

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109725197A (zh) * 2018-12-25 2019-05-07 西电通用电气自动化有限公司 一种多处理单元继电保护系统频率跟踪方法及系统
CN110297127A (zh) * 2019-05-28 2019-10-01 许昌许继软件技术有限公司 一种交流信号滤波方法及装置
CN110187174A (zh) * 2019-06-28 2019-08-30 北京智芯微电子科技有限公司 快速定点频偏测量方法
CN111257634A (zh) * 2020-03-26 2020-06-09 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 一种电网频率实时高精度测量新方法
CN111257634B (zh) * 2020-03-26 2022-01-07 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 一种电网频率实时高精度测量新方法

Also Published As

Publication number Publication date
CN107045081B (zh) 2020-02-11

Similar Documents

Publication Publication Date Title
CN107045081A (zh) 一种保护装置中实时频率的测量方法
US8532230B2 (en) Instantaneous phasor determination for poly-phase electrical grids
JP5806273B2 (ja) インバータの直流リンクコンデンサの容量推定装置
CN104360137B (zh) 一种适用于动态电压恢复器的电压暂降检测方法
CN103207308B (zh) 避雷器阻性电流和容性电流暂态值的测量方法
CN102735938A (zh) 一种电网电压基波正序相角的快速检测方法
CN103983847B (zh) 一种同步相量测量中基于rls的自适应频率跟踪测量方法
JP6416072B2 (ja) 同期フェーザ測定装置およびパルス生成装置
CN105680854B (zh) 一种锁相环及锁相方法
CN109067393A (zh) 一种电力系统的锁相方法、装置及设备
US10677834B2 (en) Distance protection of electric power delivery systems using time domain and frequency domain
CN102221639A (zh) 正负序电流实时检测的方法
CN104502707A (zh) 一种基于三次样条插值的电力系统同步相量测量方法
CN110007183B (zh) 逆变器在线开路故障检测方法
CN104333027A (zh) 预测型换相失败预防方法、装置和设备
CN105021872A (zh) 一种电网电压畸变不对称状态下电流各分量有功成分和无功成分的检测方法
CN103207307B (zh) 金属氧化物避雷器泄漏电流组成成分测量方法
CN105006828A (zh) 一种并网变流器负序电流补偿方法
CN107271772B (zh) 一种高精度且抗噪声干扰的电网频率快速检测方法
CN109031029B (zh) 一种基于单相接地故障的标准电压获取方法及装置
CN104333029B (zh) 预测型换相失败预防方法与装置
Cacciato et al. Multi-reference frame based PLL for single phase systems in voltage distorted grids
WO2021082036A1 (zh) 一种电力系统频率测量方法、母线电压校正方法及装置
CN104062545A (zh) 新型短路故障源快速定位方法
CN104280601B (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
CB02 Change of applicant information
CB02 Change of applicant information

Address after: New energy industrial park in Xuchang County of Xuchang City, Henan province 461000 China Electric Valley XJ Group

Applicant after: Xuchang smart relay Limited by Share Ltd

Address before: New energy industrial park in Xuchang County of Xuchang City, Henan province 461000 China Electric Valley XJ Group

Applicant before: The Henan intelligent Science and Technology Co., Ltd. that permitted to continue

CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Zhang Hongtao

Inventor after: Kou Wenqiang

Inventor after: Wang Zhouzhou

Inventor after: Yang Yang

Inventor after: Li Xuyong

Inventor after: Huang Zhenzhen

Inventor after: Wang Zhankui

Inventor after: Zhang Ruijie

Inventor after: Wang Xiangkun

Inventor after: Chen Yongsheng

Inventor after: Wang Qiwen

Inventor before: Yang Yang

Inventor before: Wang Zhouzhou

Inventor before: Li Xuyong

Inventor before: Huang Zhenzhen

Inventor before: Wang Zhankui

Inventor before: Zhang Ruijie

Inventor before: Wang Xiangkun

Inventor before: Chen Yongsheng

Inventor before: Wang Qiwen

Inventor before: Kou Wenqiang

GR01 Patent grant
GR01 Patent grant