核电站蒸汽发生器再循环水质量流量实时估计方法

doi:10.16183/j.cnki.jsjtu.2020.355

核电站蒸汽发生器再循环水质量流量实时估计方法

惠久武¹^,², 凌君¹, 栾振华¹, 王改霞¹, 董贺², 袁景淇^,²

1.中广核工程有限公司核电安全监测技术与设备国家重点实验室, 广东深圳 518172

2.上海交通大学自动化系, 上海 200240

A Real-Time Mass Flow Rate Estimation Method of Recirculation Water in Steam Generator of Nuclear Power Plants

HUI Jiuwu¹^,², LING Jun¹, LUAN Zhenhua¹, WANG Gaixia¹, DONG He², YUAN Jingqi^,²

1. State Key Laboratory of Nuclear Power Safety Monitoring Technology and Equipment, China Nuclear Power Engineering Co., Ltd., Shenzhen 518172, Guangdong, China

2. Department of Automation, Shanghai Jiao Tong University, Shanghai 200240, China

通讯作者: 袁景淇,男,教授,博士生导师,电话(Tel.):021-34204055;E-mail:jqyuan@sjtu.edu.cn.

责任编辑: 石易文

收稿日期: 2020-10-29

基金资助:

核电安全监控技术与装备国家重点实验室开放课题(K-2019.426)

Received: 2020-10-29

作者简介 About authors

惠久武(1996-),男,山东省日照市人,博士生,主要从事核电站关键设备机理建模及先进控制策略研究.

摘要

根据蒸汽发生器的实际结构,将蒸汽发生器划分为热段、冷段、汽水分离器和汽室,热段和冷段又划分为下降通道和上升通道.基于质量、能量、动量守恒定律,依次建立蒸汽发生器热段、冷段、汽水分离器和汽室分布的参数模型.结合工质物性参数数据库和数据采集系统 (DCS) 实时测量数据实时求解模型,从而获得蒸汽发生器再循环水的质量流量估计值.采用某CPR1000核电机组运行数据的验证结果表明,由所提方法解算得到的蒸汽发生器顶部出口饱和蒸汽温度、压力和质量流量与实测值吻合较好,从而间接证实了蒸汽发生器再循环水质量流量估计的正确性.

关键词： 蒸汽发生器; 分布参数模型; 再循环水; 质量流量; 实时估计

Abstract

According to its actual structure, the steam generator is divided into the hot leg, the cold leg, the steam-water separator, and the steam chamber, where the hot leg and the cold leg are further divided into a descending channel and an ascending channel. Based on the mass, energy, and momentum conservation laws, the distribution parameter models for the hot leg, the cold leg, the steam-water separator, and the steam chamber are established. Combining the physical parameter database of the working medium and the data collection system (DCS) real-time measurement data, the model is solved in real time and the estimated mass flow rate of recirculation water is obtained. The verification results using the operating data of a CPR1000 nuclear power unit show that the calculated temperature, the pressure, and the mass flow rate of the outlet saturated steam of the steam generator agree well with the measured values, which indirectly verifies the correctness of the estimation approach of recirculation water proposed in this paper.

Keywords： steam generator; distributed parameter model; recirculation water; mass flow rate; real-time estimation

PDF (1569KB) 元数据多维度评价相关文章导出 EndNote| Ris| Bibtex 收藏本文

本文引用格式

惠久武, 凌君, 栾振华, 王改霞, 董贺, 袁景淇. 核电站蒸汽发生器再循环水质量流量实时估计方法[J]. 上海交通大学学报（自然版）, 2022, 56(1): 21-27 doi:10.16183/j.cnki.jsjtu.2020.355

HUI Jiuwu, LING Jun, LUAN Zhenhua, WANG Gaixia, DONG He, YUAN Jingqi. A Real-Time Mass Flow Rate Estimation Method of Recirculation Water in Steam Generator of Nuclear Power Plants[J]. Journal of shanghai Jiaotong University, 2022, 56(1): 21-27 doi:10.16183/j.cnki.jsjtu.2020.355

核能作为一种清洁、经济、高效的能源,近年来获得了平稳的发展^[1].截至2019年底,我国现役核电机组47台,位居世界第3位,装机容量 44645 MW占全国总装机容量的2.3%,发电量则占全国总发电量的4.2%.

