4.3 伽马谱定量分析

4.3 伽马谱定量分析

本节摘要：定量伽马谱分析的流水线共五步：认峰、扣本底、算净面积、除以效率与分支比、除以计数时间。本节把五步走通一遍，算出真实量级的活度结果，讲清效率曲线的驼峰形状从哪里来，并给出最低可探测活度的算法——它决定"没检出"三个字的底气。这一节是从"看谱"到"出数"的最后一公里，环境监测、活化分析、废物分类都用同一条流水线。

4.2 节解决了"峰在哪"，本节解决"峰有多高、代表多少活度"。方法本身朴素得令人安心：数峰下的格子，再除以一连串已知系数。

五步流水线

第一步认峰：在能量刻度上核对候选核素的全部特征峰，单峰定核素是外行做法，铯-137 要 661.7 千电子伏一支独大且无伴生异常。第二步扣本底：峰两侧取本底区，线性内插扣除连续本底。第三步净面积：峰区内计数减去本底贡献。第四步系数链：除以探测效率、除以该能量伽马的分支比。第五步换活度：除以计数时间，必要时再除样品质量。整条链写成一行式：活度等于净面积除以效率、分支比、时间。

# 演算 4.3-1：铯-137 全能峰定量的完整一算
net = 98000.0        # 661.7 keV 全能峰净面积（已扣本底）
eff = 0.032          # 该几何条件下 661.7 keV 的绝对效率（3.2%）
I_g = 0.851          # 每次衰变放出该伽马的分支比
t_live = 600.0       # 活时间（秒）
A = net / (eff * I_g * t_live)
print(f"样品活度 = {net:.0f}/({eff}×{I_g}×{t_live}) = {A:.0f} 贝可")
# 若本底区有 1200 计数：净面积与本底不确定度见 4.4 节
print(f"合 {A / 1e3:.1f} 千贝可；对应质量可由比活度换算")
# 输出：
# 样品活度 = 98000/(0.032×0.851×600.0) = 5998 贝可
# 合 6.0 千贝可；对应质量可由比活度换算
​

一行式里每个除数都是一个学过的概念：效率是 4.2 谱仪的本事，分支比是 3.2 讲过的发货折扣，时间是 3.1 衰变定律的钟。核化学的知识在这里真正首尾相接。

效率曲线：驼峰的物理

探测效率随能量不是单调的，而是先升后降的驼峰：低能端窗户吸收与自吸收压住效率，高能端伽马穿透力太强干脆穿过晶体，峰在两百千电子伏上下最高。定量分析必须用与样品同几何条件测出的效率曲线，换几何（样品高度、容器壁厚）效率就作废——环境实验室为不同样品杯各测一条效率曲线，就是这个原因。常用标准源是混合伽马核素（镅-241、铯-137、钴-60 等多点联测），对数坐标下拟合成光滑曲线。

# 演算 4.3-2：效率曲线的形状直觉（相对值示意）
# 绝对效率受晶体与几何影响，这里只演示趋势：窗吸收-峰值-穿透三段
energies = [60, 150, 300, 662, 1332, 2500]          # 千电子伏
eff_rel = [0.35, 0.95, 1.00, 0.62, 0.30, 0.12]      # 相对驼峰值
for E, e in zip(energies, eff_rel):
    print(f"{E:4d} keV: 相对效率 {e:.2f}")
# 输出：
#   60 keV: 相对效率 0.35
#  150 keV: 相对效率 0.95
#  300 keV: 相对效率 1.00
#  662 keV: 相对效率 0.62
# 1332 keV: 相对效率 0.30
# 2500 keV: 相对效率 0.12
​

图 4-3 全能峰定量五步流水线

"没检出"也要有数学：判别下限

低水平测量里最常写的三个字是"未检出"，但未检出不等于零——它等于"低于我的判别下限"。下限由本底涨落决定：本底计数 B 越大、计数时间越短，下限越高。库里公式给出：判别下限约等于 2.71 加 4.65 乘根号 B（净计数口径）。想让报告里的下限值更低，办法只有三个：延长时间、改善几何、加强屏蔽压低本底——每一条都是真金白银的实验决策。

# 演算 4.3-3：判别下限与"未检出"的表述规范
import math
B = 800.0                          # 峰区本底计数
LD = 2.71 + 4.65 * math.sqrt(B)    # 净计数判别下限（库里公式）
eff, I_g, t_live = 0.032, 0.851, 600.0
A_LD = LD / (eff * I_g * t_live)
print(f"本底 {B:.0f} 计数 → 判别下限 {LD:.0f} 净计数")
print(f"折合活度下限 ≈ {A_LD:.0f} 贝可")
print(f"样品净面积若为 90 计数（低于下限）：报告写未检出，附下限值 {A_LD:.0f} 贝可")
# 输出：
# 本底 800 计数 → 判别下限 134 净计数
# 折合活度下限 ≈ 8 贝可
# 样品净面积若为 90 计数（低于下限）：报告写未检出，附下限值 8 贝可
​

精密测量再补两笔：符合相加与自吸收

级联伽马（比如钴-60 的两连跳）在探头附近有小小的"近水楼台"：两支伽马被同一个晶体先后收下的概率随几何靠近而上升，导致全能峰面积偏离效率曲线的预测——符合相加修正，样品越贴近窗越要补。另一笔是自吸收：含水样品自身会吃掉一部分低能伽马，同样几何、不同密度的基质效率并不同，环境水样常按"加同等密度模拟基质的标准源"或做传输修正来摆平。这两笔修正单看都不大（百分之几到十几），但痕量测量的结论常常就吊在百分之几上——4.4 的统计误差之外，它们是谱定量的系统账。

实战案例：一管沉降水的完整定量

背景：环境站收到大沉降后的一管降水样品，需要给出碘-131 与铯-137 的活度浓度。

操作：样品浓缩后装入标准几何杯，上锗谱仪计数四万秒。谱上 364.5 千电子伏（碘-131）与 661.7 千电子伏（铯-137）各有一峰；分别取峰区扣本底算净面积，查同几何效率曲线取值，除以各自分支比与活时间，再除以水样体积。

结果：碘-131 得 2.4 贝可每升（衰变修正到采样日），铯-137 得 0.9 贝可每升，两值都高于判别下限一个量级，结论可信。

解读：碘-131 的存在本身就是事故新鲜释放的指纹——八天半衰期让它在环境里只能存活两三个月（6.5 节细说）；铯-137 则可能是历史沉降的再悬浮。同一张谱读出了时间与来源两个维度，这正是谱学定量的魅力。

变式练习：若把计数时间从四万秒压到四千秒，判别下限抬高根号十倍，碘-131 的 2.4 贝可每升还算不算"检出"？动手用 4.3-3 的公式改一改——你会看到定量报告里"计数时间"一栏绝不是形式。

峰面积是链条留在谱仪里的脚印，五步流水线把脚印换算成活度；从本节起，"测量出一个数"对你而言应当意味着"同时报出它的不确定度"——那是 4.4 的主题。