CN114856726B - 一种饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法 - Google Patents
一种饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法 Download PDFInfo
- Publication number
- CN114856726B CN114856726B CN202210554966.4A CN202210554966A CN114856726B CN 114856726 B CN114856726 B CN 114856726B CN 202210554966 A CN202210554966 A CN 202210554966A CN 114856726 B CN114856726 B CN 114856726B
- Authority
- CN
- China
- Prior art keywords
- steam
- saturated
- humidity
- specific enthalpy
- mass flow
- 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
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
- F01D—NON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
- F01D21/00—Shutting-down of machines or engines, e.g. in emergency; Regulating, controlling, or safety means not otherwise provided for
- F01D21/003—Arrangements for testing or measuring
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
- F01D—NON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
- F01D21/00—Shutting-down of machines or engines, e.g. in emergency; Regulating, controlling, or safety means not otherwise provided for
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
- F01D—NON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
- F01D21/00—Shutting-down of machines or engines, e.g. in emergency; Regulating, controlling, or safety means not otherwise provided for
- F01D21/14—Shutting-down of machines or engines, e.g. in emergency; Regulating, controlling, or safety means not otherwise provided for responsive to other specific conditions
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
- F01D—NON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
- F01D25/00—Component parts, details, or accessories, not provided for in, or of interest apart from, other groups
- F01D25/30—Exhaust heads, chambers, or the like
-
- 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
- Y02E30/00—Energy generation of nuclear origin
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Control Of Turbines (AREA)
Abstract
本发明提供了一种饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法,包括:S1,获取汽轮机低压缸和低压给水加热器的设备结构参数及该汽轮机低压缸和低压给水加热器在给定时刻下的运行数据;S2,计算各级组的干基内效率和汽水分离效率;S3,计算饱和蒸汽在级组内膨胀后的比焓;S4,计算汽水分离后的饱和蒸汽质量流量和比焓;S5,计算抽汽后饱和蒸汽的质量流量;S6,将前一级组出口蒸汽的比焓和质量流量作为下一级组进口的比焓和质量流量,依次计算各级组出口的饱和蒸汽比焓及质量流量,得到乏汽湿度;S7,计算饱和蒸汽汽轮机乏汽容积流量。本发明可用于在线计算饱和蒸汽汽轮机乏汽的容积流量及湿度,解决了该变量无法实时测量的难题。
Description
技术领域
本发明涉及核电站运行优化控制技术领域,具体的,涉及一种饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法。
背景技术
饱和蒸汽汽轮机是一种以饱和蒸汽作为工质的汽轮机,饱和蒸汽是干饱和蒸汽和饱和水的混合物。在核电机组中,除了低压缸的前几级工作在过热蒸汽区域,高压缸的全部级和低压缸的末几级均工作在饱和蒸汽区域。
饱和蒸汽汽轮机乏汽容积流量及湿度可用于汽轮机运行状态的监测,目的是提高核电机组运行的安全性和经济性。若乏汽容积流量过大,则会导致低压缸末级叶片的震颤,影响机组运行的安全性;还将导致汽轮机余速损失的增加,降低机组的经济性。乏汽湿度过高则会引发叶片冲蚀,降低设备的使用寿命;同时也会加剧汽轮机的湿汽损失,从而降低汽轮机的内效率。
因饱和蒸汽汽轮机乏汽的容积流量巨大(1000MW机组的乏汽容积流量约为3000m3/s),且处于气液两相流状态,故难以通过流量计直接测量乏汽的质量流量或容积流量。同时,乏汽湿度的在线测量也有很大难度:同位素追踪方法受半衰期影响,无法长期使用;光学方法受液滴粒径影响大,精度较低;热力学方法无法实现在线连续测量。至于软测量方法,因饱和蒸汽汽轮机的各级抽汽均为湿饱和蒸汽,各级抽汽的质量流量和湿度(或比焓)均为未知变量,故现有的面向火电机组过热蒸汽汽轮机的软测量方法均不适用于饱和蒸汽汽轮机。
经过对现有技术检索后发现,申请号为CN202110229455.0,申请日为2021年03月02日的中国发明专利《汽轮机末级叶片容积流量估计方法、颤振预警方法及***和装置》与本发明相关,该专利提供的方法基于汽轮机各状态变量的标幺值,在变工况下对汽轮机乏汽容积流量进行修正,可用于火电厂汽轮机末级叶片容积流量的估算。但该方法包含众多线性化假设,如主蒸汽压力、抽汽压力、蒸汽膨胀比等,而这些线性化假设未经验证,且真实情况是非线性的,据此假设进行的计算准确性存疑。更为重要的是,该专利提供的方法仅适用于过热蒸汽,而饱和蒸汽汽轮机各级抽汽均为饱和蒸汽。
事实上,现有公开报道的文献和专利均未涉及饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法,这一空白有待填补。
发明内容
针对现有技术的不足,本发明提供了一种饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法。
本发明的第一方面,提供一种饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法,包括:
S1:获取汽轮机低压缸和低压给水加热器的设备结构参数,汽轮机的设计参数,以及该汽轮机低压缸和低压给水加热器在给定时刻下的运行数据;
S2:根据S1的汽轮机设计参数计算各级组的干基内效率和汽水分离效率;
S3:根据S1得到的低压缸进口工质比焓和压力以及抽汽口工质压力和S2的干基内效率,计算饱和蒸汽在级组内膨胀后的比焓;
S4:根据S2得到的级组的汽水分离效率和S3得到的饱和蒸汽在级组内膨胀后的比焓,计算汽水分离后的饱和蒸汽质量流量、比焓和湿度,将汽水分离后的饱和蒸汽比焓作为级组出口饱和蒸汽比焓;
S5:根据S1运行数据中的低压给水加热器进出口给水质量流量和比焓和S4计算得到的级组出口饱和蒸汽比焓,计算抽汽后饱和蒸汽的质量流量,将抽汽的饱和蒸汽质量流量作为级组出口饱和蒸汽的质量流量;
S6:将前一级组出口蒸汽的比焓和质量流量作为下一级组进口的比焓和质量流量,重复S3-S5,依次计算各级组出口的饱和蒸汽比焓及质量流量,直至计算完所有的级组;得到饱和蒸汽汽轮机乏汽的比焓、湿度和质量流量。
S7:根据S6中乏汽比焓和质量流量计算得到饱和蒸汽汽轮机乏汽容积流量。
本发明第二方面,提供一种可读介质,所述可读介质上存储有指令,该指令在电子设备上执行时使电子设备执行所述的饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法。
本发明第三方面,提供一种电子设备,包括存储器、处理器及存储在存储器上并能在处理器上运行的计算机程序,所述处理器执行所述程序时用于执行所述的饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法。
与现有技术相比,本发明实施例具有以下至少一种有益效果:
本发明提供的饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法,可用于在线确定饱和蒸汽汽轮机乏汽的容积流量及湿度,解决了乏汽容积流量及湿度无法实时测量的难题;得到的容积流量及湿度可用于汽轮机低压缸运行状态的监测,有提升核电厂饱和蒸汽汽轮机安全性和经济性的应用潜力。
本发明提供的饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法,适用于核电站饱和蒸汽汽轮机和其它应用领域的饱和蒸汽汽轮机。它与目前应用于火电汽轮机的相关技术的根本区别在于,火电汽轮机的抽汽处于过热状态,比焓可根据温度、压力计算得到;而饱和蒸汽汽轮机抽汽的比焓(湿度)和质量流量均为未知变量。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显。
图1为本发明一实施例中低压缸、低压给水加热***、凝汽器结构示意图;
图2为本发明一实施例中一种饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法的流程图;
图3为本发明一实施例中某核电机组输出功率变化图;
图4为本发明一实施例中饱和蒸汽汽轮机乏汽容积流量计算结果;
图5为本发明一实施例中饱和蒸汽汽轮机乏汽湿度计算结果;
图6为本发明一实施例中乏汽放热量与循环水吸热量对比。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进。这些都属于本发明的保护范围。
本发明针对汽轮机划分为若干串联的级组,根据级组进口饱和蒸汽的质量流量、比焓和湿度,计算得到级组出口饱和蒸汽的质量流量、比焓和湿度;从饱和蒸汽汽轮机进口依次计算各个级组,得到汽轮机乏汽的质量流量和比焓,最终得到饱和蒸汽汽轮机乏汽容积流量及湿度。
图1为本发明一实施例中低压缸、低压给水加热***、凝汽器结构示意图,该实施例应用于某1000MW压水堆核电机组;再热蒸汽进入低压缸做功后成为4级抽汽和乏汽,汽轮机每两级抽汽中间的若干级被视作一个级组,前一级组出口的蒸汽除进入抽汽管道外,剩余部分作为下一级组的入口蒸汽。
参照图2所示,为本发明一实施例中饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法的流程图。具体的,本实施例中饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法包括如下步骤:
S1:获取汽轮机低压缸和低压给水加热器的设备结构参数,汽轮机的设计参数,以及该汽轮机低压缸和低压给水加热器在给定时刻下的运行数据;
S2:根据S1的汽轮机设计参数计算各级组的干基内效率和汽水分离效率;
S3:根据S1的低压缸进口工质比焓和压力以及抽汽口工质压力和S2的干基内效率,计算饱和蒸汽在级组内膨胀后的比焓;
S4:根据S2得到的级组的汽水分离效率和S3得到的饱和蒸汽在级组内膨胀后的比焓,计算汽水分离后的饱和蒸汽质量流量、比焓和湿度,将汽水分离后的饱和蒸汽比焓作为级组出口饱和蒸汽比焓;
S5:根据S1运行数据中的低压给水加热器进出口给水质量流量和比焓和S4得到的级组出口饱和蒸汽比焓,计算抽汽后饱和蒸汽的质量流量,将抽汽的饱和蒸汽质量流量作为级组出口饱和蒸汽的质量流量;
S6:将前一级组出口蒸汽的比焓和质量流量作为下一级组进口的比焓和质量流量,重复S3-S5,依次计算各级组出口的饱和蒸汽比焓及质量流量,直至计算完所有的级组;将最后一个级组出口饱和蒸汽的比焓hout、湿度yout和质量流量Dout分别作为饱和蒸汽汽轮机乏汽的比焓hexh、湿度yexh和质量流量Dexh;
S7:根据S6中乏汽比焓和质量流量计算得到饱和蒸汽汽轮机乏汽容积流量。
作为一优选实施例,S1中,所述运行数据可通过现场的DCS获取,包括:汽轮机进口工质质量流量、比焓和湿度,各级组进、出口的抽汽压力、低压给水加热器进出口给水质量流量和比焓。
设备结构参数包括:汽轮机的级数和抽汽口位置、各给水加热器的外形尺寸和换热管的长度、管径,可根据汽轮机和给水加热器的规程得到;
汽轮机的设计参数包括:汽轮机各级组进、出口工质的温度、压力、比焓和湿度的设计值,上述设计值可根据汽轮机规程得到。
上述的设备结构参数和汽轮机的设计参数,根据产品手册查阅或者厂家提供。
作为一优选实施例,S2中,根据汽轮机设计参数计算各级组的干基内效率和汽水分离效率:
其中:
式中,ηdry是级组的干基内效率,无量纲;ηref是级组在设计工况的内效率,即汽轮机的最大内效率,无量纲;α是级组的鲍曼因子,无量纲;yin ref和yout ref分别是级组在设计工况下的进、出口蒸汽的湿度,无量纲;hin ref和hout ref分别是级组在设计工况下的进、出口蒸汽的比焓,kJ/kg;pin ref和pout ref分别是级组在设计工况下的进、出口蒸汽的压力,MPa;ρin ref和ρout ref分别是级组在设计工况下的进、出口蒸汽的密度,kg/m3;κ是级组的对数平均绝热指数,无量纲。
本优选实施例中,干基内效率采用上述公式计算,较现有方法具有实施方便、计算结果准确的优点。
作为一优选实施例,所述级组的汽水分离效率采用以下公式计算得到:
式中,ηsp是级组的汽水分离效率,无量纲;Ddw ref是设计工况下汽轮机级间疏水的质量流量,kg/s;Din ref是设计工况下级组进口工质的质量流量,kg/s。
本优选实施例中,采用上述的汽水分离效率公式,在级组内部汽水分离装置的汽水分离过程,解决了这一过程无测点的难题。
作为一优选实施例,S3中,所述饱和蒸汽在级组内膨胀后的比焓,计算包括以下步骤:
(a)计算级组进口工质的湿度yin和比熵sin:
sin=(1-yin)sV(pin)+yinsL(pin)
式中,yin是级组进口工质的湿度,无量纲;hin是级组进口工质的比焓,kJ/kg;pin是级组进口工质压力,MPa;sin是级组进口工质的比熵,kJ/(kg*K);hL(pin)和hV(pin)分别是级组进口工质压力对应的饱和水、饱和干蒸汽的比焓,kJ/kg;sL(pin)和sV(pin)分别是级组进口工质压力对应的饱和水、饱和干蒸汽的比熵,kJ/(kg*K)。
(b)根据(a)中计算得到的级组进口工质比熵sin,计算在等熵膨胀时级组出口蒸汽的湿度yout,s和比焓hout,s:
hexp,s=(1-yexp,s)hV(pout)+yexp,shL(pout)
式中,yexp,s是等熵膨胀时级组出口工质的湿度,无量纲;hexp,s是等熵膨胀时级组出口工质的比焓,kJ/kg;pout是级组出口工质压力,MPa;sL(pout)和sV(pout)分别是级组出口工质压力对应的饱和水、饱和干蒸汽的比熵,kJ/(kg*K);hL(pout)和hV(pout)分别是级组出口工质压力对应的饱和水、饱和干蒸汽的比焓,kJ/kg。
(c)计算级组的内效率:
ηexp=ηdry(1-αyin)
式中,ηexp是变工况下的级组内效率,无量纲;ηdry是级组的干基内效率,无量纲;yin是级组进口蒸汽湿度,无量纲。
(d)根据(b)得到的等熵膨胀时级组出口蒸汽的比焓hexp,s和(c)中的级组内效率ηexp,计算蒸汽在级组内膨胀后的比焓及湿度:
hexp=hin-ηexp(hin-hexp,s)
式中,hexp是蒸汽在级组内膨胀后的比焓,kJ/kg;yexp是蒸汽在级组内膨胀后的湿度,无量纲;hL(pout)和hV(pout)分别是级组出口工质压力对应的饱和水、饱和干蒸汽的比焓,kJ/kg。
(e)根据(d)中得到的级组出口湿度yexp修正级组内效率ηexp:
(f)重复(d)得到蒸汽在级组内膨胀后的比焓hexp及湿度yexp。
上述步骤,(c)中ηexp和(d)中hexp会相互影响,因此需要用(d)的结果重新计算ηexp,即(e)步骤,再用重新计算得到的ηexp来再算hexp。这样得到的结果更为准确,符合实际的汽轮机的情况。
本优选实施例中,采用上述步骤(a)-(f)计算饱和蒸汽在级组内膨胀后的比焓的过程中,考虑了蒸汽湿度对级组内效率的影响,与现有技术相比,结果的准确性更高。
作为一优选实施例,S4中,所述汽水分离后的饱和蒸汽比焓采用以下公式计算得到:
式中,hsep是汽水分离后的饱和蒸汽比焓,kJ/kg;Din是级组进口工质质量流量,kg/s;hexp是蒸汽在级组内膨胀后的比焓,kJ/kg;ηsp是级组的汽水分离效率,无量纲;yexp是蒸汽在级组内膨胀后的湿度,无量纲;pout是级组出口工质压力,MPa;hL(pout)是级组出口工质压力对应的饱和水的比焓,kJ/kg。
本优选实施例中,采用上述的汽水分离后的饱和蒸汽比焓计算公式,考虑了汽轮机级组内部除湿装置对蒸汽比焓的影响,弥补了现有技术的空缺。
作为一优选实施例,所述汽水分离后的饱和蒸汽湿度采用以下公式计算得到:
式中,ysep是汽水分离后的饱和蒸汽湿度,无量纲;hV(pout)分别是级组出口工质压力对应的饱和干蒸汽的比焓,kJ/kg。
本优选实施例中,采用上述的汽水分离后的饱和蒸汽湿度计算公式,考虑了汽轮机级组内部除湿装置对饱和蒸汽湿度的影响,弥补了现有技术的空缺。
作为一优选实施例,所述汽水分离后的饱和蒸汽质量流量采用以下公式计算得到:
Dsep=(1-ηspyexp)Din
式中,Dsep是汽水分离后的饱和蒸汽质量流量,kg/s;Din是级组进口饱和蒸汽质量流量,kg/s。
本优选实施例中,采用上述的汽水分离后的饱和蒸汽质量流量计算公式,解决了汽轮机的质量流量无法测量的难题。
作为一优选实施例,将汽水分离后的饱和蒸汽比焓hsep作为级组出口饱和蒸汽的比焓hout;将汽水分离后的饱和蒸汽湿度ysep作为级组出口饱和蒸汽的湿度yout。
作为一优选实施例,S5中,抽汽后饱和蒸汽的质量流量采用以下公式计算得到:
Dext=Dsep-Des
式中,Dext是抽汽后的饱和蒸汽的质量流量,kg/s;Dsep是汽水分离后的饱和蒸汽质量流量,kg/s;Des是抽汽质量流量,kg/s。
本优选实施例中,采用上述的抽汽后饱和蒸汽的质量流量计算公式,解决了汽轮机抽汽后的饱和蒸汽质量流量无法测量的难题。
作为一优选实施例,所述抽汽质量流量Des通过求解以下方程组得到:
式中,VV是给水加热器内蒸汽体积,m3;ρV是饱和干蒸汽密度,kg/m3;hL和hV分别是饱和水、饱和干蒸汽比焓,kJ/kg;Des是抽汽质量流量,kg/s;Ddw,in是上级给水加热器疏水的质量流量,kg/s;hdw,in是上级给水加热器疏水的比焓,kJ/kg;Dcon是凝结水质量流量,kg/s;QV是蒸汽凝结段和疏水冷却段的热流量,kW;SV是给水加热器蒸汽凝结段的换热面积,m2;Tes和Tfw分别是抽汽温度和给水温度,K;Kfw和Kes分别是给水和蒸汽的表面传热系数,kW/(m2*K);wfw是给水的流速,m/s;ρfw是给水的密度,kg/m3,λfw和λL,es分别是给水的热导率和抽汽压力对应的饱和水热导率,kW/(m*K);cp,fw是给水的定压比热容,kW/(kg*K);μfw和μL,es分别是给水的动力粘度和抽汽压力对应的饱和水动力粘度,Pa·s;di是换热管内径,m;g是重力加速度,取9.8m/s2;ρL,es是抽汽压力对应的饱和水密度,kg/m3;do是给水加热器的换热管外径,m。
将抽汽后的饱和蒸汽质量流量Dext作为级组出口饱和蒸汽的质量流量Dout。
作为一优选实施例,S7中,所述饱和蒸汽汽轮机乏汽容积流量采用以下公式计算得到:
式中,Vexh是饱和蒸汽汽轮机乏汽容积流量,m3/s;Dexh是S6得到的饱和蒸汽汽轮机乏汽质量流量,kg/s;yexh是S6得到的饱和蒸汽汽轮机乏汽湿度,无量纲;pexh是饱和蒸汽汽轮机乏汽压力,MPa;ρL(pexh)和ρV(pexh)分别是蒸汽汽轮机乏汽压力对应和饱和水密度和饱和干蒸汽密度,kg/m3。
作为一优选实施例,所述饱和蒸汽汽轮机乏汽湿度即为S6中的yexh。
在本发明的基于饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法的实施例中,首先要获取核电机组的DCS实测数据。其中,机组输出功率见图3。
图4为本发明一实施例中饱和蒸汽汽轮机乏汽容积流量的计算结果,本发明基于饱和蒸汽汽轮机的现有测点数据确定乏汽容积流量,可解决乏汽容积流量无测点的难题。
图5为本发明一实施例中饱和蒸汽汽轮机乏汽湿度的计算结果,本发明可解决饱和蒸汽汽轮机乏汽湿度无法在线测量的难题。
图6为本发明一实施例中乏汽放热量与循环水吸热量对比,其中,循环水吸热量根据循环水泵电流和循环水进出口温度计算得到。两者的平均相对误差为1.70%,其一致性可间接验证乏汽湿度和容积流量在线确定的正确性。
基于相同技术构思,在本发明另一实施例中,还提供一种可读介质,可读介质上存储有指令,该指令在电子设备上执行时使电子设备执行上述任一项实施例的饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法。
基于相同技术构思,在本发明另一实施例中,还提供一种电子设备,包括存储器、处理器及存储在存储器上并能在处理器上运行的计算机程序,处理器执行程序时用于执行上述任一项实施例的饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法。
为了更好地说明本发明的实施过程,下面以图1所示汽轮机中最具代表性的低压缸前流道第3级组(该级组包含了第2级抽汽口至第1级抽汽口之间的全部级)为例,展示某时刻下饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法:
S1,通过现场的DCS获取运行数据,包括:汽轮机进口工质质量流量1180kg/s,比焓2856kJ/kg;级组进、出口的抽汽压力分别为和0.07426MPa和0.01588MPa;对应的低压给水加热器进口给水质量流量为915.1kg/s,比焓为175.9kJ/kg。
额外地,在本发明的S6中,将前一级组出口蒸汽的比焓和质量流量作为下一级组进口的比焓和质量流量,故有:级组进口工质的比焓为2465kJ/kg;级组进口工质的质量流量为1027kg/s;“前一级组出口蒸汽的比焓和质量流量”可将本发明的S1-S5实施于前一级组得到;
查阅汽轮机运行规程得到设计参数,包括:级组进、出口工质的温度设计值分别为366.9K和329.1K,压力设计值分别为0.08086MPa和0.01750MPa、比焓设计值分别为2548kJ/kg和2390kJ/kg,以及湿度的设计值分别为5.20%和9.05%;
S2,根据汽轮机设计参数计算各级组的干基内效率和汽水分离效率:
其中:
式中,ηref是级组在设计工况的内效率,计算结果为85.05%;κ是级组的对数平均绝热指数,计算结果为1.099;α是级组的鲍曼因子,取值为1.2;yin ref和yout ref分别是级组在设计工况下的进、出口蒸汽的湿度,分别为5.20%和9.05%;hin ref和hout ref分别是级组在设计工况下的进、出口蒸汽的比焓,分别为2548kJ/kg和2355kJ/kg;pin ref和pout ref分别是级组在设计工况下的进、出口蒸汽的压力,分别为0.08086MPa和0.01750MPa;ρin ref和ρout ref分别是级组在设计工况下的进、出口蒸汽的密度,分别为0.5105kg/m3和0.1268kg/m3;
得到级组的干基内效率ηdry为0.9310%。
级组的汽水分离效率采用以下公式计算得到:
式中,Ddw ref是设计工况下汽轮机级间疏水的质量流量,为60.48kg/s;Din ref是设计工况下级组进口工质的质量流量,890.9kg/s。
得到级组的汽水分离效率ηsp为16.11%。
S3,饱和蒸汽在级组内膨胀后的比焓,计算包括以下步骤:
(a)计算级组进口工质的湿度yin和比熵sin:
sin=(1-yin)sV(pin)+yinsL(pin)
式中,hin是级组进口工质的比焓,为2465kJ/kg;pin是级组进口工质压力,为0.01588MPa;hL(pin)和hV(pin)分别是级组进口工质压力对应的饱和水、饱和干蒸汽的比焓,分别为383.2kJ/kg和2662kJ/kg;sL(pin)和sV(pin)分别是级组进口工质压力对应的饱和水、饱和干蒸汽的比熵,分别为1.210kJ/(kg*K)和7.459kJ/(kg*K)。
得到级组进口工质的湿度yin为8.647%;级组进口工质比熵sin为6.919kJ/(kg*K)。
(b)根据(a)中计算得到的级组进口工质比熵sin,计算在等熵膨胀时级组出口蒸汽的湿度yout,s和比焓hout,s:
hexp,s=(1-yexp,s)hV(pout)+yexp,shL(pout)
式中,pout是级组出口工质压力,为0.01588MPa;sL(pout)和sV(pout)分别是级组出口工质压力对应的饱和水、饱和干蒸汽的比熵,分别为0.7700kJ/(kg*K)和7.987kJ/(kg*K);hL(pout)和hV(pout)分别是级组出口工质压力对应的饱和水、饱和干蒸汽的比焓,分别为230.9kJ/kg和2600kJ/kg。
得到等熵膨胀时级组出口工质的湿度yexp,s为14.81%;等熵膨胀时级组出口工质的比焓hexp,s为2250kJ/kg。
(c)计算级组的内效率:
ηexp=ηdry(1-αyin)
式中,dry是级组的干基内效率,在S2中得到,为93.10%;yin是级组进口蒸汽湿度,为8.650%。
得到变工况下的级组内效率ηexp为83.44%。
(d)根据(b)得到的等熵膨胀时级组出口蒸汽的比焓hexp,s和(c)中的级组内效率ηexp,计算蒸汽在级组内膨胀后的比焓及湿度:
hexp=hin-ηexp(hin-hexp,s)
得到蒸汽在级组内膨胀后的比焓hexp为2285kJ/kg;蒸汽在级组内膨胀后的湿度yexp为13.30%。
(e)根据(d)中得到的级组出口湿度yexp重新计算级组内效率ηexp:
得到级组内效率ηexp为80.84%。
(f)利用(e)中的级组内效率,重复(d),得到蒸汽在级组内膨胀后的比焓hexp及湿度yexp。
得到蒸汽在级组内膨胀后的比焓hexp为2291kJ/kg;蒸汽在级组内膨胀后的湿度yexp为13.09%。
S4,汽水分离后的饱和蒸汽比焓采用以下公式计算得到:
式中,Din是级组进口工质质量流量,1027kg/s;ηsp是级组的汽水分离效率,在S2中得到,为16.11%;
得到汽水分离后的饱和蒸汽比焓hsep为2335kJ/kg。
汽水分离后的饱和蒸汽湿度采用以下公式计算得到:
得到汽水分离后的饱和蒸汽湿度ysep为11.21%。
汽水分离后的饱和蒸汽质量流量采用以下公式计算得到:
Dsep=(1-ηspyexp)Din
得到汽水分离后的饱和蒸汽质量流量Dsep为1005kg/s。
将汽水分离后的饱和蒸汽比焓hsep作为级组出口饱和蒸汽的比焓hout;将汽水分离后的饱和蒸汽湿度ysep作为级组出口饱和蒸汽的湿度yout。
S5,抽汽后饱和蒸汽的质量流量采用以下公式计算得到:
Dext=Dsep-Des
得到抽汽后的饱和蒸汽的质量流量Dext为980.8kg/s;
其中,Des是抽汽质量流量,为24.40kg/s,通过求解以下方程组得到:
式中,VV是给水加热器内蒸汽体积,为60.28m3;ρV是饱和干蒸汽密度,为0.2864kg/m3;hL和hV分别是饱和水、饱和干蒸汽比焓,分别为230.7kJ/kg和2601kJ/kg;Ddw,in是上级给水加热器疏水的质量流量,该给水加热器无上级疏水,故为0kg/s;Dcon是凝结水质量流量,46.04kg/s;QV是蒸汽凝结段和疏水冷却段的热流量,计算结果为1.761×105kW;SV是给水加热器蒸汽凝结段的换热面积,为1600m2;Tes是抽汽温度,为328.3K;Tfw是给水温度,为314.8K至317.7K的温度分布;Kfw和Kes分别是给水和蒸汽的表面传热系数,分别为10.73kW/(m2*K)至11.78kW/(m2*K)和10.72kW/(m2*K)至80.47kW/(m2*K)的分布;wfw是给水的流速,为2.243m/s至2.255m/s的分布;ρfw是给水的密度,为991.5kg/m3至986.2kg/m3的分布;λfw是给水的热导率,为0.6346kW/(m*K)至0.6480kW/(m*K)的分布;λL,es分别是抽汽压力对应的饱和水热导率,为0.6480kW/(m*K);cp,fw是给水的定压比热容,为4.169kW/(kg*K)至4.172kW/(kg*K)的分布;μfw和μL,es分别是给水的动力粘度和抽汽压力对应的饱和水动力粘度,为6.158×10-4Pa·s至5.033×10-4Pa·s的分布;μL,es是抽汽压力对应的饱和水动力粘度,为5.032×10-4Pa·s;ρL,es是抽汽压力对应的饱和水密度,987.3kg/m3;g是重力加速度,取9.8m/s2;di是换热管内径,为0.0130m;do是给水加热器的换热管外径,为0.0160m;
将抽汽后的饱和蒸汽质量流量Dext作为级组出口饱和蒸汽的质量流量Dout。
S6,将前一级组出口蒸汽的比焓和质量流量作为下一级组进口的比焓和质量流量,重复S3-S5,依次计算各级组出口的饱和蒸汽比焓及质量流量,直至计算完所有的级组;将最后一个级组出口饱和蒸汽的比焓hout、湿度yout和质量流量Dout分别作为饱和蒸汽汽轮机乏汽的比焓hexh、湿度yexh和质量流量Dexh;得到饱和蒸汽汽轮机乏汽的比焓hexh为2289kJ/kg;湿度yexh为11.57%;质量流量Dexh为980.8kg/s;
S7,饱和蒸汽汽轮机乏汽容积流量采用以下公式计算得到:
式中,pexh是饱和蒸汽汽轮机乏汽压力,0.006210MPa;ρL(pexh)和ρV(pexh)分别是蒸汽汽轮机乏汽压力对应和饱和水密度和饱和干蒸汽密度,分别为993.3kg/m3和0.04352kg/m3。
本实施例中,得到饱和蒸汽汽轮机乏汽容积流量Vexh为2608m3/s;饱和蒸汽汽轮机乏汽湿度即为S6中的yexh,为11.57%。
本发明上述各实施例中的饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法,通过各步骤中的参数设计,提高了参数的准确性,可用于在线计算饱和蒸汽汽轮机乏汽的容积流量及湿度,解决了该变量无法实时测量的难题,对饱和蒸汽汽轮机运行状态的监测具有重要意义。
本发明实施例得到的结果可以直接用于对核电饱和蒸汽汽轮机的实时在线监测,具体的,得到上述饱和蒸汽汽轮机乏汽容积流量和饱和蒸汽汽轮机乏汽湿度之后,可以根据设定的监测阈值进行对比,超过设定的监测阈值则进行预警,实现实时在线监测,对于核电饱和蒸汽汽轮机的运行监测起到重要辅助作用。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变形或修改,这并不影响本发明的实质内容。
Claims (9)
1.一种饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法,其特征在于,包括:
S1:获取汽轮机低压缸和低压给水加热器的设备结构参数,汽轮机的设计参数,以及该汽轮机低压缸和低压给水加热器在给定时刻下的运行数据;所述运行数据包括:低压缸进口工质质量流量、比焓、湿度和压力,各级组进、出口的抽汽压力、低压给水加热器进出口给水质量流量和比焓;
S2:根据S1的汽轮机的设计参数计算各级组的干基内效率和汽水分离效率;
S3:根据S1的低压缸进口工质比焓和压力以及各级组进、出口的抽汽压力和S2的干基内效率,计算饱和蒸汽在级组内膨胀后的比焓;
S4:根据S2得到的级组的汽水分离效率和S3得到的饱和蒸汽在级组内膨胀后的比焓,计算汽水分离后的饱和蒸汽质量流量、比焓和湿度,将汽水分离后的饱和蒸汽比焓作为级组出口饱和蒸汽比焓;
S5:根据S1运行数据中的低压给水加热器进出口给水质量流量和比焓和S4得到的级组出口饱和蒸汽比焓,计算抽汽后饱和蒸汽的质量流量,将抽汽后饱和蒸汽的质量流量作为级组出口饱和蒸汽的质量流量;
S6:将前一级组出口蒸汽的比焓和质量流量作为下一级组进口的比焓和质量流量,重复S3-S5,依次计算各级组出口的饱和蒸汽比焓及质量流量,直至计算完所有的级组;将最后一个级组出口饱和蒸汽的比焓hout、湿度yout和质量流量Dout分别作为饱和蒸汽汽轮机乏汽的比焓hexh、湿度yexh和质量流量Dexh;
S7:根据S6中乏汽的比焓和质量流量计算得到饱和蒸汽汽轮机乏汽容积流量。
2.根据权利要求1所述的饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法,其特征在于,S2中,根据汽轮机的设计参数计算各级组的干基内效率和汽水分离效率,采用以下公式:
其中:
式中,ηdry是级组的干基内效率,无量纲;ηref是级组在设计工况的内效率,无量纲;α是级组的鲍曼因子,无量纲;yin ref和yout ref分别是级组在设计工况下的进、出口工质的湿度,无量纲;hin ref和hout ref分别是级组在设计工况下的进、出口工质的比焓,kJ/kg;pin ref和pout ref分别是级组在设计工况下的进、出口工质的压力,MPa;ρin ref和ρout ref分别是级组在设计工况下的进、出口工质的密度,kg/m3;κ是级组的对数平均绝热指数,无量纲;
所述级组的汽水分离效率采用以下公式计算得到:
式中,ηsp是级组的汽水分离效率,无量纲;Ddw ref是设计工况下汽轮机级间疏水的质量流量,kg/s;Din ref是设计工况下级组进口工质的质量流量,kg/s。
3.根据权利要求1所述的饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法,其特征在于,S3中,所述计算饱和蒸汽在级组内膨胀后的比焓,包括:
(a)计算级组进口工质的湿度yin和比熵sin:
sin=(1-yin)sV(pin)+yinsL(pin)
式中,yin是级组进口工质的湿度,无量纲;hin是级组进口工质的比焓,kJ/kg;pin是级组进口工质压力,MPa;sin是级组进口工质比熵,kJ/(kg*K);hL(pin)和hV(pin)分别是级组进口工质压力对应的饱和水、饱和干蒸汽的比焓,kJ/kg;sL(pin)和sV(pin)分别是级组进口工质压力对应的饱和水、饱和干蒸汽的比熵,kJ/(kg*K);
(b)根据(a)中计算得到的级组进口工质比熵sin,计算在等熵膨胀时出口工质的湿度yexp,s和比焓hexp,s:
hexp,s=(1-yexp,s)hV(pout)+yexp,shL(pout)
式中,yexp,s是等熵膨胀时级组出口工质的湿度,无量纲;hexp,s是等熵膨胀时级组出口工质的比焓,kJ/kg;pout是级组出口工质压力,MPa;sL(pout)和sV(pout)分别是级组出口工质压力对应的饱和水、饱和干蒸汽的比熵,kJ/(kg*K);hL(pout)和hV(pout)分别是级组出口工质压力对应的饱和水、饱和干蒸汽的比焓,kJ/kg;
(c)计算级组的内效率:
ηexp=ηdry(1-αyin)
式中,ηexp是变工况下的级组内效率,无量纲;ηdry是级组的干基内效率,无量纲;yin是级组进口工质的湿度,无量纲;
(d)根据(b)得到的等熵膨胀时级组出口工质的比焓hexp,s和(c)中的级组内效率ηexp,计算蒸汽在级组内膨胀后的比焓及湿度:
hexp=hin-ηexp(hin-hexp,s)
式中,hexp是蒸汽在级组内膨胀后的比焓,kJ/kg;yexp是蒸汽在级组内膨胀后的湿度,无量纲;hL(pout)和hV(pout)分别是级组出口工质压力对应的饱和水、饱和干蒸汽的比焓,kJ/kg;
(e)根据(d)中得到的级组出口湿度yexp重新计算级组内效率ηexp:
(f)将(e)得到的级组内效率ηexp代回(d),重复(d)得到蒸汽在级组内膨胀后的比焓hexp及湿度yexp。
4.根据权利要求1所述的饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法,其特征在于,S4中,所述汽水分离后的饱和蒸汽比焓,采用以下公式计算得到:
式中,hsep是汽水分离后的饱和蒸汽比焓,kJ/kg;Din是级组进口工质质量流量,kg/s;hexp是饱和蒸汽在级组内膨胀后的比焓,kJ/kg;ηsp是级组的汽水分离效率,无量纲;yexp是蒸汽在级组内膨胀后的湿度,无量纲;pout是级组出口工质压力,MPa;hL(pout)是级组出口工质压力对应的饱和水的比焓,kJ/kg;
所述汽水分离后的饱和蒸汽湿度采用以下公式计算得到:
式中,ysep是汽水分离后的饱和蒸汽湿度,无量纲;hV(pout)是级组出口工质压力对应的饱和干蒸汽的比焓,kJ/kg;hL(pout)是级组出口工质压力对应的饱和水的比焓;
所述汽水分离后的饱和蒸汽质量流量,采用以下公式计算得到:
Dsep=(1-ηspyexp)Din
式中,Dsep是汽水分离后的饱和蒸汽质量流量,kg/s;Din是级组进口饱和蒸汽质量流量,kg/s;
将汽水分离后的饱和蒸汽比焓hsep作为级组出口饱和蒸汽的比焓hout;将汽水分离后的饱和蒸汽湿度ysep作为级组出口饱和蒸汽的湿度yout。
5.根据权利要求1所述的饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法,其特征在于,S5中,抽汽后饱和蒸汽的质量流量采用以下公式计算得到:
Dext=Dsep-Des
式中,Dext是抽汽后的饱和蒸汽的质量流量,kg/s;Dsep是汽水分离后的饱和蒸汽质量流量,kg/s;Des是抽汽质量流量,kg/s;
将抽汽后的饱和蒸汽的质量流量Dext作为级组出口饱和蒸汽的质量流量Dout。
6.根据权利要求5所述的饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法,其特征在于,所述抽汽质量流量Des通过求解以下方程组得到:
式中,VV是给水加热器内蒸汽体积,m3;ρV是饱和干蒸汽密度,kg/m3;hL和hV分别是饱和水、饱和干蒸汽比焓,kJ/kg;Des是抽汽质量流量,kg/s;Ddw,in是上级给水加热器疏水的质量流量,kg/s;hdw,in是上级给水加热器疏水的比焓,kJ/kg;Dcon是凝结水质量流量,kg/s;QV是蒸汽凝结段和疏水冷却段的热流量,kW;SV是给水加热器蒸汽凝结段的换热面积,m2;Tes和Tfw分别是抽汽温度和给水温度,K;Kfw和Kes分别是给水和蒸汽的表面传热系数,kW/(m2*K);wfw是给水的流速,m/s;ρfw是给水的密度,kg/m3,λfw和λL,es分别是给水的热导率和抽汽压力对应的饱和水热导率,kW/(m*K);cp,fw是给水的定压比热容,kW/(kg*K);μfw和μL,es分别是给水的动力粘度和抽汽压力对应的饱和水动力粘度,Pa·s;di是换热管内径,m;g是重力加速度,取9.8m/s2;ρL,es是抽汽压力对应的饱和水密度,kg/m3;do是给水加热器的换热管外径,m。
8.一种可读介质,其特征在于,所述可读介质上存储有指令,该指令在电子设备上执行时使电子设备执行权利要求1至7中任一项所述的饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法。
9.一种电子设备,包括存储器、处理器及存储在存储器上并能在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时用于执行权利要求1至7中任一项所述的饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210554966.4A CN114856726B (zh) | 2022-05-20 | 2022-05-20 | 一种饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210554966.4A CN114856726B (zh) | 2022-05-20 | 2022-05-20 | 一种饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114856726A CN114856726A (zh) | 2022-08-05 |
CN114856726B true CN114856726B (zh) | 2023-06-20 |
Family
ID=82640094
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210554966.4A Active CN114856726B (zh) | 2022-05-20 | 2022-05-20 | 一种饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114856726B (zh) |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2999122B2 (ja) * | 1994-06-24 | 2000-01-17 | 三菱重工業株式会社 | 複合プラントの制御装置 |
JP4155495B2 (ja) * | 2002-04-10 | 2008-09-24 | 株式会社東芝 | 原子炉隔離時冷却設備 |
CN103117101B (zh) * | 2013-01-19 | 2015-09-30 | 哈尔滨工程大学 | 用于一体化反应堆的启停辅助装置及一体化反应堆的冷启动方法 |
CN110160113A (zh) * | 2019-04-29 | 2019-08-23 | 中广核工程有限公司 | 一种核电厂供暖***及其供暖方法 |
CN112127958A (zh) * | 2020-08-28 | 2020-12-25 | 中核核电运行管理有限公司 | 一种确定核电汽轮机抽汽参数的装置及方法 |
-
2022
- 2022-05-20 CN CN202210554966.4A patent/CN114856726B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN114856726A (zh) | 2022-08-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100437015C (zh) | 汽轮机通流间隙变化在线监测方法 | |
CN108691585B (zh) | 一种凝汽式汽轮机低压缸效率的计算方法 | |
Pattanayak et al. | Thermal performance assessment of steam surface condenser | |
CN101825502B (zh) | 汽机带疏水冷却器的加热器出水及疏水温度测算方法 | |
CN111079302A (zh) | 一种低压缸效率测算***及方法 | |
CN109211439A (zh) | 一种汽轮机低压缸排汽焓值在线监测***及方法 | |
JP3965275B2 (ja) | 火力発電プラントの熱効率診断方法および装置 | |
CN211454603U (zh) | 一种低压缸效率测算*** | |
CN107543733A (zh) | 一种凝汽器热负荷在线监测***及测量方法 | |
CN105303037A (zh) | 一种引起凝汽器真空恶化的临界漏空气量的计算方法 | |
US5327772A (en) | Steam quality sensor | |
CN114856726B (zh) | 一种饱和蒸汽汽轮机乏汽容积流量及湿度实时确定方法 | |
CN104483152B (zh) | 非再热回热复合循环机组的热耗率测定方法 | |
JP4108097B2 (ja) | 原子力発電プラント熱効率診断システム、原子力発電プラント熱効率診断プログラムおよび原子力発電プラント熱効率診断方法 | |
CN101825503B (zh) | 汽机表面式加热器的出水及疏水温度测算方法 | |
CN112127957B (zh) | 一种火电厂汽轮机主蒸汽流量的测量方法 | |
CN108663216A (zh) | 一种汽轮机低压缸效率的直接测量方法 | |
CN112127958A (zh) | 一种确定核电汽轮机抽汽参数的装置及方法 | |
JP4924762B2 (ja) | 蒸気管の損失計測システム及び計測方法 | |
CN106894853B (zh) | 凝汽式汽轮机冷端节能诊断方法 | |
CN207623031U (zh) | 一种凝汽器热负荷在线监测*** | |
CN101852658B (zh) | 汽机带蒸冷和疏冷器加热器出水及疏水温度测算方法 | |
CN101832545B (zh) | 汽机带蒸汽冷却器的加热器出水及疏水温度测算方法 | |
CN113221477B (zh) | 一种确定循环水流量的热平衡计算方法 | |
CN112395817A (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 |