开发了一种用于卵石床高温气冷反应堆(pebble bed HTGR)的多尺度DLOFC源项分析框架

时间:2026年3月8日
来源:Annals of Nuclear Energy

编辑推荐:

高温气冷堆(HTR-PM)压力气冷堆强制冷却丧失(DLOFC)事故的放射性源项分析中,提出了一种多尺度耦合框架,整合中子物理、热水力、燃料性能及放射性迁移模型,通过精细化节点划分和统计燃料球操作历史,实现事故全过程中核岛内放射性核素分布的定量表征。结果显示,传统模型低估了氪-85(3.2×10⁻⁸)、氙-133(1.4×10⁻⁷)和铯-137(6.7×10⁻⁸)的放射性释放,新方法环境释放分数降低约一个数量级。

广告
   X   

成浩曹|陈俊毅|沈绍宁|梁静刚|李川|曹建柱
清华大学核能与新能源技术研究院,北京,100084,中国

摘要

高温气冷反应堆(HTGR)中的失压强制冷却(DLOFC)事故对人类健康和环境存在潜在的放射性危害。然而,现有的源项模型通常基于简化的假设,这在事故后果评估中引入了显著的不确定性。为了提高模拟的可靠性,本研究开发了一个多尺度DLOFC源项分析框架,该框架能够捕捉事故在时间和空间维度上的物理过程。该框架整合了一套高保真的、相互连接的模型,以定量描述所有关键事故阶段和反应堆区域中的放射性核素传输和分布。核心源项使用多物理场耦合模型进行模拟,而稳态条件下燃料颗粒的操作历史的统计集合定义了初始事故状态,从而能够对整个事故过程中的全堆释放进行分析。在减压之前,核岛内的放射性核素分布通过一维四区传输-沉积模型确定,该模型结合了主回路的详细节点化处理。此外,该框架还结合了FRG解吸模型和操作数据,以模拟不同事故阶段下放射性核素在堆芯、主回路、安全壳和环境中的迁移。以HTR-PM作为参考案例,进行了全面的DLOFC源项分析。结果表明,典型核素的环境释放分数分别为:长寿命惰性气体Kr-85为3.2×108,短寿命惰性气体Xe-133为1.4×107,长寿命金属核素Cs-137为6.7×108

引言

