成都理工 核反应堆物理分析9 动力学

13:43:451/41核反应堆物理分析核反应堆动力学13:43:452/41核反应堆安全运行的基础在于成功的控制中子通量密度或反应堆功率在各种情况下随时间的变化。第七章,我们讲了燃料和裂变产物同位素成分随时间的变化以及它们对keff的影响。由于这些量随时间的变化很缓慢，所以很容易对其进行控制并使反应堆维持在一定功率下运行。13:43:453/41反应堆启动、停堆或功率调节时的控制棒的移动等情况下,反应堆的keff发生迅速变化。此时反应堆成为超临界或次临界，而中子通量密度随时间急剧变化。这种变化以秒为单位来量度。了解这种中子通量密度在偏离临界状态下的瞬态变化特征，对反应堆的控制和安全运行是及其重要的。13:43:454/41本章讨论由于有效增殖因子或反应性的迅速变化所引起的反应堆内中子密度随时间的瞬态变化特性。中子密度的瞬态变化会引起反应堆内功率、温度等的瞬态变化，而这些参数的变化，又会引起反应性的变化，从而又引起中子密度的变化。即在反应堆的瞬态过程中，存在反应性反馈效应，为了简单起见本章暂不涉及反应性的反馈效应。13:43:455/419.1缓发中子的作用在以前的反应堆临界计算只需考虑中子产生与消失保持平衡即可，认为所有裂变中子为瞬发中子。然而，反应堆动力学研究功率或中子通量密度的瞬态时间特征，就必须考虑缓发中子的产生相对于裂变时刻的延迟。缓发中子占总裂变中子份额很小(235U0.65%)，但缓发时间很长，它对反应堆动态特征有重要的影响。13:43:456/41为了说明这一问题，假设所有裂变中子为瞬发中子，则堆芯内中子密度的变化率为：13:43:457/41假设在t0时，keff=1，在t=0时keff有一阶跃变化，上式积分得定义反应堆周期（T）为：反应堆内中子密度变化e倍所需要的时间。则有：反应堆内中子密度随时间的变化13:43:458/41不考虑缓发中子时,压水堆内中子的平均寿命就等于瞬发中子的寿命,既l≈10-4秒。引入=0.001的反应性扰动，反应堆周期：T=0.1秒。1秒内，堆功率将增大e10（22000）倍。核反应无法控制！！13:43:459/41上面结论是基于所有裂变中子为瞬发的假设，忽略了缓发中子。虽然，缓发中子所占比例很少，但其缓发时间却非常长（表1-6）。所以缓发中子对中子寿命所起的作用不能忽视。根据表1-6第i组缓发中子寿命应等于ti+l，这里ti是第i组缓发中子先驱核的平均寿命。考虑缓发中子后裂变中子的平均寿命为：13:43:4510/41利用表1-6的数据可以计算出考虑缓发中子时压水堆的中子平均寿命为l≈10-1秒，引入=0.001的反应性扰动，反应堆周期：T=100秒。1秒内，堆功率将增大e0.01（1%）倍。核反应完全可以控制！！13:43:4511/41缓发中子虽然份额很少，但缓发时间较长，缓发效应大大增加了两代中子之间的平均时间间隔，从而延迟了中子密度的变化率。所以缓发中子效应在研究瞬态过程和反应堆控制时不可忽略。反应堆控制正是利用缓发中子的作用才得以实现的。13:43:4512/419.2点堆中子动力学方程核反应堆动力学方程我们从单群中子扩散方程出发推导反应堆动力学方程。在研究功率的瞬态特征时要考虑每组缓发中子产生的数据延迟及其效应，以及它对动态过程的影响。根据（3-34）单群中子扩散方程为：13:43:4513/41如认为所有裂变中子为瞬发中子则：若考虑缓发中子的效应：13:43:4514/41考虑缓发中子的效应单群中子扩散方程为：13:43:4515/41每一个先驱核只放出一个缓发中子，就是堆内每秒每单位时间体积内第i组缓发中子先驱核产生率。先驱核通过衰变而消失，消失率等于，考虑缓发中子效应的联立方程为：以上两方程便是反应堆的动力学方程。可以推广到多群情况。反应堆的动力学方程描述瞬态过程堆芯中子通量密度随空间和时间的变化。一般只能用数值方法求解。13:43:4516/41点堆动力学方程在动态过程中,空间效应不是主要的,我们感兴趣的是中子通量密度随时间的变化,我们可以用点堆模型来处理.假定不同时刻中子通量密度在空间中的分布形状不变，也就是说堆内各点中子密度随时间的变化涨落是同步的，堆内中子好像没有线度尺寸一样，可以把它看作一个集总参数的系统来处理，在此基础上得到的模型称为点堆模型。13:43:4517/41假定中子通量密度和先驱核浓度可以用空间形状因子,与时间相关的幅函数和的乘积来表示:若堆芯偏离临界状态不远,并且先驱核的浓度分布具有与中子通量密度分布相同的分布函数,将以上表达式代入反应堆的动力学方程可得。),(trCi)(rgi)(tn)(tCi13:43:4518/4113:43:4519/41定义中子每代时间以上方程可以改写为：13:43:4520/41这就是考虑缓发中子效应后的反应堆动态方程，通常称为点堆动力学方程。这是耦合的一阶微分方程组。keff或反应性是时间函数。如考虑功率和温度对反应性的反馈作用，keff或反应性还是中子通量密度函数，所以以上点堆动力学方程是非线形方程。13:43:4521/41点堆模型的适用范围假定不同时刻中子通量密度在空间中的分布形状不变，也就是说堆内各点中子密度随时间的变化涨落是同步的，堆内中子好像没有线度尺寸一样，可以把它看作一个集总参数的系统来处理。点堆模型的主要限制在于它不能描述与空间有关的动力学效应，如反应性的局部扰动和过渡过程中中子通量密度空间分布随时间的快速畸变。13:43:4522/41例如，对于一个均匀平板反应堆，分为三个活性区，在t=0时刻，I区引入一个正反应性阶跃，反应堆超临界；随后在0.01s内反应性线性下降到-。平板状堆芯当反应性局部阶跃变化时快中子通量密度空间分布的计算结果13:43:4523/41点堆模型的缺陷：多数物理量不考虑空间效应，或中子通量密度分布所用的形状因子是取未受扰动时的形状函数，且不随时间变化。由于实际中子通量密度的显著畸变和倾斜，用点堆模型求得的中子通量密度将大大低于实际的峰值。点堆模型仅适用于局部扰动不大，或者空间效应不太重要的情况！13:43:4524/419.3阶跃扰动时点堆模型动态方程的解在运行过程中,反应堆的提升功率、降低功率和停堆过程分别对应于反应堆引入正和负反应性的操作。在不考虑反馈效应，点堆方程为一阶线性常系数微分方程组，假定其解的形式为：反应性的阶跃变化13:43:4525/41代入点堆方程（9-18）和（9-19）并整理可得反应性方程：或13:43:4526/41它是一个关于的7次代数方程,在给定的反应堆特性参数下,由它可以确定出7个可能的值。但求解直接该方程却非常困难。可以用图解法研究方程的根的分布却非常方便。0时：有6个负根和1个正根。0时：有7个负根。13:43:4527/41在反应性阶跃变化的情况下，点堆模型动态方程（9-18）和（9-19）是线性的，所以方程的一般解由的所有7个解所形成的线性组合给出，即：系数根据初始条件确定。用图解法确定反应性方程的根13:43:4528/41作为反应堆keff阶跃变化的例子，设在初始状态突然引入=0.001的阶跃变化，中子的寿命l≈10-4秒。由方程（9-24）确定出7个值，其中为正值，其余6个均为负值。中子密度时间响应为13:43:4529/41时方程（9-25）中只有为正值，其余为负值。反应堆周期将按增长。无缓发时T=0.1秒时所有项均指数衰减。秒550182.0/11T阶段扰动下相对中子密度水平随时间变化的曲线13:43:4530/419.4反应堆周期13:43:4531/41中子密度按指数规律变化e倍所需的时间称为反应堆周期：当引入反应性为正时，Ｔ为正，中子密度随时间指数增长当引入反应性为负时，Ｔ为负，中子密度随时间指数衰减13:43:4532/41堆内中子通量密度增长一倍所需的时间称为反应堆倍周期或倍增周期。用TD表示：所以反应堆周期的符号和大小直接反映堆内中子增减变化速率，特别是在启动或功率提升过程中，对反应堆周期的监督十分重要。周期过小，可能导致反应堆失控。一般周期限制在30s以上。反应堆还有周期保护系统。2)/exp(/)(0TTnTnddTTTd693.02ln13:43:4533/41反应堆周期可以由反应性方程（9-23）或（9-24）确定，以代入则可得：或13:43:4534/41反应性方程（9-23）或（9-24）和方程（9-40）、（9-41）都是表示反应堆周期和反应性之间的方程式。给定反应性便可确定周期Ｔ。反之，它是由测量得到的反应堆周期来确定反应性的理论依据。式中右端第一项是表示瞬发中子的作用，而第二项则表示缓发中子对反应堆周期的影响。若=0，则两式变为这便和（9-3）不考虑缓发中子时的反应堆周期相等。13:43:4535/419.4.2不同反应性引入时反应堆的响应特性13:43:4536/4113:43:4537/41时，反应堆达到临界尚需缓发中子作出贡献，因而反应堆特性在很大程度上由先驱核衰变的时间决定，称为缓发临界。时，称为瞬发超临界，此时即使不考虑缓发中子，有效增殖因子也会大于1，只靠瞬发中子就能使链式反应不断进行下去，缓发中子在决定周期方面不起作用。反应堆功率以瞬发中子决定的极短周期快速增长。13:43:4538/41利用（9-13）式可以求得瞬发临界的条件，不考虑缓发中子（9-13）式变为：临界时中子密度不随时间变化，所以瞬发中子临界条件为：上式说明去处缓发中子部分，即仅依靠瞬发中子达到临界，称为瞬发临界。13:43:4539/4113:43:4540/41当为很大的负反应性时稳定周期将接近于，即约等于80s。如果由于引入大的负反应性而突然停堆，则中子通量密度迅速下降，而在短时间内瞬变项衰减之后，中子通量密度将按指数规律下降，其周期约为80s，即大约每184s功率下降一个量级。13:43:4541/41反应堆停堆时，中子通量密度需下降10个数量级以上，其关闭时间要求至少30min，反应堆设计需要考虑。有些反应堆，即使停堆以后很长时间，还会由积累的裂变产物衰变放出射线，然后通过（,n）反应继续产生中子。这样的过程会延长停堆的时间。