尉鹏翔

(山西省水文水资源勘测局太原分局,山西太原 030002)

近年来,人们普遍注意到由于超量开采地下水和水质污染,地下水生态环境受到破坏,且日趋严重。应用科学有效的计算方法寻找处理地下水污染问题的方法,对于地下水的污染防治有着极其重要的作用[1]。目前地下水系统数值模拟方法主要有有限差分法(FDM)、有限单元法(FEM)、边界元法(BEM)和有限分析法(FAM)等,可以用来研究海水入侵、地面沉降、地下水位降低、地下水过量开采和地下水污染等问题。随着计算机技术的快速发展,数值计算方法在地下水资源分析评价中得到广泛应用。笔者拟借助Visual Modflow商业化软件研究北京某污染场区地下水溶质运移问题。

1 研究区概况

所研究污染场区位于北京市区东南部,地貌上属于永定河冲洪积扇的中下部,场区地形基本平坦,地面以下30m深度范围内一般存在2个含水层,局部区域有上层滞水,单层含水层厚度一般小于10m,仅部分地区大于10m。本文只研究埋深在10～20m的含水层和埋深在25m左右的承压含水层以及它们之间的黏土弱透水层,鉴于此,模型在空间上概化为3层:潜水含水层、弱透水层和承压含水层。

场区目前潜水流向由北向南,北部边界是侧向补给边界,主要接受侧向补给,南部边界为整个含水层系统的微量流出边界,由于地下水整体流向由北向南,故模拟区东、西两侧可以流线为边界,将它们设为隔水边界。承压含水层水流向由西北向东南,西北部为侧向补给边界,主要接受侧向补给,东南部为侧向流出边界,含水层底部大多为透水性较差的黏土地层,对系统内部影响很小,定义为隔水边界。……