CN112417680B - 蒸汽发生器二回路工质质量气含率分布估计方法及系统 - Google Patents

蒸汽发生器二回路工质质量气含率分布估计方法及系统 Download PDF

Info

Publication number
CN112417680B
CN112417680B CN202011308067.3A CN202011308067A CN112417680B CN 112417680 B CN112417680 B CN 112417680B CN 202011308067 A CN202011308067 A CN 202011308067A CN 112417680 B CN112417680 B CN 112417680B
Authority
CN
China
Prior art keywords
working medium
channel
section
inverted
loop
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
CN202011308067.3A
Other languages
English (en)
Other versions
CN112417680A (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 CN202011308067.3A priority Critical patent/CN112417680B/zh
Publication of CN112417680A publication Critical patent/CN112417680A/zh
Application granted granted Critical
Publication of CN112417680B publication Critical patent/CN112417680B/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
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F22—STEAM GENERATION
    • F22B—METHODS OF STEAM GENERATION; STEAM BOILERS
    • F22B1/00—Methods of steam generation characterised by form of heating method
    • F22B1/02—Methods of steam generation characterised by form of heating method by exploitation of the heat content of hot heat carriers
    • F22B1/08—Methods of steam generation characterised by form of heating method by exploitation of the heat content of hot heat carriers the heat carrier being steam
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • Thermal Sciences (AREA)
  • Fluid Mechanics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Sustainable Energy (AREA)
  • Sustainable Development (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Mechanical Engineering (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Engine Equipment That Uses Special Cycles (AREA)

Abstract

本发明提供了一种蒸汽发生器二回路工质质量气含率分布估计方法及系统，获取蒸汽发生器的实时运行数据；建立下降通道模型，得到当前时刻下降通道底部出口液相工质的流量、温度及压力；计算一回路冷却剂与倒U型管金属壁之间以及倒U型管金属壁与二回路工质之间的传热系数；建立一回路冷却剂模型，得到倒U型管金属壁的温度分布；建立上升通道模型，得到沿倒U型管高度的二回路工质的流速、温度、压力分布，计算沿倒U型管高度的二回路工质质量气含率分布。本发明用于蒸汽发生器二次侧虚假水位甄别，改进水位安全限控制，为蒸汽发生器运行优化及监测提供支撑。

Description

蒸汽发生器二回路工质质量气含率分布估计方法及系统
技术领域
本发明属于核电站运行优化控制技术领域，具体地，涉及一种蒸汽发生器二回路工质质量气含率分布估计方法及系统，同时提供了一种相应的终端和存储介质。
背景技术
蒸汽发生器是核电站的枢纽设备之一，它既是一回路设备，又是二回路设备，是连接一二回路之间的纽带。在核反应堆中，核裂变产生的热量由一回路冷却剂带出，通过蒸汽发生器倒U型管将热量传递给二回路工质，使过冷状态水转变成饱和蒸汽。饱和蒸汽再流入汽轮机中做功，通过同轴发电机转换为电能。
蒸汽发生器内部结构极为复杂，内部有4474根倒U型管，以及一级旋叶式汽水分离器和二级蒸汽干燥器。在蒸汽发生器内部存在着复杂的传热过程。例如在蒸汽发生器二回路侧，倒U型管向工质的热量传递包括单相对流换热、过冷沸腾以及饱和沸腾对流换热。在沸腾对流换热过程中，二回路工质局部汽化，形成气液两相流。气泡产生、成长、脱离壁面区域的过程强烈扰动二回路水位及传热阻力。而在蒸汽发生器一回路侧，冷却剂向倒U型管的热量传递为单相对流换热。由于蒸汽发生器系统的非线性、非对称性、时滞以及两相流换热过程的复杂性，目前国内外的相关研究以对蒸汽发生器集总参数建模和稳态性能仿真为主，对蒸汽发生器内部工质动态研究较少，故研究成果并不能用于改进蒸汽发生器结构设计和运行优化，也不利于提高蒸汽发生器实时液位控制系统的控制品质。
经过对现有技术检索后发现，授权公告号为CN110020476B，授权公告日为2020年06月26日的中国发明专利《一种反应堆U型管式蒸汽发生器全三维耦合模型建立方法》，提出了一种建立反应堆U型管式蒸汽发生器全三维耦合模型方法，包括管束区域的几何模型简化、管侧及壳侧计算域的计算节点划分、管侧能量源项处理、壳侧能量源项处理、管壳一次侧和二次侧相同坐标位置通过管壁耦合处能量数据交换处理以及一次侧和二次侧各节点与周围节点之间能量数据的交换以及能量源项的相互迭代过程。该专利主要贡献在于提出了一种蒸汽发生器全三维耦合建模方法，但未提出蒸汽发生器二回路工质质量气含率分布计算方法。
综上，现有公开报道均未涉及蒸汽发生器二回路工质质量气含率分布估计问题，这一空缺有待填补。
发明内容
针对现有技术的不足，本发明的目的是提供一种蒸汽发生器二回路工质质量气含率分布估计方法及系统，同时提供了一种相应的终端和存储介质。
本发明是通过以下技术方案实现的。
根据本发明的一个方面，提供了一种蒸汽发生器二回路工质质量气含率分布估计方法，将蒸汽发生器划分为热段、冷段和汽水分离器，包括：
获取给定时刻下蒸汽发生器的实时运行数据；
利用获取的所述蒸汽发生器的实时运行数据，分别建立热段、冷段的下降通道模型，得到当前时刻下降通道底部出口液相工质的温度、压力及质量流量；
利用获取的所述蒸汽发生器的实时运行数据，计算一回路冷却剂与倒U型管金属壁之间传热系数以及倒U型管金属壁与二回路工质之间的传热系数，其中一回路冷却剂与倒U型管金属壁、预热段倒U型管金属壁与二回路工质之间的传热系数采用迪图斯-贝尔特公式计算，沸腾段倒U型管金属壁与二回路工质之间的传热系数采用Chen公式计算；
利用获取的所述蒸汽发生器的实时运行数据以及得到的所述一回路冷却剂与倒U型管金属壁之间传热系数，分别建立热段、冷段的一回路冷却剂模型，得到倒U型管金属壁的温度分布；
利用获取的所述蒸汽发生器的实时运行数据、所述倒U型管金属壁与二回路工质之间的传热系数、所述倒U型管金属壁的温度分布以及所述下降通道底部出口液相工质的温度、压力及质量流量，分别建立热段、冷段的上升通道模型，利用下降通道底部出口液相工质的温度、压力及质量流量作为输入，得到当前时刻沿倒U型管高度的二回路工质的流速、温度以及压力分布，进而计算当前时刻沿倒U型管高度的二回路工质质量气含率分布。
进一步的，所述方法还包括：利用获取的所述蒸汽发生器的实时运行数据以及上升通道顶部出口气液混合物工质的流速、温度以及压力，建立汽水分离器模型，计算得到汽水分离器出口气相工质、液相工质的温度、压力和质量流量，所述汽水分离器出口液相工质是所述下降通道模型入口再循环水。
根据本发明的另一个方面，提供了一种蒸汽发生器二回路工质质量气含率分布估计系统，包括：
数据获取模块，该模块用于获取给定时刻下蒸汽发生器的实时运行数据；
下降通道模型模块，该模块利用获取的所述蒸汽发生器的实时运行数据，建立下降通道模型，得到当前时刻下降通道底部出口液相工质的温度、压力及质量流量；
传热系数计算模块，该模块利用获取的所述蒸汽发生器的实时运行数据，计算一回路冷却剂与倒U型管金属壁之间传热系数以及倒U型管金属壁与二回路工质之间的传热系数；
一回路冷却剂模型模块，该模块利用获取的所述蒸汽发生器的相关测点实时运行数据以及得到的所述一回路冷却剂与倒U型管金属壁之间传热系数，建立一回路冷却剂模型，得到倒U型管金属壁的温度分布；
上升通道模型模块，该模块利用获取的所述蒸汽发生器的实时运行数据、所述倒U型管金属壁与二回路工质之间的传热系数、所述倒U型管金属壁的温度分布以及所述下降通道底部出口液相工质的温度、压力及质量流量，建立上升通道模型，利用下降通道底部出口液相工质的温度、压力及质量流量作为输入，得到当前时刻沿倒U型管高度的二回路工质的流速、温度以及压力分布，进而计算当前时刻沿倒U型管高度的二回路工质质量气含率分布；
汽水分离器模型模块，该模块利用获取的所述蒸汽发生器的实时运行数据以及上升通道顶部出口气液混合物工质的流速、温度以及压力，建立汽水分离器模型，计算得到汽水分离器出口气相工质、液相工质的温度、压力和质量流量。
与现有技术相比，本发明实施例具有如下至少一种有益效果：
本发明提供的蒸汽发生器二回路工质质量气含率分布估计方法及系统以及相应的终端和存储介质，实现了全工况蒸汽发生器质量气含率分布实时估计，可用于蒸汽发生器二次侧虚假水位的甄别，改进水位安全限控制，为蒸汽发生器运行优化及监测提供支撑条件，有助于提高核电站运行的安全性与经济性。
本发明提供的蒸汽发生器二回路工质质量气含率分布估计方法及系统以及相应的终端和存储介质，可描述蒸汽发生器内部二回路工质热工水力特性动态变化过程，据此可估计二回路上升通道工质平均质量气含率，并可进一步应用于虚假水位甄别和液位安全限控制策略的改进。
本发明提供的蒸汽发生器二回路工质质量气含率分布估计方法及系统以及相应的终端和存储介质，可用于蒸汽发生器二次侧虚假水位的甄别，改进水位安全限控制，为蒸汽发生器运行优化及监测提供支撑条件。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：
图1为本发明一实施例中蒸汽发生器质量气含率分布实时估计方法流程图；
图2为本发明一优选实施例中蒸汽发生器质量气含率分布实时估计方法工作示意图；
图3为本发明一优选实施例中蒸汽发生器简化结构示意图；
图4为本发明一优选实施例中核电站机组输出负荷变化图；
图5为本发明一优选实施例中热段质量气含率分布计算结果；
图6为本发明一优选实施例中冷段质量气含率分布计算结果；
图7为本发明一实施例中蒸汽发生器质量气含率分布实时估计系统组成模块示意图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进。这些都属于本发明的保护范围。
图1为本发明一实施例中蒸汽发生器二回路工质质量气含率分布估计方法流程图。
如图1所示，该实施例提供的蒸汽发生器二回路工质质量气含率分布估计方法，可以包括以下步骤：
S100，获取给定时刻下蒸汽发生器的实时运行数据；
S200，利用获取的蒸汽发生器的实时运行数据，建立下降通道模型，得到当前时刻下降通道底部出口液相工质的温度、压力及质量流量；
S300，利用获取的蒸汽发生器的实时运行数据，计算一回路冷却剂与倒U型管金属壁之间传热系数以及倒U型管金属壁与二回路工质之间的传热系数；
S400，利用获取的蒸汽发生器的实时运行数据以及得到的一回路冷却剂与倒U型管金属壁之间传热系数，建立一回路冷却剂模型，得到倒U型管金属壁的温度分布；
S500，利用获取的蒸汽发生器的实时运行数据、得到的倒U型管金属壁与二回路工质之间的传热系数、得到的倒U型管金属壁的温度分布以及得到的下降通道底部出口液相工质的温度、压力及质量流量以及，建立上升通道模型，利用得到的下降通道底部出口液相工质的温度、压力及质量流量作为输入，得到当前时刻沿倒U型管高度的二回路工质的流速、温度以及压力分布，进而计算当前时刻沿倒U型管高度的二回路工质质量气含率分布。
在该实施例的一具体实例中，给定时刻下蒸汽发生器的相关测点实时运行数据，优选包括：
-机组负荷；
-给水温度、压力以及质量流量；
-饱和蒸汽温度、压力以及质量流量；
-一回路冷却剂进出口温度、压力以及质量流量；
-水位高度。
其中，上升通道顶部出口气液混合物工质的流速、温度以及压力，通过解算上升通道模型得到。该算法可以采用现有技术实现，比如龙格库塔法等。
图2为本发明一优选实施例中蒸汽发生器二回路工质质量气含率分布估计方法工作示意图。
如图2所示，该优选实施例所提供的蒸汽发生器二回路工质质量气含率分布估计方法，可以包括以下步骤：
步骤一，获取给定时刻下蒸汽发生器的相关测点实时运行数据；
步骤二，利用获取的蒸汽发生器的相关测点实时运行数据，建立下降通道模型，得到当前时刻下降通道底部出口液相工质的温度、压力及质量流量；
步骤三，利用获取的蒸汽发生器的相关测点实时运行数据，计算一回路冷却剂与倒U型管金属壁之间传热系数以及倒U型管金属壁与二回路工质之间的传热系数；
步骤四，利用获取的蒸汽发生器的相关测点实时运行数据以及得到的一回路冷却剂与倒U型管金属壁之间传热系数，建立一回路冷却剂模型，得到倒U型管金属壁的温度分布；
步骤五，利用获取的蒸汽发生器的相关测点实时运行数据、得到的倒U型管金属壁与二回路工质之间的传热系数、得到的倒U型管金属壁的温度分布以及得到的下降通道底部出口液相工质的流量、温度、压力，建立上升通道模型，利用下降通道底部出口液相工质的温度、压力及质量流量作为输入，得到当前时刻沿倒U型管高度的二回路工质的流速、温度以及压力分布，进而计算当前时刻沿倒U型管高度的二回路工质质量气含率分布；
步骤六，利用获取的蒸汽发生器的相关测点实时运行数据以及通过对上升通道模型从下向上解算得到的上升通道顶部出口气液混合物工质的流速、温度以及压力，建立汽水分离器模型，计算得到汽水分离器出口气相工质、液相工质的温度、压力和质量流量，其中液相工质为再循环水。
作为一优选实施例，本优选实施例所提供的方法，在步骤一之前或步骤一之后，还可以包括如下步骤：
步骤零，将蒸汽发生器划分为热段、冷段和汽水分离器，其中热段和冷段又分别划分为下降通道和上升通道：下降通道是指外壳与内部套筒之间工质流经的空间，工质向下流动；上升通道是指内部套筒与倒U型管壁之间工质流经的空间，工质向上流动。
作为一优选实施例，在步骤一中，从运行机组现场分布式控制系统(DCS)实时数据库中获取给定时刻的运行数据，包括：机组负荷；给水温度、压力、质量流量；饱和蒸汽温度、压力、质量流量；一回路冷却剂进出口温度、压力、质量流量以及水位高度等。
作为一优选实施例，在步骤二中，利用获取的蒸汽发生器的相关测点实时运行数据，根据工质的质量、能量、动量守恒关系，建立下降通道模型，得到当前时刻下降通道底部出口液相工质的温度、压力及质量流量。
作为一优选实施例，在步骤零中，在蒸汽发生器的上升通道中，根据二回路工质状态将上升通道划分为预热区和沸腾区。预热区和沸腾区分界面的划分依据为：
hRC(t,z)＝hsw(t,z) (1)
式中，hRC(t,z)是上升通道当前时刻t和高度z的二回路工质的比焓，kJ/kg；hsw(t,z)是当前时刻t和高度z的二回路工质饱和状态比焓，kJ/kg。
作为一优选实施例，在步骤二中，下降通道入口液相工质中占比
Figure BDA0002788892000000061
的给水流入热段，占比
Figure BDA0002788892000000062
的给水流入冷段，占比
Figure BDA0002788892000000063
的再循环水流入热段，占比
Figure BDA0002788892000000064
的再循环水流入冷段。
作为一优选实施例，
Figure BDA0002788892000000071
的取值范围为：70～90；
Figure BDA0002788892000000072
的取值范围为：40～60。
根据动量、质量和能量守恒关系，分别建立蒸汽发生器热段下降通道模型和冷段下降通道模型；
其中，建立热段下降通道模型如式(2)～(4)所示：
Figure BDA0002788892000000073
Figure BDA0002788892000000074
Figure BDA0002788892000000075
式中，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 BDA0002788892000000076
Figure BDA0002788892000000077
Figure BDA0002788892000000081
式中，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型管金属壁之间传热系数以及倒U型管金属壁与二回路工质之间的传热系数。
热段和冷段一回路冷却剂与倒U型管金属壁之间传热系数KHL,PS和KCL,PS以及热段和冷段预热区倒U型管金属壁与二回路工质之间的传热系数KHL,RC,PR和KCL,RC,PR，统一表示为K，采用迪图斯-贝尔特公式计算：
K＝0.023Rew 0.8Prw 0.3λw/dHL,MT (8)
式中，式中，Rew是对应的热段或冷段一回路或二回路工质雷诺数；Prw是对应的热段或冷段一回路或二回路工质普朗特数；λw是对应的热段或冷段一回路或二回路工质热导率；dHL,MT是倒U型管内径；
对热段和冷段沸腾区倒U型管金属壁与二回路工质之间的传热系数K*,RC,BR采用式(9)～(14)Chen公式计算，其中，在热段的K*,RC,BR表示为KHL,RC,BR，在冷段的K*,RC,BR表示为KCL,RC,BR：
K*,RC,BR＝Kcht+Kbht (9)
Figure BDA0002788892000000082
Figure BDA0002788892000000091
Figure BDA0002788892000000092
Figure BDA0002788892000000093
Figure BDA0002788892000000094
式中，Kcht、Kbht分别是对流传热部分的传热系数和泡核沸腾传热部分的传热系数；CP,w是工质定压比热容；hfs是沸腾区液相工质汽化潜热；σ沸腾区液相工质表面张力系数；ΔTMT是沸腾区倒U型管金属壁过热度；ΔPMT是沸腾区饱和蒸汽压差；x是质量气含率；ρw是上升通道液相工质密度；ρs是上升通道饱和蒸汽密度；μw是上升通道液相工质粘性系数；μs是上升通道饱和蒸汽粘性系数；dHL,MT是热段倒U型管内径；G是工质质量流量；Xtt和S是中间变量。
作为一优选实施例，在步骤四中，利用获取的蒸汽发生器的实时运行数据，根据工质的质量、能量、动量守恒关系，建立一回路冷却剂模型，分别建立蒸汽发生器热段一回路冷却剂模型和冷段一回路冷却剂模型；
其中，
建立热段一回路冷却剂模型如式(15)～(18)所示：
Figure BDA0002788892000000095
Figure BDA0002788892000000096
Figure BDA0002788892000000097
Figure BDA0002788892000000098
式中，ρ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 BDA0002788892000000101
Figure BDA0002788892000000102
Figure BDA0002788892000000103
Figure BDA0002788892000000104
式中，ρ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 BDA0002788892000000105
Figure BDA0002788892000000106
Figure BDA0002788892000000107
Figure BDA0002788892000000108
Figure BDA0002788892000000111
Figure BDA0002788892000000112
Figure BDA0002788892000000113
式中，ρ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 BDA0002788892000000121
Figure BDA0002788892000000122
Figure BDA0002788892000000123
Figure BDA0002788892000000124
式中，ρ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 BDA0002788892000000131
式中，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 BDA0002788892000000132
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；
通过对汽水分离器模型求解，得到汽水分离器出口饱和水以及饱和蒸汽的温度、压力、质量流量，其中液相工质为再循环水。
在上述各实施例基础上，还可以进一步包括S600，该步骤利用获取的蒸汽发生器的实时运行数据以及上升通道顶部出口气液混合物工质的流速、温度以及压力，建立汽水分离器模型，计算得到汽水分离器出口气相工质、液相工质的温度、压力和质量流量。其中，液相工质即为再循环水，即S200中下降通道模型入口再循环水。
该优选实施例提供的蒸汽发生器质量气含率分布实时估计方法，利用从运行机组现场DCS分布式控制系统实时测量数据库中获取给定时刻下机组负荷，给水温度、压力、质量流量，饱和蒸汽温度、压力、质量流量，一回路冷却剂进出口温度、压力、质量流量以及水位高度等数据，结合工质物性参数数据库和蒸汽发生器结构参数库，解算热段模型、冷段模型和汽水分离器模型，输出热段和冷段质量气含率分布。
如图3所示，为蒸汽发生器简化结构示意图。根据蒸汽发生器真实结构，进行简化划分为热段、冷段和汽水分离器。热段和冷段的二回路又可以根据工质是否达到饱和状态划分为预热区和沸腾区。
该优选实施例提供的蒸汽发生器二回路工质质量气含率分布估计方法，如图4所示，获取实施例核电站机组2019年6月10日不同负荷下的蒸汽发生器DCS实测数据。
图5是热段二回路工质沿z方向的质量气含率分布估计结果，图6是冷段二回路工质沿z方向的质量气含率分布估计结果。由图5和图6可知，在沸腾区任一高度z的质量气含率随着核电站机组负荷升高而升高，反之亦反。
在本发明另一实施例中，还提供了一种蒸汽发生器二回路工质质量气含率分布估计系统，如图7所示，可以包括：数据获取模块、下降通道模型模块、传热系数计算模块、一回路冷却剂模型模块、上升通道模型模块以及汽水分离器模型模块。其中：数据获取模块，该模块用于获取给定时刻下蒸汽发生器的相关测点实时运行数据；下降通道模型模块，该模块利用获取的蒸汽发生器的相关测点实时运行数据，建立下降通道模型，得到当前时刻下降通道底部出口液相工质的温度、压力及质量流量；传热系数计算模块，该模块利用获取的蒸汽发生器的相关测点实时运行数据，计算一回路冷却剂与倒U型管金属壁之间传热系数以及倒U型管金属壁与二回路工质之间的传热系数；一回路冷却剂模型模块，该模块利用获取的蒸汽发生器的相关测点实时运行数据以及得到的一回路冷却剂与倒U型管金属壁之间传热系数，建立一回路冷却剂模型，得到倒U型管金属壁的温度分布；上升通道模型模块，该模块利用获取的蒸汽发生器的相关测点实时运行数据、得到的倒U型管金属壁与二回路工质之间的传热系数、得到的倒U型管金属壁的温度分布以及得到的下降通道底部出口液相工质的温度、压力及质量流量以及，建立上升通道模型，利用下降通道底部出口液相工质的温度、压力及质量流量作为输入，得到当前时刻沿倒U型管高度的二回路工质的流速、温度以及压力分布，进而计算当前时刻沿倒U型管高度的二回路工质质量气含率分布。当然，可选地，还可以包括汽水分离器模型模块，该模块利用获取的蒸汽发生器的相关测点实时运行数据以及上升通道顶部出口气液混合物工质的流速、温度以及压力，建立汽水分离器模型，计算得到汽水分离器出口气相工质、液相工质的温度、压力和质量流量。
本发明上述实施例提供了一种蒸汽发生器二回路工质质量气含率分布估计方法，给出了一种机理建模与DCS(运行机组现场分布式控制系统)实时测量数据相结合的蒸汽发生器二回路工质质量气含率分布估计(辨识)方法，该方法以蒸汽发生器为研究对象，根据蒸汽发生器具体结构，将蒸汽发生器划分为热段、冷段和汽水分离器。基于质量、能量、动量动态衡算建立蒸汽发生器机理模型，完成模型验证，最后基于模型和DCS测量数据实时计算二回路工质质量气含率分布，质量气含率定义为单位时间内流过某一截面的两相流体总质量中气相所占的比例份额。
本发明上述实施例提供的蒸汽发生器二回路工质质量气含率分布估计方法及系统以及相应的终端和存储介质，其中，工质质量气含率指单位时间内，流过某一截面的水-水蒸气两相流流体总质量中气相所占的比例份额。本发明上述实施例获取给定时刻下蒸汽发生器相关测点的实时测量数据；将蒸汽发生器划分为热段、冷段和汽水分离器，其中热段和冷段又划分为下降通道和上升通道：下降通道是指蒸汽发生器外壳与内部套筒之间二回路工质流经的空间，上升通道是指内部套筒与倒U型管管壁之间二回路工质流经的空间；求解下降通道模型得到当前时刻下降通道底部出口二回路工质的比焓、流量、温度、压力以及密度；求解上升通道模型得到当前时刻沿倒U型管高度的二回路工质的比焓、流量、温度、压力以及密度分布，进而计算出当前时刻沿倒U型管高度的上升通道二回路工质质量气含率分布；求解汽水分离器模型计算得到汽水分离器出口工质温度、压力和流量；可描述蒸汽发生器内部二回路工质热工水力特性动态变化过程，据此可估计二回路上升通道工质平均质量气含率，并可进一步应用于虚假水位甄别和液位安全限控制策略的改进。本发明上述实施例提供的技术方案，实现了全工况蒸汽发生器质量气含率分布实时估计，可用于蒸汽发生器二次侧虚假水位的甄别，改进水位安全限控制，为蒸汽发生器运行优化及监测提供支撑条件，有助于提高核电站运行的安全性与经济性。
需要说明的是，本发明提供的方法中的步骤，可以利用系统中对应的模块、装置、单元等予以实现，本领域技术人员可以参照方法的技术方案实现系统的组成，即，方法中的实施例可理解为构建系统的优选例，在此不予赘述。
本领域技术人员知道，除了以纯计算机可读程序代码方式实现本发明提供的系统及其各个装置以外，完全可以通过将方法步骤进行逻辑编程来使得本发明提供的系统及其各个装置以逻辑门、开关、专用集成电路、可编程逻辑控制器以及嵌入式微控制器等的形式来实现相同功能。所以，本发明提供的系统及其各项装置可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构；也可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变形或修改，这并不影响本发明的实质内容。

Claims (11)

1.一种蒸汽发生器二回路工质质量气含率分布估计方法，将蒸汽发生器划分为热段、冷段和汽水分离器，其特征在于，包括：
获取给定时刻下蒸汽发生器的实时运行数据；
利用获取的所述蒸汽发生器的实时运行数据，分别建立热段、冷段的下降通道模型，得到当前时刻下降通道底部出口液相工质的温度、压力及质量流量；
利用获取的所述蒸汽发生器的实时运行数据，计算一回路冷却剂与倒U型管金属壁之间传热系数以及倒U型管金属壁与二回路工质之间的传热系数，其中一回路冷却剂与倒U型管金属壁和预热段倒U型管金属壁与二回路工质之间的传热系数采用迪图斯-贝尔特公式计算，沸腾段倒U型管金属壁与二回路工质之间的传热系数采用Chen公式计算；其中，对热段和冷段沸腾区的所述倒U型管金属壁与二回路工质之间的传热系数KHL,RC,BR和KCL,RC,BR，统一表示为K*,RC,BR，所述Chen公式采用式(9)～(14)计算如下：
K*,RC,BR＝Kcht+Kbht (9)
Figure FDA0003842039250000011
Figure FDA0003842039250000012
Figure FDA0003842039250000013
Figure FDA0003842039250000014
Figure FDA0003842039250000015
式中，Kcht、Kbht分别是对流传热部分的传热系数和泡核沸腾传热部分的传热系数；CP,w是工质定压比热容；hfs是沸腾区液相工质汽化潜热；σ沸腾区液相工质表面张力系数；ΔTMT是沸腾区倒U型管金属壁过热度；ΔPMT是沸腾区饱和蒸汽压差；x是质量气含率；ρw是上升通道液相工质密度；ρs是上升通道饱和蒸汽密度；μw是上升通道液相工质粘性系数；μs是上升通道饱和蒸汽粘性系数；dHL,MT是热段倒U型管内径；G是工质质量流量；Xtt和S是中间变量；
利用获取的所述蒸汽发生器的实时运行数据以及得到的所述一回路冷却剂与倒U型管金属壁之间传热系数，分别建立热段、冷段的一回路冷却剂模型，得到倒U型管金属壁的温度分布；
利用获取的所述蒸汽发生器的实时运行数据、所述倒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所述的蒸汽发生器二回路工质质量气含率分布估计方法，其特征在于，根据蒸汽发生器的实时运行数据中蒸汽发生器下降通道入口液相工质的动量、质量和能量守恒关系，分别建立蒸汽发生器热段下降通道模型和冷段下降通道模型；
其中：
建立的热段下降通道模型如式(2)～(4)所示：
Figure FDA0003842039250000031
Figure FDA0003842039250000032
Figure FDA0003842039250000033
式中，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 FDA0003842039250000034
Figure FDA0003842039250000035
Figure FDA0003842039250000041
式中，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.根据权利要求4所述的蒸汽发生器二回路工质质量气含率分布估计方法，其特征在于，所述蒸汽发生器下降通道入口液相工质中，占比
Figure FDA0003842039250000042
的给水流入热段，占比
Figure FDA0003842039250000043
的给水流入冷段，占比
Figure FDA0003842039250000044
的再循环水流入热段，占比
Figure FDA0003842039250000045
的再循环水流入冷段。
6.根据权利要求1所述的蒸汽发生器二回路工质质量气含率分布估计方法，其特征在于，热段和冷段的所述一回路冷却剂与倒U型管金属壁之间传热系数KHL,PS和KCL,PS以及热段和冷段预热区的所述倒U型管金属壁与二回路工质之间的传热系数KHL,RC,PR和KCL,RC,PR，统一表示为K，采用迪图斯-贝尔特公式计算如下：
K＝0.023Rew 0.8Prw 0.3λw/dHL,MT (8)
式中，Rew是对应的热段或冷段一回路或二回路工质雷诺数；Prw是对应的热段或冷段一回路或二回路工质普朗特数；λw是对应的热段或冷段一回路或二回路工质热导率；dHL,MT是倒U型管内径。
7.根据权利要求1所述的蒸汽发生器二回路工质质量气含率分布估计方法，其特征在于，根据蒸汽发生器的实时运行数据以及一回路冷却剂的动量、质量和能量守恒关系，分别建立蒸汽发生器热段一回路冷却剂模型和冷段一回路冷却剂模型；
其中，
建立的热段一回路冷却剂模型如式(15)～(18)所示：
Figure FDA0003842039250000046
Figure FDA0003842039250000051
Figure FDA0003842039250000052
Figure FDA0003842039250000053
式中，ρHL,PS是热段一回路冷却剂密度；WHL,PS是热段一回路冷却剂的流速；CP,HL,PS是热段一回路冷却剂定压比热容；THL,PS是热段一回路冷却剂温度；KHL,PS是热段一回路冷却剂通过倒U型管金属壁向二回路工质传热的传热系数；dHL,MT是热段倒U型管内径；PHL,PS是热段一回路冷却剂压力；
通过对热段一回路冷却剂模型求解，得到热段倒U型管金属壁的温度分布；
建立的冷段一回路冷却剂模型如式(19)～(22)所示：
Figure FDA0003842039250000054
Figure FDA0003842039250000055
Figure FDA0003842039250000056
Figure FDA0003842039250000057
式中，ρCL,PS是冷段一回路冷却剂密度；WCL,PS是冷段一回路冷却剂的流速；CP,CL,PS是冷段一回路冷却剂定压比热容；TCL,PS是冷段一回路冷却剂温度；KCL,PS是冷段一回路冷却剂通过倒U型管金属壁向二回路工质传热的传热系数；dCL,MT是冷段倒U型管内径；PCL,PS是冷段一回路冷却剂压力；
通过对冷段一回路冷却剂模型求解，得到冷段倒U型管金属壁的温度分布。
8.根据权利要求1所述的蒸汽发生器二回路工质质量气含率分布估计方法，其特征在于，根据蒸汽发生器的实时运行数据以及蒸汽发生器上升通道工质的动量、质量和能量守恒关系，分别建立蒸汽发生器热段上升通道模型和冷段上升通道模型；
其中：
建立的热段上升通道模型如式(23)～(30)所示：
Figure FDA0003842039250000058
Figure FDA0003842039250000059
Figure FDA0003842039250000061
Figure FDA0003842039250000062
Figure FDA0003842039250000063
Figure FDA0003842039250000064
Figure FDA0003842039250000065
Figure FDA0003842039250000066
式中，ρ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 FDA0003842039250000071
Figure FDA0003842039250000072
Figure FDA0003842039250000073
Figure FDA0003842039250000074
Figure FDA0003842039250000075
Figure FDA0003842039250000076
Figure FDA0003842039250000077
式中，ρ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型管高度的二回路工质的流速、温度、压力以及传热系数分布。
9.根据权利要求1所述的蒸汽发生器二回路工质质量气含率分布估计方法，其特征在于，所述计算当前时刻沿倒U型管高度的二回路工质质量气含率分布，包括：
Figure FDA0003842039250000081
式中，hBR是沸腾区气液混合相工质比焓；hss是沸腾区饱和蒸汽比焓；hsw是沸腾区饱和水比焓；xBR是沸腾区工质质量气含率；hBR,hsw和hsw根据沸腾区气液混合相工质温度和压力通过工质物性参数数据库计算得到。
10.根据权利要求1所述的蒸汽发生器二回路工质质量气含率分布估计方法，其特征在于，还包括：利用获取的所述蒸汽发生器的实时运行数据以及上升通道顶部出口气液混合物工质的流速、温度以及压力，建立汽水分离器模型，计算得到汽水分离器出口气相工质、液相工质的温度、压力和质量流量，所述汽水分离器出口液相工质是所述下降通道模型入口再循环水。
11.一种蒸汽发生器二回路工质质量气含率分布估计系统，其特征在于，包括：
数据获取模块，该模块用于获取给定时刻下蒸汽发生器的实时运行数据；
下降通道模型模块，该模块利用获取的所述蒸汽发生器的实时运行数据，建立下降通道模型，得到当前时刻下降通道底部出口液相工质的温度、压力及质量流量；
传热系数计算模块，该模块利用获取的所述蒸汽发生器的实时运行数据，计算一回路冷却剂与倒U型管金属壁之间传热系数以及倒U型管金属壁与二回路工质之间的传热系数；
一回路冷却剂模型模块，该模块利用获取的所述蒸汽发生器的相关测点实时运行数据以及得到的所述一回路冷却剂与倒U型管金属壁之间传热系数，建立一回路冷却剂模型，得到倒U型管金属壁的温度分布；
上升通道模型模块，该模块利用获取的所述蒸汽发生器的实时运行数据、所述倒U型管金属壁与二回路工质之间的传热系数、所述倒U型管金属壁的温度分布以及所述下降通道底部出口液相工质的温度、压力及质量流量，建立上升通道模型，利用下降通道底部出口液相工质的温度、压力及质量流量作为输入，得到当前时刻沿倒U型管高度的二回路工质的流速、温度以及压力分布，进而计算当前时刻沿倒U型管高度的二回路工质质量气含率分布；
汽水分离器模型模块，该模块利用获取的所述蒸汽发生器的实时运行数据以及上升通道顶部出口气液混合物工质的流速、温度以及压力，建立汽水分离器模型，计算得到汽水分离器出口气相工质、液相工质的温度、压力和质量流量。
CN202011308067.3A 2020-11-19 2020-11-19 蒸汽发生器二回路工质质量气含率分布估计方法及系统 Active CN112417680B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011308067.3A CN112417680B (zh) 2020-11-19 2020-11-19 蒸汽发生器二回路工质质量气含率分布估计方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011308067.3A CN112417680B (zh) 2020-11-19 2020-11-19 蒸汽发生器二回路工质质量气含率分布估计方法及系统

Publications (2)

Publication Number Publication Date
CN112417680A CN112417680A (zh) 2021-02-26
CN112417680B true CN112417680B (zh) 2022-11-22

Family

ID=74773250

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011308067.3A Active CN112417680B (zh) 2020-11-19 2020-11-19 蒸汽发生器二回路工质质量气含率分布估计方法及系统

Country Status (1)

Country Link
CN (1) CN112417680B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103616136A (zh) * 2013-11-21 2014-03-05 华北电力大学（保定） 一种热力管道阀门泄漏状态监测方法
CN103776502A (zh) * 2014-01-20 2014-05-07 上海交通大学 火电机组中低压缸入口热再热蒸汽质量流量实时计量方法
CN107289436A (zh) * 2016-03-31 2017-10-24 华北电力大学 一种具有提高蒸汽发生器换热效率的支撑导流结构

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103616136A (zh) * 2013-11-21 2014-03-05 华北电力大学（保定） 一种热力管道阀门泄漏状态监测方法
CN103776502A (zh) * 2014-01-20 2014-05-07 上海交通大学 火电机组中低压缸入口热再热蒸汽质量流量实时计量方法
CN107289436A (zh) * 2016-03-31 2017-10-24 华北电力大学 一种具有提高蒸汽发生器换热效率的支撑导流结构

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"On switching H1 controllers for nuclear steam generator water level:A multiple parameter-dependent Lyapunov functions approach";Ke Hu etc.;《Annals of Nuclear Energy》;20081231;全文 *
"蒸汽发生器二回路过冷沸腾的数值模拟";张小英，等;《华南理工大学学报（自然科学版）》;20130531;第41卷(第5期);全文 *

Also Published As

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

Similar Documents

Publication Publication Date Title
CN104834773B (zh) 一种直管式直流蒸汽发生器换热性能的仿真方法
Qi et al. Numerical simulation of the heat transfer of superheater tubes in power plants considering oxide scale
Prusek et al. Deposit models for tube support plate flow blockage in Steam Generators
Zhu et al. Dynamic modeling, validation and analysis of direct air-cooling condenser with integration to the coal-fired power plant for flexible operation
Norouzi et al. Heat recovery steam generator: Constructal thermoeconomic optimization
CN110059388B (zh) 核电厂冷凝器沿工质流动方向的参数分布建模仿真方法
Ding et al. Development of a model for thermal-hydraulic analysis of helically coiled tube once-through steam generator based on Modelica
CN112417680B (zh) 蒸汽发生器二回路工质质量气含率分布估计方法及系统
CN112417676B (zh) 基于核电蒸汽发生器分布参数模型的关键变量估计方法
Deghal Cheridi et al. Modeling and simulation of a natural circulation water-tube steam boiler
Ferruzza et al. A method to account for transient performance requirements in the design of steam generators for concentrated solar power applications
CN112417781B (zh) 核电蒸汽发生器出口饱和蒸汽质量流量估计方法及系统
CN112417780B (zh) 蒸汽发生器二回路再循环水质量流量估计方法及系统
CN112699523B (zh) 蒸汽发生器二回路工质汽化起始高度估计方法及装置
CN112417681B (zh) 一种蒸汽发生器一二次侧对流换热系数分布的估计方法
CN112380713B (zh) 一种蒸汽发生器倒u型管金属壁温度分布估计方法
CN112417782B (zh) 一种蒸汽发生器二回路工质循环倍率估计方法
CN112446097A (zh) 一种蒸汽发生器体积和负荷的多目标优化方法
Van De Graaf et al. Scaling laws and design aspects of a natural-circulation-cooled simulated boiling water reactor fuel assembly
Kumar et al. Optimization of heat transfer coefficient during condensation of refrigerant inside plain horizontal tube using teaching-learning based optimization algorithm
Xu et al. Validation of RELAP5/MOD3. 4 for Flashing-Induced Instabilities in a Natural Circulation Loop
Zhang et al. Experimental study on the pressure drop characteristics of supercritical CO
Yao et al. Refined Modeling and Simulation of Condensing Equipment in Nuclear Power Plant
Li et al. Investigation on thermohydraulic characteristics of novel superheated steam supply system for low temperature nuclear heating reactor
Jiang et al. Thermal hydraulic modeling of a natural circulation loop

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