靖剑平　刘雅宁　贾　斌　高新力　孙　微　左嘉旭　张春明

1（环境保护部核与辐射安全中心　北京　100082）2（中国科学院过程工程研究所　北京　100190）

10 MWt固态燃料熔盐堆控制棒失控抽出事故分析

靖剑平1刘雅宁2贾斌1高新力1孙微1左嘉旭1张春明1

1（环境保护部核与辐射安全中心北京100082）2（中国科学院过程工程研究所北京100190）

钍基熔盐堆核能系统项目是中科院先导科技专项之一，其战略性目标是研发第四代熔盐冷却裂变反应堆核能系统。基于10 MWt固态燃料熔盐堆的系统设计，开发了适用于球床式反应堆系统的安全分析软件，并以高温气冷堆为对象对程序计算结果的准确性进行了验证。基于该软件程序，对固态燃料球床堆（Thorium Molten Salt Reactor-Solid Fuel, TMSR-SF）控制棒失控抽出事故进行了分析计算，研究了不同停堆限值及各停堆信号对事故的影响。计算结果表明，超功率停堆限值越高，出口温度限值越大，信号延迟时间越长，反应堆停堆越晚，堆芯功率和燃料最高温度越高。在TMSR-SF控制棒失控抽出事故下，燃料最高温度不超过860°C，远低于1600°C的熔化温度限值。

固态燃料熔盐堆，安全分析，控制棒失控抽出事故，程序开发

核能具有绿色、高效、低碳排放和可规模生产的突出优势，目前，核能在世界范围内得到了较大重视，大力发展核能已成为我国能源中长期发展规划的重点。由于钍基熔盐堆具有本征安全性、核燃料长期稳定供应、核废物最小化、物理防核扩散、多用途与灵活性等特点，已经成为第四代反应堆核能系统的6种候选堆型之一，受到了国际上的广泛关注［1−2］。2011年中国科学院设置了4个A类先导专项，其中“未来先进核裂变能”项目包括“钍基熔盐（Thorium Molten Salt Reactor, TMSR）核能系统”、“加速器驱动的次临界堆（Accelerator Driven sub Critical Reactor, ADS）嬗变系统”两大内容。在TMSR项目中，分为固态燃料球床堆（Thorium Molten Salt Reactor-Solid Fuel, TMSR-SF）和液态燃料熔盐堆（Thorium Molten Salt Reactor-Liquid Fuel, TMSR-LF）［3−4］。

控制棒失控抽出事故是一个非常重要的反应性引入事故。反应堆控制系统失效或操作人员失误造成一根控制棒失控提出，堆功率迅速上升，引起堆芯燃料温度上升，影响反应堆安全。因此，开展TMSR-SF控制棒失控抽出事故分析具有重要意义。

本文借鉴了高温气冷堆安全分析程序所采用的物理模型和数学方法，结合FLiBe熔盐的物性参数和经验公式，开发了可用于TMSR-SF反应性事故分析的瞬态分析程序，并开展了TMSR-SF控制棒失控抽出事故分析和敏感性分析，为TMSR-SF的设计和安全分析提供支持。

1　TMSR-SF简介

TMSR-SF堆芯功率为10 MWt，活性区为圆柱形固定球床堆芯，一次装料11 043颗燃料元件，燃料元件为全陶瓷包覆颗粒球形燃料元件。堆芯活性区体积1.95m3，堆内包括16根控制棒、一个中子源通道和三个实验测量通道［5−6］。反应堆一回路熔盐为FliBe （LiF-BeF2），二回路熔盐为FliNaK（LiF-NaF-KF）。反应堆一回路由熔盐泵、换热器、反应堆堆芯以及连接管道组成。在正常工况下，堆芯热传递给一回路中的熔盐，通过主换热器，将热量传递给二回路的熔盐，冷却后的一回路熔盐经主熔盐泵送回堆芯，如此循环往复。一回路熔盐进口温度为600°C，出口温度为628°C，质量流速为41.3kg·s−1。TMSR-SF系统如图1所示［7］。

图1　TMSR-SF系统示意图Fig.1　Schematic of TMSR-SF.

2　物理模型和数值方法

2.1物理模型

2.1.1中子动力学模型

采用的反应堆系统分析及事故分析程序过程中经常采用的点堆模型来计算堆芯物理参数及堆芯热功率，忽略空间有关的动力学效应，考虑了6组缓发中子的影响。模型如下：

式中：n（t）为中子数量；β为缓发中子的总份额；ρ（t）为反应性；Λ为中子代时间；λ为缓发中子衰变常数；C（t）为缓发中子先驱核的浓度；i表示第i组缓发中子。

2.2.2流体流动传热模型

