CN112417780B - 蒸汽发生器二回路再循环水质量流量估计方法及系统 - Google Patents

蒸汽发生器二回路再循环水质量流量估计方法及系统 Download PDF

Info

Publication number
CN112417780B
CN112417780B CN202011308103.6A CN202011308103A CN112417780B CN 112417780 B CN112417780 B CN 112417780B CN 202011308103 A CN202011308103 A CN 202011308103A CN 112417780 B CN112417780 B CN 112417780B
Authority
CN
China
Prior art keywords
working medium
section
channel
phase working
pressure
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
Application number
CN202011308103.6A
Other languages
English (en)
Other versions
CN112417780A (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.)
Shanghai Jiaotong University
Original Assignee
Shanghai Jiaotong University
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 Shanghai Jiaotong University filed Critical Shanghai Jiaotong University
Priority to CN202011308103.6A priority Critical patent/CN112417780B/zh
Publication of CN112417780A publication Critical patent/CN112417780A/zh
Application granted granted Critical
Publication of CN112417780B publication Critical patent/CN112417780B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08—Thermal analysis or thermal optimisation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Abstract

本发明提供了一种蒸汽发生器二回路再循环水质量流量估计方法及系统，获取给定时刻下蒸汽发生器的实时运行数据；计算一回路冷却剂与倒U型管金属壁之间的传热系数以及倒U型管金属壁与二回路工质之间的传热系数；建立下降通道模型，得到当前时刻下降通道底部出口液相工质的流量、温度及压力；建立一回路冷却剂模型，得到倒U型管金属壁的温度分布；建立上升通道模型，得到当前时刻沿倒U型管高度的二回路工质的流速、温度、压力分布；建立汽水分离器模型，得到汽水分离器出口气相工质、液相工质的温度、压力和质量流量，进而计算当前时刻蒸汽发生器的二回路再循环水质量流量。

Description

