一种圆柱体铀溶液贮槽内核临界事故裂变次数估算的方法技术

本发明专利技术涉及一种圆柱体铀溶液贮槽内核临界事故裂变次数估算的方法。采用本发明专利技术所提供的圆柱体铀溶液贮槽内核临界事故裂变次数估算的方法，可以根据铀溶液贮槽几何结构参数、铀溶液贮槽内容物的材料成分参数、以及工艺流程参数，利用三维蒙卡中子输运程序建模并计算得到发生临界事故时加入系统的反应性速率ρ

【技术实现步骤摘要】
一种圆柱体铀溶液贮槽内核临界事故裂变次数估算的方法


[0001]本专利技术属于核与辐射安全
，涉及一种圆柱体铀溶液贮槽内核临界事故裂变次数估算的方法。

技术介绍

[0002]GB15146.2
‑
2008中给出与化学组成无关的
233
U、
235
U和
239
Pu的次临界质量限值分别为500g、700g和450g，当核燃料循环设施中处理的易裂变物质的量超过了该规定的量时就存在发生核临界事故的风险。世界各国对于核临界事故的计算和模拟开展了大量研究和实验，例如法国先后建立了CRAC和SILENE实验装置，日本也建造了TRACY实验装置开展一系列的临界实验。在此基础上一些国家也开发了具有动力学计算功能的用于估算核临界事故相关参数的计算程序。
[0003]如果可以明确的定义发生核临界事故时的事故情景和参数，使用这些程序可以较为准确的计算出核临界事故的裂变次数、功率、温度、压力等参数的变化。但是一则用于估算核临界事故相关参数的计算程序并不普遍，二则对于事故情景的定义很多情况下无法明确。因而国际上也有一些经验公式用于估算裂变次数。但是用经验公式估算的裂变次数往往较为包络和保守，用经验公式的值进行设计可能会使事故缓解的代价较大。
[0004]虽然用于核临界事故模拟和分析的动力学计算程序并不普遍，但是稳态的中子输运程序却是广泛使用的，这些程序可以用于计算系统的一些特性参数，比如系统的有效缓发中子份额，系统有效中子增殖因子等，这些参数也是动力学程序进行核临界事故裂变次数的计算的基础参数。

技术实现思路

