放射性核素迁移的包气带-饱和带耦合模拟研究

摘要/Abstract

摘要： 为了使地下水环评结果更加真实,通过试验和拟三维方法数值模拟开展包气带-饱和带耦合模拟工作。为进一步说明包气带对核素迁移的迟滞(延迟和滞留)作用,开展了不考虑包气带的模拟并与耦合模拟结果对比。结果表明:(1)由于包气带对核素迁移的迟滞作用,吸附性较强的Se-79、Pd-107、Cs-135和Sn-126迁移至潜水面处的浓度远小于泄漏的放射性液体核素浓度,所致公众食入年待积有效剂量已小于剂量限值0.1 mSv,只有吸附性较弱的C-14所致公众食入年待积有效剂量大于0.1 mSv,需进行耦合模拟;(2)C-14在对流作用下不断向下游迁移,在弥散作用下污染范围逐渐增大,在地下水稀释及衰变、吸附作用下浓度逐渐降低。通过耦合模拟得到的地下水中放射性核素所致公众食入年待积有效剂量小于剂量限值0.1 mSv;(3)包气带对核素迁移的迟滞作用使地下水中C-14的浓度大大降低、迁移距离明显变小:地下水中C-14峰值浓度在考虑包气带时为2.15 Bq/L,不考虑包气带时为9.4×10⁴ Bq/L;事故发生后第1、5和10年,考虑包气带时C-14迁移距离分别为1.5 m、4 m和4.2 m,不考虑包气带时分别为3 m、12 m和17 m。

关键词: 拟三维, 放射性核素, 包气带, 地下水, 耦合模拟

Abstract: In order to obtain more real environmental impact assessment results for groundwater,a simulation study on vadose-saturation zone coupling was carried out by using quasi-three-dimensional method based on experiment and numerical simulation. To further explain the hysteresis effect of vadose medium on nuclides migration, simulation without considering the vadose medium was carried out and the results was compared with that of the coupling simulation. The results show that: (1) For strong adsorbed Se-79, Pd-107, Cs-135 and Sn-126, their concentration arriving at water table is far less than the leaked liquid due to the hysteresis effect. And their annual committed effective dose for public will be less than the dose limit of 0.1 mSv. For weakly adsorbed C-14, the annual committed effective dose for public will be more than 0.1 mSv and it requires coupling simulation. (2) C-14 migrates to downstream by convection, and pollution range increases through dispersion. At the same time the C-14 concentration decreases by dilution, decay and adsorption, therefore the annual committed effective dose for public will be less than 0.1 mSv. (3) The concentration of C-14 in groundwater decreases greatly and the migration distance decreases obviously because of the hysteresis effect of vadose medium on nuclides migration. The peak concentration of C-14 in groundwater is 2.15 Bq/L taking vadose zone into consideration; while leaving out the vadose zone, the peak concentration of C-14 will be 9.4×10⁴ Bq/L.The migration distance of C-14 is 1.5 m, 4 m and 4.2 m respectively in 1, 5 and 10 years after the accident, taking vadose zone into consideration; while leaving out the vadose zone, the migration distance of C-14 will be 3 m, 12 m and 17 m respectively.

Key words: quasi-three-dimensional, radionuclide, vadose zone, groundwater, coupling simulation

中图分类号:

P641

张雪, 胡继华, 施晓文, 王小元. 放射性核素迁移的包气带-饱和带耦合模拟研究[J]. 辐射防护, 2021, 41(6): 514-522.

ZHANG Xue, HU Jihua, SHI Xiaowen, WANG Xiaoyuan. Study on vadose-saturation zone coupling simulation of radionuclide migration[J]. RADIATION PROTECTION, 2021, 41(6): 514-522.

参考文献

[1] 李书绅, 王志明, 郭择德,等.核素在非饱和黄土中迁移研究[M].北京:原子能出版社,2003.
LI Shushen, WANG Zhiming, GUO Zede, et al. Study on the migration of nuclides in unsaturated loess[M]. Beijing:Atomic Energy Press, 2003.
[2] Havard P L,Prasher S O, Bonnell R B,et al. LINKFLOW, a water flow computer model for water table management: Parti. Model development [J]. Transaction of the ASAE,1995,38(2):481-488.
[3] Facchi A,Ortuani B,Maggi D,et al. Coupled SVAT-groundwater model for water resources simulation in irrigated alluvial plains [J]. Environmental Modelling & Software, 2004, 19(11): 1053-1063.
[4] Niswonger R G, Prudic D E. Modeling variably saturated flow using kinematic waves in MODFLOW[M]//Hogan J F, Phillips F M, Scanlon B R,eds.Groundwater Recharge in a Desert Environment-The Southwestern United States: Volume 9, Water Science and Application Series. American Geophysical Union, 2004: 101-112.
[5] Niswonger R G, Prudic D E,Regan R S. Documentation of the Unsaturated-Zone Flow (UZF1) package for modeling unsaturated flow between the land surface and the water table with MODFLOW-2005 [R]. U.S. Geological Survey Techniques and Methods, 2006:6-A19-62.
[6] Twarakavi N K C, Simunek J, Sophia S. Evaluating interactions between groundwater and vadose zone using the HYDRUS-based flow package for MODFLOW [J]. Vadose Zone J, 2008,7(2): 757-768.
[7] Yakirevich A,Borisov V, Sorek S. A quasi three-dimensional model for flow and transport in unsaturated and saturated zones: 1. Implementation of the quasi two-dimensional case [J]. Adv Water Resour, 1998,21(8): 679-689.
[8] Kuznetsov M,Yakirevich A, Pachepsky Y A, et al. Quasi 3D modeling of water flow in vadose zone and groundwater [J]. J Hydrol, 2012, 450-451: 140-149.
[9] 林琳, 杨金忠, 史良胜,等.区域饱和-非饱和地下水流运动数值模拟[J].武汉大学学报, 2005,38(6):53-57,71.
LIN Lin, YANG Jinzhong, SHI Liangsheng, et al. Saturated and unsaturated groundwater flow numerical simulation in large scale zone[J].Engineering Journal of Wuhan University, 2005,38(6):53-57,71.
[10] 朱焱.区域拟三维饱和-非饱和水流与溶质运移模型研究与应用[D].武汉大学,2013.
ZHU Yan. Study on Quasi-3D water flow and solute transport model in regional scales and its application[D]. Wuhan University, 2013.
[11] 查元源.饱和-非饱和水流运动高效数值算法研究及应用[D].武汉大学,2014.
ZHA Yuanyuan. Research on cost-effective algorithm for unsaturated-saturated flow and its application[D]. Wuhan University,2014.
[12] 于洋,韩鹏,林谧.半干旱地区包气带-饱和带水盐耦合模拟[J].干旱区研究,2017,(4):741-747.
YU Yang, HAN Peng, LIN Mi. Water-salt coupled simulation in vadose zone and saturated zone in semiarid area[J]. Arid Zone Research,2017,(4):741-747.
[13] 张雅楠.基于HYDRUS-MODFLOW地表水与地下水耦合模型研究[D].郑州大学,2019.
ZHANG Yanan. Study on coupling model of surface water and groundwater based on HYDRUS-MODFLOW[D]. Zhengzhou University,2019.
[14] EPRI. Groundwater monitoring guidance for nuclear power plants[R]. 2006.