第七章 反应性控制

1、第七章：反应堆控制核反应堆工程原理反应性随时间的变化反应性随时间的变化反应性随时间的变化本章主要问题包括：核燃料同位素成分的变化和燃耗；裂变产物同位素的生成与消耗；反应堆启动和停堆后135Xe和149Sm中毒随时间的变化；反应性随时间的变化；堆芯寿期、燃耗深度以及核燃料的转换与循环等问题。核燃料中重同位素成分随时间的变化7.1.1 核燃料中重同位素的燃耗方程核燃料中重同位素的燃耗方程p在反应堆的运行过程中，核燃料中的易裂变同位素(铀和钚)不断地燃耗根据粗赂的估计，一个电功率为1000MW的核电厂每天大约要消耗3kg左右的235U(或239Pu)；p另一方面，可转换材料(如238U或232Th)

2、俘获中子后又可以转换成易裂变同位索(如239Pu或233U)；p同时，由于裂变将产生300多种裂变产物；p因此，核燃料中各种重同位素的核密度将随反应堆运行时间不断地变化。这种变化与所采用的燃料循环的类型有关。核燃料中重同位素成分随时间的变化PuCmCmAmAmAmPuPuPuPuPuNpUPuNpUU238242%75244243242%25241244242241240239239238238237236235铀钚燃料循环中重同位素燃耗链核燃料中重同位素成分随时间的变化应该指出，上图中所列出的燃耗链是经过了简化的，只保留了工程计算中有重要意义的一些核素，略去了一些半衰期较短或者吸收截面较小的

3、中间元素的作用。上图中从238U到240Pu的实际链如下图所示。Np：镎Am：镅Cm：锔图7-1核燃料中重同位素成分随时间的变化可以看出，这里我们忽略了239U，240U和240Np的作用。因为239U的吸收截面很小而半衰期又很短，所以绝大部分的239U都通过-衰变成了239Np，而生成的240U是极少数的。同理，对其他的一些核素链也作了类似的简化，这样就简化成了7-1图中的链。核燃料中重同位素成分随时间的变化Pa：镤图7-2核燃料中重同位素成分随时间的变化对于裂变产物连，情况更为复杂。把由裂变反应直接产生的裂变碎片以及随后由这些碎片经过放射性衰变形成的各种同位素统称为。它大约包括300多种放

4、射性及稳定的同位素。因此，要分别计算它们的浓度变化及其对反应性的影响是非常复杂和耗时的工作。在工程中，计算时一般只需要选取其中吸收截面大或裂变产额较大的一些主要同位素，如135Xe，149Sm等(吸收截面都大于104b)，单独进行计算。核燃料中重同位素成分随时间的变化对于其他的裂变产物，按其截面的大小及浓度随时间变化特性归并成两组“假想的集总裂变产物”(FP)：一组是吸收截面相对大一点，其浓度随运行时间的增加而缓慢的趋于饱和的，称之为慢饱和裂变产物(SSFP)；另一组是截面很小的非饱和裂变产物(NSFP)。核燃料中重同位素成分随时间的变化对于这两组假想的裂变产物的产额及截面作如下处理：令其裂变

5、产额=i，而假想吸收截面则可根据经验数据或用该组裂变产物的吸收截面对其裂变产额进行加权平均而近似得到。例如，对于SSFP，。组的所有裂变产物求和是对并入截面和裂变产额。求和素的热中子吸收种慢饱和裂变产物同位分别是第和其中，SSFPiiiaiiiiaiSSFPa核燃料中重同位素成分随时间的变化7.1.2 核燃料中重同位素的燃耗方程核燃料中重同位素的燃耗方程若要准确的计算堆运行过程中燃料内成分的变化，首先必须建立这些同位素的燃耗方程，在图7-1和7-2的燃耗链和裂变产物链中，除241Am链和Pm-Sm-Eu链外，其他大部分都是按单独裂变产物处理。 241Am链或Pm-Sm-Eu链，计算时可以予以线

6、性化处理以简化计算。如下图所示，同位素E的产生可转化成右边两个路径获得，计算时分别当作两个独立的链和核素处理，计算后再把结果相加。即EDBAECBAEDCBA转化为核燃料中重同位素成分随时间的变化经过线性化后，可以对7-1和7-3图中的每个核写出其浓度的燃耗方程。为此对图中的核素依序给与一个编号，如下表所示。编号12345678核素235U236U237Np238Pu238U238Np239Pu240Pu核燃料中重同位素成分随时间的变化这样，对其中每个核素，其燃耗方程如下:核素的产额。易裂变核裂变时对是是衰变常数，其中：或其中：iitrNtrFtrFtrNtrtrNdttrdNiiGgiigi

