CN112699523B - 蒸汽发生器二回路工质汽化起始高度估计方法及装置 - Google Patents
蒸汽发生器二回路工质汽化起始高度估计方法及装置 Download PDFInfo
- Publication number
- CN112699523B CN112699523B CN202011303312.1A CN202011303312A CN112699523B CN 112699523 B CN112699523 B CN 112699523B CN 202011303312 A CN202011303312 A CN 202011303312A CN 112699523 B CN112699523 B CN 112699523B
- 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
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
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Engine Equipment That Uses Special Cycles (AREA)
Abstract
本发明提供了一种蒸汽发生器二回路工质汽化起始高度估计方法及装置,其中:获取给定时刻下蒸汽发生器的实时运行数据;计算一回路冷却剂与倒U型管金属壁之间的传热系数以及倒U型管金属壁与二回路工质之间的传热系数;建立下降通道模型,得到当前时刻下降通道底部出口液相工质的流量、温度及压力;建立一回路冷却剂模型,得到倒U型管金属壁的温度分布;建立上升通道模型,得到当前时刻沿倒U型管高度的二回路工质的流速、温度、压力以及汽化起始高度;建立汽水分离器模型,计算得到汽水分离器出口气相工质和液相工质的温度、压力和质量流量。本发明开发了蒸汽发生器二回路不可测量的工质汽化起始高度的实时估计方法。
Description
技术领域
本发明涉及核电站运行优化控制技术领域,具体地,涉及一种蒸汽发生器二回路工质汽化起始高度估计方法及装置。
背景技术
立式倒U型自然循环蒸汽发生器是核电站的主要设备之一。由于蒸汽发生器结构极其复杂,一回路冷却剂在4474根倒U型管内流动,通过复杂的传热过程,将热量传递给二回路过冷水,使其逐渐气化变成饱和蒸汽,输入到汽轮机中带动发电机做功。在蒸汽发生器二回路侧,倒U型管向工质的热量传递包括单相对流换热、过冷沸腾以及饱和沸腾对流换热。在沸腾对流换热过程中,二回路工质局部汽化,形成气液两相流。气泡产生、成长和脱离壁面区域的过程强烈扰动二回路水位及传热阻力。而在蒸汽发生器一回路侧,冷却剂向倒U型管的热量传递为单相对流换热。由于蒸汽发生器系统的非线性、非对称性、时滞以及两相流换热过程的复杂性,目前国内外的相关研究以对蒸汽发生器集总参数建模和稳态性能仿真为主,对蒸汽发生器内部工质动态研究较少,故研究成果并不能用于改进蒸汽发生器结构设计和运行优化,也不利于提高蒸汽发生器实时液位控制系统的控制品质。
经过对现有技术检索后发现,授权公告号为CN201910173085.6,授权公告日为2019 年06月11日的中国发明专利《一种用于蒸汽发生器的内部结构的温度场建模方法》,提出了一种建立反应堆U型管式蒸汽发生器内部结构的温度场建模方法,根据蒸汽发生器的内部结构建立初始有限元模型,并通过优化得到最佳的边界条件。该专利主要贡献在于研究了蒸汽发生器内部工质的温度场分布,但未提及蒸汽发生器二回路工质汽化起始高度计算方法。
综上,现有公开报道均未涉及核电蒸汽发生器二回路工质汽化起始高度实时估计问题,这一空缺有待填补。
发明内容
针对现有技术的不足,本发明的目的是提供一种核电蒸汽发生器二回路工质汽化起始高度实时估计方法和装置,是一种机理建模与DCS(运行机组现场分布式控制系统)实时测量数据相结合的方法及装置。
本发明的第一方面,提供了一种蒸汽发生器二回路工质汽化起始高度估计方法,将蒸汽发生器划分为热段、冷段和汽水分离器,包括:
S1,获取给定时刻下蒸汽发生器的实时运行数据;
S2,利用S1得到的所述实时运行数据分别建立热段、冷段的下降通道模型,通过所述下降通道模型得到当前时刻下降通道底部出口液相工质的温度、压力及质量流量;
S3,利用S1得到的所述实时运行数据计算一回路冷却剂与倒U型管金属壁之间传热系数以及倒U型管金属壁与二回路工质之间的传热系数;
S4,利用S1得到的所述实时运行数据以及S3得到的所述一回路冷却剂与倒U 型管金属壁之间传热系数,分别建立热段、冷段的一回路冷却剂模型,计算得到倒 U型管金属壁的温度分布;
S5,利用S1得到的所述实时运行数据、S3得到的所述倒U型管金属壁与二回路工质之间的传热系数、S4得到的所述倒U型管金属壁的温度分布以及S2得到的所述当前时刻下降通道底部出口液相工质的温度、压力及质量流量,分别建立热段、冷段的上升通道模型,通过所述上升通道模型计算得到当前时刻沿倒U型管高度的二回路工质的流速、温度、压力以及汽化起始高度;其中所述当前时刻沿倒U型管高度的二回路工质的流速、温度、压力包括当前时刻上升通道顶部出口气液混合物工质的流速、温度以及压力;
S6,利用S1得到的所述实时运行数据以及S5得到的所述上升通道顶部出口气液混合物工质的流速、温度以及压力,建立汽水分离器模型,通过所述汽水分离器模型计算得到汽水分离器出口液相工质的温度、压力和质量流量,所述汽水分离器出口液相工质是所述下降通道模型入口再循环水;
S7,利用S6得到的所述汽水分离器出口液相工质的温度、压力和质量流量,重新执行步骤S2~S6若干次,得到最终的汽化起始高度估计值。
本发明的第二方面,提供了一种蒸汽发生器二回路工质汽化起始高度估计装置,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时用于执行上述的蒸汽发生器二回路工质汽化起始高度估计方法。
与现有技术相比,本发明实施例具有以下至少一种有益效果:
本发明实施例提供的蒸汽发生器二回路工质汽化起始高度估计方法,通过获取蒸汽发生器的实时运行数据;建立下降通道模型、一回路冷却剂模型、上升通道模型、汽水分离器模型,最终可以得出蒸汽发生器二回路工质汽化起始高度估计值,该值在工艺上无法通过仪表实际测量,解决了现有蒸汽发生器二回路的工质汽化起始高度无法实时测量的难题。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为本发明一优选实施例蒸汽发生器二回路工质汽化起始高度实时估计方法示意图;
图2为本发明一优选实施例蒸汽发生器简化结构示意图;
图3为本发明一优选实施例中核电站机组输出负荷变化图;
图4为本发明一优选实施例中蒸汽发生器二回路工质汽化起始高度估计结果。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进。这些都属于本发明的保护范围。
本发明实施例提供了一种核电蒸汽发生器二回路工质汽化起始高度实时估计方法,以蒸汽发生器为研究对象,根据蒸汽发生器具体结构,将蒸汽发生器划分为热段、冷段和汽水分离器。此部分均为现有结构,不做详述。本发明实施例基于质量、能量、动量动态衡算建立蒸汽发生器机理模型,完成模型验证,最后基于模型和DCS 测量数据实时计算二回路工质汽化起始高度。
具体的,参照图1所示,本实施例中的核电蒸汽发生器二回路工质汽化起始高度实时估计方法,包括以下步骤:
步骤一,获取给定时刻下蒸汽发生器的实时运行数据;
本步骤中,从运行机组现场分布式控制系统(DCS)实时数据库中获取给定时刻的运行数据,包括:机组负荷;给水温度、压力、质量流量;饱和蒸汽温度、压力、质量流量;一回路冷却剂进出口温度、压力、质量流量以及水位高度等。
步骤二,利用获取的蒸汽发生器的实时运行数据,建立下降通道模型,得到当前时刻下降通道底部出口液相工质的温度、压力及质量流量;具体的,可以根据工质的质量、能量、动量守恒关系,建立下降通道模型,得到当前时刻下降通道底部出口液相工质的温度、压力及质量流量。
步骤三,利用获取的蒸汽发生器的实时运行数据,计算一回路冷却剂与倒U型管金属壁之间传热系数以及倒U型管金属壁与二回路工质之间的传热系数;
步骤四,利用获取的蒸汽发生器的实时运行数据,根据工质的质量、能量、动量守恒关系,建立一回路冷却剂模型,得到倒U型管金属壁的温度分布;
步骤五,利用获取的蒸汽发生器的实时运行数据以及得到的下降通道底部出口液相工质的流量、温度、压力,结合工质的质量、能量、动量守恒关系,建立上升通道模型,得到当前时刻沿倒U型管高度的二回路工质的流速、温度、压力以及汽化起始高度;此处得到的参数影响步骤六,步骤六结果影响步骤二,最终影响到汽化起始高度。
步骤六,利用获取的蒸汽发生器的实时运行数据以及得到的上升通道顶部出口气液混合物工质的流速、温度以及压力,建立汽水分离器模型,计算得到汽水分离器出口液相工质的温度、压力和质量流量,汽水分离器出口液相工质即为下降通道模型入口再循环水;
步骤七,利用算得的汽水分离器出口液相工质参数重新计算步骤二至六若干次,得到最终的汽化起始高度估计值。
作为一优选实施例,参照图2所示,一般蒸汽发生器划分为热段、冷段和汽水分离器,其中热段和冷段又分别划分为下降通道和上升通道:下降通道是指外壳与内部套筒之间工质流经的空间,工质向下流动;上升通道是指内部套筒与倒U型管壁之间工质流经的空间,工质向上流动。在蒸汽发生器上升通道模型中,根据二回路工质状态将上升通道划分为预热区和沸腾区。预热区和沸腾区分界面的划分依据为:
hRC(t,z)=hsw(t,z) (1)
式中,hRC(t,z)是上升通道当前时刻t和高度z的二回路工质的比焓,kJ/kg;hsw(t,z) 是当前时刻t和高度z的二回路工质饱和状态比焓,kJ/kg。当在z0高度时,hRC(t,z0) ≥hsw(t,z0),则z0即是汽化起始高度。
作为一优选实施例,在步骤二中,下降通道入口液相工质中占比的给水流入热段,占比的给水流入冷段,占比的再循环水流入热段,占比的再循环水流入冷段。的取值范围可以为70~90;的取值范围可以为40~60。根据动量、质量和能量守恒关系,建立热段下降通道模型为式(2)~(4):
式中,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):
式中,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,上述4个参数简记为K,均分别采用迪图斯- 贝尔特公式计算:
K=0.023Rew 0.8Prw 0.3λw/dHL,MT (8)
式中,Rew是对应的热段或冷段一回路或二回路工质雷诺数;Prw是对应的热段或冷段一回路或二回路工质普朗特数;λw是对应的热段或冷段一回路或二回路工质热导率,kW/(m·K);dHL,MT是倒U型管内径,m;
将热段和冷段沸腾区倒U型管金属壁与二回路工质之间的传热系数KHL,RC,BR和KCL,RC,BR简记为K*,两参数均分别采用式(9)~(14)计算:
K=Kcht+Kbht (9)
式中,Kcht、Kbht分别是对流传热部分的传热系数和泡核沸腾传热部分的传热系数,kW/(m2·K);CP,w是工质定压比热容,kJ/(kg·K);hfs是沸腾区液相工质汽化潜热,kJ/kg;σ沸腾区液相工质表面张力系数;△TMT是沸腾区倒U型管金属壁过热度,K;△PMT是沸腾区饱和蒸汽压差,MPa;x是质量气含率;ρw是上升通道液相工质密度;ρs是上升通道饱和蒸汽密度;μw是上升通道液相工质粘性系数;μs是上升通道饱和蒸汽粘性系数;dHL,MT是热段倒U型管内径;G是工质质量流量;Xtt和S 是中间变量。
作为一优选实施例,在步骤四中,考虑热段一回路冷却剂重力压降,根据动量、质量和能量守恒关系,建立热段一回路冷却剂模型为式(15)~(18):
式中,ρ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):
式中,ρ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):
式中,ρ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):
式中,ρ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型管高度的二回路工质质量气含率分布:
式中,hBR是沸腾区气液混合相工质比焓,kJ/kg;hss是沸腾区饱和蒸汽比焓, kJ/kg;hsw是沸腾区饱和水比焓,kJ/kg;xBR是沸腾区工质质量气含率,%。hBR,hss和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)
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)是在当前时刻下高度为z的二回路工质饱和状态的比焓,kJ/kg。若在z0高度时,hRC(t,z0)≥hsw(t,z0),则z0即是汽化起始高度。
步骤三、下降通道入口液相工质中占比80%的给水流入热段,占比20%的给水流入冷段,占比50%的再循环水流入热段,占比50%的再循环水流入冷段。求解热段下降通道模型,式(2)~(4):
式中,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):
式中,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,上述4个参数简记为K,均分别采用迪图斯-贝尔特公式计算:
K=0.023Rew 0.8Prw 0.3λw/dHL,MT (8)
式中,Rew是工质雷诺数;Prw是工质普朗特数;λw是工质热导率,kW/(m·K); dHL,MT是倒U型管内径,m;
将热段和冷段沸腾区倒U型管金属壁与二回路工质之间的传热系数KHL,RC,BR和KCL,RC,BR简记为K*,两参数均分别采用式(9)~(14)计算:
K=Kcht+Kbht (9)
式中,Kcht、Kbht分别是对流传热部分的传热系数和泡核沸腾传热部分的传热系数,kW/(m2·K);CP,w是工质定压比热容,kJ/(kg·K);hfs是沸腾区液相工质汽化潜热,kJ/kg;σ沸腾区液相工质表面张力系数;△TMT是沸腾区倒U型管金属壁过热度,K;△PMT是沸腾区饱和蒸汽压差,MPa;x是质量气含率;ρw是上升通道液相工质密度;ρs是上升通道饱和蒸汽密度;μw是上升通道液相工质粘性系数;μs是上升通道饱和蒸汽粘性系数;dHL,MT是热段倒U型管内径;G是工质质量流量;Xtt和S 是中间变量。
步骤五、考虑热段一回路冷却剂重力压降,根据动量、质量和能量守恒关系,求解热段一回路冷却剂模型,式(15)~(18):
式中,ρ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):
式中,ρ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):
式中,ρ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):
式中,ρ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型管高度的二回路工质质量气含率分布:
式中,hBR是沸腾区气液混合相工质比焓,kJ/kg;hss是沸腾区饱和蒸汽比焓, kJ/kg;hsw是沸腾区饱和水比焓,kJ/kg;xBR是沸腾区工质质量气含率,%。hBR,hss和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)
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所示,为本仿真实例获取核电站机组不同负荷下的蒸汽发生器DCS实测数据。图4是蒸汽发生器二回路工质汽化起始高度估计结果。由图4可知,在沸腾区工质汽化起始高度随着核电站机组负荷降低而升高,反之亦反。
本发明上述实施例提供的核电蒸汽发生器二回路工质汽化起始高度实时估计方法及系统,其中,再循环水指流出汽水分离器的饱和水。本发明上述实施例获取给定时刻下蒸汽发生器相关测点的实时测量数据;将蒸汽发生器划分为热段、冷段和汽水分离器,其中热段和冷段又划分为下降通道和上升通道:下降通道是指蒸汽发生器外壳与内部套筒之间二回路工质流经的空间,上升通道是指内部套筒与倒U 型管管壁之间二回路工质流经的空间;求解下降通道模型得到当前时刻下降通道底部出口二回路工质的比焓、流量、温度、压力以及密度;求解上升通道模型得到当前时刻沿倒U型管高度的二回路工质的、流量、温度、压力以及汽化起始高度;求解汽水分离器模型计算得到汽水分离器出口工质温度、压力和流量。本发明上述实施例提供的技术方案,实现了全工况蒸汽发生器二回路工质汽化起始高度实时估计,可用于了解二回路工质汽化过程,为蒸汽发生器运行优化及监测提供支撑条件,有助于提高核电站运行的安全性与经济性。
需要说明的是,本发明提供的方法中的步骤,可以利用系统中对应的模块、装置、单元等予以实现,本领域技术人员可以参照方法的技术方案实现系统的组成,即,方法中的实施例可理解为构建系统的优选例,在此不予赘述。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变形或修改,这并不影响本发明的实质内容。
Claims (4)
1.一种蒸汽发生器二回路工质汽化起始高度的估计方法,将蒸汽发生器划分为热段、冷段和汽水分离器,其特征在于,包括:
S1,获取给定时刻下蒸汽发生器的实时运行数据;
S2,利用S1得到的所述实时运行数据分别建立热段、冷段的下降通道模型,通过所述下降通道模型得到当前时刻下降通道底部出口液相工质的温度、压力及质量流量;
S3,利用S1得到的所述实时运行数据计算一回路冷却剂与倒U型管金属壁之间传热系数以及倒U型管金属壁与二回路工质之间的传热系数;
S4,利用S1得到的所述实时运行数据以及S3得到的所述一回路冷却剂与倒U型管金属壁之间传热系数,分别建立热段、冷段的一回路冷却剂模型,计算得到倒U型管金属壁的温度分布;
S5,利用S1得到的所述实时运行数据、S3得到的所述倒U型管金属壁与二回路工质之间的传热系数、S4得到的所述倒U型管金属壁的温度分布以及S2得到的所述当前时刻下降通道底部出口液相工质的温度、压力及质量流量,分别建立热段、冷段的上升通道模型,通过所述上升通道模型计算得到当前时刻沿倒U型管高度的二回路工质的流速、温度、压力以及汽化起始高度;其中所述当前时刻沿倒U型管高度的二回路工质的流速、温度、压力包括当前时刻上升通道顶部出口气液混合物工质的流速、温度以及压力;
S6,利用S1得到的所述实时运行数据以及S5得到的所述上升通道顶部出口气液混合物工质的流速、温度以及压力,建立汽水分离器模型,通过所述汽水分离器模型计算得到汽水分离器出口液相工质的温度、压力和质量流量,所述汽水分离器出口液相工质是所述下降通道模型入口再循环水;
S7,利用S6得到的所述汽水分离器出口液相工质的温度、压力和质量流量,重新执行步骤S2~S6若干次,得到最终的汽化起始高度估计值;
式中,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):
式中,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是冷段下降通道当量直径;
通过对冷段下降通道模型求解,得到当前时刻冷段下降通道底部出口液相工质的温度、压力以及质量流量;
设热段一回路冷却剂与倒U型管金属壁之间传热系数为KHL,PS,冷段一回路冷却剂与倒U型管金属壁之间传热系数为KCL,PS,热段预热区倒U型管金属壁与二回路工质之间的传热系数为KHL,RC,PR,冷段预热区倒U型管金属壁与二回路工质之间的传热系数为KCL,RC,PR;
将这4个参数KHL,PS、KCL,PS、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型管金属壁与二回路工质之间的传热系数为KHL,RC,BR,冷段沸腾区倒U型管金属壁与二回路工质之间的传热系数为KCL,RC,BR,将这两个参数KHL,RC,BR、KCL,RC,BR简记为K*,两参数均分别采用式(9)~(14)计算:
K*=Kcht+Kbht (9)
式中,Kcht、Kbht分别是对流传热部分的传热系数和泡核沸腾传热部分的传热系数;CP,w是一回路冷却剂定压比热容;hfs是沸腾区液相工质汽化潜热;σ沸腾区液相工质表面张力系数;△TMT是沸腾区倒U型管金属壁过热度;△PMT是沸腾区饱和蒸汽压差;x是质量气含率;ρw是上升通道液相工质密度;ρs是上升通道饱和蒸汽密度;μw是上升通道液相工质粘性系数;μs是上升通道饱和蒸汽粘性系数;dHL,MT是热段倒U型管内径;G是工质质量流量;Xtt和S是中间变量;
考虑热段一回路冷却剂重力压降,根据动量、质量和能量守恒关系,建立热段一回路冷却剂模型为式(15)~(18):
式中,ρHL,PS是热段一回路冷却剂密度;WHL,PS是热段一回路冷却剂的流速;CP,HL,PS是热段一回路冷却剂定压比热容;THL,PS是热段一回路冷却剂温度;KHL,PS是热段一回路冷却剂通过倒U型管金属壁向二回路工质传热的传热系数;n是倒U型管根数;dHL,MT是热段倒U型管内径;PHL,PS是热段一回路冷却剂压力;g是重力加速度;
通过对热段一回路冷却剂模型求解,得到热段倒U型管金属壁的温度分布;
建立冷段一回路冷却剂模型为式(19)~(22):
式中,ρCL,PS是冷段一回路冷却剂密度;WCL,PS是冷段一回路冷却剂的流速;CP,CL,PS是冷段一回路冷却剂定压比热容;TCL,PS是冷段一回路冷却剂温度;KCL,PS是冷段一回路冷却剂通过倒U型管金属壁向二回路工质传热的传热系数;dCL,MT是冷段倒U型管内径;PCL,PS是冷段一回路冷却剂压力;
通过对冷段一回路冷却剂模型求解,得到冷段倒U型管金属壁的温度分布;
考虑热段上升通道工质重力压降、摩擦压降和加速压降,根据动量、质量和能量守恒关系,建立热段上升通道模型为式(23)~(30):
式中,ρ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):
式中,ρ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是冷段上升通道预热区液相工质压力;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型管高度的二回路工质的流速、温度、压力以及汽化起始高度;
为了汽水分离器模型,计算当前时刻沿倒U型管高度的二回路工质质量气含率分布,其中:
式中,xBR是沸腾区工质质量气含率;hBR是沸腾区气液混合相工质比焓;hss是沸腾区饱和蒸汽比焓;hsw是沸腾区饱和水比焓;hBR,hss和hsw根据沸腾区气液混合相工质温度和压力通过工质物性参数数据库计算得到;
根据上述公式,对应采用热段、冷段的对应参数,计算得到热段上升通道沸腾区出口气液混合相工质质量气含率xHL,RC,BR,out和冷段上升通道沸腾区出口气液混合相工质质量气含率xCL,RC,BR,out;
所述汽水分离器模型为式(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)
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是汽水分离器出口饱和蒸汽质量流量;GHL,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是冷段上升通道沸腾区出口气液混合相工质密度;
通过对汽水分离器模型求解,得到汽水分离器出口饱和水以及饱和蒸汽的温度、压力、质量流量,此处汽水分离器出口饱和水即为再循环水,将饱和水参数代回S2,通过迭代若干次,得到蒸汽发生器二回路工质汽化起始高度。
2.根据权利要求1所述的蒸汽发生器二回路工质汽化起始高度估计方法,其特征在于,所述给定时刻下蒸汽发生器的实时运行数据,包括:
-机组负荷;
-给水温度、压力以及质量流量;
-饱和蒸汽温度、压力以及质量流量;
-一回路冷却剂进出口温度、压力以及质量流量;
-水位高度。
3.根据权利要求1所述的蒸汽发生器二回路工质汽化起始高度估计方法,其特征在于,在所述蒸汽发生器的上升通道中,根据二回路工质状态将上升通道划分为预热区和沸腾区;其中,预热区和沸腾区分界面的划分依据为:
hRC(t,z)=hsw(t,z) (1)
式中,hRC(t,z)是上升通道在当前时刻t下高度为z的二回路工质的比焓;hsw(t,z)是在当前时刻t下高度为z的二回路工质饱和状态的比焓;若在z0高度时,hRC(t,z0)≥hsw(t,z0),则z0即是汽化起始高度。
4.一种蒸汽发生器二回路工质汽化起始高度估计装置,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时用于执行上述权利要求1-3任一项所述的蒸汽发生器二回路工质汽化起始高度估计方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011303312.1A CN112699523B (zh) | 2020-11-19 | 2020-11-19 | 蒸汽发生器二回路工质汽化起始高度估计方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011303312.1A CN112699523B (zh) | 2020-11-19 | 2020-11-19 | 蒸汽发生器二回路工质汽化起始高度估计方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112699523A CN112699523A (zh) | 2021-04-23 |
CN112699523B true CN112699523B (zh) | 2022-04-26 |
Family
ID=75505936
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011303312.1A Active CN112699523B (zh) | 2020-11-19 | 2020-11-19 | 蒸汽发生器二回路工质汽化起始高度估计方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112699523B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU46556U1 (ru) * | 2004-11-09 | 2005-07-10 | Замятин Сергей Аркадьевич | Прямоточный вертикальный парогенератор |
CN102778073A (zh) * | 2012-08-10 | 2012-11-14 | 中石化广州工程有限公司 | 强化气分装置利用余热余压回收丙烯的制冷装置及工艺 |
DE102013021831A1 (de) * | 2013-12-21 | 2015-06-25 | Helmut Pfitzner | Verfahren und Vorrichtung zur Umsetzung eines hocheffizienten thermodynamischen Kreisprozesses (Dampfkreislauf) |
CN105137957A (zh) * | 2015-06-23 | 2015-12-09 | 黄红林 | 一种核电蒸汽发生器控制方法 |
CN107784156A (zh) * | 2017-09-07 | 2018-03-09 | 中国船舶重工集团公司第七〇九研究所 | 核动力装置蒸汽排放系统参数的计算方法 |
CN108469744A (zh) * | 2018-02-11 | 2018-08-31 | 东南大学 | 一种建立核电机组蒸汽发生器机理模型的方法及其系统 |
CN109783936A (zh) * | 2019-01-15 | 2019-05-21 | 国网宁夏电力有限公司电力科学研究院 | 压水堆核电站核岛变工况蓄热增量的计算方法 |
JP2019516057A (ja) * | 2016-10-12 | 2019-06-13 | 李華玉 | シングル作業物質の蒸気連合サイクルと連合サイクル蒸気動力装置 |
CN110887027A (zh) * | 2019-11-19 | 2020-03-17 | 中广核工程有限公司 | 一种立式蒸汽发生器热力性能评估系统和方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108875212A (zh) * | 2018-06-19 | 2018-11-23 | 哈尔滨工程大学 | 核动力装置主冷却剂系统热工水力多尺度分区仿真方法 |
CN110310710B (zh) * | 2019-06-21 | 2020-06-16 | 西安交通大学 | 获得燃料元件矩形窄缝通道内干涸型临界热流密度的方法 |
-
2020
- 2020-11-19 CN CN202011303312.1A patent/CN112699523B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU46556U1 (ru) * | 2004-11-09 | 2005-07-10 | Замятин Сергей Аркадьевич | Прямоточный вертикальный парогенератор |
CN102778073A (zh) * | 2012-08-10 | 2012-11-14 | 中石化广州工程有限公司 | 强化气分装置利用余热余压回收丙烯的制冷装置及工艺 |
DE102013021831A1 (de) * | 2013-12-21 | 2015-06-25 | Helmut Pfitzner | Verfahren und Vorrichtung zur Umsetzung eines hocheffizienten thermodynamischen Kreisprozesses (Dampfkreislauf) |
CN105137957A (zh) * | 2015-06-23 | 2015-12-09 | 黄红林 | 一种核电蒸汽发生器控制方法 |
JP2019516057A (ja) * | 2016-10-12 | 2019-06-13 | 李華玉 | シングル作業物質の蒸気連合サイクルと連合サイクル蒸気動力装置 |
CN107784156A (zh) * | 2017-09-07 | 2018-03-09 | 中国船舶重工集团公司第七〇九研究所 | 核动力装置蒸汽排放系统参数的计算方法 |
CN108469744A (zh) * | 2018-02-11 | 2018-08-31 | 东南大学 | 一种建立核电机组蒸汽发生器机理模型的方法及其系统 |
CN109783936A (zh) * | 2019-01-15 | 2019-05-21 | 国网宁夏电力有限公司电力科学研究院 | 压水堆核电站核岛变工况蓄热增量的计算方法 |
CN110887027A (zh) * | 2019-11-19 | 2020-03-17 | 中广核工程有限公司 | 一种立式蒸汽发生器热力性能评估系统和方法 |
Non-Patent Citations (2)
Title |
---|
压水堆核电站U型管蒸汽发生器建模与仿真;缠阿芳;《中国优秀博硕士学位论文全文数据库(硕士) 工程科技Ⅱ辑》;20130715(第07期);全文 * |
核电机组蒸汽发生器建模与状态空间;魏志伟;《中国优秀博硕士学位论文全文数据库(硕士) 工程科技Ⅱ辑》;20190515(第05期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112699523A (zh) | 2021-04-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Franco et al. | Experimental analysis of closed loop two phase thermosyphon (CLTPT) for energy systems | |
Colombo et al. | A scheme of correlation for frictional pressure drop in steam–water two-phase flow in helicoidal tubes | |
Chung et al. | Thermal hydraulic calculation in a passive residual heat removal system of the SMART-P plant for forced and natural convection conditions | |
Kim et al. | Heat transfer characteristics and operation limit of pressurized hybrid heat pipe for small modular reactors | |
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? | |
Bilde et al. | Mathematical model of natural circulation biomass boilers during start-up | |
CN112417676B (zh) | 基于核电蒸汽发生器分布参数模型的关键变量估计方法 | |
CN112699523B (zh) | 蒸汽发生器二回路工质汽化起始高度估计方法及装置 | |
CN112417780B (zh) | 蒸汽发生器二回路再循环水质量流量估计方法及系统 | |
Tao et al. | Numerical and experimental study on thermal-hydraulic behaviors in the secondary side of passive residual heat removal heat exchanger | |
CN112417781B (zh) | 核电蒸汽发生器出口饱和蒸汽质量流量估计方法及系统 | |
CN112417681B (zh) | 一种蒸汽发生器一二次侧对流换热系数分布的估计方法 | |
CN112380713B (zh) | 一种蒸汽发生器倒u型管金属壁温度分布估计方法 | |
CN112417680B (zh) | 蒸汽发生器二回路工质质量气含率分布估计方法及系统 | |
Ghazanfari et al. | Drift flux modeling of the VVER-1000 horizontal nuclear steam generator | |
CN112417782B (zh) | 一种蒸汽发生器二回路工质循环倍率估计方法 | |
Li et al. | Evaporation flow patterns and heat transfer in copper and stainless steel three-dimensional dimpled tubes | |
Shao et al. | Numerical investigation on the thermal-hydraulic performance of separated structure steam generator with different tube arrangements | |
Sun et al. | Transient analysis and dynamic modeling of the steam generator water level for nuclear power plants | |
DOLGANOV et al. | Advanced Natural Circulation Model in Evaporation Circuits of Thermosiphon Waste Heat Boilers | |
Alipchenkov et al. | Progress in system thermohydraulic code HYDRA-IBRAE/LM models development for fast reactor simulation | |
David et al. | Heat exchanger heat transfer coefficient and CFD modelling | |
CN114781109A (zh) | 核电厂蒸汽发生器离散化模型构建方法、终端及介质 | |
Jin et al. | Thermohydraulic Characteristics of Supercritical CO2 Flowing in Printed Circuit Heat Exchanger With Large Eddy Simulation | |
Wuth | Comparative analysis of the Žukauskas method and data from the OSU MASLWR test facility steam generator |
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 |