蒸汽发生器二回路再循环水质量流量估计方法及系统
技术领域
本发明属于核电站运行优化控制技术领域，具体地，给出了一种机理建模与DCS(运行机组现场分布式控制系统)实时测量数据相结合的核电蒸汽发生器二回路再循环水质量流量实时估计(辨识)方法及系统。
背景技术
蒸汽发生器是核电站的主要设备之一，它既是连接核电站一、二回路之间的枢纽，同时又起到了隔离一二回路的作用。在核反应堆中，核裂变产生的热量由一回路冷却剂带出，通过蒸汽发生器倒U型管将热量传递给二回路工质，使过冷状态水转变成饱和蒸汽。饱和蒸汽再流入汽轮机中做功，通过同轴发电机转换为电能。
自然循环蒸汽发生器内部结构复杂，在二次侧，由给水和再循环水混合成的混合水沿着倒U型管束冲刷倒U型管，被逐渐加热成饱和蒸汽。在蒸汽发生器内部存在着复杂的传热过程。例如在蒸汽发生器二回路侧，倒U型管向工质的热量传递包括单相对流换热、过冷沸腾以及饱和沸腾对流换热。在沸腾对流换热过程中，二回路工质局部汽化，形成气液两相流。气泡产生、成长、脱离壁面区域的过程强烈扰动二回路水位及传热阻力。而在蒸汽发生器一回路侧，冷却剂向倒U型管的热量传递为单相对流换热。由于蒸汽发生器系统的非线性、非对称性、时滞以及两相流换热过程的复杂性，目前国内外的相关研究以对蒸汽发生器集总参数建模和稳态性能仿真为主，对蒸汽发生器内部工质动态研究较少，故研究成果并不能用于改进蒸汽发生器结构设计和运行优化，也不利于提高蒸汽发生器实时液位控制系统的控制品质。
经过对现有技术检索后发现，授权公告号为CN201911135338.7，授权公告日为2020年03月17日的中国发明专利《一种立式蒸汽发生器热力性能评估系统和方法》，提出了一种建立反应堆U型管式蒸汽发生器热力评估系统，根据每一不同区域的温差传热或压降损失建立壳侧不同区域的壳侧下降通道熵产模型、壳侧管束上升区域熵产模型、汽水分离混合区域熵产熵产模型以及壳侧管束入口区域熵产模型。该专利主要贡献在于研究了蒸汽发生器内部工质的换热全过程，但未提及蒸汽发生器二回路再循环水质量流量计算方法。
综上，现有公开报道均未涉及核电蒸汽发生器二回路再循环水质量流量实时估计问题，这一空缺有待填补。
发明内容
针对现有技术的不足，本发明的目的是提供一种核电蒸汽发生器二回路再循环水质量流量实时估计方法及系统，实现了对现有蒸汽发生器二回路不可测量的再循环水质量流量的实时估计。
本发明第一方面，提供一种蒸汽发生器二回路再循环水质量流量估计方法，包括：
获取给定时刻下蒸汽发生器的实时运行数据；
利用获取的蒸汽发生器的实时运行数据，建立下降通道模型，求解下降通道模型，得到当前时刻下降通道底部出口液相工质的温度、压力及质量流量；
利用获取的蒸汽发生器的实时运行数据，建立一回路冷却剂模型，求解一回路冷却剂模型，得到冷却剂温度、压力、流速分布、一次侧传热系数分布以及倒U型管金属壁的温度分布；
利用获取的蒸汽发生器的实时运行数据、所述倒U型管金属壁的温度分布以及得到的下降通道底部出口液相工质的温度、压力及质量流量，建立上升通道模型，求解上升通道模型，得到当前时刻沿倒U型管高度的二回路工质的流速、温度、压力分布以及二次侧传热系数分布；所述当前时刻沿倒U型管高度的二回路工质的流速、温度以及压力分布包括上升通道顶部出口气液混合物工质的流速、温度以及压力；
利用获取的蒸汽发生器的实时运行数据以及得到的上升通道顶部出口气液混合物工质的流速、温度以及压力，建立汽水分离器模型，其中上升通道顶部出口气液混合物工质的流速、温度以及压力作为汽水分离器模型输入，解算汽水分离器模型得到汽水分离器出口液相工质的质量流量，即为当前时刻蒸汽发生器的二回路再循环水质量流量。
本发明第二方面，提供一种蒸汽发生器二回路再循环水质量流量估计装置，包括存储器、处理器及存储在存储器上并能在处理器上运行的计算机程序，所述处理器执行所述程序时用于执行上述的蒸汽发生器二回路再循环水质量流量估计方法。
与现有技术相比，本发明实施例具有以下至少一种有益效果：
本发明提供的上述提供一种蒸汽发生器二回路再循环水质量流量估计方法和装置，实现了蒸汽发生器的二回路再循环水质量流量估计，解决了现有技术的难题，填补了空白。
本发明提供的上述提供一种蒸汽发生器二回路再循环水质量流量估计方法和装置，具有实时性，所开发的蒸汽发生器的二回路再循环水质量流量实时估计平台刷新周期为秒级。
本发明提供的上述提供一种蒸汽发生器二回路再循环水质量流量估计方法和装置，可用于蒸汽发生器液位控制优化。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：
图1为本发明一优选实施例蒸汽发生器二回路再循环水质量流量实时估计方法示意图；
图2为本发明一优选实施例蒸汽发生器简化结构示意图；
图3为本发明一优选实施例中核电站机组输出负荷变化图；
图4为本发明一优选实施例中蒸汽发生器二回路再循环水质量流量估计结果；
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进。这些都属于本发明的保护范围。
本发明一实施例提供了一种核电蒸汽发生器二回路再循环水质量流量实时估计方法，以蒸汽发生器为研究对象，根据蒸汽发生器具体结构，将蒸汽发生器划分为热段、冷段和汽水分离器，其中热段和冷段又分别划分为下降通道和上升通道：下降通道是指外壳与内部套筒之间工质流经的空间，工质向下流动；上升通道是指内部套筒与倒U型管壁之间工质流经的空间，工质向上流动。基于质量、能量、动量动态衡算建立蒸汽发生器机理模型，完成模型验证，最后基于模型和DCS测量数据实时计算二回路再循环水质量流量。
参照图1所示，本实施例提供的核电蒸汽发生器二回路再循环水质量流量实时估计方法，包括以下步骤：
步骤一，获取给定时刻下蒸汽发生器的相关测点实时运行数据；
具体的，从运行机组现场分布式控制系统(DCS)实时数据库中获取给定时刻的运行数据，包括：机组负荷；给水温度、压力、质量流量；饱和蒸汽温度、压力、质量流量；一回路冷却剂进出口温度、压力、质量流量以及水位高度等。
步骤二，将蒸汽发生器划分为热段、冷段和汽水分离器，其中热段和冷段又分别划分为下降通道和上升通道：下降通道是指外壳与内部套筒之间工质流经的空间，工质向下流动；上升通道是指内部套筒与倒U型管壁之间工质流经的空间，工质向上流动；
步骤三，利用获取的蒸汽发生器的相关测点实时运行数据，根据工质的质量、能量、动量守恒关系，建立下降通道模型，得到当前时刻下降通道底部出口液相工质的温度、压力及质量流量。
步骤四，利用获取的蒸汽发生器的相关测点实时运行数据，计算一回路冷却剂与倒U型管金属壁之间传热系数以及倒U型管金属壁与二回路工质之间的传热系数；
步骤五，利用获取的蒸汽发生器的相关测点实时运行数据，根据工质的质量、能量、动量守恒关系，建立一回路冷却剂模型，得到倒U型管金属壁的温度分布；
步骤六，利用获取的蒸汽发生器的相关测点实时运行数据以及得到的下降通道底部出口液相工质的流量、温度、压力，结合工质的质量、能量、动量守恒关系，建立上升通道模型，得到当前时刻沿倒U型管高度的二回路工质的流速、温度以及压力分布；
步骤七，利用获取的蒸汽发生器的相关测点实时运行数据以及得到的上升通道顶部出口气液混合物工质的流速、温度以及压力，建立汽水分离器模型，计算得到汽水分离器出口质量流量，即为蒸汽发生器的二回路再循环水质量流量。
作为一优选实施例，在步骤二中，在蒸汽发生器上升通道模型中，根据二回路工质状态将上升通道划分为预热区和沸腾区。预热区和沸腾区分界面的划分依据为：
hRC(t,z)＝hsw(t,z) (1)
式中，hRC(t,z)是上升通道当前时刻t和高度z的二回路工质的比焓，kJ/kg；hsw(t,z)是当前时刻t和高度z的二回路工质饱和状态比焓，kJ/kg。
作为一优选实施例，在步骤三中，下降通道入口液相工质中占比
Figure BDA0002788904890000041
的给水流入热段，占比/>
Figure BDA0002788904890000051
的给水流入冷段，占比/>
Figure BDA0002788904890000052
的再循环水流入热段，占比/>
Figure BDA0002788904890000053
的再循环水流入冷段。根据动量、质量和能量守恒关系，建立热段下降通道模型为式(2)～(4)：/>
Figure BDA0002788904890000054
Figure BDA0002788904890000055
Figure BDA0002788904890000056
式中，MHL,DC是热段下降通道液相工质质量，kg；ρHL,DC是热段下降通道底部出口液相工质密度，kg/m3；AHL,DC是热段下降通道的横截面积，m2；H是下降通道的水位高度，m；Gfw是给水质量流量，kg/s；Grw是再循环水质量流量，kg/s；GHL,DC,out是热段下降通道底部出口液相工质质量流量，kg/s；CP,HL,DC是热段下降通道液相工质的定压比热容，kJ/(kg·K)；THL,DC是热段下降通道底部出口液相工质温度，K；hHL,DC是热段下降通道液相工质比焓，kJ/kg，可根据热段下降通道液相工质温度和压力通过工质物性参数数据库计算得到；hfw是给水比焓，kJ/kg，可根据给水温度和压力通过工质物性参数数据库计算得到；hrw是再循环水比焓，kJ/kg，可根据再循环水温度和压力通过工质物性参数数据库计算得到；hHL,DC，out是热段下降通道底部出口液相工质比焓，kJ/kg，可根据热段下降通道底部出口液相工质温度和压力通过工质物性参数数据库计算得到；PHL,DC热段下降通道底部出口液相工质压力，MPa；GHL,DC是热段下降通道液相工质质量流量，kg/s；fHL,DC是热段下降通道摩擦因子；De,HL,DC是热段下降通道当量直径，m；g是重力加速度，m/s2；
通过对热段下降通道模型求解，得到当前时刻热段下降通道底部出口液相工质的温度、压力以及质量流量。
建立冷段下降通道模型为式(5)～(7)：
Figure BDA0002788904890000057
Figure BDA0002788904890000058
Figure BDA0002788904890000061
式中，MCL,DC是冷段下降通道液相工质质量，kg；ρCL,DC是冷段下降通道液相工质密度，kg/m3；ACL,DC是冷段下降通道的横截面积，m2；GCL,DC,out是冷段下降通道底部出口液相工质质量流量，kg/s；CP,CL,DC是冷段下降通道液相工质的定压比热容，kJ/(kg·K)；TCL,DC是冷段下降通道液相工质温度，K；hCL,DC是冷段下降通道液相工质比焓，kJ/kg，可根据冷段下降通道液相工质温度和压力通过工质物性参数数据库计算得到；hCL,DC,out是冷段下降通道底部出口液相工质比焓，kJ/kg，可根据冷段下降通道底部出口液相工质温度和压力通过工质物性参数数据库计算得到；PCL,DC冷段下降通道液相工质压力，MPa；GCL,DC是冷段下降通道液相工质质量流量，kg/s；fCL,DC是冷段下降通道摩擦因子；De,CL,DC是冷段下降通道当量直径，m；
通过对冷段下降通道模型求解，得到当前时刻冷段下降通道底部出口液相工质的温度、压力以及质量流量。
作为一优选实施例，在步骤四中，热段和冷段一回路冷却剂与倒U型管金属壁之间传热系数KHL,PS和KCL,PS以及热段和冷段预热区倒U型管金属壁与二回路工质之间的传热系数KHL,RC,PR和KCL,RC,PR，采用迪图斯-贝尔公式计算：
K＝0.023Rew 0.8Prw 0.3λw/dHL,MT (8)
式中，Rew是工质雷诺数；Prw是工质普朗特数；λw是工质热导率；dHL,MT是倒U型管内径；
对热段和冷段沸腾区倒U型管金属壁与二回路工质之间的传热系数K*,RC,BR采用式(9)～(14)计算，其中，在热段的K*,RC,BR用KHL,RC,BR替换，在冷段的K*,RC,BR用KCL,RC,BR替换：
K*,RC,BR＝Kcht+Kbht (9)
Figure BDA0002788904890000062
Figure BDA0002788904890000063
Figure BDA0002788904890000071
Figure BDA0002788904890000072
Figure BDA0002788904890000073
式中，Kcht、Kbht分别是对流传热部分的传热系数和泡核沸腾传热部分的传热系数；CP,w是工质定压比热容；hfs是沸腾区液相工质汽化潜热；σ沸腾区液相工质表面张力系数；ΔTMT是沸腾区倒U型管金属壁过热度；ΔPMT是沸腾区饱和蒸汽压差；x是质量气含率；ρw是上升通道液相工质密度；ρs是上升通道饱和蒸汽密度；μw是上升通道液相工质粘性系数；μs是上升通道饱和蒸汽粘性系数；dHL,MT是热段倒U型管内径；G是工质质量流量；Xtt和S是中间变量。
作为一优选实施例，在步骤五中，考虑热段一回路冷却剂重力压降，根据动量、质量和能量守恒关系，建立热段一回路冷却剂模型为式(15)～(18)：
Figure BDA0002788904890000074
/>
Figure BDA0002788904890000075
Figure BDA0002788904890000076
Figure BDA0002788904890000077
式中，ρHL,PS是热段一回路冷却剂密度，kg/m3；WHL,PS是热段一回路冷却剂的流速，m/s；CP,HL,PS是热段一回路冷却剂定压比热容，kJ/(kg·K)；THL,PS是热段一回路冷却剂温度，K；KHL,PS是热段一回路冷却剂通过倒U型管金属壁向二回路工质传热的传热系数，kW/(m2·K)；dHL,MT是热段倒U型管内径，m；PHL,PS是热段一回路冷却剂压力，MPa。
通过对热段一回路冷却剂模型求解，得到热段倒U型管金属壁的温度分布。
建立冷段一回路冷却剂模型为式(19)～(22)：
Figure BDA0002788904890000081
Figure BDA0002788904890000082
Figure BDA0002788904890000083
Figure BDA0002788904890000084
式中，ρCL,PS是冷段一回路冷却剂密度，kg/m3；WCL,PS是冷段一回路冷却剂的流速，m/s；CP,CL,PS是冷段一回路冷却剂定压比热容，kJ/(kg·K)；TCL,PS是冷段一回路冷却剂温度，K；KCL,PS是冷段一回路冷却剂通过倒U型管金属壁向二回路工质传热的传热系数，kW/(m2·K)；dCL,MT是冷段倒U型管内径，m；PCL,PS是冷段一回路冷却剂压力，MPa。
通过对冷段一回路冷却剂模型求解，得到冷段倒U型管金属壁的温度分布。
作为一优选实施例，在步骤六中，考虑热段上升通道工质重力压降、摩擦压降和加速压降，根据动量、质量和能量守恒关系，建立热段上升通道模型为式(23)～(30)：
Figure BDA0002788904890000085
Figure BDA0002788904890000086
Figure BDA0002788904890000087
Figure BDA0002788904890000088
Figure BDA0002788904890000089
/>
Figure BDA00027889048900000810
Figure BDA0002788904890000091
Figure BDA0002788904890000092
式中，ρHL,RC是热段上升通道工质密度，kg/m3；WHL,RC是热段上升通道工质的流速，m/s；ρHL,MT是热段倒U型管金属壁密度，kg/m3；CP,HL,MT是热段倒U型管金属壁的定压比热容，kJ/(kg·K)；THL,MT是热段倒U型管金属壁温度，K；n是倒U型管根数；KHL,RC,PR是热段上升通道预热区二回路工质与倒U型管金属壁间传热系数，kW/(m2·K)；dHL,MT是热段倒U型管内径，m；THL,RC，PR是热段上升通道预热区液相工质温度，K；ρHL,RC，PR是热段上升通道预热区液相工质密度，kg/m3；CP,HL,RC,PR是热段上升通道预热区液相工质定压比热容，kJ/(kg·K)；WHL,RC，PR是热段上升通道预热区液相工质的流速，m/s；KHL,RC,BR是热段上升通道沸腾区二回路工质与倒U型管金属壁间传热系数，kW/(m2·K)；THL,RC，BR是热段上升通道沸腾区气液混合相工质温度，K；ρHL,RC，BR是热段上升通道沸腾区气液混合相工质密度，kg/m3；CP,HL,RC,BR是热段上升通道沸腾区气液混合相工质定压比热容，kJ/(kg·K)；WHL,RC，BR是热段上升通道沸腾区气液混合相工质的流速，m/s；PHL,RC,PR是热段上升通道预热区液相工质压力，MPa；GHL,RC,PR是热段上升通道预热区液相工质质量流量，kg/s；fHL,RC,PR是热段上升通道预热区摩擦因子；De,HL,RC,PR是热段上升通道预热区当量直径，m；ξHL,RC,PR是热段上升通道预热区局部阻力系数；PHL,RC,BR是热段上升通道沸腾区气液混合相工质压力，MPa；GHL,RC,BR是热段上升通道沸腾区气液混合相工质质量流量，kg/s；fHL,RC,BR是热段上升通道沸腾区摩擦因子；De,HL,RC,BR是热段上升通道沸腾区当量直径，m；φ是两相倍乘因子；ξHL,RC,BR是热段上升通道沸腾区局部阻力系数；x是质量气含率，％；ρw是上升通道液相工质密度，kg/m3；ρs是上升通道饱和蒸汽密度，kg/m3；μw是上升通道液相工质粘性系数；μs是上升通道饱和蒸汽粘性系数；
通过对热段上升通道模型求解，得到热段当前时刻沿倒U型管高度的二回路工质的流速、温度、压力以及传热系数分布。
建立冷段上升通道模型为式(31)～(37)：
Figure BDA0002788904890000093
Figure BDA0002788904890000101
/>
Figure BDA0002788904890000102
Figure BDA0002788904890000103
Figure BDA0002788904890000104
Figure BDA0002788904890000105
Figure BDA0002788904890000106
式中，ρCL,RC是冷段上升通道工质密度，kg/m3；WCL,RC是冷段上升通道工质的流速，m/s；ρCL,MT是冷段倒U型管金属壁密度，kg/m3；CP,CL,MT是冷段倒U型管金属壁的定压比热容，kJ/(kg·K)；TCL,MT是冷段倒U型管金属壁温度，K；KCL,RC,PR是冷段上升通道预热区二回路工质与倒U型管金属壁间传热系数，kW/(m2·K)；dCL,MT是冷段倒U型管内径，m；TCL,RC，PR是冷段上升通道预热区液相工质温度，K；ρCL,RC，PR是冷段上升通道预热区液相工质密度，kg/m3；CP,CL,RC,PR是冷段上升通道预热区液相工质定压比热容，kJ/(kg·K)；WCL,RC,PR是冷段上升通道预热区液相工质的流速，m/s；KCL,RC,BR是冷段上升通道沸腾区二回路工质与倒U型管金属壁间传热系数，kW/(m2·K)；TCL,RC,BR是冷段上升通道沸腾区气液混合相工质温度，K；ρCL,RC,BR是冷段上升通道沸腾区气液混合相工质密度，kg/m3；CP,CL,RC,BR是冷段上升通道沸腾区气液混合相工质定压比热容，kJ/(kg·K)；WCL,RC,BR是冷段上升通道沸腾区气液混合相工质的流速，m/s；PCL,RC,PR是冷段上升通道预热区液相工质压力，MPa；GCL,RC,PR是冷段上升通道预热区液相工质质量流量，kg/s；fCL,RC,PR是冷段上升通道预热区摩擦因子；De,CL,RC,PR是冷段上升通道预热区当量直径，m；ξCL,RC,PR是冷段上升通道预热区局部阻力系数；PCL,RC,BR是冷段上升通道沸腾区气液混合相工质压力，MPa；GCL,RC,BR是冷段上升通道沸腾区气液混合相工质质量流量，kg/s；fCL,RC,BR是冷段上升通道沸腾区摩擦因子；De,CL,RC,BR是冷段上升通道沸腾区当量直径，m；ξCL,RC,BR是冷段上升通道沸腾区局部阻力系数；
通过对冷段上升通道模型求解，得到冷段当前时刻沿倒U型管高度的二回路工质的流速、温度、压力以及传热系数分布。
进一步，计算当前时刻沿倒U型管高度的二回路工质质量气含率分布：
Figure BDA0002788904890000111
式中，hBR是沸腾区气液混合相工质比焓，kJ/kg；hss是沸腾区饱和蒸汽比焓，kJ/kg；hsw是沸腾区饱和水比焓，kJ/kg；xBR是沸腾区工质质量气含率，％。hBR,hsw和hsw可根据沸腾区气液混合相工质温度和压力通过工质物性参数数据库计算得到。
作为一优选实施例，在步骤六中，所述汽水分离器模型为式(39)～(46)：
Gss,SP,out＝(xHL,RC,BR,outGHL,RC,BR,out+xCL,RC,BR,outGCL,RC,BR,out)×η (39)
Gsw,SP,out＝(1-xHL,RC,BR,out×η)GHL,RC,BR,out+(1-xCL,RC,BR,out×η)GCL,RC,BR,out (40)
Figure BDA0002788904890000112
GSP,in＝GHL,RC,BR,out+GCL,RC,BR,out (42)
PSP,in＝PHL,RC,BR,out＝PCL,RC,BR,out (43)
TSP,in＝THL,RC,BR,out＝TCL,RC,BR,out (44)
PSP,out＝Pss,SP,out＝Tsw,SP,out (45)
TSP,in＝Tss,SP,out＝Tsw,SP,out (46)
式中，Gss,SP,out是汽水分离器出口饱和蒸汽质量流量，kg/s；xHL,RC,BR,out是热段上升通道沸腾区出口气液混合相工质质量气含率，％；GHL,RC,BR,out是热段上升通道沸腾区出口气液混合相工质质量流量，kg/s；xCL,RC,BR,out是冷段上升通道沸腾区出口气液混合相工质质量气含率，％；GCL,RC,BR,out是冷段上升通道沸腾区出口气液混合相工质质量流量，kg/s；η是汽水分离器效率，％；Gsw,SP,out是汽水分离器出口饱和水质量流量，kg/s；PSP,out是汽水分离器出口工质压力，MPa；Pss,SP,out是汽水分离器出口饱和蒸汽压力，MPa；Psw,SP,out是汽水分离器出口饱和水压力，MPa；PSP,in是汽水分离器入口气液混合相工质压力，MPa；TSP,in是汽水分离器入口气液混合相工质温度，K；TCL,RC,BR,out是冷段上升通道沸腾区出口气液混合相工质温度，K；THL,RC,BR,out是热段上升通道沸腾区出口气液混合相工质温度，K；ξSP是汽水分离器局部阻力系数；GSP,in是汽水分离器入口气液混合相工质质量流量，kg/s；ρSP,in是汽水分离器入口气液混合相工质密度，kg/m3；PHL,RC,BR,out是热段上升通道沸腾区出口气液混合相工质压力，MPa；PCL,RC,BR,out是冷段上升通道沸腾区出口气液混合相工质压力，MPa；Tss,SP,out是汽水分离器出口饱和蒸汽温度，K；Tsw,SP,out是汽水分离器出口饱和水温度，K；ρHL,RC,BR,out是热段上升通道沸腾区出口气液混合相工质密度，kg/m3；ρCL,RC,BR,out是冷段上升通道沸腾区出口气液混合相工质密度，kg/m3。
通过对汽水分离器模型求解，得到汽水分离器出口饱和水以及饱和蒸汽的温度、压力、质量流量，进而计算当前时刻蒸汽发生器的二回路再循环水质量流量。
下面结合附图，对本实施例提供的核电蒸汽发生器二回路再循环水质量流量实时估计方法进一步详细描述如下。
本实施例提供的方法，涉及核电站立式U型自循环蒸汽发生器全工况二回路再循环水质量流量实时估计。图1是实施例蒸汽发生器二回路再循环水质量流量实时估计方法示意图。利用从运行机组现场DCS分布式控制系统实时测量数据库中获取给定时刻下机组负荷，给水温度、压力、质量流量，饱和蒸汽温度、压力、质量流量，一回路冷却剂进出口温度、压力、质量流量以及水位高度等数据，结合工质物性参数数据库和蒸汽发生器结构参数库，解算热段模型、冷段模型和汽水分离器模型，输出再循环水质量流量。图2是蒸汽发生器简化结构示意图。根据蒸汽发生器真实结构，进行简化划分为热段、冷段和汽水分离器。热段和冷段的二回路又可以根据工质是否达到饱和状态划分为预热区和沸腾区。
具体的，在一实施例中，上述方法包括以下步骤：
步骤一、从运行机组现场分布式控制系统(DCS)实测数据库中获取给定时刻的运行数据包括：机组负荷；给水温度、压力、质量流量；饱和蒸汽温度、压力、质量流量；一回路冷却剂进出口温度、压力、质量流量以及水位高度等；
步骤二、在蒸汽发生器上升通道模型中，根据二回路工质状态将上升通道划分为预热区和沸腾区。预热区和沸腾区分界面的划分依据为：
hRC(t,z)＝hsw(t,z) (1)
式中，hRC(t,z)是上升通道当前时刻t和高度z的二回路工质的比焓，kJ/kg；hsw(t,z)是当前时刻t和高度z的二回路工质饱和状态比焓，kJ/kg。
步骤三、下降通道入口液相工质中占比80％的给水流入热段，占比20％的给水流入冷段，占比50％的再循环水流入热段，占比50％的再循环水流入冷段。求解热段下降通道模型，式(2)～(4)：
Figure BDA0002788904890000131
Figure BDA0002788904890000132
Figure BDA0002788904890000133
式中，MHL,DC是热段下降通道液相工质质量，kg；ρHL,DC是热段下降通道底部出口液相工质密度，kg/m3；AHL,DC是热段下降通道的横截面积，m2；H是下降通道的水位高度，m；Gfw是给水质量流量，kg/s；Grw是再循环水质量流量，kg/s；GHL,DC,out是热段下降通道底部出口液相工质质量流量，kg/s；CP,HL,DC是热段下降通道液相工质的定压比热容，kJ/(kg·K)；THL,DC是热段下降通道底部出口液相工质温度，K；hHL,DC是热段下降通道液相工质比焓，kJ/kg，可根据热段下降通道液相工质温度和压力通过工质物性参数数据库计算得到；hfw是给水比焓，kJ/kg，可根据给水温度和压力通过工质物性参数数据库计算得到；hrw是再循环水比焓，kJ/kg，可根据再循环水温度和压力通过工质物性参数数据库计算得到；hHL,DC，out是热段下降通道底部出口液相工质比焓，kJ/kg，可根据热段下降通道底部出口液相工质温度和压力通过工质物性参数数据库计算得到；PHL,DC热段下降通道底部出口液相工质压力，MPa；GHL,DC是热段下降通道液相工质质量流量，kg/s；fHL,DC是热段下降通道摩擦因子；De,HL,DC是热段下降通道当量直径，m；g是重力加速度，m/s2；
通过对热段下降通道模型求解，得到当前时刻热段下降通道底部出口液相工质的温度、压力以及质量流量。
求解冷段下降通道模型，式(5)～(7)：
Figure BDA0002788904890000134
Figure BDA0002788904890000135
/>
Figure BDA0002788904890000136
式中，MCL,DC是冷段下降通道液相工质质量，kg；ρCL,DC是冷段下降通道液相工质密度，kg/m3；ACL,DC是冷段下降通道的横截面积，m2；GCL,DC,out是冷段下降通道底部出口液相工质质量流量，kg/s；CP,CL,DC是冷段下降通道液相工质的定压比热容，kJ/(kg·K)；TCL,DC是冷段下降通道液相工质温度，K；hCL,DC是冷段下降通道液相工质比焓，kJ/kg，可根据冷段下降通道液相工质温度和压力通过工质物性参数数据库计算得到；hCL,DC,out是冷段下降通道底部出口液相工质比焓，kJ/kg，可根据冷段下降通道底部出口液相工质温度和压力通过工质物性参数数据库计算得到；PCL,DC冷段下降通道液相工质压力，MPa；GCL,DC是冷段下降通道液相工质质量流量，kg/s；fCL,DC是冷段下降通道摩擦因子；De,CL,DC是冷段下降通道当量直径，m；
通过对冷段下降通道模型求解，得到当前时刻冷段下降通道底部出口液相工质的温度、压力以及质量流量。
步骤四、热段和冷段一回路冷却剂与倒U型管金属壁之间传热系数KHL,PS和KCL,PS以及热段和冷段预热区倒U型管金属壁与二回路工质之间的传热系数KHL,RC,PR和KCL,RC,PR，采用迪图斯-贝尔公式计算：
K＝0.023Rew 0.8Prw 0.3λw/dHL,MT (8)
式中，Rew是对应的热段或冷段一回路或二回路工质雷诺数；Prw是对应的热段或冷段一回路或二回路工质普朗特数；λw是对应的热段或冷段一回路或二回路工质热导率，kW/(m·K)；dHL,MT是倒U型管内径，m；
对热段和冷段沸腾区倒U型管金属壁与二回路工质之间的传热系数K*,RC,BR采用式(9)～(14)计算，其中，在热段的K*,RC,BR用KHL,RC,BR替换，在冷段的K*,RC,BR用KCL,RC,BR替换：
K*,RC,BR＝Kcht+Kbht (9)
Figure BDA0002788904890000141
Figure BDA0002788904890000142
Figure BDA0002788904890000143
Figure BDA0002788904890000151
/>
Figure BDA0002788904890000152
式中，Kcht、Kbht分别是对流传热部分的传热系数和泡核沸腾传热部分的传热系数；CP,w是工质定压比热容；hfs是沸腾区液相工质汽化潜热；σ沸腾区液相工质表面张力系数；ΔTMT是沸腾区倒U型管金属壁过热度；ΔPMT是沸腾区饱和蒸汽压差；x是质量气含率；ρw是上升通道液相工质密度；ρs是上升通道饱和蒸汽密度；μw是上升通道液相工质粘性系数；μs是上升通道饱和蒸汽粘性系数；dHL,MT是热段倒U型管内径；G是工质质量流量；Xtt和S是中间变量。
步骤五、考虑热段一回路冷却剂重力压降，根据动量、质量和能量守恒关系，求解热段一回路冷却剂模型，式(15)～(18)：
Figure BDA0002788904890000153
Figure BDA0002788904890000154
Figure BDA0002788904890000155
Figure BDA0002788904890000156
式中，ρHL,PS是热段一回路冷却剂密度，kg/m3；WHL,PS是热段一回路冷却剂的流速，m/s；CP,HL,PS是热段一回路冷却剂定压比热容，kJ/(kg·K)；THL,PS是热段一回路冷却剂温度，K；KHL,PS是热段一回路冷却剂通过倒U型管金属壁向二回路工质传热的传热系数，kW/(m2·K)；dHL,MT是热段倒U型管内径，m；PHL,PS是热段一回路冷却剂压力，MPa。
通过对热段一回路冷却剂模型求解，得到热段倒U型管金属壁的温度分布。
求解冷段一回路冷却剂模型，式(19)～(22)：
Figure BDA0002788904890000157
Figure BDA0002788904890000161
Figure BDA0002788904890000162
Figure BDA0002788904890000163
式中，ρCL,PS是冷段一回路冷却剂密度，kg/m3；WCL,PS是冷段一回路冷却剂的流速，m/s；CP,CL,PS是冷段一回路冷却剂定压比热容，kJ/(kg·K)；TCL,PS是冷段一回路冷却剂温度，K；KCL,PS是冷段一回路冷却剂通过倒U型管金属壁向二回路工质传热的传热系数，kW/(m2·K)；dCL,MT是冷段倒U型管内径，m；PCL,PS是冷段一回路冷却剂压力，MPa。
通过对冷段一回路冷却剂模型求解，得到冷段倒U型管金属壁的温度分布。
步骤六、考虑热段上升通道工质重力压降、摩擦压降和加速压降，根据动量、质量和能量守恒关系，求解热段上升通道模型，式(23)～(30)：
Figure BDA0002788904890000164
Figure BDA0002788904890000165
Figure BDA0002788904890000166
Figure BDA0002788904890000167
Figure BDA0002788904890000168
Figure BDA0002788904890000169
Figure BDA00027889048900001610
Figure BDA0002788904890000171
式中，ρHL,RC是热段上升通道工质密度，kg/m3；WHL,RC是热段上升通道工质的流速，m/s；ρHL,MT是热段倒U型管金属壁密度，kg/m3；CP,HL,MT是热段倒U型管金属壁的定压比热容，kJ/(kg·K)；THL,MT是热段倒U型管金属壁温度，K；n是倒U型管根数；KHL,RC,PR是热段上升通道预热区二回路工质与倒U型管金属壁间传热系数，kW/(m2·K)；dHL,MT是热段倒U型管内径，m；THL,RC，PR是热段上升通道预热区液相工质温度，K；ρHL,RC，PR是热段上升通道预热区液相工质密度，kg/m3；CP,HL,RC,PR是热段上升通道预热区液相工质定压比热容，kJ/(kg·K)；WHL,RC，PR是热段上升通道预热区液相工质的流速，m/s；KHL,RC,BR是热段上升通道沸腾区二回路工质与倒U型管金属壁间传热系数，kW/(m2·K)；THL,RC，BR是热段上升通道沸腾区气液混合相工质温度，K；ρHL,RC，BR是热段上升通道沸腾区气液混合相工质密度，kg/m3；CP,HL,RC,BR是热段上升通道沸腾区气液混合相工质定压比热容，kJ/(kg·K)；WHL,RC，BR是热段上升通道沸腾区气液混合相工质的流速，m/s；PHL,RC,PR是热段上升通道预热区液相工质压力，MPa；GHL,RC,PR是热段上升通道预热区液相工质质量流量，kg/s；fHL,RC,PR是热段上升通道预热区摩擦因子；De,HL,RC,PR是热段上升通道预热区当量直径，m；ξHL,RC,PR是热段上升通道预热区局部阻力系数；PHL,RC,BR是热段上升通道沸腾区气液混合相工质压力，MPa；GHL,RC,BR是热段上升通道沸腾区气液混合相工质质量流量，kg/s；fHL,RC,BR是热段上升通道沸腾区摩擦因子；De,HL,RC,BR是热段上升通道沸腾区当量直径，m；φ是两相倍乘因子；ξHL,RC,BR是热段上升通道沸腾区局部阻力系数；x是质量气含率，％；ρw是上升通道液相工质密度，kg/m3；ρs是上升通道饱和蒸汽密度，kg/m3；μw是上升通道液相工质粘性系数；μs是上升通道饱和蒸汽粘性系数；
通过对热段上升通道模型求解，得到热段当前时刻沿倒U型管高度的二回路工质的流速、温度、压力以及传热系数分布。
求解冷段上升通道模型，式(31)～(37)：
Figure BDA0002788904890000172
Figure BDA0002788904890000173
Figure BDA0002788904890000174
Figure BDA0002788904890000181
Figure BDA0002788904890000182
Figure BDA0002788904890000183
Figure BDA0002788904890000184
式中，ρCL,RC是冷段上升通道工质密度，kg/m3；WCL,RC是冷段上升通道工质的流速，m/s；ρCL,MT是冷段倒U型管金属壁密度，kg/m3；CP,CL,MT是冷段倒U型管金属壁的定压比热容，kJ/(kg·K)；TCL,MT是冷段倒U型管金属壁温度，K；KCL,RC,PR是冷段上升通道预热区二回路工质与倒U型管金属壁间传热系数，kW/(m2·K)；dCL,MT是冷段倒U型管内径，m；TCL,RC，PR是冷段上升通道预热区液相工质温度，K；ρCL,RC，PR是冷段上升通道预热区液相工质密度，kg/m3；CP,CL,RC,PR是冷段上升通道预热区液相工质定压比热容，kJ/(kg·K)；WCL,RC,PR是冷段上升通道预热区液相工质的流速，m/s；KCL,RC,BR是冷段上升通道沸腾区二回路工质与倒U型管金属壁间传热系数，kW/(m2·K)；TCL,RC,BR是冷段上升通道沸腾区气液混合相工质温度，K；ρCL,RC,BR是冷段上升通道沸腾区气液混合相工质密度，kg/m3；CP,CL,RC,BR是冷段上升通道沸腾区气液混合相工质定压比热容，kJ/(kg·K)；WCL,RC,BR是冷段上升通道沸腾区气液混合相工质的流速，m/s；PCL,RC,PR是冷段上升通道预热区液相工质压力，MPa；GCL,RC,PR是冷段上升通道预热区液相工质质量流量，kg/s；fCL,RC,PR是冷段上升通道预热区摩擦因子；De,CL,RC,PR是冷段上升通道预热区当量直径，m；ξCL,RC,PR是冷段上升通道预热区局部阻力系数；PCL,RC,BR是冷段上升通道沸腾区气液混合相工质压力，MPa；GCL,RC,BR是冷段上升通道沸腾区气液混合相工质质量流量，kg/s；fCL,RC,BR是冷段上升通道沸腾区摩擦因子；De,CL,RC,BR是冷段上升通道沸腾区当量直径，m；ξCL,RC,BR是冷段上升通道沸腾区局部阻力系数；
通过对冷段上升通道模型求解，得到冷段当前时刻沿倒U型管高度的二回路工质的流速、温度、压力以及传热系数分布。
进一步，计算当前时刻沿倒U型管高度的二回路工质质量气含率分布：
Figure BDA0002788904890000191
式中，hBR是沸腾区气液混合相工质比焓，kJ/kg；hss是沸腾区饱和蒸汽比焓，kJ/kg；hsw是沸腾区饱和水比焓，kJ/kg；xBR是沸腾区工质质量气含率，％。hBR,hsw和hsw可根据沸腾区气液混合相工质温度和压力通过工质物性参数数据库计算得到。
步骤七，求解汽水分离器模型，式(39)～(46)：
Gss,SP,out＝(xHL,RC,BR,outGHL,RC,BR,out+xCL,RC,BR,outGCL,RC,BR,out)×η (39)
Gsw,SP,out＝(1-xHL,RC,BR,out×η)GHL,RC,BR,out+(1-xCL,RC,BR,out×η)GCL,RC,BR,out (40)
Figure BDA0002788904890000192
GSP,in＝GHL,RC,BR,out+GCL,RC,BR,out (42)
PSP,in＝PHL,RC,BR,out＝PCL,RC,BR,out (43)
TSP,in＝THL,RC,BR,out＝TCL,RC,BR,out (44)
PSP,out＝Pss,SP,out＝Tsw,SP,out (45)
TSP,in＝Tss,SP,out＝Tsw,SP,out (46)
式中，Gss,SP,out是汽水分离器出口饱和蒸汽质量流量，kg/s；xHL,RC,BR,out是热段上升通道沸腾区出口气液混合相工质质量气含率，％；GHL,RC,BR,out是热段上升通道沸腾区出口气液混合相工质质量流量，kg/s；xCL,RC,BR,out是冷段上升通道沸腾区出口气液混合相工质质量气含率，％；GCL,RC,BR,out是冷段上升通道沸腾区出口气液混合相工质质量流量，kg/s；η是汽水分离器效率，％；Gsw,SP,out是汽水分离器出口饱和水质量流量，kg/s；PSP,out是汽水分离器出口工质压力，MPa；Pss,SP,out是汽水分离器出口饱和蒸汽压力，MPa；Psw,SP,out是汽水分离器出口饱和水压力，MPa；PSP,in是汽水分离器入口气液混合相工质压力，MPa；TSP,in是汽水分离器入口气液混合相工质温度，K；TCL,RC,BR,out是冷段上升通道沸腾区出口气液混合相工质温度，K；THL,RC,BR,out是热段上升通道沸腾区出口气液混合相工质温度，K；ξSP是汽水分离器局部阻力系数；GSP,in是汽水分离器入口气液混合相工质质量流量，kg/s；ρSP,in是汽水分离器入口气液混合相工质密度，kg/m3；PHL,RC,BR,out是热段上升通道沸腾区出口气液混合相工质压力，MPa；PCL,RC,BR,out是冷段上升通道沸腾区出口气液混合相工质压力，MPa；Tss,SP,out是汽水分离器出口饱和蒸汽温度，K；Tsw,SP,out是汽水分离器出口饱和水温度，K；ρHL,RC,BR,out是热段上升通道沸腾区出口气液混合相工质密度，kg/m3；ρCL,RC,BR,out是冷段上升通道沸腾区出口气液混合相工质密度，kg/m3。
通过对汽水分离器模型求解，得到汽水分离器出口饱和水以及饱和蒸汽的温度、压力、质量流量，进而计算当前时刻蒸汽发生器的二回路再循环水质量流量。
本实施例提供的核电蒸汽发生器二回路再循环水质量流量实时估计方法，如图3所示，获取实施例核电站机组2019年6月10日不同负荷下的蒸汽发生器DCS实测数据。图4是蒸汽发生器二回路再循环水质量流量估计的仿真结果。由图4可知，在沸腾区再循环水质量流量随着核电站机组负荷降低而降低，反之亦反。
在本发明另一实施例中，提供一种蒸汽发生器二回路再循环水质量流量估计装置，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时用于执行上述的蒸汽发生器二回路再循环水质量流量估计方法。
本发明实施例利用获取的蒸汽发生器的相关测点实时运行数据以及得到的上升通道顶部出口气液混合物工质的流速、温度以及压力，建立汽水分离器模型，解算汽水分离器模型得到汽水分离器出口液相工质的质量流量，即为当前时刻蒸汽发生器的二回路再循环水质量流量。其中汽水分离器出口液相工质就是再循环水，它在实际工程中无测量值，本发明实现了核电蒸汽发生器二回路再循环水质量流量的实时估计。
本发明上述实施例提供的核电蒸汽发生器二回路再循环水质量流量实时估计方法及系统，其中，再循环水指流出汽水分离器的饱和水。本发明上述实施例获取给定时刻下蒸汽发生器相关测点的实时测量数据；将蒸汽发生器划分为热段、冷段和汽水分离器，其中热段和冷段又划分为下降通道和上升通道：下降通道是指蒸汽发生器外壳与内部套筒之间二回路工质流经的空间，上升通道是指内部套筒与倒U型管管壁之间二回路工质流经的空间；求解下降通道模型得到当前时刻下降通道底部出口二回路工质的比焓、流量、温度、压力以及密度；求解上升通道模型得到当前时刻沿倒U型管高度的二回路工质的比焓、流量、温度、压力以及密度分布；求解汽水分离器模型计算得到汽水分离器出口工质温度、压力和流量，进而计算出当前时刻蒸汽发生器二回路再循环水质量流量。本发明上述实施例提供的技术方案，实现了全工况蒸汽发生器二回路再循环水质量流量实时估计，可用于了解二回路工质流动动态特性，为蒸汽发生器运行优化及监测提供支撑条件，有助于提高核电站运行的安全性与经济性。
需要说明的是，本发明提供的方法中的步骤，可以利用系统中对应的模块、装置、单元等予以实现，本领域技术人员可以参照方法的技术方案实现系统的组成，即，方法中的实施例可理解为构建系统的优选例，在此不予赘述。
以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变形或修改，这并不影响本发明的实质内容。

