一、绪论

1.1 研究背景与意义

在全球能源转型加速与碳中和目标驱动下,核电作为稳定、低碳的基荷能源,其安全高效发展成为各国能源战略的核心议题。燃料元件作为反应堆 "能量心脏",在高温、高压、强辐照的极端环境中运行,其完整性直接决定反应堆运行安全与经济性。据国际原子能机构(IAEA)统计,全球在运核电机组中,约 35% 的非计划停堆与燃料相关问题直接相关,而燃料性能模型的精准度每提升 10%,可使燃料换料周期延长 3-6 个月,单机年经济效益增加数千万元。

燃料性能模型通过对燃料芯块肿胀、致密化、裂变气体释放及燃料 - 包壳相互作用(PCI/FCMI)等微观行为的定量化描述,建立起运行工况与燃料完整性的关联桥梁。在变功率运行、负荷跟踪等复杂工况下,FPM 能够提前预判燃料元件潜在失效风险,为反应堆运行优化、燃料设计改进及安全审评提供核心技术支撑。尤其在先进核电技术(如事故容错燃料、小型模块化反应堆)快速发展的背景下,FPM 的研发与应用已成为衡量核电技术水平的重要标志。

1.2 核心研究范畴界定

本报告聚焦核电领域燃料性能模型的全生命周期分析,核心研究范畴包括:

模型本质属性:以燃料元件微观物理化学行为为核心,融合热工、力学、材料科学的多物理场耦合模型体系;

关键行为维度:涵盖燃料芯块(肿胀、致密化、重结构、开裂)、裂变气体(形核、长大、释放)、界面作用(PCI/FCMI、间隙演化)三大核心行为域;

应用场景边界:覆盖正常运行、瞬态工况、事故工况全谱系场景,重点关注变功率运行下的完整性评估;

技术体系构成:包括基础理论模型、工程化计算程序、实验验证体系及行业应用规范的完整技术链。

二、燃料性能模型的历史背景与演进脉络

2.1 萌芽期(1950s-1960s):经验模型的初步探索

2.1.1 技术背景与驱动因素

这一时期正值全球核电技术的起步阶段,反应堆类型以实验性石墨堆、压水堆为主,燃料元件设计简单(多为金属铀燃料),运行燃耗较低(通常低于 10 GWd/MTU)。由于计算机技术限制(早期真空管计算机计算能力仅为千次 / 秒级)和实验数据匮乏,燃料性能研究主要依赖堆内辐照试验与事后解剖分析。

1957 年美国 Shippingport 压水堆的建成运行,标志着商用核电的开端,也催生了对燃料性能预测的初步需求。此阶段的核心挑战是解决燃料元件的基本完整性问题,如早期金属铀燃料的肿胀破裂与包壳腐蚀问题。

2.1.2 代表性模型与技术特征

这一时期的燃料性能模型以经验公式为主,聚焦单一物理现象的简化描述:

肿胀模型:基于金属铀燃料实验数据,建立肿胀量与燃耗的线性经验公式(如 Uranium Swelling Model, USM),未考虑温度梯度影响;

热传导模型:采用恒定热导率假设,简化燃料温度场计算,典型代表为美国橡树岭国家实验室(ORNL)提出的一维热传导解析解;

裂变气体释放模型:基于压力平衡理论的半经验模型,释放率计算仅考虑温度单一变量。

1963 年,美国西屋公司开发的首个燃料性能计算程序 FUEL-1,标志着模型从手工计算向程序化的转变,该程序仅包含 4 个核心经验公式,可计算燃料中心温度与包壳应力的基本参数。

2.1.3 历史意义与局限

萌芽期模型实现了燃料性能从 "无法预测" 到 "初步量化" 的跨越,为早期反应堆燃料设计提供了基本工具。但受限于技术条件,存在三大核心局限:缺乏微观机理支撑、未考虑多物理场耦合效应、预测精度低(误差通常超过 30%),难以适应燃耗提升与工况复杂化的需求。

2.2 发展期(1970s-1990s):机理模型的体系构建

2.2.1 技术突破与需求升级

1970s 后,商用核电进入快速发展期,压水堆、沸水堆成为主流堆型,UO₂陶瓷燃料取代金属铀燃料成为主流。随着计算机技术的发展(1970s 晶体管计算机计算能力达百万次 / 秒级,1990s 奔腾处理器达亿次 / 秒级),多物理场耦合计算成为可能。

三里岛事故(1979 年)对燃料性能模型的发展产生了深远影响,推动模型从正常运行评估向事故工况扩展。此阶段燃料燃耗持续提升(从 15 GWd/MTU 增至 45 GWd/MTU),核心需求转变为提高燃料经济性与事故容错能力。

2.2.2 模型体系的形成与典型程序

这一时期形成了以机理模型为核心、多物理场初步耦合的模型体系,诞生了一批影响深远的工程化程序:

程序名称

开发机构

诞生时间

技术特征

应用范围

FRAPCON

美国INL

1974 年

一维柱状几何,热-力学耦合,模块化结构

PWR/BWR 正常运行

TRANSURANUS

欧洲联合团队

1982 年

多堆型适配,考虑燃料重结构效应

压水堆、快堆

FEMAXI

日本JAERI

1977 年

侧重沸水堆特性,细化空泡模型

BWR 燃料分析

ANSYS-Fuel

美国ANSYS 公司

1990 年

有限元框架,强力学分析能力

燃料力学行为

1980 年代,裂变气体释放的 "两阶段模型"(晶粒内扩散 - 晶界释放)成为主流机理模型,由 Lassmann 等人提出的该模型首次揭示了裂变气体释放的微观机制,使释放率预测误差降至 15% 以内。

