铀矿开采氡对区域环境210Pb库存贡献的上限

《Journal of Environmental Radioactivity》：An upper bound on uranium-mining radon contributions to regional environmental 210Pb inventories

　　

铀矿开采遗留物的氡释放通常根据吸入暴露进行评估。这里，一个独立的问题是，持续遗留释放产生的²¹⁰Pb是否可能扰动用于定年并解释土壤、沉积物、泥炭和其他环境档案的库存。本文界定了这种可能性，并检验了区域过量²¹⁰Pb解释所依据的自然-大气来源假设。由于衰变链质量平衡守恒的是原子而非活度，从3.8天放射性核素转换为22.3年放射性核素将相关活度源项减少λ₂₁₀/λ₂₂₂ ≈ 4.7 × 10^-4。对于一个有记录的修复期源规模为10¹⁴ Bq y^-1，即使在10⁶ km²上完全捕获所有产生的²¹⁰Pb，也仅产生约0.05 Bq m^-2 y^-1，而背景自然沉积范围为50–200 Bq m^-2 y^-1。一个互补的源比诊断表明，在代表性大陆足迹内，单个遗留物贡献约10^-4的自然氡源，而最大的历史区域合计贡献约10^-3–10^-2。因此，有记录的遗留释放不太可能可检测地扰动区域平均²¹⁰Pb沉积或由此衍生的库存。不能排除局部源-受体事件，但来自含²²⁶Ra残留物的颗粒物和侵蚀介导的输入在设施和流域尺度上是更合理的途径，可能需要明确考虑和有针对性的、基于途径的监测。

环境²¹⁰Pb（铅-210）库存广泛用于示踪大气沉积和为近现代沉积物、土壤、泥炭等档案提供年代学框架，其基础假设是过量²¹⁰Pb来自大气中的自然大陆²²²Rn（氡-222）呼出。铀矿开采遗留物（如尾矿设施、废石堆、矿石堆、矿井通风口等）是持续的人为增强²²²Rn源，单个遗留物年释放量可达10¹⁴–10¹⁵ Bq y^-1，且以区域集合体而非孤立点形式存在。由于每个释放的²²²Rn原子最终产生一个²¹⁰Pb原子，自然产生一个核心问题：持续遗留排放是否构成一个足以扰动档案所记录库存的区域性大气²¹⁰Pb源项。研究人员在《Journal of Environmental Radioactivity》上发表的论文给出了否定答案，并定量论证了该结论的宽裕安全边际。

研究人员首先通过衰变链质量平衡分析了原子守恒与活度转换的关系，得出²²²Rn释放活度与²¹⁰Pb活度源项之间的转换因子λ₂₁₀/λ₂₂₂≈4.7×10^-4。随后，基于文献中有记录的源项数据（包括德国Wismut复合体修复期约10¹⁴ Bq y^-1、1990年代初约10¹⁵ Bq y^-1，匈牙利Pécs尾矿池覆盖前2.5–5.0×10¹⁴ Bq y^-1，以及美国活性尾矿区历史合计10¹⁵–10¹⁶ Bq y^-1），进行了完全捕获上限计算和区域源比诊断，评估了遗留源对区域²¹⁰Pb沉积的潜在贡献，并探讨了达到显著扰动所需的优先耦合因子。研究还指出，颗粒物和侵蚀介导的²¹⁰Pb输入是更值得关注的途径。

主要关键技术方法包括：
- 衰变链质量平衡分析：通过原子守恒推导，将²²²Rn释放活度乘以衰变常数比λ₂₁₀/λ₂₂₂≈4.7×10^-4转换为²¹⁰Pb活度源项。
- 完全捕获上限计算：假设给定接收面积内所有²¹⁰Pb原子被捕获，计算空间平均沉积通量的上限。
- 区域源比诊断：定义遗留源与自然源释放活度比δ，并引入优先耦合因子E，定量评估达到显著贡献所需的条件。
- 源项数据来源：文献中监测报告和历史估计，包括Wismut（德国，2024年修复期监测数据及1990年代初历史估计）、Pécs尾矿池（匈牙利，基于通量和面积推导）、美国活性尾矿区（EPA 1983年合计估计）。

研究结果按论文小标题总结如下：

