本文主要研究了原子能院两侧宽5km、下游10km的区域上包气带和潜水层中的人工放射性核素~3H迁移的浓度的时空分布。使用美国国家橡树岭实验室开发的FEMWATER／FEMWASTE程序对模型求解，并根据该区现场的试验资料或与该区介质特性相似的文献资料来选取程序的计算参数。FEMWATER的水流方程主要是Richards方程式及其改进形式；FEMWASTE的溶质运移机制有对流、水力学弥散和分子扩散、放射性衰变、吸附、源汇、生物降解，在计算时加进区域质量平衡，还可处理含水层多层之间的渗透问题。该程序采用迦辽金有限元法求解，在时间、空间、初始条件、源汇、边界的处理上具有很大的灵活性。在对研究区进行矩形剖分的基础上，分别采用一维柱渗透模型和二维水平模型来模拟~3H穿透包气带和在潜水层中的迁移过程，并根据该区域的水文地质条件、水位和浓度的监测资料等来设置模型的边界和初始条件。假设~3H以某个浓度进入，分别计算了包气带和潜水层中各个时间点和各个节点的浓度值。然后将包气带的计算结果绘制成浓度穿透曲线和沿程分布曲线，将潜水层的计算结果绘制成区域的相对浓度等值线图。在这些分布图的基础上，对~3H在包气带和潜水层的迁移行为进行了分析。模拟结果表明，该区内包气带~3H随水下渗的速度很小，为0.001 m／d，流速分布均匀，仅有一小段层面上有溶质的上下运动，方向基本一致。潜水层内流场基本由北偏南流动，流速为1～10m／d，等水位线由西北至东南等间隔均匀下降。~3H相对浓度分布图的分析表明，~3H在包气带中的迁移分为三个阶段，第一阶段是穿透阶段，浓度随距地面的深度增加而减小，在计算开始后的第840天穿透包气带底层；第二阶段是均衡阶段，~3H的浓度分布趋于均衡，不同深度上的浓度值相差变小；第三阶段~3H在底部沉积，这时底部的零通量边界不再适用。~3H在潜水层中随着水流向下游的迁移，分别模拟了点源、线源、面源三种不同的情况。结果都显示相对浓度在纵向上变化明显，横向上变化不大，这是由于纵向弥散系数比横向弥散系数大得多。并且从模拟的区域上看，~3H的迁移基本限于大石河的东侧，出现明显的~3H浓度的狭长带。当初始浓度相同时，源的形状不同，引起的迁移的快慢也不同，范围也有所差异。