李林子-GMS在承钢双大尾矿库氟污染模型中的应用.pdf

GMS 在承钢双大尾矿库氟污染模型中的应用 李林子 南京大学水科学系 摘要：本文以承钢双大尾矿库区地下水为研究对象，在系统阐述了GMS6.0 中构建溶质运移 模型各个模块（MODFLOW、MT3DMS 和反求参数模块）的基本原理和一般方法的基础上， 在GMS 中建立研究区的地下水流模型，并对水流模型进行识别与检验，继而建立了地下水 的氟污染运移模型。模拟计算结果表明，研究区地下水所赋存的第四系松散层对污染物代谢 能力很强，若维持现状，5 年以后整个研究区地下水中的氟离子含量将基本满足国家生活饮 用水卫生标准。 关键词：GMS 氟污染 水流模型 模型识别 溶质运移模型 引言 随着经济的迅速发展和人民生活水平的提高，一方面，地下水污染问题日益 加剧，另一方面，人们对关乎自身健康的地下水水质的要求日益提高。2007 年， 承德市白庙子村村民和承德钢铁集团公司就针对地下水氟污染的侵权赔偿发生 了法律纠纷，河北省高级人民法院委托山东正元建设工程有限责任公司对此进行 相应的司法鉴定，笔者也于2007 年10 月和12 月随导师参与了该地区地下水污 染问题的调查取证和野外试验工作。 在地下水污染研究的实际应用中，由于受限于人力、财力资源，通常只能获 得有限的地下水水位和污染物浓度观测资料。于是，利用数值模型对地下水流和 溶质运移问题进行模拟的方法以其有效性、灵活性和相对廉价性逐渐成为地下水 [1] 研究领域中一种不可或缺的重要方法,并受到越来越大的重视和广泛的应用 。近 20年来，依赖于人机交互、计算机图形学和科学可视化等计算机技术的迅猛发展， [2] 具有可视化功能的地下水模拟软件迅速发展。 目前，国际上最具影响力的几个 地下水模拟软件有MODFLOW、MT3D、PEST、FEFLOW、Visual Groundwater 、GMS 等，其中，GMS 以其良好的使用界面,强大的前处理、后处理功能及优良的三维 [1,3] 可视效果已成为当今功能最完善的地下水模拟软件 。 本文在承德市白庙子村地区水文地质调查和野外试验资料的基础上，利用地 下水模拟软件GMS6.0 以该区地下水的氟污染运移为对象进行数值模拟研究，建 立了研究区的地下水水流模型和溶质运移模型，为预测其氟含量的变化范围和变 化梯次提供依据。 1 研究区概况 1.1 自然地理概况 研究区位于承德市西郊双滦区，地处冀北山地、燕山山脉中段。白庙子村坐 落于双滦区双塔山镇南部的牦牛河与滦河交汇的滦河东岸，北靠原承钢集团双大 矿业尾矿库，西北为啤酒厂和输送机厂，东为砂砾岩低山区。研究区系暖温带大 陆季风性气候，冬季长寒冷干燥，夏季短炎热多雨。区内昼夜温差大，最高气温 41.5 ℃，最低气温-23.3℃，多年平均气温8.8 ℃。降水量年内分配不均，7-8 月为 雨季即洪水期，历年24 小时最大降水量达151mm，6 至8 月的降水量约占全年 降水量的70% 。据资料统计，1960 年至2006 年，多年平均降水量为519.25mm ； 1970 年至2006 年多年平均蒸发量为1534.9mm。区内主要河流为滦河与牦牛河。 滦河发源于丰宁大滩界牌梁，呈蛇曲形东西流向，为常年性河流，从西部流入本 区，河谷宽400-2000m （白庙子段宽1500m 左右），多年平均径流量7.19 亿立方 米。牦牛河发源于大庙镇北梁村，向南流经本区的主要村镇，在双塔山镇南注入 滦河。牦牛河为季节性河流，河谷较窄，宽约 200-400m ；该河水量小，干旱季 节经常干涸无水，但每当夏季大雨时节，河水暴涨，汇入滦河。 1.2 水文地质条件 研究区区域位于中朝准地台（Ⅰ2 ）的Ⅱ级构造单元——燕山台褶带（Ⅱ2 ） 北部，北属承德拱断束（Ⅲ26 ）。区内褶皱运动微弱，构造环境相对稳定。主要 构造为近东西向构造，其形成时间较早，而后被北东向、北西向断裂切割。研究 区主要出露岩层为侏罗系后域组一段（H1j2）紫红色砾岩、砂砾岩和第四系松散 层。紫红色砾岩、砂砾