蒙特卡罗模拟小孔问题中一种提高探测器计数效率的方法.pdfVIP

第一篇蒙特卡罗方法 蒙特卡罗模拟小孔问题中一种提高探测器 计数效率的方法 潘流俊，邓力，许海燕，王瑞宏 (北京应用物理与计算数学研究所，北京]oooss) [铺要】 在使用指向概率方法计算定向发射强度、石油测井和物探等一些问题中，常常遇到只有当粒子在系统 中的一个小区域发生碰撞时才对探测器有贡献的情况。为了获得更多粒子在特定小区域内的碰撞信息，引入“虚 粒子”的概念．当粒子从源出发或在系统中发生碰撞时，定义并抽取指向特定小区域的虚粒子，计算虚粒子无碰 撞到达特定小区域边界的概率及虚粒子到达小区域边界的位置和时间。到达特定小区域边界的虚粒子在小区域内 输运、碰撞并用指向概率方法对探测器计数。计算一个简单的例子结果表明，这一方法是有效的。 ￡关键词】 蒙特卡罗方法； 指向概率： 定向发射强度： 石油测井； 虚粒子 【中图分类号】0242 【文献标识硼 A 0引言 在使j{j蒙特卡罗方法中的指向概率方法计算实际问题时，尤其在有关探测散射粒子信息的问题中，经 常会遇到只有当粒子在系统中的一个小区域发生碰撞时才对探测器有贡献的情况。 例如，在图1所示的，，．散射岩性密度测井仪问题【3】巾，从源位置出发的光子经过岩层散射到达探测器。 在使用指向概率方法计算时，只有当光子在探测器准直孔所对应的相对很小的区域内发生碰撞才会对探测 器的计数有贡献。 在计算定向发射强度【I】【21时也常常用到指向概率方法。如图2所示，要计算沿着垂直于计数平面穿过 计数平面上某些网格的粒子数。通常的计算办法【l】是，粒子从源出发或者发生碰撞时计算粒子沿着垂直于 计数平面单位方向飞行的概率，计算沿该方向飞行的粒子无碰撞到达计数平面的概率并计数。当计数网格 相对较小时，只有小区域内出发的粒子或发生的碰撞才对计数有贡献。例如，要计算图2中计数平面上一 点附近小圆内的定向发射强度，只有以该小圆为底，以垂直计数平面方向为对称轴的柱形区域内出发的粒 子或发生的碰撞才对计数有贡献。 图I岩性密度测井示意图 图2定向发射强度计算示意图 ofwell directionalemission Diagramofdesignated strength Fig．1Diagramlogging Fig．2 在用蒙特卡罗方法模拟康普顿散射检测【41，小孔成像及其他一些带有准直孔的模犁时也会遇到类似的 l作者俺介I潘漉俊(1979一)，男，Miei．助理研究员。从事蒙特卡罗方法及其应用研究，北京8009信箱16分箱．北京100088· 114 第十届全国蒙特卡罗方法及其应用学术会论文 问题。这类问题往往需要计算详细的能量谱、时问谱或者平面上按网格划分的分布，计算量大．本文针对 这类问题提出一个改善计算效率的方法，通过一个简单的例子说明了这一途径并进行了数值实验。 1计算方法 为了叙述方便，用图3所示的一个例子来说明方法。左边是系统，中间是理想屏蔽体，屏蔽体上开一 小孔，在小孔后面放置一个点探测器。显然，在使用指向概率方法时，只有锥形区域Q内出发的粒子或发 生的碰撞才对探测器计数有贡献。 僻蔽体 U 闭3模型示意图 ofmodel Fig．3Diagram 区域Q相对_丁：整个系统是一个小区域，到达Q内的粒子数不多，碰撞也就不多，探测器计数效率不 高。为了更好地模拟散射粒子对探测器的贡献，在每次源抽样或粒子发生碰撞时抽取能到达区域Q的虚粒 予，虚粒子在区域Q内发生碰撞并计数。为了区分，把跟踪粒子输运过程的粒子叫做模拟粒子，对模拟粒 予不再做计数。每个样本跟踪的具体算法如下：