首页 / 专利库 / 物理 / 临界状态 / 核反应堆达临界状态模拟方法

核反应堆临界状态模拟方法

阅读:825发布:2020-05-12

专利汇可以提供核反应堆临界状态模拟方法专利检索,专利查询,专利分析的服务。并且本 发明 属于 核反应堆 领域,特别是一种核反应堆 临界状态 模拟方法,包括以下步骤:(1) 硼 化稀释模型得到硼浓度变化量;(2)棒位控 制模 型的叠步棒位与单棒换算式得到叠步棒位变化量;(3)硼浓度变化量和叠步棒位变化量代入 反应性 模型得到反应性变化量;(4)由反应堆模型得到 中子 数;(5)当改变硼浓度和 控制棒 棒位时中子数不变,则反应堆达临界状态。使用本发明的方法:能够动态直观实时地表现出中子数与反应性的关系;由于采用基于计算机数值解法,对于反应性任何形式的变化,中子数均能被非常精确的计算出;所需资源少,只需一台计算机即可实现;可以灵活地更改参数,以适应不同结构的反应堆的达临界特性。,下面是核反应堆临界状态模拟方法专利的具体信息内容。

1.一种核反应堆临界状态模拟方法,包括以下步骤:
(1)对反应堆进行化或稀释操作,得到硼化或稀释引入的硼浓度c1、硼化或稀释体积v、反应堆的容积v0,将参数c1、v、v0代入硼化稀释方程,得到反应堆的硼浓度c变化量dc;
所述的硼化稀释方程为:
c1dv=v0dc+cdv
式中:
c——反应堆的硼浓度
-3
c1——由硼化或稀释引入的硼浓度,当为硼化状态时,c1的值为7×10 ,当为稀释状态时为0
v0——反应堆的水容积,该值与反应堆结构有关,对同一反应堆来讲,在反应堆运行时该值近似不变
v——硼化或稀释体积,该值取决于反应堆现有硼浓度和其与所要达到的硼浓度的差,无具体范围,取值范围可以很大,若不考虑时间性经济性,理论值可以在零至无穷大(2)对反应堆进行控制棒位的提升或下降操作,当进行提棒控制时,增加了叠步棒位数l的值;插棒控制时,减少了叠步棒位数l的值,由棒位控制模型的叠步棒位与单棒的换算关系式,计算得到叠步棒位l的变化量dl;
所述的叠步棒位与单棒的换算关系式为:
当反应堆有i组控制棒,最高棒位为a步,最低棒位为b步时,所述的棒位控制模型的叠步棒位与单棒的换算关系式如下式:
式中:
l-叠步棒位数
Ni-i棒棒位数
x-重叠步数
a-最高棒位步数
b-最低棒位步数
(3)将步骤(1)中得到的硼浓度变化量dc和步骤(2)中得到的叠步棒位变化量dl代入反应性计算模型,得到反应性的变化量dρ;
所述的反应性计算模型为:
dρ=k1dc+k2dl
式中:
dρ-反应性的变化量
k1-硼浓度变化量影响反应性变化量的比例系数
k2-棒位变化量影响反应性变化量的比例系数
dc-反应堆的硼浓度变化量
dl-叠步棒位变化量
所述的硼浓度变化量影响反应性变化量的比例系数k1可以为常 数或非线性系数;k1-5
为常数时其范围为-70×10 ~0;k1的值取决于反应堆硼浓度、反应堆冷却剂温度、反应堆结构、核燃料的燃耗、控制棒的位置;k1为非线性系数时,使用二次曲线拟合法或三次曲线拟合法对系数k1进行非线性化处理;
所述的棒位变化量影响反应性变化量的比例系数k2均可以为常数或非线性系数;k2为-5
常数时其范围为0~10×10 ,k2的值取决于反应堆硼浓度、反应堆冷却剂温度、反应堆结构、核燃料的燃耗、控制棒的位置;k2为非线性系数时,使用二次曲线拟合法或三次曲线拟合法对系数k2进行非线性化处理;
(4)设定初始中子数n、缓发中子份额β、衰变常数λ、平均代时间l、外加中子源强q、反应性ρ,反应性ρ是由反应堆模型初始设定的反应性ρ0与反应性计算模型得到的反应性变化量dρ相加构成的,即:ρ=ρ0+dρ;
将上述所有设定的参数和反应性ρ带入反应堆模型得到单位时间内中子数的变化量 和先驱核数的变化率
所述的反应堆模型是按照点堆动学方程建立的,反应堆模型为:
n——中子数
i——由于核反应会产生多种先驱核类型,各先驱核类型的衰变时间不同,存在6种主要类型的先驱 核,故i=1,2,…6
βi——缓发中子份额β是由先驱核衰变产生的中子数,由于存在6种主要类型的先驱核,故有βi
β——缓发中子份额,缓发中子份额是由先驱核衰变产生的中子数占核反应产生所有的中子数的比例,一般为0.0065,
λ——衰变常数,其范围为0.01~10每秒
l——平均代时间,其范围为0.0849秒
q——外加中子源强,取决于反应堆类型,无确定值,有的反应堆无外加中子源强,即取q=0
ρ——反应性
ci——先驱核数,它取决于反应堆功率,理论值可以在零至无穷大,因实际反应堆功率有限,先驱核数ci不可能无穷大
——单位时间内中子数的变化量
——所有中子均为快中子时产生的中子数
——因缓发而被扣除的中子总数
——各先驱核产生的缓发中子总数
——先驱核数的变化率
n——中子数
(5)步骤(1)中的硼浓度c变化时,硼浓度变化量dc随之变化,引起反应性的变化量dρ变化,从而改变了反应堆模型中反应性ρ的 值,因此中子数n发生变化;
步骤(2)中控制棒棒位改变时棒位变化量dl随之变化,引起反应性的变化量dρ变化,从而改变了反应堆模型中反应性ρ的值,因此中子数n发生变化;
当改变硼浓度和控制棒棒位而中子数n不发生变化,即单位时间内中子数的变化量时,核反应堆达临界状态;
当反应堆达临界状态时硼化稀释操作停止。

