核技术应用中的蒙卡计算问题

第四部分：例题

EGS4软件包给出了七个例题，供用户使用 与程序检验。
EGS4软件包的使用方法 使用EGS4软件包，需要编写用户程序。其中包 括主程序（及其所需要的若干子程序），和EGS4 程序包所要求的两个子程序：HOWFAR子程序和 AUSGAB子程序。

3.GEANT4程序包


Geant4是由欧洲核子中心主导开发的一套用于 Monte Carlo模拟程序包. 来自于美国、俄罗斯、 日本、加拿大等国家10多个实验室的100多名科学 家参与了Geant4程序的研制工作. 它包括了实验装置的描述,粒子在材料和磁场 中的运动,以及粒子与物质相互作用的物理过程模 型等一整套工具包.
六．蒙特卡罗方法在粒子辐射效应、 抗辐射加固中的应用

在核聚变堆、空间卫星及飞行器中，其关键 部分的元器件、晶体受到高能粒子（中子、光子、 电子、质子等）辐射是造成各种元器件失灵、失 控乃至事故的主要原因。因此，计算各种元器件 材料的辐射损伤，对于材料的辐照筛选，设施的 抗粒子辐射加固和安全运行有重要的意义。
一．引言


在核技术应用领域的研究中，由于受实验条件的 限制，有很多问题需要借助于理论计算来完成。对于 几何结构复杂且粒子与核的反应机制复杂的问题，一 般数值方法无法求解，而蒙特卡罗方法能够比较逼真 地描述事物的特点及物理实验过程，解决数值方法难 以解决的问题，因而该方法广泛地应用于核技术应用 研究中。 本文重点介绍蒙特卡罗方法 [1] 在几个核技术领域中的 应用，以及所使用的蒙特卡罗模拟方法；蒙特卡罗方 法的软件；最后，给出三个应用实例。
六．蒙特卡罗方法在粒子辐射效应、 抗辐射加固中的应用

高能粒子辐射到晶体后，通过与核反应，产生 多种粒子（如中子、光子、电子、质子、α粒子等） 和反冲核，而其中的每一种粒子、反冲核通过各种 反应道又可生成各种粒子与核。这种复杂的偶合输 运问题不仅需要考虑粒子（不带电与带电）与核的 反应，还需要考虑原子碰撞。
实例一
1995年10月在莫斯科举行的“第七届国际聚 变堆材料会议”上报告了该种实验设计、计算方 法与计算结果，引起了较大的反响，认为这一成 果在国际上处于先进位置，实验、计算与计算方 法首次使用。

实例二、 高纯Ge探测器在不同几何上测量的 “有效作用深度”的研究与计算
该项工作与实验工作一起，获得了部级科技 进步二等奖，以及国家级科学进步三等奖。
光子-电子-核-核的多分支方法
e
r r e e e r r e e
e
4.粒子—粒子—核—核偶合输运的多分支方法
该方法实现多种粒子的偶合输运问题。粒 子—粒子偶合输运过程与光子—电子偶合输运 的多分支方法类似，由于粒子的种类不同，需 要加入描述粒子类型的存储单元。这个方法是 在粒子—粒子偶合输运的基础上嵌入核—核偶 合输运过程。
二．蒙特卡罗方法简介


蒙特卡罗方法又称计算机模拟方法、随机 抽样技巧或统计实验方法。半个多世纪以来， 由于科学技术的发展和计算机的出现和发展， 这种方法作为一种独立的方法被提出来，并首 先在核武器的实验与研制中得到了应用。 蒙特卡罗方法是一种计算方法，但与一般 数值计算方法有很大区别，它能够比较逼真地 描述事物的特点及物理实验过程，解决数值方 法难以解决的问题，因而该方法的应用领域是 十分广泛的。
七．蒙特卡罗模拟方法
1. 对于单个粒子（中子或光子），其模拟方法 通常有：直接模拟法，加权法，统计估计法、 指数变换法等。
2.两种粒子（比如：光子与电子）偶合输运问题


对于两种粒子（比如：光子与电子）偶合输 运问题，使用“字典编辑式”的多分支方法。 该方法是：在光子的输运过程中，一旦有光子、 电子产生，将光子存储，跟踪电子。模拟的电子 历史结束后，依照后进先出的原则取出光子，进 行模拟。直至所有的次级光子和电子的历史全部 结束，一个由源出发的光子的历史结束 。
三．蒙特卡罗方法在探测器系列参数 计算中的应用

