从原子核到电网
你已经掌握了基础构件:原子核、裂变、结合能以及 E=mc²。
本模块探讨下一个问题:我们如何真正利用这些:安全、可靠地运行数十年?
核电站的核心,本质上是一种极其受控的烧水方式。关键在于“受控”一词。历史上每一次反应堆事故,都可追溯至控制失效的那一刻:无论是源于物理规律、工程缺陷,还是人为决策。
我们将从业链反应的数学原理出发,经过燃料循环与冷却剂水力学,深入探讨导致 SL-1、切尔诺贝利及三哩岛事故的特定物理失效机制。
这是社区大学水平的核工程课程。请准备好面对数字、方程以及实际的逻辑推理。
你目前了解什么?
在开始之前,让我们先进行一下校准。
中子生命周期
每个中子都有自己的故事
由裂变产生的中子在反应堆中穿行,最终会发生以下四种情况之一:引发另一次裂变、在不引发裂变的情况下被吸收、泄漏出反应堆,或者衰变(罕见情况:中子半衰期约为 10 分钟,对于反应堆物理而言速度太慢,影响微乎其微)。
一代中子数量与上一代中子数量的比值是增殖因子 k。
- k < 1:次临界:链式反应逐渐停止
- k = 1:临界:链式反应以恒定功率自我维持
- k > 1:超临界:功率正在增加
正常运行的反应堆恰好处于 k = 1 的状态。启动过程中的反应堆 k 值略高于 1。停机意味着 k 值被降至远低于 1。
为了理解什么因素控制 k,我们使用无限大反应堆(无泄漏)的四因子公式:
k∞ = η × ε × p × f
每个因子代表中子生命周期中的一个阶段。我们将逐一讲解。
四因子公式
k∞ = η × ε × p × f
η (eta),增殖因子:每个热中子被燃料吸收后产生的快中子平均数量。对于 U-235,η ≈ 2.07。对于 Pu-239,η ≈ 2.11。这是收益因子,即每次裂变能给我们提供多少新中子?
ε (epsilon),快裂变因子:用于考虑 U-238 中的快裂变。由 U-235 裂变产生的快中子在减速之前,可以引发丰度极高的 U-238 发生裂变。对于典型的压水堆(LWR)燃料组件,ε ≈ 1.03–1.07。它总是大于 1,是一个小的增益。
p:共振逃逸概率:指中子从快能级减速到热能级而未被 U-238 共振峰捕获的概率。U-238 在超热能区的特定能量(共振峰)处具有巨大的中子俘获截面。在典型的压水堆(LWR)中,p ≈ 0.75–0.80。这是最大的损失项。
f:热利用因子:指热中子被燃料吸收(而非被慢化剂、结构材料或控制棒吸收)的比例。f = Σ_fuel / Σ_total。在典型的压水堆(LWR)中,若无控制棒插入,f ≈ 0.71–0.75。
示例:η=2.07, ε=1.04, p=0.77, f=0.73 → k∞ = 2.07 × 1.04 × 0.77 × 0.73 ≈ 1.21
这意味着在无限大反应堆中,该燃料将处于高度超临界状态。实际反应堆是有限的:泄漏会使 k 降至 k∞ 以下。
理解四因子公式
反应堆操作员注意到,将控制棒插入更深会降低反应堆功率。控制棒由中子吸收材料(硼或铪)制成,并插入燃料区域。
六因子公式与泄漏
真实反应堆是有限的
四因子公式假设反应堆是无限的:没有中子逃逸。真实的反应堆具有边界,靠近表面的中子可能会流走并丢失。
六因子公式增加了两个非泄漏概率:
k_eff = η × ε × p × f × P_FNL × P_TNL
- P_FNL:快中子非泄漏概率:快中子在慢化之前不泄漏出去的概率。在大型轻水堆中通常为 0.97。
- P_TNL:热中子非泄漏概率:热中子在被吸收之前不泄漏出去的概率。在大型轻水堆中通常为 0.99。
泄漏是小反应堆更难达到临界状态的原因。小反应堆具有较高的表面积与体积比:按比例来说,更多的中子会到达边界并逃逸。
几何曲率 B² 量化了泄漏倾向。球体具有最低的表面积与体积比,因此在给定体积下具有最低的 B²:这解释了为什么弹芯是球形的(在给定质量下最大化 k_eff)。
在一座大型商用压水堆(1000 MWe)中,运行初期(BOL)在无控制棒插入的情况下,k∞ ≈ 1.2,但中子泄漏和控制棒的作用使运行期间的有效增殖因数 k_eff 恰好降至 1.000。
瞬发中子与缓发中子
为什么反应堆是可控制的
当 U-235 发生裂变时,大部分中子会立即出现:这些是瞬发中子,在裂变发生后 10⁻¹⁴ 秒内发射。所有裂变中子中约有 99.35% 是瞬发中子。
剩余的 0.65% 是缓发中子,由某些裂变产物在衰变过程中于数秒至数分钟后发射。平均延迟时间约为 13 秒,但不同中子组的延迟时间范围从 0.2 秒到 55 秒不等。
正是这一微小的缓发中子份额(对于 U-235,β = 0.0065)使得反应堆具有可控性。
瞬发临界是指仅靠瞬发中子即可使 k_eff ≥ 1 的状态:即无需依赖缓发中子份额。这是一种灾难性场景。在瞬发临界状态下,反应堆周期(功率增加 e 倍所需的时间)会从分钟级骤降至毫秒级。没有任何机械系统能够做出足够快的响应。
正常临界状态(k_eff = 1.000)依赖缓发中子来维持链式反应。有效中子代时间 ℓ_eff ≈ β/λ ≈ 0.0065/0.08 ≈ 0.08 秒:这一速度足够慢,使得机械控制棒能够调节功率。
瞬发临界的条件是:k_eff ≥ 1 + β,即对于铀-235,k_eff ≥ 1.0065。
我们将此称为超临界度 ρ ≥ β:反应堆处于“瞬发超临界”状态。
SL-1 事故(1961 年)以及 1986 年测试期间的切尔诺贝利 RBK-1000 反应堆均达到了瞬发临界状态。两者都在不到一秒内自我毁灭。
为什么缓发中子能救我们
反应堆周期与 Inhour 方程
反应度的测量
反应度 ρ 定义为 ρ = (k-1)/k。在临界状态下,ρ = 0。次临界:ρ < 0。超临界:ρ > 0。
单位美元 ($) 将反应度归一化为延迟中子份额:1$ = β ≈ 0.0065(针对 U-235)。当 ρ = 1$ = β 时,发生瞬发临界。
1 分 = 0.01美元。
反应堆周期 T 是功率增加 e 倍(≈2.718)所需的时间。小的正反应度插入会导致长周期(稳定、可控)。当接近瞬发临界时,周期会迅速缩短至零(不稳定)。
中时方程(Inhour equation)将反应度与反应堆周期联系起来。“Inhour”意为“逆小时”。该方程为:
ρ = (ℓ/T) + Σᵢ [ βᵢ / (1 + λᵢT)]
其中 βᵢ 和 λᵢ 分别是每组延迟中子的产额分数和衰变常数(对于 U-235 有 6 组),ℓ 是瞬发中子寿命。
对于小的正反应度(ρ << β),该方程给出 T ≈ β/(ρ·λ̄):反应堆周期较长且可控。
当 ρ → β(接近瞬发临界)时,T → 0:周期崩溃,功率呈爆炸性上升。
实际意义:启动需要正反应度。操作员需监视反应堆周期表。启动期间 30-60 秒的周期是正常的。若周期低于 10 秒,将触发 SCRAM(紧急停堆)。
为什么我们需要减缓中子速度
快中子与热中子
裂变产生的中子是快的:动能约为 1–2 MeV。U-235 在 1 MeV 能量下的裂变截面约为 1 靶恩(10⁻²⁴ cm²)。
将中子减速至热能量(室温下约 0.025 eV)& U-235 的裂变截面跃升至约 585 靶恩:高出近 600 倍。
这解释了为什么热中子反应堆(LWR、CANDU、AGR)使用慢化剂:一种能将中子从 MeV 减速至 eV 且不会吸收过多中子的材料。
热化通过弹性散射碰撞发生。每次碰撞都会将部分中子动能传递给靶核。每次碰撞的最大能量传递为:
ΔE/E = 4A/(1+A)²
其中 A 是靶核的原子质量。对于氢(A=1):ΔE/E = 1.0,中子可以在一次碰撞中传递全部能量。对于碳(A=12):ΔE/E = 0.28。对于铀(A=238):ΔE/E = 0.017,几乎无法减速。
这解释了为什么氢(在水中)是如此高效的慢化剂:它可以在约 18 次碰撞内使中子热化。碳(石墨)则需要约 114 次碰撞。但氢也会吸收中子(下文将详细讨论)。
慢化剂对比:H₂O vs. D₂O vs. 石墨
慢化剂的权衡
一种好的慢化剂必须:
1. 具有低原子质量(每次碰撞能量传递效率高)
2. 具有低的中子吸收截面(不要“偷走”你正在减速的中子)
对于普通氢而言,这两项要求是相互矛盾的。
轻水 (H₂O)
- 减速能力:非常高 (ξΣₛ ≈ 1.35 cm⁻¹)
- 吸收截面 (H):0.33 靶恩:显著
- 减速比 (ξΣₛ/Σₐ) ≈ 62
- 结果:虽然是一种优秀的慢化剂,但吸收的中子数量足以迫使您必须使用浓缩铀(3–5% U-235)来补偿。天然铀(0.71% U-235)无法提供足够的过剩中子来克服 H₂O 的吸收。
重水 (D₂O)
- 减速能力:低于 H₂O (ξΣₛ ≈ 0.18 cm⁻¹):需要更多次碰撞
- 吸收截面 (D):0.0005 靶恩:比 H 低 660 倍
- 慢化比 ≈ 5,500
- 结果:D₂O 几乎不吸收中子。因此可以使用天然铀(U-235 含量为 0.71%)。这解释了为什么 CANDU 反应堆使用天然铀燃料。
石墨 (C)
- 慢化能力:中等 (ξΣₛ ≈ 0.064 cm⁻¹)
- 吸收截面 (C):0.0035 靶恩:较低,但高于 D₂O
- 慢化比 ≈ 170
- 结果:可以使用天然铀或轻微富集铀。RBMK、Magnox 和 AGR 反应堆使用石墨。切尔诺贝利反应堆即为石墨慢化反应堆。
钠 (Na):非热中子慢化剂
- 钠冷快堆有意避免中子热化。直接使用快中子。不需要也不希望有慢化剂。快中子能谱允许增殖新的裂变材料(由 U-238 生成 Pu-239)。
CANDU 的优势
快堆:无需慢化剂
为什么钠冷快堆不需要慢化剂
快堆(钠冷快堆 SFR、铅冷快堆 LFR)有意维持快中子能谱。冷却剂(液态钠或铅)具有高原子质量且散射截面低:它不会使中子热化。
为什么要运行快堆?有两个原因:
1. 增殖:与热中子堆相比,快中子能更高效地将可裂变 U-238 转化为可裂变 Pu-239。快堆的增殖比(每消耗一个裂变原子所生成的新裂变原子数)可以超过 1.0,增殖堆产生的燃料多于其消耗的燃料。U-238 占天然铀的 99.3%,如果我们能实现增殖,这将是一个近乎取之不尽的燃料来源。
2. 嬗变:快中子可以裂变长寿命锕系元素(Am-241、Np-237、Cm-244),这些元素是乏核燃料中主要的长期辐射危害来源。在快堆中燃烧这些元素可将高放射性废物的寿命从 >100,000 年缩短至约 1,000 年。
权衡之处在于:钠与水及空气具有化学活性(钠火风险),快中子能谱意味着较低的裂变截面(单位中子效率较低),且工程设计更为复杂。
从矿山到燃料组件
燃料循环的前端
1. 开采:铀矿石通常含有 0.1–0.5%(按质量计)的铀。开采方式包括露天或地下开采,或原地浸出(ISL),即通过化学溶液在地下溶解铀。
2. 磨矿:矿石经过破碎和化学处理,生产出黄饼(U₃O₈):铀含量约为 85%(按质量计)。磨矿尾矿具有轻微放射性,需要谨慎处置。
3. 转化:黄饼被转化为六氟化铀(UF₆):在适度温度下呈气态。UF₆ 是浓缩过程的工作流体。反应过程:U₃O₈ + HF → UF₄ → UF₆。
4. 浓缩:天然铀由 99.3% 的 U-238 和 0.71% 的 U-235 组成。大多数反应堆需要 3–5% 的 U-235。两种商业化工艺:
气体扩散法:UF₆ 气体被泵送通过数千个多孔屏障。由于 U-235 比 U-238 略轻,因此 ²³⁵UF₆ 在每个级次的扩散速度比 ²³⁸UF₆ 快 1.004 倍。这需要数百个级次的级联以及巨大的电能(每分离单位约 2,400 千瓦时)。该工艺现已基本淘汰。
气体离心法:UF₆ 以 50,000–70,000 RPM 的速度旋转。较重的 ²³⁸UF₆ 富集在外壁;较轻的 ²³⁵UF₆ 富集在中心。每级分离因子约为 1.3(扩散法为 1.004)。耗电量约为扩散法的 1/50。现代标准方法。
浓缩度以分离功单位(SWU)衡量。从天然铀生产 1 公斤 5% 浓缩铀大约需要 8 SWU。
5. 燃料制造:浓缩 UF₆ 被转化为二氧化铀(UO₂)粉末,压制成陶瓷燃料芯块(直径约 1 厘米,高 1 厘米),在 1700°C 下烧结,堆叠进锆合金(Zircaloy)管中并密封:这些就是燃料棒。燃料棒被组装成燃料组件(例如,压水堆组件为 17×17 = 289 根燃料棒)。典型的 1000 MWe 压水堆拥有约 193 个燃料组件,总计约 80 吨铀。
浓缩水平及应用:
- 天然铀(0.71%):CANDU 堆、Magnox 堆
- 低浓缩铀(LEU,<20%):商业发电,轻水堆(LWR)使用 3–5%
- 高浓缩铀(HEU,≥20%):舰用反应堆(≥90%)、研究堆
- 武器级:≥90% U-235
离心机与气体扩散法
乏燃料与后处理
燃料循环的后端
在反应堆中运行 3–4 年后,乏燃料在物理上处于高温状态,具有极强的放射性,并且仍含有大量易裂变材料:
- 约 94% 为 U-238(U-235 已耗尽)
- 约 1% 为 U-235(仍具易裂变特性)
- 约 1% 为 Pu-239、Pu-240、Pu-241(由 U-238 俘获中子生成)
- 约 4% 为裂变产物(Cs-137、Sr-90、I-131 及其他约 200 种)
- <0.1% 为次级锕系元素(Am、Np、Cm)
一次性循环:美国政策:乏燃料首先储存在湿式乏燃料池(水屏蔽辐射并移除衰变热)中 5–10 年,然后转移到干式罐式储存。不进行后处理。高放射性废物(HLW)计划进行永久地质处置(位于尤卡山,目前处于停滞状态)。
PUREX 后处理(法国、英国、日本、俄罗斯):乏燃料溶解在硝酸中。溶剂萃取(煤油中的磷酸三丁酯)选择性提取铀和钚,留下裂变产物。回收的铀(后处理铀,RepU)可以重新浓缩。钚与贫化铀混合制成MOX 燃料(混合氧化物,约 5–7% PuO₂)。MOX 可将燃料资源延长约 10–20%。
武器级与反应堆级钚:
天然铀在反应堆中会产生 Pu-239。如果在反应堆中停留时间足够长,Pu-239 的中子俘获会产生 Pu-240。反应堆级钚(通常 Pu-240 含量 >18%)对武器制造存在问题,因为 Pu-240 具有较高的自发裂变率:它会导致枪式装置中的预爆炸(哑弹)。武器级钚需要较短的辐照时间(<3 个月)以限制 Pu-240 的积累。商业动力反应堆(燃料循环长达 18 个月以上)产生无法用于武器的反应堆级钚。这是一次性燃料循环中故意设置的扩散屏障。
微分与积分棒价值
控制棒的价值有多大?
棒价值是由插入控制棒引起的反应性变化。它不是恒定的:它取决于控制棒相对于中子通量分布的插入位置。
微分棒价值 (Δρ/Δx):指在特定位置处,每单位控制棒插入深度引起的反应性变化。它在通量最高的位置达到峰值:即堆芯中心。在顶部和底部(低通量区域),其值较低。
积分棒价值:指控制棒从完全抽出状态插入到某一特定深度时产生的总反应性变化。它形成一条S形曲线:在顶部(低通量)变化缓慢,在中心(通量峰值)变化迅速,在底部变化缓慢。
控制棒弹出事故:如果控制棒突然从堆芯中弹出(例如由于棒驱动机构故障),会在毫秒级时间内产生巨大的正反应性插入。其幅度取决于该棒的价值(从几毫分米到几美元,具体取决于棒的位置)。如果弹出棒的价值超过瞬发临界阈值(1$),则会发生瞬发临界偏离。
棒遮蔽 / 棒-棒相互作用:插入一根控制棒会降低局部通量,从而降低附近控制棒的价值。操作员在规划棒阵时必须考虑这种相互作用。
控制棒材料:硼-10(σₐ = 3,840 靶,在 0.025 eV 下),铪(σₐ = 102 靶,适中但燃耗缓慢,适用于长寿命控制棒),银-铟-镉合金(用于压水堆,Ag 提供快速响应,In 和 Cd 在燃耗过程中维持价值)。
氙中毒:无形的杀手
氙-135:已知最强的中子吸收体
氙-135 的热中子吸收截面为 260 万靶:远超任何核素。作为对比,铀-235 的裂变截面为 585 靶。氙-135 的单原子吸收能力约为其 4,400 倍。
生成:Xe-135 主要来源于 I-135(碘)的衰变,而 I-135 直接由裂变产生。仅有约 0.3% 的 Xe-135 直接来自裂变;约 95% 通过以下衰变链产生:
Te-135 → I-135(半衰期 6.6 小时)→ Xe-135(半衰期 9.2 小时)→ Cs-135
移除:Xe-135 通过两种过程被移除:(1) 放射性衰变(半衰期 9.2 小时),& (2) 中子吸收(被中子通量“烧尽”)。在高功率运行时,中子吸收是主要的移除机制。
碘坑(氙瞬变):
在稳态运行期间,Xe-135 的生成与移除达到平衡(在典型 PWR 中,氙价值 ≈ -2,500 pcm)。
当反应堆停堆时,Xe-135 的中子吸收停止。但 I-135 在随后的数小时内继续衰变为新的 Xe-135。停堆后 6–8 小时内,Xe-135 浓度上升:即碘坑。
如果剩余反应性不足,这可能导致反应堆暂时无法重启(无法克服氙中毒)。
与切尔诺贝利的关联:1986 年 4 月 26 日,切尔诺贝利 4 号机组的测试因电网需求延迟了约 9 小时。在此期间,氙不断积累。为了继续测试,操作员不得不抽出几乎全部控制棒以克服氙中毒。这使反应堆几乎没有停堆裕度:这是事故的关键前提条件。
为什么氙使反应堆在停堆后变得危险
钐中毒
Sm-149:长期存在的毒物
钐-149 是第二重要的反应堆毒物。其热中子吸收截面约为 41,000 靶恩(barns)。
生成链:Nd-149 → Pm-149(半衰期 53 小时)→ Sm-149(稳定)
与氙不同,Sm-149 是稳定的:它不会衰变消失。只能通过中子吸收将其移除。在稳态功率下,Sm-149 达到平衡浓度,代表约 -700 pcm 的反应性。
在停堆时:中子燃耗停止,但 Pm-149 继续衰变为 Sm-149。由于 Sm-149 是稳定的,它在停堆后约 100 小时内不断积累:增加约 -600 pcm 的负反应性。
在重启时:中子通量会燃耗掉过量的 Sm-149。钐中毒比氙中毒(没有等效的碘坑)要轻微,但必须在长期反应性管理中予以考虑。
综合来看,在停堆峰值时,氙和钐总共代表约 -3,000 至 -3,500 pcm 的反应性负担:在重启时,这必须通过控制棒提升或化学调节(PWR 中的硼酸)来平衡。
什么是反应性系数?
安全与不安全反应堆的区别
反应性系数是指每单位物理参数(温度、空泡份额、功率)变化所引起的反应性变化量。
负系数:随着功率增加,反应性降低:反应堆具有自限性。这是一种固有安全的设计。
正系数:随着功率增加,反应性增加:反应堆会放大扰动。这是一种潜在不稳定的设计。
反应性系数的符号决定了反应堆是固有安全的,还是需要主动干预以防止失控。这是反应堆设计中最重要的单一安全参数。
多普勒展宽:最重要的安全机制
多普勒反应性系数
多普勒展宽是一种量子力学效应:随着燃料温度的升高,U-238 原子核的热运动使其吸收中子的共振峰变宽。
在超热中子能区(1 eV 至 10 keV),U-238 具有巨大的共振吸收峰。在低温下,这些峰很窄:中子必须具有非常精确的能量才能被吸收。随着温度升高,展宽后的峰会吸收更宽能区内的中子。
对 p(共振逃逸概率)的影响:随着燃料温度升高 → U-238 共振峰变宽 → 在慢化过程中被俘获的中子增多 → p 减小 → k 减小 → 功率降低。
多普勒系数(α_D) 对于 U-235/U-238 燃料通常为 -1 至 -3 pcm/°C。这是一个强负值。
为何这是主要的安全机制:它即时起作用(温度变化速度等于热流速度:毫秒至秒级)。只要燃料中存在 U-238,它就始终存在。它不依赖于任何主动系统或操作员操作。它不会失效。
在任何反应性偏离(功率突然上升)中,多普勒效应会立即启动 & 在任何机械系统能够响应之前提供负反馈。这解释了为什么现代压水堆(LWR)燃料(燃料基体中 U-238 含量超过 95%)具有固有的自限性。
武器注记:纯 U-235 或 Pu-239 金属几乎没有多普勒反馈。这是武器使用高浓缩材料的原因之一:使动力堆安全的多普勒安全机制也会限制武器当量。
空泡系数:区分压水堆(LWR)与RBMK的关键
空泡系数与切尔诺贝利物理机制
空泡系数(α_v)是指每单位空泡份额(冷却剂中沸腾成蒸汽气泡的比例)变化所引起的反应性变化。
在轻水堆(PWR或BWR)中:
水同时充当冷却剂和慢化剂。如果水沸腾(形成空泡),慢化作用会减弱。慢化作用减弱 → 热中子数量减少 → 裂变减少 → 功率下降。此外,水会吸收部分中子:水减少意味着寄生吸收减少,这具有轻微的正面影响,但慢化作用的损失占主导地位。
结果:轻水堆中的空泡系数为负值(通常为 -100 至 -200 pcm/% 空泡)。冷却剂流失会自动降低功率。
在RBMK-1000(切尔诺贝利反应堆)中:
RBMK 使用石墨作为慢化剂,仅用水作为冷却剂。如果水沸腾:
- 慢化作用不变(石墨慢化剂不会改变)
- 水中的中子吸收减少(寄生吸收减少)
- 净效应:在低功率下具有正空泡系数
- 随着功率上升,水沸腾加剧,正空泡系数增加更多反应性,从而进一步提高功率:形成一个正反馈回路。
RBMK 中的正空泡系数幅度:在低功率且插入控制棒较少的情况下,α_v ≈ +4 至 +5 pcm/% 空泡。苏联设计师对此心知肚明,但向电厂操作员隐瞒了这一事实。
1986 年 4 月 26 日:切尔诺贝利 4 号机组在低功率(约 200 MWt,而额定功率为 3,200 MWt)下运行,大部分控制棒被抽出以克服氙中毒。在这种配置下:正空泡系数达到最大,控制棒价值最小,功率受氙抑制。当测试序列导致反应堆功率激增时,沸腾加剧,空泡系数增加了反应性,功率上升更快,导致更多沸腾:不稳定的正反馈。反应堆达到瞬发临界,并在约 3 秒内自我毁灭。
为什么 RBMK 在低功率下运行不稳定
慢化剂温度系数与功率系数
其他关键系数
慢化剂温度系数 (MTC):慢化剂温度每变化一度引起的反应性变化。在压水堆 (PWR) 中:随着水温升高,其密度降低 → 单位体积内的慢化剂减少 → 热化作用减弱 → 热中子数量减少 → k 值降低。在轻水堆 (LWR) 中,MTC 为负值(通常为 -20 至 -80 pcm/°C)。这是一项强制性的安全规范:美国核管理委员会 (US NRC) 法规要求 MTC 在任何时候都必须 ≤ 0。
燃料温度系数 (FTC):主要由多普勒展宽驱动(如上所述)。在轻水堆燃料中,该系数始终为强负值。
功率系数:所有来源引起的总反应性反馈,对应于单位功率变化。在设计良好的轻水堆中:该系数为强负值。功率上升 → 燃料温度升高(多普勒反馈)→ 慢化剂受热并产生空泡(MTC 及空泡反馈)→ 反应性降低 → 功率稳定。
综合效应:轻水堆反应堆具有固有的自调节特性。如果操作员什么都不做,会发现反应堆最终稳定在一个反馈机制使 k = 1.000 的功率水平上。这不是偶然:这是一个刻意的设计要求。
所有系数均为负值的反应堆,绝不会因热反馈事件而进入瞬发临界状态。轻水堆发生瞬发临界需要外部插入的正反应性大于瞬发临界阈值(>β ≈ 0.0065)。实际上,这意味着控制棒弹出或硼浓度快速稀释:这两种情况都在设计基准中进行了明确分析。
热移除:从燃料到冷却剂
保持燃料冷却
裂变主要通过裂变碎片的动能(~83%)和瞬发伽马辐射(~3%)产生热量,这些热量几乎全部沉积在燃料芯块内。裂变产物的β衰变(~4%)和伽马衰变(~4%)会随时间推移增加热量:这就是衰变热,它在反应堆停堆后仍会持续存在。
衰变热大致遵循“12法则”:停堆1分钟后,衰变热约为运行功率的1%。1小时后:~0.4%。1天后:~0.2%。1周后:~0.07%。一台3,000 MWt反应堆在停堆1分钟后的衰变热约为30 MWt:如果失去冷却,足以熔化堆芯。这解释了为什么紧急堆冷却系统(ECCS)如此关键。
热流路径:燃料芯块 → 燃料棒包壳(锆合金) → 冷却剂水 → 蒸汽发生器(压水堆)或直接产生蒸汽(沸水堆)
温度分布:在满功率运行时,压水堆(PWR)燃料中心线温度可达~900–1,200°C。锆合金包壳表面:~300–350°C。冷却剂主体:~290–325°C。从芯块中心到冷却剂的陡峭梯度意味着功率的微小增加会导致燃料温度的大幅上升:& 巨大的多普勒反馈。
关键热学极限:燃料中心线温度必须保持在UO₂熔点(~2,865°C)以下。包壳温度必须保持在锆合金氧化阈值(~1,200°C)以下,超过此温度锆会与蒸汽发生放热反应:Zr + 2H₂O → ZrO₂ + 2H₂。正是这一反应产生了在福岛1、3、4号机组发生爆炸的氢气。
偏离核态沸腾(DNB)
临界热流密度极限
在压水堆(PWR)中,冷却剂在约155巴的压力下保持液态(沸点约345°C)。微小的蒸汽气泡在包壳表面成核,并被流动带走,即核态沸腾,这实际上是一种极佳的热传递方式。
如果局部热流密度超过临界值(临界热流密度,CHF),气泡会聚合成包裹在燃料棒周围的连续蒸汽膜。这层蒸汽膜是绝缘体。燃料产生的热流无法通过蒸汽被带走:包壳温度会迅速飙升。这就是偏离核态沸腾(DNB)或临界热流密度超限。
DNB的后果:如果没有迅速恢复流动,包壳温度将升至1,200°C,此时锆合金开始氧化,随后趋向熔化(约1,850°C)。燃料颗粒会飞散,裂变产物释放到冷却剂中。
MDNBR(最小DNB比):局部临界热流密度与实际热流密度之比,在堆芯最受限的位置进行评估。在任何时候都维持MDNBR ≥ 1.3的安全限值(即1.3倍的DNB裕度)。此限值限制了反应堆的最大功率和流量条件。
两相流:在沸水堆(BWR)中,整体沸腾是有意为之的:堆芯在双相流(水+蒸汽)状态下运行。BWR中的等效限值是临界功率比(CPR)或最小临界功率比(MCPR) ≥ 1.2。
堆芯温度分布:轴向热流密度遵循轴向中子通量分布(在新鲜堆芯中大致为截断余弦形)。峰值通量(及最高DNB风险)位于堆芯中平面。径向峰值位于中心组件。 热通道因子(Fq或F∆H)量化了峰值局部功率比堆芯平均值高多少:在PWR中通常为2.5–3.0。
为何DNB设定了关键安全限值
PWR 和 BWR:主导设计
轻水反应堆
轻水反应堆(LWR)占全球商业核电装机容量的约 85%。
压水反应堆 (PWR)
- 一回路:水在约 155 巴(15.5 MPa)、约 290–325°C 下运行:压力高于沸点,保持液态
- 热交换器:蒸汽发生器将热量从一回路传递至二回路
- 二回路:水在约 60 巴下运行,产生约 280°C 的蒸汽以驱动汽轮机
- 优势:主回路放射性水从不接触汽轮机。维护更容易。
- 功率:每台机组 900–1,700 MWe。热效率约 33%。
- 示例:西屋 AP1000、法国 EPR、俄罗斯 VVER
沸水堆 (BWR)
- 直接循环:水在反应堆压力容器内沸腾,压力约 75 巴(约 290°C)。蒸汽直接输送至汽轮机。
- 无需蒸汽发生器:结构更简单,对压力容器压力要求较低
- 汽轮机具有轻微放射性(蒸汽中夹带的裂变气体):需要屏蔽及远程维护
- 除控制棒外,通过调节再循环流量来控制功率(流量增加 → 空隙率降低 → 慢化作用增强 → 功率升高)
- 被动安全性:较低的压力意味着储存的能量较少,ECCS(紧急核心冷却系统)设计更简单
- 热效率约为 33%,与压水堆(PWR)相似
- 示例:GE BWR/6、ABWR、ESBWR
VVER(Vodo-Vodyanoi Energetichesky Reaktor,俄语“水-水能量反应堆”):苏联/俄罗斯压水堆设计。采用水平式蒸汽发生器,而西方压水堆采用垂直式。燃料组件几何形状为六边形,而西方压水堆为方形。现代 VVER(VVER-1200)符合西方安全标准。
CANDU 和 RBMK:压力管设计
压力容器之外的替代方案
CANDU(Canada Deuterium Uranium,加拿大氘铀反应堆)
- 水平压力管内含燃料和冷却剂(高压 D₂O),周围是位于 calandria(慢化剂容器)内的低压 D₂O 慢化剂
- 在线换料:在反应堆满功率运行期间更换燃料,无需停堆。每台燃料装换机可单独访问每一根压力管。这使得容量因子达到100%,无需因换料而停堆(压水堆必须停堆约18个月进行换料)
- 天然铀燃料(UO₂):无需浓缩。CANDU的中子经济性使其成为可能。
- 还可接受MOX燃料、钍燃料及乏轻水堆燃料(回收)
- 所有反应性系数均为负值:固有稳定性
- 示例:CANDU-6(700 MWe)、ACR-1000(采用轻水冷却剂的先进设计)
RBMK-1000(Reaktor Bolshoy Moshchnosti Kanalnyy:大功率管式反应堆)
- 苏联设计:石墨慢化剂,垂直压力管内的轻水冷却剂
- 大型(1,000–1,500 MWe),低浓缩铀,在线换料
- 致命的物理缺陷:在低功率且控制棒抽出状态下存在正空泡系数(在反应性系数部分中有详细描述)
- 额外的设计缺陷:石墨端头效应,控制棒带有石墨端头。当控制棒从完全抽出状态开始插入时,首先会将堆芯底部的水排出(移除了寄生吸收),随后吸收体部分才进入活性区。因此,插入控制棒以执行紧急停堆(SCRAM)操作时,最初会引入一个短暂的正反应性脉冲,这与预期的效果完全相反。
- 这两个缺陷的结合导致了切尔诺贝利灾难。
- 所有幸存的RBMK反应堆均经过改造,以降低正空泡系数并重新设计控制棒。它们仍然是独特的苏联设计,西方没有同类机型。
第四代反应堆概念
超越现有机队
第四代核能系统国际论坛(GIF)确定了六种反应堆概念用于开发,目标是在2030年及以后部署:
熔盐堆 (MSR):燃料溶解在熔融氟化物盐(LiF-BeF₂ 或 NaF-ZrF₄)中。没有固体燃料,也没有会熔化的燃料包壳。具有被动排空至冻结塞的特性,如果失去电力,冻结塞会熔化,盐液排入次临界几何形状。在大气压力下运行(约 650°C)。可实现钍增殖。
液态氟化物钍堆 (LFTR):一种使用 Th-232/U-233 增殖循环的特定 MSR 设计。钍的丰度约为铀的 3 倍。U-233 由 Th-232 产生(Th + n → Pa-233 → U-233)。LFTR 产生的长寿命锕系废物极少。倡导社区热情高涨;但工程挑战(高温腐蚀、盐化学控制)依然显著。
钠冷快堆 (SFR):液态钠冷却剂,快中子能谱,具有增殖或锕系嬗变的潜力。挑战:钠与水及空气反应(需要惰性气氛)。现有实例:BN-800(俄罗斯)、Superphénix(法国,已退役)、Monju(日本,事故后关闭)。EBR-II(美国)在 1986 年通过故意诱导的流量丧失证明了被动安全性,反应堆在无需 SCRAM 的情况下自行安全停堆。
铅冷快堆 (LFR):铅或铅铋冷却剂。铅不与水或空气反应(与钠不同)。沸点高(1,740°C),无需加压。可能实现自然循环冷却。挑战:铅非常重,且在高温下对钢具有腐蚀性。俄罗斯潜艇反应堆曾使用 Pb-Bi 冷却剂。
超临界水堆 (SCWR):水处于其临界点以上(374°C,221 bar),单相,焓值极高。热效率可能达到 ~44%,而当前 LWR 约为 ~33%。结合了 BWR 的简单性与高效率。在超临界条件下存在重大的材料挑战。
超高温堆 (VHTR):氦气冷却,石墨慢化,出口温度 700–950°C。可通过热化学循环实现制氢。TRISO 燃料颗粒(陶瓷包覆微球)即使在无主动冷却的情况下也能保留裂变产物。实例:HTR-PM(中国,2023 年运行)。
选择反应堆类型
朗肯循环
将热能转化为功
核电站本质上是一种蒸汽动力厂。卡诺效率定理设定了效率的上限:
η_Carnot = 1 - T_cold/T_hot(温度单位为开尔文)
压水堆(PWR)蒸汽工况:T_hot ≈ 280–290°C(553–563 K),T_cold ≈ 30–40°C(303–313 K)
η_Carnot = 1 - 308/558 ≈ 0.45(45%)
实际热效率约为 33%:差距源于实际循环中的不可逆性(涡轮损失、泵功、传热温差、蒸汽中的水分)。
朗肯循环的各个阶段:
1. 给水泵:将过冷液态水加压至锅炉压力(输入功较小)
2. 蒸汽发生器 / 锅炉:来自反应堆的热量将水转化为蒸汽(输入热量较大)
3. 高压涡轮(HP):蒸汽膨胀,驱动涡轮轴旋转,压力与温度降低
4. 汽水分离器/再热器:在涡轮级之间对湿蒸汽进行干燥并重新加热
5. 低压涡轮(LP):蒸汽进一步膨胀至凝汽器压力
6. 凝汽器:蒸汽通过冷却水(河流、海洋、冷却塔)冷凝回液态
7. 给水加热器:利用从涡轮级抽出的蒸汽预热给水(再热循环:通过减少锅炉热输入和凝汽器热损失来提高循环效率)
为何核电站效率约为33%,而燃煤/联合循环燃气轮机(CCGT)为40–43%:核电站的蒸汽温度和压力显著低于现代化石燃料电厂的蒸汽。燃煤电厂可实现600°C的蒸汽温度(超临界);而压水堆(PWR)受压力容器限制及燃料温度限制,蒸汽温度仅约为280°C。较低的T_hot(热端温度)→ 较低的卡诺效率极限 → 较低的可实现效率。
为何核电站承担基荷:燃料成本几乎全部为前期投入(浓缩+制造)。可变运营成本(每MWh的燃料成本)非常低(约7美元/MWh,而天然气约为30美元/MWh)。资本成本非常高。这使得核电站成为所有可调度发电机组中边际运营成本最低的:以100%出力连续运行在经济上是最划算的。在电力市场的优先排序(merit order)中,核电站通常被优先调度。
核能效率与联合循环燃气
点动力学方程
功率随时间的变化
点动力学方程描述了中子数量(因此也是反应堆功率)随反应性变化的时变行为:
dN/dt = [(ρ - β)/ℓ]·N + Σᵢ λᵢ·Cᵢ + S
dCᵢ/dt = (βᵢ/ℓ)·N - λᵢ·Cᵢ
其中 N = 中子数量,ρ = 反应性,β = 总缓发中子份额,ℓ = 瞬发中子寿命,Cᵢ = 第 i 组缓发中子先驱核浓度,λᵢ = 第 i 组衰变常数,S = 外部中子源。
对于小的反应性引入(ρ << β),解给出稳定周期:
T ≈ β / (ρ · λ̄)
其中 λ̄ 是延迟中子的有效衰变常数(约 0.08 s⁻¹)。对于 ρ = 0.01$ = 0.0001(1 分):
T ≈ 0.0065 / (0.0001 × 0.08) ≈ 813 秒:非常稳定。
对于 ρ = 0.50$ = 0.00325:
T ≈ 0.0065 / (0.00325 × 0.08) ≈ 25 秒:仍然可控。
瞬发跃变近似:对于突然引入的反应性,中子群会在瞬发时间尺度(约 10 µs)内瞬间跃升至新水平,随后较慢的延迟中子动力学才开始主导。瞬发跃变因子为 1/(1-ρ/β)。对于 ρ = 0.50$,功率瞬间跃变因子为 1/(1-0.5) = 2,随后以 25 秒的周期上升。这解释了为什么即使是微小的反应性引入也会引起立即可见的功率响应。
反应堆启动与棒落测试
接近临界
启动程序:反应堆从次临界状态开始。控制棒被缓慢抽出。随着控制棒的抽出,k 从下方逐渐接近 1.000。
1/M 图(次临界倍增):在达到临界状态之前,监测来自启动源的能通量计数率。对于具有外部源 S 且倍增因子 M = 1/(1-k) 的次临界反应堆:
计数率 ∝ M = 1/(1-k)
绘制 1/(计数率) 与棒位的关系图,会得到一条在临界点外推至零的曲线。操作员在接近临界状态时绘制 1/M,并外推以预测临界棒位。如果 1/M 的下降速度快于预期,说明临界状态比预测的更近:操作员必须放慢操作速度。
提棒试验:将控制棒从已知位置插入堆芯。突然引入的负反应性导致功率呈指数下降。通过测量衰减速率,可以计算出棒价值。
初始衰减遵循:P(t) = P₀·exp(-t/T_negative)
其中 T_negative 取决于棒价值。棒价值越大 = 衰减越快。
逆周期仪:控制室显示反应堆周期(正 = 功率增加,负 = 功率减少)。在正常启动期间,周期保持在 30–60 秒。如果周期低于 20 秒,触发警报。如果周期低于约 10 秒,自动触发紧急停堆(SCRAM)。
临界事故(历史案例):在早期核计划中的临界事故(洛斯阿拉莫斯“龙”实验、SL-1 反应堆、日本 Tokaimura 事故)中,共同因素是反应性不受控制地增加,超过了瞬发临界阈值。在洛斯阿拉莫斯,物理学家使用裸钚半球:任何导致它们靠得过近的失误都会引发瞬发临界。Louis Slotin 在 1946 年幸存于此类事故;Harry Daghlian 在 1945 年未能幸存。
SL-1:控制棒弹出引发的瞬发临界事故(1961年)
SL-1:全球首例致死性反应堆事故
SL-1(一号低功率固定式反应堆)是美国陆军位于爱达荷国家实验室的一座小型实验反应堆。1961年1月3日,三名操作员正在进行维护工作:手动重新连接控制棒。
事故经过:中央控制棒在约0.5秒内被手动抽出约67厘米(26英寸)。这一次单棒抽出引入了约3–4个美元单位($3-4)的正反应性:远超1$的瞬发临界阈值。
物理机制:当 ρ > β = 1$ 时,达到了瞬发临界状态。点动力学方程表明,在瞬发临界状态下,稳定周期坍缩至瞬发中子寿命(约10微秒)。功率在约4毫秒内上升了约10,000倍。
能量释放:在前4毫秒内发生了约1.3 × 10¹⁷次裂变。冷却剂瞬间爆炸性汽化。蒸汽爆炸将一股水柱以约160公里/小时的速度向上推射,携带反应堆容器盖及附着的控制棒。一名操作员被控制棒刺穿并钉在天花板上。
原因:为什么单根控制棒的反应性价值高达3-4个美元单位?在SL-1反应堆中,仅三根控制棒控制整个反应堆,每根棒的反应性价值极高。仅中央棒的价值就约为5$。此外,在反应堆寿命初期,其装载了大量新鲜燃料,处于无氙状态,即最大反应性状态。
课程要点:反应堆设计应确保任何单根控制棒弹出都不会导致瞬发临界。控制棒棒值限制现已成为标准设计要求。SL-1事故直接促成了对独立停堆系统的要求以及对单根控制棒棒值的限制。
三哩岛:失水事故与操作员困惑(1979年)
三哩岛2号机组:一场系统性事故
三哩岛2号机组(压水堆,906 MWe,位于宾夕法尼亚州)于1979年3月28日发生堆芯部分熔毁。未发生瞬发临界:反应堆本身成功执行了紧急停堆(SCRAM)。该事故是一起失水事故(LOCA)与操作员失误的结合。
起始事件:稳压器上的一个卡滞开启的导阀操作式安全阀(PORV)。当压力升高时,阀门正确开启,但随后未能重新关闭。一回路冷却剂通过开启的阀门持续流失。
关键混淆点:控制面板上的一个指示灯显示PORV已收到关闭信号,但这只是一个信号指示器,而非位置指示器。阀门实际处于开启状态,但操作员误以为它已关闭。他们看到“稳压器液位上升”(水位上升是因为蒸汽空间正在填充,这是压力丧失的症状,而非水量过剩),并据此判断系统过满。因此,他们减小了应急堆芯冷却系统的注入流量。
堆芯状况:在大约2小时20分钟内,堆芯部分裸露。由于缺乏冷却,衰变热(请记住:即使在停堆状态下,衰变热也约为满功率的1%)使燃料温度升至1,200°C以上。锆合金被蒸汽氧化(Zr + 2H₂O → ZrO₂ + 2H₂)。大约45%的燃料熔化并迁移至压力容器底部。
包容成功:尽管堆芯严重受损,安全壳建筑阻止了裂变产物的显著释放。大约17居里的放射性碘和250万居里的惰性气体被释放:虽然数量可观,但远低于灾难性水平。无辐射致死案例。
课程:人因工程成为核安全中强制考量的因素。控制室进行了重新设计。关键阀门的信号指示器被位置指示器所取代。应急操作程序被重写,以基于症状(而非基于事件)的方式进行响应。美国核管理委员会进行了重组。
切尔诺贝利:正空泡系数 + 操作员越权操作(1986年)
切尔诺贝利:完美的物理风暴
切尔诺贝利核电站4号机组(RBMK-1000型,热功率3,200兆瓦)于1986年4月26日在一次安全测试中自我毁灭。这场事故是存在缺陷的反应堆设计与一系列将反应堆置于最危险配置的操作员决策共同作用的结果。
测试:汽轮机惰转测试旨在证明,惰转中的汽轮机能够提供足够的功率,以在柴油发电机启动所需的约75秒内驱动应急冷却泵。此前该测试已尝试过三次并均告失败。这是第四次尝试。
前置条件(单独看每一项都有危险;合在一起则致命):
1. 氙中毒:9小时的延迟(电网需求)导致氙积累。为了继续测试,操作员抽出了几乎所有控制棒。运行技术规范要求堆芯内至少保留15根控制棒;而在事故发生时,仅有6至8根插入。
2. 低功率运行:反应堆功率约为 200 MWt(约为额定功率的 6%)。在此功率范围内,RBMK 反应堆的空泡系数呈最强的正反馈特性。
3. 冷却剂泵全流量运行:为进行测试,额外启动了泵,导致过冷水流量增加:这抑制了沸腾,并需要抽出更多控制棒以维持功率。
4. AZ-5 控制棒设计缺陷:当控制棒从完全抽出状态完全插入时,石墨棒头会短暂增加正反应性,随后吸中子部分才进入堆芯。
事故序列:
- 测试开始。汽轮机节流阀关闭。冷却剂流量下降。水开始沸腾。
- 正空泡系数增加反应性。功率开始上升。
- 操作员意识到情况并按下 AZ-5(紧急停堆:所有控制棒插入)。
- 所有 211 根控制棒的石墨棒头同时进入堆芯,短暂增加了约 3$ 的正反应性:这与预期效果完全相反。
- 在约3秒内,功率达到估计的30,000 MWt(约为额定功率的10倍),在某些燃料通道中可能高达30,000倍。
- 瞬发超临界偏离。燃料碎裂导致蒸汽爆炸。2–3秒后,发生第二次、规模更大的爆炸(可能是更多燃料中的瞬发超临界)。
- 1,000吨重的反应堆顶盖被掀飞。石墨和燃烧的燃料散落在整个场区。
为何这在RBMK中发生且不可能在LWR中发生:
- LWR中的负空泡系数意味着沸腾会降低功率,而不是增加功率
- LWR的控制棒没有石墨端头:紧急停堆(SCRAM)总是增加负反应性
- LWR的燃料是富集的:维持功率不需要极低的控制棒插入度
对比事故分析
纵深防御
为什么反应堆拥有多重独立安全屏障
现代核安全建立在纵深防御原则之上:设置多重独立屏障,即使前序屏障失效,每一道屏障都能独立地预防或减轻事故后果。
压水堆(LWR)的五道屏障:
1. 燃料基体:UO₂ 陶瓷在极高温度下仍能保留约 97% 的裂变产物
2. 燃料包壳:锆合金管容纳燃料棒,防止裂变产物释放到冷却剂中
3. 一回路压力边界:反应堆压力容器、稳压器及一回路冷却剂管道:15 厘米厚钢材
4. 安全壳建筑:钢筋混凝土 + 钢衬里,设计用于承受内部蒸汽爆炸及外部飞机撞击
5. 隔离区:围绕场址的土地使用限制
应急系统(主动式):
- ECCS(应急堆芯冷却系统):高压和低压注入系统,在失去一次冷却剂时向堆芯注水
- SCRAM(安全控制棒斧头人:原术语字面意思即如此):所有控制棒在<2秒内插入
- 安全壳喷淋:水雾在事故后冷却并降低安全壳内的压力
被动安全(三代及以上设计:AP1000、ESBWR):
- 反应堆上方的重力供水水箱:无需泵或交流电源
- 利用水的密度差进行自然循环冷却:无需泵
- 安全壳内的被动自催化复合器(PAR):无需点火即可将 H₂ + O₂ → H₂O,防止氢爆炸
- AP1000 专为 72 小时宽限期设计,期间无需操作员干预
福岛教训:AP1000 的被动安全系统正是针对福岛事故中的失效模式而专门设计的。福岛事故中,主动紧急冷却系统(ECCS)泵因海啸淹没发电机而失去交流电源。被动系统无需外部电源。
设计一座安全反应堆
融会贯通
你现在已掌握核工程所需的完整物理工具包:四因子公式、临界性、缓发中子、慢化、燃料循环、反应性系数、热工水力学及事故分析。