由于蒸汽发生器至关重要的能量枢纽作用,近年来获得了研究人员的关注,例如一、二次侧耦合传热特性分析^[2]、传热管道破裂事故源项的计算分析^[3],蒸汽发生器热工水力数值研究^[4,5,6],蒸汽发生器水位控制^[7,8,9],蒸汽发生器换热管结构应力分析^[10],蒸汽发生器传热管泄漏程度诊断研究^[11],蒸汽发生器机理建模与仿真^{[12,13,14,15]}等.其中,蒸汽发生器机理的建模与仿真是研究热点.文献[12]采用分区与集总参数结合的方式对蒸汽发生器进行机理建模,建立了蒸汽发生器的非线性动态模型.文献[13]采用热工水力软件对大型钠冷快堆核电站蒸汽发生器进行了模型化研究.文献[14]建立了自然循环蒸汽发生器四阶动态水位机理模型,在此基础上设计了串级比例积分微分 (PID) 水位控制系统.文献[15]采用模块化机理建模方法建立了立式自然循环蒸汽发生器的非线性数学模型.

蒸汽发生器再循环水流量是机理建模所涉及的关键变量之一,但在实际工程应用中,该变量无直接测量值.为此,本文以某CPR1000核电机组倒U型自然循环蒸汽发生器为研究对象,基于质量、能量、动量守恒定律,建立蒸汽发生器的分布参数模型,再结合工质物性参数数据库和数据采集系统(DCS)实时测量数据进行模型求解,得到蒸汽发生器再循环水质量流量的实时解算值.基于上述模型化方法计算得到的蒸汽发生器顶部出口饱和蒸汽温度、压力和质量流量与实测值吻合得较好.此外,本文对蒸汽发生器运行状态进行监测,为设备结构改进和运行优化奠定了基础.

1 蒸汽发生器分布参数模型

1.1 蒸汽发生器的结构简化

立式自然循环蒸汽发生器的简化结构可以划分为热段、冷段、汽水分离器和汽室4部分,如图1所示.其中,热段和冷段又划分为下降通道和上升通道,下降通道是指外壳与内部套筒之间工质流经的空间,上升通道是指内部套筒与倒U型管壁之间工质流经的空间.此外,以下降通道底部为坐标原点,竖直向上为z轴建立坐标系.

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1 核电站蒸汽发生器结构简化图

Fig.1 Simplified structure of steam generator in nuclear power plant

根据二回路工质状态将上升通道划分为预热区和沸腾区.预热区和沸腾区分界面的划分依据为

(1)

\begin{matrix} h_{RC} (t, z) = h_{sw} (t, z) \end{matrix}

式中: $h_{RC} (t, z) 为上升通道 t 时刻沿高度 z 的二回路工质的比焓; h_{sw} (t, z) 为 t 时刻沿高度 z$ 的二回路相应温度下的工质饱和状态比焓.

1.2 分布参数模型建模

假设给水和再循环水在下降通道入口立即混合,并按设计要求,下降通道入口二回路液相工质中比例为w₁%的给水流入热段,比例为(100-w₁)%的给水流入冷段,比例为w₂%的再循环水流入热段,比例为(100-w₂)%的再循环水流入冷段.根据动量、质量和能量守恒原理,可得热段下降通道模型为

(2)

\frac{d m_{H_D}}{d t} = \frac{d ρ_{H_D} A_{H_D} H}{d t} = w_{1} % g_{m, fw} + w_{2} % g_{m, rw} - g_{m, H_D_o}

(3)

\frac{d m_{H_D} c_{H_D} T_{H_D}}{d t} = \frac{d ρ_{H_D} A_{H_D} H h_{H_D}}{d t} = w_{1} % g_{m, fw} h_{fw} + w_{2} % g_{m, rw} h_{rw} - g_{m, H_D_o} h_{H_D_o}

(4)

- \frac{d p_{H_D}}{d z} = \frac{d g_{m, H_D}}{d t} + \frac{f_{H_D} g_{m, H_D}^{2}}{2 ρ_{H_D} D_{H_D}} + ρ_{H_D} g