在应对全球气候变化的背景下,限制放射性释放是大规模部署核反应堆的关键要求(Moe, 2010)。然而,三里岛、切尔诺贝利和福岛的事故凸显了第二代和早期第三代反应堆相关的放射性释放的重大风险,从而削弱了公众对核能安全的信心。为了提高核安全,第四代反应堆(Locatelli等人,2013),特别是高温气冷反应堆(HTGR),由于其固有的安全特性、无熔化堆芯设计和被动安全系统,重新引起了全球的关注。HTGR在各国进行的测试中展示了在罕见事故场景(DBA1)中的稳健安全特性(Zhang等人,2024;Nakagawa等人,2004)。然而,某些极端事故场景(DBA2),如失压强制冷却(DLOFC)(Dou等人,2024;Zhipeng等人,2020),由于实际限制(Sieh和Bindra,2024),难以通过实验进行测试,可能带来安全风险。
在事故安全分析中,数值模拟在补充实验测试方面起着至关重要的作用。基于详细事故描述的准确模拟首先识别涉及的放射性源项,然后为这些源项提出适当的模型,最终提供可靠的放射性释放后果。DLOFC是HTGR可能发生的最严重事故。在DN65管道破裂导致的减压过程中,潜在的事故源项包括随冷却剂喷出的空气中的放射性物质、吸附在回路和氦气净化系统(HPS)表面的核素的解吸、沉积的石墨粉尘的重新悬浮以及由于堆芯加热而产生的额外释放。这些放射性物质通过泄漏口流入安全壳,并由于安全壳的减压而随后释放到环境中,造成放射性危害。
大量研究集中在DLOFC事件期间的瞬态热流体响应上(Zheng等人,2018;Husnayani等人,2021;Sato等人,2014;Guan等人,2022;Xu等人,2024),一些研究考察了来自堆芯的放射性释放(Berens等人,2024)或应用保守模型进行环境排放评估(Wang等人,2020;Ao Liu等人,2024)。然而,对核岛内放射性物质行为的多尺度分析仍然相对有限。Wagner使用了反应堆安全分析代码MELOCR(Wagner等人,2023),该代码整合了中子物理、热流体、化学和源项,分析了PBMR-400在DLOFC场景下的系统动态特性。MELOCR的强大集成分析能力使其能够模拟堆芯温度分布的变化、回路内的自然循环的建立以及由于腐蚀导致的石墨结构的退化。然而,在关注源项分析时,均质化燃料颗粒的建模忽略了燃料颗粒随机操作历史对加热阶段核素释放的影响。AVR的辐照后燃料性能测试结果不能直接代表新反应堆设计中使用的TRISO颗粒。此外,使用简化的化学吸附系数来评估空气中核素在回路表面的沉积引入了额外的不确定性。这些源项模型与实际条件之间的差异可能导致事故后果评估的不准确。
DLOFC事故期间的放射性释放可以分为两个阶段:早期减压释放和加热阶段的长期额外释放。前者源于反应堆稳态运行期间在主回路中积累的放射性物质因减压而喷射。后者源于由于堆芯加热而滞留在燃料中的裂变产物(FPs)的加速向外扩散。这两个源项模型的准确性对于模拟结果的可靠性至关重要。大量研究和项目专注于准确模拟稳态运行期间回路中核素和粉尘的传输行为,旨在评估放射性热点设备的维护安全性和辐射屏蔽的设计。例如,我们之前关于主回路源项特性的研究(Cao等人,2025),以及LOOP(Chen等人,2018)、FIST(Li,2019)、RADAX(Stassen,2006)、SPATRA(Moormann,1992)、GAMMA-FP(Yoo等人,2010)和PADLOC(Hudritsch和Smith,2018)等项目,都为这一研究领域做出了贡献。然而,很少有研究关注整个回路的详细节点离散化,并将所得的放射性分布应用于事故源项计算。对于核素从堆芯的扩散和释放,一些研究对平衡堆芯的释放进行了彻底分析,例如我们小组之前的随机颗粒全堆释放模型(Cao等人,2025a),以及STACY(Xhonneux和Allelein,2014)、FIPREX-GETTER(van der Merwe和Clifford,2008)和FIRCS(Raza等人,2025a)等代码。这些模型有潜力与完整的事故进展相结合,实现颗粒级的全堆额外释放模拟。然而,相关报告仅限于Miftasani的工作(Miftasani等人,2025)和van der Merwe的研究(van der Merwe和Clifford,2008)。Miftasani基于STACY程序对PELUIt堆芯进行了多物理场耦合模拟,整合了燃料性能、中子物理和堆芯热流体。不幸的是,PELUIt使用的是OTTO燃料循环,其操作特性与MEDUL循环反应堆(如高温气冷反应堆颗粒床模块HTR-PM Zhang等人,2018)示范电厂)有显著差异,因此不能代表所有颗粒床反应堆。尽管van der Merwe模拟了DLOFC场景下的颗粒床堆芯释放,但没有进行覆盖整个核岛范围的多尺度耦合模拟。
总之,以往的研究对于DLOFC场景采用了相对粗糙的源项模型,并缺乏对核岛内源项动态的多尺度、完全耦合分析。为了克服这些限制并提供更可靠的放射性释放预测,本研究提出了一个全面的多尺度DLOFC源项分析框架。
所提出的框架结合了中子物理、堆芯热流体、燃耗、核素扩散、热传导和燃料性能分析,构建了一个颗粒级的全堆释放模型,该模型与事故进展相结合。它解决了稳态和事故条件下的时间依赖性全堆释放率。此外,该框架采用了一维四区传输-沉积模型和集总参数模型来表示核岛内的核素行为。这些模型应用于精细节点化的主回路,能够在事故发生前计算核岛内的初始放射性分布。最后,该框架结合了FRG的解吸模型和操作数据,模拟稳态运行、减压时刻以及随后的加热阶段期间核芯、主回路、安全壳和环境之间的核素传输路径,从而能够计算特定位置和时间间隔的传输活性。本研究使用HTR-PM作为参考案例,对DLOFC场景进行了全面的源项分析,并展示了相应的放射性释放后果。

部分摘录

源项模型

如第1节所述,多尺度DLOFC源项分析框架将事故分为三个时间阶段——稳态、减压阶段和加热阶段——以及四个空间区域:堆芯、主回路、安全壳和环境。图1展示了DLOFC事故的这种二维划分,突出了框架考虑的每个阶段核岛内的潜在源项、传输路径和关键物理过程。

稳态全堆释放分析

在稳态运行期间,堆芯中FPs的释放是后续模型的主要输入,因此控制了源项分析的总体准确性。为了获得统计上稳健的结果,本研究在稳态建立阶段不断向堆芯中插入8000个随机颗粒。在平衡状态下,有超过6200个颗粒留在堆芯中——超过了Xhonneux的STACY模拟中使用的5000个样本和Raza采用的2000个样本

结论

本研究开发了一个用于DLOFC事故下源项分析的多尺度框架。通过将事故进展分解为时间和空间维度上的相互关联的模型,该框架能够详细量化所有关键区域和阶段的核素传输和分布。每个模型的主要特点总结如下:
  • 1.
    稳态全堆释放模型: 该模型整合了中子学、堆芯热流体、燃料性能等

利益冲突声明

作者声明他们没有已知的可能会影响本文报告工作的竞争性财务利益或个人关系。

致谢

本工作得到了北京自然科学基金(编号:1252021)的资助

生物通微信公众号
微信
新浪微博


生物通 版权所有