放射性同位素衰变校正程序设计
作者：程昕 唐明华
单位：程昕 唐明华 北京放射医学研究所 100850
关键词：
放射性同位素衰变校正程序设计
辐射监测中，常用已知活度的放射源刻度仪器效率、检定仪器性能，但计算校正时间较麻烦。因此，为能迅速而准确地对核素进行衰变校正，需建立相应的计算程序。
一、方法
1.放射源类型
可分为已有并入库存档的标准源的新购待测的标准源。它们又可分为单能源和混合源。
2.校正日期计算方法
, http://www.100md.com
放射性同位素的衰变校正计算公式为：
式中：I是校正后活度(Bq);I0是t为零时的活度(Bq)；T1/2是半衰期(天)；t是标定至测量的间隔时间(天)。对于标准源，I0、T1/2、标定时间及测量时间都是已 知的。
t的计算方法(用于T1/2较长的核素)如下：
(1)标定或测量时间的月份为1、2月时：
tB(或tc)=365×Y+D+31(M-1)
, 百拇医药
+INT[(Y-1)/4]-INT
{(3/4)×INT
[(Y-1)/100+1]}
(2)标定或测量时间的月份为3～12月时：
tB(或tc)=365×Y+D+31(M-1)
+INT(0.4×M+2.3)
+INT(Y/4)-INT
[(3/4)×INT(Y/100+1)]
t=tc-tB
, 百拇医药
式中：Y、M、D分别为年、月、日；INT是取整函数。根据标定和测量时间的月份选择相应的计算公式，得出测量日期的天数(tc)和标定日期的天数(tB)，二者之差即为t。
3.程序设计
本程序用BASIC语言编写。主要采用“人机对话”的菜单工作方式。程序运行时，屏幕显示 菜单，提示操作者选择放射源类型，根据所输入的代码，进行下述工作：
(1)选用库存单能源时，只需选择核素的标定日期，程序就自动赋值、计算、打印；
(2)选用待测单能源时，程序提示输入核素名称、标定日期、标定活度及半衰期后，即能自 动计算活度并打印；
(3)选用库存混合源时，只要输入混合源的代码，程序就可自动赋值、计算、打印；
, 百拇医药
(4)选用待测混合源时，首先输入混合源核素的个数，然后程序提示依次输入每个核素的名 称、标定日期、标定活度及半衰期，程序将自动记录、计算结果并打印输出。
二、结果和讨论
结果以表格方式输出(见表1、表2)
本程序的特点是采用人机对话的菜单工作方式和纠错两个方面。菜单显示、代码输入 ，便于用户使用，即使是不懂程序的人也不会感到吃力；当不慎输入错误字符时，程序可自动识别，让用户重新选择。若想对放射性核素进行全面、系统的管理，还可在此基础上建立数据库，完善后的库即可推广使用。
表1 单能源γ核素半衰期校正前后活 度
核素名称
半衰期
, 百拇医药
(年)
标定日期
测量日期
活度(Bq)
标定
测量
137Cs
30.174
1981.12.1
1997.12.12
68 450
47 356
, 百拇医药
表2 混合源γ核素半衰期校正 前后活度
核素名称
(No.1)
半衰期
标定日期
测量日期
活度(Bq)*
标定
测量
241Am
433a
1991.11.1
, 百拇医药
1997.7.1
5700
5649
109Cd
453d
1992.5.23
1997.7.1
4000
231
137Cs
30.174a
1992.9.6
, 百拇医药
1997.7.1
4700
4208
54Mn
312.5d
1992.4.16
1997.7.1
4770
70
57Co
270d
1992.6.16
, http://www.100md.com
1997.7.1
5810
52
65Zn
244d
1992.6.22
1997.7.1
3460
19
总活度
10 227(Bq)
(收稿：1998-01-04 修回：1998-06-09), http://www.100md.com
单位：程昕 唐明华 北京放射医学研究所 100850
关键词：
放射性同位素衰变校正程序设计
辐射监测中，常用已知活度的放射源刻度仪器效率、检定仪器性能，但计算校正时间较麻烦。因此，为能迅速而准确地对核素进行衰变校正，需建立相应的计算程序。
一、方法
1.放射源类型
可分为已有并入库存档的标准源的新购待测的标准源。它们又可分为单能源和混合源。
2.校正日期计算方法
, http://www.100md.com
放射性同位素的衰变校正计算公式为：
式中：I是校正后活度(Bq);I0是t为零时的活度(Bq)；T1/2是半衰期(天)；t是标定至测量的间隔时间(天)。对于标准源，I0、T1/2、标定时间及测量时间都是已 知的。
t的计算方法(用于T1/2较长的核素)如下：
(1)标定或测量时间的月份为1、2月时：
tB(或tc)=365×Y+D+31(M-1)
, 百拇医药
+INT[(Y-1)/4]-INT
{(3/4)×INT
[(Y-1)/100+1]}
(2)标定或测量时间的月份为3～12月时：
tB(或tc)=365×Y+D+31(M-1)
+INT(0.4×M+2.3)
+INT(Y/4)-INT
[(3/4)×INT(Y/100+1)]
t=tc-tB
, 百拇医药
式中：Y、M、D分别为年、月、日；INT是取整函数。根据标定和测量时间的月份选择相应的计算公式，得出测量日期的天数(tc)和标定日期的天数(tB)，二者之差即为t。
3.程序设计
本程序用BASIC语言编写。主要采用“人机对话”的菜单工作方式。程序运行时，屏幕显示 菜单，提示操作者选择放射源类型，根据所输入的代码，进行下述工作：
(1)选用库存单能源时，只需选择核素的标定日期，程序就自动赋值、计算、打印；
(2)选用待测单能源时，程序提示输入核素名称、标定日期、标定活度及半衰期后，即能自 动计算活度并打印；
(3)选用库存混合源时，只要输入混合源的代码，程序就可自动赋值、计算、打印；
, 百拇医药
(4)选用待测混合源时，首先输入混合源核素的个数，然后程序提示依次输入每个核素的名 称、标定日期、标定活度及半衰期，程序将自动记录、计算结果并打印输出。
二、结果和讨论
结果以表格方式输出(见表1、表2)
本程序的特点是采用人机对话的菜单工作方式和纠错两个方面。菜单显示、代码输入 ，便于用户使用，即使是不懂程序的人也不会感到吃力；当不慎输入错误字符时，程序可自动识别，让用户重新选择。若想对放射性核素进行全面、系统的管理，还可在此基础上建立数据库，完善后的库即可推广使用。
表1 单能源γ核素半衰期校正前后活 度
核素名称
半衰期
, 百拇医药
(年)
标定日期
测量日期
活度(Bq)
标定
测量
137Cs
30.174
1981.12.1
1997.12.12
68 450
47 356
, 百拇医药
表2 混合源γ核素半衰期校正 前后活度
核素名称
(No.1)
半衰期
标定日期
测量日期
活度(Bq)*
标定
测量
241Am
433a
1991.11.1
, 百拇医药
1997.7.1
5700
5649
109Cd
453d
1992.5.23
1997.7.1
4000
231
137Cs
30.174a
1992.9.6
, 百拇医药
1997.7.1
4700
4208
54Mn
312.5d
1992.4.16
1997.7.1
4770
70
57Co
270d
1992.6.16
, http://www.100md.com
1997.7.1
5810
52
65Zn
244d
1992.6.22
1997.7.1
3460
19
总活度
10 227(Bq)
(收稿：1998-01-04 修回：1998-06-09), http://www.100md.com