TMSR-SF采用FLiBe熔盐作为其冷却剂，在正常及事故工况下FLiBe熔盐发生沸腾相变的可能很小，因此将其作为单相流体处理。一维单相流体瞬态流动传热的基本控制方程如下。

质量守恒方程：

动量守恒方程：

能量守恒方程：

式中：ρ为流体密度；u为流体速度；p为流体压力；h为流体焓值；q为壁面传热量；g为重力加速度；Δpf为流体流动时壁面摩擦压降；Δpr为流体流动时局部压降；z代表流体流动距离；t为时间。

2.2.3辅助计算方程

在求解用以模拟上述瞬态流动传热过程的方程时，还应补充大量的辅助方程或本构方程，其中包括壁面摩擦力方程、壁面传热计算方程等。

1） 壁面摩擦力方程

球床堆芯处壁面摩擦力计算采用德国安全导则KTA3102.3规定的公式计算。具体方程如下：

固态球床熔盐摩擦力计算采用ERGUN公式：

2） 壁面传热计算方程

堆芯燃料元件与熔盐冷却剂间的传热过程的计算中采用WAKAO公式：

2.2数值方法

程序中用于流体瞬态传热流动过程、壁面传热过程计算的基本方程均为偏微分方程，无法直接求解，因此采用有限差分方式进行求解。对TMSR-SF的系统特性模拟采用了一维方法，在各瞬态参数（温度、压力等）对空间相离散过程中仅需要将其在流动方向进行离散。在对空间相离散过程中，采用了交错网格技术。交错网格技术中流体物性状态参数和流动参数存放位置不同，相互之间交错半个控制体，即将压力、温度、焓值等物性存放在控制体中心处，而将流体的速度存放在控制体边界处，如图2所示。对于在控制体边界处的流体物性状态参数采用迎风格式，即采用边界上游控制体的流体物性参数。

图2　网格划分Fig.2　Grid partitioning.

1） 质量守恒方程

2） 动量守恒方程

3） 能量守恒方程

3　程序验证

由于熔盐堆和高温气冷堆具有相同的燃料元件、相似的堆芯结构、相似的运行温度条件，因此选择高温气冷堆为对象对所开发的程序进行验证，以确定程序的计算精度和正确性。验证计算选择的工况为反应性引入事故，对比的计算结果来自高温气冷堆安全分析报告结果［8］。

图3（a）和（b）分别给出了开发程序计算与安全分析报告中得到的高温气冷堆一根控制棒在100%功率运行工况下失控提升事故中堆芯功率、温度反应性（Moderator temperature coefficient, MTC）反馈随时间变化规律。由图3可知，程序计算结果与安全分析报告结果符合性较好，变化规律也一致，从而验证了程序用于分析TMSR-SF控制棒失控抽出事故的可行性。

图3　功率（a）、反应性（b）随时间变化Fig.3　Power （a） and MTC （b） vs. time.

4　控制棒失控抽出事故分析

4.1工况划分和基本假设

目前TMSR-SF尚处于设计阶段，公开资料及文献中关于TMSR-SF反应堆保护系统动作定值的描述比较少，为更加全面、更加深入地研究和探讨TMSR-SF反应堆的安全特性，在控制棒失控抽出事故分析中，考虑了各反应堆保护系统动作定值、系统设备特性参数（如探测信号响应延迟等）对其事故的影响，具体工况划分见表1。

表1　工况划分Table 1　Condition definition.

针对TMSR-SF控制棒失控抽出事故进行了分析计算所采用的主要假设如下：

1） 初始功率取额定功率，即100%功率；

2） 一根控制棒失控提出堆芯引入的总反应性分别为0.01；

3） 反射层控制棒以1cm·s−1的速度提出，反应性线性引入；

4） 考虑最长的停堆延迟时间和反应性当量最大的控制棒卡在堆外，控制棒下落引入的总的反应性为0.1；

5） 控制棒落棒时间为5s。

4.2结果分析

4.2.1停堆功率限值的影响

如表1所示，分别假设反应堆中子注量率/反应堆核功率超过额定值20%和30%后延迟1 s触发反应堆紧急停堆。

由图4（a）和（b）可知，控制棒失控抽出事故下反应堆核功率很快达到停堆限值，并触发反应堆紧急停堆。由于控制棒下落速度较快，反应堆功率迅速降低，反应堆在极短时间内实现停堆，有效确保了反应堆的安全性。由于从事故发生至触发反应堆紧急停堆之间的时间较短，反应堆燃料最高温度峰值经历短暂上升后迅速降低，其峰值分别为692.5°C和697.6°C，远低于1600°C的燃料熔化温度限值。从工况1和2的计算结果对比可以看出，反应堆超功率停堆限值越高，停堆触发时间越晚，堆芯出口冷却剂温度越高，反应堆燃料最高温度峰值越大。