2.2.3 关键技术进展与行业影响

此阶段实现了三大技术突破:一是建立了 UO₂燃料的热 - 力学耦合基础模型,包括 Mazur 模型(1973)描述的芯块肿胀机理;二是开发了燃料 - 包壳相互作用(FCMI)的力学分析框架,引入接触压力计算模型;三是形成了初步的实验验证体系,如美国 HALDEN 反应堆的燃料辐照考验平台。

这些进展支撑了燃料燃耗从 30 GWd/MTU 向 50 GWd/MTU 的跨越,使核电经济性显著提升。同时,模型开始纳入监管体系,美国 NRC 将 FRAPCON 结果作为燃料许可的重要依据。

2.3 精细化期(2000s-2010s):多物理场耦合与程序升级

2.3.4 技术驱动与挑战

21 世纪以来,核电面临 "高燃耗、长周期、高安全" 的三重需求,燃料燃耗目标突破 60 GWd/MTU,同时福岛核事故(2011 年)推动模型向事故容错能力强化方向发展。计算机技术的飞跃(高性能计算集群峰值达千万亿次 / 秒)为多物理场精细化耦合提供了支撑。

高燃耗带来的新问题(如高燃耗结构 HBS 形成、包壳氢脆)对模型提出了新挑战,传统一维模型已难以满足精度需求,二维 / 三维耦合模型成为发展方向。

2.3.5 模型技术的核心突破

这一时期的模型实现了从 "一维到多维"、"单场到多场"、"宏观到宏微观结合" 的三大跨越:

多维度扩展:从传统一维柱状模型(FRAPCON 系列)发展到二维轴对称模型(如 TRANSURANUS 的 2D 模块)和三维有限元模型(如 BISON 的 3D 计算);

多物理场深化:实现热 - 力 - 辐照 - 化学多场强耦合,如引入氧势场对裂变气体释放的影响模型;

微观机理融入:将分子动力学模拟结果引入宏观模型,如晶粒尺度的裂变气体扩散模型。

美国爱达荷国家实验室(INL)2009 年发布的 BISON 程序,基于 MOOSE 开源有限元平台,首次实现了燃料性能的全三维多物理场耦合计算,标志着精细化模型时代的到来。

2.3.6 国际合作与标准体系构建

OECD/NEA 和 IAEA 主导的国际合作项目推动了模型的标准化与比对验证:

•1990 年代启动的 FUMEX 系列项目(FUMEX-I 至 FUMEX-III),组织全球 20 余个机构开展燃料性能模型国际比对;

•IAEA 的 Coordinated Research Project(CRP)项目,建立了涵盖 PWR、BWR、VVER 的统一验证数据库;

•美国 NRC 发布 NUREG/CR-6534 报告,建立了燃料性能模型的验证与确认(V&V)标准流程。

2.4 创新期(2010s 至今):先进模型与智能融合

2.4.1 技术变革与需求升级

福岛核事故(2011 年)后,事故容错燃料(ATF)成为研究热点,推动模型向极端工况适应性发展。同时,第四代核电技术(快堆、高温气冷堆)与小型模块化反应堆(SMRs)的研发,对模型提出了多燃料类型(MOX、金属燃料、TRISO)与灵活运行工况的适配需求。

燃耗目标持续提升至 75-90 GWd/MTU,高燃耗结构(HBS)带来的材料性能劣化问题成为模型核心挑战。人工智能与大数据技术的兴起,为模型优化与不确定性量化提供了新途径。

2.4.2 前沿模型技术特征

ATF 专用模型:针对 SiC 包壳、Cr 涂层、U₃Si₂芯块等新型材料,开发专用物性模型与界面行为模型,如西安交通大学基于 MOOSE 平台开发的 BEEs 程序,实现了多层 SiC 包壳的精细化建模;

多尺度耦合模型:建立从纳米级(裂变碎片轨迹)、微观级(晶粒演化)到宏观级(元件响应)的跨尺度模型,如 MIT 开发的 GRSIS 模型,实现金属燃料气泡演化的多尺度描述;

AI 融合模型:引入物理知情神经网络(PINNs)与代理模型,如 BISON 程序集成的多项式混沌展开不确定性量化方法,将计算效率提升 10 倍以上;

数字化孪生集成:开发与反应堆数字孪生衔接的实时预测模型,支持运行状态动态评估与寿命预测。

2.4.3 全球技术格局演变

当前形成了三大技术阵营:

美国阵营:以 INL 的 BISON 程序为核心,侧重开源化、多堆型适配与 AI 融合;

欧洲阵营:以 TRANSURANUS 为代表,强调多燃料类型(尤其是 MOX 燃料)的成熟验证;

亚洲阵营:中国(BEEs、COPERNIC)、日本(FEMAXI-8)侧重自主化与先进堆适配,韩国(PyFRAP)强化 AI 技术应用。

三、燃料性能模型的具体实现方式

3.1 模型的理论基础与核心构成

燃料性能模型是融合热传导、固体力学、材料科学、核化学的多学科理论体系,其核心构成包括三大模块:

3.1.1 燃料芯块行为模型

燃料芯块作为能量产生的核心载体,其行为模型涵盖四个关键过程:

1.热物性模型核心参数包括热导率、比定压热容、热膨胀系数,均需考虑温度、燃耗、辐照的耦合影响:

•热导率模型:采用 Maxwell-Eucken 方程描述孔隙率影响,引入辐照劣化因子,典型公式为:k(T,B)=k_0(_T)⋅f(ϕ)⋅g(P)其中 k₀为未辐照热导率,f (φ) 为辐照修正因子,g (P) 为孔隙率修正因子;