Claims (10)

1.一种蒸汽发生器二回路再循环水质量流量估计方法，其特征在于，包括：
获取给定时刻下蒸汽发生器的实时运行数据；
利用获取的蒸汽发生器的实时运行数据，建立下降通道模型，求解下降通道模型，得到当前时刻下降通道底部出口液相工质的温度、压力及质量流量；
利用获取的蒸汽发生器的实时运行数据，建立一回路冷却剂模型，求解一回路冷却剂模型，得到冷却剂温度、压力、流速分布、一次侧传热系数分布以及倒U型管金属壁的温度分布；
利用获取的蒸汽发生器的实时运行数据、所述倒U型管金属壁的温度分布以及得到的下降通道底部出口液相工质的温度、压力及质量流量，建立上升通道模型，求解上升通道模型，得到当前时刻沿倒U型管高度的二回路工质的流速、温度、压力分布以及二次侧传热系数分布；所述当前时刻沿倒U型管高度的二回路工质的流速、温度以及压力分布包括上升通道顶部出口气液混合物工质的流速、温度以及压力；
利用获取的蒸汽发生器的实时运行数据以及得到的上升通道顶部出口气液混合物工质的流速、温度以及压力，建立汽水分离器模型，其中上升通道顶部出口气液混合物工质的流速、温度以及压力作为汽水分离器模型输入，解算汽水分离器模型得到汽水分离器出口液相工质的质量流量，即为当前时刻蒸汽发生器的二回路再循环水质量流量。
2.根据权利要求1所述的蒸汽发生器二回路再循环水质量流量估计方法，其特征在于，所述给定时刻下蒸汽发生器的相关测点实时运行数据，包括：
-机组负荷；
-给水温度、压力以及质量流量；
-饱和蒸汽温度、压力以及质量流量；
-一回路冷却剂进出口温度、压力以及质量流量；
-水位高度。
3.根据权利要求1所述的蒸汽发生器二回路再循环水质量流量估计方法，其特征在于，在所述上升通道中，根据二回路工质状态将上升通道划分为预热区和沸腾区；其中，预热区和沸腾区分界面的划分依据为：
hRC(t,z)＝hsw(t,z) (1)
式中，hRC(t,z)是上升通道当前时刻t和高度z的二回路工质的比焓；hsw(t,z)是当前时刻t和高度z的二回路工质饱和状态比焓。
4.根据权利要求1所述的蒸汽发生器二回路再循环水质量流量估计方法，其特征在于，所述下降通道入口液相工质中，占比
Figure FDA0003980540470000021
的给水流入热段，占比
Figure FDA0003980540470000022
的给水流入冷段，占比
Figure FDA0003980540470000023
的再循环水流入热段，占比
Figure FDA0003980540470000024
的再循环水流入冷段，
Figure FDA0003980540470000025
的取值范围为：70～90；
Figure FDA0003980540470000026
的取值范围为：40～60；根据动量、质量和能量守恒关系，建立热段下降通道模型为式(2)～(4)：
Figure FDA0003980540470000027
Figure FDA0003980540470000028
Figure FDA0003980540470000029
式中，MHL,DC是热段下降通道液相工质质量；ρHL,DC是热段下降通道底部出口液相工质密度；AHL,DC是热段下降通道的横截面积；H是下降通道的水位高度；Gfw是给水质量流量；Grw是再循环水质量流量；GHL,DC,out是热段下降通道底部出口液相工质质量流量；CP,HL,DC是热段下降通道液相工质的定压比热容；THL,DC是热段下降通道底部出口液相工质温度；hHL,DC是热段下降通道液相工质比焓，根据热段下降通道液相工质温度和压力通过工质物性参数数据库计算得到；hfw是给水比焓，根据给水温度和压力通过工质物性参数数据库计算得到；hrw是再循环水比焓，根据再循环水温度和压力通过工质物性参数数据库计算得到；hHL,DC，out是热段下降通道底部出口液相工质比焓，根据热段下降通道底部出口液相工质温度和压力通过工质物性参数数据库计算得到；PHL,DC热段下降通道底部出口液相工质压力；GHL,DC是热段下降通道液相工质质量流量；fHL,DC是热段下降通道摩擦因子；De,HL,DC是热段下降通道当量直径；g是重力加速度；
通过对热段下降通道模型求解，得到当前时刻热段下降通道底部出口液相工质的温度、压力以及质量流量；
进一步，建立冷段下降通道模型为式(5)～(7)：
Figure FDA0003980540470000031
Figure FDA0003980540470000032
Figure FDA0003980540470000033
式中，MCL,DC是冷段下降通道液相工质质量；ρCL,DC是冷段下降通道液相工质密度；ACL,DC是冷段下降通道的横截面积；GCL,DC,out是冷段下降通道底部出口液相工质质量流量；CP,CL,DC是冷段下降通道液相工质的定压比热容；TCL,DC是冷段下降通道液相工质温度；hCL,DC是冷段下降通道液相工质比焓，根据冷段下降通道液相工质温度和压力通过工质物性参数数据库计算得到；hCL,DC,out是冷段下降通道底部出口液相工质比焓，根据冷段下降通道底部出口液相工质温度和压力通过工质物性参数数据库计算得到；PCL,DC冷段下降通道液相工质压力；GCL,DC是冷段下降通道液相工质质量流量；fCL,DC是冷段下降通道摩擦因子；De,CL,DC是冷段下降通道当量直径；
通过对冷段下降通道模型求解，得到当前时刻冷段下降通道底部出口液相工质的温度、压力以及质量流量。
5.根据权利要求1所述的蒸汽发生器二回路再循环水质量流量估计方法，其特征在于，热段和冷段一回路冷却剂与倒U型管金属壁之间传热系数KHL,PS和KCL,PS以及热段和冷段预热区倒U型管金属壁与二回路工质之间的传热系数KHL,RC,PR和KCL,RC,PR，采用迪图斯-贝尔公式计算：
K＝0.023Rew 0.8Prw 0.3λw/dHL,MT (8)
式中，Rew是对应的热段或冷段一回路或二回路工质雷诺数；Prw是对应的热段或冷段一回路或二回路工质普朗特数；λw是对应的热段或冷段一回路或二回路工质热导率；dHL,MT是倒U型管内径；
对热段和冷段沸腾区倒U型管金属壁与二回路工质之间的传热系数K*,RC,BR采用式(9)～(14)计算，其中，在热段的K*,RC,BR用KHL,RC,BR替换，在冷段的K*,RC,BR用KCL,RC,BR替换：
K*,RC,BR＝Kcht+Kbht (9)
Figure FDA0003980540470000041
Figure FDA0003980540470000042
Figure FDA0003980540470000043
Figure FDA0003980540470000044
Figure FDA0003980540470000045
式中，Kcht、Kbht分别是对流传热部分的传热系数和泡核沸腾传热部分的传热系数；CP,w是工质定压比热容；hfs是沸腾区液相工质汽化潜热；σ沸腾区液相工质表面张力系数；△TMT是沸腾区倒U型管金属壁过热度；△PMT是沸腾区饱和蒸汽压差；x是质量气含率；ρw是上升通道液相工质密度；ρs是上升通道饱和蒸汽密度；μw是上升通道液相工质粘性系数；μs是上升通道饱和蒸汽粘性系数；dHL,MT是热段倒U型管内径；G是工质质量流量；Xtt和S是中间变量。
6.根据权利要求1所述的蒸汽发生器二回路再循环水质量流量估计方法，其特征在于，考虑热段一回路冷却剂重力压降，根据动量、质量和能量守恒关系，建立热段一回路冷却剂模型为式(15)～(18)：
Figure FDA0003980540470000046
Figure FDA0003980540470000047
Figure FDA0003980540470000048
Figure FDA0003980540470000049
式中，ρHL,PS是热段一回路冷却剂密度；WHL,PS是热段一回路冷却剂的流速；CP,HL,PS是热段一回路冷却剂定压比热容；THL,PS是热段一回路冷却剂温度；KHL,PS是热段一回路冷却剂通过倒U型管金属壁向二回路工质传热的传热系数；dHL,MT是热段倒U型管内径；PHL,PS是热段一回路冷却剂压力；
通过对热段一回路冷却剂模型求解，得到热段倒U型管金属壁的温度分布；
建立冷段一回路冷却剂模型为式(19)～(22)：
Figure FDA0003980540470000051
Figure FDA0003980540470000052
Figure FDA0003980540470000053
Figure FDA0003980540470000054
式中，ρCL,PS是冷段一回路冷却剂密度；WCL,PS是冷段一回路冷却剂的流速；CP,CL,PS是冷段一回路冷却剂定压比热容；TCL,PS是冷段一回路冷却剂温度；KCL,PS是冷段一回路冷却剂通过倒U型管金属壁向二回路工质传热的传热系数；dCL,MT是冷段倒U型管内径；PCL,PS是冷段一回路冷却剂压力；
通过对冷段一回路冷却剂模型求解，得到冷段倒U型管金属壁的温度分布。
7.根据权利要求1所述的蒸汽发生器二回路再循环水质量流量估计方法，其特征在于，考虑热段上升通道工质重力压降、摩擦压降和加速压降，根据动量、质量和能量守恒关系，建立热段上升通道模型为式(23)～(30)：
Figure FDA0003980540470000055
Figure FDA0003980540470000056
Figure FDA0003980540470000057
Figure FDA0003980540470000058
Figure FDA0003980540470000059
Figure FDA0003980540470000061
Figure FDA0003980540470000062
Figure FDA0003980540470000063
式中，ρHL,RC是热段上升通道工质密度；WHL,RC是热段上升通道工质的流速；ρHL,MT是热段倒U型管金属壁密度；CP,HL,MT是热段倒U型管金属壁的定压比热容；THL,MT是热段倒U型管金属壁温度；n是倒U型管根数；KHL,RC,PR是热段上升通道预热区二回路工质与倒U型管金属壁间传热系数；dHL,MT是热段倒U型管内径；THL,RC，PR是热段上升通道预热区液相工质温度；ρHL,RC，PR是热段上升通道预热区液相工质密度；CP,HL,RC,PR是热段上升通道预热区液相工质定压比热容；WHL,RC，PR是热段上升通道预热区液相工质的流速；KHL,RC,BR是热段上升通道沸腾区二回路工质与倒U型管金属壁间传热系数；THL,RC，BR是热段上升通道沸腾区气液混合相工质温度；ρHL,RC，BR是热段上升通道沸腾区气液混合相工质密度；CP,HL,RC,BR是热段上升通道沸腾区气液混合相工质定压比热容；WHL,RC，BR是热段上升通道沸腾区气液混合相工质的流速；PHL,RC,PR是热段上升通道预热区液相工质压力a；GHL,RC,PR是热段上升通道预热区液相工质质量流量；fHL,RC,PR是热段上升通道预热区摩擦因子；De,HL,RC,PR是热段上升通道预热区当量直径；ξHL,RC,PR是热段上升通道预热区局部阻力系数；PHL,RC,BR是热段上升通道沸腾区气液混合相工质压力；GHL,RC,BR是热段上升通道沸腾区气液混合相工质质量流量；fHL,RC,BR是热段上升通道沸腾区摩擦因子；De,HL,RC,BR是热段上升通道沸腾区当量直径；φ是两相倍乘因子；ξHL,RC,BR是热段上升通道沸腾区局部阻力系数；x是质量气含率；ρw是上升通道液相工质密度；ρs是上升通道饱和蒸汽密度；μw是上升通道液相工质粘性系数；μs是上升通道饱和蒸汽粘性系数；
通过对热段上升通道模型求解，得到热段当前时刻沿倒U型管高度的二回路工质的流速、温度、压力以及传热系数分布；
建立冷段上升通道模型为式(31)～(37)：
Figure FDA0003980540470000071
Figure FDA0003980540470000072
Figure FDA0003980540470000073
Figure FDA0003980540470000074
Figure FDA0003980540470000075
Figure FDA0003980540470000076
Figure FDA0003980540470000077
式中，ρCL,RC是冷段上升通道工质密度；WCL,RC是冷段上升通道工质的流速；ρCL,MT是冷段倒U型管金属壁密度；CP,CL,MT是冷段倒U型管金属壁的定压比热容；TCL,MT是冷段倒U型管金属壁温度；KCL,RC,PR是冷段上升通道预热区二回路工质与倒U型管金属壁间传热系数；dCL,MT是冷段倒U型管内径；TCL,RC，PR是冷段上升通道预热区液相工质温度；ρCL,RC，PR是冷段上升通道预热区液相工质密度；CP,CL,RC,PR是冷段上升通道预热区液相工质定压比热容；WCL,RC,PR是冷段上升通道预热区液相工质的流速；KCL,RC,BR是冷段上升通道沸腾区二回路工质与倒U型管金属壁间传热系数；TCL,RC,BR是冷段上升通道沸腾区气液混合相工质温度；ρCL,RC,BR是冷段上升通道沸腾区气液混合相工质密度；CP,CL,RC,BR是冷段上升通道沸腾区气液混合相工质定压比热容；WCL,RC,BR是冷段上升通道沸腾区气液混合相工质的流速；PCL,RC,PR是冷段上升通道预热区液相工质压力a；GCL,RC,PR是冷段上升通道预热区液相工质质量流量；fCL,RC,PR是冷段上升通道预热区摩擦因子；De,CL,RC,PR是冷段上升通道预热区当量直径；ξCL,RC,PR是冷段上升通道预热区局部阻力系数；PCL,RC,BR是冷段上升通道沸腾区气液混合相工质压力；GCL,RC,BR是冷段上升通道沸腾区气液混合相工质质量流量；fCL,RC,BR是冷段上升通道沸腾区摩擦因子；De,CL,RC,BR是冷段上升通道沸腾区当量直径；ξCL,RC,BR是冷段上升通道沸腾区局部阻力系数；
通过对冷段上升通道模型求解，得到冷段当前时刻沿倒U型管高度的二回路工质的流速、温度、压力以及传热系数分布。
8.根据权利要求1所述的蒸汽发生器二回路再循环水质量流量估计方法，其特征在于，计算当前时刻沿倒U型管高度的二回路工质质量气含率分布：
Figure FDA0003980540470000081
式中，hBR是沸腾区气液混合相工质比焓；hss是沸腾区饱和蒸汽比焓；hsw是沸腾区饱和水比焓；xBR是沸腾区工质质量气含率；hBR,hsw和hsw根据沸腾区气液混合相工质温度和压力通过工质物性参数数据库计算得到。
9.根据权利要求1所述的蒸汽发生器二回路再循环水质量流量估计方法，其特征在于，所述汽水分离器模型为式(39)～(46)：
Gss,SP,out＝(xHL,RC,BR,outGHL,RC,BR,out+xCL,RC,BR,outGCL,RC,BR,out)×η (39)
Gsw,SP,out＝(1-xHL,RC,BR,out×η)GHL,RC,BR,out+(1-xCL,RC,BR,out×η)GCL,RC,BR,out (40)
Figure FDA0003980540470000082
GSP,in＝GHL,RC,BR,out+GCL,RC,BR,out (42)
PSP,in＝PHL,RC,BR,out＝PCL,RC,BR,out (43)
TSP,in＝THL,RC,BR,out＝TCL,RC,BR,out (44)
PSP,out＝Pss,SP,out＝Tsw,SP,out (45)
TSP,in＝Tss,SP,out＝Tsw,SP,out (46)
式中，Gss,SP,out是汽水分离器出口饱和蒸汽质量流量；xHL,RC,BR,out是热段上升通道沸腾区出口气液混合相工质质量气含率；GHL,RC,BR,out是热段上升通道沸腾区出口气液混合相工质质量流量；xCL,RC,BR,out是冷段上升通道沸腾区出口气液混合相工质质量气含率；GCL,RC,BR,out是冷段上升通道沸腾区出口气液混合相工质质量流量；η是汽水分离器效率；Gsw,SP,out是汽水分离器出口饱和水质量流量；PSP,out是汽水分离器出口工质压力；Pss,SP,out是汽水分离器出口饱和蒸汽压力；Psw,SP,out是汽水分离器出口饱和水压力；PSP,in是汽水分离器入口气液混合相工质压力；TSP,in是汽水分离器入口气液混合相工质温度；TCL,RC,BR,out是冷段上升通道沸腾区出口气液混合相工质温度；THL,RC,BR,out是热段上升通道沸腾区出口气液混合相工质温度；ξSP是汽水分离器局部阻力系数；GSP,in是汽水分离器入口气液混合相工质质量流量；ρSP,in是汽水分离器入口气液混合相工质密度；PHL,RC,BR,out是热段上升通道沸腾区出口气液混合相工质压力；PCL,RC,BR,out是冷段上升通道沸腾区出口气液混合相工质压力；Tss,SP,out是汽水分离器出口饱和蒸汽温度；Tsw,SP,out是汽水分离器出口饱和水温度；ρHL,RC,BR,out是热段上升通道沸腾区出口气液混合相工质密度；ρCL,RC,BR,out是冷段上升通道沸腾区出口气液混合相工质密度；
通过对汽水分离器模型求解，得到汽水分离器出口饱和水以及饱和蒸汽的温度、压力、质量流量，其中汽水分离器出口饱和水就是再循环水。
10.一种蒸汽发生器二回路再循环水质量流量估计装置，包括存储器、处理器及存储在存储器上并能在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时用于执行上述权利要求1-9任一项所述的蒸汽发生器二回路再循环水质量流量估计方法。
CN202011308103.6A 2020-11-19 2020-11-19 蒸汽发生器二回路再循环水质量流量估计方法及系统 Active CN112417780B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011308103.6A CN112417780B (zh) 2020-11-19 2020-11-19 蒸汽发生器二回路再循环水质量流量估计方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011308103.6A CN112417780B (zh) 2020-11-19 2020-11-19 蒸汽发生器二回路再循环水质量流量估计方法及系统

