CN114239218A - 一种变工况下汽轮机组主蒸汽压力寻优计算方法 - Google Patents

一种变工况下汽轮机组主蒸汽压力寻优计算方法 Download PDF

Info

Publication number
CN114239218A
CN114239218A CN202111321388.1A CN202111321388A CN114239218A CN 114239218 A CN114239218 A CN 114239218A CN 202111321388 A CN202111321388 A CN 202111321388A CN 114239218 A CN114239218 A CN 114239218A
Authority
CN
China
Prior art keywords
pressure
steam
working condition
main steam
extraction
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
Application number
CN202111321388.1A
Other languages
English (en)
Other versions
CN114239218B (zh
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.)
Zhongnan Electric Power Test and Research Institute of China Datang Group Science and Technology Research Institute Co Ltd
Original Assignee
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 Zhongnan Electric Power Test and Research Institute of China Datang Group Science and Technology Research Institute Co Ltd filed Critical Zhongnan Electric Power Test and Research Institute of China Datang Group Science and Technology Research Institute Co Ltd
Priority to CN202111321388.1A priority Critical patent/CN114239218B/zh
Publication of CN114239218A publication Critical patent/CN114239218A/zh
Application granted granted Critical
Publication of CN114239218B publication Critical patent/CN114239218B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Control Of Turbines (AREA)

Abstract

本发明提供一种变工况下汽轮机组主蒸汽压力寻优计算方法,包括以下步骤:设定基准工况,确定发电机功率,改变主蒸汽压力;假定综合阀位变化量,计算变化后主蒸汽流量、高压缸效率、给水泵有效功率与新的抽汽压力;计算低压缸排汽量、加热器出口水温、各级抽汽流量、主蒸汽吸热量、再热蒸汽吸热量与凝汽器放热量,得到新的发电机功率;将新的发电机功率与基准工况的发电机功率进行比较,判断是否结束转到下一步;将变工况计算得出的不同热耗率进行比较,得到最优主蒸汽压力。本发明方法无需多次重复的热力性能试验,减少了传统滑压优化试验过程中边界参数变化、运行方式改变、系统隔离等因素的影响,为后续汽轮机组滑压优化提供新思路。

Description