•熔点模型:考虑燃耗导致的熔点降低,采用经验公式:Tₘ(B) = 3110 - 32×(B/10)(单位:K,B 为燃耗 GWd/MTU)。

2.肿胀与致密化模型

•致密化模型:描述低燃耗阶段(<10 GWd/MTU)空位扩散导致的密度增加,采用扩散方程描述空位迁移;

•肿胀模型:分为固体裂变产物肿胀与气态裂变产物肿胀,前者与燃耗呈线性关系(每 1% 原子燃耗肿胀 0.3%),后者基于气泡形核长大理论,如 Forsberg-Massih 模型将气泡分为封闭气泡与开放气泡分别计算。

3.重结构与开裂模型

•重结构模型:基于物质迁移理论,描述高温梯度下的气孔迁移与柱状晶形成,采用热蒸发 - 冷凝方程模拟;

•开裂模型:考虑热应力与机械应力耦合作用,采用断裂力学准则(如最大主应力准则)判断裂纹萌生与扩展。

3.1.2 裂变气体行为模型

裂变气体(主要为 Kr、Xe)占裂变产物的 15% 左右,其释放行为直接影响燃料棒内压与包壳应力,模型涵盖三个关键阶段:

1.气泡形核模型基于经典形核理论,考虑辐照缺陷(空位、位错)的催化作用,形核率 N 与过饱和度 S 的关系为:N\=C⋅exp(−Δ_G_∗/kT)其中 ΔG * 为临界形核自由能,C 为常数。

2.气泡长大模型考虑气体原子扩散与气泡合并两种机制,高温区(>1600K)以扩散主导,低温区以合并主导。对于金属燃料,GRSIS 模型将气泡分为三类(小封闭气泡、大封闭气泡、开放气泡)分别计算长大速率。

3.释放模型采用 "晶粒内扩散 - 晶界传输 - 释放" 的两阶段模型,释放率 F 的典型计算公式为:F\=1−exp(−λt)其中 λ 为释放系数,与温度、燃耗正相关。当温度超过 1973K 时,释放率可达 70%-95%。

3.1.3 燃料 - 包壳相互作用模型

燃料 - 包壳相互作用(FCMI)与 pellet-cladding interaction(PCI)是影响燃料完整性的关键界面行为:

1.间隙演化模型描述燃料与包壳间间隙的闭合过程,考虑燃料肿胀、包壳蠕变与热膨胀的耦合作用,间隙宽度δ(t) 的计算方程为:dtdδ\=Δ_Lf_−Δ_Lc_−δ⋅_ϵcreep_其中ΔL_f 为燃料膨胀量,ΔL_c 为包壳膨胀量,ε_creep 为包壳蠕变率。

2.PCI 模型分为机械相互作用与化学相互作用:

•机械 PCI:燃料芯块膨胀导致的包壳弹性变形与塑性屈服,采用弹塑性力学模型计算接触压力;

•化学 PCI:裂变产物(Cs、Te)与包壳的应力腐蚀开裂,引入腐蚀动力学与断裂力学耦合模型。

3.包壳行为模型涵盖包壳热蠕变(采用 Power-law 模型)、辐照蠕变(Hoppe 模型)、氧化(抛物线速率定律)等行为,对于 Zr 合金包壳,氧化层厚度 δ_ox 的计算为:δox_2=_K⋅exp(−Q/RT)⋅_t_其中 K 为常数,Q 为活化能,R 为气体常数。

3.2 典型工程化计算程序解析

工程化计算程序是燃料性能模型的具体载体,目前全球主流程序超过 20 种,其中 FRAPCON、TRANSURANUS、BISON 构成三大标杆体系。

3.2.1 美国 FRAP 系列程序

由美国 INL 开发的 FRAP(Fuel Rod Analysis Program)系列是全球应用最广泛的程序族,包括 FRAPCON(稳态)、FRAPTRAN(瞬态)、FRAPTOP(事故)三个核心模块:

技术特征

•几何建模:采用一维柱状几何,将燃料棒分为多个轴向切片,每个切片含燃料、间隙、包壳三层;

•核心模型:

○热传导:考虑径向温度梯度与辐照热导率劣化;

○肿胀模型:采用 Cumberland-Glassman 模型,区分致密化与肿胀阶段;

○FCMI:基于弹性 - 塑性 - 蠕变耦合模型计算接触压力;

•数据基础:基于 MATPRO 数据库,包含 2000 + 组燃料与包壳物性数据。

应用场景:主要用于 PWR/BWR 燃料正常运行与瞬态工况分析,NRC 将其作为燃料许可的基准程序。FRAPCON-4.0 版本可支持最高 62 GWd/MTU 的燃耗计算,预测误差控制在 10% 以内。

3.2.2 欧洲 TRANSURANUS 程序

由欧洲联合研究团队(含比利时 SCK・CEN、德国 KIT 等)开发,是目前功能最全面的燃料性能程序之一:

技术特征

•多维度建模:支持 0D(整体平均)、1D(径向)、2D(径向 - 轴向)三种建模模式,2D 模式可模拟芯块开裂与重结构的轴向分布;

•多燃料适配:可处理 UO₂、MOX、金属燃料、TRISO 等多种燃料类型,在 MOX 燃料分析方面积累了丰富验证数据;

•模块化设计:包含热工、力学、裂变气体、材料劣化等 12 个独立模块,支持用户自定义模型扩展。

应用场景:覆盖 PWR、BWR、VVER、快堆等多种堆型,参与了 OECD/NEA 的 PRIMO MOX 燃料基准测试,在欧洲核电运营商中应用率超过 80%。