Publications (2)

Publication Number Publication Date
CN112417780A CN112417780A (zh) 2021-02-26
CN112417780B true CN112417780B (zh) 2023-03-24

Family

ID=74773195

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011308103.6A Active CN112417780B (zh) 2020-11-19 2020-11-19 蒸汽发生器二回路再循环水质量流量估计方法及系统

Country Status (1)

Country Link
CN (1) CN112417780B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5307833A (en) * 1992-10-26 1994-05-03 Texaco Inc. Method and apparatus for automatically transferring and measuring wet steam between priority and secondary users
CN107784156A (zh) * 2017-09-07 2018-03-09 中国船舶重工集团公司第七〇九研究所 核动力装置蒸汽排放系统参数的计算方法
CN109799142A (zh) * 2018-12-20 2019-05-24 中广核核电运营有限公司 核电站水压试验循环加热系统及循环加热方法
CN111276266A (zh) * 2020-02-21 2020-06-12 中广核研究院有限公司 一种两环路小型堆分组运行的控制系统及其控制方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5307833A (en) * 1992-10-26 1994-05-03 Texaco Inc. Method and apparatus for automatically transferring and measuring wet steam between priority and secondary users
CN107784156A (zh) * 2017-09-07 2018-03-09 中国船舶重工集团公司第七〇九研究所 核动力装置蒸汽排放系统参数的计算方法
CN109799142A (zh) * 2018-12-20 2019-05-24 中广核核电运营有限公司 核电站水压试验循环加热系统及循环加热方法
CN111276266A (zh) * 2020-02-21 2020-06-12 中广核研究院有限公司 一种两环路小型堆分组运行的控制系统及其控制方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Extended state observer-based adaptive dynamic sliding mode control;Jiuwu Hui等;《Annals of Nuclear Energy》;20200831;全文 *
不同运行工况下蒸汽发生器热工水力稳态特性数值研究;杨元龙等;《原子能科学技术》;20130520(第05期);全文 *
核电站蒸汽发生器的建模及仿真;张腾等;《仪器仪表与分析监测》;20111128(第04期);全文 *

