沙坑回灌补给地下水的数值模拟分析

2014-03-22田竹君肖长来田浩然

东北水利水电 2014年1期

关键词：沙坑渗透系数潜水

杜超，田竹君，肖长来，田浩然

(1.松辽流域水资源保护局松辽水环境科学研究所，吉林长春 130021；2.吉林大学环境与资源学院，吉林长春 130021)

沙坑回灌补给地下水的数值模拟分析

杜超1，田竹君1，肖长来2，田浩然1

(1.松辽流域水资源保护局松辽水环境科学研究所，吉林长春 130021；2.吉林大学环境与资源学院，吉林长春 130021)

文章对一般性的沙坑回灌补给地下水模型进行了概化，建立了水文地质概念模型和数值模型，采用 GMS 软件对沙坑回灌补给地下水的效果进行了模拟分析。结果表明，沙坑回灌补给地下水的累计渗漏量与时间呈线性关系；渗漏速率在渗漏初期大于渗漏后期，经过一定时间后趋于稳定；单个沙坑回灌补给引起地下水明显上升的范围较小，效果有限。

沙坑回灌；地下水；渗漏量；渗漏速率；数值模拟

由于社会经济的不断发展，在中国的很多地区尤其是北方地区，地表水资源已经远远不能满足人类的生活生产需要，许多城市供水与农业灌溉用水走上了盲目大规模开采地下水这条道路，造成了地下水资源不断减少，地下水位持续下降等不良后果。随着观念的转变，人们对地下水资源的保护越来越重视，很多地方都开始着手建立地下水回灌工程。在实施工程前，应对回灌可能产生的效果进行可行性分析，而数值模拟技术正是进行分析的有效手段。下面采用数值模拟技术对沙坑回灌补给地下水的效果进行分析，讨论回灌过程中沙坑渗漏量、渗漏率、影响半径随时间的变化关系，得出了沙坑回灌补给地下水的一般性规律。

1 沙坑回灌模型

为了探讨沙坑回灌补给地下水的一般性规律，对沙坑回灌模型做如图 1 的概化，矩形沙坑深 10 m，BC=ED=10 m，CD=20 m。沙坑底水平，沙坑水位为 3 m，未渗漏前地下水初始水位距沙坑底部 2.5m。将土层分为两种方案来分别模拟：一种是土层为层状土，上层为粗砂，下层为砂砾石；另一种是土层为均质砂砾石。对渗透系数取不同的值来模拟不同时刻沙坑内水流的渗流情况。由于是做一般性探讨，所选参数相对较大。

方案一：上部为粗砂，渗透系数 K=50 m/d；下部为砂砾石层，渗透系数 K=150m/d，给水度 μ=0.3。

方案二：土层均质砂砾卵石渗透系数 K=200 m/d，给水度 μ=0.3。

图1 沙坑回灌模型示意图

2 水文地质概念模型

2.1 含水层概化

对方案一，含水层渗透性在垂向上呈层状分布，概化为非均质各向同性潜水含水层；对方案二，含水层渗透性在垂向上不发生变化，概化为均质各向同性潜水含水层。

2.2 水力特征概化

沙坑回灌模型中地下水流各要素随时间发生变化，概化为非稳定流；地下水流主要在平面方向上流动，潜水面随时间发生动态变化，概化为潜水二维流。

2.3 边界条件概化

1）侧向边界。由于沙坑中水流不断补给地下水，将很快影响到距沙坑较近的地方。若边界条件选择离沙坑过近，则边界的水头及流量会很快发生变化，从而影响模型模拟结果。因此，边界采用远离沙坑中心的零流量边界，边界距沙坑中心 5 km。

2）垂向边界。潜水面为潜水含水层的上部水量交换边界，潜水含水层底板为隔水边界。当沙坑中的水流与地下水取得水力联系后，由于沙坑中水头固定不变，此时即为潜水含水层的定水头边界。

3 数学模型的建立与求解

3.1 数学模型的建立

均质、潜水二维地下水非稳定流数学模型如下：

式中：K——潜水含水层渗透系数，m/d；μ——潜水含水层给水度；H——潜水水位，m；B——潜水含水层底板标高，m；Qi——源汇项补给/排泄强度，m/d；H0——初始潜水水位，m；H1——研究区一类边界处潜水水位，m；q——研究区二类边界单宽流量，m3/(d·m)；x，y——坐标，m；D——计算区范围；Г1，Г2——研究区一类、二类边界。

3.2 模型离散

1）空间离散。以沙坑定水头矩形边界 4 个顶点为基准点，对沙坑定水头边界附近水力梯度变化较大的区域进行网格加密。将计算区剖分为 174 列、174 行，共 30 276 个网格单元。最小网格单元长 0.25 m，宽 0.25 m，面积 0.625 m2。网格偏离系数取 1.1。

2）时间离散。模拟期为 3 年，以观测长时间条件下沙坑回灌补给地下水的稳定情况。模拟初期地下水水位变化较快，变幅较大，故模型运行的前 5 d以1 d为一个应力期，1 h 为一个时间步长；第 6~30 d 以 1 d 为一个应力期，4 h 为一个时间步长；第 31~365 d 以 1 d 为一个时间步长；第 366~1 095 d 以 5 d 为一个时间步长。

3.3 初始水位

数学模型在沙坑水与地下水取得水力联系并使坑底地层达到饱和后才适用。显然沙坑开始注水后会迅速渗透至潜水面并达到饱和。在计算时,假设沙坑中水位由无水状态突然上升到设计水位，坑底地层达到饱和，给定边值条件，作为稳定流问题可以求得它的解 H x，（）y ，把这个解近似作为沙坑充水后非稳定流问题的初值 H0x，（）y。