3.2.3 美国 BISON 程序

作为新一代开源燃料性能程序,BISON 基于 MOOSE 多物理场有限元平台开发,代表了当前技术发展方向:

技术特征

•全三维建模:采用有限元方法实现三维几何建模,可精确模拟燃料开裂与包壳变形的空间分布;

•多物理场强耦合:实现热 - 力 - 辐照 - 化学场的全耦合求解,支持非线性材料行为与接触非线性;

•开源扩展性:代码完全开源(GitHub 托管),支持用户添加自定义模型,已集成 ATF 专用模块;

•AI 融合能力:集成代理模型与不确定性量化工具,计算效率较传统程序提升 5-10 倍。

应用场景:重点用于先进燃料(ATF、金属燃料)与先进堆(快堆、SMRs)的性能分析,美国 TerraPower 的 Natrium 快堆项目已将其作为核心燃料分析工具。

3.3 不同堆型的模型适配与技术差异

燃料性能模型需根据堆型特性(冷却剂类型、中子能谱、运行工况)进行针对性适配,主要堆型的技术差异如下:

3.3.1 压水堆(PWR)模型适配

PWR 采用水作为冷却剂,运行压力约 15.5 MPa,燃料中心温度通常低于 2400 K,模型适配重点包括:

•热工模型:考虑冷却剂沸腾对传热的影响,引入 DNB(偏离泡核沸腾)校正因子;

•FCMI 模型:重点模拟稳态运行下的间隙闭合过程,采用弹性 - 塑性蠕变模型计算包壳应力;

•高燃耗适配:针对 75 GWd/MTU 以上燃耗,强化 HBS 结构的裂变气体释放模型,如 FRAPCON-4.0 新增的 HBS 修正模块。

典型应用案例:美国西屋 AP1000 机组采用 FRAPCON-4.0 与 BISON 联合分析,将燃料换料周期从 18 个月延长至 24 个月。

3.3.2 沸水堆(BWR)模型适配

BWR 运行压力较低(约 7 MPa),燃料棒直接与饱和蒸汽接触,模型需重点考虑:

•水侧腐蚀模型:强化包壳腐蚀与氢吸收计算,采用修正的氧化动力学模型;

•空泡效应:考虑堆芯空泡率对中子通量与燃料温度的影响,建立空泡 - 温度耦合模型;

•功率瞬变:针对 BWR 负荷跟踪特性,优化瞬态热传导与 PCI 模型,如 FEMAXI-8 的瞬态响应模块。

3.3.3 快堆(Fast Reactor)模型适配

快堆采用液态金属(钠、铅铋)冷却,中子能谱硬,燃料以金属合金(U-Pu-Zr)为主,模型适配难点包括:

•金属燃料肿胀:采用 GRSIS 多尺度模型,模拟气泡形核、长大与连通过程,预测肿胀率误差小于 3%;

•包壳肿胀:针对 ODS 钢等先进材料,开发空洞肿胀模型,考虑快中子辐照(>100 dpa)的影响;

•在线换料:建立燃料轴向燃耗梯度模型,支持不停堆换料的性能评估。

俄罗斯 BN-800 快堆采用自主开发的 PERMA 程序,成功实现 U-Pu-Zr 燃料 60 GWd/MTU 燃耗的安全运行。

3.3.4 高温气冷堆(HTGR)模型适配

HTGR 采用 TRISO 包覆颗粒燃料,氦气冷却,运行温度高达 1000℃以上,模型核心是颗粒破损概率预测:

•多层包覆层模型:模拟 Buffer 层、致密热解炭层的辐照收缩与蠕变,采用修正的 PANAMA 方法计算应力分布;

•破损概率模型:考虑裂变气体内压、热应力、辐照损伤的耦合作用,建立统计性破损概率模型;

•燃料球流动:针对 pebble-bed 堆型,开发燃料球流动与燃耗分布耦合模型。

中国 HTR-PM 项目采用自主 FCAS 程序,TRISO 颗粒破损率预测值低于 1×10⁻⁵,满足安全要求。

3.4 模型的实验验证体系与标准化流程

燃料性能模型的可靠性依赖于完善的实验验证体系,国际通行 "基础实验 - 堆内考验 - 基准比对" 三级验证流程。

3.4.1 基础实验体系

针对单一物理行为开展的实验室模拟实验,主要包括:

热物性实验:采用激光闪射法测量燃料热导率,差示扫描量热法(DSC)测量比热容;

力学性能实验:通过拉伸试验、蠕变试验获取材料应力 - 应变曲线,辐照后力学性能采用辐照后检验(PIE)获取;

裂变气体实验:采用热重分析(TGA)测量释放率,透射电镜(TEM)观察气泡微观分布。

美国 INL 的燃料性能实验平台可开展温度高达 2800 K、压力 20 MPa 的模拟实验,为模型提供基础数据支撑。

3.4.2 堆内辐照考验

在实际反应堆环境中开展的燃料性能测试,是验证模型的核心环节:

仪器化燃料棒:内置热电偶、压力传感器、位移传感器,实时监测运行参数;

典型考验项目:包括高燃耗考验(如美国的 LTA 程序,燃耗达 75 GWd/MTU)、瞬态考验(如法国 CABRI 反应堆的 RIA 测试)、ATF 专用考验(如美国 ATF-1 项目);

国际共享数据库:IAEA 的 IFPE 数据库包含全球 30 余个国家的 1000 + 组辐照考验数据。

3.4.3 国际基准比对

通过国际合作开展模型预测与实验结果的比对,验证模型的通用性与精度:

OECD/NEA 基准测试:如 PRIMO MOX 燃料基准,组织 15 个国家 20 个程序开展比对,裂变气体释放预测误差最优达 5%;

IAEA CRP 项目:FUMEX-III 项目涵盖 PWR、BWR、VVER 三种堆型的 12 个基准题,重点验证高燃耗模型;

NRC 标准考题:针对 LOCA、RIA 事故工况,建立标准化验证考题,要求模型预测误差小于 15%。

3.4.4 标准化流程与质量保证

国际标准化组织(ISO)与 IAEA 建立了燃料性能模型的质量保证体系:

V&V 流程:遵循 ISO 17025 与 NUREG/CR-6850 标准,包括模型验证(与实验比对)、确认(与工程需求比对)、不确定性量化三个环节;

参数管理:采用结构化数据库管理模型输入参数,建立参数溯源与不确定性分析机制;

程序认证:NRC 对用于安全审评的程序实施认证制度,BISON 程序于 2023 年通过认证,成为首个获认证的开源程序。

四、燃料性能模型的主要争议点与各方立场

4.1 技术层面的核心争议

4.1.1 高燃耗下裂变气体释放模型的适用性争议

随着燃耗从 60 GWd/MTU 提升至 75-90 GWd/MTU,燃料芯块形成高燃耗结构(HBS),其特征为亚微米级晶粒与贯通性孔隙,导致裂变气体释放率显著升高(从 <10% 增至> 30%)。现有模型在描述这一现象时存在显著分歧:

争议焦点

•释放机制之争:传统模型基于 "晶粒内扩散 - 晶界释放" 两阶段理论,但 HBS 结构下是否出现 "体积扩散主导" 的新机制存在争议;

•模型参数之争:FRAPCON 采用的经验修正因子(HBS 因子 f=1.5-2.0)与 TRANSURANUS 的机理模型预测结果偏差可达 20%-30%;

•验证数据之争:高燃耗辐照数据匮乏,不同反应堆的 HBS 形成条件差异导致数据可比性差。

典型案例:2018 年 OECD/NEA 组织的高燃耗基准测试中,12 个参与程序对同一工况的裂变气体释放预测值范围为 18%-42%,离散度远超低燃耗工况(<10%)。

4.1.2 PCI/FCMI 模型的机理与预测精度争议

燃料 - 包壳相互作用是燃料失效的主要诱因(约占燃料失效的 40%),但其复杂的物理化学耦合机制导致模型预测精度有限:

争议焦点

•化学 PCI 机理之争:对于 Cs-Te 化合物导致的应力腐蚀开裂,存在 "液态金属脆化" 与 "氧化膜破裂" 两种主导机理假设,分别对应不同的模型构建逻辑;

•瞬态 PCI 预测之争:变功率运行下,燃料芯块快速膨胀导致的包壳塑性变形,传统准静态模型预测误差可达 40% 以上;

•界面参数之争:燃料 - 包壳接触换热系数的取值存在分歧,FRAPCON 采用经验公式,BISON 采用接触压力耦合模型,结果差异显著。

行业影响:为保证安全,核电运营商通常采用保守的功率变化速率限制(如≤5% Pn/min),导致机组负荷跟踪能力受限,经济性降低。

4.1.3 事故工况下模型的有效性争议

事故工况(如 LOCA、RIA)下燃料行为的极端性与复杂性,使模型面临 "预测可靠性" 的核心挑战:

争议焦点

•高温材料物性之争:UO₂在 1800-3100 K 区间的热物性(尤其是熔化前后的热导率)数据稀缺,不同模型采用的插值方法导致温度预测偏差达 300 K 以上;

•燃料熔化与重定位之争:LOCA 事故下燃料熔池形成与迁移的模拟,存在 "宏观流动模型" 与 "颗粒流模型" 的技术路线分歧;

•验证局限性之争:事故工况无法通过实际反应堆验证,仅能依赖模拟实验(如 CORA、PHEBUS FP),存在 "实验尺度效应" 导致的模型外推误差。

监管关切:美国 NRC 在 2022 年的安全审评报告中指出,现有模型对严重事故下燃料失效的预测不确定性高达 100%-200%,需进一步强化验证。

4.1.4 先进燃料材料的模型描述争议

事故容错燃料(ATF)等新型材料的涌现,使传统模型面临适应性挑战:

争议焦点

•SiC 包壳模型之争:SiC/SiC 复合材料的辐照蠕变与断裂行为,现有模型(基于陶瓷材料经验)预测误差达 50% 以上;

•金属燃料肿胀之争:UMo 燃料的辐照肿胀机制与 UO₂存在本质差异,GRSIS 模型与 FEAST-Metal 模型的预测结果偏差显著;

•界面行为之争:ATF 燃料的涂层 - 基体界面(如 Cr 涂层 - Zr 基体)结合强度退化模型,缺乏足够的实验数据支撑。

发展瓶颈:ATF 燃料的模型不确定性导致其商业化进程延迟,美国计划将 ATF 大规模应用时间从 2025 年推迟至 2030 年。

4.2 利益相关方的立场差异

不同利益相关方基于自身目标与诉求,对燃料性能模型的技术路线与应用标准存在显著立场差异,形成了 "安全 - 经济 - 技术" 的三角博弈关系。

4.2.1 监管机构:以安全保守为核心原则

以美国 NRC、法国 ASN、中国国家核安全局为代表的监管机构,其核心目标是保障反应堆运行安全,立场具有显著的保守性:

关键立场

•模型验证要求:强调 "实验数据充分性",要求模型在全工况范围内均有实验验证,对高燃耗、ATF 等新型工况持谨慎态度;