图4　功率（a）和燃料最高温度（b）随时间变化Fig.4　Power （a） and maximum fuel temperature （b） vs. time.

4.2.2不同停堆信号的影响

对于停堆信号的影响，主要研究了由反应堆高功率限值引起的停堆与反应堆出口温度高引起的停堆（具体见表1中工况1和3）。由于温度信号的探测和响应时间相对较长，计算中假设反应堆在出口熔盐温度达到整定值后8s触发紧急停堆。

由图5（a）和（b）可知，堆芯功率和燃料最高温度都呈现先升高后降低的趋势，工况1的燃料最高温度为692.5°C，工况3的燃料最高温度846°C，低于燃料熔化温度限值。由于控制棒失控抽出事故下反应堆出口熔盐温度达到停堆限值的时间相对较晚，工况3比工况1晚停堆42s，从而导致工况3中反应堆总的反应性较工况1大，使得堆芯功率和燃料最高温度峰值也相对较大。

图5　功率（a）和燃料最高温度（b）随时间变化Fig.5　Power （a） and maximum fuel temperature （b） vs. time

4.2.3反应堆出口温度限值的影响

分别选择反应堆出口温度限值为645 °C（工况3）和650 °C（工况4）时进行分析。从图6（a）和（b）可以发现，两种工况下堆芯功率和燃料最高温度变化规律基本一致，均为快速升高达到峰值后迅速下降，两种工况下的燃料温度峰值低于熔化限值。通过对比工况3和4的计算结果可以看出，停堆温度限值越大，停堆触发时间越晚，堆芯出口冷却剂温度峰值越高，反应堆燃料最高温度峰值也越大。

图6　功率（a）和燃料最高温度（b）随时间变化Fig.6　Power （a） and maximum fuel temperature （b） vs. time

4.2.4信号触发延迟时间的影响

对于延迟时间的影响，分别考虑反应堆温度限值达到延迟时间是8s和5s时的情况。由图7（a）和（b）可知，两种工况下堆芯功率和燃料最高温度变化规律一致，只是随着延迟时间的缩短，堆芯功率和燃料最高温度略有降低。两种工况下的燃料温度峰值低于熔化限值。

图7　功率（a）和燃料最高温度（b）随时间变化Fig.7　Power （a） and maximum fuel temperature （b） vs. time

5　结语

针对TMSR-SF系统，开发了适用于球床式反应堆系统的安全分析软件，并以高温气冷堆为对象对程序计算结果的准确性进行了验证。利用程序对TMSR-SF控制棒失控抽出事故进行了分析计算，研究了其瞬态特性及各停堆信号对其安全性的影响。计算结果表明：

1） TMSR-SF控制棒失控抽出事故下，反应堆功率迅速升高，出口温度和燃料最高温度也随之上升，随着反应堆停堆，功率下降，出口温度和燃料最高温度也迅速下降，最终达到稳定。事故过程中，燃料最高温度远低于1600°C的熔化温度限值。

2） 超功率停堆限值越高，出口温度限值越大，信号延迟时间越长，反应堆停堆越晚，堆芯功率和燃料峰值温度越高。

3） 与超功率停堆限值相比，出口温度限值达到时间更晚，堆芯功率和燃料峰值温度更高。

1 Ingersoll D T, Forsberg C W, Ott L J, et al. Status of preceonceptual design of the Advanced High-Temperature Reactor （AHTR）［R］. Knoxville, Tennessee: Oak Ridge National Laboratory, ORNL/TM-2004/14, 2004

2 Varma V K, Holcomb D E, Peretz F J, et al. AHTR mechanical, structural, and neutronic precocenptual design［R］. Knoxville, Tennessee: Oak Ridge National Laboratory, ORNL/TM-2012/320, 2012

3 TMSR研究中心. 2 MWt固态钍基熔盐实验堆概念设计报告（上）［R］. 上海: 中国科学院上海应用物理研究所, 2013

TMSR Research Center. The conceptual design report 2 MWt solid thorium-based molten salt experiment reactor （Upper V01）［R］. Shanghai: Shanghai Institute of Applied Physics, Chinese Academy of Sciences, 2013

4 田金, 夏晓彬, 彭超, 等. 10-MWt固态钍基熔盐堆乏燃料贮存系统临界安全影响分析［J］. 核技术, 2015, 38（5）: 050602. DOI: 10.11889/j.0253-3219.2015.hjs.38.050602

TIAN Jin, XIA Xiaobin, PENG Chao, et al. Impact analysis of criticality safety for 10-MWt solid thorium-based molten salt reactor spent nuclear fuel storage system［J］. Nuclear Techniques, 2015, 38（5）: 050602. DOI: 10.11889/j.0253-3219.2015.hjs.38.050602

