CN113217119A - 一种汽轮机调速系统稳定性判断方法 - Google Patents

一种汽轮机调速系统稳定性判断方法 Download PDF

Info

Publication number
CN113217119A
CN113217119A CN202110511140.5A CN202110511140A CN113217119A CN 113217119 A CN113217119 A CN 113217119A CN 202110511140 A CN202110511140 A CN 202110511140A CN 113217119 A CN113217119 A CN 113217119A
Authority
CN
China
Prior art keywords
speed regulating
steam turbine
regulating
steam
regulating 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.)
Withdrawn
Application number
CN202110511140.5A
Other languages
English (en)
Inventor
李冰天
梁卓
邢志刚
杨森
杨旭辉
王龙川
杜新江
原峥
张许阳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Datang Sanmenxia Electric Power Co ltd
Zhongnan Electric Power Test and Research Institute of China Datang Group Science and Technology Research Institute Co Ltd
Original Assignee
Datang Sanmenxia Electric Power Co ltd
Zhongnan Electric Power Test and Research Institute of China Datang Group Science and Technology Research Institute Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Datang Sanmenxia Electric Power Co ltd, Zhongnan Electric Power Test and Research Institute of China Datang Group Science and Technology Research Institute Co Ltd filed Critical Datang Sanmenxia Electric Power Co ltd
Priority to CN202110511140.5A priority Critical patent/CN113217119A/zh
Publication of CN113217119A publication Critical patent/CN113217119A/zh
Withdrawn legal-status Critical Current

Links

Images

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F01MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
    • F01DNON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
    • F01D19/00Starting of machines or engines; Regulating, controlling, or safety means in connection therewith
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F01MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
    • F01DNON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
    • F01D17/00Regulating or controlling by varying flow
    • F01D17/10Final actuators
    • F01D17/105Final actuators by passing part of the fluid
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05DINDEXING SCHEME FOR ASPECTS RELATING TO NON-POSITIVE-DISPLACEMENT MACHINES OR ENGINES, GAS-TURBINES OR JET-PROPULSION PLANTS
    • F05D2260/00Function
    • F05D2260/80Diagnostics
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05DINDEXING SCHEME FOR ASPECTS RELATING TO NON-POSITIVE-DISPLACEMENT MACHINES OR ENGINES, GAS-TURBINES OR JET-PROPULSION PLANTS
    • F05D2260/00Function
    • F05D2260/81Modelling or simulation
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05DINDEXING SCHEME FOR ASPECTS RELATING TO NON-POSITIVE-DISPLACEMENT MACHINES OR ENGINES, GAS-TURBINES OR JET-PROPULSION PLANTS
    • F05D2270/00Control
    • F05D2270/70Type of control algorithm
    • F05D2270/706Type of control algorithm proportional-integral-differential
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05DINDEXING SCHEME FOR ASPECTS RELATING TO NON-POSITIVE-DISPLACEMENT MACHINES OR ENGINES, GAS-TURBINES OR JET-PROPULSION PLANTS
    • F05D2270/00Control
    • F05D2270/70Type of control algorithm
    • F05D2270/71Type of control algorithm synthesized, i.e. parameter computed by a mathematical model

Landscapes

  • Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Control Of Turbines (AREA)

Abstract

本发明涉及一种汽轮机调速系统稳定性判断方法,技术方案是通过用一个增益系数β进行替换,在汽轮机及调速系统经典模型的油动机惯性环节(Tydj)与高压缸前汽室环节(TCH)之间加入了高调门管理函数f(cv)与高调门流量特性g(cv)以及高调门管理函数线性因子β环节,其中f(cv)、g(cv)与β可以进行切换,最终通过计算调速系统阻尼比ζ、负阻尼临界值βf、过阻尼临界值βg,同时建立汽轮机调速系统是否处于稳定状态的判断标准,对汽轮机调速系统稳定进行判断,并给出不稳定的原因,其判断准确率高,大大减少了汽轮机调速系统因为摆动问题产生的不利影响,提高了系统一次调频的合格率,其使用方便,效果好,是汽轮机调速系统稳定性判断方法上的创新。

Description

