建筑几何结构对室内γ剂量率影响的敏感性分析

《Buildings》：Sensitivity Analysis of the Influence of Building Geometric Structure on Indoor Gamma Dose Rate

　　

摘要：为探究实际应用中建筑材料对室内辐射环境的影响特征，并定量评估建筑设计参数对室内γ剂量率的调控作用，研究人员采用蒙特卡罗（Monte Carlo）方法建立了标准房间内空气吸收剂量率的三维计算模型，并分析了室内剂量场的分布。结果表明，标准封闭房间内的空气吸收剂量率呈现中心低、边缘高的梯度分布特征，由于多个近场源的叠加，角落和边缘区域出现剂量峰值。对建筑几何参数的敏感性研究表明，当房间高度固定时，中心剂量率随房间体积呈非单调变化：随着房间体积从2.8 m3扩大至约25.2 m3，中心剂量率下降，而当进一步扩大至280 m3时，与标准房间相比则逐渐升高。在固定建筑面积条件下，中心剂量率随房间高度增加呈单调递减趋势。引入门、窗等真实非对称开口后，剂量场呈现出显著的非对称分布特征。研究人员分析了窗户尺寸和位置对剂量的影响，并建立了一个典型的偏置门与居中窗的住宅布局。研究发现，在开口延伸区域形成了低剂量通道，导致室内极端剂量点偏离几何中心。

一、研究背景与问题提出

近年来，全球建筑工程需求显著增长，水泥作为主要建筑材料每年用量达数亿吨，而人们大部分时间生活在各类水泥建筑物中。因此，水泥原料中的废料及其产品带来的天然放射性危害不容忽视。建筑材料中的天然放射性主要来源于原材料（如岩石、土壤、工业废渣）中广泛存在的天然放射性核素，其中²²⁶Ra（铀系）、²³²Th（钍系）和⁴⁰K因对公众辐射剂量贡献最大而被列为重点关注对象。这些核素主要通过外照射（人体受到体外强穿透γ射线的照射）和内照射两种途径影响人体健康。建筑材料中的⁴⁰K、铀系子体（如²¹⁴Pb和²¹⁴Bi）以及钍系子体（如²²⁸Ac和²⁰⁸Tl）是室内γ辐射的主要来源，可对室内人员造成持续全身照射。

在将建筑材料放射性源转化为实际人体受照剂量的评估体系中，早期研究多依赖特定的标准几何假设和经验公式。已有研究表明，恰当考虑建筑材料几何形状会强烈影响室内剂量率。欧洲委员会辐射防护112号报告（RP112）提出了标准房间模型（尺寸5 m×4 m×2.8 m，墙厚20 cm）和活度浓度筛选指数，但此类简单指标只能提供粗略的宏观安全估计，难以准确表征真实居住环境中辐射的空间异质性。在解析计算方法方面，Markkanen基于点核法和γ射线积累因子开发了解析计算模型，但传统Markkanen模型在处理复杂空间几何拓扑和非均匀材料方面存在局限。后续学者从几何适应性、参数优化和能谱扩展等方面进行了改进，但目前主流模型多为Markkanen模型的改进型，解析计算方法相对复杂。虽然蒙特卡罗（Monte Carlo）模拟在建模真实建筑时受限于几何简化、材料均质化和计算成本等问题，但MCNP和OpenMC等软件发展迅速、功能强大，因此有必要利用Monte Carlo方法建立相关模型，对放射性剂量分布进行快速准确的模拟研究。

二、研究目的与意义

本研究聚焦于室内空气吸收剂量率，将其作为人员剂量评估的物理基础。旨在建立基于Monte Carlo方法的室内γ剂量率计算模型，揭示标准房间内的三维剂量场分布规律，并系统分析房间尺寸、高度以及门窗开口等建筑几何参数对室内γ剂量率的敏感性影响，从而为辐射防护和房间布局优化提供定量依据。

三、关键技术方法

