评MCNP程序在核测井模拟中的应用研究.pdfVIP

盯 评MCNP程序在核测井模拟中的应用 裴鹿成 (中国原子能科学研究院) 摘要 MCNP是世界上最著名的解决中子一光r一电了辐射输运问题的超级蒙特卡罗方法程序，其几何能 力与截面能力之强，应用范围之广，含蒙特 I=罗技巧之多，是任何其他蒙特P罗方法程序所无法比拟的。 MCNP程序广泛应用士许多科学与工业领域，核测井模拟是其最重要的应用领域之一 对于有的核测井问 题用MCNP程序可以很好地直接得到解决，有的则需要用户在MCNP程序的基础上适当地加入针对性强 的蒙特卡罗技巧予以解决、本文的目的是、评述MCNP程序在核测井方面的应用中，应采用哪些技巧，尚 存在哪此问题，还有哪些核测井问题直接用MCNP程序解决有困难，不适当地加人针对性强的蒙特 K罗技 巧难以很好解决，甚至根本无法解决 引 u=U7 MCNP是美国LosAlamos国家实验室众多科学家的共同成果，其中包括著名科学家 Fermi,VonNeumann,Lilam和Metropolis等人早期的关于蒙特 卜罗力法的开创性工作，也 包括许多后继者.Kahn,Cashwell和E二二tt等做出的重要贡献、 MCNP程序可以处理任意二维材料结构的粒子输运问题，几何区域的界面可以是平面、 任意二次曲面和椭圆旋转环面，几何区域中的材料可以由任意多种元素组成。 MCNP程序所采用的截面形式，按能量区分可以是连续方式的，也可以是多群方式的。 考虑了ENDF截面库给出的所有中子反应类型，对 f热中子，可以选用自由气体或S (a, 司，两种模型中的任意一种。对于光子，考虑了相千与非相干散射，并处理了光电吸收之 后可能有的荧光发射及电子对产生后的就地韧致辐射光。对于电子的考虑也是非常全面的。 MCNP程序凝聚了作者们所掌握的几乎全部重要的蒙特卡罗技巧，其中包括计算点通 量的下次事件估计 (NEE)技巧，多一次碰撞通量估计 (OMCFE)技巧，多群伴随技巧; 计算体通量的径迹长度技巧;解一般问题的指数变换技均，DXTRAN技巧，分裂与赌技巧; 解截面微扰问题的相关抽样技巧，等等 MCNP程序除r为美国政府的核武器计划做出了重要贡献外，还广泛应用于许多科学 与工业部门，在核测井模拟中的应用则是其最币要的领域之 。本文的目的是，评述MC- NP程序在核测井问题的应用中，应采用哪些技巧，尚存在哪此问题，还有哪些核测井问题 直接用MCNP程序解决有困难，不适当地加入针对性强的蒙特卡罗技巧难以很好解决;甚 至根本无法解决 核测井模拟中普遍存在的深穿透问题二’一“〕 在核测井问题中，除自然粒子测井外，大体上都属于用人工制造粒子源，粒子由源出发 13 后经过在地层中的输运到达探测器，通过探测器对到达粒子的响应了解地层情况。由于粒子 由源出发后到达探测器安在地层中经过漫长的输运方能到达，因此，属于 “探穿透问题 ”。 蒙特卡罗方法解深穿透间题所遇到的最大困难是，当到达探测器的粒子所占份预太少时计算 结果比实际结果偏低，所占份额越少情况越是如此 考虑如下简单积分的计算问题: 。一丁;;)二，)(一;J，二);(一:1;(卜工))dX (1) 其中凡与艺、为任意正实数，该积分的计算实际1一相当于这样一个简单的平板屏蔽计算问 题:屏蔽的厚度为1，粒子垂直人射，在屏蔽内碰撞前后的总截面分别为岛和马，发生碰 撞时全部散射，方向不变，第二次发生碰撞时全部被吸收，粒子穿过屏蔽的通量正好等于积 分 (飞)。 根据统计估计方法的一般原理，对于 中的蒙特卡罗估计是: 、‘一壳笼1贬Xp(一艺1“一‘。)) (2) 式中 N— 粒子的抽样数; 二。— 粒子首次碰撞位置的第，，次抽样值」 为了表明蒙特 卜罗方法解深穿透问题时存在计算结果偏低的情况，按式 :)2 进行 了具 体计算，积分的实际结果与计算结果同时被列于表 1巾。由于这里的屏蔽厚度为 1，所以若 以粒