采用有限差分代码MODFLOW SURFACT(Hydrogeologic, Inc., 1996)对概念模型进行了模拟,该代码与原始的MODFLOW(McDonald和Harbaugh, 1988)具有相似的输入要求。
本节描述了用于模拟概念模型的MODFLOW-SURFACT和MODFLOW文件所需的输入参数。
在开发和模拟ETA排水场槽模型时,使用了MODFLOW-SURFACT的BASIC、BCF4、Recharge、Drain、ATO4、RCH4、PCG4和NRB1包。
输入文件是使用Groundwater Vistas版本3.45的图形用户界面(GUI)预处理器创建的(Rumbaugh等,2002)。
3.3.1 Basic Package
Basic(基础)包是MODFLOW代码读取的任何模型的第一个包。Basic包包含有关模型几何、选择的模块、时间离散化和初始压力头数据的信息。
3.3.1.1 Model Geometry
模型几何包括有关模型中的行数、列数和层数的信息。对于非饱和区域模型,行和列的宽度应尽可能小,以提高结果的精度和准确性,但紧密的空间离散化会导致大型数据文件,从而减慢模拟速度。在考虑了0.5英尺、1英尺和2英尺的单元宽度后,决定使用1英尺的单元宽度是模型的最佳选择。因此,选择了一个60行60列的网格用于模型。
模型中顶层的层厚度范围从0.335英尺到1英尺,而底层的层厚度范围从0.5英尺到3英尺。这是因为模型表面附近的层在模拟期间由于模型顶表面的充水和蒸散发边界条件而受到了显著的湿度运动的影响。显著的湿度运动导致了收敛困难和