CN107064630B - 一种电力系统频率测量方法及装置 - Google Patents
一种电力系统频率测量方法及装置 Download PDFInfo
- Publication number
- CN107064630B CN107064630B CN201710208460.7A CN201710208460A CN107064630B CN 107064630 B CN107064630 B CN 107064630B CN 201710208460 A CN201710208460 A CN 201710208460A CN 107064630 B CN107064630 B CN 107064630B
- Authority
- CN
- China
- Prior art keywords
- line voltage
- frequency
- voltage
- fundamental phasors
- moment
- 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.)
- Expired - Fee Related
Links
- 238000000691 measurement method Methods 0.000 title claims abstract description 13
- 238000005070 sampling Methods 0.000 claims description 10
- 230000035772 mutation Effects 0.000 claims description 4
- 230000005611 electricity Effects 0.000 claims 1
- 238000000034 method Methods 0.000 description 25
- 238000005259 measurement Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 4
- 238000001914 filtration Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
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/06—Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage by converting frequency into an amplitude of current or voltage
Abstract
本发明涉及一种电力系统频率测量方法及装置，根据电力系统中获取的母线电压采样值，计算依次相差1/P周波的三个时刻的线电压，根据三个时刻的线电压的实部、虚部，计算出系统当前的频率。该测量方法原理简单、容易实现，能够准确测量电力系统的频率，并且不受过零点和谐波的影响。
Description
技术领域
本发明属于电力系统继电保护及控制测量技术领域，具体涉及一种电力系统频率测量方法及装置。
背景技术
频率是电力系统的重要电气参数之一，对电力系统频率测量算法的研究在理论和实际中都有重要的意义。电力系统正常运行时，频率是衡量电能质量的一个重要指标，电力系统出现故障时，频率是判断系统故障的重要依据。
目前，电力系统中的设备都是按照额定频率设计和制造的，只有在额定频率附近运行时，才能发挥最好的性能。系统频率过大的变动，对用户和发电厂的运行都将产生不利影响。随着负荷和发电厂输出功率的变化，频率也将发生相应的变化。因此，频率可以作为预测电力系统是否将失去稳定的判据，从而为切负荷和切机提供可行的判据。
现有的频率估计方法主要分为两类：
1、测量电压一周波的时间，计算频率；
2、利用数字滤波的方法，进行频率估计。
前者受过零点和谐波的影响较大；后者将时域信号转化到频域，寻找最大值，并进行插值，为提高频率估计精度，需要较长的数据窗，插值过程中又需要确定修正系数，增加了装置的计算量和计算时间。
发明内容
本发明的目的是提供一种电力系统频率测量方法及装置，用于解决现有数字滤波方法计算频率时计算工作量大、计算时间长的问题。
为解决上述技术问题，本发明提出一种电力系统频率测量方法，包括以下方法方案：
方法方案一，包括以下步骤：
1)根据电力系统中获取的母线电压采样值，通过傅立叶变换计算依次相差1/P周波的三个时刻的线电压基波相量，分别为当前时刻的线电压基波相量、1/P周波前的线电压基波相量、2/P周波前的线电压基波相量，P≥3；
2)根据所述三个时刻的线电压基波相量的实部、虚部，计算系统当前的频率。
方法方案二，在方法方案一的基础上，所述系统当前频率的计算式如下：
式中，f为系统当前的频率，Ts为采样间隔时间，k为1/P周波，和分别为所述当前时刻的线电压基波相量的虚部和实部，和分别为所述1/P周波前的线电压基波相量的虚部和实部，和分别为所述2/P周波前的线电压基波相量的虚部和实部。
方法方案三，在方法方案一的基础上，计算出所述三个时刻的线电压基波相量后，比较所述三个时刻的线电压基波相量的幅值，当所述三个时刻的线电压基波相量中最小的电压幅值小于设定的第一电压门槛值时，设定系统当前的频率为50赫兹。
方法方案四，在方法方案一的基础上，计算所述系统当前的频率后，还包括输出频率的步骤：当1/P周波内的线电压的最大瞬时值的幅值小于所述第一电压门槛值时，判定线电压发生突然无压，输出历史频率；当所述当前时刻的线电压的基波幅值与一周前线电压的基波幅值之差的绝对值大于设定的第二电压门槛值时，判定线电压发生突然变化，输出历史频率；当线电压没有发生突变时，输出所述系统当前的频率。
为解决上述技术问题，本发明提出一种电力系统频率测量装置，包括以下装置方案：
装置方案一，包括以下单元：
电压计算单元：用于根据电力系统中获取的母线电压采样值，通过傅立叶变换计算依次相差1/P周波的三个时刻的线电压基波相量，分别为当前时刻的线电压基波相量、1/P周波前的线电压基波相量、2/P周波前的线电压基波相量，P≥3；
频率计算单元：用于根据所述三个时刻的线电压的实部、虚部，计算系统当前的频率。
装置方案二，在装置方案一的基础上，所述系统当前频率的计算式如下：
式中，f为系统当前的频率，Ts为采样间隔时间，k为1/P周波，和分别为所述当前时刻的线电压基波相量的虚部和实部，和分别为所述1/P周波前的线电压基波相量的虚部和实部，和分别为所述2/P周波前的线电压基波相量的虚部和实部。
装置方案三，在装置方案一的基础上，还包括比较单元：用于比较所述三个时刻的线电压基波相量的幅值，当所述三个时刻的线电压基波相量中最小的电压幅值小于设定的第一电压门槛值时，设定系统当前的频率为50赫兹。
装置方案四，在装置方案一的基础上，还包括频率输出单元：用于当1/P周波内的线电压的最大瞬时值的幅值小于所述第一电压门槛值时，判定线电压发生突然无压，输出历史频率；当所述当前时刻的线电压的基波幅值与一周前线电压的基波幅值之差的绝对值大于设定的第二电压门槛值时，判定线电压发生突然变化，输出历史频率；当线电压没有发生突变时，输出所述系统当前的频率。
本发明的有益效果是：根据电力系统中获取的母线电压采样值，计算依次相差1/P周波的三个时刻的线电压，根据三个时刻的线电压的实部、虚部，计算出系统当前的频率。该测量方法原理简单、容易实现，能够准确测量电力系统的频率，并且不受过零点和谐波的影响。
附图说明
图1是本发明方法在主接线中的接入位置示意图；
图2是基于线电压的电力系统频率测量方法流程图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步的说明。
本发明的一种电力系统频率测量方法的实施例：
根据电力系统中获取的母线电压采样值，计算依次相差1/P周波的三个时刻的线电压：当前时刻的线电压、1/P周波前的线电压、2/P周波前的线电压：
其中，n为采样序列，N为50Hz频率下一周波的采样点数，和代表不同的相，和分别为相和相的当前时刻线电压的基波相量，和分别为相和相的1/P周波前线电压的基波相量，和分别为相和相的2/P周波前线电压的基波相量，为上述三个时刻线电压的基波相量。
计算当前时刻的频率：
为防止发生PT断线情况，本发明使用线电压进行频率测量，假设线电压为计算信号全周傅氏的实部和虚部，
由于滤波器为偶对称，为奇对称，所以，
式中，|Fc(ω)|和|Fs(ω)|为偶、奇滤波器在ω处的增益，α(ω)为在ω处的相位变化；
利用三角函数变换关系，进行如下变换，
将g2k(ω)和gk(ω)进行比值可得当前系统频率：
其中，Ts为采样间隔时间，为的虚部，为的实部，的虚部，为的实部，为的虚部，为的实部。
上述电力系统的频率计算是线电压在正常水平的情况得出的，当线电压出现以下异常情况时，进行相应的防错处理方法：
(1)比较三个时刻的线电压的幅值，求出和中最小的幅值，即
判断使用的三个时刻的线电压是否存在无压的情况，Uvalve为无压门槛，判据如下：
Umin＜Uvalve1
式中，Umin为上述最小的幅值，Uvalve1为设定的第一电压门槛值，当上述最小的幅值小于设定的第一电压门槛值时，设定系统当前的频率为50Hz。
(2)判别线电压是否发生突然无压和突然变化。突然无压的判据如下：
umax＜Uvalve1
式中，N1为新计算频率下的一周波采样点数，即信号由50Hz频率发生变化后的一周波的采样点数，为线电压的瞬时值，umax为最近周波线电压瞬时值绝对值最大采样点，当其小于第一电压门槛值时，判定线电压发生突然无压。
突然变化的判据如下：
式中，为一周前线电压的基波幅值，为当前时刻线电压的基波幅值，Uvalve2为设定的第二电压门槛值，上式表示若当前时刻的线电压的基波幅值与一周前线电压的基波幅值之差的绝对值大于设定的第二电压门槛值，判定线电压发生突然变化。本发明只有在线电压没有发生突然无压和突然变化时，才按照前述电力系统频率的计算公式进行当前时刻频率的计算；本发明测量电力系统频率的方法是实时采集线电压、并进行实时的频率计算，当线电压发生突然无压或突然变化时，将上次测量的频率结构作为当前时刻测量的频率值。
本发明首先使用DFT的方法计算当前时刻线电压，1/P周波前线电压，2/P周波前线电压，并记录三个时刻线电压的实部、虚部和幅值最小值；然后使用三个时刻线电压的实部和虚部，利用本发明的频率计算方法得到当前时刻的系统频率。本发明的频率测量方法基于全周滤波后的数据，不依赖于过零点，不要求高采样率，运算速度快，能够实时准确的测量各采样时刻的频率。
本发明的一种电力系统频率测量装置的实施例：
包括电压计算单元和频率计算单元。电压计算单元用于根据电力系统中获取的母线电压采样值，计算依次相差1/P周波的三个时刻的线电压：当前时刻的线电压、1/P周波前的线电压、2/P周波前的线电压，P≥3。
频率计算单元用于根据所述三个时刻的线电压的实部、虚部，计算系统当前的频率。
如图1所示的双端系统，将本发明的频率测量装置安装在母线B1时，通过互感器TV1获取母线电压采样值，采样值送入频率测量装置后，进行电力系统频率的测量，测量方法的流程如图2所示。
由于本实施例中所指的电力系统频率测量装置，实际上是基于本发明方法流程的一种计算机解决方案，即一种软件构架，上述装置即为与方法流程相对应的处理进程。由于对上述方法的介绍已经足够清楚完整，故不再详细进行描述。
Claims (6)
1.一种电力系统频率测量方法，其特征在于，包括以下步骤：
1)根据电力系统中获取的母线电压采样值，通过傅立叶变换计算依次相差1/P周波的三个时刻的线电压基波相量，分别为当前时刻的线电压基波相量、1/P周波前的线电压基波相量、2/P周波前的线电压基波相量，P≥3；
2)根据所述三个时刻的线电压基波相量的实部、虚部，计算系统当前的频率；
所述系统当前频率的计算式如下：
式中，f为系统当前的频率，Ts为采样间隔时间，k为1/P周波，和分别为所述当前时刻的线电压基波相量的虚部和实部，和分别为所述1/P周波前的线电压基波相量的虚部和实部，和分别为所述2/P周波前的线电压基波相量的虚部和实部。
2.根据权利要求1所述的电力系统频率测量方法，其特征在于，计算出所述三个时刻的线电压基波相量后，还包括比较步骤：比较所述三个时刻的线电压基波相量的幅值，当所述三个时刻的线电压基波相量中最小的电压幅值小于设定的第一电压门槛值时，设定系统当前的频率为50赫兹。
3.根据权利要求1所述的电力系统频率测量方法，其特征在于，计算所述系统当前的频率后，还包括输出频率的步骤：当1/P周波内的线电压基波相量的最大瞬时值的幅值小于设定的第一电压门槛值时，判定线电压发生突然无压，输出历史频率；当所述当前时刻的线电压的基波幅值与一周前线电压的基波幅值之差的绝对值大于设定的第二电压门槛值时，判定线电压发生突然变化，输出历史频率；当线电压没有发生突变时，输出所述系统当前的频率。
4.一种电力系统频率测量装置，其特征在于，包括以下单元：
电压计算单元：用于根据电力系统中获取的母线电压采样值，通过傅立叶变换计算依次相差1/P周波的三个时刻的线电压基波相量，分别为当前时刻的线电压基波相量、1/P周波前的线电压基波相量、2/P周波前的线电压基波相量，P≥3；
频率计算单元：用于根据所述三个时刻的线电压基波相量的实部、虚部，计算系统当前的频率；
所述系统当前频率的计算式如下：
式中，f为系统当前的频率，Ts为采样间隔时间，k为1/P周波，和分别为所述当前时刻的线电压基波相量的虚部和实部，和分别为所述1/P周波前的线电压基波相量的虚部和实部，和分别为所述2/P周波前的线电压基波相量的虚部和实部。
5.根据权利要求4所述的电力系统频率测量装置，其特征在于，还包括比较单元：用于比较所述三个时刻的线电压基波相量的幅值，当所述三个时刻的线电压基波相量中最小的电压幅值小于设定的第一电压门槛值时，设定系统当前的频率为50赫兹。
6.根据权利要求4所述的电力系统频率测量装置，其特征在于，还包括频率输出单元：用于当1/P周波内的线电压基波相量的最大瞬时值的幅值小于设定的第一电压门槛值时，判定线电压发生突然无压，输出历史频率；当所述当前时刻的线电压的基波幅值与一周前线电压的基波幅值之差的绝对值大于设定的第二电压门槛值时，判定线电压发生突然变化，输出历史频率；当线电压没有发生突变时，输出所述系统当前的频率。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710208460.7A CN107064630B (zh) | 2017-03-31 | 2017-03-31 | 一种电力系统频率测量方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710208460.7A CN107064630B (zh) | 2017-03-31 | 2017-03-31 | 一种电力系统频率测量方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107064630A CN107064630A (zh) | 2017-08-18 |
| CN107064630B true CN107064630B (zh) | 2019-11-12 |
Family
ID=59602110
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710208460.7A Expired - Fee Related CN107064630B (zh) | 2017-03-31 | 2017-03-31 | 一种电力系统频率测量方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107064630B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108169558A (zh) * | 2017-12-28 | 2018-06-15 | 长园深瑞继保自动化有限公司 | 电力系统实时测频方法 |
| CN109725197B (zh) * | 2018-12-25 | 2021-05-11 | 西电通用电气自动化有限公司 | 一种多处理单元继电保护系统频率跟踪方法及系统 |
| CN110824247A (zh) * | 2019-10-31 | 2020-02-21 | 许昌许继软件技术有限公司 | 一种电力系统频率测量方法、母线电压校正方法及装置 |
| CN112098719B (zh) * | 2020-08-06 | 2023-07-14 | 许昌许继软件技术有限公司 | 一种频率变化的正弦波信号傅氏计算方法及继电保护装置 |
| CN112051444B (zh) * | 2020-08-18 | 2023-02-10 | 许继集团有限公司 | 一种基于抛物线插值迭代的抗谐波频率校准方法及装置 |
| CN112905958B (zh) * | 2021-01-27 | 2024-04-19 | 南京国电南自电网自动化有限公司 | 基于测控装置的短时数据窗遥测数据状态辨识方法及系统 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH03128465A (ja) * | 1989-07-24 | 1991-05-31 | Sony Corp | 周波数検出回路 |
| CN102305889A (zh) * | 2011-07-18 | 2012-01-04 | 山东大学 | 双模式220v用户侧电网相量测量及发送装置及方法 |
| CN102495282A (zh) * | 2011-12-27 | 2012-06-13 | 国电南瑞科技股份有限公司 | 一种交流电压频率的测量方法 |
| CN104049144A (zh) * | 2014-06-27 | 2014-09-17 | 福州大学 | 一种滤除衰减直流分量的同步相量测量的实现方法 |
| CN104076201A (zh) * | 2014-07-15 | 2014-10-01 | 国家电网公司 | 不受非周期噪声影响的变电站母线电压瞬态频率测量方法 |
| CN105372492A (zh) * | 2014-08-31 | 2016-03-02 | 盛吉高科（北京）科技有限公司 | 基于三条dft复数谱线的信号频率测量方法 |
| CN106405230A (zh) * | 2016-08-23 | 2017-02-15 | 北京德威特继保自动化科技股份有限公司 | 频率测量方法和装置 |
-
2017
- 2017-03-31 CN CN201710208460.7A patent/CN107064630B/zh not_active Expired - Fee Related
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH03128465A (ja) * | 1989-07-24 | 1991-05-31 | Sony Corp | 周波数検出回路 |
| CN102305889A (zh) * | 2011-07-18 | 2012-01-04 | 山东大学 | 双模式220v用户侧电网相量测量及发送装置及方法 |
| CN102495282A (zh) * | 2011-12-27 | 2012-06-13 | 国电南瑞科技股份有限公司 | 一种交流电压频率的测量方法 |
| CN104049144A (zh) * | 2014-06-27 | 2014-09-17 | 福州大学 | 一种滤除衰减直流分量的同步相量测量的实现方法 |
| CN104076201A (zh) * | 2014-07-15 | 2014-10-01 | 国家电网公司 | 不受非周期噪声影响的变电站母线电压瞬态频率测量方法 |
| CN105372492A (zh) * | 2014-08-31 | 2016-03-02 | 盛吉高科（北京）科技有限公司 | 基于三条dft复数谱线的信号频率测量方法 |
| CN106405230A (zh) * | 2016-08-23 | 2017-02-15 | 北京德威特继保自动化科技股份有限公司 | 频率测量方法和装置 |
Non-Patent Citations (2)
| Title |
|---|
| 一种基于傅式变换算法的系统频率变化求取方法;王拔群等;《电力系统保护与控制》;20091201;第37卷(第23期);第54-57页 * |
| 基于基相量循环移位的基波相位和频率高精度快速检测方法;陈明凯等;《电工电能新技术》;20070131;第26卷(第1期);第38-42页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107064630A (zh) | 2017-08-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107064630B (zh) | 一种电力系统频率测量方法及装置 | |
| CN106199183B (zh) | 一种实现次同步振荡在线辨识告警的pmu和方法 | |
| CN106841778B (zh) | 基于pmu实现的次同步和超同步谐波参数的处理方法 | |
| CN107345984B (zh) | 一种基于信号识别的自适应同步相量测量方法 | |
| CN109283429B (zh) | 一种基于正负序阻抗相等原理的配电网故障测距方法 | |
| WO2009048964A1 (en) | Real-time power system oscillation detection using modal analysis | |
| JP2021148778A (ja) | 非接地配電システムのための過渡状態に基づいた故障位置特定方法 | |
| CN105372616B (zh) | 一种cvt状态在线测试系统及在线评估方法 | |
| CN104242267A (zh) | 一种风力发电送出输电线路距离保护方法 | |
| CN109284933B (zh) | 一种基于数理统计的电子式互感器状态评估系统及方法 | |
| CN102338827A (zh) | 用于电网波形畸变分析和电力谐波参数自动监测方法 | |
| CN102253282A (zh) | 纳托尔窗函数连续频谱内插电力谐波参数获取方法 | |
| CN102621388A (zh) | 基于同步时域信号的输电线路集中参数在线确定方法 | |
| CN103066591A (zh) | 一种基于实时测量的电网参数偏差识别方法 | |
| CN102798758B (zh) | 并联电容器组串联电抗率的测量方法和系统 | |
| JP5618758B2 (ja) | 電力系統の短絡容量監視方法およびそのシステム | |
| CN108414838A (zh) | 一种逆变器并联系统线路阻抗测量方法 | |
| CN108896852B (zh) | 一种公共接入点短路容量在线测量方法及系统 | |
| CN107179476B (zh) | 一种配网故障测距方法 | |
| CN102004191A (zh) | 消弧线圈并列运行系统中计算线路电容的方法 | |
| KR101664010B1 (ko) | 시각 동기화 데이터를 이용한 전력계통의 선로정수 추정방법 | |
| CN111638401A (zh) | 一种电容器在线监测系统及方法 | |
| JP5971692B2 (ja) | 電力系統の短絡容量監視方法およびそのシステム | |
| CN106156489B (zh) | 一种电力系统同步相量快速计算方法 | |
| CN109283385A (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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20191112 |