核电厂发生事故的概率极低,但是一旦发生事故,释放出的放射性物质可能对附近的公众和安全工作人员造成严重的辐射损伤,因此需要及时评估放射性核素在大气中的扩散情况,以估计放射性物质对附近人员的照射剂量。利用估计的辐射剂量值可以评估辐射对电厂附近生物的损害以及潜在危险,从而制定合理的公众紧急疏散计划及其他应急计划。本文使用MATLAB程序开发了基于高斯烟羽模型的大气扩散模型计算代码,用于估计放射性核素空气浓度、总有效剂量当量(TEDE)、地面沉积浓度和地面照射剂量率。代码包含多个输入参数/变量,如放射性核素混合物及其活度、气象条件、释放高度和接收器高度、地表类型、大气稳定性等级、地表粗糙度、逆温层高度、烟羽的浮力和动量抬升、呼吸速率、烟羽和环境温度、烟羽可吸入部分和不可吸入部分的沉积速度。以K2/K3核电厂附近的地形和气象信息作为输入数据,源项参考福岛核电站事故相关数据。代码综合考虑了放射性核素在大气中的运输、沉积和扩散,计算了放射性活度时间积分空气浓度、地面沉积浓度、地面照射剂量率、浸没剂量值、吸入剂量值、地面照射剂量值以及总有效剂量当量。此外,代码还可输出人体器官剂量信息以及放射性混合物中每种核素导致的甲状腺和骨髓剂量的信息。结果表明,需要采取就地避难及分发碘片(ITB)等手段以减轻1-131对甲状腺造成的辐射损伤,而在放射源释放点下风向1.4km外骨髓的剂量在IAEA规定的安全标准1Sv内。代码的所有结果都被导出到了 Microsoft Excel文件中,据此可绘制随下风向距离变化的剂量曲线。将代码计算结果与已得到广泛使用的标准参考代码HotSpot进行比较,结果表明代码结果与HotSpot相近,根据代码结果绘制的曲线与HotSpot得到的曲线一致。这表明,代码有实用价值并可用于进一步的类似研究。误差幅度小于2%。