数值模拟在地下水环境影响评价中应用.docVIP

数值模拟在地下水环境影响评价中应用 　　摘 要：本文利用地下水数值模拟软件Visual Modflow，在详细研究某电厂周边水文地质条件的基础上，建立了研究区的地下水流动模型和地下水溶质运移模型，对电厂事故工况下泄露到地下水中轻质柴油的运移规律进行了模拟研究，就其对地下水的影响进行了分析评价。结果表明：模拟期内柴油对含水层造成了污染，并且存在超标现象，但没有对保护目标造成明显影响。结合数值模拟过程，提出了数值模拟在地下水评价应用中的一些注意事项，可为今后的地下水评价提供参考。 　　关键词：数值模拟；地下水；环境影响评价 　　中图分类号：X820.3 　　随着我国经济社会的发展，人们对环境质量的要求越来越高[1]，而同时，我国的环境污染问题，尤其是地下水环境污染的问题越来越严重。作为重要的环境因子，地下水已经成为城市和工农业用水的主要水源[2]。在经济社会可持续发展的前提下，防治地下水环境污染已引起人们的重视，全面、准确地评价建设项目对地下水环境的影响对于保护地下水资源意义重大[3]。环保部2011年6月发布了《环境影响评价技术导则―地下水环境》（HJ610-2011）[4]（以下简称《导则》）正式作为我国建设项目地下水环境影响评价工作的规范和指导，这就对地下水环境影响评价提出了更高的要求。 　　随着科学技术的发展，数值法越来越成为解决地下水问题的重要方法，数值模拟技术以其方便、灵活和高效的特点在水文地质领域得到广泛应用[5-7]。而《导则》要求，对于地下水环境影响一级评价，必须采用数值法进行预测和评价，对于地下水环境影响二级评价，可选择采用数值法或者解析法进行预测和评价。可见，数值法在地下水环境影响评价中占有重要地位。然而，由于自《导则》开始实施以来的时间较短，数值模拟在地下水环境影响评价中的应用还存在诸多方面的问题。本文拟通过建立某电厂厂址区域的地下水数值模型，利用Visual Modflow软件模拟事故工况下污染物在地下水中的迁移规律，进而就本工程对该区域地下水环境的影响进行预测和评价。在此基础上，对数值模拟在地下水环境影响评价应用中的一些问题进行分析，并给出应用中应注意的事项。 　　1 研究区概况 　　研究区位于桂林市西南部，属于低山丘陵地貌，位于低山丘陵冲积扇的中部和下部，地形北高南低，紧邻相思江。研究区的东、北部连绵土丘形成的孤形地下水和地表水分水岭，与西部大溪河，南部相思江一起形成一个典型的水文地质单元区。区内由第四系残坡积物及冲积层形成潜水含水系统，是可能受到工程污染的主要含水层。该层底板埋深一般在10～33m，在靠近丘陵地带为由粉质粘土混角砾组成的残坡积物，在靠近相思江和大溪河部分主要淤泥、砂砾及粉质粘土组成的冲积层。 　　研究区的计算面积为3.41km2，区内，潜水含水层分布广、厚度大，在常温常压下符合达西定律。地下水系统的输入和输出随时间、空间变化，故地下水为非稳定流；参数随空间的变化体现了系统的非均质性，同时具有明显的方向性，所以参数概化成各向异性。 　　2 地下水污染数值模型 　　2.1 数学模型的建立 　　2.3 地下水源汇项 　　模型源汇项包括补给项和排泄项。补给项主要包括降雨入渗量，排泄项主要为向河流排泄量和蒸发量。 　　2.4 定解条件的设定 　　研究区潜水流向由北向南，垂向上顶部潜水面接受大气降雨的补给，为一补给边界，同时地下水又通过其蒸发，是一排泄边界，将其概化为自由面边界。侧向上，研究区北部分为水岭边，可作为隔水边界；接受地下水排泄的相思江及西北部溪流伸出分水岭外，与地下水联系密切，水位变动不大，可概化为定水头边界。底部下伏泥页岩硅质岩，导水性弱，将其概化为隔水边界。 　　初始水位以丰水期（2012年5月）统测的地下水水位进行计算。以枯水期（2012年12月）统测的地下水流场作为模拟流场。 　　3 模型识别与校正及参数选取 　　参数确定是建立数值模型的关键[11]。水文地质参数分区原则主要参考地层和钻孔抽水试验资料，按含水介质的特征并考虑初始流场渗流特征进行划分。模型中水文地质参数主要为渗透系数K和给水度μ值，各水文地质参数按研究区内进行的水文地质试验的结果给出，并结合岩性特征和经验值给定初始值，通过模型模拟校正，最终获得模拟所需的水文地质参数。 　　经识别、验证，各含水层水文地质等效参数分区见图2，各区参数值见表1，观测值和计算值拟合程度见图3。 　　根据当地多处进行的弥散试验，弥散系数基本跟渗透系数成正比，推测评价区的纵向弥散系数为2.50 m2/d，横向弥散系数为0.25 m2/d，垂向弥散系数为0.025 m2/d。同时设定有效孔隙度为0.4～0.5。 　　4 污染物源强的确定 　　本文考虑厂址区可能出现的污染事故点对地下水造成污染的因素较复杂，在设计可能出现的事故情景时，重点考