[0005]针对现有技术中所存在的缺陷，本专利技术的目的在于提供一种圆柱体铀溶液贮槽内核临界事故裂变次数估算的方法，用以替代需要使用动力学计算程序进行裂变次数计算或仅使用经验公式进行估算的方式，实现了仅利用三维蒙卡中子输运程序的计算结果结合前期实验数据拟合得到的参数变化趋势来估算核临界事故的裂变次数。
[0006]为实现此目的，本专利技术提供一种圆柱体铀溶液贮槽内核临界事故裂变次数估算的方法，所述方法包括如下步骤：
[0007](1)确定所述圆柱体铀溶液贮槽的几何结构参数、所述圆柱体铀溶液贮槽内容物的材料成分参数、以及工艺流程参数；
[0008](2)根据步骤(1)中所确定的参数建模，计算搜索贮槽内料液系统的临界液位高度，并计算贮槽内料液系统的中子动力学参数、贮槽内料液不同高度时的系统有效中子增殖因子k
eff
和对应的系统反应性ρ；然后通过计算确定所述圆柱体铀溶液贮槽内发生临界事故时加入系统的反应性速率ρ
t
；
[0009](3)根据所述加入系统的反应性速率ρ
t
与裂变次数体密度N
V
的对应关系，计算得到所述圆柱体铀溶液贮槽内发生临界事故时单位体积内的裂变次数；再结合发生临界事故时
所述圆柱体铀溶液的总体积V，计算得到所述圆柱体铀溶液贮槽内发生临界事故时的总裂变次数N。
[0010]进一步，所述步骤(1)中，所述圆柱体铀溶液贮槽的几何结构参数为圆柱体铀溶液贮槽的直径D。
[0011]进一步，所述步骤(1)中，所述圆柱体铀溶液贮槽内容物的材料成分参数包括铀溶液中铀的浓度C、U
‑
235的富集度E。
[0012]进一步，所述步骤(1)中，所述工艺流程参数为料液添加速度上限l。
[0013]进一步，所述步骤(2)中，所述建模是采用三维蒙卡中子输运程序建立计算模型。
[0014]进一步，所述步骤(2)中，所述贮槽内料液系统的中子动力学参数包括系统的有效缓发中子份额β
eff
、衰变常数λ
i
。
[0015]进一步，根据贮槽内料液不同高度时的系统有效中子增殖因子k
eff
和对应的系统反应性ρ＝1
‑
1/k
eff
，得到到达临界附近处时，所述系统反应性ρ随高度H的变化率dρ；根据所述圆柱型溶液贮槽的直径D、料液添加速度上限l，计算得到贮槽内料液高度H变化的最大速率s；然后根据下列关系式：
[0016][0017]计算系统临界时加入系统的反应性速率ρ
t
。
[0018]进一步，所述步骤(3)中，所述加入系统的反应性速率ρ
t
与裂变次数体密度N
V
的对应关系是由前期实验数据拟合得到的关于所述裂变次数体密度N
V
与加入系统的反应性速率ρ
t
的经验公式：
[0019]N
V
＝(1.81ρ
‑
0.645)
×
10
15
[0020]将步骤(2)中得到的所述加入系统的反应性速率ρ
t
代入计算，得到裂变次数体密度N
V
。
[0021]进一步，所述裂变次数体密度N
V
与加入系统的反应性速率ρ
t
的经验公式中，所述加入系统的反应性速率ρ
t
的参数变化范围为0.5
‑
3.0$/s。
[0022]进一步，所述加入系统的反应性速率ρ
t
的参数变化范围为0.8～1.7$/s。
[0023]本专利技术的有益效果在于，采用本专利技术所提供的圆柱体铀溶液贮槽内核临界事故裂变次数估算的方法，可以利用常用的稳态中子输运程序相对合理的估算含有硝酸铀溶液的圆柱体贮槽内发生核临界事故时的裂变次数。采用本专利技术的方法，根据铀溶液贮槽几何结构参数、铀溶液贮槽内容物的材料成分参数、以及工艺流程参数来确定发生核临界事故时裂变次数及事故大小，相比于简单的经验公式，其估算值更为合理。同时，计算仅需利用常用的中子输运程序，无需开发和使用专用于核临界事故模拟的中子动力学计算程序，对核临界事故的估算更为便捷。
附图说明
[0024]图1为本专利技术所述的圆柱体铀溶液贮槽内核临界事故裂变次数估算的方法示意图。
[0025]图2为本专利技术实施方式中圆柱体铀溶液贮槽内核临界事故裂变次数估算的方法示意图。
[0026]图3为本专利技术实施方式中溶液贮槽内系统反应性ρ随贮槽内料液高度H的变化拟合曲线图。
[0027]图4为本专利技术所述的裂变次数体密度随加入系统的反应性速率变化的拟合曲线图。
具体实施方式
[0028]下面结合附图和具体实施方式对本专利技术作进一步描述。
[0029]一种装载了浓度为300gU/L的硝酸铀酰溶液、直径为80cm的圆柱体贮槽核临界事故裂变次数估算的方法，包括如下步骤：
[0030]1)首先利用三维蒙卡中子输运程序进行建模，计算搜索该系统的临界液位高度，并计算系统的缓发中子有效份额β
ieff
和衰变常数λ
i
等中子动力学参数、系统有效中子增殖因子k
e本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种圆柱体铀溶液贮槽内核临界事故裂变次数估算的方法，其特征在于，所述方法包括如下步骤：(1)确定所述圆柱体铀溶液贮槽的几何结构参数、所述圆柱体铀溶液贮槽内容物的材料成分参数、以及工艺流程参数；(2)根据步骤(1)中所确定的参数建模，计算搜索贮槽内料液系统的临界液位高度，并计算贮槽内料液系统的中子动力学参数、贮槽内料液不同高度时的系统有效中子增殖因子k
eff
和对应的系统反应性ρ；然后通过计算确定所述圆柱体铀溶液贮槽内发生临界事故时加入系统的反应性速率ρ
t
；(3)根据所述加入系统的反应性速率ρ
t
与裂变次数体密度N
V
的对应关系，计算得到所述圆柱体铀溶液贮槽内发生临界事故时单位体积内的裂变次数；再结合发生临界事故时所述圆柱体铀溶液的总体积V，计算得到所述圆柱体铀溶液贮槽内发生临界事故时的总裂变次数N。2.根据权利要求1所述的一种圆柱体铀溶液贮槽内核临界事故裂变次数估算的方法，其特征在于，所述步骤(1)中，所述圆柱体铀溶液贮槽的几何结构参数为圆柱体铀溶液贮槽的直径D。3.根据权利要求1所述的一种圆柱体铀溶液贮槽内核临界事故裂变次数估算的方法，其特征在于，所述步骤(1)中，所述圆柱体铀溶液贮槽内容物的材料成分参数包括铀溶液中铀的浓度C、U
‑
235的富集度E。4.根据权利要求2所述的一种圆柱体铀溶液贮槽内核临界事故裂变次数估算的方法，其特征在于，所述步骤(1)中，所述工艺流程参数为料液添加速度上限l。5.根据权利要求1所述的一种圆柱体铀溶液贮槽内核临界事故裂变次数估算的方法，其特征在于，所述步骤(2)中，所述建模是采用三维蒙卡中子输运程序建立计算模型。6.根据权利要求4所述的一种圆柱体铀溶液贮槽内核临界事故裂变次数估算的方法，其特征在于，所述步骤(2)中，所述贮槽内料液系统的中子动力学参数包括系统的有效缓发...

【专利技术属性】
技术研发人员：易璇，于淼，霍小东，邵增，刘国明，杨海峰，郭治鹏，李想，胡小利，陈添，
申请(专利权)人：中国核电工程有限公司，
类型：发明
国别省市：