7、gfiiiGggigiiiiGggigaiiii,1,11,111,11,表示由于同位素i-1的吸收中子或由于衰变而导致同位素i的产生率表示由于同位素吸收中子和衰变而引起的总消失率表示由于裂变反应引起的产生率核燃料中重同位素成分随时间的变化是一个非线性问题。而严格来讲，方程间分布，因成分的核子密度及其空度的分布又取决于燃料通量密有关，而且反过来中子与中子通量密度位素成分杂的是，不仅燃料的同算比较困难。同时更复方程，计是一个变系数的偏微分的函数，因此方程和时间间量密度和核密度都是空由于在反应堆内中子通27,27trtrNtri核燃料中重同位素成分随时间的变化在实际计算中，通常采用一些近似方法来解

8、决上述问题。首先，把堆芯划分成若干个子区，称为燃耗区。在每个燃耗区内，中子通量密度和核子密度随空间位置变化不大，可以认为等于常数，或可以用它们在该区的平均值近似代替。这样，在给定的燃耗区内中子通量密度和核子密度就不再是空间 的函数；其次，把运行时间t也分成许多的时间间隔，每一时间间隔(tn-1,tn)称为燃耗时间步长。由于运行的反应堆内堆芯成分的变化并不很快，中子通量密度的空间分布形状随时间的变化很缓慢，所以时间步长可以取到长达几个星期或更长。在每个时间步长中，可以近似认为中子通量密度不随时间变化而等于常数。这样，就在每个燃耗步长每消去了中子通量密度函数对自变量t的依赖关系。r核燃料中重同位素

9、成分随时间的变化做了上述假设后，对于给定的燃耗区，在给定的燃耗步长内，燃耗方程便可以简化为常系数的常微分方程组 1,111,11,iiiGggigaiaiiaiiiiiiiIIIFNtNdttdN或其中：核燃料中重同位素成分随时间的变化同时，方程(7-2)是一组耦合方程组，为了解析求解的方便，可以进一步加以线性化。为此，令 为燃耗步长。，来代替子数即用该步长内的平均核ttNNdttNtNiitttiinn1核燃料中重同位素成分随时间的变化经过这样的线性化后，方程(7-2)和(7-5)的产生项Fi便等于常数。在燃耗步长t内，i，i和Fi都与时间无关。而且对于非裂变产物(图7-1)中的燃耗链核素F

10、i=0。核燃料中重同位素成分随时间的变化7.1.3 燃耗方程的求解燃耗方程的求解对于上述线性化近似下，方程(7-5)为一常系数的一阶微分方程，可以用解析方法或常规数值方法求解。核燃料中重同位素成分随时间的变化1解析方法设在燃耗时间步长t=tn+1-tn内，令=t-tn，tn为该燃耗步长的起始时间。对常微分方程组，自i=1方程可以直接写出其解的形式为C11是由初始条件确定的常数。 1111111;1FDDeCN核燃料中重同位素成分随时间的变化对于i=2的方程，N2()的解由方程的齐次部分的通解和特解两部分组成，C21，C22和D2为待定常数。 22122211DeCeCN核燃料中重同位素成分随时

11、间的变化这样一直顺推下去，可以写出(7-5)的通解，其形式如下：为确定系数Cij，可将通解带入方程(7-5)，得到 iijijiDeCNj1iiijijiiijjiiijjijFDeCDeCeCjjj1111, 111核燃料中重同位素成分随时间的变化这样，根据方程两端对应阶次项的系数相等的原则，可得FDDijCCFDDCCCjiiiijijiiijijiiiiijijiijij, 11, 11, 11, 111;,0;因而：核燃料中重同位素成分随时间的变化对于i=j，(7-11)为恒等式，因此Cii必须用初始条件来确定。令t=tn(即=0)时刻同位素i的和密度为Ni(0)，从(7-10)得 i

12、ijijiiiiijijiDCNCFNCDCN111111110;00因而：核燃料中重同位素成分随时间的变化这样就给出了线性化后燃耗方程(7-5)在00：超临界； 1，反应堆运行时调节keff使其为1，停堆时调节keff使其小于1。反应堆冷态停堆情况下(假使全部停堆系统全部移出堆芯)反应性大于0的部分称作剩余反应性。 核燃料中重同位素成分随时间的变化7.1.2 裂变产物裂变产物135Xe和和149Sm的中毒的中毒在下面我们从简单的单群四因子模型来讨论一下裂变产物对反应性的影响，首先写出在单群近似下的有效增殖因数表达式其中F，M表示燃料和慢化剂。pkMaFaFf核燃料中重同位素成分随时间的变化假