•不确定性控制:要求模型不确定性量化结果纳入安全裕量计算,通常采用 2 倍以上的安全系数;

•程序认证标准:实施严格的程序认证制度,优先认可经长期工程验证的成熟程序(如 FRAPCON),对新型开源程序(如 BISON)认证周期长达 5-8 年;

•事故工况侧重:重点关注模型在极端事故下的预测能力,将 LOCA/RIA 工况的燃料失效预测作为审评核心。

典型政策:美国 NRC 在 NUREG-0933 报告中明确,高燃耗燃料的失效阈值从 170 cal/g 下调至 30 cal/g,显著提高了安全裕量。

4.2.2 核电运营商:以经济性与实用性为导向

以 EDF、Exelon、中核为代表的核电运营商,追求燃料经济性最大化与运行灵活性,立场更注重模型的工程实用性:

关键立场

•高燃耗支持:积极推动模型向更高燃耗(>75 GWd/MTU)适配,要求模型降低保守性,释放安全裕量;

•负荷跟踪需求:强调模型对变功率运行的精准预测,以支持机组参与电网调峰(功率变化速率≥10% Pn/min);

•计算效率优先:倾向于使用计算速度快的成熟程序(如 FRAPCON),对高精度但高耗时的三维模型持谨慎接受态度;

•成本控制导向:要求模型优化燃料设计参数(如芯块密度、间隙尺寸),以降低燃料制造成本与换料频率。

典型实践:法国 EDF 通过优化 TRANSURANUS 模型参数,将压水堆燃料换料周期从 12 个月延长至 18 个月,年运营成本降低约 15%。

4.2.3 燃料供应商:以技术创新与市场竞争为核心

以西屋、法马通、中核燃料为代表的燃料供应商,致力于通过技术创新提升市场竞争力,立场侧重模型的技术先进性:

关键立场

•先进燃料研发:大力投入 ATF、MOX 等新型燃料的模型开发,以形成技术壁垒;

•精细化设计:推动三维多物理场模型应用,以优化燃料设计(如环形芯块、功能梯度材料);

•数据独占性:倾向于建立自有数据库与模型体系,对开源模型的参与度较低;

•国际标准参与:积极参与 OECD/NEA、IAEA 的模型比对与标准制定,以主导技术话语权。

典型案例:西屋公司开发的 AP1000 先进燃料组件,基于自主改进的 FRAPCON 模型,实现了 72 GWd/MTU 的高燃耗设计,占据全球 25% 的燃料市场份额。

4.2.4 研究机构:以技术突破与科学严谨为目标

以 INL、KIT、西安交通大学为代表的研究机构,追求模型的科学严谨性与技术前沿性,立场相对中立但更具前瞻性:

关键立场

•机理研究优先:注重微观机理模型的开发,推动多尺度、多物理场耦合模型创新;

•开源合作倡导:积极开发开源程序(如 BISON、BEEs),促进全球技术合作与数据共享;

•前沿技术融合:探索 AI、数字孪生等新技术与传统模型的融合,提升模型性能;

•基础数据积累:投入大量资源开展基础实验,建立高精度数据库(如 INL 的 MATPRO 数据库)。

典型成果:西安交通大学基于 MOOSE 平台开发的 BEEs 程序,解决了传统程序无法模拟多层 SiC 包壳的技术瓶颈,为 ATF 研发提供了关键工具。

4.3 国际标准与规范的差异争议

燃料性能模型的应用依赖于统一的标准与规范,但由于技术路线、安全文化、监管体系的差异,国际上尚未形成全球统一的标准体系,存在显著的区域差异。

4.3.1 中美欧三大区域标准体系对比

对比维度

美国体系(NRC)

欧洲体系(EURATOM)

中国体系(国家能源局)

核心标准

NUREG 系列报告

EURATOM 技术导则

NB/T 系列标准

基准程序

FRAPCON/FRAPTRAN

TRANSURANUS

自主程序+ FRAPCON

燃耗限值

75 GWd/MTU(许可上限)

70 GWd/MTU(推荐值)

65 GWd/MTU(当前限值)

验证要求

需通过NRC 程序认证

需参与OECD/NEA 比对

需通过中核集团验证

不确定性处理

采用概率安全评估(PSA)

采用确定性安全裕量

混合概率- 确定性方法

ATF 适配

2023 年发布 ATF 审评指南

2024 年启动 ATF 标准制定

2025 年计划发布 ATF 标准

4.3.2 关键争议领域

1.模型认可标准之争美国采用 "程序认证 + 个案审评" 模式,对模型的技术细节与验证数据要求极为严格;欧洲采用 "通用导则 + 自主验证" 模式,给予运营商更大的灵活性;中国采用 "引进消化 + 自主创新" 模式,在参考国际标准的同时强化自主程序的应用。

2.数据共享与知识产权之争美国 NRC 要求模型输入数据需公开可追溯,而欧洲燃料供应商倾向于保护核心数据的知识产权,导致部分关键模型参数(如 MOX 燃料的裂变气体释放系数)无法共享。中国在自主数据库建设过程中,面临 "国际数据封锁 + 国内数据积累不足" 的双重挑战。

3.先进堆型标准空白之争对于快堆、高温气冷堆等先进堆型,国际上尚无成熟的模型标准体系。美国依托 BISON 程序主导快堆模型标准制定,欧洲在高温气冷堆 TRISO 燃料模型方面占据优势,中国则在小型模块化反应堆(SMRs)模型标准制定中积极发声。

4.监管互认障碍之争由于标准差异,不同国家对燃料性能模型的审评结果难以互认,导致燃料组件的跨境取证成本高昂。例如,中国自主燃料组件进入欧洲市场,需同时满足 EURATOM 导则与法国 ASN 的额外要求,认证周期长达 3-5 年。

