汾河二库蓄水对兰村水源地影响的数值模拟研究(2)

　　2.2 模型

　　根据模拟区地下水系统水文概念模型，建立了计算区内水位分布数学模型。潜水和中深层承压水系统地下水运动的数学模型分别为:

　　潜水系统：

　　                   （1）

　　承压水系统：

　　                        （2）

　　式中: ——浅层地下水子系统的平均渗透系数(LT^-1)；

　　 ——中深层地下水子系统导水系数(L²T^-1)；

　　 ——潜水系统水位标高（L）；

　　 ——中深层地下水子系统水位标高（L）；

　　 ——潜水含水层的给水度（无量纲）；

　　 ——承压含水层的贮水系数（无量纲）；

　　 ——潜水系统源汇项（LT^-1）；

　　 ——t时间时（x，y，z）边界流量（LT^-1）；

　　 ——弱透水层系统顶部边界标高（L）；

　　D——浅层地下水子系统和中深层地下水子系统同时确定的平面渗流域（L²）；

　　 ——D的已知流量边界；

　　 ——二类边界外法线向量。

　　2.3 模型剖分

　　研究区的总面积为162km²，在平面上将其剖分为209行，184列，在垂向上分为两层，其网格为矩形，共计76912个单元格。其中，有效单元格45465个研究区的总面积为162km²，在平面上将其剖分为209行，184列，在垂向上分为两层，其网格为矩形，共计76912个单元格。其中，有效单元格45465个，无效单元格31447个，每个单元格的平面面积约为0.007km²。