硼浓度由下式进行计算
式中,B是一回路硼的
质量,单位g;M是一回路硼溶液的质量,单位t;BC 是硼浓度,单位mg/kg。因此
B=M·BC (2)
B10=xB=xM·BC (3)
式中,B10是一回路硼-10的质量,单位g;x是一回路硼-10的富集度(同位素 质量含量)。
如果不考虑硼的核消耗效应,并假设一回路排出的冷却剂硼浓度与一回 路的硼浓度相同,则在硼化稀释过程中,一回路硼浓度的变化量
式中,m是注入(或排出)一回路硼化稀释
水的质量,单位t;BCi表示注入一回 路硼化稀释水的硼浓度,单位mg/kg。对方程(4)进行积分,则有
上式(5)就是常用的一回路硼化稀释公式。式中,BCb表示硼化稀释前 一回路冷却剂的硼浓度,单位mg/kg;BCe表示硼化稀释后一回路冷却剂的硼浓 度,单位mg/kg。
常用的一回路硼化稀释公式(5)在工程应用上具有足够的
精度,但由于没 有考虑硼的核消耗效应,它只反映出一回路硼化稀释对硼浓度的影响(或硼质 量的影响),而没有反映出反应堆功率对硼浓度的影响(或硼质量的影响),也 就无法反映出10B丰度变化的影响。
另外,反应堆一回路可溶硼硼-10丰度变化较小且与功率运行史有密切关 系,因此传统的
堆芯反应性跟踪与核设计中不考虑它的变化。由此在实际堆 芯跟踪过程中,会有以下不利后果:
(1)无法真实地反映出理论计算临界硼浓度与实际跟踪的差别,对核设 计评价不够真实;
(2)反应堆的后备反应性的大小估计产生偏差,进而影响反应堆停堆日 期的准确预测;
(3)对硼微分价值随燃耗的变化估计不足,进而影响寿期末慢化剂
温度 系数的测量结果的准确性和可信度;
(4)产生“难以解释”的现象。
本发明的目的是:提供一种可以计算反应堆一回路中可溶硼硼-10丰度、 10B核消耗量以及10B核消耗产生的
能量的一种计算方法。
本发明的技术方案是:一种反应堆一回路可溶硼硼-10丰度的跟踪计算方 法,它包括以下步骤:
a)读入反应堆参数(d1)与反应堆功率运行史数据(d2),进行参数初始化 及预处理;
b)根据反应堆功率运行史两个状态点分析硼浓度的数据(d2),判断出从 状态点1到状态点2是硼化操作还是稀释操作,从而确定注入一回路硼化稀 释水的硼浓度BCi;
c)根据状态点1和2的燃耗数据,把反应堆参数按燃耗区间进行插值, 得到当前燃耗下的反应堆参数;
d)假设状态点1硼-10丰度是已知的(即假设第一个状态点1硼-10丰度 为自然丰度),并同时假设状态点2硼-10丰度等于状态点1的硼-10丰度;
e)根据计算了10B核消耗的一回路硼化稀释公式和一回路
泄漏处理模型 计算第二个状态点10B富集度x2,并进行
迭代计算,直到x2收敛,所说的硼化 稀释公式为:
其中,t1为状态1点的时间,t2为状态2点的时间,x1为t1时刻的10B同位 素富集度,x2为t2时刻的10B同位素富集度,BC1为t1时刻的硼浓度,BC2为t2 时刻的硼浓度,xi为注入一回路冷却剂的硼-10富集度,M为一回路硼溶液的 质量,
m为注入(或排出)一回路硼化稀释水的质量,t为时间,
σa为10B
中子微观吸收截面,φ为堆芯活性区中子通量,V为 一回路
水体积,v为活性区水体积;
f)将10B富集度换算成丰度,并把计算结果输出;
g)为了进行下一步计算,保存x2值,即令第一个状态点10B富集度x1等 于第二个状态点10B富集度x2;
h)检查反应堆功率运行史的下一个数据,如果
数据处理未结束,则继续 循环步骤b)至步骤g)的处理;如果数据处理结束,则程序结束。
本发明的优点是:根据实测数据,能够较准确地计算出一回路系统内的 10B丰度、10B核消耗量、10B核消耗产生的能量;根据核设计硼降数据,能够 计算出一回路系统内的10B丰度变化及硼浓度变化,预计的10B核消耗量、10B 核消耗产生的能量等。同时能够解释反应堆跟踪过程中与硼-10丰度有关的现 象,如在实际跟踪过程中,参考硼浓度与修正硼浓度的差值之间的拱形曲线 以及在反应堆停堆小修前后硼浓度发生了明显的不连续等现象,并能更真实 的反映了理论计算临界硼浓度与实际跟踪的差别。
附图说明
图1是反应堆一回路可溶硼硼-10丰度的跟踪计算数据处理简要
流程图;
图2是某A循环堆芯硼浓度跟踪曲线图;
图3是对图2做了硼-10核消耗效应修正的堆芯硼浓度跟踪曲线图;
图4是某B循环堆芯硼浓度跟踪曲线;
图5是对图4做了硼-10核消耗效应修正的堆芯硼浓度跟踪曲线图;
图6是由程序abdc计算得到的某B循环硼-10丰度跟踪曲线图。
本发明是这样实现的,一种反应堆一回路可溶硼硼-10丰度的跟踪计算方 法,它包括以下步骤:
a)读入反应堆参数(d1)与反应堆功率运行史数据(d2),进行参数初始化 及预处理;
b)根据反应堆功率运行史两个状态点分析硼浓度的数据(d2),判断出从 状态点1到状态点2是硼化操作还是稀释操作,从而确定注入一回路硼化稀 释水的硼浓度BCi;
c)根据状态点1和2的燃耗数据,把反应堆参数按燃耗区间进行插值, 得到当前燃耗下的反应堆参数;
d)假设状态点1硼-10丰度是已知的(即假设第一个状态点1硼-10丰度 为自然丰度),并同时假设状态点2硼-10丰度等于状态点1的硼-10丰度;
e)根据计算了10B核消耗的一回路硼化稀释公式和一回路泄漏处理模型 计算第二个状态点10B富集度x2,并进行迭代计算,直到x2收敛,所说的硼化 稀释公式为:
其中,t1为状态1点的时间,t2为状态2点的时间,x1为t1时刻的10B同位 素富集度,x2为t2时刻的10B同位素富集度,BC1为t1时刻的硼浓度,BC2为t2 时刻的硼浓度,xi为注入一回路冷却剂的硼-10富集度,M为一回路硼溶液的 质量,
m为注入(或排出)一回路硼化稀释水的质量,t为时间,
σa为10B中子微观吸收截面,φ为堆芯活性区中子通量,V为 一回路水体积,v为活性区水体积;
f)将10B富集度换算成丰度,并把计算结果输出;
g)为了进行下一步计算,保存x2值,即令第一个状态点10B富集度x1等 于第二个状态点10B富集度x2;
h)检查反应堆功率运行史的下一个数据,如果数据处理未结束,则继续 循环步骤b)至步骤g)的处理;如果数据处理结束,则程序结束。
在步骤e)中考虑了由于硼的核消耗效应具体方法是:首先假设一回路 排出的冷却剂硼浓度与一回路的硼浓度相同,那么在一回路中,10B的质量满 足以下方程:
10B的核反应率为σaNφ,单位cm-3·s-1。σa是10B中子微观吸收截面, 10-24cm2;φ是堆芯活性区中子通量,单位n·cm-2·s-1;N是10B的核
密度,单 位cm-3,由下式计算
式中,V是一回路水体积,单位m3。
单位时间10B的消耗数量为σaNφv×106,单位s-1。v是活性区水体积, 单位m3。其折合摩尔数为
单位mol/s。则单位时间核反应消耗的 10B质量为
单位g/s。将式(7)代入,则单位时间核反应消耗的 10B质量为,单位g/s
因此式(6)可以表示为:
式中,B10 in是进入一回路硼-10的质量,单位g;B10 out是排出一回路硼-10的质量, 单位g。
根据混合均匀假设,有
则式(9)变为
式中,xi是注入一回路冷却剂的硼-10富集度。令流量
则
令
并积分,得
上式就是考虑了10B核消耗效应的一回路硼化稀释公式。已知t1时刻的10B 同位素质量含量(富集度)x1、硼浓度BC1、t2时刻硼浓度BC2,则可由上式解得 t2时刻的10B富集度
如果BCi=0,即系统在稀释,则上式简化为
如果σa=0,即系统不考虑硼-10核消耗效应,则式(14)最后可以化简成 式(5)。
最后可由下式将富集度转换成丰度
式中,y1是10B丰度;x1是10B富集度;A1是10B
原子质量数;A2是11B原子质 量数。
同时,在步骤e)中考虑了一回路泄漏处理模型计算第二个状态点10B的 富集度。具体计算如下:
在循环寿期初由于硼浓度很高,所需要的稀释水量较少,该值有可能小 于一回路泄漏水量,因此为了尽量不对反应堆的反应性产生额外的影响,还 需要配出与一回路相同硼浓度的水向一回路补充(即假设补充水的硼浓度与 一回路冷却剂硼浓度相同)。
根据两个时刻的分析硼浓度,用式(5)可计算出需要的硼化或稀释水量为 m。
已知一回路冷却剂泄漏率为Qleak,则一回路冷却剂泄漏量为mleak
mleak=Qleak(t2-t1)
如果m<mleak,(或Q<Qleak)则需要补充的水量为
mcomplement=mleak-m
式中mcomplement为补充水的质量,则需要补充水的平均流量为
式中Qcomplement为补充水的平均流量,由于式(15)计算中已经把10B的核消 耗效应考虑进去,因此补充计算不再考虑这个效应,即Tm=Qcomplement。
另外根据假设BCi=BC1=BC2,则由式(15)
如果m>=mleak,即不需要额外补充水量,因此不需要做一回路泄漏率影响 的处理。
另外由于硼-10核消耗效应除能引起硼-10丰度的变化外,同时还能产生 一定的能量。硼-10核反应式如下:
10B+n→7Li+4He+2.79MeV(6.1%)
10B+n→7Li*+4He+2.31MeV(93.9%)
7Li*→7Li+γ+0.48MeV
其一次反应释放的能量为2.79Mev
10B核消耗的质量为mB10,其质量数为AB10,则10B消耗的数量为
释放的能量为
上式(19)即为计算10B核消耗所产生能量的公式。
由于本发明的计算方法数据处理量较大,因此根据计算了10B核消耗的一 回路硼化稀释公式,编制了相应的
计算机程序abdc,程序的计算流程图如图 1所示。图中的矩形框代表处理过程,四边形代表数据,菱形代表条件判断。 该图简要的表达了反应堆一回路可溶硼硼-10丰度的跟踪计算方法和数据处 理过程。
abdc程序有如下功能:
(1)可以比较方便、自由地设定反应堆参数。
(2)读入每天的硼浓度分析值与反应堆功率数据。
(3)可以根据燃耗对反应堆参数进行插值计算,获得当前燃耗下的反应 堆参数。
(4)根据堆芯跟踪实测数据计算10B的丰度变化,每天的10B核消耗量及 累计核消耗量。
(5)根据核设计的理论硼降数据计算10B的丰度变化、硼浓度变化、每 天的10B核消耗量及累计核消耗量,得到一个包络的核设计理论硼降曲线。
(6)可以计算10B核反应释放的能量。
另外利用本发明中的反应堆一回路可溶硼硼-10丰度的跟踪计算方法还 可以解释在堆芯跟踪过程中发现的一些“难以解释”的现象。
(1)正常运行过程中的反应性跟踪
图2是某A循环堆芯硼浓度跟踪曲线图,其中曲线1是参考硼浓度;曲 线2是修正硼浓度;曲线3是修正硼浓度与参考硼浓度差值,即硼浓度偏差。 它反映了堆芯剩余反应性的状况,可以以此预测循环停堆燃耗与停堆日期。 通常把反应堆在满功率、所有
控制棒全提、平衡氙这个
临界状态作为参考状 态。由理论计算的参考状态下的临界硼浓度称为参考硼浓度;而在堆芯监督 过程中,将实际测量得到的反应堆状态修正到参考状态时的硼浓度称为修正 硼浓度。在实际跟踪过程中发现,参考硼浓度与修正硼浓度的差值,即硼浓 度偏差不是一个常数,而总是形成一个拱型曲线。这使靠硼浓度偏差来预测 循环停堆燃耗与停堆日期的方法误差太大。
图3是对图2做了硼-10核消耗效应修正的堆芯硼浓度跟踪曲线图,其中 曲线1是参考硼浓度;曲线4是做了硼-10核消耗效应修正的修正硼浓度;曲 线5是修正硼浓度与参考硼浓度差值,即硼浓度偏差。从图中可以看出,硼 浓度偏差的拱型曲线消失了。这是由于随着燃耗加深,硼-10丰度下降,硼-10 浓度相同的情况下,硼浓度就要更高。在寿期初硼-10丰度下降不大,因此 影响较小;在寿期中硼-10丰度下降较大,同时硼浓度也比较高,与自然丰度 的硼浓度相比,一回路就需要更高的硼浓度来维持临界;而寿期末硼-10丰度 虽然下降最大,但其硼浓度很低,硼-10丰度影响就越来越小。这样就产生了 拱型曲线。
(2)寿期中发生停堆小修的反应性跟踪
图4是某B循环堆芯硼浓度跟踪曲线,其中曲线6是参考硼浓度;曲线 7是修正硼浓度;曲线8是硼浓度偏差;曲线9是在堆芯燃耗7600MWd/tU时 进行了停堆小修,小修前后硼浓度偏差发生了明显的不连续。在堆芯燃耗 7600MWd/tU时进行了停堆小修,小修前后硼浓度偏差发生了明显的不连续, 令人“难以理解”。实际上这是由于停堆小修时,一回路需要较大程度的硼化, 使回路中硼-10丰度基本上接近到自然丰度的水平(如图6),从而使硼浓度下 降。这种停堆小修后硼浓度阶梯式下降的现象,一般也是在寿期中比较明显。
图5是某B循环考虑了硼-10核消耗效应的堆芯硼浓度跟踪曲线图,其中 曲线6是参考硼浓度;曲线10是做了硼-10核消耗效应修正的修正硼浓度; 曲线11是硼浓度偏差及其拟合曲线。此时硼浓度偏差基本上是连续,硼浓度 阶梯式下降没有了。
图6是由abdc程度计算得到的某B循环硼-10丰度跟踪曲线图,其中曲 线12是停堆后立即恢复临界引起的硼-10丰度变化;曲线13是停堆小修引起 的硼-10丰度变化;曲线14是功率变化引起的硼-10丰度变化。
综合以上各图可以得知,停堆小修后硼浓度阶梯式下降的现象从另外一 个方面证实了硼-10核消耗效应的存在,也反映了计算模型与abdc程度计算 的正确性。
本发明主要用于计算主回路中可含有可溶毒物的各种堆型中的反应堆一 回路10B浓度的变化,跟踪计算10B丰度的变化,10B的核消耗量,10B核消耗 所产生的能量,硼微分价值随燃耗的修正等,并能解释反应堆跟踪过程中与 硼-10丰度有关的现象。