发明名称 一种反应堆反应性测量值的外源修正方法
摘要 本发明涉及一种反应堆反应性测量值的外源修正方法,包括:假定ΛS(1)、ρ(1);对实测计数序列(tj,nj)进行修正得(tj,n′j);对n′j取自然对数,按ln(n′j)=A′+B′*t关系对t进行线性拟合,得到拟合参数A′和B′;由A′和B′以及倒时方程等得到更新的Λ(S)和ρ(1)值;使更新的ΛS(1)和ρ(1)值代替ΛS(1)、ρ(1),再次对计数序列(tj,nj)进行修正,重复以上步骤,直到ΛS(1)和ρ(1)达到预设的收敛要求。本发明提供的反应堆反应性测量值的外源修正方法,采用拟合迭代法能够很快得到收敛的ΛS和ρ,与现有基于逆动态法的修正结果相差不大,结果有效可靠。
申请公布号 CN102800372A 申请公布日期 2012.11.28
申请号 CN201210295976.7 申请日期 2012.08.17
申请人 中国原子能科学研究院 发明人 胡赟;张坚;赵郁森;周科源;杨晓燕;徐李
分类号 G21C17/104(2006.01)I;G06F19/00(2011.01)I 主分类号 G21C17/104(2006.01)I
代理机构 代理人
主权项 一种反应堆反应性测量值的外源修正方法,包括以下步骤:(Ⅰ)将选择的计数序列(tj,nj)作为拟合迭代处理的原始数据,tj表示时间序列,nj表示对应于tj的测量的计数率;(Ⅱ)假定S、ρ初始大小S(0)、ρ(0),将ΛS(0)、ρ(0)代入式(1)对计数序列(tj,nj)进行修正,得到修正后的计数序列(tj,n′j), <mrow> <msub> <msup> <mi>n</mi> <mo>&prime;</mo> </msup> <mi>j</mi> </msub> <mo>=</mo> <msub> <mi>n</mi> <mi>j</mi> </msub> <mo>+</mo> <mfrac> <mi>&Lambda;S</mi> <mi>&rho;</mi> </mfrac> <mo>=</mo> <mi>c</mi> <mo>&CenterDot;</mo> <msup> <mi>e</mi> <msub> <mi>&omega;t</mi> <mi>j</mi> </msub> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>式(1)中,c表示与反应性和初始计数n0相关的常数,ω表示反应性方程(3)对应于稳定周期的解;(Ⅲ)对修正后的计数n′j取自然对数,然后按ln(n′j)=A′+B′*t关系对时间变量t进行线性拟合,得到拟合参数A′和B′,并根据式(1)得到:A'=lnc    (2.1)B'=ω    (2.2) <mrow> <mi>c</mi> <mo>=</mo> <msub> <mi>n</mi> <mn>0</mn> </msub> <mo>+</mo> <mfrac> <mi>&Lambda;S</mi> <mi>&rho;</mi> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2.3</mn> <mo>)</mo> </mrow> </mrow>式(2.3)中,n0表示所处理的计数率序列初始值(t=0);(Ⅳ)由拟合得到的参数A′和B′,按式(2.1‑2.3)以及式(3)所示的倒时方程得到ΛS、ρ的迭代更新值,计为ΛS(1)、ρ(1), <mrow> <mi>&rho;</mi> <mo>=</mo> <mi>&Lambda;&omega;</mi> <mo>+</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <mfrac> <mrow> <msub> <mi>&beta;</mi> <mi>i</mi> </msub> <mi>&omega;</mi> </mrow> <mrow> <mi>&omega;</mi> <mo>+</mo> <msub> <mi>&lambda;</mi> <mi>i</mi> </msub> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>式(3)中,βi表示第i组有效缓发中子份额,λi表示第i组缓发中子先驱核衰变常数,M为缓发中子先驱核总的分组数;(Ⅴ)以ΛS(1)、ρ(1)代替步骤(Ⅱ)中的ΛS(0)、ρ(0),重复第(Ⅱ)~(Ⅳ)步,直到ΛS、ρ的迭代更新值达到预设的收敛要求。
地址 102413 北京市房山区北京市275信箱65分箱