第 38 卷第 4 期核技术 ol.38, No.4 05 年 4 月 NUCLEAR TECHNIQUES April 05 压水堆燃料棒裂变产物释放产生比计算方法 许锐 董冰 付亚茹 梅其良 ( 上海交通大学上海 0040) ( 上海核工程研究设计院上海 0033) 王德忠 摘要为对压水堆一回路源项进行准确分析, 在研究了美国核管制委员会 (Nuclear Regulatory Commission, NRC) 发布的计算气体裂变产物释放产生比 (Release to irth Ratio, R/) 的 ANSI/ANS-5.4 标准及法国 METEOR.5 程序计算气体裂变产物释放产生比的方法后, 分别利用几种方法计算了燃料棒在不同功率水平以及不同燃耗下的几种气体裂变产物的释放产生比, 对几种方法进行了研究分析 结果表明, 燃料棒的裂变产物释放产生比会随着其功率以及燃耗而增大, 相较于美国 ANSI/ANS-5.4-0 标准, 法国 METEOR.5 程序方法更为保守, 而 NRC 原有的 ANSI/ANS-5.4-98 标准最为保守, 计算出的释放产生比最大,ANSI/ANS-5.4-0 标准能较好地适用于压水堆核电站裂变产物释放份额的计算 关键词燃料棒, 裂变产物, 释放产生比,ANSI/ANS-5.4 标准,METEOR.5 程序方法中图分类号 TL99 OI: 0.889/.053-39.05.hs.38.040605 反应堆在运行过程中燃料芯块会释放出大量的裂变产物到燃料芯块与包壳的间隙之中, 当包壳发生破损时这些裂变产物又会通过包壳释放到一回路冷却剂中, 释放到一回路冷却剂中的裂变产物成为核电厂的主要源项 这些裂变产物按其挥发性的差异可分为气体裂变产物 ( 主要为 Kr I Xe 等惰性气体和卤素 ) 半挥发性裂变产物 (a Ru Sr 等 ) 和难挥发性裂变产物 (Rh Pd Tc 等 ) [], 由于气体裂变产物较易释放, 因而在进行源项分析时考虑的主要是气体裂变产物 在分析核材料裂变生成的气体放射性核素时, 尤需关注的是释放率的问题 [] 裂变产物的释放产生比 (Release to irth Ratio, R/) 是某种裂变产物从燃料芯块中释放到燃料芯块与包壳间隙的释放量与芯块产生的该种裂变产物产生量之比, 其值的大小决定了裂变产物从燃料芯块释放到芯块与包壳间隙中的多少, 对其进行准确计算对确定核电厂一回路源项具有重要意义 目前美国 法国等核电发达国家提出了多种计算裂变产物释放产生比的方法, 这些方法计算出的结果大多较保守, 美国核管制委员会 (Nuclear Regulatory Commission, NRC) 在 0 年发布了 ANSI/ANS-5.4-0 标准, 该标准详细介绍了 NRC 推荐的计算裂变产物释放产生比的方法, 且在原有的标准基础上进行了最优化改进和参数修正, 提高了计算的准确性 为对压水堆一回路源项进行准确分析, 确定较为准确的计 国家自然科学基金 (No.758) 上海市教育委员会科研创新项目(No.ZZ0) 资助第一作者 : 许锐, 男,990 年出生,0 年毕业于上海交通大学, 现为硕士研究生通讯作者 : 王德忠,E-mail: dzwang@stu.edu.cn 收稿日期 :04--6, 修回日期 :05-0-05 040605- 算裂变产物释放产生比的方法, 本文对 NRC 发布的 ANSI/ANS-5.4-0 标准和原来的 ANSI/ANS-5.4-98 标准以及法国 METEOR.5 程序方法进行了研究分析和对比 扩散与释放模型的建立与求解 由于 UO 核芯颗粒可近似为球形, 因此裂变产物在其中的扩散可用球坐标下的扩散方程描述 [3] 主裂变产物在 UO 颗粒中的扩散遵循菲克扩散定律 [4] : C (rc) =λc+ () t r r 式中,C 为燃料中的核素浓度,cm 3 ;t 为时间,s; 代表核素产生率,cm 3 s ;λ 代表核素衰变常数, s ;r 是燃料球的半径,cm; 为核素的扩散系数, cm s 当反应堆功率基本稳定时可以将芯块晶粒中裂变气体产物浓度视为稳定, 有 : C 0 t 于是基本方程可简化为 : ( rc) C 0 r r 令 y=rc, 则有 : () (3)
核技术 05, 38: 040605 y y r r 解以上常微分方程可得 : 0 r e r (4) y C C e r (5) 当 r=0 时,C +C =0, 即 C = C ; 当 r=a 时, C=0, 则有 : ( rc) a (6) r ra 由式 (5) (6) 可得, 燃料晶粒的裂变产物释放产生比 : C 4πa R r ra 4 π 3 a 3 C a a a a C C C e e e e 3 a a a a a a ( Ce Ce ) (7) 令 μ=a /( 无量纲化 ), 则释放产生比为 : R 3 coth (8) 对于球形晶粒半径 a 有 :a=3/s( 为晶粒体积,S 为表面积 ), 则式 (8) 可变形为 : R S [coth ] (9) 对于尺寸极小的晶粒来说其 μ 值极小, 因此式 (9) 可近似为 : R S (0) 式 (0) 是式 (9) 相对保守的近似, 其误差小于 % 裂变产物释放产生比计算方法. ANSI/ANS-5.4 标准 美国 NRC 发布了用于计算反应堆功率稳定情况下气体裂变产物释放产生比的 ANSI/ANS-5.4 标准, 其最新版本为 ANSI/ANS-5.4-0 ANSI/ANS-5.4-0 在旧的较保守的 ANSI/ANS- 5.4-98 基础上进行了估算方法的优化 参数的修正以及适用范围的扩充 (ANSI/ANS-5.4-98 标准仅适用于燃耗在 GWd tu 以下的 UO 燃料棒的计算 [5] ) 在 ANSI/ANS-5.4-0 中引入了先驱核因子的无量纲修正因子 α, 即式 (0) 被修正为 : R S () ANSI/ANS-5.4-0 中的扩散系数 是根据 [6] Turnbull 等在 98 年提出的关系式修正而来, 即 : i 7.60 exp( 35000 Ti) 5 0.5.40 F exp( 3800 Ti ) 40 0 F () 式中, i 为扩散系数,m s ;T i 为燃料温度 ;F 代表裂变率 第一项为高温时芯块的固有扩散系数 ; 第二项为因辐照引发铀的空位机制而增加的扩散系数 ; 第三项为因辐照直接引发的扩散系数 ANSI/ANS-5.4-0 还用了一个简单的经验关系式 itanza 临界准则来考虑晶间的气泡联合效应 即燃耗 u 8. MWd kg U 时, 晶间气泡开始联合的温度为 : 9800 Tlink 73 (3) ln(76 u) 当燃耗 u>8. MWd kgu 时, 晶间气泡开始联合的温度 : Tlink 434.85u 73 (4) 式中,T link 是在给定的计算节点处的温度,K 由以上关系确定的燃料节点 i 处的半衰期小于 6 h 裂变产物释放产生比为 : R S nuclidei i,nuclide i nuclide (5) 半衰期大于 6 h 但小于 60 d 的裂变产物释放产生比为 : 0.5 R nuclide 85m Kr S nuclidei i,nuclide nuclide 85m Kr i nuclide (6) 式中, 当 T i T link 时,S/=0 cm, 当 T i >T link 时, S/=650 cm,t i 为 i 节点处的燃料温度. 法国计算方法 法国 METEOR.5 程序中, 其裂变气体扩散系数由热项 无热项以及联合项三项构成 [7] 热项 : 355 0 th 7.60 exp Tk (7) 040605-
许锐等 : 压水堆燃料棒裂变产物释放产生比计算方法 无热项 : 联合项 : ath. 0 39 F (8) 3870 5 0.5 mix 6.950 F exp Tk (9) pellet P ave P P node rod ring ring m R i, i, k pellet node nuclide P i R i, (0) 式中,F 为裂变密度,m 3 s ;T k 为燃料的局部温度,K 3 整根燃料棒释放产生比计算模型 利用以上方法进行计算时先将燃料棒的径向和轴向进行网格划分使燃料棒划分为若干个环状单元, 计算时假定每个环状单元内的温度以及燃耗相等, 据此分别计算出每个环状单元的释放产生比, 然后根据每个单元处的功率水平 体积以及释放产生比进行加权平均即可得到整个燃料棒的释放产生比 对于径向节点为 k 轴向节点为 的燃料棒, 其释放产生比可表示为 : 式中,P 是功率,kW ft,ave= 棒的平均值,pellet= 轴向平均值,ring= 径向环 ; 为体积,cm 3,rod= 棒的总体积,node= 轴向节点体积,ring= 径向环体积 燃料棒的网格划分如图 所示, 整根燃料棒的释放产生比的 MATLA 计算流程如图 所示 图 燃料棒网格划分 Fig. Grid of the fuel rod. 图 计算燃料棒的释放产生比流程 Fig. Calculation process of the release to birth ratio. 040605-3
核技术 05, 38: 040605 4 结果分析 4. 不同功率水平下的释放产生比 研究 UO 燃料棒功率水平对于燃料棒裂变产物释放产生比的影响, 以 ANSI/ANS-5.4-0 中数据匹配较好的核素 85m Kr (4.48 h) 和较为常见的半衰期较长的核素 3 I (8.04 d) 为研究对象, 分别利用 ANSI/ANS-5.4-0 标准 ANSI/ANS-5.4-98 标准以及法国 METEOR.5 程序的方法计算百万千瓦级压水堆核电站的 UO 燃料棒从 0% 功率到满功率时 85m Kr 以及 3 I 的释放产生比, 得到的结果如图 3 98 标准以及法国 METEOR.5 程序的方法计算 UO 燃料棒燃耗从 GWd tu 到 60 GWd tu 时 85m Kr 以及 3 I 的释放产生比, 并将不同燃耗下计算得到的 3 [8] I 的释放产生比结果与 Halden 实验所得数据进行对比, 最终得到的结果如图 4 所示 结果表明,ANSI/ANS-5.4-0 标准的计算结果与实验结果符合最好, 在几种方法中能最准确地计算出裂变产物的释放产生比 当 UO 燃料棒燃耗增加时其裂变产物的释放产生比会增大, 在不同的燃耗下法国 METEOR.5 程序的计算方法相较于 ANSI/ANS-5.4-0 标准更保守, 在 ANSI/ANS-5.4-98 标准适用燃耗范围内 ANSI/ANS-5.4-98 标准最保守 图 3 85m Kr (a) 和 3 I (b) 在不同功率水平下的释放产生比 Fig.3 Release to birth ratio of 85m Kr (a) and 3 I (b) at different power. 将几种方法的计算结果进行对比, 结果表明, 当燃料棒功率水平增加时其裂变产物的释放产生比也增大, 在不同的功率水平下法国 METEOR.5 程序的计算方法相较于 ANSI/ANS-5.4-0 标准更保守, 功率水平较高时 ANSI/ANS-5.4-98 最保守 4. 不同燃耗下的释放产生比 为研究燃耗对于裂变产物释放产生比的影响, 分别利用 ANSI/ANS-5.4-0 标准 ANSI/ANS-5.4- 图 4 不同燃耗下 85m Kr (a) 和 3 I (b) 释放产生比 Fig.4 Release to birth ratio of 85m Kr (a) and 3 I (b) at different burn-up. 5 结语 根据以上计算结果对比分析得出如下结论 : () 不同功率水平下裂变产物释放产生比的计算表明, 当 UO 燃料棒燃耗一定时其裂变产物的释放产生比会随着功率水平而增大 () 不同燃耗下裂变产物释放产生比的计算表明, 当 UO 燃料棒处于满功率运行状态时, 其裂变产物的释放产生比会随着其燃耗而增大 (3) 几种方法计算结果与实验数据的对比表明,ANSI/ANS-5.4-0 标准的计算结果与实验结 040605-4
许锐等 : 压水堆燃料棒裂变产物释放产生比计算方法 果符合最好, 在几种方法里能最准确地计算出裂变产物的释放产生比 METEOR.5 程序的计算方法相对于 ANSI/ANS-5.4-0 标准更保守, 而 ANSI/ANS-5.4-98 标准在几种方法中最保守 ( 除低功率水平外 ) 综上所述,ANSI/ANS-5.4-0 在原来方法的基础上进行了优化修正, 不但对适用范围进行了拓展, 其计算结果相对于原有标准也更准确, 能较好地适用于压水堆核电站裂变产物释放份额的计算 参考文献 rillant G, Marchetto C, Plumecocq W. Fission product release from nuclear fuel I: physical modelling in the ASTEC code[j]. Annals of Nuclear Energy, 03, 6: 88 95 伍怀龙, 郝樊华, 唐元明. 反应堆中气体裂变产物释放问题的研究 [J]. 核技术, 007, 30(7): 633 636 WU Huailong, HAO Fanhua, TANG Yuanming. Investigation of fission gases release in nuclear reactors[j]. Nuclear Techniques, 007, 30(7): 633 636 3 房超. 半衰期为数天的惰性气体和卤素裂变产物在 UO 颗粒中的扩散与释放 [J]. 原子能科学技术, 03, 47(): 09 3 FANG Chao. iffusion and release of noble gas and halogen fission products with several days half-life in UO particle[j]. Atomic Energy Science and Technology, 03, 47(): 09 3 4 erfondern K, Nabielek H, Moormann R, et al. Fuel performance and fission product behavior in gas cooled reactors[r]. IAEA-TECOC-978, ienna: IAEA, 997 5 ANSI/ANS-5.4-98, Method for calculating the fractional release of volatile fission products from oxide fuel[s]. U. S. NRC, 98 6 ANSI/ANS-5.4-0, Method for calculating the fractional release of volatile fission products from oxide fuel[s]. U. S. NRC, 0 7 何晓军, 陈彭. 燃料元件行为分析程序 METEOR l.5 的校核 [J]. 核动力工程, 007, 8(): 5 HE Xiaoun, CHEN Peng. alidation of fuel behavior analysis code METEOR.5[J]. Nuclear Power Engineering, 007, 8(): 5 8 Turnbull J A, eyer C E. ackground and derivation of ANS-5.4 standard fission product release model[m]. U. S. NRC, 00 Methods for calculating the release to birth ratio of volatile fission products XU Rui ONG ing FU Yaru MEI Qiliang WANG ezhong (Shanghai Jiao Tong University, Shanghai 0040, China) (Shanghai Nuclear Engineering Research & esign Institute, Shanghai 0033, China) Abstract ackground: A large quantity of fission products would be produced when the nuclear reactor runs. These fission products would be released to the pellet-rod gap and then to the primary circuit while the fuel rod is broken. The fission products in the primary circuit become the main source term of the reactor. Purpose: This study aims to analysis the source term accurately by calculating the release to birth ratio (RR) of volatile fission products precisely. Methods: After researching American NRC ANSI/ANS-5.4 methods and French methods for calculating the RR of volatile fission products, we calculated the RR of volatile fission products of a fuel rod at different power and different burn-up by these methods, and made some comparison and analysis of these results. Results: Computation results showed that the RR of volatile fission products would increase with the power and burn-up of the fuel rod. Compared with the ANSI/ANS-5.4-0, the French method is more conservative whilst the ANSI/ANS-5.4-98 is the most conservative among these methods. Conclusion: ANSI/ANS-5.4-0 is the best method among these three methods and it can be used to calculate the release to birth ratio of volatile fission products with good results. Key words Fuel rod, Fission products, Release to birth ratio, ANSI/ANS-5.4, METEOR.5 CLC TL99 040605-5