(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202211150637.X (22)申请日 2022.09.21 (71)申请人 清华大学 地址 100084 北京市海淀区清华园1号 (72)发明人 佘顶　夏冰　石磊　 (74)专利代理 机构 北京三聚阳光知识产权代理 有限公司 1 1250 专利代理师 张建纲 (51)Int.Cl. G06F 30/23(2020.01) G06T 17/20(2006.01) (54)发明名称 球床反应堆倒料的模拟计算方法 (57)摘要 本发明涉及一种球床反应堆倒料的模拟计 算方法， 包括： 将球床反应堆的整个球床区(10) 划分为若干个球流网格(60)， 每个流道在球床反 应堆轴向上等体积地划分为若干层(50)， 每个流 道的每层中具有一个球流网格， 每个球流网格中 燃料球的批次数(M)小于等于球流网格中燃 料球 数目， 并且 各个流道的球流网格中燃料球的批次 数(M)按照各流道的球流网格的体积比例设定； 每经过一个倒料时间步长， 所有燃料球沿流道下 移一层进入下一层球流网格并进行燃耗计算； 将 从堆芯底部卸出的燃料球以随机抽样方式按照 各流道球流网格体积比例重新倒入堆芯顶部各 个流道， 在重新倒入各个流道之前一部分从堆芯 底部卸出的燃料球根据燃耗计算结果由新燃料 球替换。 权利要求书1页 说明书4页 附图4页 CN 115510703 A 2022.12.23 CN 115510703 A 1.一种球床反应堆倒料的模拟计算方法， 其特 征在于， 该模拟计算方法包括如下步骤： 将所述球床反应堆的整个球床区(10)划分为若干个球流网格(60)， 其中， 将所述球床 反应堆的整个球床区(10)在所述球床反应堆的径向上划分为若干个流道(40)， 每个所述流 道(40)在所述球床反应堆的轴向上等体积地划分为若干层(50)， 每个所述流道(40)的每一 层中具有一个所述球流网格(60)， 其中， 每个球流网格(60)中的燃料球的批次数(M)小于等 于该球流网格(60)中的燃料球的数目， 并且 各个流道的球流网格(60)中的燃料球的批次数 (M)按照各流道的球流网格的体积比例设定； 每经过一个倒料时间步长， 每个球流网格(60)中的各个燃料球分别下移到其它球流网 格(60)， 并进行燃耗计算； 将从堆芯底部卸出的燃料球以随机抽样方式按照各流道的球流网格体积比例重新倒 入堆芯顶部各个流道(40)， 其中， 在重新倒入 各个流道(40)之前， 一部 分从堆芯底部卸出的 燃料球根据燃耗计算结果由新燃料球替换。 2.根据权利要求1所述的球床反应堆倒料的模拟计算方法， 其特征在于， 每个球流网格 (60)中的所有燃料球沿所述流道(40)下移到当前所在流道(40)的位于该球流网格(60)下 方的球流网格(60)， 或者每个球流网格(60)中的一部分燃料球沿所述流道(40)下移到当前 所在流道(40)的位于该球流网格(60)下方的球流网格(60)且其余燃料球下移到相邻的其 它流道(40)的球流网格(6 0)。 3.根据权利要求1或2所述的球床反应堆倒料的模拟计算方法， 其特征在于， 所述球床 区(10)的各所述流道(40)的各球流网格(60)具有相同的体积与批次数(M)， 从堆芯底部卸 出的燃料球以随机抽样方式按照各流道的球流网格的体积比例均分地重新倒入堆芯顶部 的各个流道(40)。 4.根据权利要求3所述的球床反应堆倒料的模拟计算方法， 其特征在于， 各所述球流网 格(60)的批次数(M)等于所述球流网格(6 0)内的燃料球的实际数量。 5.根据权利要求1或2所述的球床反应堆倒料的模拟计算方法， 其特征在于， 所述球床 反应堆采用在线换 料。 6.根据权利要求5所述的球床反应堆倒料的模拟计算方法， 其特征在于， 所述球床反应 堆是球床高温气冷堆， 处于初装堆建立过程、 过渡过程或者平衡堆芯状态。 7.根据权利要求1、 2、 5 中任一项所述的球床反应堆倒料的模拟计算方法， 其特征在于， 部分所述燃料球由石墨球替代。权　利　要　求　书 1/1 页 2 CN 115510703 A 2球床反应堆倒料的模拟计算方 法 技术领域 [0001]本发明涉及一种核反应堆物理计算方法， 尤其是一种用于球床反应堆、 特别是高 温气冷堆倒料的模拟计算方法。 背景技术 [0002]在球床高温气冷堆中， 流道底部卸出的燃料球会被重新倒入到堆芯顶部进入下一 次燃料循环， 被重新倒入堆芯顶部之前， 这些燃料球中达到目标燃耗深度的乏燃料会被丢 弃并以新燃料球进行替换。 [0003]在压水堆以及其它大部分类型反应堆的燃耗计算当中， 燃料元件的燃耗历史是确 定的， 可以通过乏燃料核素密度测 量值来验证理论计算值的正确 性。 但对于球床 高温气冷 堆， 任意一个燃料球在堆芯内的运动轨迹与燃耗历史都是随机的， 从理论上无法 “准确”预 测乏燃料核素密度。 因此， 球床 高温气冷堆燃耗计算验证并不是基于特定计算值与实验值 之间的对比验证， 而是基于概 率分布或统计学意 义上的对比验证。 [0004]为模拟球床高温气冷堆的倒料过程， 现有确定论物理计算软件， 例如  PANGU与 VSOP， 采用了下文中结合图1至图3所述的倒料计算模型。 [0005]如图1所示， 首先将整个球床区在径向上划分为若干个流道(channel)， 流道边界 可以为直线或曲线； 每个流道在轴向上划分为若干层(layer)， 将整个球床区划分为若干个 体积相等的球流网格。 不同流道的层数可不同， 以此模拟各流道不同的球流速度。 每个球流 网格内的燃料球被划分为若干批次(batch)， 每个批次的燃料球具有相同物质成分， 批次数 N通常等于 燃料球通过次数。 [0006]如图3所示， 在每个倒料时间步， 所有燃料球沿流道下移一层， 进入下一层球流网 格(最下层燃料球被从流道底部卸出)， 并进行燃耗计算。 对各流道底部卸出 的燃料球的对 应批次进 行合并平均处理， 也就是将各个流道中的第n批次合并(1≤n≤N)置入卸料暂存盒 的N个独立空间中的一个， 从而保证卸出的燃料仍具有相同的批次数N。 然后在装料暂存盒 (同样有N  个独立空间)接收这些燃料盒之前， 对其中若干批次进行替换和排序等操作， 例 如丢弃燃耗最深的乏燃料并将其替换为新燃料。 最后， 将它们重新倒入到堆芯顶部各个流 道， 并认为各流道顶层的球流网格具有相同的燃料球组成， 也即每个流道都分别接收到了N 批次燃料球中的每一批次。 这里 的卸料暂存盒与装料暂存盒是抽象概念， 用于在计算过程 中模拟卸料与装料 过程。 [0007]在上述倒料模型中， 球流网格内批次N代表了不同燃耗深度的燃料球， 即认为该 网 格内包含N种不同燃耗深度的燃料球。 而网格内的实际燃料球数远大于批次数N， 仅用N个批 次无法准确反映出精细的燃耗深度分布。 在对堆芯底部卸出的燃料球进 行批次合并平均的 过程中， 以及假设进入各流道顶层的球流网格具有相同的燃料球组成的过程中， 进一步抹 平了燃耗分布的信息， 无法考虑倒料 过程的随机性。说　明　书 1/4 页 3 CN 115510703 A 3