**2. 质量平衡守恒什么：从氡释放到²¹⁰Pb活度**
通过衰变链分析，研究人员指出一个释放的²²²Rn原子最终产生一个²¹⁰Pb原子，但活度不守恒；活度等于原子数乘衰变常数。因此，年释放X Bq的²²²Rn对应的²¹⁰Pb活度源项为S₂₁₀ = X × (λ₂₁₀/λ₂₂₂) ≈ 4.7×10^-4 × X。对于Wismut规模（10¹⁴ Bq y^-1），S₂₁₀≈4.7×10¹⁰ Bq y^-1。该因子解释了自然大陆²²²Rn呼出（约20 mBq m^-2 s^-1）与观测到的²¹⁰Pb沉积范围（50–200 Bq m^-2 y^-1）之间的一致性。

**3. 有记录的源项和一级筛选**
结合文献中四种源项类型（监测、历史估计、推导、合计），研究人员计算了完全捕获上限。对于Wismut修复期释放（10¹⁴ Bq y^-1），在10⁶ km²面积上完全捕获所有²¹⁰Pb原子，平均沉积通量上限为47 mBq m^-2 y^-1，稳态库存上限为1.5 Bq m^-2，远低于自然沉积下限（50 Bq m^-2 y^-1）和自然库存（1600–6400 Bq m^-2）。即使考虑历史最大合计（10¹⁵–10¹⁶ Bq y^-1），完全捕获上限也仅约0.5–5 Bq m^-2 y^-1，仍低于自然范围一个数量级。

**4. 区域源比诊断及显著性所需的优先耦合**
定义源比δ = E × X / (Φ_nat × A)（其中E为优先耦合因子，Φ_nat≈6.3×10⁵ Bq m^-2 y^-1为自然²²²Rn释放通量），研究人员计算了无优先耦合（E=1）时δ值：Wismut修复期在2×10⁶ km²足迹上δ≈8×10^-5；Pécs尾矿池δ≈2–4×10^-4；历史最大合计δ≈10^-3–10^-2。要达到1%贡献（δ=0.01），Wismut需E≈1.3×10²，Pécs需E≈25–51，历史合计需E≈1–10²。文献未发现如此高的优先耦合证据。近场自由平流情景估算（1–30 km半径内平均沉积通量10–0.35 Bq m^-2 y^-1）仍低于自然范围。

**5. 更合理的途径：颗粒物和侵蚀介导的²¹⁰Pb和²¹⁰Po输入**
研究人员指出，含²²⁶Ra（镭-226）残留物（如未覆盖尾矿）可直接释放含²¹⁰Pb的颗粒物，绕过氡途径。假设残留物²¹⁰Pb活度约10 Bq g^-1，局部粉尘沉积率0.1–1 g m^-2 y^-1，可产生1–10 Bq m^-2 y^-1的²¹⁰Pb沉积，达自然下限的百分之几到五分之一。归因需结合²²⁶Ra、²¹⁰Pb、²¹⁰Po活度及粒度、矿物学示踪剂。覆盖和抑尘可控制该途径。

**6. 对环境档案和遗留场地监测的影响**
对于环境档案，自然-大气来源假设对于有记录的遗留源项是稳健的，单个遗留物安全边际三到四个数量级，历史最大合计安全边际约一到两个数量级。对于遗留场地监测，研究区分了区域档案保护与局部覆盖性能：氡释放清单对吸入评估和源头控制仍然必要，但无需额外考虑²¹⁰Pb区域档案贡献；局部监测应关注裸露残留物的颗粒物和侵蚀途径，包括²¹⁰Pb、²¹⁰Po、²²⁶Ra、粉尘沉积、覆盖状态等。

**讨论与结论**
讨论部分强调了源项历史演变、近场滞留效应、优先耦合的物理限制，以及颗粒物途径的优先级。研究结论指出：持续的²²²Rn释放产生的²¹⁰Pb活度仅为释放活度的4.7×10^-4；对于有记录的源项10¹⁴–10¹⁶ Bq y^-1，完全捕获上限仅0.05–5 Bq m^-2 y^-1，远低于自然范围50–200 Bq m^-2 y^-1；因此有记录的遗留释放不太可能产生可检测的区域平均²¹⁰Pb沉积扰动；局部源-受体事件不能排除，但更合理的途径是来自含²²⁶Ra残留物的颗粒物和侵蚀介导输入。