一种变工况下汽轮机组主蒸汽压力寻优计算方法
技术领域
本发明涉及一种变工况下汽轮机组主蒸汽压力寻优计算方法,属于火力发电厂汽轮机组运行优化领域。
背景技术
汽轮机滑压优化是指在电负荷一定时,改变主蒸汽压力导致蒸汽参数也随之变化,蒸汽参数变化时对应的调门开度也不同,则会有如下影响:1)部分进汽及节流损失发生变化,使调节级焓降偏离设计值,从而改变高压缸效率;2)若汽动给水泵出口压力下降,功耗降低,使拖动给水泵的小汽轮机耗汽减小,经济性提高;3)蒸汽参数降低会使汽轮机焓降减少,循环热效率降低,而机组的绝对内效率是相对内效率与循环热效率的乘积,故机组的绝对内效率降低。综合这些因素即可找出高压缸效率、给水泵耗功与循环热效率三者之间的最优组合,得出不同负荷下的最佳运行主蒸汽压力,提升汽轮机组的经济性能。
汽轮机热耗率是综合反映汽轮机组经济性的热力学指标,因此同一工况下对应汽轮机热耗率最低的主蒸汽压力即为最优主蒸汽压力。传统汽轮机滑压优化汽轮机最优主蒸汽压力是通过现场热力性能试验计算得到,通过比较相同负荷,不同主蒸汽压力下的汽轮机热耗率,选择汽轮机热耗率最低的主蒸汽压力作为最优主蒸汽压力,但试验期间,受运行方式改变、入炉煤质变化、给水流量测点精度、系统隔离等因素的影响,现场试验得到的热耗率误差较大,不能真实体现主蒸汽压力的变化对热耗率的影响,且机组的边界参数变化,也会对最优主蒸汽压力产生影响。
发明内容
针对现有技术的不足,本发明的目的是提供一种变工况下汽轮机组主蒸汽压力寻优计算方法,无需做多次的热力性能试验,便可直接计算得到滑压优化结论。。
为了实现上述目的,本发明所采用的技术方案是:
一种变工况下汽轮机组主蒸汽压力寻优计算方法,包括以下步骤:
步骤1:在汽轮机组确定负荷下,固定主蒸汽压力,进行汽轮机热力性能试验,按照试验标准,采集系统参数并计算汽轮机热耗率,确定发电机功率,以此次试验作为基准工况,然后改变主蒸汽压力;
步骤2:改变主蒸汽压力后假定综合阀位变化量,根据改变后的综合阀位变化量计算变化后主蒸汽流量、变化后高压缸效率、给水泵有效功率与新的抽汽压力;
步骤3:根据变工况后计算出来的抽汽压力,计算低压缸排汽量、加热器出口水温、各级抽汽流量、主蒸汽吸热量、再热蒸汽吸热量与凝汽器放热量,从而得到新的发电机功率;
步骤4:将步骤3中得到的发电机功率与步骤1中的基准工况的发电机功率进行比较,如果二者偏差绝对值除以基准工况的发电机功率的结果小于0.1%,即认为计算结束转到步骤5,否则返回到步骤2;
步骤5:将变工况计算得出的不同热耗率相互进行比较,热耗率最小工况下的主蒸汽压力即为最优主蒸汽压力。
所述步骤2中变化后主蒸汽流量、变化后高压缸效率、变化后给水泵有效功率与新的抽汽压力的计算由改变后的主蒸汽压力大小与假定的综合阀位进行计算得来。
所述变化后主蒸汽流量的计算方法如下:
根据弗留格尔公式:
Figure BDA0003345696950000021
式中:G0为基准工况主蒸汽流量,t/h;G1为变工况主蒸汽流量,t/h;p00为基准工况主蒸汽压力,MPa;p01为变工况主蒸汽压力,MPa;T00为基准工况调节级后温度,K;T01为变工况调节级后温度,K;A0为基准工况汽轮机进汽面积,m2;A1为变工况汽轮机进汽面积,m2
由于主蒸汽流量与调节级压力呈正比关系,于是用调节级压力代表主蒸汽流量,同时忽略温度的变化,对式(1)变换得到式(2);
Figure BDA0003345696950000022
式中:Ψ为综合阀位,%;K1为调节级压力与主蒸汽流量的线性系数,通过试验获取;p10为基准工况调节级压力,MPa;p11为变工况调节级压力,MPa;
当主蒸汽压力改变后,为维持发电机功率不变,综合阀位会发生改变,进而改变主蒸汽流量,对式(2)进一步化简得到式(3);
ΔG1=K2Δψp01 (3)
式中ΔG1:主蒸汽流量变化量,t/h;ΔΨ:综合阀位变化量,%;p01:变工况主蒸汽压力,MPa;K2:线性系数,由运行数据计算得到。
所述变化后高压缸效率计算方法如下:
针对试验机组进行现场试验,得出综合阀位与高压缸效率之间的关系曲线,拟合得到高压缸效率的计算公式:
ηgyg=aΨ+b (4)
式中ηgyg:高压缸效率,%;Ψ:综合阀位,%;a、b为线性系数,由现场试验拟合关系曲线得到。
所述给水泵有效功率计算方法如下:
首先计算得到给水泵出口压力,主蒸汽压力改变之后,利用给水管道压降,以基准工况给水泵出口压力加上工况改变后的给水压降得出变工况给水泵出口压力;给水管道的压降与给水流量近似呈正比关系,计算公式如下:
Figure BDA0003345696950000031
式中下标1:变工况,下标0:基准工况,Δp:给水压降,MPa;Dgs:给水流量,t/h;Kgs:压降系数,通过试验获取;
给水泵扬程计算公式如式(6):
Figure BDA0003345696950000032
给水泵有效功率计算公式如式(7):
Figure BDA0003345696950000033
给水泵效率计算公式如式(8):
Figure BDA0003345696950000034
小汽轮机输出功率计算公式如式(9):
Figure BDA0003345696950000035
小汽轮机进汽量计算公式如式(10):
Figure BDA0003345696950000036
泵组效率计算公式如式(11):
ηbz=ηxj·ηsb·ηlzq (11)
式(6)至式(11)中,H:给水泵扬程,m;pjk、pck:给水泵进、出口压力,MPa;ρ:给水泵平均密度,kg/m3;vjk、vck:给水泵进、出口流速,m/s;g:重力加速度,取9.81m/s2;ΔZ:给水泵进出口高度差,m;Dgs:给水流量,t/h;ηsb:给水泵效率,%;ηlzq:联轴器效率,%;ηxj:小汽轮机效率,%;ηbz:泵组效率,%;Pe:给水泵有效功率,kW;Pxj:小汽轮机输出功率,kW;Δhsb:给水泵焓升,kJ/kWh;Δhxj:小汽轮机理想焓降,kJ/kWh;Dxj:小汽轮机进汽量,t/h。
所述新的抽气压力计算方法如下:
对于中间各级组,在计算抽汽压力时,参照《ASME PTC6-2004汽轮机性能试验规程》中一类修正的处理方法,通过计算抽汽口后的蒸汽流量,再间接计算得到新的抽汽压力,计算公式如式(12);
Figure BDA0003345696950000041
式中Cjz:流量系数,Djz:抽汽口后的蒸汽流量,t/h;pcq:新的抽汽压力,MPa,vcq:抽汽比容,m3/kg;
级段压力变工况计算采用迭代方法进行,步骤如下:首先假设抽汽压力,然后根据变化后高压缸效率与抽汽蒸汽参数,求出抽汽焓、比容与流量,再对加热器进行热平衡计算,得到抽汽流量和抽汽口后蒸汽流量,然后利用式(12)计算得到新的抽汽压力,将新的抽汽压力与假设的抽汽压力进行对比,直至假设抽气压力与新的抽汽压力之差的绝对值除以假设的抽汽压力得到的结果小于1%,便结束迭代计算。
所述步骤3中根据变工况后计算出来的抽汽压力,根据汽轮机性能试验标准计算低压缸排汽量、加热器出口水温、各级抽汽流量、主蒸汽吸热量、再热蒸汽吸热量与凝汽器放热量,从而得到新的发电机功率,这一计算过程中为提高计算模型的准确性,模型中的其他参数设置如下:
(1)抽汽压损、加热器端差取基准工况值,变工况计算时与基准工况值相近,假定参数值不变;
(2)凝汽器的过冷度取基准工况值,变工况计算时保持不变;
(3)再热器压损取基准工况值,变工况计算时保持不变;
(4)由于中调阀维持全开,中压缸效率、低压缸ELEP效率取基准工况值,变工况计算时与基准工况值相近,假定参数值不变;
(5)低压缸排汽容积损失由变工况下的低压缸排汽流量计算得到;
(6)主蒸汽温度、再热蒸汽温度、排汽压力取基准工况值,变工况计算时保持不变。
所述步骤5中当需要计算某个负荷点的最优运行主蒸汽压力时,采用“穷举法”,在步骤1中以0.5MPa为步长改变主蒸汽压力,遍历所有压力区间,选择其中对应热耗率最小的主蒸汽压力作为该负荷点的最优主蒸汽压力。
本发明有益效果:
本发明提供一种变工况下汽轮机组主蒸汽压力寻优计算方法,可以准确计算出汽轮机组对应负荷下的最优主蒸汽压力,无需进行多次重复的热力性能试验,减少相对于传统滑压优化试验过程中边界参数变化、运行方式改变、系统隔离等因素的影响,为后续的汽轮机组滑压优化提供新的思路。
附图说明
图1为综合阀位与高缸效率的关系曲线图。
图2为给水压降与主给水流量的关系曲线图。
其中,Δp1/Δp0:变工况给水压降/基准工况给水压降,MPa;Dgs1/Dgs0:变工况给水流量/基准工况给水流量。
图3为最优主蒸汽压力变工况计算流程图。
具体实施方式
以下结合实施例对本发明的具体实施方式作进一步详细说明。
实施例
一种变工况下汽轮机组主蒸汽压力寻优计算方法,包括以下步骤:
步骤1:在机组确定负荷下,固定主蒸汽压力,进行汽轮机热力性能试验,按照试验标准,采集系统参数并计算汽轮机热耗率,确定发电机功率,以此次试验作为基准工况,然后改变主蒸汽压力。
步骤2:改变主蒸汽压力后假定综合阀位变化量,根据改变后的综合阀位变化量计算变化后主蒸汽流量、变化后高压缸效率、给水泵耗功与新的抽汽压力。
步骤3:根据变工况后计算出来的抽汽压力,计算低压缸排汽量、加热器出口水温、各级抽汽流量、主蒸汽吸热量、再热蒸汽吸热量与凝汽器放热量,从而得到新的发电机功率。
步骤4:将步骤3中得到的发电机功率与步骤1中的基准工况的发电机功率进行比较,如果二者偏差绝对值除以基准工况的发电机功率的结果小于0.1%,即认为计算结束转到步骤5,否则返回到步骤2。
步骤5:将变工况计算得出的不同热耗率相互进行比较,热耗率最小工况下的主蒸汽压力即为最优主蒸汽压力。
其中,步骤2中变化后主蒸汽流量、变化后高压缸效率、变化后给水泵耗功与新的抽汽压力指标的计算由改变后的主蒸汽压力大小与假定的综合阀位进行计算得来,
(1)变化后主蒸汽流量计算方法如下:
根据弗留格尔公式:
Figure BDA0003345696950000051
式中:G0为基准工况主蒸汽流量,t/h;G1为变工况主蒸汽流量,t/h;p00为基准工况主蒸汽压力,MPa;p01为变工况主蒸汽压力,MPa;T00为基准工况调节级后温度,K;T01为变工况调节级后温度,K;A0为基准工况汽轮机进汽面积,m2;A1为变工况汽轮机进汽面积,m2
对于高压缸来说,其进汽面积随高调门的开度变化而变化,而高调门开度则由综合阀位指令通过阀门管理函数生成,故可用综合阀位表征高压缸的进汽面积。
由于主蒸汽流量与调节级压力呈正比关系,可用调节级压力代表主蒸汽流量,同时忽略温度的变化,对式(1)稍作变换得到式(2)。
Figure BDA0003345696950000061
式中:Ψ为综合阀位,%;K1为调节级压力与主蒸汽流量的线性系数,通过试验获取;p10为基准工况调节级压力,MPa;p11为变工况调节级压力,MPa。
当主蒸汽压力改变后,为维持发电机功率不变,综合阀位会发生改变,进而改变主蒸汽流量,对式(2)进一步化简得到式(3)。
ΔG1=K2Δψp01 (3)
式中ΔG1:主蒸汽流量变化量,t/h,ΔΨ:综合阀位变化量,%,p01:变工况主蒸汽压力,MPa,K2:线性系数,可由运行数据计算得到。
(2)变化后高压缸效率计算方法如下:
针对试验机组进行现场试验,可得出综合阀位与高压缸效率之间的关系曲线,拟合得到公式:
ηgyg=aΨ+b (4)
式中ηgyg:高压缸效率,%;Ψ:综合阀位,%;a、b:线性系数,由现场试验拟合关系曲线得到。
需要注意的是,汽轮机高压缸效率与高调门的重叠度有关,若高调门的重叠度发生了变化,相应的函数关系也将发生变化。
(3)计算给水泵耗功需要提前计算得到给水泵出口压力。主蒸汽压力改变之后,可利用给水管道压降,以基准工况给水泵出口压力加上工况改变后的给水压降得出变工况给水泵出口压力。给水管道的压降与给水流量近似呈正比关系,计算公式如下:
Figure BDA0003345696950000062
式中下标1:变工况,下标0:基准工况,Δp:给水压降,MPa;Dgs:给水流量,t/h。Kgs:压降系数,通过试验获取。
给水泵耗功(即给水泵有效功率)计算方法如下:
给水泵扬程可由式(6)计算:
Figure BDA0003345696950000071
给水泵有效功率计算公式如式(7):
Figure BDA0003345696950000072
给水泵效率计算公式如式(8):
Figure BDA0003345696950000073
小汽轮机输出功率计算公式如式(9):
Figure BDA0003345696950000074
小汽轮机进汽量计算公式如式(10):
Figure BDA0003345696950000075
泵组效率计算公式如式(11):
ηbz=ηxj·ηsb·ηlzq (11)
式(6)至式(11)中,H:给水泵扬程,m;pjk、pck:给水泵进、出口压力,MPa;ρ:给水泵平均密度,kg/m3;vjk、vck:给水泵进、出口流速,m/s;g:重力加速度,取9.81m/s2;ΔZ:给水泵进出口高度差,m;Dgs:给水流量,t/h;ηsb:给水泵效率,%;ηlzq:联轴器效率,%;ηxj:小汽轮机效率,%;ηbz:泵组效率,%;Pe:给水泵有效功率,kW;Pxj:小汽轮机输出功率,kW;Δhsb:给水泵焓升,kJ/kWh;Δhxj:小汽轮机理想焓降,kJ/kWh;Dxj:小汽轮机进口流量(即小汽轮机进汽量),t/h。
(4)级段压力变工况(即新的抽气压力)计算方法如下:
对于中间各级组,在计算抽汽压力时,参照《ASME PTC6-2004汽轮机性能试验规程》中关于一类修正的处理方法,即通过计算抽汽口后的蒸汽流量,再间接计算得到新的抽汽压力,计算公式如式(12)。
Figure BDA0003345696950000076
式中Cjz:流量系数,Djz:抽汽口后的蒸汽流量,t/h;pcq:新的抽汽压力,MPa,vcq:抽汽比容,m3/kg。
Cjz实际上表征了汽轮机各级段的通流面积,其在各个负荷段几乎不变,可根据运行数据预先得到各级段的Cjz
级段压力变工况计算采用迭代方法进行,步骤如下:首先假设抽汽压力,然后根据变化后高压缸效率与抽汽蒸汽参数(抽汽压力、抽汽温度),求出抽汽焓、比容与流量,再对加热器进行热平衡计算,得到抽汽流量和抽汽口后蒸汽流量,然后利用式(12)计算得到新的抽汽压力,将新的抽汽压力与假设的抽汽压力进行对比,直至假设抽气压力与新的抽汽压力之差的绝对值除以假设的抽汽压力得到的结果小于1%,便结束迭代计算。
其中,步骤3中根据变工况后计算出来的抽汽压力,可根据汽轮机性能试验标准计算低压缸排汽量、加热器出口水温、各级抽汽流量、主蒸汽吸热量、再热蒸汽吸热量与凝汽器放热量,从而得到新的发电机功率,这一计算过程中为提高计算模型的准确性,模型中的其他参数设置如下:
(1)抽汽压损、加热器端差取基准工况值,变工况计算时与基准工况值近似,假定参数值不变;
(2)凝汽器的过冷度取基准工况值,变工况计算时保持不变;
(3)再热器压损取基准工况值,变工况计算时保持不变;
(4)由于中调阀维持全开,中压缸效率、低压缸ELEP效率取基准工况值,变工况计算时与基准工况值近似,假定参数值不变;
(5)低压缸排汽容积损失由变工况下的低压缸排汽流量计算得到;
(6)主蒸汽温度、再热蒸汽温度、排汽压力取基准工况值,变工况计算时保持不变。
其中,步骤5中当需要计算某个负荷点的最优运行主蒸汽压力时,采用“穷举法”,在步骤1中以0.5MPa为步长改变主蒸汽压力,遍历所有压力区间,选择其中对应热耗率最小的主蒸汽压力作为该负荷点的最优主蒸汽压力。
应用例
以实施例所述变工况下汽轮机组主蒸汽压力寻优计算方法,进行如下实验:
步骤1:以某超临界600MW机组为例,采用该机组投产后不同发电机功率下得到试验工况作为基准工况。然后选择要进行优化的工况,确定发电机功率,然后改变主蒸汽压力,在原工况压力基础上以0.5MPa为步长,增加或减少3MPa,共进行6次计算。
步骤2:改变主蒸汽压力后,假设综合阀位变化量,根据式(3)求出变化后主蒸汽流量;采用现场试验的方法,获取综合阀位与高压缸效率的变化关系如图1所示。拟合得到综合阀位与高压缸效率的函数关系为式(13),根据关系式求出变化后的高压缸效率。
ηgyg=0.2162ψ+69.5 (13)
式中ηgyg:高压缸效率,%;Ψ:综合阀位,%;
根据机组运行数据得出该机组的给水压降与主给水流量关系曲线如图2所示,由图可求得该机组的压降系数Kgs为1.139。再根据压降计算公式计算出给水泵的出口压力,各负荷下的给水流量和锅炉压降如表1所示,进一步根据式(6)至式(11)求出给水泵耗功、小汽轮机进汽量。小汽轮机对应的抽汽压力由接下来的级段压力变工况计算得出。
表1给水流量与给水压降关系
给水流量/t/h 1761 1179 838 704 531
给水压降/MPa 3.01 1.85 1.2 0.98 0.67
级段压力变工况计算采用迭代方法进行,步骤如下:首先假设抽汽压力,然后根据高压缸效率与抽汽蒸汽参数,求出抽汽焓、比容与流量,再对加热器进行热平衡计算,得到抽汽流量和抽汽口后蒸汽流量,然后利用式(12)计算得到新的抽汽压力,将新的抽汽压力与假设的抽汽压力进行对比,直至假设压力与新的抽汽压力之差的绝对值除以假设的抽汽压力得到的结果小于1%,便结束迭代计算得到变工况下的新抽汽压力。
步骤3:根据变工况后计算出来的抽汽压力,计算低压缸排汽量、加热器出口水温、各段抽汽流量、主蒸汽吸热量、再热蒸汽吸热量与凝汽器放热量,从而得到新的发电机功率。
步骤4:将新的发电机功率与原发电机功率对比,当二者偏差绝对值除以基准工况的发电机功率的结果小于0.1%时,即认为计算结束。否则返回到步骤2。
步骤5:当需要计算某个负荷点的最优运行主蒸汽压力时,采用“穷举法”,以0.5MPa为步长,遍历所有压力区间,选择其中对应热耗率最小的主蒸汽压力作为该负荷点的最优主蒸汽压力。变工况计算流程如图3所示,各个工况下的计算结果如表2所示:
表2某超临界600MW机组最优主蒸汽压力计算结果
Figure BDA0003345696950000091
由表2可知,该超临界机型采用变工况计算得到的最优初压相比设计值普遍要高,在滑压运行时应适当提高主蒸汽压力。

