CN111239489A - 一种结合prony与fft算法的次同步振荡分析方法 - Google Patents
一种结合prony与fft算法的次同步振荡分析方法 Download PDFInfo
- Publication number
- CN111239489A CN111239489A CN201811441047.6A CN201811441047A CN111239489A CN 111239489 A CN111239489 A CN 111239489A CN 201811441047 A CN201811441047 A CN 201811441047A CN 111239489 A CN111239489 A CN 111239489A
- Authority
- CN
- China
- Prior art keywords
- oscillation
- subsynchronous oscillation
- analysis
- prony
- subsynchronous
- 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/16—Spectrum analysis; Fourier analysis
Landscapes
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开一种结合PRONY与FFT算法的次同步振荡分析方法,包括如下步骤:(1)从次同步振荡监测主站获取暂态录波数据,根据三相瞬时电压与瞬时电流计算出监测设备的瞬时功率;(2)使用PRONY算法对瞬时功率数据进行振荡模式分析,定位出振荡的起止时间,得到的次同步振荡特征;(3)根据步骤(2)中定位出振荡的起止时间,对此段时间内的瞬时电流进行FFT,精确计算出次同步振荡时简谐波电流的频率、幅值、分量等信息;(4)将分析结果进行综合统计,生成事件简报,并将结果展示。此种方法可用于对大量暂态录波数据进行次同步振荡分析,提高了次同步振荡自动分析的效率,实现了次同步振荡的快速、准确辨识。
Description
技术领域
本发明属于电力系统自动化领域,特别涉及一种对大量暂态录波数据进行次同步振荡分析的方法。
背景技术
近年来,随着电力系统规模的不断扩大,大规模新能源电厂并网运行,以及大量电力电子设备的应用,使得电力系统运行面临的新挑战。由风电场诱发次同步振荡现象在新疆哈密地区多次发生,波及到了周边的火电机组,对电网安全稳定运行产生了不良影响。
次同步振荡是一类严重的系统稳定性问题,不但会使系统产生振荡现象,而且极易造成汽轮发电机组的大轴损毁,目前次同步振荡的原因和机理、监测和控制仍在研究中。当前,次同步振荡分析一般采传统的信号分析方法,分析手段较少,效率不高,无法满足对大量暂态录波数据的分析需求。
传统的信号分析方法具有局限性,比如FFT算法只能给出信号的频率分布情况,不能给出频率随时间变化的情况,对于非平稳信号分析缺乏对时间的定位功能。
本发明提出一种结合PRONY与FFT算法的次同步振荡分析方法,弥补了以上不足,满足对大量暂态录波数据的分析需求,提高了次同步振荡自动分析的效率,实现了次同步振荡的快速、准确辨识。
发明内容
本发明的主要目的,在于提供一种结合PRONY与FFT算法的次同步振荡分析方法,其提高了次同步振荡自动分析的效率,实现了次同步振荡的快速、准确辨识。
为了达成上述目的,本发明的解决方案是:
一种结合PRONY与FFT算法的次同步振荡分析方法,包括如下步骤:
步骤(1)、从次同步振荡监测主站获取暂态录波数据,根据三相瞬时电压与瞬时电流计算出监测设备的瞬时功率;
步骤(2)、使用PRONY算法对瞬时功率数据进行振荡模式分析,定位出振荡的起止时间,得到次同步振荡中功率振荡的特征量;
步骤(3)、根据步骤(2)中定位出的振荡起止时间,对此段时间内的瞬时电流进行FFT,精确计算出次同步振荡时简谐波电流的频率、幅值、分量等信息。
进一步地,上述方法还包括
步骤(4)、将分析结果进行综合统计,生成事件简报,并将结果展示。
进一步地,所述步骤(2)具体包括:
步骤(21)、读入暂态录波数据,以固定时间窗为分析单位,形成连续的时间段;
步骤(22)、使用PRONY算法对固定时间窗内的瞬时功率进行振荡模式分析,定位出振荡的开始时间T_start,结束时间T_end,得到次同步振荡频率F_p、幅值Amp_p。
所述步骤(3)具体包括:
根据步骤(2)中定位出的振荡起止时间,在振荡的时间内即开始时间T_start到结束时间T_end,对监测设备的瞬时电流进行FFT运算,精确计算出次同步振荡时简谐波电流的频率F_i、幅值Amp_i、分量Cpn_i等信息。
进一步地,所述步骤(3)还包括:
由于振荡功率频率与简谐波电流频率存在互补的关系,结合PRONY计算出的振荡频率F_p与FFT计算出的简谐波电流的频率F_i,综合辨识出次同步振荡的特性。
PRONY算法是用一组指数项的线性组合来拟和等间距采样数据的方法,可以从数据中分析出信号的幅值、相位、阻尼因子、频率等信息。其最大的优点是既可以对仿真结果进行分析,又可以对实时测量数据进行分析。电网发生次同步振荡时,功率振荡模式明显,使用PRONY进行分析比较适合。对于振荡时的电流、电压的简谐波分量,偏向于研究其频率特性,使用快速傅里叶变换算法(FFT)比较合适,快速傅里叶变换算法(FFT)可以分析信号在各个频率成分中的分布情况,从而分析出简谐波的频率和分量。
本发明的有益效果为:采用PRONY与FFT算法结合的分析方法,可用于对大量暂态录波数据进行次同步振荡分析,提高了次同步振荡自动分析的效率,实现了次同步振荡的快速、准确辨识。
附图说明
图1是本发明的方法流程图。
具体实施方式
以下将结合附图及具体实施例,对本发明的实现过程进行详细说明。
参考图1所示,本发明提供一种结合PRONY与FFT算法的次同步振荡分析方法,包括如下步骤:
步骤(1)、从次同步振荡监测主站获取暂态录波数据,根据三相瞬时电压与瞬时电流计算出监测设备的瞬时功率;
步骤(2)、使用PRONY算法对瞬时功率数据进行振荡模式分析,定位出振荡的起止时间,得到的次同步振荡中功率振荡的特征量;
具体包括:
步骤(21)、读入暂态录波数据,以固定时间窗为分析单位,形成连续的时间段;
步骤(22)、使用PRONY算法对固定时间窗内的瞬时功率进行振荡模式分析,定位出振荡的开始时间T_start,结束时间T_end,得到次同步振荡频率F_p、幅值Amp_p。
步骤(3)、根据步骤(2)中定位出振荡的起止时间,对此段时间内的瞬时电流进行FFT,精确计算出次同步振荡时简谐波电流的频率、幅值、分量等信息;
具体包括
(31)在振荡的时间内(开始时间T_start,结束时间T_end)对监测设备的瞬时电流进行FFT,精确计算出次同步振荡时简谐波电流的频率F_i、幅值Amp_i、分量Cpn_i等信息;
(32)由于振荡功率频率与简谐波电流频率存在互补的关系,结合PRONY计算出的振荡频率F_p与FFT计算出的简谐波电流的频率F_i,综合辨识出次同步振荡的特性。
步骤(4)、将分析结果进行综合统计,生成事件简报,并将结果展示。
以上实施例仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明保护范围之内。
Claims (4)
1.一种结合PRONY与FFT算法的次同步振荡分析方法,其特征在于,包含如下步骤:
步骤(1)、从次同步振荡监测主站获取暂态录波数据,根据三相瞬时电压与瞬时电流计算出监测设备的瞬时功率;
步骤(2)、使用PRONY算法对瞬时功率数据进行振荡模式分析,定位出振荡的起止时间,得到次同步振荡特征;
步骤(3)、根据步骤(2)中定位出的振荡起止时间,对此段时间内的瞬时电流进行FFT,精确计算出次同步振荡时简谐波电流的频率、幅值、分量信息。
2.如所述权利1所述结合PRONY与FFT算法的次同步振荡分析方法,其特征在于,还包括
步骤(4)、将分析结果进行综合统计,生成事件简报,并将结果展示。
3.如所述权利1所述结合PRONY与FFT算法的次同步振荡分析方法,其特征在于,所述步骤(2)具体包括:
步骤(21)、读入暂态录波数据,以固定时间窗为分析单位,形成连续的时间段;
步骤(22)、使用PRONY算法对固定时间窗内的瞬时功率进行振荡模式分析,定位出振荡的开始时间T_start,结束时间T_end,得到次同步振荡频率F_p、幅值Amp_p;
所述步骤(3)具体包括:
根据步骤(2)中定位出的振荡起止时间,在振荡的时间内即开始时间T_start到结束时间T_end,对监测设备的瞬时电流进行FFT运算,计算出次同步振荡时简谐波电流的频率F_i、幅值Amp_i、分量Cpn_i信息。
4.如所述权利1所述结合PRONY与FFT算法的次同步振荡分析方法,其特征在于,所述步骤(3)还包括:
由于振荡功率频率与简谐波电流频率存在互补的关系,结合PRONY计算出的振荡频率与FFT计算出的简谐波电流的频率,相互验证,相互补充,综合辨识出次同步振荡的特性。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811441047.6A CN111239489B (zh) | 2018-11-29 | 2018-11-29 | 一种结合prony与fft算法的次同步振荡分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811441047.6A CN111239489B (zh) | 2018-11-29 | 2018-11-29 | 一种结合prony与fft算法的次同步振荡分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111239489A true CN111239489A (zh) | 2020-06-05 |
CN111239489B CN111239489B (zh) | 2022-02-18 |
Family
ID=70879233
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811441047.6A Active CN111239489B (zh) | 2018-11-29 | 2018-11-29 | 一种结合prony与fft算法的次同步振荡分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111239489B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112688321A (zh) * | 2021-01-11 | 2021-04-20 | 河海大学 | 大电网次/超同步振荡路径获取方法及系统 |
CN114301055A (zh) * | 2022-02-17 | 2022-04-08 | 河海大学 | 基于宽频量测的电力系统间谐波潮流获取方法与系统 |
CN115083133A (zh) * | 2022-06-16 | 2022-09-20 | 西安热工研究院有限公司 | 一种立式水轮发电机组振荡和失步检测预警方法及系统 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1093171A (zh) * | 1994-01-31 | 1994-10-05 | 清华大学 | 监测电力系统次同步振荡的方法及监测仪 |
EP0683398A1 (en) * | 1993-02-02 | 1995-11-22 | HIRATA, Yoshimutsu | Non-harmonic analysis of waveform data and synthesizing processing system |
CN102122823A (zh) * | 2011-02-23 | 2011-07-13 | 南京南瑞继保电气有限公司 | 电力系统振荡扰动源定位的方法 |
CN104635109A (zh) * | 2015-01-23 | 2015-05-20 | 云南电网有限责任公司电力科学研究院 | 一种从故障录波数据中直接计算提取电网暂态非周期分量幅值和时间常数的高精度方法 |
CN104852392A (zh) * | 2015-04-24 | 2015-08-19 | 神华国华(北京)电力研究院有限公司 | 一种基于Prony算法的次同步振荡模态衰减系数计算方法 |
CN105652087A (zh) * | 2016-01-11 | 2016-06-08 | 南京南瑞继保电气有限公司 | 带间谐波分析和连续采样值录波功能的同步相量测量装置 |
CN107966611A (zh) * | 2017-11-24 | 2018-04-27 | 广东电网有限责任公司电力调度控制中心 | 一种基于矢量匹配的电力系统次超同步谐波检测方法 |
CN108490254A (zh) * | 2018-01-29 | 2018-09-04 | 国电南瑞科技股份有限公司 | 一种基于多算法融合的宽带多频信号测量方法 |
CN108845230A (zh) * | 2018-06-22 | 2018-11-20 | 国网陕西省电力公司电力科学研究院 | 一种次同步振荡随机时变模态辨识方法 |
-
2018
- 2018-11-29 CN CN201811441047.6A patent/CN111239489B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0683398A1 (en) * | 1993-02-02 | 1995-11-22 | HIRATA, Yoshimutsu | Non-harmonic analysis of waveform data and synthesizing processing system |
CN1093171A (zh) * | 1994-01-31 | 1994-10-05 | 清华大学 | 监测电力系统次同步振荡的方法及监测仪 |
CN102122823A (zh) * | 2011-02-23 | 2011-07-13 | 南京南瑞继保电气有限公司 | 电力系统振荡扰动源定位的方法 |
CN104635109A (zh) * | 2015-01-23 | 2015-05-20 | 云南电网有限责任公司电力科学研究院 | 一种从故障录波数据中直接计算提取电网暂态非周期分量幅值和时间常数的高精度方法 |
CN104852392A (zh) * | 2015-04-24 | 2015-08-19 | 神华国华(北京)电力研究院有限公司 | 一种基于Prony算法的次同步振荡模态衰减系数计算方法 |
CN105652087A (zh) * | 2016-01-11 | 2016-06-08 | 南京南瑞继保电气有限公司 | 带间谐波分析和连续采样值录波功能的同步相量测量装置 |
CN107966611A (zh) * | 2017-11-24 | 2018-04-27 | 广东电网有限责任公司电力调度控制中心 | 一种基于矢量匹配的电力系统次超同步谐波检测方法 |
CN108490254A (zh) * | 2018-01-29 | 2018-09-04 | 国电南瑞科技股份有限公司 | 一种基于多算法融合的宽带多频信号测量方法 |
CN108845230A (zh) * | 2018-06-22 | 2018-11-20 | 国网陕西省电力公司电力科学研究院 | 一种次同步振荡随机时变模态辨识方法 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112688321A (zh) * | 2021-01-11 | 2021-04-20 | 河海大学 | 大电网次/超同步振荡路径获取方法及系统 |
CN112688321B (zh) * | 2021-01-11 | 2023-05-23 | 河海大学 | 大电网次/超同步振荡路径获取方法及系统 |
CN114301055A (zh) * | 2022-02-17 | 2022-04-08 | 河海大学 | 基于宽频量测的电力系统间谐波潮流获取方法与系统 |
CN115083133A (zh) * | 2022-06-16 | 2022-09-20 | 西安热工研究院有限公司 | 一种立式水轮发电机组振荡和失步检测预警方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN111239489B (zh) | 2022-02-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111239489B (zh) | 一种结合prony与fft算法的次同步振荡分析方法 | |
Moreno et al. | Propagation of flicker in electric power networks due to wind energy conversions systems | |
US9903896B2 (en) | Method for testing dynamic model parameters of wind power plant | |
CN108196187B (zh) | 双馈风力发电机转子绕组三相不对称负载故障诊断方法 | |
CN111082455A (zh) | 一种风电场并网性能监测方法及系统 | |
CN108252873A (zh) | 一种风力发电机组在线数据监测及其性能评估的系统 | |
Ma et al. | Evaluation and design tools for the reliability of wind power converter system | |
CN103701389B (zh) | 一种双馈感应发电机控制器参数辨识方法 | |
CN114865653A (zh) | 电网振荡源分布式定位方法、系统、设备及存储介质 | |
CN109061297A (zh) | 一种次/超同步振荡在线监测方法和装置 | |
CN112018784B (zh) | 一种基于同步相量测量数据的次同步谐振溯源方法 | |
CN113219278A (zh) | 新能源场站一次调频测试系统及方法 | |
CN110879321A (zh) | 一种基于实际电网扰动波形的风电场avc子站性能测试方法 | |
CN115864449B (zh) | 新能源接入系统振荡监测方法、装置、系统、设备及介质 | |
Wang et al. | Reliability assessment of wind farm active power based on sequential monte-carlo method | |
Zang et al. | Variable parameter Kalman filter based dynamic harmonic state estimation for power systems with wind energy integration | |
CN111413578B (zh) | 一种次同步振荡的实时监测预警方法 | |
CN104201712B (zh) | 一种基于空间平均风速的风电功率实时预测计算方法 | |
Kumar et al. | Improvements in Synchronous Generator Parameter Tuning Using PMU Data | |
Lv et al. | Research on detection method of power quality for wind power connected to grid | |
CN110865259A (zh) | 一种风电场电能质量评估方法和装置 | |
CN203981778U (zh) | 一种基于dsp的电力参数测量装置 | |
Fan et al. | Research about the harmonic/inter-harmonic characteristic of grid connected with wind plant based on wide-frequency measurement data | |
Wanzheng et al. | Simulation evaluation of fast frequency response capacity of new energy power stations | |
Yang et al. | An improved aggregated model of wind farm with fixed-speed wind turbines |
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 |