13、设裂变产物对逃脱共振俘获概率p和不泄漏概率不产生重大的影响，那么有裂变产物积累时的有效增值因数将等于观吸收截面。是裂变产物的热中子宏其中：PaPaMaFaFfpk核燃料中重同位素成分随时间的变化这样，根据反应性的定义，可以导出裂变产物所引起的反应性的变化为a为不存在裂变产物时芯部的热中子宏观吸收截面。这种由于裂变产物的存在，吸收中子而引起的反应性变化称为裂变产物中毒中毒。aPaMaFaPakkk应该指出，由于裂变产物的分布是非均匀的，并且采用均匀裸堆的单群四因子模型也是很粗略的，因此采用上式来计算裂变产物中毒的误差是比较大的。在反应堆的实际设计中，一般都是采用数值方法直接对裂变产物中毒进行计算

14、。核燃料中重同位素成分随时间的变化在裂变产物中，对热中子反应堆来说，有两种同位素特别重要：和。一方面是因为它们具有非常大的热中子吸收截面和裂变产额，因而其浓度在反应堆启动后便迅速增长，不久便趋于饱和，对反应性有效大的影响；另一方面，由于放射性的衰变使它们的浓度在工况变化时发生迅速的变化。 这些将使在反应堆的启动、停堆及功率升降时反应性在较短时间内发生较大的变化，给运行造成困难。因此，对于热中子反应堆来讲，有必要对这两个同位素因此，对于热中子反应堆来讲，有必要对这两个同位素的中毒情况进行单独的研究和计算。的中毒情况进行单独的研究和计算。核燃料中重同位素成分随时间的变化7.2.1 135Xe中毒中

15、毒在热中子反应堆中， 135Xe是所有裂变产物中最重要的一种同位素，一方面因为它的热中子吸收截面非常大，如图所示，在中子能量为0.025eV时， 135Xe微观截面达2.7106b左右。在热能范围内它的平均吸收截面大约为3106b，因此在热中子反应堆中，必须认真考虑135Xe中毒所带来的影响。但是在高能区，135Xe的吸收截面随中子能量的增加而显著的下降，因此在快中子反应堆中，氙中毒的影响是比较小的。核燃料中重同位素成分随时间的变化另一方面，135Xe的裂变产额比较大，虽然235U核裂变时，135Xe的直接产额仅为0.00228，但是它的先驱核的直接裂变产额去很高，它们经过-衰变后形成135X

16、e，这样135Xe的中体产额，可达到6%以上。下图给出了质量数为135的裂变产物的衰变链。核燃料中重同位素成分随时间的变化由图可以看出，135Sb(锑)和135Te(碲)的半衰期都非常短，因此可以忽略它们在中间过程中的作用，把135Sb和135Te的裂变产额和135I的直接裂变产额之和作为135I的裂变产额，即其中， 是135I的直接裂变产额。由于135I的热中子截面仅为8b，它的半衰期也只有6.7h，在热中子通量密度为1014cm-2s-1的时候， ，即135I由吸收中子引起的损失项远小于它衰变引起的损失项，因此可以忽略135I对热中子的吸收，认为135I全部都衰变成135Xe，这样就可以简

17、化上图，得到简化后的135Xe衰变图。ITeSbII410IIa核燃料中重同位素成分随时间的变化核燃料中重同位素成分随时间的变化以单群为例，根据上图，可以写出135I和135Xe的浓度随时间变化的方程式 给出。的数值由表和，其中，27XeIXeIXeXeaXeIIfXeXeIIfIItNtNdttdNtNdttdN核燃料中重同位素成分随时间的变化由于135Xe具有很大的吸收截面和短的半衰期，因而在反应堆启动后，135Xe的浓度将很快的增加并趋近饱和，而停堆后又将很快的衰变，这些变化将使反应性在短时间内发生较大的变化，给反应堆的运行带来很多问题。下面将讨论反应堆在启动、停堆以及功率变化时的氙中毒

18、情况。核燃料中重同位素成分随时间的变化1反应度起动时反应度起动时135Xe中毒中毒对于一个新的反应堆，135I和135Xe的初始浓度都是0。若反应堆在t=0时刻开始启动，并且很快就达到了满功率，那么就可以近似的认为在t=0时刻中子通量密度瞬时达到了额定值，并且一直保持不变。利用(7-24)和(7-25)式，并采用下列的初始条件： 000XeINN核燃料中重同位素成分随时间的变化核燃料中重同位素成分随时间的变化核燃料中重同位素成分随时间的变化核燃料中重同位素成分随时间的变化核燃料中重同位素成分随时间的变化核燃料中重同位素成分随时间的变化动力反应堆在额定功率运行时的热中子通量密度一般都比较高，因此