说明书全文

核反应堆临界状态模拟方法

技术领域

[0001] 本发明涉及一种核反应堆模拟方法,特别是一种核反应堆达临界状态模拟方法。

背景技术

[0002] 核反应堆达到临界,即将反应堆的中子数稳定在某一个值不变化,可以使核反应过程受控,这样可以利用这种稳定的核反应所释放出的能量来发电。
[0003] 使核反应堆达到临界状态是任何一座核反应堆操作人员所必需掌握的技能,现有的核反应堆达到临界状态是通过大型模拟机模拟操作实现的。使用大型模拟机模拟核反应堆达到临界状态存在许多缺点:对人物力要求高,需专人配合;需要日常维护;缺乏灵活性,不能随时更改参数用于理论研究;大型模拟机地点固定、数量有限且模拟周期长。

发明内容

[0004] 本发明的目的在于提供一种核反应堆达临界状态模拟方法。
[0005] 一种核反应堆达临界状态模拟方法,包括以下步骤:
[0006] (1)对反应堆进行化或稀释操作,得到硼化或稀释引入的硼浓度c1、硼化或稀释体积v、反应堆的容积v0,将参数c1、v、v0代入硼化稀释方程,得到反应堆的硼浓度c变化量dc;
[0007] 所述的硼化稀释方程为:
[0008] c1dv=v0dc+cdv
[0009] 式中:
[0010] c——反应堆的硼浓度
[0011] c1——由硼化或稀释引入的硼浓度,当为硼化状态时,c1的值为7×10-3,当为稀释状态时为0
[0012] v0——反应堆的水容积,该值与反应堆结构有关,对同一反应堆来讲,在反应堆运行时该值近似不变
[0013] v——硼化或稀释体积,该值取决于反应堆现有硼浓度和其与所要达到的硼浓度的差,无具体范围,取值范围可以很大,若不考虑时间性经济性,理论值可以在零至无穷大[0014] (2)对反应堆进行控制棒位的提升或下降操作,当进行提棒控制时,增加了叠步棒位数l的值;插棒控制时,减少了叠步棒位数l的值,由棒位控制模型的叠步棒位与单棒的换算关系式,计算得到叠步棒位l的变化量dl;
[0015] 所述的叠步棒位与单棒的换算关系式为:
[0016] 当反应堆有i组控制棒,最高棒位为a步,最低棒位为b步时,所述的棒位控制模型的叠步棒位与单棒的换算关系式如下式:
[0017]
[0018] 式中:
[0019] l-叠步棒位数
[0020] Ni-i棒棒位数
[0021] x-重叠步数
[0022] a-最高棒位步数
[0023] b-最低棒位步数。
[0024] (3)将步骤(1)中得到的硼浓度变化量dc和步骤(2)中得到的叠步棒位变化量dl代入反应性计算模型,得到反应性的变化量dρ;
[0025] 所述的反应性计算模型为:
[0026] dρ=k1dc+k2dl
[0027] 式中:
[0028] dρ-反应性的变化量
[0029] k1-硼浓度变化量影响反应性变化量的比例系数
[0030] k2-棒位变化量影响反应性变化量的比例系数
[0031] dc-反应堆的硼浓度变化量
[0032] dl-叠步棒位变化量;
[0033] 所述的硼浓度变化量影响反应性变化量的比例系数k1可以为常数或非线性系数;-5
k1为常数时其范围为-70×10 ~0;k1的值取决于反应堆硼浓度、反应堆冷却剂温度、反应堆结构、核燃料的燃耗、控制棒的位置;k1为非线性系数时,使用二次曲线拟合法或三次曲线拟合法对系数k1进行非线性化处理。
[0034] 所述的棒位变化量影响反应性变化量的比例系数k2均可以为常数或非线性系数;-5
k2为常数时其范围为0~10×10 ,k2的值取决于反应堆硼浓度、反应堆冷却剂温度、反应堆结构、核燃料的燃耗、控制棒的位置;k2为非线性系数时,使用二次曲线拟合法或三次曲线拟合法对系数k2进行非线性化处理。
[0035] (4)设定初始中子数n、缓发中子份额β、衰变常数λ、平均代时间l、外加中子源强q、反应性ρ,反应性ρ是由反应堆模型初始设定的反应性ρ0与反应性计算模型得到的反应性变化量dρ相加构成的,即:ρ=ρ0+dρ;
[0036] 将上述所有设定的参数和反应性ρ带入反应堆模型得到单位时间内中子数的变化量 和先驱核数的变化率
[0037] 所述的反应堆模型是按照点堆动力学方程建立的,反应堆模型为:
[0038]
[0039] n——中子数
[0040] i——由于核反应会产生多种先驱核类型,各先驱核类型的衰变时间不同,存在6种主要类型的先驱核,故i=1,2,…6
[0041] βi——缓发中子份额β是由先驱核衰变产生的中子数,由于存在6种主要类型的先驱核,故有βi
[0042] β——缓发中子份额,缓发中子份额是由先驱核衰变产生的中子数占核反应产生所有的中子数的比例,
[0043] 一般为0.0065,
[0044] λ——衰变常数,其范围为0.01~10每秒
[0045] l——平均代时间,其范围为0.0849秒
[0046] q——外加中子源强,取决于反应堆类型,无确定值,有的反应堆无外加中子源强,即取q=0
[0047] ρ——反应性
[0048] ci——先驱核数,它取决于反应堆功率,理论值可以在零至无穷大,因实际反应堆功率有限,先驱核数ci不可能无穷大
[0049] ——单位时间内中子数的变化量
[0050] ——所有中子均为快中子时产生的中子数
[0051] ——因缓发而被扣除的中子总数
[0052] ——各先驱核产生的缓发中子总数
[0053] ——先驱核数的变化率
[0054] n——中子数
[0055] (5)步骤(1)中的硼浓度c变化时,硼浓度变化量dc随之变化,引起反应性的变化量dρ变化,从而改变了反应堆模型中反应性ρ的值,因此中子数n发生变化;
[0056] 步骤(2)中控制棒棒位改变时棒位变化量dl随之变化,引起反应性的变化量dρ变化,从而改变了反应堆模型中反应性ρ的值,因此中子数n发生变化;
[0057] 当改变硼浓度和控制棒棒位而中子数n不发生变化,即单位时间内中子数的变化量 时,核反应堆达临界状态;
[0058] 当反应堆达临界状态时硼化稀释操作停止。
[0059] 使用本发明的方法:能够动态直观实时地表现出中子数与反应性的关系;由于采用基于计算机的数值解法,对于反应性任何形式的变化,中子数均能被非常精确的计算出;所需资源少,只需一台计算机即可实现;可以灵活地更改参数,以适应不同结构的反应堆的达临界特性。

