CN106099952B - 一种调速***引发低频振荡的机理识别方法 - Google Patents
一种调速***引发低频振荡的机理识别方法 Download PDFInfo
- Publication number
- CN106099952B CN106099952B CN201610515402.4A CN201610515402A CN106099952B CN 106099952 B CN106099952 B CN 106099952B CN 201610515402 A CN201610515402 A CN 201610515402A CN 106099952 B CN106099952 B CN 106099952B
- Authority
- CN
- China
- Prior art keywords
- oscillation
- frequency
- signal
- criterion
- governing system
- 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.)
- Active
Links
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/24—Arrangements for preventing or reducing oscillations of power in networks
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了属于电力***安全稳定领域的一种调速***引发低频振荡的机理识别方法。该方法通过对电网有关参数和调速***相关信号的测量,对***进行初步模态分析,结合对电网功率振荡信号频率和波形的分析,形成低频振荡产生机理的判据。结合低频振荡产生机理的判据和电网功率信号与调速***信号的振荡频率的一致性计算,形成判断调速***低频振荡机理的方法。
Description
技术领域
本发明属于电力***安全稳定领域,特别涉及电力***低频振荡机理分析。尤其涉及汽轮发电机调速***引发低频振荡的机理识别方法,适用于生产实践。
背景技术
互联电网可以提高电力***运行的经济性,然而也导致电力***动态稳定问题非常突出,国内外多次发生过低频振荡现象。现代快速、高放大倍数励磁***是引起负阻尼的主要原因,加装电力***稳定器PSS可以对这类低频振荡现象进行比较有效的抑制。随着调速***的快速发展,调速***动态对电力***动态稳定性的影响日益突出,***中多次出现由调速***引发的低频振荡现象,严重影响的***的安全稳定运行。
现代大型汽轮机调速***通常采用功频电液控制,其响应速度比以往机械液压式调速***大大提高,已经能够对电力***的动态产生较大影响。调速***包括转速测量机构、控制器、电液转换器、油动机以及阀门控制等,某个环节参数设置不合理或者出现故障均可能导致调节阀门摆动,进而造成电网低频功率振荡。
目前,电力***低频振荡的机理主要包括负阻尼机理和强迫振荡机理。负阻尼机理是解释自由振荡稳定性的成熟理论,而强迫振荡是指***在外施扰动源激励作用下产生的振荡。对于汽轮发电机调速***引起的振荡,究竟是负阻尼机理还是强迫振荡机理并没有完善的判断方法。本发明专利通过对调速***和电网有关参数的测量,形成判断调速***低频振荡机理的方法。
发明内容
本发明针对调速***引发的低频振荡,提出了一种调速***引发低频振荡的机理识别方法,其特征在于,所述方法包括
步骤1、计算互联电网的各振荡模式,根据特征值求各振荡模式的振荡频率与阻尼比;计算参与因子,参与因子Pki表示第i个模式中第k个变量的参与程度;
步骤2、实测信号,信号包括电网功率、发电机转速、调节阀门指令、等效阀位、阀位开度以及调节级压力;
步骤3、提取实测信号的振荡频率及其衰减因子;
步骤4、基于上述步骤和判据判断调速***低频振荡的机理。
所述步骤4的具体过程为
步骤401、对***进行初步模态分析;
根据广域测量***同步相量测量装置的监测数据,若是监测到至少两个区域的机组群发生低频振荡现象,则可判断***振荡模式为区域振荡模式,振荡产生机理为负阻尼振荡;若是仅监测到单个发电厂内的机组发生低频振荡现象,则可初步判断***振荡模式为本地振荡模式,振荡产生机理须进一步判断;
步骤402、判据F1:电网功率信号振荡频率;
由步骤1求得互联电网本地振荡模式的振荡频率,其中频率最小值为fmin,最大值为fmax;当电网功率信号的振荡频率fPe在[fmin,fmax]范围以外时,判据F1=0,振荡产生机理为强迫振荡;当fPe在[fmin,fmax]范围以内时,判据F1=1,振荡产生机理须进一步判断;
步骤403、判据F2:电网功率信号振荡波形;
在F1=1的基础上,起振阶段为增幅振荡而稳态阶段为等幅振荡时,判据F2=0,振荡产生机理为强迫振荡;起振阶段与稳态阶段均为增幅振荡时,判据F2=1,振荡产生机理为负阻尼振荡;
步骤404、判据F3:电网功率信号与调速***信号的振荡频率的一致性;
计算电网实测信号与调速***实测信号的振荡频率方差:
其中,fPe为电网功率信号的振荡频率,为调节阀门指令的振荡频率,fω为发电机转速信号的振荡频率,fcv为等效阀位信号的振荡频率,fgv为阀位开度信号的振荡频率,fPm为调节级压力信号的振荡频率;μ为所有实测信号的平均振荡频率;
当D(f)≤0.005时,判据F3=1,判断该强迫振荡是由调速***引起;当D(f)>0.005时,判据F3=0,判断该强迫振荡不是由调速***引起。
有益效果
本发明方法结合***结构、运行参数以及广域测量装置对***各信号的检测,形成识别***低频振荡产生机理的判据,能够判断***低频振荡产生的原因,从而能够更有针对性地采取控制措施抑制低频振荡。
附图说明
图1为调速***引发低频振荡的机理识别方法原理图;
图2为电网功率信号的在线检测图;
图3为调节阀门指令信号的在线检测图;
图4为等效阀位信号的在线检测图;
图5为阀位开度信号的在线检测图;
图6为调节级压力信号的在线检测图;
图7为发电机转速信号的在线检测图。
具体实施方式
本发明提出了一种判断调速***低频振荡机理的方法,图1为调速***引发低频振荡的机理识别方法原理图,具体步骤为
(1)计算互联电网的各振荡模式、参与因子;
(a)计算互联电网的振荡模式;
线性***的状态空间方程表示如下:
其中Δx是***的状态向量,Δu是***的输入向量,Δy为输出向量。矩阵A与电网***的发电机参数和***网络参数有关,通过求取矩阵A的特征值和特征向量,可以获得***模态。
***的特征方程与特征向量求解如下:
det(A-λI)=0 (2)
Aφi=λiφi (3)
由特征方程可以求得状态矩阵的特征值,其中每一对共轭复数特征值σ±jω对应***的一个振荡模式,并可根据特征值求得振荡的频率f与阻尼比ξ如下:
(b)计算参与因子;
由于式(2)中矩阵A为非对角阵,变量之间存在耦合,为解耦作如下变换,取Δx=ΦΔz,其中Φ为矩阵A对应的右模态矩阵。取Δz=ψΔx,其中ψ为矩阵A对应的左模态矩阵。则参与矩阵P=[P1 P2…Pn]结合左模态矩阵与右模态矩阵,作为状态变量与模式之间的一种度量。
其中
参与因子Pki表示第i个模式中第k个变量的参与程度。
(2)电网功率、发电机转速、调节阀门指令、等效阀位、阀位开度、调节级压力等信号的检测;
Prony算法是针对等间距采样点,假设模型是一些具有任意幅值、相位、频率和衰减因子的指数函数的线性组合,即认为测量输入x(0),…,x(n-1)的估计值可以表示为:
其中Am为幅值,fm为频率,am为衰减因子,θm为初相,Δt为时间间隔。
算法的主要步骤为:
(a)构造样本函数矩阵;
其中
(b)确定R阵的有效秩p;
(c)根据式(9)求解a;
(d)求多项式(10)的根z;
1+a1z-1+…+apz-p=0 (10)
并根据式(11)递推出
(e)计算参数b;
(f)利用式(13)计算振幅、相位、频率和衰减因子;
(3)实测信号的振荡频率和衰减因子提取;
以调速***的电网功率信号为例,通过Prony算法对在线监测的电网功率信号进行分析后,得到一组具有不同幅值、相位、频率和衰减因子的指数函数信号的线性组合。指数信号的幅值Aj最大时对应的信号为主导振荡信号。取主导振荡信号的频率fj取为振荡频率,取主导振荡信号的衰减因子ai为振荡衰减因子,即
则主导振荡模式的阻尼比由式(16)求得
同理可得调节阀门指令信号的振荡频率fuT,发电机转速信号的振荡频率fω、等效阀位信号的振荡频率fcv、阀位开度信号的振荡频率fgv、调节级压力信号的振荡频率fPm。
(4)判断调速***低频振荡的机理;
(a)对***进行初步模态分析;
根据广域测量***(WAMS)同步相量测量装置(PMU)的监测数据,若是监测到两个或多个区域的机组群发生低频振荡现象,则可判断***振荡模式为区域振荡模式,振荡产生机理为负阻尼振荡;若是仅监测到单个发电厂内的机组发生低频振荡现象,则可初步判断***振荡模式为本地振荡模式,振荡产生机理须进一步判断。
(b)判据F1:电网功率信号振荡频率;
在初步分析***为本地振荡模式的基础上,进一步判断电网功率信号的振荡频率fPe是否在本地振荡模式的频率范围内。由步骤(1)求得互联电网的各个振荡模式的振荡频率,对本地模式振荡频率进行排序,其中本地振荡模式频率最大值为fmax,最小值为fmin。若电网功率信号的振荡频率fPe在[fmin,fmax]范围以外,判断振荡产生机理为强迫振荡;若电网功率信号的振荡频率fPe在[fmin,fmax]范围以内,振荡产生机理须进一步判断。
fPe在[fmin,fmax]范围以外时,F1=0,振荡产生机理为强迫振荡;
fPe在[fmin,fmax]范围以内时,F1=1,振荡产生机理须进一步判断。
(c)判据F2:电网功率信号振荡波形;
当电网功率信号的振荡频率fPe在本地模式范围内时,通过振荡信号的波形特征判断振荡机理。若振荡波形的起振阶段为增幅振荡,阻尼呈现为负值,而稳态阶段为等幅振荡,阻尼呈现为弱阻尼,振荡产生机理为强迫振荡;若振荡波形的起振阶段和稳态阶段均为增幅振荡,阻尼一直呈现为负值,振荡产生机理为负阻尼振荡。
起振阶段为增幅振荡而稳态阶段为等幅振荡时,F2=0,振荡产生机理为强迫振荡;
起振阶段与稳态阶段均为增幅振荡时,F2=1,振荡产生机理为负阻尼振荡。
(d)判据F3:电网功率信号与调速***相关信号的振荡频率的一致性;
在判断***功率振荡产生机理为强迫振荡的基础上,利用电网功率信号与调速***信号的振荡频率的一致性来判断该强迫振荡是否由调速***引起。
方差是在概率论和统计方差衡量随机变量或一组数据时离散程度的度量,概率论中方差用来度量随机变量和其数学期望之间的偏离程度,方差的计算如式(17)所示
其中,μ=E(X)
令电网功率信号与各调速***信号的振荡频率为变量,则频率的方差如式(18)所示,用来度量电网功率的振荡频率与调速***信号的振荡频率的一致性。
其中,μ=E(f);
方差的大小与信号振荡频率的一致程度呈负相关性,即方差越小,则数据间的波动越小,信号间的振荡频率一致性越高;方差越大,则数据间的波动越大,信号间的振荡频率一致性越低。
D(f)≤0.005时,判据F3=1,判断该强迫振荡是由调速***引起;
D(f)>0.005时,判据F3=0,判断该强迫振荡不是由调速***引起。
南方电网某电厂共有2台额定功率330MW的汽轮发电机组。2013年5月8日,机组检修后并网运行,1号机组带220MW负荷运行,2号机组带230MW负荷运行,13:50分1号机组开始进行单顺阀切换操作。阀切换过程中,1号机组发生低频功率振荡,有功功率在186MW~279MW范围内摆动,振荡持续77秒,振荡频率为0.171HZ,PMU/WAMS***记录到了机组有功振荡,其实测信号如图2所示。
(a)实测信号的振荡频率提取;
图3—图7分别为调速***的调节阀门指令、等效阀位、阀位开度、调节级压力、发电机转速等实测信号图,通过Prony算法,可依次得出各实测信号的振荡频率。其中,fuT=0.163HZ,fcv=0.166HZ,fgv=0.163HZ,fPm=0.180HZ,fω=0.164HZ。
(b)电网功率信号振荡频率;
功率振荡的频率fPe为0.171HZ,电网本地振荡模式的最小振荡频率fmin为0.8HZ,最小振荡频率fmax为1.2HZ。由于fPe在[fmin,fmax]范围以外时,判据F1=0,振荡产生机理为强迫振荡。
(c)电网功率信号与调速***信号的振荡频率的一致性;
根据公式(18),对信号的振荡频率进行方差计算,求得实测信号振荡频率的期望值
E(f)=0.168HZ
求得实测信号振荡频率的方差
D(f)=3.714×10-5
由于D(f)≤0.005,判据F3=1,因此,可以判断本次低频振荡是由调速***引起的强迫振荡。
Claims (1)
1.一种调速***引发低频振荡的机理识别方法,其特征在于,所述方法包括
步骤1、计算互联电网的各振荡模式,根据特征值求各振荡模式的振荡频率与阻尼比;计算参与因子,参与因子Pki表示第i个模式中第k个变量的参与程度;
步骤2、实测信号,信号包括电网功率、发电机转速、调节阀门指令、等效阀位、阀位开度以及调节级压力;
步骤3、提取实测信号的振荡频率及其衰减因子;
步骤4、基于上述步骤和判据判断调速***低频振荡的机理;
所述步骤4的具体过程为:
步骤401、对***进行初步模态分析;
根据广域测量***同步相量测量装置的监测数据,若是监测到至少两个区域的机组群发生低频振荡现象,则可判断***振荡模式为区域振荡模式,振荡产生机理为负阻尼振荡;若是仅监测到单个发电厂内的机组发生低频振荡现象,则可初步判断***振荡模式为本地振荡模式,振荡产生机理须进一步判断;
步骤402、判据F1:电网功率信号振荡频率;
由步骤1求得互联电网本地振荡模式的振荡频率,其中频率最小值为fmin,最大值为fmax;当电网功率信号的振荡频率fPe在[fmin,fmax]范围以外时,判据F1=0,振荡产生机理为强迫振荡;当fPe在[fmin,fmax]范围以内时,判据F1=1,振荡产生机理须进一步判断;
步骤403、判据F2:电网功率信号振荡波形;
在F1=1的基础上,起振阶段为增幅振荡而稳态阶段为等幅振荡时,判据F2=0,振荡产生机理为强迫振荡;起振阶段与稳态阶段均为增幅振荡时,判据F2=1,振荡产生机理为负阻尼振荡;
步骤404、判据F3:电网功率信号与调速***信号的振荡频率的一致性;
计算电网实测信号与调速***实测信号的振荡频率方差:
其中,fPe为电网功率信号的振荡频率,为调节阀门指令的振荡频率,fω为发电机转速信号的振荡频率,fcv为等效阀位信号的振荡频率,fgv为阀位开度信号的振荡频率,fPm为调节级压力信号的振荡频率;μ为所有实测信号的平均振荡频率;
当D(f)≤0.005时,判据F3=1,判断该强迫振荡是由调速***引起;
当D(f)>0.005时,判据F3=0,判断该强迫振荡不是由调速***引起。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610515402.4A CN106099952B (zh) | 2016-07-01 | 2016-07-01 | 一种调速***引发低频振荡的机理识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610515402.4A CN106099952B (zh) | 2016-07-01 | 2016-07-01 | 一种调速***引发低频振荡的机理识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106099952A CN106099952A (zh) | 2016-11-09 |
CN106099952B true CN106099952B (zh) | 2018-07-31 |
Family
ID=57211780
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610515402.4A Active CN106099952B (zh) | 2016-07-01 | 2016-07-01 | 一种调速***引发低频振荡的机理识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106099952B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108267689B (zh) * | 2017-12-19 | 2020-07-14 | 中国神华能源股份有限公司 | 发电机组的调速***的阻尼极性判别方法、*** |
CN112751347B (zh) * | 2020-12-30 | 2023-05-09 | 国网河北省电力有限公司电力科学研究院 | 一种汽轮机调节***参与的低频振荡控制方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102055201A (zh) * | 2010-12-09 | 2011-05-11 | 北京四方继保自动化股份有限公司 | 基于微扰动信号振荡模式辨识的电力***低频振荡机理分析方法 |
CN102928695A (zh) * | 2012-10-18 | 2013-02-13 | 中国电力科学研究院 | 基于直线法判别负阻尼振荡与强迫振荡的方法 |
CN102944798A (zh) * | 2012-11-29 | 2013-02-27 | 武汉华中电力电网技术有限公司 | 一种负阻尼低频振荡与强迫功率振荡判别方法 |
CN104538952A (zh) * | 2014-12-08 | 2015-04-22 | 广东电网有限责任公司电力科学研究院 | 汽轮机调速***低频振荡抑制方法及装置 |
-
2016
- 2016-07-01 CN CN201610515402.4A patent/CN106099952B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102055201A (zh) * | 2010-12-09 | 2011-05-11 | 北京四方继保自动化股份有限公司 | 基于微扰动信号振荡模式辨识的电力***低频振荡机理分析方法 |
CN102928695A (zh) * | 2012-10-18 | 2013-02-13 | 中国电力科学研究院 | 基于直线法判别负阻尼振荡与强迫振荡的方法 |
CN102944798A (zh) * | 2012-11-29 | 2013-02-27 | 武汉华中电力电网技术有限公司 | 一种负阻尼低频振荡与强迫功率振荡判别方法 |
CN104538952A (zh) * | 2014-12-08 | 2015-04-22 | 广东电网有限责任公司电力科学研究院 | 汽轮机调速***低频振荡抑制方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN106099952A (zh) | 2016-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Avendaño-Valencia et al. | Stationary and non-stationary random vibration modelling and analysis for an operating wind turbine | |
CN106525415B (zh) | 一种风电机组传动链健康状态评价***及方法 | |
CN112886638B (zh) | 一种发电机惯量在线辨识方法及*** | |
CN112834224B (zh) | 一种核电汽轮发电机健康状态评估方法及*** | |
CN106099952B (zh) | 一种调速***引发低频振荡的机理识别方法 | |
CN104993501A (zh) | 一种励磁调节器抑制低频振荡性能的在线评估方法 | |
CN109753689A (zh) | 一种电力***机电振荡模态特征参数在线辩识方法 | |
CN104578115A (zh) | 一种基于相关函数的电力***低频振荡模式辨识方法 | |
CN110492497A (zh) | 一种发电机组功率振荡综合诊断方法及*** | |
CN104133129B (zh) | 汽轮发电机组调速***的低频振荡的监测方法和*** | |
CN108241761A (zh) | 发电机部件的疲劳损伤的确定方法和设备 | |
Seppänen et al. | Long Term Characteristics of Ultra Low Frequency Oscillations in the Nordic Power System | |
CN105158687A (zh) | 一种发电机次同步振荡模态信号的测量设备和方法 | |
Verdugo et al. | Testbed for PSS tuning using synchrophasor measurements and a real-time digital simulator | |
Fabian et al. | Generation of synthetic data with low-dimensional features for condition monitoring utilizing Generative Adversarial Networks | |
Gong et al. | Synchrophasor-based online modal analysis to mitigate power system interarea oscillation | |
Senesoulin et al. | Deep neural network estimation of inter-area oscillation mode based on synchrophasor data | |
CN108267689B (zh) | 发电机组的调速***的阻尼极性判别方法、*** | |
Shi et al. | Research on fault classification of wind turbine based on IMF kurtosis and PSO-SOM-LVQ | |
Pinheiro et al. | A procedure for the estimation of frequency response using a data-driven method | |
Zhou et al. | Ambient data-driven oscillation modes extraction using output-only observer/kalman filter identification and dynamic mode decomposition | |
Niu et al. | An Analytical Model for Frequency Nadir Prediction Based on Probabilistic Production Simulation | |
Han et al. | An approach for estimating mode shape for participation of inter-area oscillation mode | |
Wang et al. | Dominant inter-area oscillation mode identification using local measurement and modal energy for large-scale power systems with high grid-tied VSCs penetration | |
Prakash et al. | Novel approach for signal selection and optimal control in emerging wide area monitored systems |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |