一种汽轮机排汽焓的检测方法转让专利

申请号 : CN201210245091.6

文献号 : CN102749156B

文献日 :

基本信息:

PDF:

法律信息:

相似专利:

发明人 : 胥建群张赟杨涛周克毅石永锋

申请人 : 东南大学

摘要 :

本发明公开了一种汽轮机排汽焓的检测方法,首先利用设计数据或实验测量数据测算得到某个特定参数在变工况时的变化规律,然后计算得到上述特定参数在某一工况下的数值,并利用测量仪表测得该工况下低压缸的部分热力数据,再基于级内损失理论做出合理假设,得到测算汽轮机排汽焓的数值模型,最后将上述特定参数的数值和热力数据代入排汽焓测算模型,即可计算得到该工况下汽轮机的排汽焓值。本发明的方法以级内损失理论为基础,将汽轮机级内各类损失划分为与容积流量相关的排汽损失、与湿度相关的湿汽损失和与理想焓降相关的其他损失,计算过程简单清晰、精度高;无需知道低压缸具体结构参数,所需测点少,测量误差对计算结果造成的影响较小。

权利要求 :

1.一种汽轮机排汽焓的检测方法,其特征在于,该方法包括以下步骤:

1)分别在至少四种不同功率运行工况下测算得到汽轮机的低压缸的理想焓降Δht及低压缸喷嘴损失、动叶损失、叶高损失、扇形损失、叶轮摩擦损失、漏汽损失之和占低压缸理想焓降的百分比b,得到各自工况下的Δht和b,所述测算Δht值和b值的方法如下:

11)利用温度压力测量仪表测量得到汽轮机低压缸进口工质压力p0和低压缸进口工质温度t0,根据p0和t0,通过IFC-67工业用水和水蒸汽热力性质模型计算低压缸进口工质焓h0和低压缸进口工质熵s0;

利用压力测量仪表测量得到汽轮机低压缸排汽压力pc,根据IFC-67工业用水和水蒸汽热力性质模型,通过s0和pc计算工质自低压缸进口等熵膨胀的理想排汽焓hct;将h0减去hct,得到低压缸的理想焓降Δht;

根据IFC-67工业用水和水蒸汽热力性质模型,通过pc计算排汽压力下的饱和蒸汽焓hg和排汽压力下的饱和水焓hl;

根据设计数据查得汽轮机排汽流量Gc和排汽焓hcu;

根据IFC-67工业用水和水蒸汽热力性质模型,通过hcu和pc计算排汽比容v;

利用汽轮机的排汽损失-排汽容积流量图,根据Gc和v查得排汽损失Δhex;将设计排汽焓hcu减去排汽损失Δhex,得到未考虑排汽损失时的排汽焓hce;

根据IFC-67工业用水和水蒸汽热力性质模型,通过hce和pc计算得到未考虑排汽损失时的排汽熵sce;

l

12)根据psq 通过IFC-67工业用水和水蒸汽热力性质模型,计算得到本次迭代Ⅰ中的l工质进入湿蒸汽区时的饱和蒸汽焓值hsq,然后代入公式 计算得l

到本次迭代Ⅰ中的工质进入湿蒸汽区时的饱和蒸汽熵ssq ;所述l为迭代Ⅰ次数,没有迭代

0 l

时l=0,psq 为步骤11)测算得到的pc,l大于0时,psq 为上一次迭代Ⅰ中得到的工质进入湿蒸汽区时的饱和压力;

l l

13)根据所述步骤12)中求得的hsq 和ssq,通过IFC-67工业用水和水蒸汽热力性质模l+1型求得工质进入湿蒸汽区时的饱和压力psq ;

l+1 l

判断是否满足|psq -psq|≤r,如是则进入步骤14),否则令l=l+1并回到步骤12);所-3述r为工程应用的精度要求,0≤r≤10 ;

l+1

14)根据步骤13)求得的满足精度要求的psq ,通过IFC-67工业用水和水蒸汽热力性质模型,计算得到工质进入湿蒸汽区时的饱和蒸汽焓值hsq;

15)将步骤14)求得的hsq代入公式 求解湿蒸汽级的理想焓降在低压缸理想焓降中所占的百分比a;

16)根据公式 求解b值;

2)根据步骤1)中得到的不同功率运行工况下的Δht值和b值,利用最小二乘方法对b值进行二次拟合,得到第一常数t1,第二常数t2和第三常数t3,从而得到b值与Δht值的2

函数关系为b=t1Δht+t2Δht+t3;

3)求取任意功率运行工况下的排汽焓hcu,具体步骤为:

31)按照步骤1)中所述方法求得任意功率运行工况下的p0、t0、h0、s0、pc、hct、Δht、hg和hl;

2

32)将所述步骤31)中求得的Δht代入b=t1Δht+t2Δht+t3中,求得当前工况下的b值;

k

33)将所述步骤1)中的任一功率运行工况作为基准工况,然后根据hce 和所述步骤31)中求得的pc,通过IFC-67工业用水和水蒸汽热力性质模型,计算得到本次迭代Ⅱ中的未考k虑排汽损失时的排汽熵sce,然后按照所述步骤12)至15)的方法,计算得到本次迭代Ⅱ中k的当前工况下的湿蒸汽级的理想焓降在低压缸理想焓降中所占的百分比a ;所述k为迭代

0 k

Ⅱ次数,没有迭代时k=0,hce 为步骤1)在基准工况下的求得的hce,k大于0时,hce 为上一次迭代Ⅱ中得到的当前工况下的未考虑排汽损失时的排汽焓;

34)将所述步骤31)中求得的hg、hl、h0、hct,步骤32)中求得的b值,以及步骤33)中k求得的a 值代入公式 求解得到当前工

k+1

况下的未考虑排汽损失时的排汽焓hce ;

k+1 k

判断是否满足|hce -hce|≤e,如是则进入步骤35),否则令k=k+1并回到步骤33);所-3述e为工程应用的精度要求,0≤e≤10 ;

35)将步骤11)中得到的基准工况下的排汽流量Gc、基准工况下的低压缸进口工质压力p0分别用Gc0和p00来表示,通过公式 求解当前工况下的排汽流量Gd,式中,p0为所述步骤31)中求得的当前工况下的低压缸进口工质压力;

k+1

36)根据IFC-67工业用水和水蒸汽热力性质模型,通过步骤34)中得到的hce 和步骤31)中测得的pc计算得到当前工况下的排汽比容v;利用汽轮机的排汽损失-排汽容积流量图,根据步骤35)中得到的Gd和当前工况下的排汽比容v查得当前工况下的排汽损失k+1Δhex,将当前工况下的排汽损失Δhex加上步骤34)所得未考虑排汽损失的排汽焓值hce ,即可求得当前工况下汽轮机的实际排汽焓值hcu。

说明书 :

一种汽轮机排汽焓的检测方法

技术领域

[0001] 本发明涉及汽轮机热力系统分析领域,尤其涉及汽轮机的排汽焓值测算技术。

背景技术

[0002] 汽轮机作为常规火电站中的关键设备,其热力性能的优劣,不仅影响发电企业的经济效益,而且直接影响全社会的能源消耗及环境状况。因此,热力发电企业需要定期进行汽轮机组热力性能试验,或者借助先进的检测技术和计算机技术在线监测机组的热力性能,来了解和掌握汽轮机组各主要设备或部件的性能状态,评估机组的热力性能、经济与安全指标。
[0003] 汽轮机组的热经济指标主要有汽耗率、热耗率、相对内效率、循环热效率、绝对电效率等。此外,为了评价和鉴定机组的安装水平、调试质量、维护、技术改造,以及新技术的推广应用,经常需要进行热力试验。无论是汽轮机组的热力性能分析,还是试验结果的计算、修正,都存在低压缸排汽焓计算这一问题。
[0004] 汽轮机在运行过程中,工质膨胀到低压缸末级或末几级时可能会进入湿蒸汽区,而处于湿蒸汽区的蒸汽焓值不能够有压力和温度查得,低压缸的效率亦无法准确计算。因此,汽轮机排汽焓的计算一直是汽轮机热力计算的难点之一。深入研究这一问题,有助于研究人员从理论上进一步了解汽轮机系统,同时能够为现场运行试验人员掌握通流部分的工作状态,获取准确的汽轮机热力性能参数提供帮助,为汽轮机组的运行优化、性能试验与技术革新等工作提供参考依据。
[0005] 目前,国内外涉及排汽焓计算的研究有很多,这些方法或计算过程复杂、所需测点过多;或依赖流量测量仪的精度,计算结果受排汽流量精度影响大;或计算结果精度不高。为了得到较为精确的计算结果,汽轮机排汽焓的计算方法应满足以下几点要求:基本理论清晰,计算过程简单,所需测点较少。本发明基于级内损失理论,提出了一种新的排汽焓值测算方法;与其他方法相比,该方法具有不需要了解结构参数、所需测点少、测量误差对计算结果造成的影响较小、计算精度高、适用于在线测算等优点,便于现场实际应用。
[0006] 常用的排汽焓计算方法主要是通过级内详细计算和整体热平衡两个方面来计算的。
[0007] 级内详细计算方法需要知道级的结构参数。机组运行一段时间后,级的结构参数可能发生变化,此时若仍用设计结构数据,计算结果将会产生误差。
[0008] 整体热平衡计算方法所需要的测点可能较多。大型汽轮机组回热系统中最后两个低压加热器通常位于凝汽器内,它们之间的凝结水管道通常没有温度测点;实际机组循环水管道直径比较大,测量其参数时需要一段足够长的直管段,但现场往往不能满足此要求,测量误差较大。这些使得系统热平衡计算结果精度不够。

