和1990年ICRP推荐接轨的内照射剂量估算方法的研究及其应用软件的开发.pdfVIP

全国职业照射个人剂量监测与评价学术研讨会 方法的研究及其应用软件的开发 陈宝维 马如维 毛 永 (中国辐射防护研究院，太原，030006) 摘要／本文介绍的是用于进行个人内照射剂量估算的新一代计算机应用软件系统。该系统参照 【cRP第56、66、67、68、69、71、78号等一系列出版物提供的代谢模型和参数，采用了适于在微机上实现的 矩阵算法，使用vc++编制而成，可在Win95或更高版本的微软操作系统下运行。沱可以将个人内照 射监测量换算为摄人量、有效剂量和其它评价量，适用于职业受照人员的常规个人内照射剂最估算。和 先前的IND0相比，本套软件采用了新的模型和参数，可以直接解任何不循环的和循环的库模型。这使 得它可用于与1990年lcRP推荐和In年的国际电离辐射防护和辐射源安全的基本安全标准接轨的个 人内照射剂量管理和估算。此外，本套软件在视窗操作系统下运行，图形界面使得它给用户的感觉更加 舒服、使用更加方便。 关键词 内照射剂量估算摄入量计算机程序 1引言 年剂量和年摄人量是目前辐射防护中主要的评价量。然而，通常个人内照射监测的测定量，主 要是全身、体(Syste血c)及器官或组织的滞留量和经尿、粪等途径的排出量。因此，将监测量换算为 评价量是辐射防护中的一个重要问题，各国研究人员为此做了大量的工作，并得出了多种由监测量 推算摄人量的方法。由于放射性核素在人体的滞留和排泄行为与其化合物类型、摄入途径、气溶胶 粒径(AMAD)以及受影响的器官或组织等因素密切相关，故这类计算是十分繁复而费时的。目前， 轨的内照射剂量估算软件删Do①。但由于内照射剂量学的不断完善和发展，INDO不能适应当前 参数，采用了适于在微机上使用的矩阵算法②及能够充分利用C／C++语言的快速特点及32位微 机性能的vc++开发工具，使得该软件界面友好、计算快速、结果精确。该软件基本保留了IND0 的功能，可以在数秒钟内，将个人内照射监测量换算为摄入量和有效剂量，可用于与1990年IcRP 推荐和1997年的国际电离辐射防护和辐射源安全的基本安全标准接轨的职业照射个人内照射剂 量管理和估算。本文主要介绍栅2000的算法、功能、特点和设计思想。 2算法的选择及矩阵算法 编制与1990年lcRP推荐接轨的内照射剂量估算软件，关键是组建并求解库模型。物质在模 型中任何一个库中含量的变化由两个相互竞争的过程来决定：由进入途径而增加和由离开途径而 减少。如果％表示某一时刻库i的物质含量，kq表示库i向库j的转移速率，则可以写出描述这个 系统的一阶线性微分方程组： 瓮=*|吨瓠 ，≠i J≠I 204· 全国职业照射个人剂量监测与评价学术研讨会 式中wf和∞分别是库i和j的含量，它们是时间的函数。 我们软件的目的，实际上是解此方程组。 一般来说，微分方程组(1)是很难求得解析解的，除非全身模型本身就是或可变化为简单的线 性库链模型，因此通常采用数值解的方法解得托(￡)。数值解的方法很多，不同作者采用了不同的 方法，它取决于使用的计算工具，要求的精度，要求的计算时间，以及所解的模型本身。我们在微机 上运算，采用过积分一步变步长的Gill方法，叠代方法和矩阵方法，试验表明在计算的复杂程度、精 度和用时等方面，各有各的优缺点。综合考虑结果，最终选择了矩阵算法。下面简述以下矩阵算 法： 定义一个n阶矩阵[A]，它的各个元素。。为： ％=颤，i=1至n，j=1至n且i≠』 (2) o“=一∑i=1至n (3) 。 ，=1 J■I 将(2)式和(3)式代入方程组(1)，可得： + 掣=耋。卉(‘)怕m(#)㈦至n(4)