研究放射性核素在地下水环境中传输的动力学行为是评价核废料处置库安全性能的一个关键性问题。以双重介质理论为基础,充分考虑了地下水变密度对核素浓度分布的影响,建立了孔隙-裂隙介质中核素随地下水迁移的耦合数值模型,其中相对密度项使地下水流控制方程具有非线性性,为此将Galerkin有限元法和Picard迭代法相结合对耦合模型进行了离散,并模拟了核素泄漏对地下水污染的动态过程。数值仿真的结果表明:污染面积随时间推移不断扩大,且纵向扩散显著;数值仿真结果没有出现弥散现象,验证了数值算法的可靠性;核素浓度的分布对水密度变化表现较为敏感,且不同核素(226Ra、210Pb和90Sr)由于各自放射性衰变特征的巨大差异,受到水密度影响的程度也有很大差别。因此,在对裂隙介质核素迁移的数值仿真分析过程中,不能忽略地下水密度的变化对其造成的影响。