-
熔盐堆作为六种第四代反应堆中唯一一个液态燃料反应堆,具有更好的本征安全性、防核扩散性、燃料循环方式灵活性和较高的钍资源利用率等特点[1]。2011年,中国科学院启动了“未来先进核裂变能—钍基熔盐堆核能系统”战略性先导科技专项(TMSR),旨在为我国核能大发展开辟一条符合中国国情的创新之路[2]。其中,钍铀燃料循环核数据的研究是必须首要解决的关键核心问题之一。
钍铀燃料循环核数据是指钍铀燃料循环核反应链上主要核素的核数据,如231Th、232Th、233Th、231Pa、232Pa、233Pa、234Pa、232U、233U、234U、235U等。结合钍基熔盐堆的特点,还需要考虑慢化剂、冷却剂、结构材料等材料的核数据。目前,现有核数据大多是面向铀钚燃料循环研制的,将该数据用于钍基熔盐堆,将会引入较高的核数据不确定度[3]。世界各国都相继开展了相关的研究,针对钍铀燃料循环核数据进行改进。改进的方法分为两种,一种是以经过筛选后的实验数据为基础,通过改进核数据评价方法,重新评价相关核素的核数据;另一种是结合反应堆的实际情况,重新开展关键核素的核数据测量并评价,比如,印度开展了一系列的钍铀燃料循环相关核素的截面实验测量并取得了重要的成果[4]。
为了满足TMSR对核数据的需求,中国科学院上海应用物理研究所对钍铀燃料循环核数据开展了一系列的研究与改进工作。首先,联合中国核数据中心(CNDC)研制并发布了钍铀燃料循环专用核数据库(CENDL-TMSR-V1)[5],用于钍基熔盐堆的临界计算以及屏蔽设计[6]。该数据库以钍基熔盐堆的物理设计为主要原则,结合宏观检验及微观评价的合理性,从国际主流数据库中挑选合适的核数据[7]。同时,根据钍基熔盐堆的特点,重新选择合适的实验,对232Th和233U的核数据进行了重新评价。其次,自主设计并建造了15 MeV电子直线加速器驱动的白光中子源装置[8],用于钍铀燃料循环关键核素核数据的测量,为其重新评价提供基础实验数据。
本文将介绍白光中子源及开展的一系列实验测量。首先阐述白光中子源的物理参数、主要特征以及数据获取系统。其次介绍已开展的实验测量,主要包括核级石墨、氧化钍、GH3535合金的总截面测量,以及钍铀转换率与核级石墨的硼当量实验测量。
-
基于15 MeV电子加速器驱动的白光中子源为紧凑型中子源,所有设备均集中在同一实验大厅中,其布局见图1。电子加速器可产生能量为15~20 MeV的电子,脉冲宽度为3 ns~3 μs(可调),脉冲频率为1~266 Hz(可调),平均脉冲电流为0.1 mA,设计最大功率为1.5 kW。
电子加速器产生脉冲电子,加速后达到15 MeV,通过三个挡块刮掉其余能量的电子后,使用二极磁铁对电子进行90°偏转,然后轰击钨靶。电子受到靶的作用而减速,发生轫致辐射反应,产生高能X射线。高能X射线与靶核发生(γ, n)反应,产生中子。初级中子的峰值能量约1 MeV,通过10 cm厚的聚乙烯对中子进行慢化,可获得需要的具有连续能谱的热中子。图2给出了通过MCNP5模拟计算得到的慢化前后的中子能谱。
中子从垂直电子束的90°方向出射,通过两根低真空管道进行传输。在每根管道头尾及中间的位置放置10 cm的聚乙烯(含硼5%)和5 cm厚的铅对中子进行准直。两根管道中间为样品区,用于放置测量样品。中子飞行距离为5.8 m。
-
白光中子源采用了钨靶作为中子产生靶,直径为60 mm,厚度为48 mm,纯度99.7%。钨靶通过特殊工艺直接焊接在铜底座上。在中子产生靶的不同位置安装了13只热偶温度计,用以监测中子产生靶的温度变化,达到间接监测束斑位置及大小的作用。整个中子产生靶系统采用了独立水循环进行冷却,放置在由铝、聚乙烯和铅组合而成的屏蔽壳中。该布局的主要优势在于结构简单,体积小,中子产生溢出率高[9]。
与其他中子源不同的是,受到建筑空间的限制,白光中子源的所有设备均集中在实验大厅中,虽然对中子产生靶、探测器进行了一系列屏蔽,但实验大厅中的环境中子/伽马本底仍然较高,而且实验大厅外墙壁附近的剂量也较高,严重影响了中子源的实验测量能力及日常运行。为了降低实验大厅中的环境本底,自主设计并安装了复合局部屏蔽体[10]。最终,中子能量在0.01~0.1 eV之间,有效中子与本底中子的比值达到约100:1,在实验大厅外的环境剂量降至本底水平。该项工作在2018年1月份通过专家组的技术验收。
-
白光中子源采用了对热中子较为敏感的6LiF/ZnS(Ag)探测器,直径为50 mm,厚度为0.5 mm,热中子探测效率为73.42%,不确定度为0.3%[11]。实验中采用了两个相同型号的中子探测器,一个用于测量通过实验样品后的透射中子,即TOF探测器;一个用于监测束流强度的变化,即监测探测器。
基于CAEN公司生产的数字波形采样器(Digitizer,型号为DT5720),自主开发了多功能数据获取软件系统(TMSR-Digitizer-DAQ)[12],并与实验靶控制软件进行耦合,实现了硬件控制、数据获取、波形在线显示、靶位自动定时更换、数据文件自动开机、自动报警等功能。该数据获取系统采用了自触发模式,即在加速器的一个脉冲宽度范围内,以加速器的Gamma Flash巨脉冲作为脉冲起始时间,期间的小信号(中子或者伽马)都会自触发并被记录,其时间作为事件终止时间。这种触发模式的优势在于可以记录任何触发阈值以上的信号,提高中子/伽马射线的计数率。
-
利用白光中子源,对钍基熔盐堆的关键材料开展了一系列的总截面实验测量,主要包括慢化剂、结构材料以及燃料熔盐中的一些主要核素。同时,也开展了一些应用实验的研究。
-
在钍基熔盐堆中,石墨作为慢化剂,GH3535合金作为结构材料,采用指定配比的含钍燃料熔盐(LiF-BeF2-ThF4-UF4)[13]。因此,基于白光中子源,主要测量了核级石墨、GH3535合金以及钍在低能区的总截面。
首先,对核级石墨进行了总截面的测量。图3给出了石墨总截面的测量结果,并与评价数据及已有实验数据进行了对比。在低能区,石墨的散射截面远大于吸收截面,可认为总截面近似于散射截面。由于美国的ENDF/B-VII.0评价数据库及EXFOR实验数据为常用的数据库,其可靠性得到了广大用户的一致认可[14],因此,将测量结果与ENDF/B-VII.0中石墨的热中子散射数据(ENDF/B-VII.0/TSL)以及EXFOR中的实验数据(EXFOR/J.A.H)进行了对比。图中显示,在0.0045~0.007 eV及0.015~0.03 eV之间,基于白光中子源测量得到的截面要显著低于评价数据,主要是受到石墨晶体结构的影响。在其他能区,测量值与评价数据基本一致。在J.A.H开展的石墨总截面测量实验中,当中子能量小于在0.003 eV时,截面都要高于评价数据。整体来看,基于白光中子源得到的截面要与ENDF/B-VII.0中的评价数据更为接近。
其次,测量了773 K时GH3535合金在热区的总截面,测量结果见图4。图中显示,在0.006 eV以上,测量截面与ENDF/B-VII.0库中评价数据基本符合;在0.006 eV以下,测量截面出现了明显的热中子散射效应的布拉格峰,但是评价数据库中并没有相应的热中子散射数据。测量结果表明在低能区GH3535合金也存在热中子散射效应,但是现有评价数据并未考虑热中子散射效应的影响,后续需开展相关的理论计算,对该数据进一步的验证,并分析对钍基熔盐堆的影响。
最后,钍是钍基熔盐堆中的重要燃料核素,其核数据的精确性非常重要。因此,利用白光中子源测量了钍的总截面。实验中使用的是二氧化钍,厚度为9 mm,扣除氧的总截面后可得到钍的总截面。图5给出了钍的总截面测量结果,并与评价数据(ENDF/B-VII.1)以及EXFOR数据库中已有实验数据(EXFOR/P.A.E)进行了对比。图中显示,在0.02 eV以上,测量结果、已有实验数据、评价数据基本一致。但是,在0.02 eV以下,评价数据呈光滑状态,已有实验数据有两个吸收峰,而测量结果呈现出明显的布拉格峰。氧化铀与氧化钍均为重金属燃料,在ENDF/B-VII.1热中子散射数据库中,氧化铀的热中子散射数据也存在类似结构[15],因此也需要考虑氧化钍的热中子散射效应,后续将进一步优化实验,降低钍总截面的测量不确定度,并结合相关理论模型,开展钍的热中子散射数据研究。
-
白光中子源除了可以进行总截面的测量,还可以开展其他的实验测量,如钍铀转换率、钍裂变率等中子学积分实验。目前主要完成了石墨的硼当量测量与氧化钍的辐照实验。
在反应堆设计中,通常以硼当量来表示核材料的纯度,硼当量是核材料纯度的重要指标之一[16]。当样品中含有吸收中子的杂质时,会引起宏观吸收截面的变化,宏观吸收截面的变化由穿透率的对数差表示(
$\Delta \ln TR = $ $ \ln (T{R_0}) - \ln (T{R_n})$ ,$\ln (T{R_n}) = \ln ({I_n}/I)$ ,I为入射中子束强度,${I_n}$ 为经过第n个几何形状相同、硼当量不同的靶后未经碰撞的中子束强度,n=0表示该靶的硼当量近似为零)。通过一系列已知硼当量石墨样品的测量,建立$\Delta \ln TR$ 与硼当量的关系,获得硼当量刻度曲线。根据该刻度曲线,通过白光中子源测量得到被测样品的$\Delta \ln TR$ ,则可计算得到被测石墨样品的硼当量。石墨硼当量刻度曲线见图6。其中,理论刻度曲线为通过公式直接计算得到的刻度曲线,模拟刻度曲线为通过MCNP5模拟计算得到的刻度曲线。结果显示,在考虑各种不确定度后,刻度曲线的斜率为6.35×10-5,相对不确定度为8.64%。根据该刻度曲线得到石墨测量样品的硼当量测量精度为9.69%,高于同条件下已有实验的硼当量测量精度[17]。为了检验钍铀燃料循环的转换性能,初步开展了232Th/233U转换率的中子学积分实验研究。基于白光中子源,对氧化钍进行辐照,利用高纯锗探测器测量其伽马射线谱(图7),然后计算232Th/233U的转换率,即233U的产生量占232Th初装量的份额。图中显示,在辐照前,可明显观测到228Ac在338.32和340.97 keV的特征峰。辐照后,可明显观测到233Pa在311.9以及340.48 keV的特征峰,表明二氧化钍经过辐照后,发生(
${\rm{n}},{\rm{\gamma}} $ )反应,生成233Th,233Th发生${\rm{\beta }}$ 衰变生成了233Pa。由于233Pa的半衰期为27 d,而在白光中子源上辐照时间仅为1 d,产生的233U非常少,未测量到233U。假设生成的233Pa全部转换为233U,则此次实验232Th/233U的转换率为3.2×10–12。 -
TMSR白光中子源为紧凑型中子源,主要集中于热区及共振区的核数据测量,已开展了石墨、钍、GH3535合金的总截面实验测量。结果显示,石墨的测量结果与ENDF/B-VII.0中的评价数据基本符合,GH3535合金在0.006 eV以上时测量结果与评价数据基本符合,但在0.006 eV以下则出现明显的布拉格峰,而钍在0.02 eV以下也同样出现了布拉格峰,表明在低能区GH3535合金和钍也存在热中子散射效应,后续需进一步验证。其次,还开展了一些应用实验,建立了白光中子源上的石墨硼当量测量方法,并对钍铀燃料循环转换性能进行了初步验证。测量结果表明,利用白光中子源可进行钍铀燃料循环关键核素的核数据测量,为核数据的验证与改进提供基础实验数据,以提高钍基熔盐堆物理设计的安全性和经济性。
-
摘要: 钍铀燃料循环核数据的精度和可靠性直接关系着钍基熔盐堆的安全性和经济性。目前大多数核数据都是基于铀钚燃料循环进行开发,若直接用在钍基熔盐堆上将会出现核设计不确定度较高的问题。为了提高钍基熔盐堆物理设计所需核数据的适用性,中国科学院上海应用物理研究所自行设计并建造了紧凑型的15 MeV电子加速器驱动的白光中子源(Photoneutron Source,PNS),用于开展钍铀燃料循环核数据的实验测量。该装置已通过技术验收,并进行了一系列关键核素的核数据测量,检验了现用核数据的可靠性,为相关核素的核数据评价与改进提供了基础实验数据。Abstract: The accuracy and reliability of Th-U fuel cycle nuclear data have an important impact on the safety and economy of thorium based molten salt reactor. However, the existing nuclear data, which are developed based on U-Pu fuel cycle, will arise high uncertainty in the nuclear design of thorium molten salt reactor. In order to improve the applicability of nuclear data for the physical design of thorium molten salt reactor, a compact photoneutron source (PNS) driven by a 15 MeV electron LINAC was designed and built up by Shanghai Institute of Applied Physics of CAS for the measurement of Th-U fuel cycle nuclear data. The PNS has passed the technical acceptance, and carried out a series of nuclear data measurement for key nuclides, which are used for verifying the reliability of existing nuclear data and providing the basic experimental data for the evaluation and improvement of nuclear data.
-
[1] GIF. A Technology Roadmap for Generation IV Nuclear Energy Systems[R]. Paris: OECD Nuclear Energy Agency, 2013. [2] 江绵恒, 徐洪杰, 戴志敏. 中国科学院院刊, 2012, 27(3): 366. doi: 10.3969/j.issn.1000-3045.2012.03.016 JIANG Mianheng, XU Hongjie, DAI Zhimin. Bulletin of the Chinese Academy of Sciences, 2012, 27(3): 366. (in Chinese) doi: 10.3969/j.issn.1000-3045.2012.03.016 [3] WARIS A, PRAMUDITYA S, PERKASA Y S, et al. International Journal of Nuclear Energy Science and Technology, 2013, 8(1): 1. [4] GANESAN S. Nuclear Data Development Related to Th-U Fuel Cycle in India[C]. Thorium Energy Conference, ThEC2013. 2013. doi: 10.1007/978-3-319-26542-1_31. [5] 核物理主题数据库[EB/OL]. [2020-02-12]. http://www.nuclear.csdb.cn. [6] 王小鹤, 胡继峰, 陈金根, 等. 原子能科学技术, 2019, 53(8): 1. doi: 10.7538/yzk.2018.youxian.0749 WANG Xiaohe, HU Jifeng, CHEN Jingen, et al. Atomic Energy Science and Technology, 2019, 53(8): 1. (in Chinese) doi: 10.7538/yzk.2018.youxian.0749 [7] 赵秋娟, 吴海成, 吴小飞, 等. 原子核物理评论, 2017, 34(3): 682. doi: 10.11804/NuclPhysRev.34.03.682 ZHAO Qiujuan, WU Haicheng, WU Xiaofei, et al. Nuclear Physics Review, 2017, 34(3): 682. (in Chinese) doi: 10.11804/NuclPhysRev.34.03.682 [8] 王宏伟, 陈金根, 蔡翔舟, 等. 核技术, 2014, 37(10): 100522. doi: 10.11889/j.0253-3219.2014.hjs.37.100522 WANG Hongwei, CHEN Jingen, CAI Xiangzhou, et al. Nuclear Techniques, 2014, 37(10): 100522. (in Chinese) doi: 10.11889/j.0253-3219.2014.hjs.37.100522 [9] LIN Zuokang, SUN Guomin, CHEN jingen, et al. Nuclear Science and Techniques, 2012, 23: 272. doi: 10.1016/j.nima.2012.07.050 [10] WANG X H, LIU L X, HU J F, et al. Progress in Nuclear Energy, 2019, 115: 74. doi: 10.1016/j.pnucene.2019.03.026 [11] 朱亮, 刘龙祥, 王宏伟, 等. 原子核物理评论, 2016, 33(3): 308. doi: 10.11804/NuclPhysRev.33.03.308 ZHU Liang, LIU Longxiang, WANG Hongwei, et al. Nuclear Physics Review, 2016, 33(3): 308. (in Chinese) doi: 10.11804/NuclPhysRev.33.03.308 [12] LIU L X, WANG H W, MA Y G, et al. Chinese Physics C, 2017, 40(5): 70. [13] JI Ruimin, LI Minghai, ZOU Yang, et al. Nuclear Science and Techniques, 2017, 28(6): 76. doi: 10.1007/s41365-017-0234-7 [14] MOSTELLER R. Nuclear Data Sheets, 2014, 118: 442. doi: 10.1016/j.nds.2014.04.101 [15] WORMALD J L, HAWARI A I. Transactions of the American Nuclear Society, 2016, 115: 1156. [16] IAEA. Guidelines for Nuclear Transfers[S]. INFCIRC/254/Rev.5/Part 1, 2002. [17] BOLEWSKIJR A, CIECHANOWSKI M, DYDEJCZYK A, et al. Nucl Instr and Meth B, 2005, 237(3-4): 602. doi: 10.1016/j.nimb.2005.03.002 -