某生活垃圾填埋场地下水污染物运移模拟_宋旌鸣.pdf

引 用 格 式：，（）：宋旌鸣，黄勇某生活垃圾填埋场地下水污染物运移模拟 甘肃科学学报，（）：某生活垃圾填埋场地下水污染物运移模拟宋旌鸣，黄勇（河海大学地球科学与工程学院，江苏 南京 ）摘要以南京溧水某生活垃圾填埋场作为研究对象，详细调查水文地质条件，概化水文地质概念模型，利用 软件对该场地进行建模和拟合调参，选取特征污染物氨氮作为模拟对象，对地下水污染物的运移规律进行了模拟。模型考虑防渗有效与防渗失效两种工况下污染物的运移规律，模拟污染物在 天、天、年、年和 年的迁移规律、污染范围及浓度情况。结果表明：在防渗有效时，污染物 年内运移的最大距离为 ，污染总面积为 ；防渗失效时，污染物运移 主 要 受 水 力 梯 度 影 响，污 染 物 年 内 运 移 的 最 大 距 离 为 ，污 染 总 面 积 为 。模拟结果旨在为垃圾填埋场的污染防治提供科学合理的参考依据。关键词垃圾填埋场；渗滤液；地下水污染；软件中图分类号：文献标志码：文章编号：（）地下水是我国水资源中重要的资源之一，我国水资源总量为 亿，其中地下水资源总量为 亿，占总水资源总量的。作为饮用水，我国约有 的人口把地下水资源作为主要的饮用水源，在农村超过 的人口饮用地下水，在我国北方地区的干旱半干旱城市、乡村，地下水成为唯一的饮用水来源。我国有 个城市通过开采地下水，将其作为城市的生活供水水源，城市比例达到全国城市总数的 左右。但是，随着社会的高速发展及人口、经济的快速增长，人们对环境的保护越来越不重视，导致地下水污染越来越严重。目前，地下水数值模拟软件的发展越来越快，多种数值模拟软件（如 、等）被应用于地下水的溶质运移模拟中，模拟污染物的运移规律等。研究以南京溧水区生活垃圾填埋场的环境修复和开发利用为背景，针对南京溧水生活垃圾填埋场封场后垃圾渗滤液的泄露可能造成的地下水环境污染，通过一系列试验取得水文地质参数进行地下水数值模拟，进而对研究区进行地下水环境影响预测，旨在为垃圾填埋场的污染防治提供科学合理的参考依据。研究区概况研究区位置研究区位于南京市中南部，即溧水区，面积为 ，是南京四大副城之一，地处苏皖两省、宁镇马 三 市 交 界 处，是华东地区重要的交通、物流枢纽，同时也是国家重要影视基地以及长三角地区的制造业基地。垃圾填埋场西临新戴运河，其西侧 处（距研究区 ）为马庄村，西侧 处（距研究区 ）是 街 集 村；南 面 是 美 人 山 公 墓 和农田。研究区地层概况场地位于溧水区栖凤路中段（永阳镇范围内）的 高速公路西侧，交通便利。场地地形总体呈西第 卷第期 年月 甘 肃 科 学 学 报 收稿日期：；修回日期：作者简介：宋旌鸣（），男，辽宁本溪人，硕士研究生，研究方向为水文地质。：通信作者：黄勇，：高东 低 之 势，勘 探 期 间 场 地 范 围 内 地 面 标 高 在 ，地貌单元属岗地。根据野外钻探、原位测试及室内土工试验成果，在勘察深度范围内揭示的地层自上而下为：（）（）杂填土：杂色，由黏性土、碎砖、碎石等建筑垃圾组成，结构松散，厚度在 。（）（）素填土：黄褐灰褐色，主要由黏性土组成，偶见碎砖、碎石及植物根茎，填土年限大于 年，结构松散，均匀性差，厚度在。（）（）杂填土：灰黑色，主要由生活垃圾组成，局部夹碎石、砖块等，厚度在。（）（）粘土夹卵砾石：棕红色，卵石粒径为 ，质量分数为 ，以亚圆形为主，排列不规则，局部为块石。硬塑，局部可塑，中压缩性，厚度在。（）（）强风化角砾岩：棕红色，呈碎块状，标贯击数大于 击，锤击易碎，属较软岩，岩体基本质量等级为级，厚度在。（）（）中风化角砾岩：棕红色，裂隙稍发育，取芯率一般为 ，岩体总体较完整，岩芯呈短柱状，锤击可碎，属较软岩，岩体基本质量等级级，厚度在。水文地质条件场地地下水主要属孔隙潜水和基岩裂隙水。孔隙潜水赋存于填土中，主要补给来源为大气降水、地表水，径流滞缓，排泄方式以蒸发为主，由于季节的不同对降水、蒸发存在影响，因此水位动态变化也会更加明显，地表水与地下水联系密切。勘探期间，测量得孔隙潜水稳定水位埋深在 ，稳定水位标高平均值为，初见水位一般比稳定水位高 。据区域水文地质资料表明，场地地下水位年变幅在 左右。基岩裂隙发育，赋存少量裂隙水，富水性不均。区域地下水补给来源主要有两种：一是垂直方向上的补给，来自大气降水入渗；二是侧向补给，降水量平均值为 ，雨季较长，主要集中在夏季。降水是地下水的主要补给源，地下水位与降水量关系密切，地下水位随降水量的增加而上升，随降水量的减少而下降。碎屑岩类裂隙含水层补给径流排泄条件受地形及第四系残积土厚度控制明显，在山坡基岩出露及松散覆盖层厚度较薄处直接接受大气降水补给。地下水排泄方式包括蒸发和地表径流。据气象资料显示，该地年平均蒸发量为 ，但地下水的蒸发量与地下水位埋深有关，蒸发量的大小与蒸发极限深度有关，研究中取 埋深，超过蒸发极限深度，蒸发的影响可以忽略不计，且实际地下水蒸发量比水面蒸发量要小得多。地下水的另一个排泄方式为向地表水塘、湖泊和河流排泄。模型建立水文地质概念模型水文地质概念模型首先要对地下水系统进行综合分析，在此前提下，对模拟区域的地质情况、含水层的实际边界条件以及水力特征等水文地质条件进行汇总、归纳和处理，把一个复杂的水文地质体转化为一个便于进行数学、物理模拟的概化模型。水文地质概念模型的建立需要考虑以下几个方面：模型需反应出研究区水文地质原型；模型的各类边界条件要与研究区的地下水流场特征相符合；模型中概化后的边界要尽可能与自然边界相对应；模型中人为的边界条件要从不利因素角度考虑。研究区东、西、北三侧被一大型河流所包围，因此将这三侧概化为第一类边界，称为定水头边界，河流水位即为水头；南侧位置较高，为地下水分水岭，当作隔水边界处理，潜水含水层底部为基岩，岩体完整性较好，可作为隔水边界，由此得出调查评价范围以及研究区的水文地质概念模型（见图）。图水文地质概念模型 数学模型非均质、各向异性、空间三维结构、非稳定地下水流系统的地下水水流模型为第 卷宋旌鸣等：某生活垃圾填埋场地下水污染物运移模拟（）（）（），（，）（，），（，），（，）（，），（，），?（，），（，），（）其中：为研究区模拟范围；为含水层贮水率，即当水头降低单位时所释出的水量（）；为含水层的水位（）；为源汇项，即水流模型中外来的补给或排出（）；、为个方向上的渗透系数（）；（，）为水位的分布（已知）（）；为地下水水流模型中的第一类边界；为地下水水流模型中的第二类边界；为渗流区域二类边界的外法线方向；为渗透系数张量（）；（，）为渗流区域二类边界上的一个流量函数（已知），它的正负值取决于渗流区域内流量的流入与流出，当处在隔水边界时相当于既没有流入也没有流出，即为。污染物控制方程可表示为 （）（）?，（，）（，），（，），（，）（，），（，），（，），（，），（）其中：为地下水流经介质的密度（）；为地下水所流经的某种介质的孔隙度，无量纲；为组分质量分数（）；?为溶质在介质中所被吸附的质量分数（）；为迟滞系数，无量纲；为时间（）；表征在一定流速下，多孔介质对溶质的弥散能力，即水动力弥散系数（）；为地下水渗流速度张量（）；为水流的源汇项，即水流模型中外来的补给或排出（）；为组分的质量浓度（）；（，）为浓度的分布（已知）；为溶解相一级反应速率（）；为吸附相反应速率（）；为研究区模拟范围；为给定浓度边界；（，）为上的浓度分布函数；为通量边界；f（，）为通量边界上已知的弥散通量函数。模型参数识别研究区的中心点选取为模型的坐标原点，、轴分别为正东、正北方向，轴为垂直方向，根据研究区地层资料，将垂直方向上分为层，由于研究区面积较大，因此离散也比较密集，即 个节点，个单元（见图）。为了验证模型的合理性和有效性，需要对所建的数值模型进行参数识别，即确保所建模型能够准确地反映出评价区水文地质条件概况。根据检测孔的地下水位，获得了整个模拟区的等水位线图（见图），将现场试验和勘查所获得的各水文参数作为模型参数的初始值，通过对计算水位与 个实测水位（见表）进行比较，反复调整模型参数，获得了二者最佳拟合结果。图研究区域剖分图 图评价区地下水等水位线 甘 肃 科 学 学 报 年第期表现场地下水位调查一览表 孔号纬度经度实测值计算值 拟合后模型参数：渗透系数为 ，为 ，为 ，为，为，为 ，为 ，为 ，为 ；给水度为；水力坡度为 ；孔隙度为；纵向弥散度 ；横向弥散度 ；防渗材料渗透系数为 ；防渗材料厚度为；垃圾场氨氮质量浓度为 。运行期计算工况由于垃圾填埋场渗滤液中污染物成分复杂，在此仅选择污染物超标最为严重的氨氮进行预测分析，针对研究区垃圾填埋场对水质的影响范围，分别考虑以下两种工况 ：（）假设研究区正常运行，垃圾场防渗完好，将污染物在研究区及其在周边范围内的运移考虑在内，运行 时 间 为 年，预 测 时 段 分 别 为 天、天、年、年和 年。（）突发事故条件下，垃圾场防渗完全失效，在这种情况下，污染物的运移速度变快，渗滤液流量也因此增 大，预 测 时 间 为 年，预 测 时 段 分 别 为 天、天、年、年和 年。模拟结果分析正常工况（防渗帷幕正常工作）垃圾场的氨氮初始质量浓度为 ，正常工况下第 天时，垃圾场污染源迁移距离为；第 天时，污染源迁移距离为；第年时，污染源迁移距离为；第 年时，污染源迁移距离为；第 年时，污染源迁移距离达到最大，为 ，平面上受到污染的总面积为 。由于四周存在防渗墙且东北方向有一大坝，所以在模拟时间内污染物扩散范围随时间变化缓慢，浓度升高亦缓慢（见图）。在垂直方向上，由于没有水平防渗，其影响较大，年后污染物的影响深度已经到达含水层底部（见图）。表明在正常工况下，污染物在水平方向上的迁移距离较小，而垂直方向上由于没有水平防渗影响较大。图特征污染物氨氮水平迁移图（正常工况）（）图特征污染物氨氮垂直迁移图（正常工况）（）非正常工况（防渗帷幕完全失效）当防渗墙、大坝失效后第 天时，污染源迁移距离为 ，污染范围为 ；第 天 时，污 染 源 迁 移 距 离 为 ，污 染 范 围 为第 卷宋旌鸣等：某生活垃圾填埋场地下水污染物运移模拟 ；第年时，污染源迁移距离为 ，污染范围为 ；第 年时，污染源迁移距离为，污染范围为 ；第 年时，污 染 源 迁 移 距 离 达 到 ，污 染 范 围 为 。由于防渗墙、大坝失效，污染物迁移速度明显加快，并且迁移方向受地下水流影响较大，地下水流方向从西南流向东北，所以污染物运移在东北方向较为突出，污染范围在东北方向也逐渐增大（见图）。垂直方向上与正常工况下相差不大（见图），年后污染物的影响深度也到达含水层底部，在防渗失效条件下，污染物在很短的时间内扩散的范围很大，且质量浓度升高较快，因此防渗墙的防渗性能很重要，要定期检查防渗墙的性能，若发现防渗失效须及时采取处理措施 。图特征污染物氨氮水平迁移图（非正常工况）（）图特征污染物氨氮垂直迁移图（非正常工况）（）结论（）通过 软件对南京溧水区某生活垃圾填埋场进行地下水污染物的运移模拟，预测特征污染物氨氮在平面上的污染趋势、范围及垂直方向上的污染程度，为污染场地的治理和修复提供了参考。（）通过两种工况下的模拟和预测结果，该垃圾填埋场在防渗失效后，污染物的运移主要受水力梯度影响，因此呈现出由西南到东北方向的运移，污染物 年内运移的最远污染距离达到 ，受污染范围为 ，垂直方向上污染物的影响深度到达含水层底部。因此，须做好防渗帷幕的跟踪监测。（）本次模拟仅在垃圾填埋场附近进行，本着风险最大化原则，垂直方向上只考虑污染物在潜水含水层中的运移，在计算时并未考虑土壤的吸附、解吸作用，以及地下水中化学反应和生物降解作用，因此，模拟结果具有一定的局限性，但由于在垃圾填埋场周围的渗滤液中氨氮质量浓度较高，具有代表性，因此模拟过程中选择污染源氨氮为定质量浓度 。参考文献：熊玲，鄢贵权 浅谈我国地下水的污染现状及防治措施 贵州科学，（）：金云龙，邱锦安，刘远锋，等地下水污染物运移数值模拟研究及应用综述地下水，（）：周维博，施坰林，杨路华地下水利用北京：中国水利水电出版社，于夷，