19、可采用(7-32)式来近似地计算平衡氙中毒。但若反应堆在低于额定功率下运行，则平衡氙中毒就与运行的功率大小有关，图7-10表示在某反应堆运行过程中，平衡氙中毒与反应堆中子通量密度水平或运行功率的关系。核燃料中重同位素成分随时间的变化 2停堆后停堆后135Xe中毒中毒从图7-8可知，在反应堆运行时，135Xe的产生有两条途径，即由燃料核裂变直接产生和由135I的-衰变而产生，前者与反应堆的中子通量密度值有关。由于135Xe的裂变产额比较小，而且只要反应维运行两天以后， 135I已达到饱和浓度，这时135Xe主要是由135I的-衰变产生的。 135Xe的消失也有两条途径，即由于135Xe吸收中子和

20、135Xe的-衰变而消失，前者也与反应堆的中子通量密度有关。当反应堆中平均热中子通量密度值为0.7561013cm-2s-1时，由135Xe吸收中子和由135Xe的-衰变所引起的消亡率两者刚好相等。但在热中子反应堆中，平均热中子通量密度一般都大于这个值，因此在正常功率运行时，135Xe主要靠吸收中子而消失。核燃料中重同位素成分随时间的变化核燃料中重同位素成分随时间的变化核燃料中重同位素成分随时间的变化核燃料中重同位素成分随时间的变化核燃料中重同位素成分随时间的变化核燃料中重同位素成分随时间的变化图7-11是停堆前后135Xe浓度和剩余反应性随时间变化的示意图。从图中可知，停堆后135Xe的浓度

21、先是增加到最大值，然后逐渐地减小；剽余反应性随时间变化则与135Xe浓度的变化刚好相反，先是减小到最小值，然后又逐渐地增大，通常把这一现象称为“碘坑碘坑”，因为这一现象主要是由于停堆后135I继续衰变成135Xe，使135Xe浓度增大所引起的。从停堆时刻开始直到剩余反应性又回升到停堆时刻的值时所经历的时间称为碘坑时间，以tI表示。在碘坑时间内，若剩余反应性还大于零，则反应堆还能靠移动控制棒来启功，这段时间称为允许停堆时间，以tp表示；若剩余反应性小于或等于零，则反应堆无法启动，这段时间称为强迫停堆时间，以tf表示。核燃料中重同位素成分随时间的变化停堆后反应堆剩余反应性下降到最小值的程度称为碘坑

22、深度。碘坑深度与反应堆停堆前运行的热中子通量密度(或运行功率)值密切相关，热中子通量密度愈大，碘坑深度愈深。图7-12给出了在不同热中子通量密度水平下运行的反应堆，在停堆后氙中毒随时间的变化曲线。从图中可看出，若热中子通量密度水平低于11013cm2s-1，则停堆后撤中毒变化很微小；但若热中子通量密度大于11013cm2s-1，则停堆后氙中毒变化很显著(也即碘坑深度很深)。假使停堆前反应堆的剩余反应性不足以补偿其氙中毒，那就会出现强迫停堆现象。核燃料中重同位素成分随时间的变化核燃料中重同位素成分随时间的变化 3功率过渡时的功率过渡时的135Xe中毒中毒假设反应堆在稳定功率下运行了一段时间，已在

23、堆内建立了平衡氙浓度，而在t=0时刻突然改变它的功率，则相应热中子通量密度要从l变成2，堆心内的135I和135Xe的浓度也要发生改变。这里，在解方程组(7-24)和(7-25)时，所采用的初始条件为核燃料中重同位素成分随时间的变化核燃料中重同位素成分随时间的变化核燃料中重同位素成分随时间的变化核燃料中重同位素成分随时间的变化核燃料中重同位素成分随时间的变化核燃料中重同位素成分随时间的变化但是必须注意到，以上的过程并不会单向地无限制地发展下去。有两个因素限制着它的变化：p在中子通量密度分布进一步倾斜后，形成了中子通量密度的梯度，中于通量密度高的区域向中子通量密度低的区域有一个净的中子流，这使中

24、子通量密度趋向平坦一些；核燃料中重同位素成分随时间的变化p在中子通量密度下降的区域内，135I的产生量也会相应地减少，因而由它定变成135Xe的量也减少，这就使135Xe的浓度由原来的增加逐渐转为减小，相应地，该区的中子吸收由原来的增加逐渐转为减少，从而使该区的中子通量密度和功率密度由原来下降转为上升。而在中子通虽密度上升的区域内的情况与上述的刚好相反，该区135Xe的浓度由原来减小转为增加，中子通量密度由原来上升转为下降，见图7-14(b)。这样，中子通量密度(或功率密度)变化将沿着原来相反方向进行，并重复地循环下去。这就形成了功率密度、中子通量密度和135Xe浓度的空间振荡，简称氙振荡氙振