Claims (8)

1.一种变工况下汽轮机组主蒸汽压力寻优计算方法,其特征在于,包括以下步骤:
步骤1:在汽轮机组确定负荷下,固定主蒸汽压力,进行汽轮机热力性能试验,按照试验标准,采集系统参数并计算汽轮机热耗率,确定发电机功率,以此次试验作为基准工况,然后改变主蒸汽压力;
步骤2:改变主蒸汽压力后假定综合阀位变化量,根据改变后的综合阀位变化量计算变化后主蒸汽流量、变化后高压缸效率、给水泵有效功率与新的抽汽压力;
步骤3:根据变工况后计算出来的抽汽压力,计算低压缸排汽量、加热器出口水温、各级抽汽流量、主蒸汽吸热量、再热蒸汽吸热量与凝汽器放热量,从而得到新的发电机功率;
步骤4:将步骤3中得到的发电机功率与步骤1中的基准工况的发电机功率进行比较,如果二者偏差绝对值除以基准工况的发电机功率的结果小于0.1%,即认为计算结束转到步骤5,否则返回到步骤2;
步骤5:将变工况计算得出的不同热耗率相互进行比较,热耗率最小工况下的主蒸汽压力即为最优主蒸汽压力。
2.如权利要求1所述的方法,其特征在于,所述步骤2中变化后主蒸汽流量、变化后高压缸效率、变化后给水泵有效功率与新的抽汽压力的计算由改变后的主蒸汽压力大小与假定的综合阀位进行计算得来。
3.如权利要求2所述的方法,其特征在于,所述变化后主蒸汽流量的计算方法如下:
根据弗留格尔公式:
Figure FDA0003345696940000011
式中:G0为基准工况主蒸汽流量,t/h;G1为变工况主蒸汽流量,t/h;p00为基准工况主蒸汽压力,MPa;p01为变工况主蒸汽压力,MPa;T00为基准工况调节级后温度,K;T01为变工况调节级后温度,K;A0为基准工况汽轮机进汽面积,m2;A1为变工况汽轮机进汽面积,m2
由于主蒸汽流量与调节级压力呈正比关系,于是用调节级压力代表主蒸汽流量,同时忽略温度的变化,对式(1)变换得到式(2);
Figure FDA0003345696940000012
式中:Ψ为综合阀位,%;K1为调节级压力与主蒸汽流量的线性系数,通过试验获取;p10为基准工况调节级压力,MPa;p11为变工况调节级压力,MPa;
当主蒸汽压力改变后,为维持发电机功率不变,综合阀位会发生改变,进而改变主蒸汽流量,对式(2)进一步化简得到式(3);
ΔG1=K2Δψp01 (3)
式中ΔG1:主蒸汽流量变化量,t/h;ΔΨ:综合阀位变化量,%;p01:变工况主蒸汽压力,MPa;K2:线性系数,由运行数据计算得到。
4.如权利要求2所述的方法,其特征在于,所述变化后高压缸效率计算方法如下:
针对试验机组进行现场试验,得出综合阀位与高压缸效率之间的关系曲线,拟合得到高压缸效率的计算公式:
ηgyg=aΨ+b (4)
式中ηgyg:高压缸效率,%;Ψ:综合阀位,%;a、b为线性系数,由现场试验拟合关系曲线得到。
5.如权利要求2所述的方法,其特征在于,所述给水泵有效功率计算方法如下:
首先计算得到给水泵出口压力,主蒸汽压力改变之后,利用给水管道压降,以基准工况给水泵出口压力加上工况改变后的给水压降得出变工况给水泵出口压力;给水管道的压降与给水流量近似呈正比关系,计算公式如下:
Figure FDA0003345696940000021
式中下标1:变工况,下标0:基准工况,Δp:给水压降,MPa;Dgs:给水流量,t/h;Kgs:压降系数,通过试验获取;
给水泵扬程计算公式如式(6):
Figure FDA0003345696940000022
给水泵有效功率计算公式如式(7):
Figure FDA0003345696940000023
给水泵效率计算公式如式(8):
Figure FDA0003345696940000024
小汽轮机输出功率计算公式如式(9):
Figure FDA0003345696940000025
小汽轮机进汽量计算公式如式(10):
Figure FDA0003345696940000031
泵组效率计算公式如式(11):
ηbz=ηxj·ηsb·ηlzq (11)
式(6)至式(11)中,H:给水泵扬程,m;pjk、pck:给水泵进、出口压力,MPa;ρ:给水泵平均密度,kg/m3;vjk、vck:给水泵进、出口流速,m/s;g:重力加速度,取9.81m/s2;ΔZ:给水泵进出口高度差,m;Dgs:给水流量,t/h;ηsb:给水泵效率,%;ηlzq:联轴器效率,%;ηxj:小汽轮机效率,%;ηbz:泵组效率,%;Pe:给水泵有效功率,kW;Pxj:小汽轮机输出功率,kW;Δhsb:给水泵焓升,kJ/kWh;Δhxj:小汽轮机理想焓降,kJ/kWh;Dxj:小汽轮机进汽量,t/h。
6.如权利要求2所述的方法,其特征在于,所述新的抽气压力计算方法如下:
对于中间各级组,在计算抽汽压力时,参照《ASME PTC6-2004汽轮机性能试验规程》中一类修正的处理方法,通过计算抽汽口后的蒸汽流量,再间接计算得到新的抽汽压力,计算公式如式(12);
Figure FDA0003345696940000032
式中Cjz:流量系数,Djz:抽汽口后的蒸汽流量,t/h;pcq:新的抽汽压力,MPa,vcq:抽汽比容,m3/kg;
级段压力变工况计算采用迭代方法进行,步骤如下:首先假设抽汽压力,然后根据变化后高压缸效率与抽汽蒸汽参数,求出抽汽焓、比容与流量,再对加热器进行热平衡计算,得到抽汽流量和抽汽口后蒸汽流量,然后利用式(12)计算得到新的抽汽压力,将新的抽汽压力与假设的抽汽压力进行对比,直至假设抽气压力与新的抽汽压力之差的绝对值除以假设的抽汽压力得到的结果小于1%,便结束迭代计算。
7.如权利要求1所述的方法,其特征在于,所述步骤3中根据变工况后计算出来的抽汽压力,根据汽轮机性能试验标准计算低压缸排汽量、加热器出口水温、各级抽汽流量、主蒸汽吸热量、再热蒸汽吸热量与凝汽器放热量,从而得到新的发电机功率,这一计算过程中为提高计算模型的准确性,模型中的其他参数设置如下:
(1)抽汽压损、加热器端差取基准工况值,变工况计算时与基准工况值相近,假定参数值不变;
(2)凝汽器的过冷度取基准工况值,变工况计算时保持不变;
(3)再热器压损取基准工况值,变工况计算时保持不变;
(4)由于中调阀维持全开,中压缸效率、低压缸ELEP效率取基准工况值,变工况计算时与基准工况值相近,假定参数值不变;
(5)低压缸排汽容积损失由变工况下的低压缸排汽流量计算得到;
(6)主蒸汽温度、再热蒸汽温度、排汽压力取基准工况值,变工况计算时保持不变。
8.如权利要求1所述的方法,其特征在于,所述步骤5中当需要计算某个负荷点的最优运行主蒸汽压力时,采用“穷举法”,在步骤1中以0.5MPa为步长改变主蒸汽压力,遍历所有压力区间,选择其中对应热耗率最小的主蒸汽压力作为该负荷点的最优主蒸汽压力。
CN202111321388.1A 2021-11-09 2021-11-09 一种变工况下汽轮机组主蒸汽压力寻优计算方法 Active CN114239218B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111321388.1A CN114239218B (zh) 2021-11-09 2021-11-09 一种变工况下汽轮机组主蒸汽压力寻优计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111321388.1A CN114239218B (zh) 2021-11-09 2021-11-09 一种变工况下汽轮机组主蒸汽压力寻优计算方法

Publications (2)

Publication Number Publication Date
CN114239218A true CN114239218A (zh) 2022-03-25
CN114239218B CN114239218B (zh) 2023-02-28

Family

ID=80748807

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111321388.1A Active CN114239218B (zh) 2021-11-09 2021-11-09 一种变工况下汽轮机组主蒸汽压力寻优计算方法

Country Status (1)

Country Link
CN (1) CN114239218B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114943127A (zh) * 2022-06-28 2022-08-26 西安热工研究院有限公司 一种汽轮机顺序阀状态修正到单阀状态下的计算方法
CN115099029A (zh) * 2022-06-21 2022-09-23 上海电气集团股份有限公司 一种超临界火电机组缸效率变化对热耗率影响的计算方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104156882A (zh) * 2014-07-10 2014-11-19 华南理工大学 基于等效焓降及弗留格尔公式计算热力发电厂变工况方法
CN105201564A (zh) * 2015-08-26 2015-12-30 国网河南省电力公司电力科学研究院 一种基于主蒸汽流量的汽轮机滑压优化的控制方法
CN112444396A (zh) * 2020-11-11 2021-03-05 西安热工研究院有限公司 性能试验与全面变工况计算相结合的汽轮机滑压优化方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104156882A (zh) * 2014-07-10 2014-11-19 华南理工大学 基于等效焓降及弗留格尔公式计算热力发电厂变工况方法
CN105201564A (zh) * 2015-08-26 2015-12-30 国网河南省电力公司电力科学研究院 一种基于主蒸汽流量的汽轮机滑压优化的控制方法
CN112444396A (zh) * 2020-11-11 2021-03-05 西安热工研究院有限公司 性能试验与全面变工况计算相结合的汽轮机滑压优化方法

