在辐照条件下,溶液中会发生离子碘(I −)向元素碘(I 2)的转化现象,导致I 2从溶液中挥发出来,从而对事故源项产生重要影响,该现象在核电厂的事故源项分析中得到了极大的关注。本文对无辐照条件、辐照条件的碘化学模型进行了归纳总结,对美国核管会NUREG/CR-5950中的公式推导及反应率常数进行了介绍,可用于支持确定论事故后果分析,最后给出了总结和后续研究建议。 Radiolytic conversion of iodine ions (I −) to elemental iodine (I 2) occurs in aqueous solutions under radiation, resulting in the escape of I 2 from the solution. This phenomenon could influence the accident source term significantly so that it has to be considered seriously when performing accident source term analysis. In this paper, the iodine chemical model without irradiation condition and irradiation condition is summarized. The formula promotion and reaction rate constant in NUREG/CR-5950 of Nuclear Regulatory Commission of the United States are introduced, which can be used to support deterministic accident consequence analysis. Finally, a summary and suggestions for further research are given.
在辐照条件下,溶液中会发生离子碘(I−)向元素碘(I2)的转化现象,导致I2从溶液中挥发出来,从而对事故源项产生重要影响,该现象在核电厂的事故源项分析中得到了极大的关注。本文对无辐照条件、辐照条件的碘化学模型进行了归纳总结,对美国核管会NUREG/CR-5950中的公式推导及反应率常数进行了介绍,可用于支持确定论事故后果分析,最后给出了总结和后续研究建议。
核电厂,辐照,事故,离子碘向元素碘转化,模型,转化份额
Yaru Fu, Shengqin Gao, Qiliang Mei, Dawei Sun
Shanghai Nuclear Engineering Research & Design Institute, Shanghai
Received: Oct. 7th, 2020; accepted: Oct. 23rd, 2020; published: Oct. 30th, 2020
Radiolytic conversion of iodine ions (I−) to elemental iodine (I2) occurs in aqueous solutions under radiation, resulting in the escape of I2 from the solution. This phenomenon could influence the accident source term significantly so that it has to be considered seriously when performing accident source term analysis. In this paper, the iodine chemical model without irradiation condition and irradiation condition is summarized. The formula promotion and reaction rate constant in NUREG/CR-5950 of Nuclear Regulatory Commission of the United States are introduced, which can be used to support deterministic accident consequence analysis. Finally, a summary and suggestions for further research are given.
Keywords:Nuclear Power Plant, Radiation, Accident, Conversion of I− to I2, Models, Conversion Fraction
Copyright © 2020 by author(s) and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
事故条件下,对反应堆安全壳内碘化学进行分析是很重要的。在裂变产物中大量存在碘元素,其同位素元素将产生大剂量的辐射影响;且在一定条件下,极易挥发释放到环境中,从而对人类健康造成不可挽回的伤害。所以对于反应堆安全分析来说,准确估计安全壳内挥发性物质的化学形式、浓度及分布是十分有必要的。
在不同的反应堆系统和材料中,碘元素可以通过多种途径进行化学物理反应,其整体行为是非常复杂的。关于事故条件下水溶液中碘离子辐解分布已有大量的实验研究,研究表明,反应堆水环境中的碘离子(I‑)在特定的条件下可以转化为挥发性的元素碘(I2),元素碘挥发进入空气中成为气载放射源。为了完整的评估碘的综合影响,由于池水的pH、温度、初始碘浓度,有机物含量等条件都会在一定程度上影响该转化效应,多种条件的相互影响更使转化复杂化。
目前,研究者对碘辐射化学反应中各因素(如pH、初始碘浓度、剂量率、温度、溶解气体和杂质等)的影响开展了大量工作 [
另一方面,在综合严重事故模拟软件(COCOSYS [
根据已有的反应堆碘化学的研究中,可以建立两类碘化学系统:1) 无辐照下的碘化学系统;2) 有辐照下的放射性碘化学。有辐照下的碘化学系统在无辐照系统的基础上,将产生新的反应通道,对碘转化效应产生重要影响。但同时,有辐照下的系统分析又离不开无辐照系统的理论基础。本文将对无辐照条件、辐照条件的碘化学模型进行归纳总结,对美国核管会NUREG/CR 5950 [
碘元素的化合价在−1到+7价之间变化,在水中,主要有I− 、I2、HOI、HIO2、IO2、HIO3、HIO4 。其中+7价的化合物化学性质较为稳定,比较难被氧化和缓解,所以暂不考虑其化学变化。碘元素在无辐照水溶液中的完整行为模型可以用水解反应(1),电离(2)和歧化反应(4)等来描述 [
I 2 + H 2 O ⇌ I − + H O I + H + (1)
H O I ⇌ O I − + H + (2)
I 2 + I − ⇌ I 3 − (3)
H O I + H O I ⇌ I − + H I O 2 + H + (4a)
H O I + H I O 2 ⇌ I − + H I O 3 + H + (4b)
其中分子碘的水解反应是最重要的热反应之一,而次碘酸和次碘根离子的歧化反应进一步促进了分子碘浓度的减小,使整个反应向右进行。
对于无辐照条件下的基础完整反应关系也在图1中给出。反应(1)~(4)足够完整的来描述纯水中的碘化学反应关系。除上述基本的反应外,图1中还描述了其他的附加反应,包含 I 2 O H − , H 2 O I + 等其它碘化合物。
当加入辐照条件后,水中会产生过氧化氢(H2O2)、自由基产物(OH、H)等,其中自由基产物的浓度非常低,约10−10M或更低,而过氧化氢则会增加到与水池内碘化物相当的浓度(10−4到10−6 M)。最终达到的浓度取决于辐照剂量以及与溶液中其他物质的反应程度。因为溶液中过氧化氢的存在,在反应中引入新的氧化还原反应,包括:
IO − + H 2 O 2 ⇌ I − + O 2 + H 2 O (5)
HOI + H 2 O 2 ⇌ H + + I − + O 2 + H 2 O (6)
I − + H O I ⇌ I 2 O H − (7a)
I 2 O H − + H 2 O 2 ⇌ 2 I − + O 2 + H + + H 2 O (7b)
IO 3 − + H 2 O 2 → I O 2 − + H 2 O + O 2 (7c)
IO 3 − + H + + H 2 O 2 → H I O 2 + H 2 O + O 2 (7d)
HIO 2 + H 2 O 2 → H I O + O 2 + H 2 O (7e)
HIO 2 + I − + H + → 2 H I O + H 2 O (7f)
HIO + I − + H + ⇌ I 2 + H 2 O (7g)
氧化反应包括:
IO 3 − + H 2 O 2 → I O 4 − + H 2 O (8a)
HIO 2 + H 2 O 2 → I O 3 − + H + + H 2 O (8b)
OH + HOI → IO + H 2 O (8c)
I + HOI → IO + I − + H + (8d)
I 2 − + HOI → IO + 2 I − + H + (8e)
辐照的溶液中,整体的碘化学关系图见图2。在众多反应关系中,重点关注次碘酸和双氧水之间的还原反应(5)和(6),此反应与氧化反应(4)形成竞争反应,抑制了HIO3 的形成,并很大程度上影响了分子碘的浓度大小,进而影响了碘离子向分子碘的转化份额。
图1. 无辐照条件下溶液中碘化合物关系图
图2. 辐照条件下溶液中碘化合物关系图
Abel [
[ H + ] 2 × [ I − ] 2 [ I 2 ] = a + b × [ H + ] (9)
式中, [ ] 代表该括号内粒子的浓度(mol/L),a、b为与温度有关的常数(25℃时, a = 6.05 × 10 − 14 , b = 1.47 × 10 − 9 )。
已知溶液的pH值时,H+的浓度可由下式计算得到:
[ H + ] = 10 − p H (10)
其推导过程从I−与H2O2的各反应平衡中得到,如下:
① 反应率描述方程
对于氧化反应:
I − + H 2 O 2 → IO − + H 2 O (11)
有 r 1 = k 1 × [ I − ] × [ H 2 O 2 ] (12)
对于还原反应:
IO − + H 2 O 2 ⇄ k 5 I − + O 2 + H 2 O (13)
HOI + H 2 O 2 ⇄ k 6 H + + I − + O 2 + H 2 O (14)
有 r 2 = k 6 × [ HIO ] × [ H 2 O 2 ] + k 5 × [ IO − ] × [ H 2 O 2 ] (15)
式中, r 1 为氧化反应(11)的反应速率, r 2 为还原反应(13)和(14)的总反应速率, k 1 为反应(11)的反应率常数, k 5 、 k 6 分别为反应(13)和(14)的反应率常数。
② 稳定状态意味着氧化反应和还原反应的反应率相同,即 r 1 = r 2 ,所以有:
k 1 × [ I − ] × [ H 2 O 2 ] = k 6 × [ HIO ] × [ H 2 O 2 ] + k 5 × [ IO − ] × [ H 2 O 2 ] (16)
等式两边除去相同的双氧水浓度,有:
k 1 [ I − ] = k 6 [ HIO ] + k 5 [ IO − ] (17)
③ 次碘酸浓度的推导
根据反应关系式(1)和(2)可以得到:
[ HIO ] = K 1 × I 2 [ H + ] × [ I − ] (18)
[ IO − ] = K 2 × [ HIO ] [ H + ] = K 2 × K 1 × [ I 2 ] [ H + ] 2 × [ I − ] (19)
其中,
K 1 为反应(2)的平衡常数;
K 2 为反应(1)的平衡常数。
④ 将式(18)和(91)带入式(17),得到:
k 1 × [ I − ] = k 6 × K 1 × [ I 2 ] [ H + ] × [ I − ] + k 5 × K 2 × K 1 × [ I 2 ] [ H + ] 2 × [ I − ] = K 1 × [ I 2 ] [ H + ] 2 × [ I − ] × ( k 6 × [ H + ] + k 5 × K 2 ) (20)
进一步简化得到方程(9)的形式
[ H + ] 2 × [ I − ] 2 [ I 2 ] = K 1 k 1 × ( k 6 × [ H + ] + k 5 × K 2 ) (21)
⑤ 根据方程形式得到:
a = K 1 × K 2 × k 5 k 1 (22)
b = K 1 × k 6 k 1 (23)
这些反应平衡常数的值在25℃下较好获得,且有较多研究围绕类常数进行了测量,表1给出NURER/CR-5950 [
方程 | 参数 | 数值(t = 25˚C) | 作者 | 年份 |
---|---|---|---|---|
I − + H 2 O 2 → IO − + H 2 O | k 1 | 0.012 L/mol∙s | H. A. Liebhafsky | 1932 |
IO − + H 2 O 2 → I − + O 2 + H 2 O | k 5 | (6.6 ± 2.0) × 107 L/mol∙s | H. Shiraishi | 1992 |
HOI + H 2 O 2 → H + + I − + O 2 + H 2 O | k 6 | 37 L/mol∙s | H. A. Liebhafsky | 1932 |
I 2 + H 2 O ⇌ HOI + I − | K 1 | 4.77 × 10−13 | D. A. Palmer | 1985 |
HIO ⇌ IO − + H + | K 2 | 2.31 × 10−11 | M. Eigen | 1962 |
a = K 1 × K 2 × k 5 k 1 | a | (6.05 ± 1.83) × 10−14 | C. F. Weber | 1992 |
b = K 1 × k 6 k 1 | b | 1.47 × 10−9 | C. F. Weber | 1992 |
表1. NUREG-5950中引用的各反应常数
本文对无辐照条件、辐照条件的碘化学模型进行了归纳总结,对美国核管会NUREG/CR 5950中的公式推导及反应率常数进行了介绍,并给出了总结和后续建议,可用于支持确定论事故后果分析。
国际上各主要事故分析软件中考虑的离子碘向元素碘转化模型所考虑的化学反应是有所不同,但重点关注的反应是一致的;各数学模型中的反应率常数或平衡常数是有一定的差别的,后续可重点对反应的反应率常数或平衡常数进行进一步的研究,重点关注温度的影响,用于模型改进。
付亚茹,高圣钦,梅其良,孙大威. 离子碘向元素碘转化行为模型研究Study on Models for the Conversion Behavior of Iodine Ions to Elemental Iodine in Aqueous Solutions[J]. 核科学与技术, 2020, 08(04): 237-243. https://doi.org/10.12677/NST.2020.84028