五、燃料性能模型的未来发展方向

5.1 技术发展趋势:多维度创新驱动模型升级

5.1.1 多尺度多物理场耦合的深度融合

未来模型将实现从 "宏观唯象" 向 "宏微观耦合" 的根本性转变,构建涵盖量子力学、分子动力学、连续介质力学的跨尺度模型体系:

微观机理嵌入:通过分子动力学模拟获取裂变气体扩散系数、材料界面结合能等基础参数,嵌入宏观模型,如采用第一性原理计算 SiC 包壳的辐照损伤演化;

介观模型衔接:开发晶粒尺度模型,描述 HBS 结构形成与裂变气体迁移的介观过程,建立介观 - 宏观耦合界面;

多物理场强耦合:实现热 - 力 - 辐照 - 化学 - 中子场的全耦合求解,消除传统 "场分离" 求解带来的误差,如 BISON 程序正在开发的中子通量 - 热导率 - 肿胀率强耦合模块。

技术目标:到 2030 年,实现高燃耗工况下裂变气体释放预测误差≤10%,FCMI 应力预测误差≤15%。

5.1.2 人工智能与传统模型的深度融合

AI 技术将成为燃料性能模型优化的核心驱动力,形成 "物理机理 + 数据驱动" 的混合模型体系:

物理知情神经网络(PINNs):将热传导方程、力学平衡方程等物理约束嵌入神经网络,解决小数据量下的模型泛化问题,如用于 ATF 材料的未知物性预测;

代理模型加速:针对三维精细模型计算耗时问题,开发高精度代理模型,将计算效率提升 100-1000 倍,支撑实时在线预测;

不确定性量化(UQ):采用深度学习与蒙特卡洛模拟结合的方法,实现模型输入 - 输出不确定性的高效量化,如 BISON 程序集成的多项式混沌展开 UQ 模块;

智能参数校准:利用强化学习算法自动优化模型参数,实现模型与实验数据的自适应匹配。

应用场景:SMRs 的实时负荷跟踪、燃料数字孪生的在线评估、ATF 燃料的快速设计优化。

5.1.3 面向先进燃料与堆型的模型创新

随着 ATF、金属燃料、TRISO 燃料等新型燃料与快堆、SMRs、高温气冷堆等先进堆型的发展,模型需针对性创新:

ATF 专用模型:开发 SiC/SiC 复合材料的辐照蠕变、断裂模型,Cr 涂层的剥落与腐蚀模型,U₃Si₂芯块的肿胀与裂变气体释放模型;

金属燃料模型:优化 GRSIS 模型,实现 U-Pu-Zr 合金的肿胀、重熔与裂变气体释放的精准预测,支撑快堆闭式循环发展;

TRISO 燃料模型:建立考虑辐照损伤累积的多层包覆层应力演化模型,将破损概率预测精度提升至 1×10⁻⁶量级;

SMRs 适配模型:开发灵活功率运行下的瞬态响应模型,支持负荷跟踪与频繁启停工况的性能评估。

5.1.4 数字化孪生与模型的深度集成

燃料性能模型将成为反应堆数字孪生的核心组成部分,实现从 "离线分析" 向 "在线智能决策" 的转变:

实时数据接口:开发与堆芯监测系统(如堆外中子探测器、燃料组件温度监测)的数据接口,实现模型输入参数的动态更新;

状态感知与预测:基于实时运行数据与模型预测结果,实现燃料状态的动态感知与剩余寿命预测;

闭环优化决策:构建 "监测 - 预测 - 优化 - 控制" 的闭环系统,自动生成最优运行策略(如功率调整、换料时机);

虚拟试验平台:基于数字孪生构建虚拟燃料试验平台,减少物理实验成本,加速新型燃料研发。

5.2 应用场景拓展:从安全评估到全生命周期优化

5.2.1 燃料设计优化的智能化

传统燃料设计采用 "经验试错" 模式,周期长达 3-5 年,未来将基于模型实现 "精准设计":

多目标优化:以燃耗最大化、失效概率最小化、成本最低化为目标,采用多目标遗传算法优化燃料参数(芯块密度、包壳厚度、间隙尺寸);

虚拟验证:通过模型开展大规模虚拟试验,减少堆内考验数量,将燃料研发周期缩短至 1-2 年;

定制化设计:针对不同反应堆的运行特性(如调峰堆、基荷堆),开发定制化燃料设计方案,如为 SMRs 设计耐频繁启停的燃料组件。

案例展望:西屋公司计划基于 AI 融合模型,实现 ATF 燃料的全参数智能化设计,研发周期从 4 年缩短至 18 个月。

5.2.2 运行优化的精细化

模型将深度融入核电运行管理,实现运行策略的动态优化:

负荷跟踪优化:基于实时模型预测,动态调整功率变化速率,在避免 PCI 失效的同时最大化调峰能力;

换料策略优化:结合堆芯物理与燃料性能模型,优化换料方案,实现堆芯功率分布均匀化与燃料利用率最大化;

故障诊断与预警:基于模型预测结果与监测数据,实现燃料失效的早期诊断与预警,如裂变气体泄漏的定位与严重程度评估;

寿期延长评估:针对超设计寿期运行的反应堆,开展燃料性能的长期评估,制定延寿运行策略。

5.2.3 燃料循环的全链条支撑

模型将覆盖燃料 "设计 - 制造 - 运行 - 退役" 全生命周期,支撑闭合燃料循环发展:

制造过程优化:开发燃料制造过程模型,优化烧结温度、压制压力等参数,减少制造缺陷;