研究人员采用Monte Carlo方法，使用MCNP6软件建立标准房间（500 cm×400 cm×280 cm，墙厚20 cm）的几何模型，墙体设为密度2350 kg/m³的均质混凝土，室内空气密度取1.205 kg/m³。源项定义中假设放射性核素均匀分布在墙体、天花板和地板中，并以²²⁶Ra和²³²Th代表铀系和钍系衰变链，其衰变子体处于长期平衡状态。剂量计算选用MCNP的F5点探测器获取光子注量，结合国际放射防护委员会（ICRP）第74号出版物的剂量转换系数，将光子注量转换为吸收剂量。为便于标准化比较，引入了剂量率转换因子（Dose rate conversion factor, DCF），单位为(nGy·h^-1)/(Bq·kg^-1)。所有模拟任务输运粒子数设为10⁶~10⁸，确保计算结果统计误差不超过1%。

四、研究结果

1. 标准房间模型中剂量率梯度分布的计算

研究人员首先考察了固定标准房间模型内的空间剂量分布。结果显示，各方向上的整体剂量场均呈现中心低、边缘高的非均匀分布模式。在距房间中心X方向1.5 m和Y方向1 m范围内，剂量率波动较小；当测量点靠近周围墙壁时，剂量率呈非线性加速上升趋势。Z方向因房间高度小于长宽尺寸，垂直方向空间梯度变化更为紧凑。不同核素的空间梯度存在差异，这由其释放γ射线的能谱特性决定。⁴⁰K发射高能单能γ光子（1460.8 keV），²³²Th衰变系列的主要剂量贡献来自其子体²⁰⁸Tl发射的高能光子（2614.5 keV），这些高能光子在建筑材料中平均自由程长、穿透能力强，当测量点靠近源平面时，高能光子直射增强效应尤为敏感，导致更高的峰值相对偏差。而²²⁶Ra衰变系列发射的γ射线主要集中在中低能区，在建筑材料中衰减截面大、自吸收严重，因此靠近墙壁时低能射线的通量增幅小于高能射线。三维网络切片显示，封闭房间内空气吸收剂量率范围为59.4 nGy/h至61.2 nGy/h，靠近角落和边缘区域因同时暴露于多个相互正交的面源，近场内多个高强度辐射源的空间矢量叠加使边缘区域剂量率达到局部峰值。

2. 房间尺寸的敏感性分析

研究人员在固定净空高度H=2.8 m条件下，构建了一系列方形平面房间模型，水平边长从1 m增至10 m（对应体积2.8 m³至280 m³）。以标准房间体积56 m³为基准，计算了各尺寸模型几何中心剂量率的百分比偏差。结果显示，三种核素的剂量率相对偏差曲线随房间体积呈先降后升的非单调变化特征，这是辐射输运中的几何衰减效应与源强增加效应动态博弈的结果。在小尺度空间中，侧墙离中心探测点很近（仅0.5 m），短辐射距离导致中心接收较高辐射通量；随着房间平面尺寸初步扩大，侧墙向外推离，射线发生显著几何发散，侧墙对中心点的剂量贡献迅速下降，此时尽管天花板和地板体积同步增大，但新增放射源位于远离中心的边缘区域，其低效剂量贡献不足以弥补侧墙后退造成的损失，宏观上表现为总剂量率快速下降。随着平面尺寸继续扩大，侧墙后退减弱了几何衰减效应的影响，在固定墙厚下，房间平面扩张直接导致天花板和地板面积扩大及四面墙总长度增加，即围护结构总体积和质量增加，放射性核素总活度增大，增加的放射源强度向室内发射更多光子，从而使总剂量率回升。

3. 房间高度的敏感性分析

研究人员保持建筑面积5 m×4 m不变，将房间高度H从2.0 m逐步增加至4.0 m（步长0.2 m）。以标准房间模型H=2.8 m为基准，观察点设于几何中心(0,0,0)。结果表明，对于给定建筑面积，三种核素产生的吸收剂量率相对偏差均与房间高度呈显著负相关趋势。在2.0 m低矮空间中，⁴⁰K、²³²Th和²²⁶Ra贡献的剂量率分别比标准2.8 m层高增加约0.86%、0.62%和0.26%；而在4.0 m高度空间中，三者剂量率分别下降约0.78%、0.58%和0.51%。驱动这一负相关趋势的物理机制除源体积变化外，还包括空间立体角再分配和辐射源自吸收的协同效应。虽然房间高度H增加导致墙体体积增大、核素总活度上升，但这种几何拉伸同时引发了平均辐射距离的增加和斜入射引发的墙体自吸收加剧两方面衰减效应。侧墙新增立体角主要分布在天花板和地板附近的远端边缘区域，这些远端侧墙发出的γ射线到达房间中心时出射倾角较大，需穿透更长的等效物理厚度才能逸出墙体进入空气，产生一定的自吸收衰减。侧墙新增立体角的辐射贡献效率和源强增加的正效应弱于斜入射引发墙体自吸收和平均辐射距离增加带来的负效应，因此整体呈递减规律。