发明内容

[0009] 技术问题:本发明的目的在于提供一种测算模型简单、计算精度高、测量成本低的汽轮机排汽焓值测算方法。
[0010] 技术方案:本发明的汽轮机排汽焓的检测方法,包括以下步骤:
[0011] 1)分别在至少四种不同功率运行工况下测算得到汽轮机的低压缸的理想焓降Δht及低压缸喷嘴损失、动叶损失、叶高损失、扇形损失、叶轮摩擦损失、漏汽损失之和占低压缸理想焓降的百分比b,得到各自工况下的Δht和b,所述测算Δht值和b值的方法如下:
[0012] 11)利用温度压力测量仪表测量得到汽轮机低压缸进口工质压力p0和低压缸进口工质温度t0,根据p0和t0,通过国际公式化委员会(IFC)提出的IFC-67工业用水和水蒸汽热力性质模型计算低压缸进口工质焓h0和低压缸进口工质熵s0;
[0013] 利用压力测量仪表测量得到汽轮机低压缸排汽压力pc,
[0014] 根据IFC-67工业用水和水蒸汽热力性质模型,通过s0和pc计算工质自低压缸进口等熵膨胀的理想排汽焓hct;将h0减去hct,得到低压缸的理想焓降Δht;
[0015] 根据IFC-67工业用水和水蒸汽热力性质模型,通过pc计算排汽压力下的饱和蒸汽焓hg和排汽压力下的饱和水焓hl;
[0016] 根据设计数据查得汽轮机排汽流量Gc和排汽焓hcu;
[0017] 根据IFC-67工业用水和水蒸汽热力性质模型,通过hcu和pc计算排汽比容v;
[0018] 利用汽轮机的排汽损失-排汽容积流量图,根据Gc和v查得排汽损失Δhex;将设计排汽焓hcu减去排汽损失Δhex,得到未考虑排汽损失时的排汽焓hce;
[0019] 根据IFC-67工业用水和水蒸汽热力性质模型,通过hce和pc计算得到未考虑排汽损失时的排汽熵sce;
[0020] 12)根据psql通过IFC-67工业用水和水蒸汽热力性质模型,计算得到本次迭代Ⅰ中的工质进入湿蒸汽区时的饱和蒸汽焓值hsql,然后代入公式 计l
算得到本次迭代Ⅰ中的工质进入湿蒸汽区时的饱和蒸汽熵ssq ;所述l为迭代Ⅰ次数,没有
0 l
迭代时l=0,psq 为步骤11)测算得到的pc,l大于0时,psq 为上一次迭代Ⅰ中得到的工质进入湿蒸汽区时的饱和压力,;
[0021] 其中,公式 的推导过程为:
[0022] 根据现场实际运行情况,做出合理假设:近似认为低压缸中的蒸汽膨胀线为一条直线,得到关系式 式中:hsq为工质进入湿蒸汽区时的饱和蒸汽焓值,ssq为工质进入湿蒸汽区时的饱和蒸汽熵;然后根据该关系式推导得到上述公式。
l l
[0023] 13)根据所述步骤12)中求得的hsq 和ssq,通过IFC-67工业用水和水蒸汽热力性l+1质模型求得工质进入湿蒸汽区时的饱和压力psq ;
l+1 l
[0024] 判断是否满足|psq -psq|≤r,如是则进入步骤14),否则令l=l+1并回到步骤-312);所述r为工程应用的精度要求,0≤r≤10 ;
l+1
[0025] 14)根据步骤13)求得的满足精度要求的psq ,通过IFC-67工业用水和水蒸汽热力性质模型,计算得到工质进入湿蒸汽区时的饱和蒸汽焓值hsq。
[0026] 15)将步骤14)求得的hsq代入公式 求解湿蒸汽级的理想焓降在低压缸理想焓降中所占的百分比a;
[0027] 16)根据公式 求解b值;
[0028] 其中,公式 的推导过程为:
[0029] 根据现场实际运行情况,做出合理假设:当级组结构参数、压比、反动度保持不变时,喷嘴损失、动叶损失、叶高损失、扇形损失、叶轮摩擦损失、漏汽损失占理想焓降的比例不变;在工况变化不大时,撞击损失可忽略。
[0030] 将湿汽损失与排汽损失、喷嘴损失、动叶损失、叶高损失、扇形损失、叶轮摩擦损失、漏汽损失分开考虑,得到未考虑排汽损失的排汽焓计算数值模型xc×hg+(1-xc)×hl=hce,式中:xc为排汽干度,可
以推导得到上述公式。
[0031] 2)根据步骤1)中得到的不同功率运行工况下的Δht值和b值,利用最小二乘方法对b值进行二次拟合,得到第一常数t1,第二常数t2和第三常数t3,从而得到b值与Δht2
值的函数关系为b=t1Δht+t2Δht+t3;
[0032] 3)求取任意功率运行工况下的排汽焓hcu,具体步骤为:
[0033] 31)按照步骤1)中所述方法求得任意功率运行工况下的p0、t0、h0、s0、pc、hct、Δht、hg和hl;
[0034] 32)将所述步骤31)中求得的Δht代入b=t1Δht2+t2Δht+t3中,求得当前工况下的b值;
[0035] 33)将所述步骤1)中的任一功率运行工况作为基准工况,然后根据hcek和所述步骤31)中求得的pc,通过IFC-67工业用水和水蒸汽热力性质模型,计算得到本次迭代Ⅱ中k的未考虑排汽损失时的排汽熵sce,然后按照步骤12)至15)的方法,计算得到本次迭代Ⅱk
中的当前工况下的湿蒸汽级的理想焓降在低压缸理想焓降中所占的百分比a ;k为迭代Ⅱ
0 k
次数,没有迭代时k=0,hce 为步骤1)在基准工况下的求得的hce,k大于0时,hce 为上一次迭代Ⅱ中得到的当前工况下的未考虑排汽损失时的排汽焓,;
[0036] 34) 根 据 步 骤 16) 中 的 排 汽 焓 计 算 数 值 模 型xc×hg+(1-xc)×hl=hce,推导得到未考虑排汽损
失的排汽焓计算公式 将步骤31)中求得
k
的hg、hl、h0、hct,步骤32)中求得的b值,以及步骤33)中求得的a 值代入上述公式,求解k+1
得到当前工况下的未考虑排汽损失时的排汽焓hce ;
[0037] 判断是否满足|hcek+1-hcek|≤e,如是则进入步骤35),否则令k=k+1并回到步骤-333);所述e为工程应用的精度要求,0≤e≤10 ;
[0038] 35)将步骤11)中得到的基准工况下的排汽流量Gc、基准工况下的低压缸进口工质压力p0分别用Gc0和p00来表示,通过公式 求解当前工况下的排汽流量Gd,式中,p0为所述步骤31)中求得的当前工况下的低压缸进口工质压力;
k+1
[0039] 36)根据IFC-67工业用水和水蒸汽热力性质模型,通过步骤34)中得到的hce和步骤31)中测得的pc计算得到当前工况下的排汽比容v;利用汽轮机的排汽损失-排汽容积流量图,根据步骤35)中得到的Gd和当前工况下的排汽比容v查得当前工况下的排汽损失Δhex,将当前工况下的排汽损失Δhex加上步骤34)所得未考虑排汽损失的排汽焓值k+1hce ,即可求得当前工况下汽轮机的实际排汽焓值hcu。
[0040] 有益效果:相比于传统的排汽焓计算方法,本发明提出的新的汽轮机排汽焓测算方法具有以下优点:
[0041] 1.本发明基于级内损失理论,将汽轮机级内各类损失划分为与容积流量相关的排汽损失、与湿度相关的湿汽损失和与理想焓降相关的其他损失。测算过程简单清晰,无需了解结构参数,便于现场实际应用。
[0042] 2.本发明所需测点少,测量成本低:仅需要测量低压缸进口压力、温度;最后一个抽汽口的压力;背压。同时使得测量误差对结果造成的影响较小。
[0043] 3.本发明在测算过程中,排汽流量仅用于排汽损失的计算,因此最后计算结果不依赖排汽流量的精度。进一步提高了测算结果的精度。
[0044] 4.本发明可用于现场对汽轮机排汽焓的在线测算,为准确获取汽轮机的实时热力性能提供帮助。

