Nuclear Science and Technology
Vol.05 No.01(2017), Article ID:19504,12
pages
10.12677/NST.2017.51001
Effects on keff and Neutron Flux of MA Transmutation in Fast Reactor
Jin Cai, Bin Liu, Haoyang Yu, Wenxin Zhang
School of Nuclear Science and Engineering, North China Electric Power University, Beijing
Received: Dec. 19th, 2016; accepted: Jan. 6th, 2017; published: Jan. 12th, 2017
Copyright © 2017 by authors and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
ABSTRACT
The radiation hazard of minor actinides (MA nuclides) in spent fuel cannot eliminate until tens of thousands of years’ decay. Partitioning and Transmutation technology (P&T) is a method to eliminate MA nuclides’ radioactive hazard. By loading into the operating nuclear reactor, MA nuclides separated from the spent fuel are transmuted to short-lived nuclides or stable nuclides. Fast reactor is efficient and mature neutron source for transmuting MA nuclides. We build a simulation model of fast reactor by MCNP code to study the effects on effective multiplication factor (keff) caused by the addition of MA nuclides and to compare different loading amounts and different nuclides include 237Np, 241Am, 243Am, 244Cm, 245Cm and the mixture of MA nuclides above (MA mixture). Results show that the addition of 237Np, 241Am, 243Am and MA mixture can cause a decrease of keff, while the addition of 244Cm and 245Cm boost the keff instead, and the effects on effective multiplication factor (keff) of the reactor core enhance with the raising of MA loading amount, which threatens the safety of the reactor, thus, we consider the minimum amount as the best loading amount. Simulation results of neutron flux show that the minimum loading amount of MA barely affects the neutron flux, which is a advantage to transmute MA nuclides in fast reactor.
Keywords:MA Nuclides, Transmutation, Fast Reactor, keff
MA核素的嬗变对快堆keff和中子通量的影响
蔡进,刘滨,于浩洋,张文鑫
华北电力大学核科学与工程学院,北京
收稿日期:2016年12月19日;录用日期:2017年1月6日;发布日期:2017年1月12日
摘 要
乏燃料中大量的次锕系元素(MA核素)要经过上万年的衰变才能消除其放射性危害。分离和嬗变技术(P&T)将乏燃料中的MA核素分离出来,再将MA核素放入反应堆中嬗变为短寿命核素或稳定核素,从而达到消除其放射性危害的目的。快堆是效率较高且技术较成熟的MA嬗变装置。为研究和比较237Np,241Am,243Am,244Cm,245Cm和混合MA核素的添加和上述核素装载量的不同对堆芯有效增殖因子(keff)的影响,采用MCNP程序搭建快堆模型进行模拟计算。计算结果显示,将237Np,241Am,243Am或混合MA核素以嬗变棒的形式替换燃料棒加入堆芯时,堆芯keff均下降,而将244Cm和245Cm嬗变棒加入堆芯时,keff反而上升,且装载量越大,keff变化越多,这将威胁到反应堆的安全性,因此我们认为方案一为最佳装载量。中子通量密度的计算结果显示,MA核素的添加几乎对堆芯的中子通量没有影响,这是使用快堆嬗变MA核素的一个优点。
关键词 :MA核素,嬗变,快堆,keff
1. 引言
能否安全合理的处置乏燃料是限制核电发展的一个重要因素,而长寿命高放射性废物的处置是乏燃料后处理所面临的主要困难 [1] 。长寿命高放射性废物中的长寿命裂变产物(Long-lived Fission Products,简称LLFP)和次锕系核素(Minor Actinides,简称MA,主要包括:237Np、241Am、243Am、244Cm和245Cm)主导着乏燃料的长期放射毒性 [2] [3] ,LLFP和MA核素要经过上万年的衰变,其放射性水平才能达到天然铀对环境的放射性水平。一般的乏燃料后处理方法如“一次性通过”直接将乏燃料包装后长期埋存或后处理核燃料循环(RFC)将长寿命高放核素从乏燃料中分离出来,经玻璃固化后在将其埋存。这两种方法都直接将长寿命高放核素埋藏在环境中,对埋藏地存在长期的潜在放射性危害。采用分离和嬗变(P&T)方法能有效处置长寿命高放核素 [4] [5] 。P&T将乏燃料中的钚、铀等可循环利用的核素以及长寿命高放核素分离出来,再将长寿命高放核素放入运行中的反应堆中嬗变,如此循环,直到他们嬗变为短寿命核素或稳定核素。热堆、快堆和超临界堆都可用于嬗变MA核素 [6] 。快堆具备高中子通量和高中子能量等特性,能够有效地嬗变MA核素 [7] [8] ,研究MA核素在快堆中的嬗变特性具有重要意义。本文采用嬗变棒的形式,将MA核素制成嬗变棒放入快堆中进行嬗变,研究不同MA核素和不同装载量对堆芯keff的影响,并以此为依据选择出最佳添加方式,并计算最佳方式下,MA核素的添加对堆芯中子通量密度的影响。
2. 快堆嬗变MA核素
如表1所示,MA核素在热堆中的裂变截面和俘获截面都比在快堆中的大,但MA核素只有发生裂变反应时才能转换成短寿命的核素,这样才能有效地消除其放射性危害。MA核素在快堆的俘获裂变比要比在热堆中小得多,即MA核素在快堆中与中子发生反应时,裂变反应所占的份额比在热堆中更高,这说明快堆能更有效地将MA核素嬗变为短寿命核素或稳定核素。已有研究表明,一座快堆可嬗变4~10
Table 1. Cross section and capture-to-fission ratio of main nuclides in thermal reactor and fast reactor [9] [13]
表1. 堆芯主要核素在热堆和快堆中的核反应截面和俘获裂变比 [9] [13]
座等功率规模的热堆产生的MA核素 [10] 。
如图1所示,主要MA核素的裂变截面都比238U大,MA核素在高能区的裂变截面随着中子能量的增大而迅速增加。所以向快堆中添加MA核素后,低能区的中子数目将会减少,高能区中子数目增加,快堆的中子能谱将会硬化 [11] 。
如图1所示,主要MA核素的裂变截面都比238U大,MA核素在高能区的裂变截面随着中子能量的增大而迅速增加。所以向快堆中添加MA核素后,低能区的中子数目将会减少,高能区中子数目增加,快堆的中子能谱将会硬化 [11] 。
3. 快堆嬗变堆芯设计
本文利用MCNP程序建立快堆嬗变模型来研究MA核素加入快堆后keff和中子通量密度分布的变化。
快堆堆芯由燃料组件、再生区组件、控制棒、屏蔽组件和作为冷却剂的液态钠组成,采用交替排列燃料组件和再生区组件的非均匀结构 [12] ,其堆芯结构如图2所示。这种结构能产生更高的增殖比并减小钠空洞系数,但要求有较高的易裂变燃料总投量,堆芯燃料只包含一种富集度。采用标准三角形栅格,将燃料棒排列成燃料组件。燃料组件和再生区组件布置如图3所示。燃料采用混合UO2-PO2燃料,MOX燃料的成分含量如表2所示。在燃料芯块上下方装有天然UO2组成的再生区芯块进行增殖,燃料棒结构如图4所示。
4. 快堆嬗变MA核素
4.1. MA核素嬗变方案
为研究MA核素的添加对快堆堆芯的影响,我们将MA核素嬗变材料制成嬗变棒,其尺寸结构与燃料棒一致,将嬗变棒插入燃料组件,在结构上替代一些燃料棒。这种结构使MA核素嬗变材料与核燃料完全隔离,有助于MA核素嬗变材料从反应堆堆芯卸载后的再加工 [6] 。
Figure 1. Fission cross section of fuel nuclides and MA nuclides
图1. 燃料和MA核素的裂变截面
Figure 2. Configuration of the core
图2. 堆芯结构
为比较不同MA核素对堆芯影响的区别,将237Np、241Am、243Am、244Cm、245Cm和上述MA核素的混合分别制成嬗变棒加入堆芯,混合MA核素嬗变棒中各核素的比例按照压水堆乏燃料中MA核素的比例混合,如表3所示。嬗变棒的结构如图5所示。
同时,为比较不同装载量对堆芯影响的不同,设计五种添加方案:方案一至五分别将每个燃料组件中的1、7、19、37、61根燃料棒换为嬗变棒,如图6所示。
4.2. MA核素装载量
堆芯中有156组燃料组件,6组交替组件,模拟计算中将交替组件置为燃料组件,各材料的密度如表4所示,由此可计算出各种方案下MOX燃料和嬗变材料的装载量以及它们的质量比,见表5。
4.3. MA核素的添加对堆芯keff的影响
使用MCNP模拟嬗变棒加入堆芯的情况,在使用重复结构排列组件时,如图6所示,在嬗变棒的位
(a) 燃料组件 (b) 再生区组件
Figure 3. Radial schematics of a fuel assembly and a blanket assembly
图3. 燃料组件和再生区组件布置
Figure 4. Axial schematics of a fuel rod
图4. 燃料棒结构
Table 2. Constituent of MOX fuel [1]
表2. MOX燃料成分 [1]
置将宇称值(即U卡的值)置为各种嬗变棒的宇称值,这样一来,便可模拟图6所示五种方案,得出不同方案下的keff值。
表6列出了以五种方案分别向临界的快堆堆芯中添加237Np、241Am、243Am、244Cm、245Cm和混合
Figure 5. Axial schematics of a transmutation rod
图5. 嬗变棒结构
Figure 6. Five different transmutation rod distribution patterns in fuel assembly
图6. 五种方案中嬗变棒的分布
Table 3. The ratio of each MA nuclide in the spent fuel of PWR
表3. 压水堆乏燃料中MA核素的比例
Table 4. The density of core material
表4. 堆芯材料密度
Table 5. The loading amount and the ratio of loading amount of MA nuclides and MOX
表5. MA核素和MOX燃料装载量(kg)和质量比(MA/MOX)
Table 6. keff after adding different MA nuclides to the fast reactor
表6. 加入不同MA核素嬗变棒后快堆的keff
MA核素嬗变棒后keff的变化,计算结果的标准偏差均保持在0.0003左右。从表中可以看出,不同的MA核素对keff的影响不同。加入237Np、241Am、243Am和混合MA核素嬗变棒将使keff有不同程度的降低,而加入244Cm和245Cm嬗变棒将使keff升高。
图7显示了不同MA核素对keff的影响程度的差异以及随着MA核素装载量的增加,keff的变化趋势。在降低keff方面,随着装载量的增加,243Am引起keff降低的能力最大,241Am次之,再其次是混合MA,237Np使keff轻微地降低。随着241Am、243Am和混合MA的装载量的增加,keff呈单调降低趋势,即241Am、243Am和混合MA的装载量越大,堆芯keff降低得越多。而随着237Np的装载量的增加,keff先降低,随后又开始上升。这说明存在一个临界值,当装载量超过该临界值时,237Np装载量的增加将不会继续降低堆芯keff。在升高keff方面,随着装载量的增加,245Cm使keff剧烈的升高,相比而言,244Cm只使keff轻微地增加。
将燃料棒替换成MA核素嬗变棒之后,MA核素以及它们的俘获产物和裂变产物对反应堆内中子平
Figure 7. The variation of keff of fast reactor along with the loading amount
图7. 快堆堆芯keff随MA核素装载量的变化
衡的影响将改变堆芯的keff。快堆没有慢化过程,堆芯内中子的平衡主要取决于裂变、快中子增殖、辐射俘获和泄漏过程,裂变反应产生中子,辐射俘获和泄漏时中子减少。由keff的定义:
在同一堆芯中,快中子增殖系数和不泄漏概率不变,故不同核素对堆芯keff的影响取决于所加核素以及它们的裂变产物和俘获产物的裂变和俘获过程,即在同一堆芯中,keff的值与俘获裂变比呈负相关。
堆芯引入MA核素嬗变材料后,部分燃料棒被嬗变棒替换,裂变材料减少。MA核素在快中子能区不仅能发生裂变反应,还能发生俘获反应。若MA核素发生俘获反应则消耗中子,使中子通量降低,反应性减小。若发生裂变反应,则放出裂变中子,补偿因裂变减少而引起的中子通量的降低和反应性的减小。表2列出了堆芯主要核素在快堆的裂变截面和俘获裂变比数值。MOX燃料中易裂变核素如235U和239Pu的裂变截面都较大,裂变性能较好。根据表3所列MOX燃料的成分含量和表2所列混合MA中核素比例可算出MOX燃料和混合MA的计算裂变截面和俘获裂变比,列在表7中。从表7可以看出237Np、241Am、243Am和混合MA的裂变截面均小于MOX燃料,俘获裂变比大于MOX燃料,故将这些MA核素嬗变棒代替部分燃料棒加入快堆堆芯后,keff下降。245Cm的裂变截面大于239Pu和235U的裂变截面,俘获裂变比也比239Pu和235U小,它在反应堆中的裂变性能比易裂变燃料还要好,故将245Cm嬗变棒加入堆芯后,keff剧烈上升。虽然244Cm的有效裂变中子数小于MOX燃料,但它裂变截面比MOX燃料要高,且其俘获产物245Cm的裂变性能很好,所以将244Cm嬗变棒加入反应堆后,keff也将上升。并且,裂变截面的大小由高到低分别为245Cm、244Cm、237Np、混合MA、241Am 、243Am,与俘获裂变比大小由低至高
的顺序一致,也这与keff变化量一致。
已有研究表明,以嬗变棒的形式向压水堆中添加1%的混合MA核素,keff降低至0.98634 [6] 。从表5和表6可以看出,同样的向快堆中添加2.39%混合MA核素嬗变棒,keff只降到0.99789,这说明MA核素的添加对快堆堆芯性能的影响比热堆小,利用快堆嬗变MA核素能取得较理想的结果。
通过修改MCNP程序材料卡中MOX燃料中各核素的比例调整MOX燃料中PuO2的含量,使反应堆重新达到临界,即在加入MA核素嬗变棒的情况下使keff重新回到1。结果如表8所示。
由于239Pu的缓发中子份额小于235U,如果燃料中239Pu的含量太高会明显降低堆芯的缓发中子份额,导致反应堆的控制困难,故MOX燃料中PuO2的含量不能超过40%。由表8可知,无论使哪种方式下的反应堆重回临界都不用使MOX燃料中PuO2的含量超过40%。但即使使239Pu和235U的含量为零,也不能使方案三、四、五的堆芯重回临界。
4.4. MA核素的添加对堆芯中子通量密度的影响
堆芯加入MA核素会影响堆芯物理性能,且装载量越大,对堆芯的扰动越大。综合队堆芯的影响和MA核素分离技术的难度,我们认为方案一下的混合MA是最佳添加方式。接下来讨论在最佳添加方式下,MA核素的添加对堆芯能谱和径向中子通量密度的影响。
MCNP程序计算中子通量时模拟的是一个中子引起的相对的情况,要得到真实的计数,必须结合堆芯功率对计算结果归一化。堆芯平均每次裂变释放的可利用能量为198 MeV,归一化因子即单位时间内系统每释放1 J所需要的裂变反应次数:
该系统热功率水平为970 MeV,每次裂变产生的平均中子数为2.5942,若keff = 1,则每提供2.5942个源中子就会发生一次裂变。则该系统的归一化因子为:
Table 7. The fission cross section and capture-to-fission ratio of main nuclides of transmutation fast reactor core
表7. 快堆嬗变堆芯主要核素的裂变截面和俘获裂变比
Table 8. Mass faction of PuO2 in MOX fuel when reactor core returning to criticality
表8. 堆芯重回临界时MOX燃料中PuO2的含量(%)
即MCNP程序对该系统进行计数是使用的归一化源强常数 [2] 。
对MCNP计算结果处理过后的堆芯能谱如图8所示。加入MA核素前后堆芯径向中子通量密度的分布情况如图9所示。蓝框部分放大如图。
由图8可知,堆芯加入MA核素后,由于MA核素在低能区俘获截面很大,所以低能区中子减少明显,而其他能区基本没有变化。由图9可看出,堆芯径向中子通量在燃料和再生区交替排列的区域分布比较平缓,加入MA核素后变化较小。说明此方案对快堆堆芯中子通量的影响较小。
综上所述,采用方案一将MA核素混合嬗变棒加入堆芯嬗变对堆芯keff、中子通量密度分布等物理性能产生的影响较小,可以取得理想的嬗变效果。
5. 结论
本文比较了五种不同MA核素以及它们的混合物以嬗变棒形式加入堆芯时,在不同装载量下对堆芯keff的影响。计算结果显示,237Np、241Am、243Am和混合MA将使keff降低。随着装载量的增加,243Am引起keff降低的幅度最大,241Am次之,再其次是混合MA,237Np使keff轻微地降低。随着241Am、243Am和混合MA的装载量的增加,keff呈单调降低趋势,而随着237Np的装载量的增加,keff先降低,随后又开始上升。通过增加MOX燃料中PuO2的含量,均可使堆芯重返临界。244Cm和245Cm的添加将使keff升高,且随着装载量的增加,245Cm使keff剧烈的升高,相比而言,244Cm只使keff轻微地增加。降低MOX燃料中PuO2的含量可使244Cm的五种方案和245Cm的方案一和方案二的堆芯重返临界,但即使使239Pu和235U
Figure 8. The neutron flux of reactor core with and without MA nuclides
图8. 堆芯加入MA核素前后的中子通量
Figure 9. The radial distribution of neutron flux of reactor core with and without MA nuclides
图9. 堆芯加入MA核素前后的径向中子通量
的含量为零,也不能使245Cm的方案三、四、五的堆芯重回临界。综合MA核素对堆芯的影响和MA核素分离技术的难度,我们认为方案一下的混合MA是最佳添加方式。在该添加方式下,MA核素的添加对堆芯能谱和径向中子通量密度的影响都很小,且该方式下MA核素的装载量可达55.644 kg,嬗变效果可观。
文章引用
蔡进,刘滨,于浩洋,张文鑫. MA核素的嬗变对快堆keff和中子通量的影响
Effects on keff and Neutron Flux of MA Transmutation in Fast Reactor[J]. 核科学与技术, 2017, 05(01): 1-12. http://dx.doi.org/10.12677/NST.2017.51001
参考文献 (References)
- 1. 周培德. MOX燃料模块快堆嬗变研究[D]: [博士学位论文]. 北京: 中国原子能科学研究院, 2000.
- 2. Artisyuk, V., Saito, M. and Stankovsky, A. (2005) Challenge of Transmutation of Long-Lived Nuclides. Progress in Nuclear Energy, 47, 327-338.
- 3. Tucek, K. (2004) Neutronic and Burnup Studies of Accelerator Driven Systems Dedicated to Nuclear Waste Trans-mutation. PhD Thesis, Department of Physics, Royal Institute of Technology, Stockholm.
- 4. Steinberg, M. (1964) Neutron Burning of Long-Lived Fission Products For Waste Disposal. BNL, New York.
- 5. 罗璋琳. 核废料核素价值研究(续前) [J]. 原子能科学技术, 2004, 38(4): 334-338.
- 6. Liu, B., Wang, K., Tu, J., Liu, F., Huang, L.M. and Hu, W.C. (2014) Transmutation of Minor Actinides in the Pressurized Water Reactors. Annals of Nuclear Energy, 64, 86-92. https://doi.org/10.1016/j.anucene.2013.09.042
- 7. Iwasaki, T. and Hirakawa, N. (1994) Neutron Economy of Transmutation of TRU in Thermal and Fast Neutron Fields. Journal of Nuclear Science and Technology, 31, 1255-1264. https://doi.org/10.1080/18811248.1994.9735288
- 8. Wakabayashi, T. (2002) Transmutation Characteristics of MA and LLFP in a Fast Reactor. Progress in Nuclear Energy, 40, 457-463. https://doi.org/10.1016/S0149-1970(02)00038-0
- 9. Liu, B., Hu, W., Wang, K., Huang, L.M., Ouyang, X., Tu, J. and Zhu, Y. (2013) Transmutation of MA in the High Flux Thermal Reactor. Journal of Nuclear Materials, 437, 95-101. https://doi.org/10.1016/j.jnucmat.2013.01.348
- 10. 徐銤. 中国实验快堆[M]. 北京: 高技术通讯, 1995.
- 11. 胡赟, 王侃, 徐銤. 钠冷氧化物燃料快堆嬗变MA研究[J]. 核动力工程, 2010, 31(1): 18-22.
- 12. 苏著亭, 叶长源, 闫凤文. 钠冷快增殖堆[M]. 北京: 原子能出版社, 1991.
- 13. Berthou, V., Degueldre, C. and Magill, J. (2003) Trans-mutation Characteristics in Thermal and Fast Neutron Spectra: Application to Americium. Journal of Nuclear Materials, 320, 156-162. https://doi.org/10.1016/S0022-3115(03)00183-1