一种汽轮机调速系统稳定性判断方法
技术领域
本发明涉及汽轮机调速系统稳定性判断,特别是一种考虑高调门管理函数的汽轮机调速系统稳定性判断方法。
背景技术
汽轮机调速系统的稳定性决定了汽轮机的调节品质。衡量调速系统的稳定性主要从一次调频性能、功率或转速闭环控制性能、配汽方式切换是否平稳、调节汽阀及EH油系统是否稳定等方面考虑。在导致调速系统稳定性下降的原因中,调速系统摆动问题最为突出。调速系统摆动不仅会导致调节汽阀、油动机等设备频繁动作,同时也有可能导致功率、转速出现大幅波动,并进一步诱发电力系统低频振荡事故。
汽轮机及其调速系统模型、参数实测与建模是电力系统稳定分析计算的基础之一,通过电力汽轮机及其调速系统模型,可以分析电力系统低频振荡、功率波动等故障。
如图1所示,为汽轮机及调速系统经典模型图,图1中DB—转速死区带;Kδ—速度不等率的倒数;K1—一次调频前馈系数;Tydj—油动机时间常数,s;LVDT—线性位移传感器;Pe—功率指令;Pm—实测功率;PGV—高调门开度;TCH为高压缸前汽室时间常数;TR为再热器时间常数;TIH为中压缸时间常数;FL为低压缸做功比例;FI为中压缸做功比例;FH为高压缸做功比例;Tcl为功率测量时间常数;Pg为输出功率;Pref为功率指令;p0为主蒸汽压力。
当转速w超越一次调频死区以后,经转速不等率折线函数,形成DEH功率指令,叠加在当前功率指令Pref上,然后与实际功率Pg对比后,计算二者偏差,然后经PID控制器,形成指令,驱动油动机动作(惯性环节Tydj),油动机动作后开启高调门,经高压缸前汽室(惯性环节TCH),进入高压缸做功,再经再热器(惯性环节TR),进入中压缸做功,最后经中低压缸连通管(惯性环节TIH),进入低压缸做功,高、中低压缸做功比例分别用FH、FI、FL代表,三者相加为1,最后形成实际输出功率,模型中的功率采用标幺值。
该模型忽略了汽轮机高调门管理函数这一环节,因高调门的开度与流量并非线性关系,在DEH逻辑中采用高调门管理函数进行修正,使综合阀位指令与蒸汽流量呈线性关系。但从实际运行情况来看,高调门管理函数常出现与高调门流量特性不匹配的情况,导致调速系统稳定性下降,进而导致机组一次调频能力差、功率波动、调速系统摆动时有发生。因此,经典的汽轮机及其调速系统模型未考虑高调门管理函数的影响,这会导致依据该模型得到的分析结果与实际情况不符。
因此,其改进和创新势在必行。
发明内容
针对上述情况,为克服现有技术之缺陷,本发明之目的就是提供一种汽轮机调速系统稳定性判断方法,可有效解决汽轮机调速系统稳定性分析及故障诊断的问题。
本发明解决的技术方案是:
一种考虑高调门管理函数的汽轮机调速系统稳定性判断方法,包括以下步骤:
步骤1:在电厂DCS系统(distributed control systems)中提取相关参数,导出机组负荷处于30%~100%之间、综合阀位指令处于50%~100%之间的历史数据,时间间隔为1s,具体包括以下参数:发电机功率、主蒸汽流量、主蒸汽压力、主蒸汽温度、调节级压力、调节级蒸汽温度、综合阀位指令、1号高调门指令、1号高调门开度、2号高调门指令、2号高调门开度、3号高调门指令、3号高调门开度、4号高调门指令、4号高调门开度;根据提取的参数,计算机组在30%~100%负荷的高调门管理函数线性因子β值,具体计算方法如下:
(1)高调门流量特性g(cv)
从调节级后至高压缸排汽末级,应用弗留格尔公式,得到下式:
Figure BDA0003060349370000021
式中:G0、G1分别为变工况前、后主蒸汽流量,kg/s;p1代表调节级压力,p10、p11分别为变工况前、后调节级压力,MPa;T10、T11分别为变工况前、后调节级蒸汽温度,K;A为通流面积,m2
主蒸汽压力为p0,采用p1/p0来表示高调门的流量,代表了经主蒸汽压力修正后的蒸汽流量G',G'再经过标幺化处理后即可得到高调门的流量特性公式为:
Figure BDA0003060349370000022
式中,G'x表示在x阀位处的高调门流量的标幺值,(p1/p0)1、(p1/p0)0、(p1/p0)x分别表示高调门在全开、全关和x开度下的流量;
(2)高调门管理函数线性因子β
高调门管理函数线性因子β代表汽轮机高调门管理函数的线性程度,它可用修正后蒸汽流量G'标幺值变化量与综合阀位变化量的比值代表,计算式为
Figure BDA0003060349370000023
式中ΔG'、ΔF分别为修正后流量标幺值的变化量和综合阀位指令的变化量,%;p01、p00分别为变工况、设计工况下的主蒸汽压力,MPa;G'xf1、G'xf2、G'xf3、G'xf4表示4个高调门在x阀位的标幺流量;Z12、Z3、Z4表示1/2号高调门、3号高调门、4号高调门的流量分配系数,计算式如下:
Figure BDA0003060349370000031
Figure BDA0003060349370000032
Z4=1-Z12-Z3
式中w1、w2、w3、w4表示1、2、3、4号高调门对应的调节级喷嘴数量;
步骤2:汽轮机机组在75%额定负荷运行,退出汽轮机功率闭环控制,退出一次调频功能,高调门单阀运行,在机组DEH逻辑中,强制综合阀位指令阶跃5%,采集相关参数,参数包括发电机功率、主蒸汽压力、调节级压力、高压缸排汽压力、再热压力、中压缸排汽压力、综合阀位指令、1号高调门指令、1号高调门开度,利用《同步发电机原动机及其调节系统参数实测与建模导则》中的汽轮机及调速系统经典模型,辨识得到汽轮机及其调速系统模型中的参数:TCH,TCH为高压缸前汽室时间常数;
步骤3:利用机组停机的机会,在机组DEH逻辑中,依次强制开启1号高调门、2号高调门、3号高调门、4号高调门,记录相关参数,参数包括:高调门开度、高调门指令、综合阀位指令,利用《同步发电机原动机及其调节系统参数实测与建模导则》中的汽轮机及调速系统经典模型,辨识得到汽轮机及其调速系统模型中的油动机时间常数Tydj
步骤4:计算得到机组的ζ、βf、βg
(1)调速系统阻尼比ζ
调速系统阻尼比ζ为:
Figure BDA0003060349370000033
(4)负阻尼临界值βf、过阻尼临界值βg
将ζ=0、1时对应的β值分别定义为负阻尼临界值和过阻尼临界值,用βf、βg表示,计算式分别为:
Figure BDA0003060349370000034
Figure BDA0003060349370000041
式中:Tyjd为油动机时间常数;TCH为高压缸前汽室时间常数;KP为DEH功率闭环PID控制器的比例系数;KI为DEH功率闭环PID控制器的积分系数;
步骤5:汽轮机调速系统稳定性判断
根据β、ζ、βf、βg的计算结果,判断汽轮机调速系统是否处于稳定状态,判断条件如下:
ζ>0,汽轮机调速系统处于稳定状态;
ζ≤0,汽轮机调速系统处于不稳定状态,调速系统将出现摆动;
β≥βf,汽轮机调速系统处于不稳定状态,调速系统将出现摆动;
βg≤β≤βf,汽轮机调速系统处于稳定状态,响应速度正常;
β≤βg,汽轮机调速系统处于稳定状态,但响应速度过慢。
当汽轮机调速系统处于非稳定状态时,调速系统将出现摆动的现象有调速系统摆动、功率振荡、转速波动、高调门摆动,分别有以下两种情况进行调整:
A、若ζ≤0,且βg≤β≤βf,这表明汽轮机调速系统处于不稳定状态,其原因是由于PID控制器中KI、KP参数设置不当造成的,需要重新整定PID参数;
B、若ζ≤0,且β≥βf,这表明汽轮机调速系统处于不稳定状态,其原因是由于高调门管理函数线性度差造成的,需要重新优化高调门管理函数,使β处于βg与βf之间。
本发明方法简单有效,理想情况下,高调门管理函数f(cv)与流量特性g(cv)二者是反函数的关系,而实际运行机组f(cv)与g(cv)非反函数关系,需要具体计算,由于高调门管理函数与高调门流量特性的修正,在汽轮机及其调速系统模型中主要是以增益形式体现的,因此本发明用一个增益系数β进行替换,在汽轮机及调速系统经典模型的油动机惯性环节(Tydj)与高压缸前汽室环节(TCH)之间加入了高调门管理函数f(cv)与高调门流量特性g(cv)以及高调门管理函数线性因子β环节,其中f(cv)、g(cv)与β可以进行切换,最终通过计算调速系统阻尼比ζ、负阻尼临界值βf、过阻尼临界值βg,同时建立汽轮机调速系统是否处于稳定状态的判断标准,对汽轮机调速系统稳定进行判断,并给出不稳定的原因,其判断准确率高,大大减少了汽轮机调速系统因为摆动问题产生的不利影响,提高了系统一次调频的合格率,其使用方便,效果好,是汽轮机调速系统稳定性判断方法上的创新。
附图说明
图1为本发明汽轮机及调速系统经典模型图。
图2为本发明高调门流量修正模型示意图。
图3为含高调门管理函数的汽轮机及调速系统模型图。
图4为简化后含高调门管理函数的汽轮机及其调速系统模型图。
图中w为实际转速;DB为一次调频死区带;Kδ为转速不等率倒数;PID为功率闭环控制器;Tydj为油动机时间常数;f(cv)为高调门管理函数;g(cv)为阀门流量特性;β为高调门管理函数线性因子;TCH为高压缸时间常数;TR为再热器时间常数;TIH为中压缸时间常数;FL为低压缸做功比例;FI为中压缸做功比例;FH为高压缸做功比例;Tcl为功率测量时间常数;Pg为输出功率;Pref为功率指令;p0为主蒸汽压力。
具体实施方式
以下结合附图对本发明的具体实施方式作进一步详细说明。
由图1-4给出,本发明一种考虑高调门管理函数的汽轮机调速系统稳定性判断方法,包括以下步骤:
步骤1:建立包含高调门管理函数的汽轮机及调速系统模型,即建立高调门流量修正模型,如图2所示,理想情况下,高调门管理函数f(cv)与流量特性g(cv)二者是反函数的关系,而实际运行机组f(cv)与g(cv)非反函数关系,需要具体计算,由于高调门管理函数与高调门流量特性的修正,在汽轮机及其调速系统模型中主要是以增益形式体现的,因此可用一个增益系数β,代替图3的模型,β定义为高调门管理函数线性因子,基于上述,建立如图4所示的包含高调门管理函数的汽轮机及其调速系统模型,该模型与经典汽轮机及其调速系统模型的主要区别是,在经典模型的油动机惯性环节(Tydj)与高压缸前汽室环节(TCH)之间加入了高调门管理函数f(cv)与高调门流量特性g(cv)及线性因子β环节,其中f(cv)、g(cv)与β可以切换。
在电厂DCS系统(distributed control systems)中提取相关参数,导出机组负荷处于30%~100%之间、综合阀位指令处于50%~100%之间的历史数据,时间间隔为1s,具体包括以下参数:发电机功率、主蒸汽流量、主蒸汽压力、主蒸汽温度、调节级压力、调节级蒸汽温度、综合阀位指令、1号高调门指令、1号高调门开度、2号高调门指令、2号高调门开度、3号高调门指令、3号高调门开度、4号高调门指令、4号高调门开度;根据提取的参数,计算机组在30%~100%负荷的高调门管理函数线性因子β值,具体计算方法如下:
(1)高调门流量特性g(cv)
从调节级后至高压缸排汽末级,应用弗留格尔公式,得到下式:
Figure BDA0003060349370000051
式中:G0、G1分别为变工况前、后主蒸汽流量,kg/s;p1代表调节级压力,p10、p11分别为变工况前、后调节级压力,MPa;T10、T11分别为变工况前、后调节级蒸汽温度,K;A为通流面积,m2
主蒸汽压力为p0,忽略蒸汽温度的变化,根据上式,蒸汽流量与调节级压力接近正比关系,因此可用调节级压力代表主蒸汽流量,调节级压力不仅受高调门开度影响,同时也受主蒸汽压力p0的影响;采用p1/p0来表示高调门的流量,代表了经主蒸汽压力修正后的蒸汽流量G',G'再经过标幺化处理后即可得到高调门的流量特性公式为:
Figure BDA0003060349370000061
式中,G'x表示在x阀位处的高调门流量的标幺值,(p1/p0)1、(p1/p0)0、(p1/p0)x分别表示高调门在全开、全关和x开度下的流量;
(2)高调门管理函数线性因子β
高调门管理函数线性因子β代表汽轮机高调门管理函数的线性程度,它可用修正后蒸汽流量G'标幺值变化量与综合阀位变化量的比值代表,计算式为
Figure BDA0003060349370000062
式中ΔG'、ΔF分别为修正后流量标幺值的变化量和综合阀位指令的变化量,%;p01、p00分别为变工况、设计工况下的主蒸汽压力,MPa;G'xf1、G'xf2、G'xf3、G'xf4表示4个高调门在x阀位的标幺流量;Z12、Z3、Z4表示1/2号高调门、3号高调门、4号高调门的流量分配系数,计算式如下:
Figure BDA0003060349370000063
Figure BDA0003060349370000064
Z4=1-Z12-Z3
式中w1、w2、w3、w4表示1、2、3、4号高调门对应的调节级喷嘴数量;可通过机组设计资料,得到调节级喷嘴数量w1、w2、w3、w4
汽轮机机组在75%额定负荷运行,退出汽轮机功率闭环控制,退出一次调频功能,高调门顺序阀运行,在机组DEH逻辑中,强制高调门动作,依次将1号高调门、2号高调门、3号高调门、4号高调门由0%开度,逐渐开启至100%开度,然后再由100%开度逐渐全关,采集相关参数:发电机功率、主蒸汽压力、调节级压力、综合阀位指令、1号高调门指令、1号高调门开度、2号高调门指令、2号高调门开度、3号高调门指令、3号高调门开度、4号高调门指令、4号高调门开度,依照步骤1中的高调门流量特性g(cv)的计算方法,就可以分别得到各高调门的流量特性;高调门的流量特性g(cv)与高调门管理函数f(cv)共同决定了β,当β超出正常范围后,可依据各高调门的流量特性和管理函数作进一步的分析,比如判断哪一个高调门造成β的异常,进行调整;
步骤2:汽轮机机组在75%额定负荷运行,退出汽轮机功率闭环控制,退出一次调频功能,高调门单阀运行,在机组DEH逻辑中,强制综合阀位指令阶跃5%,采集相关参数,参数包括发电机功率、主蒸汽压力、调节级压力、高压缸排汽压力、再热压力、中压缸排汽压力、综合阀位指令、1号高调门指令、1号高调门开度,利用《同步发电机原动机及其调节系统参数实测与建模导则》中的汽轮机及调速系统经典模型,辨识得到汽轮机及其调速系统模型中以下参数:TCH、TR、TIH、Tcl,其中,TCH为高压缸前汽室时间常数,TR为再热器时间常数,TIH为中压缸时间常数,Tcl为功率测量时间常数;TR、TIH、Tcl在调速系统稳定性判断中未涉及到,其主要作用是使调速系统模型参数更完整;
步骤3:利用机组停机的机会,在机组DEH逻辑中,依次强制开启1号高调门、2号高调门、3号高调门、4号高调门,记录相关参数,参数包括:高调门开度、高调门指令、综合阀位指令,利用《同步发电机原动机及其调节系统参数实测与建模导则》中的汽轮机及调速系统经典模型,辨识得到汽轮机及其调速系统模型中的油动机时间常数Tydj
查询机组DEH逻辑,得到调速系统PID控制启中的KI、KP,其中KI为DEH功率闭环PID控制器的积分系数,KP为DEH功率闭环PID控制器的比例系数;同时查询机组DEH逻辑,得到汽轮机高调门管理函数f(cv);根据该机组汽轮机热平衡图,通过热力性能计算得到FI、FH、FL;其中,FL为低压缸做功比例;FI为中压缸做功比例;FH为高压缸做功比例;
高调门管理函数f(cv)和高调门的流量特性g(cv)共同决定了β,当β超出正常范围后,可依据各高调门的流量特性和管理函数作进一步的分析,比如判断哪一个高调门造成β的异常,进行调整;FI、FH、FL是高、中、低压缸的做功比例,当调速系统稳定性下降后(比如功率出现了波动),依据做功比例,进一步分析高、中、低压缸分别对调速系统稳定性下降的权重,为电厂人员提供分析意见;所述步骤3、4均为本领域的常规技术手段,通过机组现有的系统即可得到步骤中的参数;
步骤4:计算得到机组的ζ、βf、βg
(1)调速系统阻尼比ζ
如图4所示,依据建立的简化版汽轮机及其调速系统模型,计算其闭环控制时的传递函数为:
Figure BDA0003060349370000081
根据劳斯判据,汽轮机及其调速系统稳定的判断条件为:
Figure BDA0003060349370000082
然后引入调速系统阻尼比ζ表示闭环控制下系统的阻尼特性,调速系统阻尼比ζ为:
Figure BDA0003060349370000083
(4)负阻尼临界值βf、过阻尼临界值βg
将ζ=0、1时对应的β值分别定义为负阻尼临界值和过阻尼临界值,用βf、βg表示,计算式分别为:
Figure BDA0003060349370000084
Figure BDA0003060349370000085
式中:Tyjd为油动机时间常数;TCH为高压缸前汽室时间常数;KP为DEH功率闭环PID控制器的比例系数;KI为DEH功率闭环PID控制器的积分系数;
步骤5:汽轮机调速系统稳定性判断
根据β、ζ、βf、βg的计算结果,判断汽轮机调速系统是否处于稳定状态,判断条件如下:
ζ>0,汽轮机调速系统处于稳定状态;
ζ≤0,汽轮机调速系统处于不稳定状态,调速系统将出现摆动;
β≥βf,汽轮机调速系统处于不稳定状态,调速系统将出现摆动;
βg≤β≤βf,汽轮机调速系统处于稳定状态,响应速度正常;
β≤βg,汽轮机调速系统处于稳定状态,但响应速度过慢。
当汽轮机调速系统处于非稳定状态时,调速系统将出现摆动的现象有调速系统摆动、功率振荡、转速波动、高调门摆动,分别有以下两种情况进行调整:
A、若ζ≤0,且βg≤β≤βf,这表明汽轮机调速系统处于不稳定状态,其原因是由于PID控制器中KI、KP参数设置不当造成的,需要重新整定PID参数;
B、若ζ≤0,且β≥βf,这表明汽轮机调速系统处于不稳定状态,其原因是由于高调门管理函数线性度差造成的,需要重新优化高调门管理函数,使β处于βg与βf之间。
本发明经实际应用均取得了相同或相近似的效果,具体应用例如下:
某机组一次调频考核合格率常低于50%且高负荷出现调速系统、功率及转速的持续性摆动问题,利用本发明方法,首先建立了该机组的汽轮机及其调速系统模型,然后通过现场测试、采集历史数据、查询设计资料等手段,计算得到该模型的相关参数。表1为线性因子β在综合阀位指令为50%~100%之间的计算结果。
表1线性因子β计算结果
Figure BDA0003060349370000091
该机组正常运行期间DEH功率闭环投入,KP=0.3,KI=10,Tydj=0.11s,TCH=0.32s。计算得到机组的βf=1.85,βg=0.096。由表1可见:该机组综合阀位在72%~75%区间时,β≈βg,导致机组进入过阻尼区域,响应速率下降,从而降低了一次调频调节能力,因此β值偏小是机组一次调频合格率较低的原因之一;而综合阀位在80%~90%区间时,β>βf,调速系统稳定性下降,故在一次调频动作或升降负荷时易出现功率波动现象。
随后对该机组进行高调门管理函数的重新优化,该机组上述问题得到妥善的处理,一次调频合格率明显提高,汽轮机调速系统未再出现摆动问题。优化后的线性因子β如表2。
表2线性因子β优化结果
Figure BDA0003060349370000092