具体实施方式

[0060] 实施例1
[0061] (1)对反应堆进行硼化稀释操作,得到硼化或稀释引入的硼浓度c1、硼化或稀释体积v、反应堆的水容积v0,将参数c1、v、v0代入硼化稀释方程,硼化稀释方程如下:
[0062] c1dv=v0dc+cdv
[0063] 式中:
[0064] c——反应堆的硼浓度
[0065] c1——由硼化或稀释引入的硼浓度,当为硼化状态时,c1的值为7×10-3,当为稀释状态时为0
[0066] v0——反应堆的水容积,该值与反应堆结构有关,对同一反应堆来讲,在反应堆运行时该值近似不变
[0067] v——硼化或稀释体积,该值取决于反应堆现有硼浓度和其与所要达到的硼浓度的差,无具体范围,取值范围可以很大,若不考虑时间性经济性,理论值可以在零至无穷大[0068] 由上式得到反应堆的硼浓度c及反应堆的硼浓度的变化量dc,当反应堆达临界状态时硼化稀释操作停止。
[0069] (2)进行控制棒位的提升或下降,当进行提棒控制时,增加了叠步棒位数l的值;插棒控制时,减少了叠步棒位数l的值;反应堆有4组控制棒,其最高棒位为225步,最低棒位为5步,叠步棒位与单棒的换算关系式如下:
[0070]
[0071] 式中:
[0072] l-叠步棒位数
[0073] N1-1棒棒位数
[0074] N2-2棒棒位数
[0075] N3-3棒棒位数
[0076] N4-4棒棒位数
[0077] x-重叠步数
[0078] 从而由上式计算出1、2、3、4棒、叠步棒位l的数值,由叠步棒位l的数值变化得出叠步棒位l的变化量dl。
[0079] (3)将步骤(1)中得到的硼浓度变化量dc和步骤(2)中得到的棒位变化量dl代入反应性计算模型;由于达临界过程中,反应堆内冷却剂温度变化量微小,因此在考虑对反应性影响因素时,忽略了因冷却剂温度变化对反应性的影响,考虑了棒位与硼浓度变化时对反应性的影响;反应性计算模型如下式:
[0080] dρ=k1dc+k2dl
[0081] 式中:
[0082] dρ-反应性的变化量
[0083] k1-硼浓度变化量影响反应性变化量的比例系数
[0084] k2-棒位变化量影响反应性变化量的比例系数
[0085] dc-反应堆的硼浓度变化量
[0086] dl-棒位变化量
[0087] 由上式得出反应性的变化量dρ。
[0088] 将棒位与硼浓度对反应性的影响线性化,即上式中k1、k2系数为常数,k1的范围-5 -5为-70×10 ~0,k2的范围为0~10×10 。如需要更为精确的数值,而且对于不同的反应堆类型或同一反应堆不同时期k1、k2的值不同,k1、k2的值取决于反应堆硼浓度、反应堆冷却剂温度、反应堆结构、核燃料的燃耗、控制棒的位置。使用二次曲线拟合法或三次曲线拟合法对系数k1、k2进行非线性化处理。
[0089] (4)设定中子数n、缓发中子份额β、衰变常数λ、平均代时间l、外加中子源强q、反应性ρ,反应性ρ是由反应堆模型初始设定的反应性ρ0与反应性计算模型得到的反应性变化量dρ相加构成的,即:ρ=ρ0+dρ,其中ρ0在一次模拟过程中作为初始设定是常数,对于不同次的模拟过程ρ0为不同的常数;将上述所有设定的参数和反应性ρ带入反应堆模型,反应堆模型是按照点堆动力学方程建立的,反应堆模型如下式:
[0090]
[0091] 式中:
[0092] n——中子数
[0093] i——由于核反应会产生多种先驱核类型,各先驱核类型的衰变时间不同,存在6种主要类型的先驱核,故i=1,2,…6
[0094] βi——缓发中子份额β是由先驱核衰变产生的中子数,由于存在6种主要类型的先驱核,故有βi
[0095] β——缓发中子份额,缓发中子份额是由先驱核衰变产生的中子数占核反应产生所有的中子数的比例,一般为0.0065,
[0096] λ——衰变常数,其范围为0.01~10每秒
[0097] l——平均代时间,其范围为0.0849秒
[0098] q——外加中子源强,取决于反应堆类型,无确定值;有的反应堆无外加中子源强,即取q=0
[0099] ρ——反应性
[0100] ci——先驱核数,它取决于反应堆的功率,理论值可以在零至无穷大,,因实际反应堆功率有限,先驱核数ci不可能无穷大
[0101] ——单位时间内中子数的变化量
[0102] ——所有中子均为快中子时产生的中子数
[0103] ——因缓发而被扣除的中子总数
[0104] ——各先驱核产生的缓发中子总数
[0105] ——先驱核数的变化率
[0106] 由上式得出单位时间内中子数的变化量 和先驱核数的变化率
[0107] (5)步骤(1)中的硼浓度的变化时硼浓度变化量dc随之变化,引起反应性的变化量dρ变化,从而改变了反应堆模型中反应性ρ的值,因此中子数n发生变化;步骤(2)中控制棒棒位改变时棒位变化量dl随之变化,引起反应性的变化量dρ变化,从而改变了反应堆模型中反应性ρ的值,因此中子数n发生变化。
[0108] 假设在时刻t时中子数为n,那么在时刻t+dt时,由反应堆模型计算出的中子数为n′,n′=n时核反应堆达临界状态;即当改变硼浓度和控制棒棒位时,单位时间内中子数的变化量 则核反应堆达临界状态。
[0109] 当反应堆达临界状态时,由于中子数n不变,即 而缓发中子份额β为常数,故缓发中子数 不变,又由于先驱核衰变常数λ不变,故作为缓发中子的先驱核数ci不变,因此
[0110] 将单位时间内中子数的变化量 带入中子周期T的公式:
[0111]
[0112] 式中:
[0113] T——反应堆中子周期
[0114] n——中子数
[0115] 单位时间内中子数的变化量 为零时核反应堆达临界状态,由上可知核反应堆达临界状态时中子周期为无穷大。
[0116] 实施例2
[0117] 与实施例1的不同在于步骤(2)中的反应堆有3组控制棒,其最高棒位为225步,最低棒位为5步,叠步棒位与单棒的换算关系式如下:
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