25、荡。这种振荡可能是稳定的，也可能是不稳定的，这取决于反应堆的中于通量密度水平和它的物理特性。图7-15形象地给出了反应堆内由于氙振荡引起的功率振荡图，氙振荡的周期大约是15-30 h。核燃料中重同位素成分随时间的变化核燃料中重同位素成分随时间的变化核燃料中重同位素成分随时间的变化核燃料中重同位素成分随时间的变化核燃料中重同位素成分随时间的变化核燃料中重同位素成分随时间的变化核燃料中重同位素成分随时间的变化核燃料中重同位素成分随时间的变化核燃料中重同位素成分随时间的变化核燃料中重同位素成分随时间的变化反应性随时间的变化与燃耗深度反应性随时间的变化与燃耗深度反应性随时间的变化与燃耗深度反应性随时间

26、的变化与燃耗深度反应性随时间的变化与燃耗深度反应性随时间的变化与燃耗深度在图7-19中给出某压水堆的考虑平衡氙中毒和最大氙中毒后的有效增殖系数随燃耗深度变化的计算曲线。注意，这里的运行时间是指有效满功率天(EFPD)，它表示反应堆在满功率下运行的天数。考虑到在堆芯寿期末，反应堆运行时控制调节需要一定的反应性，因此堆芯寿期未的有效增殖系数应稍大于1，例如图中取k=1.005，由此便可确定出堆芯的寿期。 TCdtPtPEFPDT00反应性随时间的变化与燃耗深度反应性随时间的变化与燃耗深度反应性随时间的变化与燃耗深度反应性随时间的变化与燃耗深度反应性随时间的变化与燃耗深度反应性随时间的变化与燃耗深度

27、反应性随时间的变化与燃耗深度反应性随时间的变化与燃耗深度核燃料的转换与循环核燃料的转换与循环核燃料的转换与循环核燃料的转换与循环核燃料的转换与循环 从上式中可知，转换比是空间和时间的函数。但是我们感兴趣的主要是堆芯或某个区域的平均转换比，因此对整个维芯或区域进行体积积分，有核燃料的转换与循环核燃料的转换与循环进一步考察转换(或增殖)过程的物理特征和实现增殖的条件。易裂变核每吸收一个中子的中子产额为，显然，为维持链式反应必须的一个中子被易裂变核吸收，余下的再扣除被其它材料所吸收和泄漏损失后剩余的中子，加上可转换材料的快裂变中子数才是被可转换材料吸收而用于转换过程。因此，根据中子平衡和CR的定义有

28、核燃料的转换与循环核燃料的转换与循环核燃料的转换与循环讨论影响增殖的一些因素。首先，核燃料和可转换材料以外的其它物质(冷却剂、结构材料和裂变产物等)的吸收(A)将使转换比(或增殖比)降低。由于在热能区城内这些材料的中子吸收截面(尤其是裂变产物识)比较大，因而它们对热中子反应堆的影响最为严重，并使转换比显著降低。对于快中子反应堆，结构材料以及裂变产物的有害吸收就比较小，影响增殖较小。这从表7-4中A的数值也可以看出。泄漏损失同样使转换比(或增殖比)降低。对于热中子反应堆来说，这个数值并不太大。但是对于快中子反应堆，由于芯部体积较小，L的数值相当可观，因此为了减少中子的泄漏损失，通常在快中子反应堆

29、的芯部外面加一层由可转换材料构成的“再生区”，用来吸收泄漏出堆芯的中子，以提高增殖比。再生区通常很厚，因而采用了再生区以后，快堆的泄漏损失便可减小到极小的程度。238U和232Th的快裂变份额F对于热中子反应堆来说是很小的，而对于快中子堆则可以达到0.20左右。表7-4中列出了典型的热中子堆和快中子堆的一些数据。综上所述，可以看到快堆具有作为增殖堆的许多有利的条件。核燃料的转换与循环核燃料的转换与循环核燃料的转换与循环核燃料的转换与循环核燃料的转换与循环最后应该指出，上述的倍增时间，并没有考虑到增殖的裂变材料在反应堆内的停留时间、燃料元件后处理和元件制造所需的时间以及燃料运输所需的时间。实际的