Claims (2)

1.一种考虑高调门管理函数的汽轮机调速系统稳定性判断方法,其特征在于,包括以下步骤:
步骤1:在电厂DCS系统中提取相关参数,导出机组负荷处于30%~100%之间、综合阀位指令处于50%~100%之间的历史数据,时间间隔为1s,具体包括以下参数:发电机功率、主蒸汽流量、主蒸汽压力、主蒸汽温度、调节级压力、调节级蒸汽温度、综合阀位指令、1号高调门指令、1号高调门开度、2号高调门指令、2号高调门开度、3号高调门指令、3号高调门开度、4号高调门指令、4号高调门开度;根据提取的参数,计算机组在30%~100%负荷的高调门管理函数线性因子β值,具体计算方法如下:
(1)高调门流量特性g(cv)
从调节级后至高压缸排汽末级,应用弗留格尔公式,得到下式:
Figure FDA0003060349360000011
式中:G0、G1分别为变工况前、后主蒸汽流量,kg/s;p1代表调节级压力,p10、p11分别为变工况前、后调节级压力,MPa;T10、T11分别为变工况前、后调节级蒸汽温度,K;A为通流面积,m2
主蒸汽压力为p0,采用p1/p0来表示高调门的流量,代表了经主蒸汽压力修正后的蒸汽流量G',G'再经过标幺化处理后即可得到高调门的流量特性公式为:
Figure FDA0003060349360000012
式中,G'x表示在x阀位处的高调门流量的标幺值,(p1/p0)1、(p1/p0)0、(p1/p0)x分别表示高调门在全开、全关和x开度下的流量;
(2)高调门管理函数线性因子β
高调门管理函数线性因子β代表汽轮机高调门管理函数的线性程度,它可用修正后蒸汽流量G'标幺值变化量与综合阀位变化量的比值代表,计算式为
Figure FDA0003060349360000013
式中ΔG'、ΔF分别为修正后流量标幺值的变化量和综合阀位指令的变化量,%;p01、p00分别为变工况、设计工况下的主蒸汽压力,MPa;G'xf1、G'xf2、G'xf3、G'xf4表示4个高调门在x阀位的标幺流量;Z12、Z3、Z4表示1/2号高调门、3号高调门、4号高调门的流量分配系数,计算式如下:
Figure FDA0003060349360000021
Figure FDA0003060349360000022
Z4=1-Z12-Z3
式中w1、w2、w3、w4表示1、2、3、4号高调门对应的调节级喷嘴数量;
步骤2:汽轮机机组在75%额定负荷运行,退出汽轮机功率闭环控制,退出一次调频功能,高调门单阀运行,在机组DEH逻辑中,强制综合阀位指令阶跃5%,采集相关参数,参数包括发电机功率、主蒸汽压力、调节级压力、高压缸排汽压力、再热压力、中压缸排汽压力、综合阀位指令、1号高调门指令、1号高调门开度,利用《同步发电机原动机及其调节系统参数实测与建模导则》中的汽轮机及调速系统经典模型,辨识得到汽轮机及其调速系统模型中的参数TCH,TCH为高压缸前汽室时间常数
步骤3:利用机组停机的机会,在机组DEH逻辑中,依次强制开启1号高调门、2号高调门、3号高调门、4号高调门,记录相关参数,参数包括:高调门开度、高调门指令、综合阀位指令,利用《同步发电机原动机及其调节系统参数实测与建模导则》中的汽轮机及调速系统经典模型,辨识得到汽轮机及其调速系统模型中的油动机时间常数Tydj
步骤4:计算得到机组的ζ、βf、βg
(1)调速系统阻尼比ζ
调速系统阻尼比ζ为:
Figure FDA0003060349360000023
(4)负阻尼临界值βf、过阻尼临界值βg
将ζ=0、1时对应的β值分别定义为负阻尼临界值和过阻尼临界值,用βf、βg表示,计算式分别为:
Figure FDA0003060349360000024
Figure FDA0003060349360000025
式中:Tyjd为油动机时间常数;TCH为高压缸前汽室时间常数;KP为DEH功率闭环PID控制器的比例系数;KI为DEH功率闭环PID控制器的积分系数;
步骤5:汽轮机调速系统稳定性判断
根据β、ζ、βf、βg的计算结果,判断汽轮机调速系统是否处于稳定状态,判断条件如下:
ζ>0,汽轮机调速系统处于稳定状态;
ζ≤0,汽轮机调速系统处于不稳定状态,调速系统将出现摆动;
β≥βf,汽轮机调速系统处于不稳定状态,调速系统将出现摆动;
βg≤β≤βf,汽轮机调速系统处于稳定状态,响应速度正常;
β≤βg,汽轮机调速系统处于稳定状态,但响应速度过慢。
2.根据权利要求1所述的汽轮机调速系统稳定性判断方法,其特征在于,当汽轮机调速系统处于非稳定状态时,调速系统将出现摆动的现象有调速系统摆动、功率振荡、转速波动、高调门摆动,分别有以下两种情况进行调整:
A、若ζ≤0,且βg≤β≤βf,这表明汽轮机调速系统处于不稳定状态,其原因是由于PID控制器中KI、KP参数设置不当造成的,需要重新整定PID参数;
B、若ζ≤0,且β≥βf,这表明汽轮机调速系统处于不稳定状态,其原因是由于高调门管理函数线性度差造成的,需要重新优化高调门管理函数,使β处于βg与βf之间。
CN202110511140.5A 2021-05-11 2021-05-11 一种汽轮机调速系统稳定性判断方法 Withdrawn CN113217119A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110511140.5A CN113217119A (zh) 2021-05-11 2021-05-11 一种汽轮机调速系统稳定性判断方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110511140.5A CN113217119A (zh) 2021-05-11 2021-05-11 一种汽轮机调速系统稳定性判断方法