式中:m_{H_D}为热段下降通道液相工质质量;ρ_{H_D}为热段下降通道液相工质密度;A_{H_D}为热段下降通道的横截面积;H为下降通道的水位高度;g_m_,fw为给水质量流量; g_m_,rw为再循环水质量流量;g_m_{,H_D_o}为热段下降通道出口液相工质质量流量; c_{H_D}为热段下降通道液相工质的定压比热容;T_{H_D}为热段下降通道液相工质温度;h_fw为给水比焓;h_rw为再循环水比焓;h_{H_D}为热段下降通道液相工质比焓;h_{H_D_o}为热段下降通道出口液相工质比焓;p_{H_D}为热段下降通道液相工质压力; g_m_{,H_D}为热段下降通道液相工质质量流量;f_{H_D}为热段下降通道摩擦因子;D_{H_D}为热段下降通道当量直径;g为重力加速度.

考虑热段上升通道二回路工质重力压降、摩擦压降和加速压降,根据动量、质量和能量守恒原理,所得热段上升通道模型为

(5)

\frac{\partial ρ_{H_R}}{\partial t} + \frac{\partial ρ_{H_R} v_{H_R}}{\partial z} = 0

(6)

\frac{\partial ρ_{H_M} c_{H_M} T_{H_M}}{\partial t} = n K_{H_R_P} π d_{H_M} (T_{H_R_P} - T_{H_M})

(7)

\frac{\partial ρ_{H_R_P} c_{H_R_P} T_{H_R_P}}{\partial t} + \frac{\partial ρ_{H_R_P} c_{H_R_P} T_{H_R_P} v_{H_R_P}}{\partial z} = n K_{H_R_P} π d_{H_M} (T_{H_M} - T_{H_R_P})

(8)

\frac{\partial ρ_{H_M} c_{H_M} T_{H_M}}{\partial t} = n K_{H_R_B} π d_{H_M} (T_{H_R_B} - T_{H_M})

(9)

\frac{\partial ρ_{H_R_B} c_{H_R_B} T_{H_R_B}}{\partial t} + \frac{\partial ρ_{H_R_B} c_{H_R_B} T_{H_R_B} v_{H_R_B}}{\partial z} = n K_{H_R_B} π d_{H_M} (T_{H_M} - T_{H_R_B})

(10)

- \frac{\partial p_{H_R_P}}{\partial z} = \frac{\partial g_{m, H_R_P}}{\partial t} + \frac{\partial}{\partial z} (\frac{g_{m, H_R_P}^{2}}{ρ_{H_R_P}}) + \frac{f_{H_R_P} g_{m, H_R_P}^{2}}{2 ρ_{H_R_P} D_{H_R_P}} ϕ^{2} + ρ_{H_R_P} g + ξ_{H_R_P} \frac{g_{m, H_R_P}^{2}}{ρ_{H_R_P}}

(11)

- \frac{\partial p_{H_R_B}}{\partial z} = \frac{\partial g_{m, H_R_B}}{\partial t} + \frac{\partial}{\partial z} (\frac{g_{m, H_R_B}^{2}}{ρ_{H_R_B}}) + \frac{f_{H_R_B} g_{m, H_R_B}^{2}}{2 ρ_{H_R_B} D_{H_R_B}} ϕ^{2} + ρ_{H_R_B} g + ξ_{H_R_B} \frac{g_{m, H_R_B}^{2}}{ρ_{H_R_B}}

(12)

ϕ^{2} = [1 + x (\frac{ρ_{w}}{ρ_{s}} - 1)] [1 + x (\frac{μ_{w}}{μ_{s}} - 1)]