附图说明

[0045] 图1为实例中汽轮机低压缸热力过程示意图。
[0046] 图2为本发明的逻辑流程框图。

具体实施方式

[0047] 以下通过不同工况下的实例进一步说明本发明。
[0048] 实施例1:根据四种不同工况检测上汽某300MW汽轮机任意工况下的排汽焓值。
[0049] 1)分别在四种不同功率运行工况:额定工况(THA工况)、75%额定工况(75%THA工况)、60%额定工况(60%THA工况)和50%额定工况(50%THA工况)下测算得到汽轮机的低压缸的理想焓降Δht及低压缸喷嘴损失、动叶损失、叶高损失、扇形损失、叶轮摩擦损失、漏汽损失之和占低压缸理想焓降的百分比b,得到各自工况下的Δht和b,所述测算Δht值和b值的方法如下:
[0050] 11)利用温度压力测量仪表测量得到汽轮机低压缸进口工质压力p0和低压缸进口工质温度t0,根据p0和t0,通过国际公式化委员会(IFC)提出的IFC-67工业用水和水蒸汽热力性质模型计算低压缸进口工质焓h0和低压缸进口工质熵s0;
[0051] 利用压力测量仪表测量得到汽轮机低压缸排汽压力pc,
[0052] 根据IFC-67工业用水和水蒸汽热力性质模型,通过s0和pc计算工质自低压缸进口等熵膨胀的理想排汽焓hct;将h0减去hct,得到低压缸的理想焓降Δht;
[0053] 根据IFC-67工业用水和水蒸汽热力性质模型,通过pc计算排汽压力下的饱和蒸汽焓hg和排汽压力下的饱和水焓hl;
[0054] 根据设计数据查得汽轮机排汽流量Gc和排汽焓hcu;
[0055] 根据IFC-67工业用水和水蒸汽热力性质模型,通过hcu和pc计算排汽比容v;
[0056] 利用汽轮机的排汽损失-排汽容积流量图,根据Gc和v查得排汽损失Δhex;将设计排汽焓hcu减去排汽损失Δhex,得到未考虑排汽损失时的排汽焓hce;
[0057] 根据IFC-67工业用水和水蒸汽热力性质模型,通过hce和pc计算得到未考虑排汽损失时的排汽熵sce;
[0058] 上述四种不同工况下的各热力参数如表1所示。
[0059] 表1四种不同工况下测算得到的各热力参数
[0060]
[0061] 12)根据psql通过IFC-67工业用水和水蒸汽热力性质模型,计算得到本次迭代Ⅰl中的工质进入湿蒸汽区时的饱和蒸汽焓值hsq,然后代入公式 计
l
算得到本次迭代Ⅰ中的工质进入湿蒸汽区时的饱和蒸汽熵ssq ;l为迭代Ⅰ次数,没有迭代
0 l
时l=0,psq 为步骤11)测算得到的pc,l大于0时,psq 为上一次迭代Ⅰ中得到的工质进入湿蒸汽区时的饱和压力,;
[0062] 13)根据所述步骤12)中求得的hsql和ssql,通过IFC-67工业用水和水蒸汽热力性l+1质模型求得工质进入湿蒸汽区时的饱和压力psq ;
[0063] 判断是否满足|psql+1-psql|≤r,如是则进入步骤14),否则令l=l+1并回到步骤-312);所述r为工程应用的精度要求,0≤r≤10 ;
[0064] 14)根据步骤13)求得的满足精度要求的psql+1,通过IFC-67工业用水和水蒸汽热力性质模型,计算得到工质进入湿蒸汽区时的饱和蒸汽焓值hsq。
[0065] 经过多次迭代计算,误差r小于10-5,认为已经收敛,得到四种工况下工质进入湿蒸汽区时的饱和蒸汽压力psq分别为0.0586MPa、0.0409MPa、0.0298MPa和0.0227MPa;四种工况下的工质进入湿蒸汽区时的饱和蒸汽焓值hsq分别为2652.597kJ/kg、2637.732kJ/kg、2625.194kJ/kg和2614.706kJ/kg;
[0066] 15)将步骤14)求得的hsq代入公式 求解湿蒸汽级的理想焓降在低压缸理想焓降中所占的百分比a;计算得到四种工况下的a值分别为0.430、0.381、0.332和
0.286。
[0067] 16)根据公式 求解b值;得到四种工况下的b值分别为0.0944、0.0791、0.0792和0.0853。
[0068] 2)根据步骤1)中得到的不同功率运行工况下的Δht值和b值,利用最小二乘方法-6对b值进行二次拟合,得到第一常数t1为6.974×10 ,第二常数t2为-0.0105,第三常数t3-6 2
为4.004,从而得到b值与Δht值的函数关系为b=6.974×10 Δht-0.0105Δht+4.004;
[0069] 3)求取任意功率运行工况,如阀门全开工况(VWO工况)下的排汽焓hcu,具体步骤为:
[0070] 31)按照步骤1)中所述方法求得VWO工况下的p0、t0、h0、s0、pc、hct、Δht、hg和hl分别为0.646MPa、296.109℃、3053.000kJ/kg、7.324kJ/(kg·K)、0.005MPa、2242.755kJ/kg、810.245kJ/kg、2564.015kJ/kg和143.388kJ/kg。
[0071] 32)将所述步骤31)中求得的Δht代入b=6.974×10-6Δht2-0.0105Δht+4.004中,求得当前工况下的b值为0.10395;
[0072] 33)将所述步骤1)中的任一功率运行工况,如额定工况THA工况作为基准工况,k然后根据hce 和所述步骤31)中求得的pc,通过IFC-67工业用水和水蒸汽热力性质模型,k
计算得到本次迭代Ⅱ中的未考虑排汽损失时的排汽熵sce,然后按照步骤12)至13)的方法,计算得到本次迭代Ⅱ中的当前工况下的湿蒸汽级的理想焓降在低压缸理想焓降中所占k 0
的百分比a ;k为迭代Ⅱ次数,没有迭代时k=0,hce 为步骤1)在基准工况下的求得的hce,kk
大于0时,hce 为上一次迭代Ⅱ中得到的当前工况下的未考虑排汽损失时的排汽焓,;
[0073] 34)将所述步骤31)中求得的hg、hl、h0、hct,步骤32)中求得的b值,以及步骤33)k中求得的a 值代入公式 求解得到当前
k+1
工况下的未考虑排汽损失时的排汽焓hce ;
[0074] 判断是否满足|hcek+1-hcek|≤e,如是则进入步骤35),否则令k=k+1并回到步骤-333);所述e为工程应用的精度要求,0≤e≤10 ;
[0075] 经过多次迭代计算,误差e小于10-5,认为已经收敛,得到ak值为0.442,当前VWOk+1工况下的未考虑排汽损失时的排汽焓hce 为2326.810kJ/kg。
[0076] 35)将步骤11)中得到的THA工况下的排汽流量Gc、THA工况下的低压缸进口工质压力p0分别用Gc0和p00来表示,通过公式 求解当前VWO工况下的排汽流量Gd,式中,p0为所述步骤31)中求得的当前VWO工况下的低压缸进口工质压力;计算得到VWO工况下的排汽流量Gc为663.735t/h。
[0077] 36)根据IFC-67工业用水和水蒸汽热力性质模型,通过步骤34)中得到的hcek+1和3
步骤31)中测得的pc计算得到当前VWO工况下的排汽比容v为23.871m/kg;利用汽轮机的排汽损失-排汽容积流量图,根据步骤35)中得到的Gd和当前VWO工况下的排汽比容v查得当前VWO工况下的排汽损失Δhex为16.264kJ/kg,将当前VWO工况下的排汽损失Δhexk+1
加上步骤34)所得未考虑排汽损失的排汽焓值hce ,即可求得当前VWO工况下汽轮机的实际排汽焓值hcu为2343.074kJ/kg。
[0078] 实施例2:根据五种不同工况检测上汽某300MW汽轮机任意工况下的排汽焓值。
[0079] 1)分别在五种不同功率运行工况:额定工况(THA工况)、最大连续出力(TMCR工况)、75%额定工况(75%THA工况)、60%额定工况(60%THA工况)和50%额定工况(50%THA工况)下测算得到汽轮机的低压缸的理想焓降Δht及低压缸喷嘴损失、动叶损失、叶高损失、扇形损失、叶轮摩擦损失、漏汽损失之和占低压缸理想焓降的百分比b,得到各自工况下的Δht和b,所述测算Δht值和b值的方法如下:
[0080] 11)利用温度压力测量仪表测量得到汽轮机低压缸进口工质压力p0和低压缸进口工质温度t0,根据p0和t0,通过国际公式化委员会(IFC)提出的IFC-67工业用水和水蒸汽热力性质模型计算低压缸进口工质焓h0和低压缸进口工质熵s0;
[0081] 利用压力测量仪表测量得到汽轮机低压缸排汽压力pc,
[0082] 根据IFC-67工业用水和水蒸汽热力性质模型,通过s0和pc计算工质自低压缸进口等熵膨胀的理想排汽焓hct;将h0减去hct,得到低压缸的理想焓降Δht;
[0083] 根据IFC-67工业用水和水蒸汽热力性质模型,通过pc计算排汽压力下的饱和蒸汽焓hg和排汽压力下的饱和水焓hl;
[0084] 根据设计数据查得汽轮机排汽流量Gc和排汽焓hcu;
[0085] 根据IFC-67工业用水和水蒸汽热力性质模型,通过hcu和pc计算排汽比容v;
[0086] 利用汽轮机的排汽损失-排汽容积流量图,根据Gc和v查得排汽损失Δhex;将设计排汽焓hcu减去排汽损失Δhex,得到未考虑排汽损失时的排汽焓hce;
[0087] 根据IFC-67工业用水和水蒸汽热力性质模型,通过hce和pc计算得到未考虑排汽损失时的排汽熵sce;
[0088] 上述五种不同工况下的各热力参数如表2所示。
[0089] 表2五种不同工况下测算得到的各热力参数
[0090]
[0091] 12)根据psql通过IFC-67工业用水和水蒸汽热力性质模型,计算得到本次迭代Ⅰ中的工质进入湿蒸汽区时的饱和蒸汽焓值hsql,然后代入公式 计l
算得到本次迭代Ⅰ中的工质进入湿蒸汽区时的饱和蒸汽熵ssq ;l为迭代Ⅰ次数,没有迭代
0 l
时l=0,psq 为步骤11)测算得到的pc,l大于0时,psq 为上一次迭代Ⅰ中得到的工质进入湿蒸汽区时的饱和压力,;
l l
[0092] 13)根据所述步骤12)中求得的hsq 和ssq,通过IFC-67工业用水和水蒸汽热力性l+1质模型求得工质进入湿蒸汽区时的饱和压力psq ;
l+1 l
[0093] 判断是否满足|psq -psq|≤r,如是则进入步骤14),否则令l=l+1并回到步骤-312);所述r为工程应用的精度要求,0≤r≤10 ;
l+1
[0094] 14)根据步骤13)求得的满足精度要求的psq ,通过IFC-67工业用水和水蒸汽热力性质模型,计算得到工质进入湿蒸汽区时的饱和蒸汽焓值hsq。-5
[0095] 经过多次迭代计算,误差r小于10 ,认为已经收敛,得到五种工况下工质进入湿蒸汽区时的饱和蒸汽压力psq分别为0.0586MPa、0.0619MPa、0.0409MPa、0.0298MPa和0.0227MPa;五种工况下的工质进入湿蒸汽区时的饱和蒸汽焓值hsq分别为2652.597kJ/kg、
2654.895kJ/kg、2637.732kJ/kg、2625.194kJ/kg和2614.706kJ/kg;
[0096] 15)将步骤14)求得的hsq代入公式 求解湿蒸汽级的理想焓降在低压缸理想焓降中所占的百分比a;计算得到五种工况下的a值分别为0.430、0.437、0.381、0.332和0.286。
[0097] 16)根据公式 求解b值;得到五种工况下的b值分别为0.0944、0.0994、0.0791、0.0792和0.0853。
[0098] 2)根据步骤1)中得到的不同功率运行工况下的Δht值和b值,利用最小二乘方法-6对b值进行二次拟合,得到第一常数t1为7.067×10 ,第二常数t2为-0.0106,第三常数t3-6 2
为4.056,从而得到b值与Δht值的函数关系为b=7.067×10 Δht-0.0106Δht+4.056;
[0099] 3)求取任意功率运行工况,如阀门全开工况(VWO工况)下的排汽焓hcu,具体步骤为:
[0100] 31)按照步骤1)中所述方法求得VWO工况下的p0、t0、h0、s0、pc、hct、Δht、hg和hl分别为0.646MPa、296.109℃、3053.000kJ/kg、7.324kJ/(kg·K)、0.005MPa、2242.755kJ/kg、810.245kJ/kg、2564.015kJ/kg和143.388kJ/kg。
[0101] 32)将所述步骤31)中求得的Δht代入b=7.067×10-6Δht2-0.0106Δht+4.056中,求得当前工况下的b值为0.10365;
[0102] 剩余步骤与实例1中步骤33)~步骤36)相同,求解得到当前VWO工况下的未考虑k+1排汽损失时的排汽焓hce 为2327.277kJ/kg;VWO工况下的排汽流量Gc为663.735t/h;VWO
3
工况下的排汽比容v为23.892m/kg;查得当前VWO工况下的排汽损失Δhex为16.028kJ/kg;最终求得当前VWO工况下汽轮机的实际排汽焓值hcu为2343.305kJ/kg。
[0103] 实施例3:根据六种不同工况检测上汽某300MW汽轮机任意工况下的排汽焓值。
[0104] 1)分别在六种不同功率运行工况:额定工况(THA工况)、最大连续出力(TMCR工况)、夏季工况(TRL工况)、75%额定工况(75%THA工况)、60%额定工况(60%THA工况)和50%额定工况(50%THA工况)下测算得到汽轮机的低压缸的理想焓降Δht及低压缸喷嘴损失、动叶损失、叶高损失、扇形损失、叶轮摩擦损失、漏汽损失之和占低压缸理想焓降的百分比b,得到各自工况下的Δht和b,所述测算Δht值和b值的方法如下:
[0105] 11)利用温度压力测量仪表测量得到汽轮机低压缸进口工质压力p0和低压缸进口工质温度t0,根据p0和t0,通过国际公式化委员会(IFC)提出的IFC-67工业用水和水蒸汽热力性质模型计算低压缸进口工质焓h0和低压缸进口工质熵s0;
[0106] 利用压力测量仪表测量得到汽轮机低压缸排汽压力pc,
[0107] 根据IFC-67工业用水和水蒸汽热力性质模型,通过s0和pc计算工质自低压缸进口等熵膨胀的理想排汽焓hct;将h0减去hct,得到低压缸的理想焓降Δht;
[0108] 根据IFC-67工业用水和水蒸汽热力性质模型,通过pc计算排汽压力下的饱和蒸汽焓hg和排汽压力下的饱和水焓hl;
[0109] 根据设计数据查得汽轮机排汽流量Gc和排汽焓hcu;
[0110] 根据IFC-67工业用水和水蒸汽热力性质模型,通过hcu和pc计算排汽比容v;
[0111] 利用汽轮机的排汽损失-排汽容积流量图,根据Gc和v查得排汽损失Δhex;将设计排汽焓hcu减去排汽损失Δhex,得到未考虑排汽损失时的排汽焓hce;
[0112] 根据IFC-67工业用水和水蒸汽热力性质模型,通过hce和pc计算得到未考虑排汽损失时的排汽熵sce;
[0113] 上述六种不同工况下的各热力参数如表3所示。
[0114] 表3六种不同工况下测算得到的各热力参数
[0115]
[0116] 12)根据psql通过IFC-67工业用水和水蒸汽热力性质模型,计算得到本次迭代Ⅰl中的工质进入湿蒸汽区时的饱和蒸汽焓值hsq,然后代入公式 计
l
算得到本次迭代Ⅰ中的工质进入湿蒸汽区时的饱和蒸汽熵ssq ;l为迭代Ⅰ次数,没有迭代
0 l
时l=0,psq 为步骤11)测算得到的pc,l大于0时,psq 为上一次迭代Ⅰ中得到的工质进入湿蒸汽区时的饱和压力,;
[0117] 13)根据所述步骤12)中求得的hsql和ssql,通过IFC-67工业用水和水蒸汽热力性l+1质模型求得工质进入湿蒸汽区时的饱和压力psq ;
[0118] 判断是否满足|psql+1-psql|≤r,如是则进入步骤14),否则令l=l+1并回到步骤-312);所述r为工程应用的精度要求,0≤r≤10 ;
[0119] 14)根据步骤13)求得的满足精度要求的psql+1,通过IFC-67工业用水和水蒸汽热力性质模型,计算得到工质进入湿蒸汽区时的饱和蒸汽焓值hsq。
[0120] 经过多次迭代计算,误差r小于10-5,认为已经收敛,得到六种工况下工质进入湿蒸汽区时的饱和蒸汽压力psq分别为0.0586MPa、0.0619MPa、0.0683MPa、0.0409MPa、0.0298MPa和0.0227MPa;六种工况下的工质进入湿蒸汽区时的饱和蒸汽焓值hsq分别为 2652.597kJ/kg、2654.895kJ/kg、2659.039kJ/kg、2637.732kJ/kg、2625.194kJ/kg 和
2614.706kJ/kg;
[0121] 15)将步骤14)求得的hsq代入公式 求解湿蒸汽级的理想焓降在低压缸理想焓降中所占的百分比a;计算得到六种工况下的a值分别为0.430、0.437、0.376、
0.381、0.332和0.286。
[0122] 16)根据公式 求解b值;得到六种工况下的b值分别为0.0944、0.0994、0.0889、0.0791、0.0792和0.0853。
[0123] 2)根据步骤1)中得到的不同功率运行工况下的Δht值和b值,利用最小二乘方法-6对b值进行二次拟合,得到第一常数t1为7.105×10 ,第二常数t2为-0.0107,第三常数t3-6 2
为4.077,从而得到b值与Δht值的函数关系为b=7.105×10 Δht-0.0107Δht+4.077;
[0124] 3)求取任意功率运行工况,如阀门全开工况(VWO工况)下的排汽焓hcu,具体步骤为:
[0125] 31)按照步骤1)中所述方法求得VWO工况下的p0、t0、h0、s0、pc、hct、Δht、hg和hl分别为0.646MPa、296.109℃、3053.000kJ/kg、7.324kJ/(kg·K)、0.005MPa、2242.755kJ/kg、810.245kJ/kg、2564.015kJ/kg和143.388kJ/kg。-6 2
[0126] 32)将所述步骤31)中求得的Δht代入b=7.105×10 Δht-0.0107Δht+4.077中,求得当前工况下的b值为0.10347;
[0127] 剩余步骤与实例1中步骤33)~步骤36)相同,求解得到当前VWO工况下的未考虑k+1排汽损失时的排汽焓hce 为2327.631kJ/kg;VWO工况下的排汽流量Gc为663.735t/h;VWO
3
工况下的排汽比容v为23.659m/kg;查得当前VWO工况下的排汽损失Δhex为15.841kJ/kg;最终求得当前VWO工况下汽轮机的实际排汽焓值hcu为2343.427kJ/kg。
[0128] 通过上述三个实例可以看出,本发明公布的测算方法能够根据不同个数的工况检测汽轮机任意工况下的排汽焓值,工况的数目直接影响拟合得到的b值与Δht值的函数关系,从而影响最终的检测结果。
[0129] 利用本发明公布的测算方法对实例进行排汽焓计算的结果表明,本发明测算过程简单清晰,所需测点少,在不同负荷下计算的排汽焓值都有较高的精度。
[0130] 以上所述仅为本发明的较佳实施方式,本发明的保护范围并不以上述实施方式为限,但凡本领域普通技术人员根据本发明所揭示内容所作的等效修饰或变化,皆应纳入权利要求书中记载的保护范围内。