Publications (1)

Publication Number Publication Date
CN113217119A true CN113217119A (zh) 2021-08-06

Family

ID=77094671

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110511140.5A Withdrawn CN113217119A (zh) 2021-05-11 2021-05-11 一种汽轮机调速系统稳定性判断方法

Country Status (1)

Country Link
CN (1) CN113217119A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114398843A (zh) * 2022-01-18 2022-04-26 中国大唐集团科学技术研究院有限公司中南电力试验研究院 一种适用于多种地形的三维尾流风速分布计算方法
CN114430181A (zh) * 2022-02-25 2022-05-03 华北电力科学研究院有限责任公司 燃煤发电机组的调频仿真方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CH511365A (de) * 1968-04-19 1971-08-15 Westinghouse Electric Corp Regelvorrichtung für Dampfturbinen
CN103382860A (zh) * 2013-05-24 2013-11-06 华电国际电力股份有限公司山东分公司 汽轮机发电供热系统控制方法
CN109578088A (zh) * 2018-11-16 2019-04-05 华南理工大学 一种汽轮机dcs系统下的deh控制系统及控制方法
CN109888797A (zh) * 2019-03-09 2019-06-14 四川大学 基于态势感知的含风电接入送端系统频率控制方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CH511365A (de) * 1968-04-19 1971-08-15 Westinghouse Electric Corp Regelvorrichtung für Dampfturbinen
CN103382860A (zh) * 2013-05-24 2013-11-06 华电国际电力股份有限公司山东分公司 汽轮机发电供热系统控制方法
CN109578088A (zh) * 2018-11-16 2019-04-05 华南理工大学 一种汽轮机dcs系统下的deh控制系统及控制方法
CN109888797A (zh) * 2019-03-09 2019-06-14 四川大学 基于态势感知的含风电接入送端系统频率控制方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李冰天等: "考虑阀门管理函数的汽轮机调速系统", 《热力发电》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114398843A (zh) * 2022-01-18 2022-04-26 中国大唐集团科学技术研究院有限公司中南电力试验研究院 一种适用于多种地形的三维尾流风速分布计算方法
CN114398843B (zh) * 2022-01-18 2024-05-31 中国大唐集团科学技术研究院有限公司中南电力试验研究院 一种适用于多种地形的三维尾流风速分布计算方法
CN114430181A (zh) * 2022-02-25 2022-05-03 华北电力科学研究院有限责任公司 燃煤发电机组的调频仿真方法及系统
CN114430181B (zh) * 2022-02-25 2024-03-29 华北电力科学研究院有限责任公司 燃煤发电机组的调频仿真方法及系统

