CN107179440A - 一种在线自适应频率变化的次同步振荡间谐波提取方法 - Google Patents
一种在线自适应频率变化的次同步振荡间谐波提取方法 Download PDFInfo
- Publication number
- CN107179440A CN107179440A CN201710338027.5A CN201710338027A CN107179440A CN 107179440 A CN107179440 A CN 107179440A CN 201710338027 A CN201710338027 A CN 201710338027A CN 107179440 A CN107179440 A CN 107179440A
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- frequency
- bandpass filter
- oscillation
- 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
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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Stabilization Of Oscillater, Synchronisation, Frequency Synthesizers (AREA)
Abstract
本发明公开了一种在线自适应频率变化的次同步振荡间谐波提取方法。该方法首先通过在线实时获取线路的相电流、相电压原始采样数据,或者通过相电流、相电压原始采样数据计算线路瞬时功率,从而得到待计算提取的输入数据,将输入数据经过窄带的梳状滤波器组进行滤波,获取每一个窄带滤波器的实时输出数据,对每一个窄带滤波器的输出实时波形峰值进行比较,得到宽频带幅值的极值点,极值点对应的滤波器的中心频率即为振荡的主要振荡频率,通过对该滤波器的输出数据通过过零比较数据拟合,确定精确的主要振荡频率。该算法能够应用于次同步振荡监测在线监测装置上,实现振荡信号的在线快速辨识。
Description
技术领域
本发明属于电力系统动态监测领域,具体涉及一种在线提取次同步、超同步间谐波的方法。
背景技术
近年来,随着特高压电网快速发展,风电、光伏等新能源大量并网,远距离跨区输电规模持续增长,电网电力电子化特征凸显,电网特性发生深刻变化。各种电力电子设备或装置,如静止同步补偿器(STATCOM)、柔性直流输电换流器得到了广泛的应用。然而,一旦这些新型电力设备接入电网运行,在特定的控制模式或运行方式下,易形成持续次同步和/或超同步谐波注入系统,激发近区机组轴系扭振,造成火电机组跳闸、直流闭锁等连锁反应,危及电网的安全稳定运行和影响供电质量。
随着电力电子化系统宽频带稳定性问题日益突出,对系统进行实时在线的宽频监测非常必要,目前在电网上应用的监测手段,主要是针对工频信号或者工频的整次谐波进行的,在这些频点上,常规的监测手段非常有效。而对于非整次谐波间谐波的监测与分析手段,相对来说还不完整,对于一个信号的分析,主要是分析信号两个域空间的信号特征:时域和频域特征,离线傅里叶变换在数据窗足够长的情况下可以有效的对信号进行频域辨识,其频率分辨率和数据窗的选择息息相关,因此为了获取比较准确的频率信息,往往采用10s或者更长的数据窗进行分析,响应时间太长。在线的信号提取方法有无限冲击响应IIR滤波、有限冲击响应FIR滤波等等数字滤波技术,前提条件是频率信息必须已知并且信号频率要单一,否则频率混叠问题无法解决。
本发明在综合无限冲击响应IIR滤波的特点以及频率变化适应的需求,设计了一种由窄带带通滤波器组成的梳状滤波滤波器组,保证每一个带通滤波器输出频率的唯一性,通过与过零监测、曲线拟合方法相结合适应,从而在线适应频率变化,能够较准确的获取间谐波信号的幅值、频率的特征信息,此方法可以作为在线监测或者保护装置来使用。
发明内容
本发明的目的在于提出一种在线自适应频率变化的次同步振荡间谐波提取方法,该方法可以用来在线监测电网中的次同步和/或超同步谐波电压、谐波电流与谐波功率的特征并能够适应谐波频率的变化。本发明采用以下技术方案:
一种在线自适应频率变化的次同步振荡间谐波提取方法,其特征在于,所述谐波提取方法包括以下步骤:
(1)以至少1kHz的采样率采集某一条线路的三相相电压ua、ub、uc,三相相电流ia、ib、ic,或者根据三相相电压、三相相电流计算线路的瞬时功率p;
(2)设置梳状带通滤波器组,在要分析的频带内上实现频率细分,带通滤波器组的通带宽度为Δfp;
(3)把三相相电压ua、ub、uc,三相相电流ia、ib、ic或者线路瞬时功率p分别通过梳状带通滤波器组,获取各个带通滤波器组输出三相电压upai、upbi、upci,三相相电流ipai、ipbi、ipci或者瞬时功率ppi,其中i表示梳状滤波器组中各个带通滤波器的序号;
(4)求取各个带通滤波器输出upai、upbi、upci,ipai、ipbi、ipci或者ppi在同一时间窗的幅值序列[Mu1,Mu2,...,Mun]、[Mi1,Mi2,...,Min]或者[Mp1,Mp2,...,Mpn],其中:n是否为带通滤波器的个数;在幅值序列[Mu1,Mu2,...,Mun]、[Mi1,Mi2,...,Min]或者[Mp1,Mp2,...,Mpn]中求取该幅值序列中的极值点,形成极值点序列[Msso1,Msso2,...,Mssom],从而根据极值点对应的带通滤波器中心频率得到各个振荡模式的初选频率序列[fsso1,fsso2,...,fssom],其中:m为同一时间窗幅值序列中的极值点个数。
(5)对初选频率序列[fsso1,fsso2,...,fssom]对应带通滤波器输出的实时数据选取大于一个周波的数据进行过零检测,得到复选后的振荡频率[f’sso1,f’sso2,...,f’ssom];
(6)得到的幅值极值点[Msso1,Msso2,...,Mssom]和复选后的振荡频率[f’sso1,f’sso2,...,f’ssom]数据即为提取出来的次同步振荡间谐波特征量。
本发明进一步包括以下优选方案:
在步骤(2)中梳状带通滤波器组需要同时满足以下三个条件:
2.1带通滤波器组的通带宽度Δfp≤1Hz,同时带通滤波器组通带范围为5~45Hz和55~95Hz;
2.2梳状带通滤波器组各带通滤波器的通带范围无缝对接;
其中,n为梳状滤波器组包含的带通滤波器个数;
2.3各带通滤波器对工频信号的衰减要在40dB以上。
在步骤(4)中,极值点的判断是每个带通滤波器输出幅值和相邻的两个带通滤波器输出幅值做比较,以线路相电流为例,假设第i个带通滤波器输出幅值为Mii,该带通滤波器的中心频率为f0i,则第i-1和i+2个滤波器的输出幅值为Mi(i-1)和Mi(i+1),则:
Mii、Mi(i-1)、Mi(i+1)分别等于一时刻相电流第i、i-1、i+1个带通滤波器输出选取的数据窗的最大值的绝对值;
如果Mii≥Mi(i-1)并且Mii≥Mi(i+1),则判断Mii是极值点,记作Mssoi;
振荡频率的初选频率fsso=f0i;
其中:fs为采样频率。
对于电压的幅值序列[Mu1,Mu2,...,Mun]或功率的幅值序列[Mp1,Mp2,...,Mpn],其极值点的判断方式和电流的幅值序列相同。
在步骤(5)中,对极值点对应的带通滤波器实时波形数据采用过零监测+数据拟合的方法得到复选后的振荡频率频率,以线路电流信号为例,采样间隔时间第一个上升沿过零点之前的采样点序号为nl,对应的采样值为vl;第一个上升沿过零点之后的采样点序号为nl+1,对应的采样值为vl+1;第二个上升沿过零点之前的采样点序号为nm,对应的采样值为vm;第二个上升沿过零点之前的采样点序号为nm+1,对应的采样值为vm+1;采取一阶线性拟合的方式或者高阶多项式拟合的方式对电流信号进行拟合,获取过零点时刻;
以一阶线性拟合的方式为例,过零点到nl+1的时间ΔT1为:
同理,第二个过零点到nl+1的时间ΔT2为:
复选后的振荡频率[f’sso1,f’sso2,...,f’ssom],第i个振荡频率的表达式为:
本发明具有以下有益的技术效果:
根据本发明在线自适应频率变化的次同步振荡间谐波提取方法,可以在监测或者保护装置上实现在线提取振荡的相关信息特征,从而为监测预警及保护判据的实施提供相关依据。
附图说明
图1是本发明在线自适应频率变化的次同步振荡间谐波提取方法的流程示意图;
图2是梳状滤波器组的幅频特性图;
图3是梳状滤波器组输出极值点选择;
图4是过零比较与线性插值示意图。
具体实施方式
下面结合说明书附图对本发明的技术方案的具体实施作进一步详细说明。
本发明提出的方法在嵌入式开发的装置上实现,用来对次同步振荡信号进行在线监测。本发明针对次同步振荡间谐波的提取流程如图1所示。
步骤1:待测线路三相相电压、三相相电流或者三相瞬时功率的采集
在本实现中,装置以1200Hz采样率采集某一条线路的三相相电压ua、ub、uc,三相相电流ia、ib、ic,以及根据三相相电压、三相相电流计算线路的瞬时功率p,计算瞬时功率的公式为:
p=uaia+ubib+ucic
步骤2:梳状滤波器组的设计
(1)在本实现中,为了解决多个振荡模式同时存在情况下的带通滤波器输出信号频率唯一性问题,取Δfp=1Hz,同时滤波器组通带范围为5~45Hz,55~95Hz,兼顾次同步频率和超同步频率信号,所以梳状滤波器组一共由80个窄带带通滤波器组成,窄带带通滤波器选择巴特沃斯4阶带通滤波器来实现,80个窄带带通滤波器的幅频特性如图2所示,通过这些滤波器的设计,可以尽可能实现每个滤波器输出信号频率信号的唯一性。
(2)各个带通滤波器组的通带范围要实现无缝对接;
(3)带通滤波器的设计要实现对工频信号的足够衰减,从而提高次同步/超同步信号的辨识能力,通常带通滤波器对工频信号的衰减要40dB以上。
步骤3:梳状滤波器组的计算输出
每个带通滤波器的离散化格式为:
其中N(z),D(z)为算子z的多项式,即:
N(z)=bmzm+bm-1z(m-1)+...+b1z+b0;
D(z)=amzm+am-lz(m-l)+...+alz+a0;
以中心频率为5.5Hz为例,通带宽度为1Hz,m统一取4,则其存储系数为:
[10.5811141289120e-003,-71.8451679981765e-003,199.556371100672e-003,
-271.311026351323e-003,127.721330104613e-003,
-271.423985741603e-003,191.059425476604e-003,-66.4348295925944e-003,
9.47571062195649e-003,8.45678532126587e-003,];
前5个系数为分子系数,后5个系数为分母系数,都按照m的降幂排列。
步骤4:初选频率的自适应选择过程
极值点的判断主要是每个带通滤波器输出幅值和相邻的两个滤波器输出幅值做比较,如图3所示,以电流信号为例,假设滤波器输入信号含有19.6Hz信号的次同步谐波,采样频率fs=1200Hz,假设第i个带通滤波器输出幅值为Mii=1,带通滤波器的中心频率为f0i=20Hz,Δfp=1Hz,第i-1和i+2个滤波器的输出幅值为Mi(i-1)=0.95和Mi(i+1)=0.90,其中:
Mii、Mi(i-1)、Mi(i+1)分别等于一时刻第i、i-1、i+1个带通滤波器输出选取的数据窗的最大值的绝对值;
因为Mii≥Mi(i-1)并且Mii≥Mi(i+1),则判断Mii是极值点,则初选频率fsso=f0i=20Hz为振荡频率之一。
步骤5:精确频率的获取
过零比较及插值处理过程如图4所示,在本实现中,以20Hz滤波器输出为初选信号,经过滤波器后信号的幅值为1,则过零点到nl+1的时间ΔT1为:
同理,第二个过零点到nl+1的时间ΔT2为:
信号的准确频率为:
从计算的结果来看,最终的计算出来的频率和输入的次同步谐波信号相差很小。
步骤6:把计算出来的频率信息、幅值信息以及20Hz为中心频率的滤波器输出数据可以作为后续监测和保护应用功能使用。
本发明申请人结合说明书附图对本发明的实施示例做了详细的说明与描述,但是本领域技术人员应该理解,以上实施示例仅为本发明的优选实施方案,详尽的说明只是为了帮助读者更好地理解本发明精神,而并非对本发明保护范围的限制,相反,任何基于本发明的发明精神所作的任何改进或修饰都应当落在本发明的保护范围之内。
Claims (4)
1.一种在线自适应频率变化的次同步振荡间谐波提取方法,其特征在于,所述谐波提取方法包括以下步骤:
(1)以至少1kHz的采样率采集某一条线路的三相相电压ua、ub、uc,三相相电流ia、ib、ic,或者根据三相相电压、三相相电流计算线路的瞬时功率p;
(2)设置梳状带通滤波器组,在要分析的频带内上实现频率细分,带通滤波器组的通带宽度为Δfp;
(3)把三相相电压ua、ub、uc,三相相电流ia、ib、ic或者线路瞬时功率p分别通过梳状带通滤波器组,获取各个带通滤波器组输出三相电压upai、upbi、upci,三相相电流ipai、ipbi、ipci或者瞬时功率ppi,其中i表示梳状滤波器组中各个带通滤波器的序号;
(4)求取各个带通滤波器输出upai、upbi、upci,ipai、ipbi、ipci或者ppi在同一时间窗的幅值序列[Mu1,Mu2,...,Mun]、[Mi1,Mi2,...,Min]或者[Mp1,Mp2,...,Mpn],其中:n为带通滤波器的个数;在幅值序列[Mu1,Mu2,...,Mun]、[Mi1,Mi2,...,Min]或者[Mp1,Mp2,...,Mpn]中求取该幅值序列中的极值点,形成极值点序列[Msso1,Msso2,...,Mssom],从而根据极值点对应的带通滤波器中心频率得到各个振荡模式的初选频率序列[fsso1,fsso2,...,fssom],其中:m为同一时间窗幅值序列中的极值点个数;
(5)对初选频率序列[fsso1,fsso2,...,fssom]对应带通滤波器输出的实时数据选取大于一个周波的数据进行过零检测,得到复选后的振荡频率[f’sso1,f’sso2,...,f’ssom];
(6)得到的幅值极值点[Msso1,Msso2,...,Mssom]和复选后的振荡频率[f’sso1,f’sso2,...,f’ssom]数据即为提取出来的次同步振荡间谐波特征量。
2.根据权利要求1所述的在线自适应频率变化的次同步振荡间谐波提取方法,其特征在于:
在步骤(2)中梳状带通滤波器组需要同时满足以下三个条件:
2.1带通滤波器组的通带宽度Δfp≤1Hz,同时带通滤波器组通带范围为5~45Hz和65~95Hz;
2.2梳状带通滤波器组各带通滤波器的通带范围无缝对接;
其中,n为梳状滤波器组包含的带通滤波器个数;
2.3各带通滤波器对工频信号的衰减要在40dB以上。
3.根据权利要求1所述的在线自适应频率变化的次同步振荡间谐波提取方法,其特征在于:
在步骤(4)中,极值点的判断是每个带通滤波器输出幅值和相邻的两个带通滤波器输出幅值做比较:
针对电流的幅值序列[Mi1,Mi2,...,Min],假设第i个带通滤波器输出幅值为Mii,该带通滤波器的中心频率为f0i,则第i-1和i+2个滤波器的输出幅值为Mi(i-1)和Mi(i+1),则:
Mii、Mi(i-1)、Mi(i+1)分别等于一时刻相电流第i、i-1、i+1个带通滤波器输出选取的数据窗的最大值的绝对值;
如果Mii≥Mi(i-1)并且Mii≥Mi(i+1),则判断Mii是极值点,记作Mssoi;
振荡频率的初选频率fsso=f0i;
其中:fs为采样频率;
对于电压的幅值序列[Mu1,Mu2,...,Mun]或功率的幅值序列[Mp1,Mp2,...,Mpn],其极值点的判断方式和电流的幅值序列相同。
4.根据权利要求1所述的在线自适应频率变化的次同步振荡间谐波提取方法,其特征在于:
在步骤(5)中,对极值点对应的带通滤波器实时波形数据采用过零监测+数据拟合的方法得到复选后的振荡频率频率,采样间隔时间第一个上升沿过零点之前的采样点序号为nl,对应的采样值为vl;第一个上升沿过零点之后的采样点序号为nl+1,对应的采样值为vl+1;第二个上升沿过零点之前的采样点序号为nm,对应的采样值为vm;第二个上升沿过零点之前的采样点序号为nm+1,对应的采样值为vm+1;采取一阶线性拟合的方式或者高阶多项式拟合的方式对实时数据进行拟合,获取过零点时刻;
过零点到nl+1的时间ΔT1为:
<mrow>
<msub>
<mi>&Delta;T</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>|</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>|</mo>
</mrow>
<mrow>
<mo>(</mo>
<mo>|</mo>
<msub>
<mi>v</mi>
<mi>l</mi>
</msub>
<mo>|</mo>
<mo>+</mo>
<mo>|</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>|</mo>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>n</mi>
<mrow>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>n</mi>
<mi>l</mi>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mo>|</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>|</mo>
</mrow>
<mrow>
<mo>(</mo>
<mo>|</mo>
<msub>
<mi>v</mi>
<mi>l</mi>
</msub>
<mo>|</mo>
<mo>+</mo>
<mo>|</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>|</mo>
<mo>)</mo>
</mrow>
</mfrac>
<msub>
<mi>T</mi>
<mi>s</mi>
</msub>
</mrow>
同理,第二个过零点到nl+1的时间ΔT2为:
<mrow>
<msub>
<mi>&Delta;T</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>|</mo>
<msub>
<mi>v</mi>
<mi>m</mi>
</msub>
<mo>|</mo>
</mrow>
<mrow>
<mo>(</mo>
<mo>|</mo>
<msub>
<mi>v</mi>
<mi>m</mi>
</msub>
<mo>|</mo>
<mo>+</mo>
<mo>|</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>m</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>|</mo>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>n</mi>
<mrow>
<mi>m</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>n</mi>
<mi>m</mi>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mo>|</mo>
<msub>
<mi>v</mi>
<mi>m</mi>
</msub>
<mo>|</mo>
</mrow>
<mrow>
<mo>(</mo>
<mo>|</mo>
<msub>
<mi>v</mi>
<mi>m</mi>
</msub>
<mo>|</mo>
<mo>+</mo>
<mo>|</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>m</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>|</mo>
<mo>)</mo>
</mrow>
</mfrac>
<msub>
<mi>T</mi>
<mi>s</mi>
</msub>
</mrow>
复选后的振荡频率[f’sso1,f’sso2,...,f’ssom],第i个振荡频率的表达式为:
<mrow>
<msubsup>
<mi>f</mi>
<mrow>
<mi>s</mi>
<mi>s</mi>
<mi>o</mi>
<mi>i</mi>
</mrow>
<mo>,</mo>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msub>
<mi>&Delta;T</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&Delta;T</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>n</mi>
<mi>m</mi>
</msub>
<mo>-</mo>
<msub>
<mi>n</mi>
<mrow>
<mi>l</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>T</mi>
<mi>S</mi>
</msub>
</mrow>
</mfrac>
<mo>.</mo>
</mrow>
2
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710338027.5A CN107179440B (zh) | 2017-05-15 | 2017-05-15 | 一种在线自适应频率变化的次同步振荡间谐波提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710338027.5A CN107179440B (zh) | 2017-05-15 | 2017-05-15 | 一种在线自适应频率变化的次同步振荡间谐波提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107179440A true CN107179440A (zh) | 2017-09-19 |
CN107179440B CN107179440B (zh) | 2019-12-10 |
Family
ID=59831161
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710338027.5A Active CN107179440B (zh) | 2017-05-15 | 2017-05-15 | 一种在线自适应频率变化的次同步振荡间谐波提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107179440B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107576842A (zh) * | 2017-10-12 | 2018-01-12 | 惠州市成惠电子科技有限公司 | 一种宽频带同步采样方法 |
CN108490254A (zh) * | 2018-01-29 | 2018-09-04 | 国电南瑞科技股份有限公司 | 一种基于多算法融合的宽带多频信号测量方法 |
CN109301842A (zh) * | 2018-10-23 | 2019-02-01 | 国网吉林省电力有限公司 | 一种基于负阻尼贡献的风电场次同步振荡切除方法 |
CN109546663A (zh) * | 2017-09-21 | 2019-03-29 | 南京南瑞继保电气有限公司 | 一种次同步振荡抑制系统阻尼控制信号调整方法 |
CN111398676A (zh) * | 2020-04-08 | 2020-07-10 | 电子科技大学 | 一种通过次/超同步振荡频率判定振荡传播路径的方法 |
CN111896796A (zh) * | 2020-08-03 | 2020-11-06 | 珠海格力电器股份有限公司 | 一种pfc电路的电流振荡检测装置、方法和电源系统 |
CN115065055A (zh) * | 2022-07-07 | 2022-09-16 | 中广核新能源安徽有限公司 | 一种风电汇集站间谐波控制方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6073370A (ja) * | 1983-09-30 | 1985-04-25 | Anritsu Corp | 周波数選択性のあるレベル測定装置 |
CN103487652A (zh) * | 2013-09-03 | 2014-01-01 | 电子科技大学 | 一种频率自适应实时分次谐波检测方法 |
CN105223418A (zh) * | 2015-09-22 | 2016-01-06 | 清华大学 | 次同步和超同步谐波相量的测量方法及测量装置 |
CN105470970A (zh) * | 2015-11-16 | 2016-04-06 | 清华大学 | 电力系统中次同步和超同步振荡的保护方法及系统 |
CN105606895A (zh) * | 2016-01-07 | 2016-05-25 | 国家电网公司 | 电力系统次同步振荡成份的在线检测及滤除方法 |
CN106154035A (zh) * | 2016-06-20 | 2016-11-23 | 哈尔滨工业大学 | 一种快速谐波及间谐波检测方法 |
CN106483375A (zh) * | 2016-09-22 | 2017-03-08 | 电子科技大学 | 一种多频率分次谐波检测方法 |
CN106645919A (zh) * | 2016-12-26 | 2017-05-10 | 北京四方继保自动化股份有限公司 | 基于三相瞬时功率的电网全谱功率振荡相量同步测量方法 |
-
2017
- 2017-05-15 CN CN201710338027.5A patent/CN107179440B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6073370A (ja) * | 1983-09-30 | 1985-04-25 | Anritsu Corp | 周波数選択性のあるレベル測定装置 |
CN103487652A (zh) * | 2013-09-03 | 2014-01-01 | 电子科技大学 | 一种频率自适应实时分次谐波检测方法 |
CN105223418A (zh) * | 2015-09-22 | 2016-01-06 | 清华大学 | 次同步和超同步谐波相量的测量方法及测量装置 |
CN105470970A (zh) * | 2015-11-16 | 2016-04-06 | 清华大学 | 电力系统中次同步和超同步振荡的保护方法及系统 |
CN105606895A (zh) * | 2016-01-07 | 2016-05-25 | 国家电网公司 | 电力系统次同步振荡成份的在线检测及滤除方法 |
CN106154035A (zh) * | 2016-06-20 | 2016-11-23 | 哈尔滨工业大学 | 一种快速谐波及间谐波检测方法 |
CN106483375A (zh) * | 2016-09-22 | 2017-03-08 | 电子科技大学 | 一种多频率分次谐波检测方法 |
CN106645919A (zh) * | 2016-12-26 | 2017-05-10 | 北京四方继保自动化股份有限公司 | 基于三相瞬时功率的电网全谱功率振荡相量同步测量方法 |
Non-Patent Citations (2)
Title |
---|
王彭 等: "基于时域准同步的谐波和间谐波检测算法", 《仪器仪表学报》 * |
谢小荣 等: "电力系统次同步和超同步谐波相量的检测方法", 《电力系统自动化》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109546663A (zh) * | 2017-09-21 | 2019-03-29 | 南京南瑞继保电气有限公司 | 一种次同步振荡抑制系统阻尼控制信号调整方法 |
CN109546663B (zh) * | 2017-09-21 | 2021-12-10 | 南京南瑞继保电气有限公司 | 一种次同步振荡抑制系统阻尼控制信号调整方法 |
CN107576842A (zh) * | 2017-10-12 | 2018-01-12 | 惠州市成惠电子科技有限公司 | 一种宽频带同步采样方法 |
CN107576842B (zh) * | 2017-10-12 | 2020-03-24 | 惠州市成惠电子科技有限公司 | 一种宽频带同步采样方法 |
CN108490254A (zh) * | 2018-01-29 | 2018-09-04 | 国电南瑞科技股份有限公司 | 一种基于多算法融合的宽带多频信号测量方法 |
CN109301842A (zh) * | 2018-10-23 | 2019-02-01 | 国网吉林省电力有限公司 | 一种基于负阻尼贡献的风电场次同步振荡切除方法 |
CN109301842B (zh) * | 2018-10-23 | 2022-04-22 | 国网吉林省电力有限公司 | 一种基于负阻尼贡献的风电场次同步振荡切除方法 |
CN111398676A (zh) * | 2020-04-08 | 2020-07-10 | 电子科技大学 | 一种通过次/超同步振荡频率判定振荡传播路径的方法 |
CN111398676B (zh) * | 2020-04-08 | 2021-03-30 | 电子科技大学 | 一种通过次/超同步振荡频率判定振荡传播路径的方法 |
CN111896796A (zh) * | 2020-08-03 | 2020-11-06 | 珠海格力电器股份有限公司 | 一种pfc电路的电流振荡检测装置、方法和电源系统 |
CN111896796B (zh) * | 2020-08-03 | 2021-07-20 | 珠海格力电器股份有限公司 | 一种pfc电路的电流振荡检测装置、方法和电源系统 |
CN115065055A (zh) * | 2022-07-07 | 2022-09-16 | 中广核新能源安徽有限公司 | 一种风电汇集站间谐波控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107179440B (zh) | 2019-12-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107179440A (zh) | 一种在线自适应频率变化的次同步振荡间谐波提取方法 | |
CN105223418A (zh) | 次同步和超同步谐波相量的测量方法及测量装置 | |
CN103713237B (zh) | 一种电力系统输电线路短路故障诊断方法 | |
CN106199183B (zh) | 一种实现次同步振荡在线辨识告警的pmu和方法 | |
CN104638664B (zh) | 风电厂次同步振荡保护方法及保护装置 | |
Xie et al. | Improved synchrophasor measurement to capture sub/super‐synchronous dynamics in power systems with renewable generation | |
Gao et al. | Waveform-based method for fast and accurate identification of subsynchronous resonance events | |
CN103605016B (zh) | 一种电能质量数据处理方法和装置 | |
CN108802566B (zh) | 一种基于hht信号分析的接地选线方法 | |
CN108152586A (zh) | 一种电力系统宽频带多振荡模式分量的辨识方法 | |
CN106771594A (zh) | 一种电力系统的次/超同步谐波检测方法 | |
CN103454537A (zh) | 基于小波分析的风力发电低电压穿越检测设备及方法 | |
Wang et al. | Real-time subsynchronous control interaction monitoring using improved intrinsic time-scale decomposition | |
CN110618338B (zh) | 一种电能质量瞬态测量方法、装置和设备 | |
CN109301842B (zh) | 一种基于负阻尼贡献的风电场次同步振荡切除方法 | |
Lu et al. | A wavelet/PSO based voltage regulation scheme and suitability analysis of copper-and aluminum-rotor induction machines for distributed wind power generation | |
Damala et al. | A simple decision tree-based disturbance monitoring system for VSC-based HVDC transmission link integrating a DFIG wind farm | |
Vanfretti et al. | A PMU-based fast real-time sub-synchronous oscillation detection application | |
CN108106827B (zh) | 一种旋翼台升转速状态下放大器端子松脱识别方法 | |
Gutierrez et al. | Power quality in grid-connected wind turbines | |
CN107860969A (zh) | 风力发电机组电能信息采集装置、系统和方法 | |
CN108088546B (zh) | 一种旋翼台工作转速状态下放大器端子松脱识别方法 | |
Guajardo et al. | Error compensation in distance relays caused by wind power plants in the power grid | |
CN107144767B (zh) | 一种故障指示装置及故障信号检测方法 | |
CN102467614A (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 |