为了把实验测得的多能光子脉冲高度谱分解成 单色光子脉冲高度谱，需要对探测器进行刻度。虽 然可以通过实验的方法用光子源测得其能量响应， 但是这样的放射性核素是有限的，特别是高能光子 源更难以得到；对于低能光子源（小于2Mev）,实 验测得的结果偏高。因此，要得到对任意能量的光 子的响应函数，需要借助于理论计算来完成。
五．蒙特卡罗方法在核辐射医学计算 中的应用

关于体内剂量场分布，主要有两个途径获得： 一是基于各种测试数据和经验公式，使用近似插 值方法；另一种方法是利用核物理知识和模型， 进行理论计算。射线在人体内的反应机制是光子 和电子偶合输运过程，非常复杂，一般数值方法 难于求解。蒙特卡罗方法是解决该类问题十分有 效的方法。
2.1 蒙特卡罗方法的计算原理

当所求问题的解为某个事件的概率，或者 是某个随机变量的数学期望，或者是与概率、 数学期望有关的量时，通过随机实验的方法， 得到该事件发生的频率，或者该随机变量若干 个具体观察值的算术平均值，以频率代替概率， 或者算术平均值代替期望值得到问题的解。
2.2 蒙特卡罗方法的特点
EGS4软件包括四个部分



第一部分：PEGS4一截面数据处理程序 它的功能是给出系统介质所需的截面及介质中各 种元素的截面与各种反应的分支比。截面数据是 使用分段线性拟合方法产生的。 使用该程序时，首先需要按照规格要求填写输入 卡片（输入文件），文件名为：USERINP。DAT， 然后运行执行程序PEGS4。计算结果在输出文件 FOR012.DAT 中。运行过程在文件EGSOUT.DAT 中可以看到。
三．蒙特卡罗方法在探测器系列参数 计算中的应用

光子源有点源、面源和体源。
由于光子与探测器中的物质发生反应后，产 生次级光子和电子，而电子和正电子在输运过程 中又可产生光子。因此，这一类计算问题需要考 虑光子和电子的偶合输运过程。

四．蒙特卡罗方法在核辐射防护计算 中的应用

一个核设施（例如：加速器，实验室等）的 建立，辐射安全问题是必需考虑的问题。粒子 （中子或光子）经过屏蔽层后，有多少泄露，屏 蔽层的厚度为多少才能满足防护安全的要求，以 及如何设计核设施中的结构，即所谓最佳设计方 案问题，需要用蒙特卡罗方法模拟。这样可以节 省大量的实验时间和经费。
实例一

由计算结果可看出：当辐照温度低于700°C 时，热扩散系数随损伤剂量的增加而减小，但对 于损伤剂量为0.3dpa时，其热扩散系数比前者大， 原因不仅是由于辐照温度高（770°C），而且由 于高温弛豫时间长（13.3min）。当辐照温度 800°C时，损伤剂量为0.1dpa时，其扩散系数却 比第一、二种情况减少了许多，唯一的不同是高 温弛豫时间短（6.3min），因此，可以得到一个 重要结论：当高温弛豫时间大于6min时，能够减 少对晶体的损伤。

实例三、 I-125体内放射源剂量场分布计算 1. 物理模型
实例三


放射源装置最内层是圆柱形的银棒，中间层是 空气层，最外层是圆柱形的钛层。 γ射线均匀地分布在钛层表面，方向为各向同性。 能量服从如下能谱：
22.1 0.0058 24.9 0.0345 25.5 0.0107 27.2 0.2524 27.4 0.4641 31.0 0.1490 31.7 0.0342 35.5 0.0493

粒子-粒子-核-核偶合输运的分支方法
+ + + P D ‫ס‬
‫ס‬ ‫ס‬
He3 ∆ • • α ∆ t ∆
•
八．几个常用的蒙特卡罗模拟软件

MCNP软件包
EGS4 软件包 GEANT4程序包


1.MCNP软件包 MCNP的全名是：Monte Carlo Neutron and Photo Transport Code. 它是由美国Los Alamos 实验室研制的一个大型的多功能的蒙特卡罗程序 包。可用于计算中子、中子光子、光子电子及其 组合的输运问题，以及临界（包括次临界和超临 界）系统的本征值计算问题。是目前国内外普遍 使用的程序包。
光子-电子的偶合输运多分支方法
e e e e r r r e e e r r
e
r r
r
e
e
r
3. 光子—电子—核—核偶合输运的多分支方法

