巴斯克地区辐射监测网络毕尔巴鄂辐射站降水所致剂量率增量的计算

《Sensors》：Calculation of the Increase in Dose Rate Due to Precipitation at the Bilbao Radiological Station of the Basque Country Radiological Surveillance Network

　　

巴斯克自治区辐射监测网络持续监测环境剂量率水平以探测潜在的辐射异常。然而，在降雨事件期间，尽管没有发生辐射事故，却观察到剂量率增加，且偶尔超过预设的报警阈值。本研究提出一个预测模型，用于估算毕尔巴鄂监测站降水事件引起的剂量率增量。该模型考虑了雨滴对大气中氡子体（radon progeny）的清除作用及其随后在探测器周围水平表面上的沉积。模型公式中包含了大气传输、放射性衰变、冲刷过程（washout）和雨水累积动态。模型采用有限差分格式进行数值求解，并用监测站的实验测量数据进行验证。结果显示，在降水期间，计算剂量率增量与测量值之间具有良好一致性（相关系数、BIAS、RMSE），表明观测到的变化主要归因于短寿命氡子体的沉积。

**论文解读：巴斯克地区辐射监测网络毕尔巴鄂辐射站降水所致剂量率增量的计算**

**研究背景与问题**
环境辐射监测网络在早期探测潜在放射性释放（如事故或未经授权的排放）中发挥关键作用。然而，天然降水事件常导致环境γ剂量率暂时升高，这源于大气中氡（²²²Rn）衰变子体（如²¹⁴Pb和²¹⁴Bi）被雨滴冲刷并沉积于地面。这种自然现象可能使监测值接近甚至超过预设报警阈值，引发误报警和不必要的验证程序。巴斯克自治区辐射监测网络在毕尔巴鄂站多次观测到降水期间剂量率升高，但未发现人工放射性释放。为区分自然气象效应与潜在辐射异常，需建立物理过程模型定量解释降水引起的剂量率响应。目前，已有研究主要聚焦于观测分析或清除过程表征，缺乏可操作性预测模型。因此，本研究旨在开发并验证一个数值模型，以估算毕尔巴鄂站的降水诱导剂量率增量。

**研究内容与结论**
研究人员构建了一套基于质量平衡方程的模型，描述氡子体从大气传输到探测器周围水平表面的完整过程：大气分布与传输、降水冲刷、雨水相运输、沉积与累积、放射性衰变及γ辐射场产生。模型采用有限差分数值格式求解，并通过蒙特卡罗（MCNP）代码计算沉积核素到剂量率的转换因子。利用毕尔巴鄂站10分钟间隔的剂量率和降水实测数据，对50次独立降水事件进行回顾性验证，重点分析包含3至6个连续降水间隔的代表性案例。结果显示，模拟剂量率增量与实测值吻合良好，多数事件相关系数>0.9，均方根误差（RMSE）和偏差（BIAS）较低。模型成功再现了降水期间剂量率上升、降水后因衰变而下降的典型时间演化模式。研究表明，观测到的剂量率升高主要归因于短寿命氡子体（特别是²¹⁴Pb和²¹⁴Bi）的冲刷沉积，模型可作为减少误报警、优化监测协议的有效工具。该论文发表在《Sensors》。

**关键方法**
1. **大气传输与冲刷模型**：基于大气传输方程，忽略湍流扩散和沉降项，仅考虑放射性衰变和降水冲刷，用恒定冲刷系数描述氡子体从气相向水相的转移。
2. **雨水相与沉积过程**：建立雨水相放射性活度平衡方程，考虑雨滴下落、破碎后沉积于水平表面（探测器屋顶），并引入由流体动力学导出的排水系统模型。
3. **剂量率计算**：利用MCNP蒙特卡罗模拟计算沉积核素（²¹⁴Pb、²¹⁴Bi）发射的γ粒子通量，结合空气比释动能转换因子和有效剂量转换因子（ISO照射几何）得到环境剂量率增量。模型采用有限差分法在时间和空间上离散求解，时间步长秒级，空间步长米级。