4. 房间开口的敏感性分析

研究人员在标准房间模型的不同墙面和天花板上设置了不同尺寸和位置的开口以模拟窗户，并计算其对室内中心参考点吸收剂量率的相对影响。结果显示，窗户的存在显著降低了室内空气吸收剂量率，降低程度与窗户面积成正比。这是因为窗口面积本身不再是辐射源，直接减少了该墙的总活度；同时窗口为墙体γ射线提供了直接逃逸室外的通道，减少了其对中心点的贡献。窗户位置对剂量率降低效果有显著影响，将窗户置于墙面中心比置于角落更能有效降低室内中心点剂量率，因为房间中心探测器对其正对且最近的墙面中心区域最敏感，在墙面中心开窗占据了最大的立体角，恰好移除了对其贡献最大的辐射源。剂量率降低效果还与开窗所在墙面有关，在天花板开窗对中心点²²⁶Ra剂量率的降低效果最强（3.36%），其次是长侧墙（1.95%），最后是短侧墙（1.42%），这同样由几何因素决定。

基于真实住宅环境，研究人员根据《住宅设计规范》（GB 50096-2011）构建了包含门（0.9 m×2 m）和窗（2.0 m×1.5 m）开口的三维物理模型。对比理想封闭模型与非对称门窗模型发现，门窗开口区域直接切断了相应区域的墙体源项贡献，在包含门窗的侧墙区域，三维热图中出现了穿透墙体的深蓝色低剂量区，窗口正前方和门口延伸区域的局部最小剂量率可显著降至35 nGy/h以下。这既体现了源体积减少导致的源强削弱效应，也直观反映了开口对室内射线的逃逸效应。其次，由于缺失立体角导致空间梯度畸变，在接近带窗墙面时，截面中心明显出现与窗户尺寸对应的矩形低剂量区域，由此产生的辐射冷点沿垂直于窗面的轴向室内延伸，与门口形成低剂量梯度通道。门的偏心非对称设置进一步增强了室内剂量场的非均匀性，门口区域的低剂量冷点从地面延伸至2 m高度并偏向房间一侧，而另一侧的完整实心墙维持着较高的辐射贡献，这种非对称梯度分布导致室内极端剂量点在一定程度上偏离几何中心。

五、总结与结论

传统基于无窗标准房间的计算模型，由于将三维辐射场简化为单一几何点并忽略了门窗开口造成的立体角损失和辐射逃逸效应，往往倾向于高估真实环境中的室内剂量水平。本研究建立的三维剂量评估模型能够更准确地表征空间异质性，为室内高占用区域的安全布局提供了定量依据。研究结论如下：（1）揭示了室内辐射场的向心衰减与近场叠加的空间分布模式。在标准封闭房间中，墙体发出的γ射线受空间传输距离和空气介质衰减影响，室内剂量场呈现中心低、边缘高的梯度特征；在靠近角落和边缘的区域，由于多个正交辐射源的近场叠加，剂量率达到全局峰值。（2）阐明了建筑面积和房间高度对剂量率的几何重塑机制。对于建筑面积，剂量率先因几何衰减而下降；随着源体积扩大，后期源强增强导致剂量率逐渐上升。在固定建筑面积下，中心剂量率随房间高度单调递减，因为平均辐射距离增加和斜入射自吸收的负效应大于源体积增加的正效应。（3）门窗结构作为辐射逃逸通道，对室内剂量分布产生局部重塑。开门窗移除了相应体积的辐射源并为室内射线提供了直接逃逸路径，从而破坏了封闭剂量场的对称性，在门窗开口延伸区域形成低剂量区，并使极端剂量点偏离几何中心。本研究的发现可为辐射防护和房间布局优化提供一定参考。