该方法是：在光子—电子偶合输运的多分支方 法的基础上，嵌入核—核偶合输运的多分支方法。 即：在光子—电子偶合输运的过程中，如有反冲 核产生，即转向核—核偶合输运。当核—核偶合 输运结束后，返回到光子—电子偶合输运的过程 中去。

蒙特卡罗方法的主要应用范围包括：粒子 输运问题，统计物理。在粒子输运问题中的应 用，主要包括：实验核物理、反应堆物理以及 高能物理等方面，以及真空技术，激光技术、 医学、生物、探矿等方面。
三．蒙特卡罗方法在探测器系列参数 计算中的应用

探测器系列参数的计算包括探测器对光子的 探测效率、能量沉积谱、响应函数、全能蜂以及 与它们有关的量。
2.2 蒙特卡罗方法的特点
④
⑤
⑥
蒙特卡罗方法具有同时计算多个方案与多个 未知量的能力。 蒙特卡罗方法的计算误差容易确定，它的程 序结构简单、且易于实现。 同时，蒙特卡罗方法具有以下的缺点：收敛 速度慢，误差具有概率性，在粒子输运问题 中，计算结果与系统大小有关。
2.3 蒙特卡罗方法的主要应用范围
五．蒙特卡罗方法在核辐射医学计算 中的应用

核辐射医学在治疗疑难病症（如癌症）及重 症病（如心脏病）有着非常重要的作用，在放射 性治疗（外照射或内照射）过程中，病体所受剂 量及其分布是治疗过程中必须考虑的问题，剂量 小影响治疗效果，剂量大对病体造成伤害，因此， 当射线进入人体后各部分受到的剂量必须要进行 理论计算。
①
②
③
蒙特卡罗方法能够比较逼真地描述具有随机 性质的事物的特点及物理实验过程。它可以 部分代替物理实验，甚至可以得到物理实验 难以得到的结果。 用蒙特卡罗方法解决实际问题，可以直接从 实际问题出发，而不必从物理方程或数学表 达式出发。它有直观形象的特点。 蒙特卡罗方法受几何条件限制较小。它的收 敛速度与问题的维数无关。
第二部分：EGS4软件包

该软件包由两部分组成：SHOWER和HATCH。 SHOWER是模拟光子一电子偶合输运过程的子程 序，它在运行过程中又调用其他子程序和两个用 户子程序HOWFAR（几何输运程序）和AUSGAB （记录输出程序）；HATCH是输入介质数据子程 序。
第三部分：PRO—公共区文件 包括八个公共区，可根据需要使用其中有关 部分。
核技术应用研究中的 蒙特卡罗计算问题
报告人：许淑艳 （中国原子能科学研究院）
目录





一．引言 二．蒙特卡罗方法简介 三．蒙特卡罗方法在探测器系列参数计算中的应用 四．蒙特卡罗方法在核辐射防护计算中的应用 五．蒙特卡罗方法在核辐射医学计算中的应用 六．蒙特卡罗方法在粒子辐射效应、抗辐射加固中的应用 七．蒙特卡罗模拟方法 八．常用蒙特卡罗模拟软件 九．应用实例
3.GEANT4程序包

它的一个突出特点是它包含了非常丰富的物理 模型,并且能够在非常大的能量范围内处理粒子与物 质的相互作用.正是由于它的粒子种类多，物理模型 全，能量范围大这些特点，使得他的应用领域越来 越广泛，包括高能物理，核试验，加速器，医学， 生物科学，辐射防护等多个领域。
Βιβλιοθήκη Baidu
九．应用实例

1.MCNP软件包

MCNP—3B是1989年公开发表的版本，目前使用 的版本大多为MCNP—4A，MCNP--4B，MCNP— 4C版本。
使用MCNP软件包需要填写输入卡片：INP。计 算结果放在输出文件OUTP中。

2．EGS4 软件包

EGS（Electro一Gamma Shower）软件包是光 子一电子在任意几何中偶和输运的蒙特卡罗摸拟 通用软件包。由美国斯坦福直线加速器中心编制， EGS4软件包是较新版本。



实例一、聚变堆第一壁材料SiC中电子辐射损 伤 实例二、 高纯Ge探测器在不同几何上测量的 “有效作用深度”的研究与计算 实例三、 I-125体内放射源剂量场分布计算
实例一、聚变堆第一壁材料SiC中电子辐射 损失

计算了不同能量的入射电子在不同剂量下的 位移损伤强度，得到了入射能量为14Mev的电子， 在不同的辐照温度，不同的损伤剂量及不同的高 温弛豫时间，热扩散系数随测量温度的变化情况。