Analysis of different burnup calculation models on nuclide components of spent fuel assembly in commercial pressurized water reactor
-
摘要: 燃耗计算精度对提高乏燃料贮存效率有着重要影响,在应用燃耗信用制时,燃耗计算得到的核素成分偏差决定了乏燃料贮存的临界安全裕量。不同燃耗计算模型所得到的核素成分偏差各不相同,为提高燃耗计算精度,提出了一种装载不同燃料富集度的多组件燃耗计算模型,并使用不同燃耗计算模型分别对TMI-1反应堆NJ07OG组件中的6个样本进行了计算、对比和分析。结果表明,相比其他模型,考虑不同燃料富集度的多组件模型得到的235U、238U和239Pu等核素平均相对偏差更接近于零且6个样本的相对偏差分布更为平均。Abstract: Burnup credit has an important impact on improving the efficiency of spent fuel storage. In the burnup credit, the burnup calculation model can affect the nuclide composition deviation, and the more accurate the nuclide composition, the lower the critical safety margin for spent fuel storage. To improve the accuracy of the burnup calculation, a multi-assembly burnup calculation model loaded with different fuel enrichment is proposed in this paper. Six samples of TMI-1 reactor NJ07OG assemblies were calculated, compared and analyzed by using different burnup calculation models. The results show that the average relative deviations of 235U, 238U and 239Pu obtained from the multi-assembly burnup model with different fuel enrichment are closer to zero and the relative deviations are more evenly distributed among the six samples than that of other models.
-
Key words:
- burnup calculation model /
- SFCOMPO-2.0 /
- nuclide deviation /
- burnup credit
-
脉冲功率技术广泛应用于激光聚变、离子束、微波功率源、雷达发射机等领域,脉冲功率是通过高压脉冲电容器储存能量,闸流管或固态开关的快速导通释放能量来获取所需的高脉冲功率,而脉冲电容的高压充电电源是脉冲功率设备充电的关键[1]。高压充电电源拓扑包括电阻限流型恒压充电、LC工频谐振型充电、桥式谐振变换器等,其中,桥式谐振型中的全桥串联谐振变换器以具备恒流充电特性、负载适应性强、充电效率高、软开关、功率密度大等一系列优点得到了广泛的讨论和应用[1-8]。文献[1-3]对LC串联谐振变换器的基本原理和应用有详细介绍;文献[4]分析了采用固定占空比调制的串联谐振变换器,因工程中变压器分布电容、二极管级间电容等寄生电容的存在,实际恒流充电特性会受到影响。文献[5]基于LCC串并联谐振变换器,LCC兼有串联谐振抗负载短路能力和并联谐振抗负载开路能力的优点,更适用于脉冲电容负载;但峰值电流的减小会影响充电能力,在高重复频率应用时会受到限制;文献[9]针对双谐振拓扑,搭建样机,测试了双谐振变换器在输入电压变化时的恒流充电能力。
本文在分析串联谐振变换器的基础上,建立双谐振变换器的数学模型,根据电压电流传输特性曲线中双谐振变换器存在第二谐振点这一特征,提出一种脉冲电容的充电控制方式。该控制方式将充电过程分为两个阶段,充电阶段提高开关频率,等效为LC串联谐振,具备恒流充电、软开关等优点;充电保持阶段,降低开关频率至第二谐振点附近,通过较小的充电电流补偿脉冲电容的自放电压降,这将显著提高脉冲电容的高压稳定度。基于Matlab/simulink环境搭建仿真模型,相同参数下,通过串联谐振和双谐振变换器的仿真对比,验证了所提出的双谐振变换器的变频率调制的可行性,这对于进一步提高脉冲功率设备的脉冲高压重复稳定度提供一个可选的参考,对后期样机的研制提供理论基础。
1. LC串联谐振与双谐振电路分析
双谐振变换器的电路拓扑结构如图 1所示,其由直流侧母线电容、全桥开关管Q1~Q4、谐振电感Lr、谐振电容Cr、谐振电感Lp、变比为n的高频变压器T、续流二极管D1~D4、高频整流桥D5~D8组成,负载储能电容Cload。若不考虑谐振电感Lp,电路即为LC串联谐振变换器拓扑。
对于LC串联谐振变换器,假定负载电容较大,等效到变压器原边为n2C(若n2C≫Cr,则对谐振回路的影响可忽略),电路阻抗Z=√Lr/Cr,谐振频率为fr=1/(2π√LrCr),串联谐振充电电流平均值为I0=8UdcfsCr/n。其中,Udc为直流侧母线电压,fs为开关频率。
根据谐振频率fr与开关频率fs的关系,串联谐振有3种工作模式:(a) 当fs < fr/2,电流断续工作,软开关,实现零电流开通和关断;(b) fr/2 < fs < fr,连续工作模式,谐振回路呈容性;(c) fs>fr,连续工作模式,谐振回路呈感性。为保持高效率,一般工作在断续工作模式较多。图 2为串联谐振断续工作模式下的谐振电流波形,一个开关周期有2个谐振周期[10-12]。
相比于串联谐振,双谐振变换器在谐振电容Cr上并联谐振电感Lp。因此工作中将存在Lr、Cr、Lp的串联谐振回路和Cr、Lp组成的并联谐振回路,两个回路的谐振频率分别用fr1和fr2表示。
fr1=12π√Lp+LLpLCr (1) fr2=12π√1LpCr (2) 2. 双谐振变换器的传输特性及比较
根据等效电路,构建串联谐振电压与双谐振电路电压的传递函数,其输入输出电压传输特性分别为
Ggainl =VoVi=11+jQ(fn−1/fn) (3) Ggain2=VoVi=11+jfnQ[1+1/(k−f2n)] (4) 式中:Q为品质因数;fn=fs/fr为标幺化频率;k=Lr/Lp。串联谐振与双谐振对应的传输特性曲线如图 3所示。
由图 3可知,对于双谐振电路,当fr2 < f < fr1之间,其电压传输特性与串联谐振基本一致;因此开关频率工作在此区间时充电特性近似保持一致。当开关频率工作在谐振点fr2时,双谐振电路的电压增益为0,此时Lp、Cr并联回路的等效阻抗无穷大,回路等效为开路。
构建串联谐振和双谐振的电流传递函数,其传输特性如下
io1=ViZ=Vi[1+jQ(fn−1/fn)]R=ViZr1/Q+jZr1(fn−1/fn) (5) io2=ViZ=Vi[1+jfnQ(1+1/(k−f2n))]R=ViZr2/Q+jfnZr2[1+1/(k−f2n)] (6) 对于串联谐振和双谐振,其特征阻抗Zr1和Zr2可表示为Zr1=√Lr/Cr,Zr2=√(Lp+Lr)Lr/LpCr;其电流传输特性曲线如图 4所示。由图可知,在谐振点附近回路电流最大。当开关频率偏离谐振点时,偏离的越远,当负载变化时,电流的增益曲线基本不受影响,此时电路变现出良好的恒流特性,以及抗短路特性。而对于双谐振变换器,与电压传输特性一样,因第二个谐振点的存在,当fr2 < f < fr1之间,其电流传输特性与串联谐振基本一致;而当无线接近fr2时,电流将趋近于0,而此时变换器处于工作状态。利用这一特性,提出基于双谐振变换器的变频控制方式,即充电阶段,双谐振变换器以一定开关频率(fs < fr1/2)工作在断续模式给储能电容充电;当电压处于保持阶段时,通过频率调制,将开关频率调制到逼近fr2,用极小的电流来补偿储能电容的电压降,从而提高充电电源的精度。这对于脉冲功率设备尤其是存在大电容的固态调制器而言,对提高脉变次级的脉冲高压稳定度具有相当的优势。
3. 仿真实验验证
基于Matlab/simulink搭建仿真模型,如图 5所示,设定Udc=540 V,Lr=3.6 mH,Cr=0.6 mH,Lp=220 mH。对于串联谐振变换器,其谐振频率fr=108.346 kHz,为了实现软开关状态,实现零电流开通,电流断续工作模式,即开关频率fs < fr/2。取理想化的条件,开关频率为33 kHz,则计算得到串联谐振变换器的充电电流平均值为0.8 A;图 6为串联谐振对应的驱动、谐振电流、谐振电压的仿真结果。串联谐振采用的控制方式为充电阶段恒流给储能电容充电,因此电压曲线呈线性变化;当恒流充电到所需的高压以后,停止充电;此时,处于高压保持阶段,因电容漏电流的存在电容自放电,电容电压将以一定时间常数下降。对于双谐振,第一个谐振点fr1 =109.229 kHz,第二个谐振点fr2 =13.86 kHz,因谐振电感Lp的存在,在相同开关频率下,充电电流要略小于LC串流谐振变换器,所以相同的高压下,充电时间要略长,这与图 4所示的电流传输特性曲线相吻合。当处于高压保持阶段,因第二个谐振点的存在,则可以通过频率调制将开关频率接近该谐振点,以保持微弱的恒流充电状态从而补偿储能电容的电压降,从而保持储能电容电压的恒定,这将大幅提高后级脉冲高压的重复稳定度。图 7(b)图为该状态的放大波形,可观测到串联谐振充电的电压的阶梯增长波形,以及电压保持阶段的电压补偿效果。
图 8为双谐振变换器充电及电压保持阶段,驱动、谐振电流、谐振电容电压以及高压的波形变化曲线,可监测频率调制过程,因仿真参数较为理想,频率调制过程变化明显,变频率调制过程中储能电容保持了较好的控制精度。图 9(a)图为充电阶段,驱动、谐振电流、谐振电容电压的放大图,与图 6基本保持一致,图 9(b)为充电保持阶段,驱动、谐振电流、谐振电容电压的放大图。
4. 结论
本文针对双谐振拓扑结构变换器,根据双谐振电路存在两个谐振点的特性,提出采用频率调制的脉冲电容器充电电源充电控制方式,即充电分为两个阶段,充电阶段采用串联谐振工作模式实现快速的恒流充电,高压保持阶段降低开关频率至接近第二谐振点,实现脉冲电容自放电压降的动态补偿,从而保证充电电源的充电精度的同时,提高脉冲电容的高压稳定度,本文详细给出了双谐振充电电源的原理及工作模式。仿真结果验证了在双谐振充电电源对脉冲电容自放电压降的补偿效果,这对于更高精度脉冲电源的研制具有参考意义。
-
表 1 燃料组件几何参数表
Table 1. Fuel assembly geometric parameters
(mm) fuel pellet inner
diameterclad inner
diameterclad outer
diametercell pitch absorber rod
pellet diameterabsorber rod cladding
inner diameter9.40 9.58 10.92 14.43 8.64 9.14 absorber rod cladding
outer diameterguide tube
inner diameterguide tube
outer diameterinstrument tube
inner diameterinstrument tube
outer diameterassembly
pitch10.92 12.65 13.46 11.2 12.52 218.11 表 2 选取的核素列表
Table 2. Nuclides chosen
actinide nuclides fission products 234U, 235U, 236U, 238U 151Eu, 153Eu, 143Nd, 145Nd, 148Nd 238Pu, 239Pu, 240Pu, 241Pu, 242Pu 147Sm, 149Sm, 150Sm, 151Sm, 152Sm 237Np, 241Am, 243Am, 244Cm 155Gd -
[1] 肖雨生. 中国核电发展与乏燃料贮存及后处理的关系[J]. 电工技术, 2020(18):24-25,57 doi: 10.19768/j.cnki.dgjs.2020.18.010Xiao Yusheng. Relationship between China’s nuclear power development and spent fuel storage and reprocessing[J]. Electric Engineering, 2020(18): 24-25,57 doi: 10.19768/j.cnki.dgjs.2020.18.010 [2] 徐健, 王伟, 黄庆勇, 等. 国外核电厂乏燃料贮存方式对比研究[J]. 中国核电, 2021, 14(6):901-909Xu Jian, Wang Wei, Huang Qingyong, et al. Comparison study on spent fuel storage of nuclear power plant[J]. China Nuclear Power, 2021, 14(6): 901-909 [3] 杨海峰, 霍小东, 易璇, 等. 高密度乏燃料贮存格架临界安全设计[J]. 核动力工程, 2014, 35(s2):167-169 doi: 10.13832/j.jnpe.2014.S2.0167Yang Haifeng, Huo Xiaodong, Yi Xuan, et al. Criticality safety design for high density spent fuel storage rack[J]. Nuclear Power Engineering, 2014, 35(s2): 167-169 doi: 10.13832/j.jnpe.2014.S2.0167 [4] 倪梓宁, 于涛, 谢金森, 等. 燃耗信任制下燃耗计算对临界计算的偏差及不确定度的研究[J]. 原子能科学技术, 2022, 56(9):1915-1923Ni Zining, Yu Tao, Xie Jinsen, et al. Research on critical calculation bias and uncertainty from burnup calculation based on burnup credit[J]. Atomic Energy Science and Technology, 2022, 56(9): 1915-1923 [5] Radulescu G, Gauld I C, Ilas G, et al. Approach for validating actinide and fission product compositions for burnup credit criticality safety analyses[J]. Nuclear Technology, 2014, 188(2): 154-171. doi: 10.13182/NT13-154 [6] Michel-Sendis F, Gauld I C, Martinez J S, et al. SFCOMPO-2.0: an OECD NEA database of spent nuclear fuel isotopic assays, reactor design specifications, and operating data[J]. Annals of Nuclear Energy, 2017, 110: 779-788. doi: 10.1016/j.anucene.2017.07.022 [7] Gauld I C, Giaquinto J M, Delashmitt J S, et al. Re-evaluation of spent nuclear fuel assay data for the Three Mile Island unit 1 reactor and application to code validation[J]. Annals of Nuclear Energy, 2016, 87: 267-281. doi: 10.1016/j.anucene.2015.08.026 [8] Hermann O W, Bowman S M, Parks C V, et al. Validation of the SCALE system for PWR spent fuel isotopic composition analyses[R]. ORNL/TM-12667, 1995. [9] 汪天雄, 张滕飞, 吴海成, 等. 基于组件计算的燃耗实验基准题建模分析[J]. 核技术, 2020, 43:060003 doi: 10.11889/j.0253-3219.2020.hjs.43.060003Wang Tianxiong, Zhang Tengfei, Wu Haicheng, et al. Modeling and analysis of depletion experiment benchmark based on assembly calculation[J]. Nuclear Techniques, 2020, 43: 060003 doi: 10.11889/j.0253-3219.2020.hjs.43.060003 [10] 肖越, 吴海成, 吴小飞, 等. 燃耗数据库基准检验方法研究[J]. 原子能科学技术, 2022, 56(5):952-960Xiao Yue, Wu Haicheng, Wu Xiaofei, et al. Research on benchmarking method of burnup database[J]. Atomic Energy Science and Technology, 2022, 56(5): 952-960 [11] 杨森涵, 李云召, 邵睿智, 等. 基于NECP-Bamboo程序的商用压水堆乏燃料组件核素成分分析[J/OL]. 原子能科学技术: 1-10[2023-01-12]. https://kns.cnki.net/KCMS/detail/detail.aspx?dbcode=CAPJ&dbname=CAPJLAST&filename=YZJS20221108000.Yang Senhan, Li Yunzhao, Shao Ruizhi, et al. Nuclide composition evaluation for commercial PWR nuclear spent fuel assembly based on NECP-Bamboo[J/OL]. Atomic Energy Science and Technology: 1-10[2023-01-12]. https://kns.cnki.net/KCMS/detail/detail.aspx?dbcode=CAPJ&dbname=CAPJLAST&filename=YZJS20221108000. [12] Ade B J. SCALE/TRITON Primer: a primer for light water reactor lattice physics calculations[R]. NUREG/CR-7041, ORNL/TM-2011/21, 2012. [13] Greene N M. BONAMI: resonance self-shielding by the Bondarenko method[R]. ORNL/TM-2005/39, 2011. [14] Westfall R M, Petrie L M, Greene N M, et al. NITAWL: SCALE system module for performing resonance shielding and working library production[R]. ORNL/TM-2005/39, 2011. [15] Scaglione J M. Three Mile Island unit 1 radiochemical assay comparisons to SAS2H calculations[R]. CAL-UDC-NU-000011, 2002. [16] Gauld I C. ORIGEN-S: depletion module to calculate neutron activation, actinide transmutation, fission product generation, and radiation source terms[R]. ORNL/TM-2005/39, 2011. [17] Wimmer L B. Summary report of commercial reactor criticality data for Three Mile Island unit 1[R]. TDR-UDC-NU-000004, 2001. [18] 夏兆东, 周小平, 李晓波, 等. 田湾核电站乏燃料水池采用燃耗信任制的计算研究[J]. 原子能科学技术, 2013, 47(11):2098-2102 doi: 10.7538/yzk.2013.47.11.2098Xia Zhaodong, Zhou Xiaoping, Li Xiaobo, et al. Calculation study of TNPS spent fuel pool using burnup credit[J]. Atomic Energy Science and Technology, 2013, 47(11): 2098-2102 doi: 10.7538/yzk.2013.47.11.2098 [19] ASTM E321-2, Standard test method for atom percent fission in uranium and plutonium fuel (neodymium-148 method)[S]. 期刊类型引用(3)
1. 赵言昊,饶波,杨勇,俞豪俊,张明. 氢闸流管工作特性研究. 强激光与粒子束. 2024(05): 79-83 . 本站查看
2. 谌怡,黄子平,张篁,刘毅,丁明军,夏连胜. 氢闸流管驱动三同轴电缆Blumlein线的kHz重频脉冲功率源. 强激光与粒子束. 2024(05): 66-71 . 本站查看
3. 黄明阳,许守彦,齐欣,王生. 中国散裂中子源加速器束流引出调试和束损优化. 强激光与粒子束. 2023(12): 52-59 . 本站查看
其他类型引用(0)
-