30、倍增时间显然要比(7-67)式算出的要大得多。核燃料的转换与循环倍增时间是增殖堆增殖特性的一个重要指标。从倍增时间是增殖堆增殖特性的一个重要指标。从(7-72)式可以看出，如果要缩短倍增时间，则必须提高增殖比式可以看出，如果要缩短倍增时间，则必须提高增殖比和反应堆的单位质量易裂变材料的功率和反应堆的单位质量易裂变材料的功率(重量比功率重量比功率)输输出。对于动力快堆，倍增时间必须结合反应维的经济分出。对于动力快堆，倍增时间必须结合反应维的经济分析，综合地加以考虑。析，综合地加以考虑。温度效应与反应性控制温度效应与反应性控制温度效应与反应性控制反应性系数反应堆的反应性相对于反应堆的某一个参数的变

31、化率称为该参数的反应性系数。如：反应性相对于温度的变化率称为反应相的温度系数；相对于功率的变化率称为功率系数等。参数变化引起的反应性的变化将造成反应堆中子密度或功率变化，该变化又会引起参数的进一步变化，这样就造成了一种反馈效应。反应性系数的大小决定了反馈的强弱。反应性温度系数及其对核反应堆稳定性的影响反应性温度系数及其对核反应堆稳定性的影响反应性温度系数及其对核反应堆稳定性的影响反应性温度系数及其对核反应堆稳定性的影响反应性温度系数及其对核反应堆稳定性的影响反应性温度系数及其对核反应堆稳定性的影响为了进一步说明温度系数对反应堆稳定性的影响，图8-1表示在不同温度系数的情况下。反应性温度系数及其

32、对核反应堆稳定性的影响当反应堆内引入一个阶跃正反应性之后，反应维的功率随时间变化情况可以从图上看出，在温度系数大干零的情况下，反应堆的功率将很快地升高。当温度系数小于零且它的绝对值很小时，同时热量导出又足够快的情况下，反应堆的功率在开始时也较快地上升。但功率上升使反应堆的温度逐渐地升高，反应堆的反应性逐渐地减小。当反应堆的功率上升到某一水平，温度效应所引起的负反应性刚好等于引入的正反应性时，反应堆就在这一功率水平下稳定运行。在温度系数小于零且它的绝对值又很大，同时热量的导出不够快的情况下，反应堆的功率开始时也较快地上升。由于导热不快，所以反应推的温度增加很快，反应推的正反应性很快地就下降到零以

33、下。这时，反应堆就处于次临界状态，反应堆的功率开始下降，温度也随之下降。温度下降所引起的正反应性使反应堆的反应性开始上升。当功率下降到某一值时，反应堆的反应性刚好为零，这时，反应推就在这一功率下稳定地运行。由此可见，负温度系数对反应维的调节和运行安全由此可见，负温度系数对反应维的调节和运行安全都具有重要的意义。都具有重要的意义。燃料温度系数燃料温度系数燃料温度系数当反应堆的功率发生变化时，燃料温度立即发生变化，而慢化剂温度还来不及发生变化。这时在(8-6)式中只有I随燃料温度而变化。把(8-6)式代人(8-5)式使得到当燃料温度升高时，有效共振积分增加，即dI/dTF0。所以在以低富集铀为燃料

34、的反应堆中，燃料温度系数总是负的。图8-2给出某压水堆燃料温度系数与燃料温度的关系。燃料温度系数此外，燃料温度系数与燃料燃耗也有关系。在低富集铀为燃料的反应堆中，随着反应堆的运行，239Pu和240Pu不断地积累。240Pu对于能量靠近热能的中子有很强的共振吸收峰，它的多普勒效应使燃料负温度系数的绝对值增大。在核反应堆物理设计时，通常必须计算堆芯运行初期和运行末期在不同功率负荷情况下的燃料温度系数。慢化剂温度系数由单位慢化剂温度变化所引起的反应性变化称为慢化剂温度系数。由于热量是在燃料棒内产生，热量从燃料棒通过包壳传递到慢化剂需要一段时间，因而慢化剂的温度变化要比燃料的温度变化滞后一段时间。所

35、以，慢化剂温度系数滞后于功率的变化，故慢化剂温度系数属于缓发温度系数。慢化剂温度系数慢化剂温度系数慢化剂温度系数慢化剂负温度系数有利于反应堆功率的自动调节。例如在压水动力堆中，当外界负荷减小时，汽轮机的控制阀就自动关小一些，这就使进入堆芯的水温度升高。当慢化剂温度系数为负值时，反应堆的反应性减小，功率也随之降低，反应堆在较低功率的情况下又达到平衡。同理，当外界负荷增加时，汽轮机的控制阀自动开大一些，这就使进入堆芯的水温下降，反应堆的反应性增大，功率随之升高，反应堆在较高的功率下又达到平衡。其他反应性系数其他反应性系数其他反应性系数其他反应性系数其他反应性系数从核电厂运行的角度看，更有意义的是功