式中:ρ_{H_M}为热段倒U型管金属壁密度;ρ_{H_R}为热段上升通道工质密度;ρ_{H_R_P}为热段上升通道预热区液相工质密度;ρ_{H_R_B}为热段上升通道沸腾区气液混合相工质密度;v_{H_R}为热段上升通道工质的流速;v_{H_R_P}为热段上升通道预热区液相工质的流速;v_{H_R_B}为热段上升通道沸腾区气液混合相工质的流速;c_{H_M}为热段倒U型管金属壁的定压比热容;c_{H_R_P}为热段上升通道预热区液相工质的定压比热容;c_{H_R_B}为热段上升通道沸腾区气液混合相工质的定压比热容;D_{H_R_P}为热段上升通道预热区水力直径;D_{H_R_B}为热段上升通道沸腾区水力直径;T_{H_M}为热段倒U型管金属壁温度;T_{H_R_P}为热段上升通道预热区液相工质温度;T_{H_R_B}为热段上升通道沸腾区气液混合相工质温度;n为倒U型管数目;p_{H_R_P}为热段上升通道预热区液相工质压力;p_{H_R_B}为热段上升通道沸腾区气液混合相工质压力;g_m_{,H_R_P}为热段上升通道预热区液相工质的质量流量;g_m_{,H_R_B}为热段上升通道沸腾区气液混合相工质的质量流量;f_{H_R_P}为热段上升通道预热区摩擦因数;f_{H_R_B}为热段上升通道沸腾区摩擦因数;ξ_{H_R_P}为热段上升通道预热区局部阻力系数;ξ_{H_R_B}为热段上升通道沸腾区的局部阻力系数;d_{H_M}为热段倒U型管内径;ϕ为两相倍乘因子;x为质量气含率;ρ_w为上升通道液相工质密度;ρ_s为上升通道饱和蒸汽密度;μ_w为上升通道液相工质黏性系数;μ_s为上升通道饱和蒸汽黏性系数;K_{H_R_P}为热段上升通道预热区二回路工质与倒U型管金属壁间的传热系数;K_{H_R_B}为热段上升通道沸腾区二回路工质与倒U型管金属壁间的传热系数.

考虑热段一回路冷却剂重力压降,根据动量、质量和能量守恒原理,可得热段一回路模型为

(13)

\frac{\partial ρ_{H_P}}{\partial t} + \frac{\partial ρ_{H_P} v_{H_P}}{\partial z} = 0

(14)

\frac{\partial ρ_{H_P} c_{H_P} T_{H_P}}{\partial t} + \frac{\partial ρ_{H_P} c_{H_P} T_{H_P} v_{H_P}}{\partial z} = K_{H_P} π d_{H_M} (T_{H_M} - T_{H_P})

(15)

\frac{\partial ρ_{H_M} c_{H_M} T_{H_M}}{\partial t} = K_{H_P} π d_{H_M} (T_{H_P} - T_{H_M})

(16)

\frac{\partial p_{H_P}}{\partial z} = - ρ_{H_P} g

式中:ρ_{H_P}为热段一回路冷却剂的密度;v_{H_P}为热段一回路冷却剂的流速;c_{H_P}为热段一回路冷却剂的定压比热容; T_{H_P}为热段一回路冷却剂的温度; p_{H_P}为热段一回路冷却剂的压力;K_{H_P}为热段一回路冷却剂通过倒U型管金属壁向二回路工质传热的传热系数.

一回路冷却剂向预热区二回路工质的传热为管内强制对流换热,采用Dittus-Boelter公式计算传热系数^[16,17],则有:

(17)

\begin{matrix} K_{PR} = 0.023 R e_{w}^{0.8} P r_{w}^{0.3} λ_{w} / d_{H_M} \end{matrix}

式中:K_PR为预热区一回路冷却剂通过倒U型管金属壁向预热区二回路工质传热的传热系数;Re_w为一回路冷却剂雷诺数;Pr_w为一回路冷却剂普朗特数;λ_w为一回路冷却剂热导率.

对二回路沸腾区采用Chen公式计算传热系数^[18,19],则有:

(18)

K_{BR} = K_{cht} + K_{bht}

(19)

K_{cht} = 0.023 {[\frac{g_{m} (1 - x) d_{H_M}}{μ_{w}}]}^{0.8} P r_{w}^{0.4} \frac{λ_{w}}{d_{H_M}} F

(20)

