次临界反应性测量的空间修正及其应用综述

张 瑜,王释伟,徐琳琳,郑 正

次临界反应性测量的空间修正及其应用综述

张 瑜,王释伟,徐琳琳,郑 正

(山东核电有限公司技术部,山东海阳265116)

次临界下的反应性测量技术有着自身的特点,次临界下控制棒的动作、堆芯的次临界度以及外中子源的存在都会对堆芯中子通量的分布产生影响,因此通常情况下堆芯的次临界度只能“监视”,无法准确测量。在堆芯模拟软件发展的基础上,国外科研人员提出了次临界下点堆模型的空间修正方法,将这种方法用于动态棒价值测量(DRWM),并在此基础上进一步发展了次临界控制棒价值测量(SRWM),这些技术有的已经被国内核电站使用,但是国内对空间修正的原理及方法鲜有介绍。本文针对这种需求,总结概括了国外商用堆次临界反应性测量的基本原理与方法,并结合反应性测量仪表技术,给出了次临界反应性仪的数据处理流程,这对于推进国内商用堆次临界反应性测量的研究和实际应用具有较为重要的意义。

次临界;空间效应;点堆模型;控制棒价值

压水堆核电厂对次临界条件下反应性测量有着实质的需求,例如:堆芯装载过程中的反应性、提棒或稀释达临界时的反应性。次临界状态下的反应性测量,不仅有利于监视堆芯状态,防止意外临界;还能够在达临界的过程中进行堆芯物理相关参数的测量与验证,缩短物理试验时间,提高经济效益[1］。

国外对次临界条件下反应性的测量进行了较为广泛而深刻的研究,发展出源增殖方法(MSA/MSM)、脉冲中子源法(PNS)和噪声法(Noise)等次临界度测量方法,并进行了相关的应用研究[2-5］。上述方法,从次临界中子源的增殖方程出发,寻找源量程探测器计数与堆芯反应性的修正关系,但并没有重视商用堆中外中子源对堆芯动力学的影响。J.Wringt[6］在其论文中介绍了源调制方法(Source Modulation Method),这种方法首先基于一维单群输运的裸堆模型,得到有外中子源的静态方程,然后进行点堆近似,进一步研究了有外中子源的次临界堆芯的动力学方程。J.Wringt得到了一个重要结论是:对次临界的堆芯来说,系统的静态次临界反应性会影响系统的动力学行为。J.Wringt的研究,为商用堆次临界条件下反应性的测量提供了较为完善的理论基础。

在上述研究的基础上,Renato Y.R. Kuramoto[7］提出了发展一种用于压水堆临界预测和控制棒价值评价的数字化反应性仪表(Digital Reactivity Meter),这种仪表不仅可以在达临界的过程中连续的监视次临界度,并预测预期的临界点,还可以根据直接的次临界度的测量来评估控制棒的价值;同时Renato Y. R.Kuramoto也指出了将数字化反应性仪表直接应用于次临界反应性测量的困难:1)次临界情况下,反应堆的功率水平较低,外中子源的影响不可忽略;2)需要对源量程探测器(SR)的计数进行适当的处理以便输入数字化反应性仪表。为了解决上述问题,Renato Y.R. Kuramoto提出了一种修正的中子源增殖(NSM)方法,利用Win Naing[8］定义的三种修正因子对求解反应性方程的外中子源项进行了相关修正,再通过对探测器信号的修正,达到测量次临界反应性的目的。Renato Y.R. Kuramoto所提出的方法成为数字化反应性仪表用于次临界反应性测量的基本原理。

动态棒价值测量方法(DRWM)是商用压水堆次临界状态下反应性测量的典型代表[9］。Y.A.Chao[10］详细介绍了DRWM的测量过程、数据分析原理及空间修正影响。Y.A.Chao提出了用于次临界反应性测量的空间修正的概念,他指出经过合适的空间修正,可以补偿由于堆芯基本模型(配置)的改变(例如:控制棒插入)引起中子通量分布空间变化的影响,通过对探测器的信号应用此修正可以使其满足点堆动力学方程的基本假设。Chao还进一步指出使用3D时空动力学模拟(ANC软件)对整个测量过程进行模拟,进而求出每一状态下的堆芯次临界度,以此来修正由于测量过程中缓发中子的滞后效应对反应性测量引起的影响。

DRWM是在临界后低功率物理实验阶段(5%功率以下)测量反应性的,虽然在测量过程中控制棒的下插,导致堆芯变为次临界状态,但通常情况下,由于单根控制棒的价值较低,堆芯的次临界度较低,因此外中子源对堆芯功率分布的影响较小。Y.A.Chao[1］在此基础上进一步研究了有外中子源影响下的空间修正效应,较为详细地研究了次临界度及外中子源对堆芯动力学的影响。如图1所示,分别在HZP、ARO状态下模拟了无外中子源临界、无外中子源次临界(低于HZP、ARO临界硼浓度100 ppm)及有外中子源次临界(低于HZP、ARO临界硼浓度100 ppm)三种情形,可以看出前两种情况下堆芯的功率分布基本相同,这是由于堆芯的次临界度较小(低于HZP、ARO临界硼浓度100 ppm对应的Keff约为0.992 5),表明对于趋于临界的堆芯来说,静态的堆芯次临界度在无外中子源的情况下,对功率分布的影响较弱。对于有外中子源的次临界堆芯来说(第三种情况),即使在堆芯趋于临界的情况下,由于外中子源而引起的通量再分布效应仍然存在,而且可以预期,这种影响会随着堆芯次临界度的增加而迅速扩大。Y.A.Chao提出了利用空间修正的倒计数率(SCICR)测量控制棒价值的方法,即次临界棒价值测量方法(SRWM)。此方法最大的优势在于可以将控制棒价值的测量融合在提棒达临界或者落棒时间测量的过程中,从而缩减反应堆物理试验时间,提高经济效益。

图1 有/无、临界/次临界条件下堆芯功率分布比较[1］Fig.1 The comparison of power distributions with/without extraneous neutron sources,and under critical and subcritical core conditions.[1］

国内工作者在对国外次临界测量方法研究的基础上,针对每种方法自身的特点进行了修正和改善。其中张志伟[11］针对Greenspan的线性源增殖法(SME)在基态不处于临界状态时会放大测量误差的问题,提出了指数源增殖方法,该方法放宽了试验条件,缩短了测量过程,提高了精确性。彭飞[12］研究了脉冲中子源法中不同次临界度下的kβ/l与瞬发中子衰减常数α的实验关系式,并根据各次临界系统的栅格特性,通过改变次临界度得到关系式中的待定参数,从而得到系统在深次临界度情况下的α和反应性的新的关系式,并对新关系式的特性进行了较为详细的分析。次临界测量的空间修正方面的研究主要有:李文虎、朱国盛[13］等研究了利用断源法进行次临界负反应性测量的空间效应,指出在其所研究的试验条件下空间积分法是消除空间效应有效方法。尹延朋[14］等研究了带外源的次临界系统,阶跃加入反应性(落棒或者快速提棒)后,根据中子密度的解析解得到计算反应性的方程,并对CFBR-Ⅱ堆试验数据进行了验证核算。史永谦[15］等在试验装置上进行了源倍增方法的研究,深入研究了有外源情况下的有源次临界倍增因子Ks与外源位置和能谱的关系,同时测量了有效增殖因子Keff,验证了Ks与Keff之间的关系。郭海兵[16］等研究了临界外推过程中控制棒价值的线性修正,指出在考虑控制棒价值非线性影响的情况下,使用H-ΔK关系式做临界外推曲线具有很好的准确性和收敛速度,同时指出距离临界点较远(次临界度较深)的情况下,其准确性不如使用ΔKeff。蔡光明[17］等使用离线反应性仪计算程序(ODRM)计算了所有控制棒全插时的棒价值,开启了次临界测量在商用堆中的验证应用。总结国内工作者的研究,大多数集中在次临界测量方法的改进及修正方面并利用实验堆进行了改进验证。在此基础上,针对某一具体的次临界反应性测量方法,研究了其空间效应的影响并利用实验堆数据进行验证。

综上所述,国外科研工作者针对堆芯装料、达临界过程以及临界后的次临界反应性测量技术,特别是次临界反应性的空间修正技术进行了比较广泛而深入的研究,其中动态棒价值测量(DRWM)方法已经在商业堆中广泛应用,可以预期在未来的反应性测量方法中,这些技术有明显的应用需求和优势。国内这方面的工作基本上都在实验堆上进行,内容集中在某一方法的改进,缺乏系统性。对于商用堆国内还处于消化吸收国外技术的阶段,而且这种消化吸收更倾向于实践使用,忽略了原理及其方法的研究。本文针对这种需求,总结概括了国外商用堆次临界反应性测量的基本原理与方法,结合反应性测量仪表技术,给出了次临界反应性仪的数据处理流程,这对于推进国内商用堆次临界反应性测量的研究和实际应用具有较为重要的意义。

1 方法概括

1.1 点堆模型及其空间修正

实质上堆芯的反应性无法直接测量,反应性仪通过处理堆外探测器探测到的堆芯中子通量水平随时间的变化来推导出反应性,因此我们需要找到反应性与堆外探测器信号之间的关系。为了方便起见,先考虑没有外中子源的单群反应堆动力学方程:

上述公式具有普遍适用性。为了解上述方程,应用点堆模型假设:堆芯中子通量的分布(形状)不随时间发生变化,只是中子通量的大小(幅度)随时间变化。要满足上述条件,则要保证:

1)中子通量φ和缓发中子先驱核浓度Ci具有相同的基本形状;

2)φ和Ci在时空上是可以分离的;

3)堆芯处于稳定的状态,反应性阶跃变化引起的所有衰减项都已湮灭。

基于上述假设,在稳态下,瞬发中子和缓发中子在时间上按相同的规律发生变化,那么φ和Ci可以表示为:

将式(3)、式(4)代入式(2),可得:

再将(5)式代入式(1),则式(1)表示为:

考虑无外中子源的静态本征方程(静态情况下,缓发中子达到动态平衡,本征值是一个常数):

联合公式(6)和(7),可得:

公式(8)只有在以下两种情况下能够实现:

1)堆芯处于接近临界状态(即ω＝0, k＝1);

2)缓发中子参数βi和λi是空间常数并且或者Σf是空间常数或者1/νΣf可忽略不计。实际上反应堆的这种平衡状态几乎是不存在的,因为缓发中子参数随着不同裂变产物的分布在空间上有着明显的变化。

因此传统的反应性仪应用的情况是:

1)缓慢的插入或者提出控制棒,相应的进行稀释或者硼化操作,以补偿堆芯反应性;

2)在堆芯阶跃引入反应性后,等待足够长的时间,以使堆芯进入稳态;

3)通常测量小的反应性变化。

综上所述,点堆动力学方程应用于堆芯接近临界的稳态。对于次临界情况下的点堆动力学方程失效的情况,可以从临界外推倒计数率曲线(ICRR)的非线性得到验证。次临界情况下,控制棒动作、缓发中子效应、堆芯是否稳态以及外中子源都会产生空间效应而使点堆模型失效。

1.2 空间效应的程序修正流程

次临界反应性测量的关键问题是对次临界情况下堆芯表现出来的空间效应进行修正。近些年来,堆芯输运模拟软件的发展,使得我们可以模拟得到任何堆芯状态下探测器的信号,只要我们强迫模拟输入满足点堆模型的假设并结合对真实物理试验过程的模拟,我们就可以利用模拟的结果对实际测量过程的探测器信号及反应性仪的计算结果进行修正,这样便可以精确的求解次临界反应性。

研究发现[10］,在堆芯中子通量总量不变的情况下,控制棒在堆芯的位置不同,堆外探测器探测到的信号不相同。这是因为堆芯的不同配置(控制棒位置)导致了中子通量在堆芯分布的改变,这种由于堆芯配置变化而产生的空间效应,我们称之为静态空间效应。对静态空间效应修正的方法为:在堆芯中子通量总量不变的输入条件下,模拟控制棒插入导致的堆芯中子通量再分布效应,进一步得到堆外探测器信号,将此信号与标准状态(例如ARO)下的信号进行比较,便可以得出静态修正因子(SSF)。如果将不同堆芯配置下反应性仪测量的堆外探测器信号使用静态修正因子进行修正,那么就可以认为不同堆芯配置下中子通量的分布没有发生变化(按照标准状态分布,例如ARO),这样就将控制棒插入导致的静态空间效应转变为满足点堆动力学方程的稳态过程,便可以通过反应性仪求解出此稳态情况下的反应性。图2展示了使用MCNP和ANC软件进行上述修正的程序流程。利用MCNP软件可以求解堆芯各个节点对探测器的响应函数,响应函数可以采用伴随计算的方法求解[18-19］。根据探测器内10B与中子发生(n,α)反应的截面,确定探测器发射中子的能量。通过大量的随机数取样确定探测器发射的中子通量φw,经过伴随中子通量计算得出堆芯每一节块的中子通量φ＋,然后叠加节块的中子能谱χ,确定节块内每产生一个中子对探测器的贡献,即响应函数Ri。在总通量保持不变的前提下,利用ANC软件模拟不同堆芯配置下的中子通量分布,确定节块内的中子数。最终响应函数乘以节块中子数便可以求出探测器信号。

图2 点堆模型修正程序流程Fig.2 The process of point core spatial correction.

但是上述修正过程仍然建立在堆芯处于稳态且无外中子源的假设下,实质上任何反应性阶跃,堆芯都不可能瞬间处于稳态,而且如果要解堆芯某一刻的反应性,必须考虑缓发中子分布的滞后效应,另外中子源的存在也会对堆芯中子通量的分布产生影响。对于次临界反应性测量过程中的上述空间效应,使用ANC在“当前堆芯通量分布”输入条件下分别求解稳态/瞬态(DRWM)、无/有外中子源(SRWM)的反应性或探测器信号计算结果,这样可以使用软件模拟结果对反应性仪测量的经过静态空间修正的反应性进行二次修正,由此将软件模拟的结果与实际的物理测量过程联系起来,得到物理过程精确的堆芯反应性。

2 应用综述

2.1 低功率阶段次临界反应性测量—DRWM

2.1.1 DRWM空间修正的基本思路

在DRWM过程中控制棒不间断连续插入导致的静态空间效应可以通过静态修正因子(SSF)进行修正。反应性仪处理静态修正后的通量信号得到每一控制棒位置下稳态反应性。实际上在相应控制棒位置下,一方面堆芯不可能瞬间趋于稳态,另外缓发中子也不可能瞬间达到平衡分布状态,以上两种影响中子通量分布的效应被称为动态空间效应。在常规的控制棒价值测量过程中(例如:调硼法)其实也存在这种影响,但我们可以通过等待足够的时间让缓发中子达到平衡状态;而在DRWM过程中,由于没有反应性补偿,控制棒必须连续不间断插入,以减少测量时间,防止通量下降到次临界增殖的水平上。动态修正可以通过ANC对当前堆芯通量分布的瞬态反应性与稳态反应性进行比较实现,我们将这两者的比值称为动态修正因子(DSF)。通过静态修正因子及动态修正因子的修正,便可以准确求解出控制棒插入时刻的反应性。

2.1.2 ADRC处理DRWM数据流程图

图3 DRWM中ADRC处理数据原理及流程Fig.3 The theory and process of DRWM using ADRC

能够进行DRWM的反应性仪称为先进数字化反应性仪(ADRC),图3展示了在DRWM过程中ADRC处理数据的原理、数据流以及基本方程。ADRC接收到功率量程探测器的顶部和底部区域电流信号,电流信号经过SSF修正后进行归一化处理,SSF修正过程中的ERIF及EAIF参数来源于伴随中子通量计算(图3,公式1)。ADRC归一化后的信号使用SCM(Stiffness Confinement Method)方法[20］求解点堆动力学方程。ADRC通过求解方程得到的反应性为测量的稳态反应性,将此反应性经过DSF修正,便可以得到最终求解的反应性,DSF的计算公式见图3(公式3),其数值为ANC模拟的稳态反应性与瞬态反应性的比值。

2.1.3 DRWM测量实例

图4(A)显示了DRWM的测量过程,首先提出咬量棒,引入约60 pcm的反应性,使中子通量增长到设定值的75%水平,此时连续的以最大速度插入控制棒到堆芯底部,测量完成后全部提出控制棒,使通量再增长到设定值的75%水平,然后插入下一组测量棒,如此往复,直到所有控制棒价值测量完毕。

图4(B)显示了测量某一组控制棒时反应性和堆外探测器信号的变化曲线,虚线显示了控制棒开始插入和停止插入的时间,可见控制棒插入过程中,中子通量不断下降,在控制棒停止插入后,由于堆芯处于较深的负反应性,中子通量继续下降,在控制棒提升过程中,由于连续提棒的阶跃效应,使中子通量呈现上升趋势。

图4 DRWM测量过程以及SSF和DSFFig.4 DRWM measurement process,and SSF and DSF

图4 (C)显示了DRWM过程中静态修正因子随控制棒插入位置的变化曲线,从图中可以看出顶部探测器与底部探测器的修正因子的变化呈现相反的趋势,两者倾向于相互抵消。底部探测器的修正因子随着控制棒的插入呈现先增大后减少的趋势,而顶部探测器的修正因子则正好相反,两者在控制棒插入开始或结束时趋向于同一数值。这是可以理解的,在控制棒未插入时,由于此时堆芯仍然处于稳态过程,通量分布并没有改变,此时修正因子为1;随着控制棒的插入,堆芯上部的中子通量因控制棒的插入而降低,堆芯下部的中子通量相对增多(静态修正的思想是假设此时没有反应性引入,仅中子通量重新分布),这样需要对上部探测器信号引入一个小于1的修正量,对下部探测器信号引入一个大于1的修正量。当控制棒插入堆芯中平面以下,堆芯下部中子通量的减少量大于堆芯上部中子通量的减少量,因此需要给予与控制棒在中平面以上时相反的修正量。当控制棒完全插入时,此时堆芯上下部通量减少量趋于相同,因此上下部探测器的信号修正趋于相同的小于1的数值。

图4(D)显示了DRWM过程中动态修正因子随控制棒插入位置的变化曲线,从图中可以看出,在控制棒的两种插入速率下,动态修正因子的变化趋势基本重合,表明控制棒插入速度(0.83 s/步)还没有达到可以影响缓发中子导致功率重新分布的时间间隔。动态修正因子随控制棒的插入的变化并没有明显的趋势,但其数值波动较小(0.929～1.115);相比之下静态修正因子的波动范围较大(顶部:0.34～1,底部:0.73～1.8)。由于静态修正的过程是探测器的信号除以静态修正因子,归一化联合后再带入点堆动力学方程求解反应性,可以推导静态修正因子对反应性(控制棒价值)的影响与其本身的自然对数成正比,它与动态修正因子是两种不同的修正方式。

2.2 临界前次临界反应性测量—SRWM

2.2.1 SRWM空间修正的基本思路

临界前的次临界棒价值测量(SRWM)与低功率阶段的DRWM相比一个显著区别就是:临界前堆芯中子通量水平较低,外中子源对通量再分布的影响不可忽略。对于次临界的堆芯来说,中子通量是通过外中子源保持的。

SRWM的空间修正可以分成两部分: 1)由于堆芯配置(控制棒位置)的改变而引起的通量再分布效应,这一点与DRWM的空间修正相似,但不同的是动态空间效应可以通过等待足够长的时间来消除掉。因此SRWM堆芯配置的改变而引起的空间修正实质上是DRWM的静态修正过程。2)由于外中子源的存在对中子通量空间分布的影响,这种效应是SRWM所特有的,而且这种效应随着堆芯次临界度的增加而显著的增强,因此硼浓度的改变和控制棒的移动都会导致这种效应的变化。

外中子源的模拟可以使用ORIGEN2软件。ORIGEN2是一个广泛使用的放射性材料积累、衰变和处理的计算软件[18］。利用此软件得到堆芯外中子源随燃料组件燃耗变化的函数,然后作为输入导入ANC中进行计算(图5,公式3)。ANC模拟外中子源在堆芯中的分布并进行外中子源存在情况下的中子通量分布计算,进一步计算SRWM过程中堆外探测器信号的响应,提取次临界空间修正因子。

2.2.2 SDRC处理SRWM数据流程图

能够进行SRWM的反应性仪,我们称为次临界数字化反应性仪(SDRC),图5显示了临界前提棒过程中,测量堆芯次临界度的SDRC数据处理流程。对于堆芯配置空间修正,在“当前堆芯通量水平”下,进行无外中子源的通量分布计算,得到探测器信号以及堆芯次临界度,对探测器信号应用静态修正因子(图5,公式1)得到稳态探测器信号(满足点堆模型),进一步得到IRCC与Keff的线性关系(图5,公式2)。对于外中子源的空间修正,在“当前堆芯通量水平”下,进行有外中子源的通量分布计算,获取探测器信号并与稳态探测器信号进行比较,形成次临界空间修正因子。在实际测量过程中对测量信号应用次临界空间修正因子使堆芯行为满足点堆动力学方程,进一步利用倒计数率与次临界度的线性关系求出堆芯的次临界度。

2.2.3 SCICR测量实例

图5 SRWM次临界度测量的SDRC处理流程Fig.5 The the ory and process of SRW Musing SDRC

图6 MHI次临界空间修正验证(A,B)以及VogtleSRWM测量结果(C,D)[1］Fig.6 The verification of spatially corrected methodin MHI plant(A,B)and the SRWM measure mentresultin Vogtleplant(C,D)[1］

图6 (A,B)展示了对M HI电站二回路(A)和三回路(B)机组的提棒及稀释过程的ICRR数据进行空间修正的验证,堆芯初始状态是进行了很深的硼化并且有控制棒插入堆芯。从图中可以看出,经过空间修正后的ICRR曲线呈现较好的线性性质。当堆芯的次临界度很高时,未修正的ICRR曲线呈现较明显的非线性性质。在次临界度为0.06～0.03(A)和0.12～0.082(B)的区间内,未修正的ICRR接近于常数,这是由于此时的次临界度很深,二次中子源的输运对中子的再分布起决定作用,此时外中子源空间修正对探测器的信号起支配作用。可见当堆芯的次临界度很高时,在一定的次临界区间内,堆芯的动力学行为可以近似为外中子源在某静态次临界度下的输运行为,此时由于堆芯配置改变而引起的空间修正效应可以忽略。在次临界度为0.08～0的区间内,未修正的ICRR曲线趋近于线性分布,随着堆芯次临界度的降低,堆芯的动力学行为由两种空间效应共同决定。

图6(c,d)展示了在Vogtle Unit 2电站Cycle 12首次应用SRWM方法进行控制棒价值测量的结果,SRWM是结合在控制棒落棒时间测量规程中进行的。图c展示了控制棒从所有控制棒全插(ARI)到所有控制棒全提(ARO)过程中探测器N-31和N-31的测量的ICRR,图d展示了对测量的ICRR进行空间修正后,通过回归分析得到的修正ICRR与ANC预测的Keff之间的线性关系,测量结果发现测量的总控制棒价值与预测的总控制棒价值的偏差在1.9%,满足电站的验收准则,而且大约节约了12 h大修时间。

3 结论

本文在总结国内外次临界反应性测量技术发展的基础上,阐述了DRWM和SRWM技术中对点堆模型空间修正的基本思路,整理出空间修正程序模拟的流程及处理方法,并结合实例对堆芯配置和外中子源空间修正进行比较详细的说明,这对促进国内次临界反应性测量技术及先进反应性仪的研制和开发,对核电厂物理试验人员了解先进测量方法等方面具有一定的意义,通过本文还有助于了解未来物理试验可能的发展方向,促进国内相关技术的发展。

[1］ Y.A.Chao,H.Q.Lam.The Spatially Corrected Inverse Count Rate(SCICR)Method for Subcritical Reactivity Measurement.WCAP-16260-NP,2004.

[2］ M.Carta,G.R.Imel,C.C.Jammes,et al.Reactivity Measurements in Subcritical Core:RACE-T Experimental Activities,Fifth International Workshop on the Utilisation and Reliability of High Power Proton Accelerators SCK CEN Mol,Belgium,2007.

[3］ C.C.Jammes,R.I.George.Absolute Reactivity Calibration of Accelerator-Driven Systems after RACET Experiments,PHYSOR-2006,ANS Topical Meeting on Reactor Physics,2006.

[4］ G.R.Imel.Measurement of Sub-critical Reactivity, ADSS IAC,2005.

[5］ S.Dulla,P.Picca.Methods for the Reactivity Evaluation in Subcritical System Analysis:A Review International Conference on Mathematics and Computational Methods Applied to Nuclear Science and Engineering(M&C 2011),2011,ISBN 978-85-63688-00-2.

[6］ J.Wright,I.P'azsit.Neutron kinetics in subcritical cores with application to the source modulation method [J］．Annals of Nuclear Energy,2006,33(2):149-158．

[7］ R.Y.R.Kuramoto1,A.F.Miranda1,G.L. Valladares.Development of A Digital Reactivity Meter for Criticality Prediction and Control Rod Worth Evaluation in Pressurized Water Reactors,2009 International Nuclear Atlantic Conference-INAC 2009, ISBN:978-85-99141-03-8.

[8］ W.Naing,T.Masashi.Subcriticality Measurement of Pressurized Water Reactors by the Modified Neutron Source Multiplication Method[J］．Journal of Nuclear Science and Technology,2003,40(12),983-988．

[9］ Y.A.Chao.Dynamic rod worth measurement(DRWM) next generation rod worth measurement method,In: Proc.Topl.Mtg.Advances in Nuclear Fuel Management II,Myrtle Beach,South Carolina,USA, 1997 March 23-26:14-53,ANS

[10］ Y.A.Chao,D.M.Chapman,D.J.Hill,L.R. Grobmyer.Dynamic rod worth measurement[J］．Nuclear Technology,2000,132:403-412．

[11］ 张志伟.测量次临界度和负反应性的指数源增殖方法(SME)[J］．核科学与工程,1984,4(2):167-170．

[12］ 彭飞.次临界系统中瞬发中子衰减常数与反应性的关系[J］．核科学与工程,1987,7(2):180-183．

[13］ 李文虎,朱国盛.负反应性测量中断源法空间效应的研究[J］．原子能科学技术,1984,5:598-601．

[14］ 尹延朋,龚建,郑春.反应性阶跃方法测量次临界反应性[J］．原子能科学技术,2012,46:366-369．

[15］ 史永谦,朱庆福,夏普,李义国.反应堆物理实验中的源倍增法研究[J］．核科学与工程,2005,25(1):14-19．

[16］ 郭海兵,李润东,牛伟力.临界外推中对控制棒价值非线性的修正[J］．原子能科学技术,2013,47(1): 101-104．

[17］ 蔡光明.落棒法测量所有控制棒全插时的棒价值[J］．核科学与工程,2005,25(2):154-158．

[18］ I.X.Rootsaert.Source Range Neutron Detector Response in Core Reloading Conditions,BNEN/ 2006-2007．

[19］ S.L.Anderson,K.R.Balkey.Adjoint Flux Program for Point Beach Units 1 and 2,WCAP-10638,1984.

[20］ Y.A.Chao.A.Attard,A Resolution of the Stiffness Problem of Reactor Kinetics[J］．NuclearScience& Engineering,1985,90(1):40-46．

The Spatial Effect in Subcritical Reactivity Measurement and Its Applications

ZHANG Yu,WANG Shi-wei,XU Lin-lin,ZHENG Zheng

(Engineering Department,Shangdong Nuclear Power Company,Haiyang,Shandong Prov.265116,China)

The technology of reactivity measurement in subcritical core has special feature since the distribution of core flux will be affected by the movement of the control rod,the core subcritical reactivity and the extraneous neutron sources.Therefore,in general,the subcritical core reactivity can only be“monitored”but never accurately“measured”.With the development of core simulation codes,the foreign scientists have researched the spatially corrected methods applied to the point core in subcritical reactivity measurements.Based on these methods,the dynamic rod worth measurement (DRWM)and sub-critical rod worth measurement(SRWM)were applied in plant control rod worth measurements and some nuclear plants in China have used these methods but few people in China introduced the theory of spatially corrected methods.Based on this condition,this article summarized the theory and methods in subcritical core reactivity measurements and provided the data treating process using the reactivity measurement instruments.This will be useful to the research and application of subcritical reactivity measurement in commercial nuclear power plants in China.

subcritical;spatial effect;point core;control rod worth

O571.53

A

0258-0918(2016)01-0016-11

2014-11-28

中国电力投资集团科技项目(2013-049-ZHD-KJ-X)

张 瑜(1985—),男,山东烟台人,工程师,硕士,现从事反应堆物理试验、堆芯管理工作