**研究结果**
**3.1 三降水间隔案例**
两个降水事件（雨量分别为0.7、2.2、0.2 L m^?2和0.7、0.2、0.1 L m^?2）的模拟与实测剂量率增量对比显示，模型精确再现了上升和下降阶段，相关系数分别为0.99和0.70，RMSE分别为0.00025 μSv h^?1和0.0012 μSv h^?1，BIAS接近零。

**3.2 四降水间隔案例**
雨量0.5、0.2、0.9、0.4 L m^?2的事件中，标准模拟（正常排水）相关系数0.90，RMSE 0.0018 μSv h^?1，BIAS +0.0006 μSv h^?1；假设排水受阻的替代模拟虽改善下降阶段匹配，但总体统计指标变差，表明标准排水假设更优。

**3.3 五降水间隔案例**
雨量0.8、0.4、0.2、0.1、0.1 L m^?2的事件，模拟与实测曲线趋势一致，相关系数0.96，但BIAS ?0.0165 μSv h^?1，RMSE 0.028 μSv h^?1，显示模型在事件末期轻微低估。

**3.4 六降水间隔案例**
雨量0.3、0.2、0.7、0.5、0.1、0.1 L m^?2的事件，标准模拟相关系数0.915，RMSE 0.00233 μSv h^?1，BIAS ?0.00123 μSv h^?1；假设无排水的替代模拟相关系数升至0.972，RMSE降至0.00133 μSv h^?1，BIAS几乎为零，表明排水阻塞可能解释了部分偏差。

**讨论与结论**
讨论部分指出，模型总体成功再现了降水诱导的剂量率演化，高相关系数和低误差支持其反映主导物理机制（冲刷、沉积、衰变）。但存在若干不确定性来源：10分钟间隔内降水强度均匀假设可能引入偏差，尤其在快速变化的事件中；恒定冲刷系数未考虑雨滴谱、大气稳定度等动态参数；雨量计检测阈值可能遗漏微弱降水事件。此外，五间隔案例的较大误差暗示了云底高度估计、局部排水条件等未显式因素。本研究为概念验证性回顾验证，未来需基于更大数据集评估模型稳健性，并整合动态冲刷参数化和更精确的降水测量。

翻译结论部分：本研究提出一个数值模型，成功模拟了毕尔巴鄂辐射监测站降水事件期间环境剂量率的增加。模型将这些增加归因于短寿命氡子体的大气冲刷和后续沉积。模型与实测剂量率增量对比显示良好一致性，具有高相关系数和相对较低的BIAS及RMSE值，表明模型能可靠预测降雨事件期间的剂量率变化。尽管模型为观测到的环境剂量率增加提供了有力解释，但仍存在一些局限性：10分钟间隔内均匀降水的假设以及恒定冲刷系数的使用可能引入轻微偏差；雨量计灵敏度和排水问题对氡子体沉积的影响也会影响预测精度。未来改进可聚焦于引入更精确的降水数据、可变冲刷系数及更精细的排水系统模拟。尽管存在这些局限，该模型证明是解释降雨事件期间剂量率数据的有价值工具，可提供氡子体冲刷对剂量率贡献的定量估算，有助于区分自然辐射波动与潜在辐射事故。模型减少误报警和优化监测协议的能力使其成为辐射监测网络的重要补充。在实际应用中，该模型可集成到实时监测系统，自动预测降水事件期间的剂量率增加，提高数据解释效率和响应速度。此外，通过适应区域气象条件和冲刷特征差异，模型可应用于广泛地理区域，增强全球辐射监测的准确性和有效性。尽管本研究为回顾性验证，但所提出的建模框架旨在支持环境辐射监测系统的实时运行，通过估算降水诱导的氡子体冲刷对剂量率增量的预期贡献，帮助操作员区分自然大气效应与潜在辐射异常，从而更高效地解释监测数据。然而，当前验证基于有限代表性降水事件，应视为概念验证而非全面运行评估。未来需使用覆盖不同季节、降雨类型和气象条件的更大数据集进行额外验证，以评估模型稳健性、不确定性和可迁移性。此外，应解决站点特定运行参数的整合以及模型在实时辐射监测系统中的实施问题。