F = \{\begin{array}{l} 1.0, & \frac{1}{X_{tt}} \leq 0.1 \\ 2.35 {(\frac{1}{X_{tt}} + 0.213)}^{0.736}, & \frac{1}{X_{tt}} > 0.1 \end{array}

(21)

X_{tt} = {(\frac{1 - x}{x})}^{0.9} {(\frac{ρ_{s}}{ρ_{w}})}^{0.5} {(\frac{μ_{s}}{μ_{w}})}^{0.1}

(22)

K_{bht} = 0.00122 \frac{ρ_{w}^{0.49} c_{w}^{0.45} λ_{w}^{0.79}}{ρ_{s}^{0.24} h_{fs} σ^{0.49} μ_{w}^{0.29}} Δ T_{MT}^{0.24} Δ p_{MT}^{0.75} S

(23)

S = \{\begin{array}{l} \frac{1}{1 + 3.305 \times 10^{- 6} R e_{w}^{1.14}}, & R e_{w} \leq 3.25 \\ \frac{1}{1 + 0.42 R e_{w}^{0.78}}, & 32.5 < R e_{w} < 70 \\ \frac{1}{10}, & R e_{w} > 70 \end{array}

式中:K_BR、K_cht、K_bht分别为沸腾区、沸腾区对流传热部分和沸腾区泡核沸腾传热部分的传热系数;g_m为质量流量;σ为沸腾区液相工质表面张力系数;c_w为一回路冷却剂定压比热容;h_fs为沸腾区液相工质的汽化潜热;ΔT_MT为沸腾区倒U型管金属壁过热度;Δp_MT为沸腾区饱和蒸汽压差.

汽水分离器稳态模型为

(24)

g_{m, ss, S_o} = η (x_{H_R_B_o} g_{m, H_R_B_o} + x_{C_R_B_o} g_{m, C_R_B_o})

(25)

g_{m, sw, S_o} = (1 - η x_{H_R_B_o}) g_{m, H_R_B_o} + (1 - η x_{C_R_B_o}) g_{m, C_R_B_o}

(26)

p_{S_o} = p_{S_i} - ξ_{SP} \frac{g_{m, S_i}^{2}}{ρ_{S_i}}

(27)

g_{m, S_i} = g_{m, H_R_B_o} + g_{m, C_R_B_o}

(28)

p_{S_i} = p_{H_R_B_o} = p_{C_R_B_o}

(29)

ρ_{S_i} = ρ_{H_R_B_o} = ρ_{C_R_B_o}

式中:g_m_{,H_R_B_o}为热段上升通道沸腾区出口气液混合相工质的质量流量;g_m_{,C_R_B_o}为冷段上升通道沸腾区出口气液混合相工质的质量流量;g_m_{,ss,S_o}为汽水分离器出口饱和蒸汽的质量流量;g_m_{,sw,S_o}为汽水分离器出口饱和水的质量流量;g_m_{,S_i}为汽水分离器入口气液混合相工质质量流量;x_{H_R_B_o}为热段上升通道沸腾区出口气液混合相工质的质量气含率;x_{C_R_B_o}为冷段上升通道沸腾区出口气液混合相工质的质量气含率;η为汽水分离器效率; p_{H_R_B_o}为热段上升通道沸腾区出口气液混合相工质压力;p_{C_R_B_o}为冷段上升通道沸腾区出口气液混合相工质压力;p_{S_i}为汽水分离器入口气液混合相工质压力;p_{S_o}为汽水分离器出口饱和蒸汽压力;ρ_{H_R_B_o}为热段上升通道沸腾区出口气液混合相工质密度;ρ_{C_R_B_o}为冷段上升通道沸腾区出口气液混合相工质密度;ξ_SP为汽水分离器局部阻力系数;ρ_{S_i}为汽水分离器入口气液混合相工质密度.

汽室动态模型为

(30)

\begin{matrix} τ \frac{d g_{m, SC_o}}{d t} = g_{m, SC_i} - g_{m, SC_o} \end{matrix}

式中:g_m_{,SC_i}为汽室入口饱和蒸汽质量流量;g_m_{,SC_o}为汽室出口饱和蒸汽质量流量;τ为时间常数.

2 模型验证

模型验证数据来自某 1000 MW核电机组.该核电机组蒸汽发生器设计参数如表1所示.分布式控制系统实时数据采样时间间隔Δt=5 s.数据包括机组负荷、给水温度、给水压力、给水质量流量、饱和蒸汽温度、饱和蒸汽压力、饱和蒸汽质量流量、一回路冷却剂进出口温度、一回路冷却剂进出口压力、一回路冷却剂进出口质量流量以及一回路冷却剂进出口水位高度等.由于给水压力和一回路冷却剂压力及质量流量在核电站变负荷过程中变化较小,所以取其测量均值.给水压力取为6.7 MPa, 一回路冷却剂压力取为15.4 MPa, 一回路冷却剂体积流量取为 23790 m³/h.基于工质物性参数数据库和蒸汽发生器结构参数库,采用Runge-Kutta 法解算热段模型、冷段模型、汽水分离器模型和汽室模型,输出蒸汽发生器再循环水质量流量估计值.输出的结果作为下降通道模型的输入,进行下一步解算.

表1 某1000 MW核电机组蒸汽发生器设计参数

Tab.1 Rated parameters of steam generator in a 1000 MW nuclear power unit

参数	取值
机组额定功率/MW	1000
倒U型管数目	4474
倒U型管计算高度/m	9.36
倒U型管平均直径/mm	17.96
倒U型管厚度/mm	1.09
w₁/%	80
w₂/%	50

新窗口打开| 下载CSV

核电机组于2019年5月10日15:00—16:00的输出功率变化曲线及水位高度实测值如图2所示,其中: P为机组输出功率,其间由 1030 MW下降到690 MW.蒸汽发生器实测给水温度、一回路冷却剂进出口温度如图3所示,其中:T为温度.在该时段出口饱和蒸汽质量流量、压力和温度的模型仿真值与实测值对比如图4和5所示.其中:T_{SC_o}为汽室出口饱和蒸汽温度;p_{SC_o}为汽室出口饱和蒸汽压力.由图4和5可知,蒸汽发生器出口饱和蒸汽的温度、压力和质量流量最大相对误差分别为1.41%, 0.17% 和1.72%, 从工程上看这些误差已经很小,从而验证了所建立模型的准确性.给水质量流量实测值和再循环水质量流量的实时估计值如图6所示.由图6可知,再循环水流量随着核电机组负荷的降低逐渐下降.

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2 某1000 MW核电机组输出功率及蒸汽发生器水位高度实测值

Fig.2 Output power of a 1000 MW nuclear power unit and measured water level of steam generator

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3 蒸汽发生器工质温度测点

Fig.3 Measured temperatures of working mediums in steam generator

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4 蒸汽发生器出口饱和蒸汽质量流量

Fig.4 Mass flow rate of saturated steam at outlet of steam generator

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5 蒸汽发生器出口饱和蒸汽压力和温度

Fig.5 Pressure and temperature of saturated steam at outlet of steam generator

图6

新窗口打开| 下载原图ZIP| 生成PPT

图6 蒸汽发生器给水质量流量实测值和再循环水质量流量估计值

Fig.6 Measured mass flow rate of feedwater and estimated mass flow rate of recirculation water of steam generator

由式(1)~(30)构成的蒸汽发生器分布参数系统是机理模型,比较复杂,但对如图2所示的 3600 s 时间区间和5 s仿真步长,采用i7 PC机(4CPU,主频3.4 GHz)所需的仿真时间是994 s,即每步仿真消耗的计算时间为1.38 s,故所提模型的应用具有实时性保障.

3 结语

本文对蒸汽发生器再循环水质量流量估计问题进行了研究.首先,根据蒸汽发生器的实际几何结构对其进行结构简化,将其划分为热段、冷段、汽水分离器和汽室4部分.然后,基于质量、动量和能量守恒定律,建立所划分的热段、冷段、汽水分离器和汽室分布参数模型.最后,结合工质物性参数数据库和某核电机组分布式控制系统的实时测量数据进行模型解算.基于模型仿真得到的蒸汽发生器出口饱和蒸汽温度、压力和质量流量的计算值与实测值吻合得较好,从而验证了本文所建立模型的准确性,也间接说明了作为实时仿真输出变量之一的再循环水质量流量实时估计的准确性.

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子

[1]

刘亮

从“十三五”规划看电力转型发展

[J]. 大众用电, 2019, 34(4):3-4.