6.2 裂变产物链与衰变热

6.2 裂变产物链与衰变热

本节摘要：控制棒落下、链式反应停止，但堆芯里已经生产出的裂变产物仍在按自己的半衰期衰变，持续放出百分之几的功率——衰变热。本节给出 Way-Wigner 经验式与完整演算，解释停堆一小时的衰变热为什么足以熔毁堆芯、氙-135 尖峰为什么卡住重启窗口。这一节也是对切尔诺贝利与福岛事件最诚实的物理注脚：堆芯事故的本质是衰变热管理失败。

链式反应可以一键叫停，放射性衰变不行——它不讲纪律，只讲半衰期。堆芯停堆后的一切应急时间表，都从这句话开始。

衰变热：不熄火的余温

裂变瞬间释放的能量里约百分之六随缓发贝塔与伽马"分期到账"——这些期票的承兑方就是裂变产物链。堆运行得越久，产物库存越足，停堆后的衰变热越大。工程上的标准估算 Way-Wigner 公式给出：停堆后 t 时刻的衰变功率与满功率之比约为 0.066 乘 t 的负零点二次方幂（t 以秒计，运行时间足够长）。幂律下降很慢——停堆一天仍有约百分之零点七，停堆一年仍有千分之二。

# 演算 6.2-1：停堆衰变热的幂律（3000 MWt 电厂为例）
import math
P0 = 3000.0    # MWt
def decay_heat(t_s):
    return 0.066 * t_s ** -0.2      # P/P0，Way-Wigner 近似
for lbl, t in [("1 s", 1), ("1 h", 3600), ("1 d", 86400), ("1 a", 3.15576e7)]:
    print(f"停堆后 {lbl}: 衰变功率 ≈ {decay_heat(t) * 100:.2f}% 满功率 "
          f"= {P0 * decay_heat(t):.1f} MW")
# 输出：
# 停堆后 1 s: 衰变功率 ≈ 6.60% 满功率 = 198.0 MW
# 停堆后 1 h: 衰变功率 ≈ 1.28% 满功率 = 38.5 MW
# 停堆后 1 d: 衰变功率 ≈ 0.68% 满功率 = 20.4 MW
# 停堆后 1 a: 衰变功率 ≈ 0.21% 满功率 = 6.3 MW
​

三十八兆瓦——相当于一个小城镇的供暖功率，在没有泵送的情况下足以在数小时内把裸露的堆芯加热到锆水反应与熔化。福岛事故的物理内核正是这组数字：海啸切断了电源，衰变热的余量撑过了前几个小时，然后堆芯开始裸露。"停堆"这个词在核工程里的意思是"换一种方式发热"。

氙-135：链条给重启设的卡

3.2 节算过碘-135 到氙-135 的链：氙-135 的中子吸收截面冠绝所有核素，停堆后它失去"边生边烧"的平衡，浓度先涨后落，十余小时见顶——这就是停堆后反应性井。重启窗口被它卡住：井底附近强行提棒，反应性预算可能不够；等氙衰变过去（约一天半），功率才能重新爬坡。氙毒管理是运行教科书的第一课，物理出处全在 Bateman 方程。

# 演算 6.2-2：停堆后氙-135 存量的时间曲线（相对满功率平衡值）
import math
LI = math.log(2) / 6.57     # 碘-135 衰变常数，每小时
LX = math.log(2) / 9.14     # 氙-135 衰变常数，每小时
cap = 0.60                  # 满功率下氙的中子烧损清除率（每小时）
P = 1.0                     # 归一化的裂变产额供给速率
I0 = P / LI                 # 满功率平衡碘库存
X0 = P / (LX + cap)         # 满功率平衡氙库存（衰变与烧损共同清除）
for t in [0, 5, 9, 11, 20, 36]:
    # 停堆后烧损归零：氙初值项加碘衰变喂入项
    X = X0 * math.exp(-LX * t) + P * (math.exp(-LI * t) - math.exp(-LX * t)) / (LX - LI)
    print(f"停堆后 {t:>2d} h：氙-135 库存为满功率平衡值的 {X / X0:.2f} 倍")
# 输出：
# 停堆后  0 h：氙-135 库存为满功率平衡值的 1.00 倍
# 停堆后  5 h：氙-135 库存为满功率平衡值的 2.83 倍
# 停堆后  9 h：氙-135 库存为满功率平衡值的 3.20 倍
# 停堆后 11 h：氙-135 库存为满功率平衡值的 3.19 倍
# 停堆后 20 h：氙-135 库存为满功率平衡值的 2.46 倍
# 停堆后 36 h：氙-135 库存为满功率平衡值的 1.04 倍
​

停堆约十一小时后氙库存冲到满功率平衡值的三倍以上，三十六小时后才回落接近平常态——运行规程里"停堆两天内谨慎重启"的依据就是这条曲线。切尔诺贝利事故的当晚，操作人员正是在氙井底附近把反应性推到危险的边缘。

图 6-2 停堆后的三条时间线

实战案例：福岛一号机组的事故时间轴复盘

背景：2011 年 3 月 11 日地震后，福岛各机组自动停堆；约一小时后海啸淹没应急发电机。

操作：按本节物理走时间轴——停堆时刻衰变热约百分之六到七（机组刚换料后停堆数月的一号机组库存较低）；失去全部交流电后，堆芯冷却靠电池驱动的仪表与余热排出系统；电源耗尽后，衰变热把水位压到燃料顶部，锆合金包壳在高温下与水汽反应放出氢气。

结果：一号机组在停堆约一天内出现燃料损伤与氢气爆炸，与衰变热幂律给出的时间预算吻合。

解读：复盘的价值在"如果"：任何一路冷却水源——哪怕消防水枪——只要在衰变热曲线上方接管，事故即止步。衰变热物理不复杂，复杂的电力与组织系统在极端灾害下同时失效才是教训本体。

变式练习：把演算 6.2-1 的停堆前运行时长改成"刚启动一周"（产物的短寿命主角占比不同），Way-Wigner 偏保守还是偏危险？答案是偏保守——运行时间越长衰变热越大，公式按无限长运行取上限。工程公式的保守方向，是读公式必看的一栏。

乏燃料水池：衰变热账本最诚实的应用

乏燃料出堆后先泡在水池里，水一身兼三职：屏蔽伽马、载走衰变热、隔着几米水层把剂量率压到池边可走近。水池的冷却系统按"全炉存料加刚出堆最热一炉"设计，余量留给停循环事故。演算 6.2-1 的幂律在这里换一个用法：出堆三十天的组件衰变功率约是刚出堆的几分之一，水池分区摆放"新热夹老凉"，就是给幂律曲线排的座位表。福岛之后各国给水池加了冗余补水与监测，逻辑仍然只有一句——只要库存按半衰期衰减，冷却就要按库存的最不利时刻设计。

停堆不是链条的终点，是它换一种方式的继续；衰变热的幂律、氙毒的尖峰、乏燃料的百年，三个闹钟提醒你：核安全管的全是"停不下来的东西"。