Also Published As

Publication number Publication date
CN112417780A (zh) 2021-02-26

Similar Documents

Publication Publication Date Title
Dumont et al. Mathematical modelling and design of an advanced once-through heat recovery steam generator
He et al. Experimental investigation on turbulent heat transfer characteristics of molten salt in a shell-and-tube heat exchanger
Colorado et al. Development and experimental validation of a computational model for a helically coiled steam generator
Omar et al. Thermal Hydraulic analysis of Syrian MNSR research reactor using RELAP5/Mod3. 2 code
Walter et al. How can the heat transfer correlations for finned-tubes influence the numerical simulation of the dynamic behavior of a heat recovery steam generator?
CN112417676B (zh) 基于核电蒸汽发生器分布参数模型的关键变量估计方法
Hong et al. Heat transfer performance test of PDHRS heat exchangers of PGSFR using STELLA-1 facility
CN112417780B (zh) 蒸汽发生器二回路再循环水质量流量估计方法及系统
CN112699523B (zh) 蒸汽发生器二回路工质汽化起始高度估计方法及装置
CN112417781B (zh) 核电蒸汽发生器出口饱和蒸汽质量流量估计方法及系统
CN112417681B (zh) 一种蒸汽发生器一二次侧对流换热系数分布的估计方法
CN112380713B (zh) 一种蒸汽发生器倒u型管金属壁温度分布估计方法
Lee et al. Implementation of a new empirical model of steam condensation for the passive containment cooling system into MARS-KS code: Application to containment transient analysis
CN112417680B (zh) 蒸汽发生器二回路工质质量气含率分布估计方法及系统
CN112417782B (zh) 一种蒸汽发生器二回路工质循环倍率估计方法
Tao et al. Numerical and experimental study on thermal-hydraulic behaviors in the secondary side of passive residual heat removal heat exchanger
Ghazanfari et al. Drift flux modeling of the VVER-1000 horizontal nuclear steam generator
JANG et al. Transient analysis of natural circulation nuclear reactor REX-10
Van De Graaf et al. Scaling laws and design aspects of a natural-circulation-cooled simulated boiling water reactor fuel assembly
Jun et al. Validation of the TASS/SMR-S code for the core heat transfer model on the steady experimental conditions
Alipchenkov et al. Progress in system thermohydraulic code HYDRA-IBRAE/LM models development for fast reactor simulation
Cherubini et al. Use of the natural circulation flow map for natural circulation systems evaluation
Wuth Comparative analysis of the Žukauskas method and data from the OSU MASLWR test facility steam generator
David et al. Heat exchanger heat transfer coefficient and CFD modelling
Albdour et al. Steam Condensation Heat Transfer in the Presence of Noncondensable Gases (NCGs) in Nuclear Power Plants (NPPs): A Comprehensive Review of Fundamentals, Current Status, and Prospects for Future Research

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