Non-Patent Citations (2)

* 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
CN115099029A (zh) * 2022-06-21 2022-09-23 上海电气集团股份有限公司 一种超临界火电机组缸效率变化对热耗率影响的计算方法
CN115099029B (zh) * 2022-06-21 2024-03-12 上海电气集团股份有限公司 一种超临界火电机组缸效率变化对热耗率影响的计算方法
CN114943127A (zh) * 2022-06-28 2022-08-26 西安热工研究院有限公司 一种汽轮机顺序阀状态修正到单阀状态下的计算方法
CN114943127B (zh) * 2022-06-28 2023-07-11 西安热工研究院有限公司 一种汽轮机顺序阀状态修正到单阀状态下的计算方法

Also Published As

Publication number Publication date
CN114239218B (zh) 2023-02-28

Similar Documents

Publication Publication Date Title
CN114239218B (zh) 一种变工况下汽轮机组主蒸汽压力寻优计算方法
CN105787211B (zh) 针对燃气透平劣化的联合循环余热锅炉压力调整方法
CN105201564B (zh) 一种基于主蒸汽流量的汽轮机滑压优化的控制方法
CN109298690B (zh) 一种开式循环水冷端系统优化控制方法
US20140290264A1 (en) Control of the gas composition in a gas turbine power plant with flue gas recirculation
CN109372594B (zh) 一种二次再热汽轮机滑压运行优化方法
CN111666675B (zh) 一种设计带低温省煤器的汽轮机热力性能考核试验修正计算方法
CN111652418B (zh) 一种火电机组动态精细化复合参数滑压曲线生成方法
CN113719325A (zh) 一种汽轮机变背压特性试验的改进方法
CN112444396B (zh) 性能试验与全面变工况计算相结合的汽轮机滑压优化方法
CN111666676B (zh) 一种低温省煤器系统节能量考核试验的修正计算方法
CN113107623A (zh) 一种双低压缸汽轮机低负荷运行时提升供热抽汽参数的装置及方法
US9470112B2 (en) System and method for heat recovery and steam generation in combined cycle systems
Kler et al. An effective approach to optimizing the parameters of complex thermal power plants
CN112127957B (zh) 一种火电厂汽轮机主蒸汽流量的测量方法
CN111927576A (zh) 一种汽轮发电机组空载能耗测试方法
CN111737859B (zh) 一种改进型汽轮机组变压运行耗差定量计算模型构建方法
Kler et al. Accounting for the variable nature of heat loads in optimizing cogeneration combined heat and power plants
KR20130034523A (ko) 과열기 온도 제어 방법
CN110298492B (zh) 一种凝汽器抽真空系统节能改造收益的评价方法
CN111507011B (zh) 供热抽汽对汽轮机滑压运行影响量的修正方法
WO2017096144A1 (en) Gas turbine firing temperature control with air injection system
CN114413251B (zh) 一种调节兀型对冲直流锅炉主汽温度的系统
CN112709610B (zh) 一种零号高压加热器调整抽汽控制曲线的优化方法
CN114961890B (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