核爆放射性核素的分凝研究
合集下载
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
9 5:
,
粒子分布离散为l 计算各组质量份额‘ , [ 组, a)
由火球温度函数和粒子凝固温 度, 计算粒子冷凝时刻‘
由B t a ∞m议程计算 c 时各衰变链核素原子甄 对 比粒子凝固温度计算各链难熔份额F I
妯
"/ 例\" Z r
’ b ’ N 3d 5
 ̄eig S i 公式计 算尺度k rl n 中链j 的分布份额H0) 【
表 2 核素分布计算条件
=
F= r 所有尺度组总的当量裂变数 ; F( = i也) 尺度组 k 里质量链 i 的当量裂变 数;
= 质 量链 的裂 变产额 。 i 根据 Felg的推导 , . 粒子 中质 量链 rin i 第 i } 组
i 的当量裂变数为n :
( )=
在 核 爆火 球 冷 却 的过 程 中 , 为 主要 组 分 作 的 非放 射 性物 质 蒸 气 通过 成 核 、 冷凝 和凝 并 等
素 首 先 冷凝 到 粒子 上 , 然后 难 熔 或 易熔 元 素 相 继 冷凝 到 粒子 上 , 发性 气体 元 素 除非 它 们 很 挥
作用形成载体粒子 。同时, 放射性裂变产物在 载体粒子表 面发生冷凝。在冷凝过程 中, 由于
( £ )=∑ ∑ H x( A ) 3 c e 一 r ;() p t
l0 l I
n—l
n
C =—; L一 h 『上 。 () 4
2 分 凝 模 型
基 于放 射 性 衰 变 链 的 B t a a m n公 式 和 e
Fel g的核素 半经 验 分凝 模 型 , 建 核爆 粒 子 rin i 构
1 5 0 K时 , 7 估计 大 多数难熔 元素 , 锫和稀 土 , 如
载体 粒 子 上 分 布 不 同, 种 现 象 称 为 “ 这 分
凝” …。分凝效应使得裂 变产物在放射性粒子 上 的分布不均匀 , 使得核爆放射性沾染场 的辐 射性质在不 同地点和不同时间不同。在给定 的 核爆袭击情景下, 评估预期伤亡时 , 分凝效应有 很 大 的影 响 。
3 33 . 3 3× 1 2 0-
于气 化 土 壤 , 那么 将 有 成千 上 万 吨土 壤卷 入火
1 分凝效应
1 1 分 凝机 理 .
在 核爆 粒子形 成 过程 中并 不是 所有 的放 射 性 核 素 同时 冷凝 到粒 子 上 , 些非 常难 熔 的元 一
球 中。气化了的土壤分压将是裂变产物分压的
一
百多万倍 。较难熔的成分, 如二氧化硅 , 同沸
放射性衰变链 的 Bt a 公式构建 核爆 粒子放射性活 度分布的计算模型 。计算 了给定载体物质 和粒子 am n e 尺度分布条件下的粒子放射性活度分 布 , 果表 明挥发 性核 索呈 面分布 特征 , 结 难熔性 核素呈 体分 布特 征 。分析了粒子尺度 、 固温度 和凝 固时间对 放射性活 度分布的影响程度 , 凝 结果表明粒子尺度分 布和凝
怒 ( 扎 (
() 7
粒子质量份额
3 33 . 3 3× 1 - 0 2
3 33 . 3 3× 1 2 0-
其中:
表 3 总 B活度和质量链 8 ,5的分 布 99
粒子尺度/ m
3. 6 ×1 3—1 6 3 0×1 3 6 89 0- 2 0
.
总 B活度/ q B
5 6× 1 1 .1 02
1 O ×1 2 .0 0 ̄ 1 O ×1 2 .0 0o
1 01×1 2 . 0 ̄
I 0 3 0×1 3—7 5 4 6×1 ‘ . 2 0- 4 0一
.
3 3 3 3×1 2 .3 0- 3. 3 ×1 2 3 33 0- 3. 3 ×1 2 3 33 0- 3. 3 ×1 2 3 33 0-
郭立军 陈琳2郑毅。 , ,
(. 1海军 陆战学 院 , 广州 5 0 3 ;. 14 0 2 防化 指挥工程学 院 , 北京 12 0 ;. 0 2 53 防化研究 院 , 北京 1 2 5 2 0 0)
摘 要 : 了研究核爆粒子 中放射性核素 分凝 的影 响因素和规律 , 为 基于 Fe i r l g的半经验分凝模 型和 in
核素 不 同的挥 发性 而导致 在不 同尺 度 和性质 的
快衰变为其他元素 , 否则不冷凝到粒子上 。在 粒子未凝固之前冷凝 的核素 , 近似均匀分布在 粒子体 内。粒子凝 固之后冷凝 的核素, 主要分 布在粒子表面。 从 物 态相平 衡 的观点 来 看 , 凝 取 决 于冷 冷 凝 相物 质 的沸点 和蒸气 相 中的分 压 。火球 热力 学温度下降到 3 0 400K时 , 0~ 0 0 其总压力大 致与周 围大气 平衡 J 。以 l00 k 爆 炸为 0 t 例, 估计此时裂变产物分压 为 l “大气压。在 0
n ( 一 A)
式中: () £是链 中时刻 t n个核素的原 第 子数 ; 是链 中 0时刻第 ( n 个元素 的原 i ) 子数 ; A 是链 中第 七 k ) ( sn 个核 素的衰 变常 数。
t 时刻 , 素 的放 射性 活度 为 : 核 A ()= A () t t。
.
4. 7×1 2 5 01 4. 4×1 1 6 02 4 6 ×1 2 .9 01 4 7 ×1 2 .4 01 4 7 ×1 2 .9 0l 4. 3 ×1 2 8 01 4. 8 ×1 2 8 01 4. 2× 1 1 9 02 4 9 ×1 1 .6 02 5 01×1 2 . 01 5 0 ×l 2 .6 0l 5. O ×1 2 1 01
.
∑
^= 1
系数相关性相符 , Fe i 半径指数模型 的基 是 r lg i n 本依据。 将粒子群根据尺度大小分为若干组, : 令
=
D取值 10 m。那 么 第 后组 粒 子 中质 量 0 链i 的原 子数 为 : N ( ̄ id)= F( ) i (O 1)
k 尺度组粒子的平均直径 : 尺度组 的粒子数 ;
3 模 拟计算
给定粒子凝固时刻 的时间和温度 、 射性 放 粒子的总量和尺度分布参数 , 模型可以计算 总 的放射性活度的尺度分布和单个质量链的尺度 分布 。表 2 列出计算条件。表 3列出总的放射 性活度分布 和质量链 8 , 9 质量链 9 5的分 布。 图 2和 图 3 示 为总 的 p比活 度 和质 量链 8 、 所 9 9 5的比活度分布。
由B 她m l 蛆议程计算指定时刻t 各衰
变链的放射性活ቤተ መጻሕፍቲ ባይዱ
1个 沉 降 样 本 中 给 定 核 素 i的 裂 变 数
为 : = a 。 / () 1
由各链 的放射性和 分布份额计算 时囊 t I 总的 放 射性分 布呀指 定衰变 链的分布
式 中 , 样本 中核 素 i a是 的原 子数 , 是 核 素 的裂 变 产 额 。那 么 核 素 i 对 于核 素 的 相 分凝 系数 定义 为 :
E = —— i —L~
h=l
;
() 8
(/ 一 )  ̄ 1次幂, b= 令 i
对数 相关 性为 :
lr 9 n i5=
,
, 推出分凝系数的
∑ ( ) )
厶 ( ): 一; () 9
( 1一b)n8 +a。 ilr 9 i () 6 这 与 Feig从 放 化 分析 得 到 的核 素 分 凝 rin l
其 中的一 种 或几 种 核 素 ( 括 Ⅳ1可 以直 包 ) 接 由裂 变产 生 。除 了 Ⅳ , m 每一 种 核素 按 照衰 变 律 ep(一Af 衰 变 。 除 了 Ⅳ 每 种 核 素 按 照 x i ) ,
时刻衰变链 中难熔性核素所 占份额 ( ) 决定
73 7
整个 链 的 冷 凝 行 为 。 Feig 提 出 给 定 质 量 rin l 链 i的 放 射 性 活 度 分 布 正 比 于 粒 子 直 径 的
r J= 。 () 2
图 1 计 算 程 序 流 程
ep A t增长。衰变链中每种核素的生长 x (一 ) 和 衰减历 程 可 以用 B t n8方程 表示 : ae [ ma
通常将产额 大、 半衰期 长 的曲 r z 分 s 和 r 别作为典型挥发性和难 熔性的参考核素 , 分析 核素相 对参 考核 素 的分 凝 系数 关 系。
核素 的分凝 程度 用 相对 于不 分凝 核素 的分
。
凝系数 r 表示 。所谓不分凝 核素是指一些难 熔 的 核素 , 们 的母 体 尽 管是 一 些 易 熔 或 挥发 它 性核素 , 但是半衰期很短 , 很快就会衰变成为难
熔 核 素 , 为 它 们 不 发 生 分 凝 。 例 如 质 量 链 认
2 2 核素 分 凝 .
放射性活度分 布的计 算模 型, 计算流程如 图 1
所示 。
2 1 核 素衰 变 .
放射 性核 素 在 冷 凝 的 同时 发 生 着 衰 变 , 使
其浓度不断变化。计算分凝效应必须考虑每种 核素的衰变和生长 , 以确定每种核素在冷凝时 刻的浓度 。单个放射 性核 素按照其 半衰期 衰 减, 而裂变产物往往形成许多 由相同质量数核 素组 成 的衰变 链 , 如某 个 衰 变链 可 以表示 为 :
收 稿 日期 :02— 1— 0 2 1 0 2
点差不多、 分压低得多的裂变产物相 比, 冷凝得
早得多。
作者简 介 : 郭立军 , , 男 黑龙 江拜 泉人 , 讲师 , 士研 硕
究生, 研究 方向为核辐射探测 。
77 2
表 1 裂变产物 中各元素 的冷凝性质 ‘
1 2 分凝程 度 .
固温度的影响较大 , 子的放射性活 度分布与爆 炸当量 、 粒 比高 以及爆炸地点 的地理条件密切相关。
关键词 : ; 核爆 放射性核素 ; 分凝
中图分类号 : T 1 L8 文献标 志码 : A 文章编号 : 0 5 -94 2 1 ) 70 7 -6 2 80 3 (0 2 0 - 20 7
本 文应 用数 值模 拟 的方法 计算 核爆 粒 子 的
放射性活度分布 , 分析放 射性核素分凝作用 的
影 响 因素和 规律 。
平 衡蒸 气 压 在 l 大 气 压 以下 , O 它们 会 发 生冷
凝 。而卤素、 惰性气体等元素沸点较低 , 平衡蒸
汽压 在 l 大气 压 以上, 们 不 会 发 生 冷 0 它 凝 。各种裂变质量链 中元素 的沸点见表 l 】 。 对于地爆而言, 如果有百分之几 的爆炸能量用
第3 2卷
21 0 2年
第 7期
7月
核 电子学 与探 测技 术
Nu la e to i s& De e t n Te h oo y c e rEl cr n c tc i c n lg o
V0. 2 No 7 13 .
J l. 2 1 uy 0 2
核 爆 放 射 性 核 素 的分 凝研 究
4. 8×1 }1 4 (2
质量链 8 / q 9B
2. 3 x1 2 O 2 3 ×1 9 .2 01 2 41×1 9 . 01 2. 9 ×1 9 4 01 2. 8 ×l 5 09
质量链 9/ q 5 B
99 09 . 7 x1
1 6 3 0 ×1 3—1 0 3 0×1 3 .2 0- 2 0
Nl N2_ +N3 … N
m
() 5
Fei n 以土 壤粒 子 凝 固的 温度为 rlg t in 界, 将沸点低于该温度的核素归为挥发性核素 ,
沸点高于该温度的核素归为难熔性核素 。挥发 性核素仅仅冷凝在粒子表面 , 而难熔性 核素均
匀分 布在 粒 子体 内。任意 裂变 产物 衰变 质量链 的组 分 数 量 都 在 随 时 间变 化 , 土壤 凝 固温度 在
,
粒子分布离散为l 计算各组质量份额‘ , [ 组, a)
由火球温度函数和粒子凝固温 度, 计算粒子冷凝时刻‘
由B t a ∞m议程计算 c 时各衰变链核素原子甄 对 比粒子凝固温度计算各链难熔份额F I
妯
"/ 例\" Z r
’ b ’ N 3d 5
 ̄eig S i 公式计 算尺度k rl n 中链j 的分布份额H0) 【
表 2 核素分布计算条件
=
F= r 所有尺度组总的当量裂变数 ; F( = i也) 尺度组 k 里质量链 i 的当量裂变 数;
= 质 量链 的裂 变产额 。 i 根据 Felg的推导 , . 粒子 中质 量链 rin i 第 i } 组
i 的当量裂变数为n :
( )=
在 核 爆火 球 冷 却 的过 程 中 , 为 主要 组 分 作 的 非放 射 性物 质 蒸 气 通过 成 核 、 冷凝 和凝 并 等
素 首 先 冷凝 到 粒子 上 , 然后 难 熔 或 易熔 元 素 相 继 冷凝 到 粒子 上 , 发性 气体 元 素 除非 它 们 很 挥
作用形成载体粒子 。同时, 放射性裂变产物在 载体粒子表 面发生冷凝。在冷凝过程 中, 由于
( £ )=∑ ∑ H x( A ) 3 c e 一 r ;() p t
l0 l I
n—l
n
C =—; L一 h 『上 。 () 4
2 分 凝 模 型
基 于放 射 性 衰 变 链 的 B t a a m n公 式 和 e
Fel g的核素 半经 验 分凝 模 型 , 建 核爆 粒 子 rin i 构
1 5 0 K时 , 7 估计 大 多数难熔 元素 , 锫和稀 土 , 如
载体 粒 子 上 分 布 不 同, 种 现 象 称 为 “ 这 分
凝” …。分凝效应使得裂 变产物在放射性粒子 上 的分布不均匀 , 使得核爆放射性沾染场 的辐 射性质在不 同地点和不同时间不同。在给定 的 核爆袭击情景下, 评估预期伤亡时 , 分凝效应有 很 大 的影 响 。
3 33 . 3 3× 1 2 0-
于气 化 土 壤 , 那么 将 有 成千 上 万 吨土 壤卷 入火
1 分凝效应
1 1 分 凝机 理 .
在 核爆 粒子形 成 过程 中并 不是 所有 的放 射 性 核 素 同时 冷凝 到粒 子 上 , 些非 常难 熔 的元 一
球 中。气化了的土壤分压将是裂变产物分压的
一
百多万倍 。较难熔的成分, 如二氧化硅 , 同沸
放射性衰变链 的 Bt a 公式构建 核爆 粒子放射性活 度分布的计算模型 。计算 了给定载体物质 和粒子 am n e 尺度分布条件下的粒子放射性活度分 布 , 果表 明挥发 性核 索呈 面分布 特征 , 结 难熔性 核素呈 体分 布特 征 。分析了粒子尺度 、 固温度 和凝 固时间对 放射性活 度分布的影响程度 , 凝 结果表明粒子尺度分 布和凝
怒 ( 扎 (
() 7
粒子质量份额
3 33 . 3 3× 1 - 0 2
3 33 . 3 3× 1 2 0-
其中:
表 3 总 B活度和质量链 8 ,5的分 布 99
粒子尺度/ m
3. 6 ×1 3—1 6 3 0×1 3 6 89 0- 2 0
.
总 B活度/ q B
5 6× 1 1 .1 02
1 O ×1 2 .0 0 ̄ 1 O ×1 2 .0 0o
1 01×1 2 . 0 ̄
I 0 3 0×1 3—7 5 4 6×1 ‘ . 2 0- 4 0一
.
3 3 3 3×1 2 .3 0- 3. 3 ×1 2 3 33 0- 3. 3 ×1 2 3 33 0- 3. 3 ×1 2 3 33 0-
郭立军 陈琳2郑毅。 , ,
(. 1海军 陆战学 院 , 广州 5 0 3 ;. 14 0 2 防化 指挥工程学 院 , 北京 12 0 ;. 0 2 53 防化研究 院 , 北京 1 2 5 2 0 0)
摘 要 : 了研究核爆粒子 中放射性核素 分凝 的影 响因素和规律 , 为 基于 Fe i r l g的半经验分凝模 型和 in
核素 不 同的挥 发性 而导致 在不 同尺 度 和性质 的
快衰变为其他元素 , 否则不冷凝到粒子上 。在 粒子未凝固之前冷凝 的核素 , 近似均匀分布在 粒子体 内。粒子凝 固之后冷凝 的核素, 主要分 布在粒子表面。 从 物 态相平 衡 的观点 来 看 , 凝 取 决 于冷 冷 凝 相物 质 的沸点 和蒸气 相 中的分 压 。火球 热力 学温度下降到 3 0 400K时 , 0~ 0 0 其总压力大 致与周 围大气 平衡 J 。以 l00 k 爆 炸为 0 t 例, 估计此时裂变产物分压 为 l “大气压。在 0
n ( 一 A)
式中: () £是链 中时刻 t n个核素的原 第 子数 ; 是链 中 0时刻第 ( n 个元素 的原 i ) 子数 ; A 是链 中第 七 k ) ( sn 个核 素的衰 变常 数。
t 时刻 , 素 的放 射性 活度 为 : 核 A ()= A () t t。
.
4. 7×1 2 5 01 4. 4×1 1 6 02 4 6 ×1 2 .9 01 4 7 ×1 2 .4 01 4 7 ×1 2 .9 0l 4. 3 ×1 2 8 01 4. 8 ×1 2 8 01 4. 2× 1 1 9 02 4 9 ×1 1 .6 02 5 01×1 2 . 01 5 0 ×l 2 .6 0l 5. O ×1 2 1 01
.
∑
^= 1
系数相关性相符 , Fe i 半径指数模型 的基 是 r lg i n 本依据。 将粒子群根据尺度大小分为若干组, : 令
=
D取值 10 m。那 么 第 后组 粒 子 中质 量 0 链i 的原 子数 为 : N ( ̄ id)= F( ) i (O 1)
k 尺度组粒子的平均直径 : 尺度组 的粒子数 ;
3 模 拟计算
给定粒子凝固时刻 的时间和温度 、 射性 放 粒子的总量和尺度分布参数 , 模型可以计算 总 的放射性活度的尺度分布和单个质量链的尺度 分布 。表 2 列出计算条件。表 3列出总的放射 性活度分布 和质量链 8 , 9 质量链 9 5的分 布。 图 2和 图 3 示 为总 的 p比活 度 和质 量链 8 、 所 9 9 5的比活度分布。
由B 她m l 蛆议程计算指定时刻t 各衰
变链的放射性活ቤተ መጻሕፍቲ ባይዱ
1个 沉 降 样 本 中 给 定 核 素 i的 裂 变 数
为 : = a 。 / () 1
由各链 的放射性和 分布份额计算 时囊 t I 总的 放 射性分 布呀指 定衰变 链的分布
式 中 , 样本 中核 素 i a是 的原 子数 , 是 核 素 的裂 变 产 额 。那 么 核 素 i 对 于核 素 的 相 分凝 系数 定义 为 :
E = —— i —L~
h=l
;
() 8
(/ 一 )  ̄ 1次幂, b= 令 i
对数 相关 性为 :
lr 9 n i5=
,
, 推出分凝系数的
∑ ( ) )
厶 ( ): 一; () 9
( 1一b)n8 +a。 ilr 9 i () 6 这 与 Feig从 放 化 分析 得 到 的核 素 分 凝 rin l
其 中的一 种 或几 种 核 素 ( 括 Ⅳ1可 以直 包 ) 接 由裂 变产 生 。除 了 Ⅳ , m 每一 种 核素 按 照衰 变 律 ep(一Af 衰 变 。 除 了 Ⅳ 每 种 核 素 按 照 x i ) ,
时刻衰变链 中难熔性核素所 占份额 ( ) 决定
73 7
整个 链 的 冷 凝 行 为 。 Feig 提 出 给 定 质 量 rin l 链 i的 放 射 性 活 度 分 布 正 比 于 粒 子 直 径 的
r J= 。 () 2
图 1 计 算 程 序 流 程
ep A t增长。衰变链中每种核素的生长 x (一 ) 和 衰减历 程 可 以用 B t n8方程 表示 : ae [ ma
通常将产额 大、 半衰期 长 的曲 r z 分 s 和 r 别作为典型挥发性和难 熔性的参考核素 , 分析 核素相 对参 考核 素 的分 凝 系数 关 系。
核素 的分凝 程度 用 相对 于不 分凝 核素 的分
。
凝系数 r 表示 。所谓不分凝 核素是指一些难 熔 的 核素 , 们 的母 体 尽 管是 一 些 易 熔 或 挥发 它 性核素 , 但是半衰期很短 , 很快就会衰变成为难
熔 核 素 , 为 它 们 不 发 生 分 凝 。 例 如 质 量 链 认
2 2 核素 分 凝 .
放射性活度分 布的计 算模 型, 计算流程如 图 1
所示 。
2 1 核 素衰 变 .
放射 性核 素 在 冷 凝 的 同时 发 生 着 衰 变 , 使
其浓度不断变化。计算分凝效应必须考虑每种 核素的衰变和生长 , 以确定每种核素在冷凝时 刻的浓度 。单个放射 性核 素按照其 半衰期 衰 减, 而裂变产物往往形成许多 由相同质量数核 素组 成 的衰变 链 , 如某 个 衰 变链 可 以表示 为 :
收 稿 日期 :02— 1— 0 2 1 0 2
点差不多、 分压低得多的裂变产物相 比, 冷凝得
早得多。
作者简 介 : 郭立军 , , 男 黑龙 江拜 泉人 , 讲师 , 士研 硕
究生, 研究 方向为核辐射探测 。
77 2
表 1 裂变产物 中各元素 的冷凝性质 ‘
1 2 分凝程 度 .
固温度的影响较大 , 子的放射性活 度分布与爆 炸当量 、 粒 比高 以及爆炸地点 的地理条件密切相关。
关键词 : ; 核爆 放射性核素 ; 分凝
中图分类号 : T 1 L8 文献标 志码 : A 文章编号 : 0 5 -94 2 1 ) 70 7 -6 2 80 3 (0 2 0 - 20 7
本 文应 用数 值模 拟 的方法 计算 核爆 粒 子 的
放射性活度分布 , 分析放 射性核素分凝作用 的
影 响 因素和 规律 。
平 衡蒸 气 压 在 l 大 气 压 以下 , O 它们 会 发 生冷
凝 。而卤素、 惰性气体等元素沸点较低 , 平衡蒸
汽压 在 l 大气 压 以上, 们 不 会 发 生 冷 0 它 凝 。各种裂变质量链 中元素 的沸点见表 l 】 。 对于地爆而言, 如果有百分之几 的爆炸能量用
第3 2卷
21 0 2年
第 7期
7月
核 电子学 与探 测技 术
Nu la e to i s& De e t n Te h oo y c e rEl cr n c tc i c n lg o
V0. 2 No 7 13 .
J l. 2 1 uy 0 2
核 爆 放 射 性 核 素 的分 凝研 究
4. 8×1 }1 4 (2
质量链 8 / q 9B
2. 3 x1 2 O 2 3 ×1 9 .2 01 2 41×1 9 . 01 2. 9 ×1 9 4 01 2. 8 ×l 5 09
质量链 9/ q 5 B
99 09 . 7 x1
1 6 3 0 ×1 3—1 0 3 0×1 3 .2 0- 2 0
Nl N2_ +N3 … N
m
() 5
Fei n 以土 壤粒 子 凝 固的 温度为 rlg t in 界, 将沸点低于该温度的核素归为挥发性核素 ,
沸点高于该温度的核素归为难熔性核素 。挥发 性核素仅仅冷凝在粒子表面 , 而难熔性 核素均
匀分 布在 粒 子体 内。任意 裂变 产物 衰变 质量链 的组 分 数 量 都 在 随 时 间变 化 , 土壤 凝 固温度 在