辐照行为预测:精准预测堆内辐照行为,为燃料换料与贮存提供依据;

贮存安全评估:开发乏燃料贮存性能模型,预测长期贮存中的腐蚀与临界安全风险;

后处理支撑:为乏燃料后处理提供性能数据,如裂变气体释放量、燃料溶解特性。

5.2.4 事故应急的智能化决策

在事故工况下,模型将成为应急决策的核心支撑工具:

事故进程快速预测:开发事故工况下的快速响应模型,在事故发生后 10 分钟内给出燃料状态预测;

缓解措施优化:模拟不同缓解措施(如注水、降压、安注)对燃料性能的影响,推荐最优应急方案;

后果评估与善后:评估事故后燃料破损程度与放射性释放量,支撑环境影响评估与善后处理。

5.3 行业生态演进:标准化与国际化合作深化

5.3.1 国际标准体系的融合

未来将逐步推进全球燃料性能模型标准的融合,减少技术壁垒:

核心参数统一:推动建立全球统一的燃料物性数据库与模型参数体系,如 IAEA 正在推动的 "全球燃料性能数据库"(GFPD)项目;

验证基准一致化:制定全球统一的模型验证基准题与评估方法,提高不同程序预测结果的可比性;

审评标准互认:推动中美欧等主要核电国家的模型审评标准互认,降低跨境取证成本;

先进堆标准制定:主导或深度参与快堆、SMRs 等先进堆型的模型标准制定,抢占技术话语权。

5.3.2 开源生态与国际合作深化

开源化将成为模型发展的重要趋势,促进全球技术合作:

开源程序发展:以 BISON、MOOSE 为代表的开源程序将吸引更多国家参与开发,形成全球协作的开源生态;

数据共享机制:建立 "公开数据 + 保密数据" 的分级共享机制,在保护知识产权的同时促进数据交流;

联合研发项目:开展跨国联合研发,如中美欧联合开展 ATF 模型开发,共享实验数据与技术成果;

能力建设合作:通过 IAEA 技术合作项目,向发展中国家提供模型技术培训与支持,推广先进模型应用。

5.3.3 人才与技术体系建设

面对模型技术的快速发展,需加强人才培养与技术体系建设:

跨学科人才培养:培养融合核科学、材料科学、计算机科学、人工智能的跨学科人才,满足模型发展需求;

实验平台建设:建设高水平燃料性能实验平台,如高燃耗辐照考验平台、ATF 专用实验装置;

自主化体系完善:中国等新兴核电国家需加快自主模型程序与数据库建设,打破国际技术垄断;

行业协同机制:建立 "监管机构 - 运营商 - 供应商 - 研究机构" 的协同机制,形成技术创新合力。

5.4 挑战与应对策略

5.4.1 核心挑战

1.基础科学瓶颈:高燃耗、极端工况下的燃料微观行为机理(如 HBS 形成机制、SiC 辐照损伤)尚未完全明晰;

2.数据匮乏制约:先进燃料与事故工况的实验数据稀缺,导致模型验证不足;

3.计算能力限制:多尺度多物理场耦合模型的计算复杂度极高,现有计算能力难以满足实时应用需求;

4.标准融合困难:不同国家的技术路线与利益诉求差异,导致国际标准融合进展缓慢;

5.人才缺口问题:跨学科高端人才短缺,制约模型技术创新与应用。

5.4.2 应对策略

1.强化基础研究:加大对燃料微观机理的基础研究投入,采用原位表征技术(如原位 TEM 辐照实验)揭示核心行为机制;

2.建设共享数据库:推动国际合作建设先进燃料与堆型的实验数据库,采用区块链技术实现数据安全共享;

3.发展高效计算技术:结合 GPU 加速、云计算、量子计算等先进计算技术,提升模型计算效率;

4.推动标准协调:通过 OECD/NEA、IAEA 等国际组织,建立政府、行业、科研机构的多层次标准协调机制;

5.完善人才培养:设立跨学科专业方向,开展校企合作培养,建立全球人才交流机制。

六、结论与展望

6.1 主要结论

燃料性能模型作为核电行业的核心技术支撑,经历了从经验模型到多尺度耦合智能模型的演进历程,形成了涵盖基础理论、工程程序、实验验证的完整技术体系。其具体实现方式呈现 "多物理场耦合、多尺度衔接、工程化适配" 的特征,不同堆型(PWR、BWR、快堆、HTGR)的模型适配形成了差异化技术路线。

当前燃料性能模型面临高燃耗裂变气体释放、PCI/FCMI 机理、事故工况预测、先进材料描述四大技术争议,监管机构、运营商、供应商、研究机构基于自身诉求形成了不同立场,国际标准体系存在显著区域差异。这些争议本质上是 "安全保守与经济高效"、"技术创新与工程成熟"、"数据独占与共享合作" 的矛盾体现。

未来燃料性能模型将向多尺度多物理场深度耦合、AI 与物理机理融合、数字化孪生集成三大方向发展,应用场景从安全评估拓展至燃料设计、运行优化、全生命周期管理的全链条,行业生态将呈现标准化融合与开源合作深化的趋势。

6.2 展望

燃料性能模型的发展水平将直接决定下一代核电技术的竞争力与安全性。到 2030 年,预计将实现三大突破:

1.技术突破:形成 "宏微观耦合 + AI 驱动" 的新一代模型体系,高燃耗工况预测误差降至 10% 以内;

2.应用突破:模型深度融入数字孪生,实现燃料全生命周期的智能化管理;

3.生态突破:建立全球统一的模型标准与数据共享机制,形成开放协同的全球技术生态。

核技术论坛

阅读 分享