36、率系数的积分效应，即功率亏损，这里的亏损指的是反应堆功率升高时，向堆芯引入了负的反应性，是反应性亏损了。功率亏损PD指从零功率变化到满功率时反应性的变化P0表示满功率。从上式可以看出，反应堆功率升高时，d/dP一般为负。由于功率亏损一定得向堆芯引入一定量的正反应性来补偿功率亏损引入的负反应性，才能维持反应堆在新的功率下稳定运行，这是非常重要的一点。00PPDdPdPd温度系数的计算上面定性地分析了影响反应堆温度系数的各种因素，但温度系数的具体计算是比较复杂的。实际上需要对反应堆作不同温度T下的临界计算。计算时，首先计算在不同的燃料或慢化剂温度条件下堆芯的群常数，然后利用堆芯扩散计算程序，对反应

37、推进行临界计算，直接计算出在不同的燃料或慢化剂温度下的有效增殖系数keff(T)，求出k和T的比值，从而求得温度系数。温度系数的计算 以这种方法所计算出的结果是指在所计算的温差范围内的平均温度系数。计算的准确度与所取的温差大小有关。一般说来，所取的温差越小，计算所得的温度系数精度越高。但是，当温差很小时，有可能在这个温差下计算的k与k本身的计算误差相当，反而影响计算的准确度。虽然采用提高临界计算精度的方法能改善k计算的准确度，但由于在有效增谊系数的计算中存在着固有的计算误差，这样就限制了温度系数计算的准确度。在这种情况下，采用微扰理论方法来计算温度系数是比较合适的。反应性控制的任务和方式反应性

38、控制的任务和方式反应性控制的任务和方式反应性控制的任务反应性控制的任务反应性控制的方式反应性控制的方式根据上述控制方法，目前反应堆采用的反应性控制方式主要有如下三种： 控制棒控制； 固体可燃毒物控制，主要用于补偿部分初始过剩反应性； 化学补偿控制，主要在冷却剂中加入可溶性硼酸溶液来补偿过剩反应性。反应性控制的方式控制棒控制控制棒的一般作用和一般考虑控制棒的一般作用和一般考虑控制棒的一般作用和一般考虑控制棒的一般作用和一般考虑控制棒价值的计算控制棒价值的计算控制棒插入深度对控制棒价值的影响插入堆芯不同深度的控制棒价值通常用控制棒的积分价值和微分价值来表示p控制棒的积分价值当控制棒从一初始参考位置

39、插入到某一高度时，所引入的反应性称为这个高度上控制棒积分价值。参考位置选择堆芯顶部，则插棒向堆芯引入负反应性。随着插入深度越大，所引入的负反应性也越大。积分价值在帮为处于顶部时等于零。图8-6给出了典型的控制棒积分价值曲线，图中PCM为习惯上采用的反应性单位，1PCM=10-5；控制棒插入深度对控制棒价值的影响p控制棒的微分价值在反应堆设计和运行时，不仅需要知道控制棒完全插入时的价值，而且还需要知道控制棒在堆芯不同高度处移动单位距离所引起的反应性变化，即控制棒的微分价值，其单位常采用PCM/cm。它的表现形式如下为反应性的变化，H为棒位变化量。控制棒的微分价值是随控制棒在堆芯内的移动位置而变化

40、的。图8-7给出了典型的PWR反应堆中控制棒组的微分价值与其高度的关系。控制棒组表示一起移动的一组控制棒。棒微分价值是积分价值曲线上相应点的切线斜率。Hdzdc控制棒插入深度对控制棒价值的影响 从图8-6和图8-7可知，当控制棒位于靠近堆芯顶部和底部时，控制棒的微分价值很小并且与控制棒的移动距离呈非线性关系；当控制棒插入到中间一段区间时，控制棒的微分价值比较大并且与控制棒的移动距离基本上呈线性关系。根据这一原理，反应堆中调节棒的调节带一般都选择在堆芯的轴向中间区段。这样，调节棒移动时所引起的价值与它的插入深度呈线性关系。控制棒间的干涉效应控制棒间的干涉效应为了定性地说明相互干涉效应，我们考虑堆

41、芯中只有两根控制棒的情况。见图8-8所示。堆芯中没有控制棒插入时，径向中子通量密度分布如图8-8中虚线所示。当第一根控制棒完全插入堆芯时，径向中子通量密度分布如图中实线所示。控制棒的价值是与其所在处中子通量密度平方成正比。假如把第二根控制棒插在第一根控制棒附近的d1处，由于该处的中子通量密度比原来无控制棒时的中子通量密度下降了，因此第二根控制棒的价值比它单独插入堆芯时的价值低。如果把第二根控制棒插在离第一根控制棒较远的d2处，这时该处的中子通量密度比原来(没有第一根控制棒时)高，因此，第二根控制棒的价值比它单独插入堆芯时的价值高。同理，当第二根控制棒插入堆芯时，它也会使中子通量密度分布发生畸变

