相关抽样方法及其在煤浆系统的研究中的应用.pdfVIP

’98垒田棱仪器应用学术会议论文集 of InshMrmnt ProooeOngstIe’98Nat]coalC∞暗硼∞l*acloar 相关抽样方法及在煤浆系统研究中的应用 王汝赡姜宏宇 (清华大学工程物理系，北京，100084) 本文介绍了一种有效的方差减小技巧——相关抽样．通过网格权重窗方法和相关抽样方法的结合使用 大大提高了Monte-Carlo方法的计算效率．我们将此方法用于实际计算中。取得了很好的效果． 关键词；Monte-Carlo方法相关抽样煤浆系统 1引言 煤浆浮选对于合理利用煤能源及减少环境污染有重要的意义，煤浆在线检测系统是煤浆 浮选过程中不可缺的重要部分．由于煤浆中固、渡、气三相的存在，使得通常采用的双能(射线 法将产生较大误差．一种比较好的解决方法是在双能7射线法的基础上引入热中予探头来测 量煤浆单位体积中的含氢量．在此基础上结合密度测量则可消除气泡的影响． 在煤浆测量系统的研究中，确定在各种几何条件下固体物含量，充气量，以及温度等可变 化因素对中子计数的影响是探测部分的重点．为解决这一问题，实验方法费用太高。而且无法 完全分离各种影响因素的分别贡献。采用求解输运方程的方法，因系统的几何结构复杂结果不 够准确。Monte—Carlo方法对于复杂几何结构的处理能力很强且计算结果较为准确，缺点是收 敛速度慢。实际应用中，能否采用适合于问题特点的有效的减方差技巧往往成为能否成功进行 Monte-Carlo模拟的关键。 在我们为研究煤浆系统中子输运而开发的Monte—Carlo程序中，网格权重窗方法被作为 主要的减方差技巧，在实际应用中取得了很好的效果m”。在此基础上，我们又引入了相关抽样 方法用来计算系统中某个因素变化时对中子探测器计数的影响。该方法与权重窗方法的结合 使用，取得了比单独使用权重窗方法更好的效果． 2基本原理 相关抽样方法可以估计系统中某个因素发生变化时对探测器响应的影响．通常情况下，当 这种变化比较小时，引起的探测器响应的变化相应也比较小，这时若采用通常的计算过程，由 于统计误差的影响，往往需要达到很高的计算精度时。这种因素对计算结果的影响才能比较明 显表现出来．而采用相关抽样方法，则可在计算结果的精度不很高的情况下，充分反映出计算 结果的变化．在我们的Monte—Carlo程序中，引入了成分和密度的相关抽样．我们使用相关抽 样研究了以下两种情况： a)改变煤浆中水和煤的重量比，相应煤浆中各种元素的比例发生变化，同时煤浆的密度 也发生变化。 84 b)在煤浆中充入气泡，气泡的引入不改变煤浆中各种元素的比例，而只改变煤浆的密度． 在计算过程中，按照正常的方式计算某一个样品，这个样品称为参考样品(reference sample)，除此之外，同时计算另外一组样品，这些样品在元素比侧，密度上和参考样品有所不 sample)．在计算过程中，强迫各个比较样品中的中子沿 同，这些样品称为相关样品(correlation 与参考样品中的中子同样的径迹运动，而对比较样品的偏置影响通过对中子权重(weight)的 修正进行补偿．由于对中子跟踪过程中比较费时间的部分如核素微观截面的获得和中子几何 位置的处理对所有样品一共只需处理一次，因此和所有样品分别单独计算相比，计算效率大大 提高。 记r』为对第j个比较样品的权重修正因子，在各种情况下的处理如下： 1)发射源中子时。源中子的能量和方向的确定显然和煤浆的成分和充气量无关，所以～ 一1． 2)抽样从第扛1次碰撞点到到第i次碰撞点的距离．这个抽样过程主要依赖中子与煤浆作 用的总截面，显然总截面与煤浆组成关系密切，应进行修正．对参考样品，按正常方式，由公式 1 d=一告ln(t)}为0～1间的随机数 -R Pc』，P且分别为第，个比较样品和参考样品在d～d+△d内发生碰撞的概率 3)抽样第i次碰撞与中子发生作用的靶核以及发生的反应类型 若对参考样品用正常方法抽样出与D核发生碰撞，相应的与D核发