5 TMSR-SFl堆物理分总体. 10 MWt TMSR-SFl总体物理方案和参数［R］. 上海: 中国科学院上海应用物理研究所, 2013

Reactor Physics Department of TMSR-SFI. Physics programs and parameters for 10-MWt TMSR-SF1［R］. Shanghai: Shanghai Institute of Applied Physics, Chinese Academy of Sciences, 2013

6 焦小伟, 王凯, 何兆忠, 等. 固态熔盐堆全厂断电ATWS事故工况下的堆芯安全探讨［J］. 核技术, 2015, 38（2）: 020604. DOI: 10.11889/j.0253-3219.2015.hjs.38. 020604

JIAO Xiaowei, WANG Kai, HE Zhaozhong, et al. Core safety discussion under station blackout ATWS accident of solid fuel molten salt reactor［J］. Nuclear Techniques, 2015, 38（2）: 020604. DOI: 10.11889/j.0253-3219.2015. hjs.38.020604

7 梅牡丹, 邵世威, 何兆忠, 等. 固态钍基熔盐堆概率安全评价始发事件分析研究［J］. 核技术, 2014, 37（9）: 090601. DOI: 10.11889/j.0253-3219.2014.hjs.37.090601

MEI Mudan, SHAO Shiwei, HE Zhaozhong, et al. Research on initial event analysis for solid thorium molten salt reactor probabilistic safety assessment［J］. Nuclear Techniques, 2014, 37（9）: 090601. DOI: 10.11889/ j.0253-3219.2014.hjs.37.090601

8 华能山东石岛湾核电有限公司. 华能山东石岛湾核电厂高温气冷堆核电站示范工程初步安全分析报告［R］.北京: 中国华能集团公司, 2008

Huaneng Shidao Shandong Bay Nuclear Power Co., Ltd. Preliminary safety analysis report on the demonstration project of the high temperature gas cooled reactor nuclear power plant in Shidao Bay, Shandong, China［R］. Beijing: China Huaneng Group, 2008

Accident analyses of uncontrolled rod cluster control assembly bank withdrawal for 10-MWt thorium-based molten salt reactor-solid fuel

JING Jianping1LIU Yaning2JIA Bin1GAO Xinli1SUN Wei1ZUO Jiaxu1ZHANG Chunming1

1（Nuclear and Radiation Safety Center, Ministry of Environmental Protection, Beijing 100082, China）
2（Institute of Process Engineering, Chinese Academy of Sciences, Beijing 100190, China）

Background: Thorium-based molten salt reactor nuclear power system project is one of the “Strategic Priority Program” of Chinese Academy of Sciences, and its strategic goal is to develop the fourth generation moltensalt cooled fission reactor nuclear power system. Based on the system design of 10-MWt molten salt reactor-solid fuel, a safety analysis software is developed for the pebble bed reactor system and verified the correctness and accuracy for the high-temperature gas-cooled reactor of results calculated via this software program. Purpose: This study aims at accident analyses of uncontrolled rod cluster control assembly bank withdrawal for thorium molten salt reactor-solid fuel （TMSR-SF） by using the developed program. Methods: The TMSR-SF reactor core model is built according to the designed parameters. Analyses of the impacts of different shutdown limits and shutdown signals on the withdrawal accident of the TMSR-SF uncontrolled rod cluster control assembly bank is carried out through the safety analysis models such as neutron kinetics model, fluid flow and heat transfer model, control model, etc., in the developed program. Results: In the withdrawal accident of the uncontrolled rod cluster control assembly bank of the TMSR-SF, the maximum temperature of the reactor fuel temperature is no more than 860°C. Conclusions:Computional results show that the higher the limit value of over power shutdown is, the greater the outlet temperature is； the longer the signal delay time is, the later the reactor shutdown happens, and the higher the core power and fuel temperature is. Compared with the limit value of shutdown, the outlet temperature reaches the limit later in time sequence, the maxumum fuel temperature is much lower than the fuel pebble melting temperature limit of 1600°C under no circumstances.

TMSR-SF, Safety analysis, Uncontrolled rod cluster control assembly bank withdrawal, Application development

ZUO Jiaxu, E-mail: zuojiaxu@chinansc.cn

TL99

10.11889/j.0253-3219.2016.hjs.39.100604

中国科学院战略性先导科技专项（No.XDA02050500）资助

靖剑平，男，1983年出生，2010年于哈尔滨工业大学获博士学位，研究领域为反应堆热工水力及安全分析

左嘉旭，E-mail: zuojiaxu@chinansc.cn

Supported by Strategic Priority Research Program of Chinese Academy of Sciences （No.XDA02050500）First author: JING Jianping, male, born in 1983, graduated from Harbin Institute of Technology with a doctoral degree in 2010, focusing on the reactor thermal hydraulic and safety analysis

2016-06-02，

2016-08-17