42、，因而影响到周围控制棒的价值。事实上，这样的影响是相互的，每一根控制棒的插入都将引起其它控制棒价值的变化。控制棒间的干涉效应从8-9上可以清晰地看出这种相互干涉效应的结果。图中虚线表示单根偏心控制棒价值的两倍，实线表示两根偏心控制棒同时插入堆芯时的价值。从图中可知，在两根控制棒相距较近时，两根同时插入堆芯时的总价值比它们单独插入时价值的总和要小；在两根控制棒相距较远时，两根棒同时插入所得的价值比单独插入所得到的价值的总和要大。考虑到控制棒的相互干涉效应，通常在设计堆芯时，应使控制棒的间距大于热中子扩散长度。控制棒间的干涉效应控制棒插入不同深度对堆芯功率分布的影响控制棒插入不同深度对堆芯功率分布

43、的影响图8-10给出了控制棒束插入不同深度时的轴向中子通量密度分布。当控制棒未插入时，堆内轴向中子通量密度呈正弦对称分布；随着控制棒逐渐插入，中子通量密度的峰值逐渐向底部偏移，且峰值也变大。在动力反应堆中，通常，新堆芯的初始剩余反应性都比较大。特别是在第一个换料周期的初期，堆芯中全部核燃料都是新的，这时的剽余反应性最大。如果全部靠控制棒来补偿这些剩余反应性，那么就需要很多控制棒，而每一控制棒(或棒束)都需要一套复杂的驱动机构。这非但不经济，而且在压力容器封头上要开许多孔，结构强度也不许可。如果全部依靠增加化学补偿毒物(如硼酸)浓度来满足要求，那么硼浓度可能超过限值，从而使慢化剂温度系数出现正值

44、。尤其是在轻水反应堆中，这个问题更为突出。为了解决这个问题，可以采用控制棒、可燃毒物与化学补偿毒物三种方式的联合控制，以减少控制棒的数目。可燃毒物控制-可燃毒物的作用可燃毒物控制-可燃毒物的作用可燃毒物控制-可燃毒物的作用可燃毒物的布置及其对反应性的影响1均匀布置情况均匀布置情况可燃毒物在堆芯中可以采用均匀和非均匀的布置。为了了解可燃毒物在维芯中的分布对反应性的影响，首先分析可燃毒物与慢化剂-燃料均匀混合的情况。为了简化起见，假设堆芯中没有中子泄漏，这时燃料和可燃毒物的核密度随时间变化的方程分别为：NF，NP和NFP分别为燃料、可燃毒物和裂变产物的核密度。 tNttdtdNtNtdttdNtN

45、tdttdNFFPafFPFPPPaPFFaF,可燃毒物的布置及其对反应性的影响对上述各式对t积分，可得 。是中子注量，其中，t dttFtFtFtNtFNtNtFNtNtFPaFPaFFPFPPaPPFaFF0,exp1exp0exp0可燃毒物的布置及其对反应性的影响假设堆芯中没有中子泄漏，而且慢化剂、冷却剂和结构材料等宏观吸收截面与时间无关，那么堆芯的有效增值因数可以用下式近似表示这样，将前面所得的核子密度代入上市便可求出不同的可燃毒物吸收截面情况下有效增殖因数随时间的变化曲线(假设初始时刻的有效增殖因数都相等并等于kex)。 FPaFPKaPaPFaFFFftNtNtNtNtk,可燃毒物

46、的布置及其对反应性的影响从图中可知，在反应堆运行刚开始的一段时间内，随着时间的增加，可燃毒物消耗所引起反应性的释放率比燃料燃耗所引起反应性的下降率要快得多，因此有效增殖系数上升很快。但是，当可燃毒物大量消耗后，每单位体积中含可燃毒物的核子数减少，这时可燃毒物消耗所引起反应性的释放率小于燃料燃耗所引起反应性的下降率，因此，有效增殖系数上升到某一最大值后又开始下降。可燃毒物的布置及其对反应性的影响从图中可知有效增殖系数偏离初始值的程度与可燃毒物的吸收截面a,P有关，其值越大，有效增殖因数偏离初始值也越大，这说明可燃毒物的消耗与堆芯中剩余反应性的减小不匹配。可燃毒物的布置及其对反应性的影响我们希望随着可燃毒物的消耗，在整个堆芯寿期内，有效增殖