CN102305891A - 一种电力系统低频振荡在线监测方法 - Google Patents
一种电力系统低频振荡在线监测方法 Download PDFInfo
- Publication number
- CN102305891A CN102305891A CN201110184516A CN201110184516A CN102305891A CN 102305891 A CN102305891 A CN 102305891A CN 201110184516 A CN201110184516 A CN 201110184516A CN 201110184516 A CN201110184516 A CN 201110184516A CN 102305891 A CN102305891 A CN 102305891A
- Authority
- CN
- China
- Prior art keywords
- oscillation
- frequency
- mode
- amplitude
- tfa
- 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
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明涉及一种电力系统低频振荡在线监测方法。本发明创造性地利用时频原子复带通滤波函数在振荡频率区间作并行滤波计算,获取振荡模式个数、实际振荡频率分布和实时幅值信息。基于相邻实时幅值信息的冗余性,利用最小二乘优化估计能获取各振荡模式幅值和衰减时间常数。基于已获取的振荡频率分布、幅值和衰减时间常数可以建立低频振荡信号模型,其中只有各模式初相位和直流分量幅值未知。利用粒子群优化算法对该模型进行优化估计,获取初相位和直流分量幅值。该方法对噪声鲁棒性好,能准确辨识各复合振荡模式,有助于电力系统强非线性模式分析,便于在线监测应用。
Description
技术领域
本发明涉及一种在线监测方法,尤其涉及一种电力系统低频振荡在线监测方法。
背景技术
监测和分析电力系统在扰动下的动态行为、确定其动态特性并研究适当的对策(包括各种控制措施)是电力系统设计和运行中最重要也是最复杂的任务之一。电力系统动态特性的主要监测和控制对象是机电模式振荡的参数,即频率、阻尼和幅值等。大区电网互联初期的弱交流联系、电力市场形式下趋于稳定极限的运行方式更易于诱发各区域电网机群间的低频振荡,因此动态特性的监测和分析更加重要。
在建设坚强智能电网架构下,随着全国电网互联格局形成,大容量机组不断投运,快速励磁系统普遍使用,低频振荡问题日益严重,这不但限制了区域间的功率传输,严重时还会导致系统解列或失稳。实现低频振荡有效监测能为系统决策和广域阻尼控制提供可靠依据,对于抑制振荡,保持系统安全稳定运行有重要意义。
广域测量系统(wide area measurement system,WAMS)通过同步相量测量单元(phasor measurement unit,PMU)来实现全网情景知晓。基于WAMS实测数据监测低频振荡避免了对复杂高维非线性系统建模带来的维数灾问题。目前常用方法包括时频分析法和参数法。时频分析法有傅里叶算法、小波算法和希尔伯特-黄变换(Hilbert-huang Transform,HHT)等。傅里叶算法和小波算法对噪声鲁棒性好,能通过频谱分析提供主导振荡模式频率,但傅里叶算法难以提供时域局部化信息,较难描述信号衰减特性。小波脊算法和小波能量系数法能同时给出时频域振荡特性,但在固定振荡模式频率约束下,小波的时频域带宽难以灵活调整,导致相邻尺度频带易发生混叠,影响模式分辨率和辨识精度,且观测时窗较长,不利于分析时变信号。HHT算法通过经验模态分解(empirical mode decomposition,EMD)能有效剥离各振荡模式,并对非线性模式有较好跟踪能力,但EMD分解的物理意义尚不明确,缺乏验证固有模态分量为实际振荡模式还是虚假分量的理论依据,且没有给出算法观测时窗长度和模式分辨能力的定量分析。
参数法主要包括Prony、ESPRIT和矩阵束等。Prony在工程中得到广泛运用,但该方法对噪声较为敏感,且对非平稳振荡模式拟合能力不理想。ESPRIT和矩阵束通过对构造的Hankel矩阵进行奇异值分解,划分信噪空间并剔除噪声干扰因子,能有效提高算法对噪声鲁棒性,但算法的模式分辨能力、观测时间窗长度以及对非平稳信号的监测能力都需要进一步研究。
因此,必须研究具有较高模式分辨率,较短观测时窗长度和适用于非平稳振荡模式分析的电力系统低频振荡在线监测方法。
发明内容
一种电力系统低频振荡在线监测方法,其特征在于,包括以下步骤:
步骤1,利用时频原子复带通滤波函数在振荡频率区间对低频振荡信号进行初次并行滤波计算,获取该振荡频率区间的振荡模式个数、实际振荡频率分布和实时幅值信息;
步骤2,基于相邻实时幅值信息的冗余性,利用最小二乘优化获取各振荡模式幅值和衰减时间常数;
步骤3,根据步骤1已获取的振荡频率分布以及步骤2获取的各振荡模式幅值和衰减时间常数建立低频振荡信号模型,并根据粒子群优化算法对该模型进行优化估计,获取初相位和直流分量幅值。
本发明创造性地利用时频原子复带通滤波函数在低频带作并行复带通滤波获取各低频振荡模式频率和时域信息,并结合智能优化算法获取各振荡模式相位和直流分量。本发明在噪声干扰条件下对低频振荡波形中区间和局部振荡模式的频率、阻尼、幅值、相位和直流分量幅值参数具有较高的测量精度,具有较高的模式分辨率和较短的观测时窗长度,动态特性好,特别适用于非平稳振荡模式分析,为基于WAMS量测数据的电力系统动态特性在线监测和控制创造了条件。
在上述的低频振荡在线监测方法,所述的步骤1中,定义时频原子复带通滤波函数为对一具有低通特性的实函数η(t)进行频域调制和伸缩处理构造TFA复带通滤波函数:
式中,ε为带通宽度伸缩参数,ωk为带通中心频率调制参数,其中,k=1,2,L,N,N为所设置的TFA单元个数;上式对应的频域表达式为:
式中,Z(ε,ωk,t)为t时刻s(t)在以ωk为中心频率的带通内的滤波输出复相量,若s(t)在该频域带通内含有振荡模式分量,则Ak(t)为包含时域衰减特征的实时幅值信息,θk(t)为相位信息,对应振荡频率为:
具体操作方法如下:
步骤1.1,设置初始带通宽度伸缩参数ε0,则初始模式分辨率为2Δω0,观测时窗长度为2Δt0,其中Δω0、Δt0分别为对应于ε0的时频原子频率带宽半径和时窗半径,定义如下:
并设置所覆盖的频率区间大于低频振荡频率范围0.1Hz-2.5Hz,以及所需的TFA单元个数N;当N个TFA单元录满2Δt0时长采样数据时,获取每个单元输出的对应于2Δt0采样时刻的实时幅值信息Ak(2Δt0)和振荡频率fk(2Δt0);
步骤1.2,对完成步骤1.1的振荡模式单元的极大值进行检验,对于通过检测Ak(2Δt0)为极大值的TFA单元,极大值的TFA单元个数即为振荡模式个数R(R<N),并对该R个模式设置TFA滤波通道,并重新编号为ωi(i=1,2,L,R);否则屏蔽该单元;由此获取被测信号中含有的振荡模式个数R和振荡频率分布ωi;
步骤1.3,对完成步骤1.2的TFA单元进行重调参数:重设Δω为Δωa/2,并通过高斯窗半径计算公式获取相应的带通宽度伸缩参数ε和Δt值,据此更新各并行TFA单元参数,并发出参数重调指令。
在上述的低频振荡在线监测方法,所述的步骤2中,振荡模式幅值和衰减时间常数方法如下:
接受到重调指令后,输入新的采样点,利用更新后的TFA单元以重调后的计算数据窗进行分析,得到各振荡模式的瞬时频率fi(t)和包含衰减特征的Ai(t),其中,定义实时幅值信息Ak(t)为Mk为振荡模式幅值,αk为衰减时间常数。然后取两个不同采样值对应的Ak(t),联立解得振荡模式幅值Mk和衰减时间常数αk。
在上述的低频振荡在线监测方法,所述的步骤3中,为了进一步获取初相位i=1,2,L,R,以及直流分量幅值M0,基于已知的各振荡模式频率、振荡模式幅值和衰减系数,建立低频振荡信号模型:
因此,本发明具有如下优点:能在噪声干扰下在线、准确获取低频振荡各振荡模式频率、幅值、阻尼、相位和直流分量幅值参数,具有较高的模式分辨率和较短的观测时窗长度,动态特性好,适用于非平稳、非线性振荡模式分析。
附图说明
附图1是第k个TFA滤波单元计算示意图,
附图2是本发明的详细流程图,
附图3是有功功率振荡的现场数据波形,
附图4是本发明方法和Prony方法的测量效果对比图。
具体实施方式
下面通过实施例,并结合附图,对本发明的技术方案作进一步具体的说明。
实施例:
首先介绍一下本发明需要涉及到的参数设置原则:
低频振荡的频率范围为0.1Hz-2.5Hz,每个TFA单元的频域带宽为2Δω。为了使滤波计算能覆盖全部低频区间,在低频区间内并行设置N个如图1所示TFA单元,其中N应取不小于2π(2.5-0.1)/(2Δω)的整数,此时每个单元对应的中心频率为ωk=0.2π+2Δω(k-1),k=1,2,L,N,各单元所能达到的模式分辨率为2Δω。当Δω取较小数值以提高模式分辨率时,单元数N会增加,降低了运算效率;当减少单元数以提高运算效率时,Δω取值变大,模式分辨率会随之降低。为此,本发明先以较高模式分辨率对振荡信号作初次并行滤波计算,通过极大值检测获取实际振荡模式频率分布,再针对性地设置滤波通道,并重调ε和Δω参数,在保证模式分辨率的同时减少并行滤波单元数,并缩短计算数据窗长度,从而提高运算效率。
本发明的电力系统低频振荡在线监测方法,包括以下步骤:
步骤1,利用时频原子复带通滤波函数在振荡频率区间对低频振荡信号作进行初次并行滤波计算,获取该振荡频率区间的振荡模式个数、实际振荡频率分布和实时幅值信息;首先定义时频原子复带通滤波函数为对一具有低通特性的实函数η(t)进行频域调制和伸缩处理构造TFA复带通滤波函数:
式中,Z(ε,ωk,t)为t时刻s(t)在以ωk为中心频率的带通内的滤波输出复相量,若s(t)在该频域带通内含有振荡模式分量,则Ak(t)为包含时域衰减特征的实时幅值信息,θk(t)为相位信息,对应振荡频率为:
实际计算中通常以差分运算代替式(11)中的求导运算。将式(9)代入式(10)的频域表达式可得TFA复带通滤波的频域变换式如下所示
式中:S(ω)为s(t)的傅里叶变换。调制参数ωk决定TFA复带通的中心频率,ε决定带通有效宽度。构造不同中心频率对应的TFA单元,则ωk对应的第k个TFA单元滤波计算示意图如附图1所示
由式(9)和式(12)可知,所构造的TFA在频域上呈现以ωk为中心频率,有效带宽受ε调整的带通特性,设有效带宽半径为Δω,则该带通起始频率为ωk-Δω,截止频率为ωk+Δω,且带通为紧支撑,能滤波输出带通内模式分量,而同时阻隔带通外其他频率分量和杂散信号的干扰。
由式(9)可知,TFA的中心频率和频域带宽分别受ωk和ε两个独立参数约束,二者可以分别设置。在满足海森堡测不准原理下,TFA滤波函数的观测时窗长度2Δt也能经ε灵活调整,因此TFA单元具有柔性滤波控制特性。根据式(8)、(9)所示TFA时频域表达式,利用高斯窗半径等量定义公式,可以获取不同ε对应的频域带宽半径Δω和观测时窗半径Δt。在实际应用时,通过灵活设置ε可以获取紧凑的频域带宽和较短的观测时窗长度,从而满足对复合、时变振荡模式分析要求。
步骤1的具体操作方法如下:
步骤1.1,设置初始带通宽度伸缩参数ε0,则初始模式分辨率为2Δω0,观测时窗长度为2Δt0,其中Δω0、Δt0分别为对应于ε0的时频原子频率带宽半径和时窗半径,定义如下:
并设置所覆盖的频率区间大于低频振荡频率范围0.1Hz-2.5Hz,以及所需的TFA单元个数N;当N个TFA单元录满2Δt0时长采样数据时,获取每个单元输出的对应于2Δt0采样时刻的实时幅值信息Ak(2Δt0)和振荡频率fk(2Δt0);
由PMU装置从电网采集的功角信号先储存于相量数据集中器(phasordata concentrator,PDC),再传入下一级分析软件。本发明中设置ε0=5时,Δω0=0.03πrad/s,Δt0=5.5s,此时每个滤波单元所能达到的模式分辨率为0.03Hz,需要的并行单元数N0=81,所覆盖的频率区间为0.085-2.515Hz。当TFA单元录满2Δt0=11s时长采样数据时,每个单元输出相应的Ak(Δt0)和fk(Δt0)。需要说明的是,该步骤中单元数N0,模式分辨率2Δω0以及观测时窗长度2Δt0均可由用户根据实际测量需要通过参数ε灵活设置;
步骤1.2,对完成步骤1.1的振荡模式单元的极大值进行检验,对于通过检测Ak(2Δt0)为极大值的TFA单元,极大值的TFA单元个数即为振荡模式个数R(R<N),并对该R个模式设置TFA滤波通道,并重新编号为ωi(i=1,2,L,R);否则屏蔽该单元;由此获取被测信号中含有的振荡模式个数R和振荡频率分布ωi;
步骤1.3,对完成步骤1.2的TFA单元进行重调参数:在初次滤波环节,为了获取较高的模式分辨率,设置的频域带宽较窄,并行滤波计算量较大。获取各实际模式频率分布后,根据频差,通过改变ε值重新调整未屏蔽TFA单元的频域带宽和时窗长度。由于各时变模式频率变化范围通常较小,大量实验研究表明:若Δωa为频域上相距最近的两个振荡模式频差的绝对值,则Δω应重设为Δωa/2,此时能满足分析要求。有高斯窗半径计算公式获取相应的ε和Δt值,据此更新各并行TFA单元参数,并对外发出参数重调指令。由实际频差对单元参数重调可以使该方法在满足模式分辨率要求的同时,适当增加频域带宽,从而缩短观测时窗长度,减少并行滤波次数,提高了方法运算效率。
步骤2,基于相邻实时幅值信息的冗余性,利用最小二乘优化获取各振荡模式幅值和衰减时间常数;即接受到参数重调指令后,当参数重调判断环节接受到重调指令后,输入新的采样点,直接利用更新后的TFA单元以重调后的计算数据窗进行分析,得到各振荡模式的瞬时频率fi(t)和包含衰减特征的实时幅值信息Ai(t)。其中Ai(t)可表示为取两个不同采样值对应的Ai(t),便可联立解得振荡模式幅值Mi和衰减时间常数αi。实际上,由于重调后的计算数据窗长度Δt小于初次滤波时的Δt0值,因此,采样数据信息是冗余的,可用最小二乘法(least square,LS)对参数进行优化估计。该步骤中,最小二乘法基于:
假设曲线
使得
成立,则称曲线y*(x)为在曲线族(xi,yi)中按最小二乘原则确定的对于数据(xi,yi)的拟合曲线。由线性无关向量组j=0,1,L,n做基底构成一个Rm+1的一个子空间,记A=[φ0,φ1,Lφn],y=(y0,y1,L ym)T。向量组满足条件(16)的拟合曲线y*(x)存在且唯一,并且从方程
ATAc*=ATy; (17)
步骤3,根据步骤1已获取的振荡频率分布以及步骤2获取的各振荡模式幅值和衰减时间常数建立低频振荡信号模型,并根据粒子群优化算法对该模型进行优化估计,获取初相位和直流分量幅值,基于已知的各振荡模式频率、幅值和衰减系数,建立低频振荡信号模型:
设搜索空间为D维,总粒子数为n,第i个粒子位置表示为向量xi=(xi1,xi2 L xiD),第i个粒子迄今位置搜索到的最优位置为pbesti=(Pi1,Pi2 L PiD),整个粒子群迄今为止搜索到的最优位置为gbest=(g1,g2 L gD),第i个粒子的位置变化率(速度)为相量vi=(vi1 vi2 L viD)。
粒子的每维速度和位置按如下公式进行变化:
vid(t+1)=vid(t)+c1*r1*(pid(t)-xid(t))+c2*r2*(pgd(t)-xid(t)); (19)
xid(t+1)=xid(t)+vid(t+1);1≤i≤n;1≤d≤D; (20)
其中c1,c2为正常数,称为加速因子,c1调节粒子飞向自身最好位置方向的步长;c2调节粒子向全局最好位置飞行的步长;r1,r2为[0,1]之间的随机数。为了减少在进化过程中,粒子离开探索空间的可能性,通常,第d(l≤d≤D)维的位置变化范围限定在内,速度变化范围限定在内(即在迭代中若vid和xid)超出了边界值,将之设为边界值)。
粒子群初始位置和速度随机产生,然后按式(19)和(20)进行迭代,直至找到满意的解。
并行TFA复带通滤波方法应用于低频振荡在线监测的实现流程如附图2所示。
本实施例针对实际电网,进行了相关测试及结果验证,具体情况如下:
A区域电网一变电站(HNHJ)相关传输线路发生有功功率振荡事件,该区域电网WAMS系统完整地记录下这一过程,从相量数据集中器(PDC)中提取该传输线有功功率振荡波形如附图3所示,为了便于说明,取出时段(23:00:16-23:00:31 UTC)中的振荡数据进行分析,将该数据导入并行TFA复带通滤波分析模块,采样频率为50Hz,仅列出振荡幅度较大的模式监测结果如表1所示。利用DSA中Prony分析工具进行主导模式辨识的结果如表2所示。
表1本发明方法对主导模式辨识结果
表2Prony方法对主导模式辨识结果
由表1和表2可知,本发明方法和Prony方法均辨识出0.808Hz模式,该模式振幅较大,阻尼最弱,并呈现负阻尼的增幅振荡形式,为该振荡事件中主导振荡模式,说明系统在该局部区域阻尼不足,应及时采取抑制措施消除振荡危害。除了直流分量外,其他模式或者振荡能量较小,或者阻尼较大,能在短时间内平息,对系统危害性较小。
根据本发明方法和Prony方法辨识的全部模式信息分别对原始波形作曲线拟合,结果如附图4所示。
附图4中本发明方法拟合信噪比为54dB,拟合曲线能很好逼近原始振荡波形,而Prony方法拟合信噪比为36dB,逼近精度略差,特别在数据末端信号非线性增强时,Prony的拟合偏差更大,这是由于Prony的线性化处理较难适应于强非线性电力系统模型的缘故。因此,本发明方法对噪声鲁棒性好,能较好地跟踪非平稳振荡模式变化,且分析准确性不依赖于模型阶数选取,对低频振荡模式监测的总体性能优于传统Prony方法。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (4)
1.一种电力系统低频振荡在线监测方法,其特征在于,包括以下步骤:
步骤1,利用时频原子复带通滤波函数在振荡频率区间对低频振荡信号进行初次并行滤波计算,获取该振荡频率区间的振荡模式个数、实际振荡频率分布和实时幅值信息;
步骤2,基于相邻实时幅值信息的冗余性,利用最小二乘优化获取各振荡模式幅值和衰减因子;
步骤3,根据步骤1已获取的振荡频率分布以及步骤2获取的各振荡模式幅值和衰减时间常数建立低频振荡信号模型,并根据粒子群优化算法对该模型进行优化估计,获取初相位和直流分量幅值。
2.根据权利要求1所述的一种电力系统低频振荡在线监测方法,其特征在于,所述的步骤1中,定义时频原子复带通滤波函数为对一具有低通特性的实函数η(t)进行频域调制和伸缩处理构造TFA复带通滤波函数:
式中,Z(ε,ωk,t)为t时刻s(t)在以ωk为中心频率的带通内的滤波输出复相量,若s(t)在该频域带通内含有振荡模式分量,则Ak(t)为包含时域衰减特征的实时幅值信息,θk(t)为相位信息,对应振荡频率为:
具体操作方法如下:
步骤1.1,设置初始带通宽度伸缩参数ε0,则初始模式分辨率为2Δω0,观测时窗长度为2Δt0,其中Δω0、Δt0分别为对应于ε0的时频原子频率带宽半径和时窗半径,定义如下:
并设置所覆盖的频率区间大于低频振荡频率范围0.1Hz-2.5Hz,以及所需的TFA单元个数N;当N个TFA单元录满2Δt0时长采样数据时,获取每个单元输出的对应于2Δt0采样时刻的实时幅值信息Ak(2Δt0)和振荡频率fk(2Δt0);
步骤1.2,对完成步骤1.1的振荡模式单元的极大值进行检验,对于通过检测Ak(2Δt0)为极大值的TFA单元,极大值的TFA单元个数即为振荡模式个数R,其中R<N,并对该R个模式设置TFA滤波通道,并重新编号为ωi,其中,i=1,2,L,R;否则屏蔽该单元;由此获取被测信号中含有的振荡模式个数R和振荡频率分布ωi;
步骤1.3,对完成步骤1.2的TFA单元进行重调参数:重设Δω为Δωa/2,并通过高斯窗半径计算公式获取相应的带通宽度伸缩参数ε和Δt值,据此更新各并行TFA单元参数,并发出参数重调指令。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110184516 CN102305891B (zh) | 2011-07-04 | 2011-07-04 | 一种电力系统低频振荡在线监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110184516 CN102305891B (zh) | 2011-07-04 | 2011-07-04 | 一种电力系统低频振荡在线监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102305891A true CN102305891A (zh) | 2012-01-04 |
CN102305891B CN102305891B (zh) | 2013-09-11 |
Family
ID=45379768
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110184516 Expired - Fee Related CN102305891B (zh) | 2011-07-04 | 2011-07-04 | 一种电力系统低频振荡在线监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102305891B (zh) |
Cited By (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102801175A (zh) * | 2012-08-09 | 2012-11-28 | 国电南瑞科技股份有限公司 | 一种pss相位补偿环节时间常数计算方法 |
CN102820665A (zh) * | 2012-05-30 | 2012-12-12 | 邵俊松 | 一种快速识别风电并网系统发生次同步震荡的方法 |
CN102969713A (zh) * | 2012-12-10 | 2013-03-13 | 长沙理工大学 | 一种电力系统低频振荡模式时频分析方法 |
CN103208808A (zh) * | 2013-03-07 | 2013-07-17 | 武汉大学 | 一种电力系统次同步振荡模态辨识方法 |
CN103337866A (zh) * | 2013-07-19 | 2013-10-02 | 中国南方电网有限责任公司 | 一种从随机响应数据中辨识电力系统低频振荡参数的方法 |
CN104007312A (zh) * | 2014-06-20 | 2014-08-27 | 西南大学 | 基于低频方波信号频率预测的信号频率实时测量方法 |
CN104198809A (zh) * | 2014-06-30 | 2014-12-10 | 中国南方电网有限责任公司 | 一种电力系统多频率振荡的频率测量方法 |
CN104753075A (zh) * | 2015-03-19 | 2015-07-01 | 中国农业大学 | 互联电力系统的主导振荡模态的辨识方法及装置 |
CN105867353A (zh) * | 2016-05-25 | 2016-08-17 | 广东电网有限责任公司电力科学研究院 | 一种振荡监测方法及系统 |
CN106338651A (zh) * | 2016-08-31 | 2017-01-18 | 长沙理工大学 | 应用于电力系统低频振荡模式识别的粒子滤波分析方法 |
CN106405237A (zh) * | 2016-08-31 | 2017-02-15 | 长沙理工大学 | 一种应用于多通道电力系统信号中的低频振荡模式识别的分析方法 |
CN106610590A (zh) * | 2017-01-25 | 2017-05-03 | 北京建筑大学 | 提高电力系统稳定性的迭代辨识广域阻尼控制器构建方法 |
CN107086566A (zh) * | 2017-04-19 | 2017-08-22 | 清华大学 | 基于广域信息的lmd互联电力系统低频振荡分析方法 |
CN107271768A (zh) * | 2017-05-26 | 2017-10-20 | 东南大学 | 一种最小二乘拟合动态频率测量方法 |
CN107545253A (zh) * | 2017-09-11 | 2018-01-05 | 国网江西省电力公司电力科学研究院 | 一种电力系统低频振荡监测方法 |
CN107681658A (zh) * | 2017-09-30 | 2018-02-09 | 广东电网有限责任公司电力科学研究院 | 一种面向调度主站的电网低频振荡分析测试方法及系统 |
CN107905270A (zh) * | 2017-11-28 | 2018-04-13 | 辽宁工程技术大学 | 一种深基坑变形预测方法 |
CN107947196A (zh) * | 2017-11-16 | 2018-04-20 | 国网四川省电力公司 | 一种基于改进粒子群算法的超低频振荡抑制方法 |
CN108152586A (zh) * | 2017-12-07 | 2018-06-12 | 北京四方继保自动化股份有限公司 | 一种电力系统宽频带多振荡模式分量的辨识方法 |
CN109425784A (zh) * | 2017-08-28 | 2019-03-05 | 通用电器技术有限公司 | 用于检测和评估电网中振荡的系统和方法 |
CN109542089A (zh) * | 2018-12-21 | 2019-03-29 | 浙江大学 | 一种基于改进变分模态分解的工业过程非线性振荡检测方法 |
CN109613617A (zh) * | 2019-01-24 | 2019-04-12 | 吉林大学 | 基于磁共振响应信号参数提取的地下水探测方法及系统 |
CN111505379A (zh) * | 2020-03-06 | 2020-08-07 | 天津大学 | 一种基于奇异值分解的配电网同步相量量测数据压缩方法 |
CN111596166A (zh) * | 2020-05-12 | 2020-08-28 | 国网天津市电力公司电力科学研究院 | 一种配电网不接地系统低频振荡冲击电流快速估算方法 |
CN112385109A (zh) * | 2018-07-06 | 2021-02-19 | 乌本产权有限公司 | 用于识别低频振荡的方法和用于其的检测装置 |
CN113960363A (zh) * | 2021-09-22 | 2022-01-21 | 国网新疆电力有限公司 | 基于多分量Morlet小波的高频振荡相量测量方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101368987A (zh) * | 2008-09-27 | 2009-02-18 | 武汉大学 | 一种电力系统稳态谐波和/或间谐波测量方法 |
CN101408586A (zh) * | 2008-11-28 | 2009-04-15 | 北京四方继保自动化股份有限公司 | 基于经验模态分解的在线低频振荡检测和节点同调分群法 |
WO2009048964A1 (en) * | 2007-10-09 | 2009-04-16 | Schweitzer Engineering Laboratories, Inc. | Real-time power system oscillation detection using modal analysis |
CN101447676A (zh) * | 2008-12-01 | 2009-06-03 | 中国电力科学研究院 | 一种电力系统低频振荡分析方法 |
-
2011
- 2011-07-04 CN CN 201110184516 patent/CN102305891B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009048964A1 (en) * | 2007-10-09 | 2009-04-16 | Schweitzer Engineering Laboratories, Inc. | Real-time power system oscillation detection using modal analysis |
CN101368987A (zh) * | 2008-09-27 | 2009-02-18 | 武汉大学 | 一种电力系统稳态谐波和/或间谐波测量方法 |
CN101408586A (zh) * | 2008-11-28 | 2009-04-15 | 北京四方继保自动化股份有限公司 | 基于经验模态分解的在线低频振荡检测和节点同调分群法 |
CN101447676A (zh) * | 2008-12-01 | 2009-06-03 | 中国电力科学研究院 | 一种电力系统低频振荡分析方法 |
Cited By (44)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102820665B (zh) * | 2012-05-30 | 2015-01-07 | 江苏创航电气有限公司 | 一种快速识别风电并网系统发生次同步震荡的方法 |
CN102820665A (zh) * | 2012-05-30 | 2012-12-12 | 邵俊松 | 一种快速识别风电并网系统发生次同步震荡的方法 |
CN102801175A (zh) * | 2012-08-09 | 2012-11-28 | 国电南瑞科技股份有限公司 | 一种pss相位补偿环节时间常数计算方法 |
CN102969713A (zh) * | 2012-12-10 | 2013-03-13 | 长沙理工大学 | 一种电力系统低频振荡模式时频分析方法 |
CN102969713B (zh) * | 2012-12-10 | 2014-12-10 | 长沙理工大学 | 一种电力系统低频振荡模式时频分析方法 |
CN103208808A (zh) * | 2013-03-07 | 2013-07-17 | 武汉大学 | 一种电力系统次同步振荡模态辨识方法 |
CN103208808B (zh) * | 2013-03-07 | 2015-01-14 | 武汉大学 | 一种电力系统次同步振荡模态辨识方法 |
CN103337866A (zh) * | 2013-07-19 | 2013-10-02 | 中国南方电网有限责任公司 | 一种从随机响应数据中辨识电力系统低频振荡参数的方法 |
CN103337866B (zh) * | 2013-07-19 | 2015-04-29 | 中国南方电网有限责任公司 | 一种从随机响应数据中辨识电力系统低频振荡参数的方法 |
CN104007312A (zh) * | 2014-06-20 | 2014-08-27 | 西南大学 | 基于低频方波信号频率预测的信号频率实时测量方法 |
CN104007312B (zh) * | 2014-06-20 | 2017-01-18 | 西南大学 | 基于低频方波信号频率预测的信号频率实时测量方法 |
CN104198809A (zh) * | 2014-06-30 | 2014-12-10 | 中国南方电网有限责任公司 | 一种电力系统多频率振荡的频率测量方法 |
CN104198809B (zh) * | 2014-06-30 | 2015-05-20 | 中国南方电网有限责任公司 | 一种电力系统多频率振荡的频率测量方法 |
CN104753075A (zh) * | 2015-03-19 | 2015-07-01 | 中国农业大学 | 互联电力系统的主导振荡模态的辨识方法及装置 |
CN105867353A (zh) * | 2016-05-25 | 2016-08-17 | 广东电网有限责任公司电力科学研究院 | 一种振荡监测方法及系统 |
CN105867353B (zh) * | 2016-05-25 | 2018-06-05 | 广东电网有限责任公司电力科学研究院 | 一种振荡监测方法及系统 |
CN106338651A (zh) * | 2016-08-31 | 2017-01-18 | 长沙理工大学 | 应用于电力系统低频振荡模式识别的粒子滤波分析方法 |
CN106338651B (zh) * | 2016-08-31 | 2018-09-14 | 长沙理工大学 | 应用于电力系统低频振荡模式识别的粒子滤波分析方法 |
CN106405237B (zh) * | 2016-08-31 | 2018-10-02 | 长沙理工大学 | 一种应用于多通道电力系统信号中的低频振荡模式识别的分析方法 |
CN106405237A (zh) * | 2016-08-31 | 2017-02-15 | 长沙理工大学 | 一种应用于多通道电力系统信号中的低频振荡模式识别的分析方法 |
CN106610590A (zh) * | 2017-01-25 | 2017-05-03 | 北京建筑大学 | 提高电力系统稳定性的迭代辨识广域阻尼控制器构建方法 |
CN107086566A (zh) * | 2017-04-19 | 2017-08-22 | 清华大学 | 基于广域信息的lmd互联电力系统低频振荡分析方法 |
CN107086566B (zh) * | 2017-04-19 | 2019-11-22 | 清华大学 | 基于广域信息的lmd互联电力系统低频振荡分析方法 |
CN107271768A (zh) * | 2017-05-26 | 2017-10-20 | 东南大学 | 一种最小二乘拟合动态频率测量方法 |
CN107271768B (zh) * | 2017-05-26 | 2019-06-21 | 东南大学 | 一种最小二乘拟合动态频率测量方法 |
CN109425784A (zh) * | 2017-08-28 | 2019-03-05 | 通用电器技术有限公司 | 用于检测和评估电网中振荡的系统和方法 |
CN107545253A (zh) * | 2017-09-11 | 2018-01-05 | 国网江西省电力公司电力科学研究院 | 一种电力系统低频振荡监测方法 |
CN107545253B (zh) * | 2017-09-11 | 2020-10-20 | 国网江西省电力公司电力科学研究院 | 一种电力系统低频振荡监测方法 |
CN107681658B (zh) * | 2017-09-30 | 2020-04-21 | 广东电网有限责任公司电力科学研究院 | 一种面向调度主站的电网低频振荡分析测试方法及系统 |
CN107681658A (zh) * | 2017-09-30 | 2018-02-09 | 广东电网有限责任公司电力科学研究院 | 一种面向调度主站的电网低频振荡分析测试方法及系统 |
CN107947196B (zh) * | 2017-11-16 | 2019-09-20 | 国网四川省电力公司 | 一种基于改进粒子群算法的超低频振荡抑制方法 |
CN107947196A (zh) * | 2017-11-16 | 2018-04-20 | 国网四川省电力公司 | 一种基于改进粒子群算法的超低频振荡抑制方法 |
CN107905270A (zh) * | 2017-11-28 | 2018-04-13 | 辽宁工程技术大学 | 一种深基坑变形预测方法 |
WO2019109969A1 (zh) * | 2017-12-07 | 2019-06-13 | 北京四方继保自动化股份有限公司 | 一种电力系统宽频带多振荡模式分量的辨识方法 |
CN108152586A (zh) * | 2017-12-07 | 2018-06-12 | 北京四方继保自动化股份有限公司 | 一种电力系统宽频带多振荡模式分量的辨识方法 |
CN112385109A (zh) * | 2018-07-06 | 2021-02-19 | 乌本产权有限公司 | 用于识别低频振荡的方法和用于其的检测装置 |
CN109542089B (zh) * | 2018-12-21 | 2020-07-17 | 浙江大学 | 一种基于改进变分模态分解的工业过程非线性振荡检测方法 |
CN109542089A (zh) * | 2018-12-21 | 2019-03-29 | 浙江大学 | 一种基于改进变分模态分解的工业过程非线性振荡检测方法 |
CN109613617A (zh) * | 2019-01-24 | 2019-04-12 | 吉林大学 | 基于磁共振响应信号参数提取的地下水探测方法及系统 |
CN111505379A (zh) * | 2020-03-06 | 2020-08-07 | 天津大学 | 一种基于奇异值分解的配电网同步相量量测数据压缩方法 |
CN111505379B (zh) * | 2020-03-06 | 2021-07-16 | 天津大学 | 一种基于奇异值分解的配电网同步相量量测数据压缩方法 |
CN111596166A (zh) * | 2020-05-12 | 2020-08-28 | 国网天津市电力公司电力科学研究院 | 一种配电网不接地系统低频振荡冲击电流快速估算方法 |
CN113960363A (zh) * | 2021-09-22 | 2022-01-21 | 国网新疆电力有限公司 | 基于多分量Morlet小波的高频振荡相量测量方法及装置 |
CN113960363B (zh) * | 2021-09-22 | 2023-09-15 | 国网新疆电力有限公司 | 基于多分量Morlet小波的高频振荡相量测量方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN102305891B (zh) | 2013-09-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102305891B (zh) | 一种电力系统低频振荡在线监测方法 | |
CN103454497B (zh) | 基于改进加窗离散傅立叶变换的相位差测量方法 | |
Sirojan et al. | Embedded edge computing for real-time smart meter data analytics | |
CN101915874B (zh) | 一种基于傅立叶变换的谐波检测方法 | |
CN101975910A (zh) | 一种特高压直流输电线路故障智能分类与测距方法 | |
CN104865474A (zh) | 一种基于pmu数据实时监测低频振荡源的方法 | |
CN106841778A (zh) | 基于pmu实现的次同步和超同步谐波参数的处理方法 | |
CN103197218A (zh) | 一种高压电缆绝缘缺陷局放带电检测诊断方法 | |
CN102323480A (zh) | 一种基于Hilbert-Huang变换的电能质量分析方法 | |
CN103901273B (zh) | 电网谐波检测方法和电网谐波检测装置 | |
CN105004939A (zh) | 一种复合电能质量扰动信号定量分析方法 | |
CN106980044B (zh) | 一种适应风电接入的电力系统谐波电流估计方法 | |
CN103941088A (zh) | 一种基于三相信号的电力系统频率快速测量方法 | |
CN107576851B (zh) | 基于旋转电流相量的系统谐波阻抗测量方法 | |
CN102401870A (zh) | 基于Hilbert-huang变换和遗传算法的配电网单相接地故障测距方法 | |
Cui et al. | HVDC transmission line fault localization base on RBF neural network with wavelet packet decomposition | |
CN106250904A (zh) | 基于改进s变换的电能扰动分析仪及分类方法 | |
CN111308260A (zh) | 一种基于小波神经网络的电能质量监测和电器故障分析系统及其工作方法 | |
Samal et al. | New signal subspace approach to estimate the inter‐area oscillatory modes in power system using TLS‐ESPRIT algorithm | |
CN102323487A (zh) | 基于谐波分量的输电线路零序电容抗干扰测量方法 | |
CN205539204U (zh) | 一种变频接地阻抗测量仪 | |
CN113466548A (zh) | 一种基于相量测量技术的智能电表台区识别方法 | |
CN103913654B (zh) | 高压直流输电系统接地极故障特征的识别方法以及设备 | |
CN109581045A (zh) | 一种满足iec标准框架的间谐波功率计量方法 | |
CN111638401A (zh) | 一种电容器在线监测系统及方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130911 Termination date: 20170704 |
|
CF01 | Termination of patent right due to non-payment of annual fee |