Similar Documents

Publication Publication Date Title
CN108708775B (zh) 基于热电负荷转换模型的热电联产机组快速变负荷控制方法
CN113217119A (zh) 一种汽轮机调速系统稳定性判断方法
CN110162870B (zh) 一种基于季节的节流调节汽轮机最优滑压曲线确定方法
CN109638861B (zh) 一种超临界机组参与一次调频的控制方法及控制系统模型
CN110242370B (zh) 一种超临界中间再热式机组滑压运行时参与电网一次调频的控制方法及控制系统模型
CN108757059B (zh) 一种汽轮机在节流配汽下的定滑压曲线试验优化方法
CN110792482A (zh) 一种超超临界二次再热机组参与电网一次调频的控制系统及方法
CN112650169B (zh) 基于焓值及燃料在线热值计算的发电机组主参数控制系统
WO2023087906A1 (zh) 一种基于事件触发的模糊神经网络温度控制系统及方法
CN109882833B (zh) 一种二次再热火电机组变负荷过程蒸汽温度控制方法
CN109378833B (zh) 一种通过控制汽轮机抽汽量实现机组快速调频的方法
CN113483319A (zh) 一种火电机组凝结水泵变频节能优化控制系统及方法
CN112627923B (zh) 极端工况下基于阀门特性曲线的汽轮机转速控制方法
CN108361683B (zh) 一种全负荷段再热气温智能控制系统
CN111765447B (zh) 一种基于多变量解耦的发电锅炉主汽温控制方法和系统
CN111045321B (zh) 一种深度调峰下的协调控制带嵌入式内模控制器的方法
CN113700533B (zh) 一种双抽汽轮机解耦调节系统及其控制方法
CN110579968A (zh) 一种超超临界机组深度调峰协调系统预测控制策略
CN106765520A (zh) 一种实现供热机组最优初压运行的自动控制方法
CN116241441A (zh) 锅炉给水系统汽动给水泵自动并泵控制方法
CN109802416B (zh) 提高汽轮发电机组deh一次调频性能的方法
CN113219242B (zh) 一种百万机组额定负荷工况下的一次调频试验方法
CN113110316B (zh) 一种联合循环机组汽机一次调频控制方法
CN113898434A (zh) 阶梯档位负载箱式二氧化碳发电试验平台透平调节方法
CN113606001A (zh) 一种600mw超临界机组旁路系统及其旁路温度控制方法

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
WW01 Invention patent application withdrawn after publication
WW01 Invention patent application withdrawn after publication

Application publication date: 20210806