特征统计算法CSA用于屏蔽优化研究.pdfVIP

第十一届反应堆数值计算和粒子输运学术会议暨2006年反应堆物理会议 特征统计算法CSA用于屏蔽优化 王炳衡1，刘志宏2，施工1，胡永明2 清华大学，北京，100084 摘要iCARR堆冷源孔道组合屏蔽优化设计是通过使用CSA算法f11结合—维SN程序ANISN编写的屏蔽优化程 序，得到冷源闸门上部屏蔽层内屏蔽材料优化组合，使闸门关闭时混凝土墙外侧总剂量率达到较低剂量率水平。 关键词：特征统计算法CSA：辐射屏蔽；最优化 1引言 中国先进研究碓(CARR)堆芯周围布置有中子源孔道。其中冷源孔道如图l所示，冷中 子源位于重水反射层内，外侧依次是冷源孔道、冷源闸门。并在外侧设置60cm厚的密度为 4．89／cm3的重混凝土墙。闸门如图2所示，两个可起降闸门。闸门由上下两部分组成，上部 填充屏蔽材料，下部填充重混凝土，每一闸门开有两个长方孔道。实验时，闸门处于打开状 态，粒子通过冷源孔道、闸门孔道，到达实验大厅。平时，闸门处于关闭状态。这时冷源孔 道对着闸门上部。为了在关闭状态下，使大厅的剂量率满足监督区要求(辐射剂量率限值为 0．001～O．003mSv／h)，对闸门上部屏蔽材料和厚度选择进行了优化设计。 3 3 冷中子 闸门屏蔽层 混凝土墙 图1冷源孔道屏蔽示意图 l l 2口口 2口口 图2冷源孔道闸门径向示意图 卜上部屏蔽层2一重混凝土3一闸门孔道 第十一届反应堆数值计算和粒子输运学术会议暨2006年反应堆物理会议 2冷源孑L道组合屏蔽优化 CARR堆冷源孔道组合屏蔽优化是通过使用CSA算法【l】结合一维SN程序ANISN编写的 屏蔽优化程序，得到冷源闸门上部屏蔽层内屏蔽材料优化组合，便闸门关闭时混凝土墙外侧 总剂量率达到较低剂量率水平。(以ANISN计算得到的混凝土墙外侧的剂量率作为目标函数， 优化程序就是使目标函数极小化。) ANISN程序的计算模型如图3所示，厚度为1．18米的闸门上部屏蔽层共分24层，每层 可以在4种材料(铅，铁、聚乙烯、含硼聚乙烯)中随意选择，是典型的组合优化问题。(这 里采用分层办法是为了作离散化处理，并且综合考虑计算精度与计算时间选择分24层，前23 层每层厚度为5cm，第24层厚度为3cm。优化程序通过改变ANISN输入文件中24层材料的 输入，得到不同屏蔽材料组合的剂量率值。)冷源闸门孔道、特别是闸门顶部空气部分存在泄 漏问题，ANISN程序采用一维SN计算方法【们，不能直接处理泄漏，为此可将冷源闸门材料密 度均匀化作近似处理。另外在ANISN源强的处理上，合理的选择是以MCNP计算中取的计数 面作为ANISN的源面取平均得到源强。 1 2 3 4 5 20 21 22 23 24 源 闸门(i’’24层) 混凝土墙 特征统计算法(CSA)是一种高效的全局优化算法【11，根据具体优化问题选择概率特征信 息。这里，指定的特征就是各层的材料分布。保留一组较优化的解作为样本群体，从样本群 体中提取特征信息，用以指导产生新的优化解；并将新优化的解加进样本群体，指导下一步 的搜索，直至获得最优解群体。其基本步骤如下： (a)在域内随机地给出2N个点，计算2N个点的目标函数值。 (b)将2N个点的目标函数值按大小排序，取出其中前N个较好点作为样本群体。 (c)根据指定的特征，对样本群体作统计，得到特征分布。 (d)随机给出N个方案，扰动产生新解，保留与样本群体材料分布相逼近的N个方案 (e)计算N个新点的目标函数值，将N个新点加入到样本群体中，然后返回步骤(b) 3优化算法研究 3．1搜索方法比较 683