一种基于中子平均能量的次临界系统次临界度测量方法
【技术领域】
[0001] 本发明涉及测量次临界系统次临界度测量的技术领域,尤其涉及一种基于中子平 均能量的次临界系统次临界度测量方法。
【背景技术】
[0002] 加速器驱动次临界系统通过外源中子驱动裂变堆芯裂变产生裂变中子,同时释放 能量。由于外源中子和裂变中子能谱和分布不同,其在堆芯不同位置处的平均能量也不相 同;次临界系统中,与次临界度对应的系统有效增殖因子k rff可以表示为裂变中子在所有 中子的份额。
[0003] 由于堆芯所有中子的平均能量除与外源中子平均能量Es和裂变中子平均能量E f 相关外,还与两种中子的比例相关,因此可能通过所有中子的平均能量求得次临界系统的 次临界度。
[0004] 在已有的手段中,可以通过SuperMC等程序计算给出外源中子和裂变中子单独存 在时各自的平均能量,也可以通过替代燃料的方法测量出各自的平均能量。
【发明内容】
[0005] 本发明的目的在于:提供一种基于中子平均能量的次临界系统次临界度测量方 法,该方法避免传统测量方法由于其基于点堆假设而在次临界系统中不适用的弊端。本发 明涉及的方法直接基于次临界度的含义,通过测量堆芯中子的平均能量进行测量。
[0006] 本发明采用的技术方案为:一种基于中子平均能量的次临界系统次临界度测量方 法,所述该次临界系统为带外源中子的加速器驱动次临界系统,所述方法包括以下步骤:
[0007]步骤(101)、通过程序计算或替代燃料的实验,得到堆芯中单个或多个位置P处的 外源中子平均能量Es,_;
[0008] 步骤(102)、通过程序计算或临界实验,得到上述步骤中位置P处的裂变中子平均 能量 Ef,ave;
[0009] 步骤(103)、通过程序计算,得到外源中子相对与裂变中子的中子价值
[0010]步骤(104)、在次临界系统运行过程中,测量得到上述步骤中P处位置所有中子的 平均能量,通过以下关系式得到系统源有效增殖因子ks:
[0012] 其中,Eave为整个堆芯内所有中子的平均能量;
[0013] 步骤(105)、在得到源有效增殖因子匕之后,由以下关系式得到待测次临界度 keff:
[0014]
[0015] 其中,其所需参数Es,_、Ef,_、^未知情况下,通过如下步骤完成次临界度测量:
[0016] 步骤一、多次改变次临界系统的状态,测量得到这些状态下探测器位置处的中子 平均能量;
[0017] 步骤二、对于上述步骤一中的次临界状态,通过其它测量手段得到对应的次临界 度;
[0018] 步骤三、将上述步骤一、二中一系列次临界状态的中子平均能量与次临界度建立 刻度曲线;
[0019] 步骤四、对任一次待测临界状态,测量得到探测器位置处的中子平均能量,之后利 用上述步骤三中的刻度曲线得到系统的次临界度值。
[0020] 其中,在计算或测量中子平均能量时,统计的中子可选择所有能量的中子或是仅 选择高能量段例如大于等于7MeV中子并获得不同测量精度。
[0021] 本发明与现有技术相比的优点在于:
[0022] (1)、本发明方法利用外源中子与裂变中子平均能量不同进行测量,对于加速器驱 动次临界系统(ADS)具有天然适用性;
[0023] (2)、本发明方法原理上未做点堆假设,具有更高的测量精度,可以测量较深的次 临界度。
【附图说明】
[0024] 为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现 有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发 明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以 根据这些附图获得其他的附图。
[0025] 图1表示的是本发明一实例的测量方法流程图;
[0026] 图2表示的是本发明实例提供的铅冷快堆的堆芯布置示意图;
[0027] 图3利用所有中子的平均能量建立刻度曲线并进行测量得到的测量值与真实值 的比较;
[0028] 图4表示的是利用高能中子的平均能量建立刻度曲线并进行测量得到的测量值 与真实值的比较。
【具体实施方式】
[0029] 为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例 中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是 本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员 在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
[0030] 本发明提供的一种基于中子平均能量的次临界系统次临界度测量方法,其流程图 如图1所示。包括如下步骤:
[0031]步骤(101)、通过SuperMC等程序关闭燃料裂变截面进行计算或通过替代燃料的 实验,得到堆芯中单个或多个位置的外源中子平均能量;
[0032]步骤(102)、通过SuperMC等程序进行计算或临界实验,得到步骤(101)中同样位 置的裂变中子平均能量;
[0033]步骤(103)、通过SuperMC等程序进行计算,得到外源中子相对于裂变中子的中子 价值;
[0034]步骤(104)、在堆芯中与步骤(101)相同位置处,布置单个或多个探测器测量中子 的平均能量;
[0035]步骤(105)、利用上述步骤中得到的参数,求得系统的次临界度。
[0036] 可选的,上述步骤(101)至(103)的参数在无法得到的情况下,可以通过其它测量 方法(如脉冲源法)测量多个状态(通过调节控制棒等)的次临界度,再由步骤(104)中 测得的平均能量,反推得到上述步骤(101)至(103)的参数。
[0037] 可选的,在不实现上述步骤(101)至(103)的前提下,可以通过其它测量方法(如 脉冲源法)测量多个状态(通过调节控制棒等)的次临界度,再与步骤(104)中测得的平 均能量建立刻度曲线。之后对于待测的状态,在通过步骤(104)测得中子的平均能量后,可 以根据该刻度曲线进行测量。
[0038] 可选的,在上述步骤(101)、(102)和(104)中,可以只得到高能中子的平均能量, 并利用高能中子的平均能量进行测量,可以显著提高测量精度。
[0039] 为解决传统测量方法在次临界系统中不适用的问题,本申请提供了一种不基于点 堆方程而是基于平均能量的次临界度测量方法。该方法原理简单,易于实现,其原理为:
[0040] 从已有的资料中,可以知道系统有效增殖因子keff与源有效增殖因子k 3
具有如下 的关系:
[0042] 上式中的f表示的是外源中子相对于裂变中子的价值;
[0043] 假设某一稳态时刻,堆芯内中子中有ns个由外源直接释放的中子,有n/h裂变产 生的中子,其平均能量分别为Es,_和Ef,_。则整个堆芯内所有中子的平均能量£_为:
[0045] 由源有效增殖因子匕的定义可以知道:
[0047] 由⑵式和⑶式可以得到:
[0048] Eave= (1 -ks)Es; ave+ksEf; ave=E s; ave-ks (Es; ave-Ef;ave) (4)
[0050] 在测得中子平均能量后,可由(5)式得到源有效增殖因子匕,再由(1)式可得到系 统有效增殖因子k rff,进而可以得到系统的次临界度。
[0051] 可选的,在上述原理介绍中,可以通过只对高能中子进行分析以提高测量精度。
[0052] 可选的,如果外源中子平均能量、裂变中子平均能量Ef,_以及外源中子价值 不可知的情况下,假设:
[0053] (1)外源中子和裂变中子在次临界度改变时能谱不发生或发生很小变化,即可以 认为Es,ave和Ef,J呆持不变;
[0054] (2)外源中子价值(效率)^在各种工况下不发生变化,或退而求其次,只随次临 界度改变而改变,这实际上是忽略了外源价值在不同工况时的改变;
[0055] 在如上假设前提下,则可以通过建立中子平均能量与次临界度的刻度曲线进行测 量。
[0056] 为详细说明本发明申请的实施方式,下面在一个参考堆芯上进行了通过刻度曲线 进行测量的模拟实验:
[0057] 图2是一个铅冷次临界堆芯示意图,其中外源中子在中间组件产生,裂变中子在 燃料区域产生。
[0058] (1)首先通过同步调节控制棒得到不同的状态,得到不同状态的理论次临界度: 这里使用SuperMC计算得到的理论k eff代替其它实验手段(如脉冲源法)的测量结果;
[0059] (2)对于步骤(1)中得到的不同状态,通过"模拟"探测器一1-6#探测器进行模拟 测量,可以得到堆芯中子平均能量;
[0060] (3)使用步骤⑴和⑵中得到的理论次临界度与中子平均能量建立刻度曲线;
[0061] (4)对于次临界系统的任意不同状态(如异步调节各个控制棒的棒位),利用1-6# 探测器模拟测量中子平均能量,并利用步骤(3)中得到的刻度曲线进行测量,得到待测状 态的"测量"次临界度。
[0062] 对图2所示的模拟堆芯进行了模拟测量,期间利用程序直接计算模拟得到堆芯中 子平均能量,结果如图3所示:其中横坐标为堆芯所有中子的平均能量,纵坐标为各个不同 状态的理论次临界度。在图3中还标识出了五个待测状态相对于刻度曲线的位置。
[0063] 本申请发明还提供了进一步提高测量精度的方法:
[0064] (1)首先通过同步调节控制棒得到不同的状态,得到不同状态的理论次临界度: 这里使用SuperMC计算得到的理论k eff代替其它实验手段(如脉冲源法)的测量结果;
[0065] (2)对于步骤(1)中得到的不同状态,通过"模拟"探测器一1-6#探测器进行模拟 测量,可以得到堆芯高能中子的平均能量;
[0066] (3)使用步骤(1)和(2)中得到的理论次临界度与高能中子平均能量建立刻度曲 线;
[0067] (4)对于次临界系统的任意不同状态(如异步调节各个控制棒的棒位),利用1-6# 探测器模拟测量高能中子平均能量,并利用步骤(3)中得到的刻度曲线进行测量,得到待 测状态的"测量"次临界度。
[0068] 对图2所示的模拟堆芯进行了模拟测量,结果如图4所示:其中横坐标为高能中子 (7MeV以上)的平均能量,纵坐标为各个不同状态的理论次临界度。在图4中还标识出了五 个待测状态相对于刻度曲线的位置。
【主权项】
1. 一种基于中子平均能量的次临界系统次临界度测量方法,所述该次临界系统为带外 源中子的加速器驱动次临界系统,其特征在于:所述方法包括以下步骤: 步骤(101)、通过程序计算或替代燃料的实验,得到堆芯中单个或多个位置P处的外源 中子平均能量Es,_; 步骤(102)、通过程序计算或临界实验,得到上述步骤中位置P处的裂变中子平均能量 Ef,ave; 步骤(103)、通过程序计算,得到外源中子相对与裂变中子的中子价值^; 步骤(104)、在次临界系统运行过程中,测量得到上述步骤中P处位置所有中子的平均 能量,通过以下关系式得到系统源有效增殖因子ks:其中,Eav^整个堆芯内所有中子的平均能量; 步骤(105)、在得到源有效增殖因子ks之后,由以下关系式得到待测次临界度k rff:2. 根据权利要求1所述的基于中子平均能量的次临界系统次临界度测量方法,其特征 在于:其所需参数Es,_、EtaTC、^未知情况下,通过如下步骤完成次临界度测量: 步骤一、多次改变次临界系统的状态,测量得到这些状态下探测器位置处的中子平均 能量; 步骤二、对于上述步骤一中的次临界状态,通过其它测量手段得到对应的次临界度; 步骤三、将上述步骤一、二中一系列次临界状态的中子平均能量与次临界度建立刻度 曲线; 步骤四、对任一次待测临界状态,测量得到探测器位置处的中子平均能量,之后利用上 述步骤三中的刻度曲线得到系统的次临界度值。3. 根据权利要求1或2所述的基于中子平均能量的次临界系统次临界度测量方法,其 特征在于:在计算或测量中子平均能量时,统计的中子可选择所有能量的中子或是仅选择 高能量段即大于等于7MeV中子并获得不同测量精度。
【专利摘要】本发明提供一种基于中子平均能量的次临界系统次临界度测量方法,其包括:在堆芯布置单个或多个探测器,测量探测器位置所有中子的平均能量;通过SuperMC等程序或替代燃料等实验方法,得到探测器位置处外源中子和裂变中子各自的平均能量;通过SuperMC等程序计算外源中子相对裂变中子的中子价值;利用得到的外源中子平均能量、裂变中子平均能量和外源中子价值,从探测得到的所有中子平均能量计算得到待测的次临界度;也可以在测量之前得到中子平均能量与次临界度的刻度曲线,对于待测状态,在测得所有中子平均能量之后,利用刻度曲线得到待测状态的次临界度;在上述过程中的测量或计算中子平均能量过程,只测量或计算高能部分中子,可以显著提高测量精度。
【IPC分类】G01T3/00
【公开号】CN104898155
【申请号】CN201510274557
【发明人】杨英坤, 常博, 刘超, 孙艳婷
【申请人】中国科学院合肥物质科学研究院
【公开日】2015年9月9日
【申请日】2015年5月26日