古山坑水源地三维流三维流分析.docxVIP

古山坑水源地三维流三维流分析 近年来，随着经济的加快、工业生产的扩大，该地区水资源的需求增加。人们对环境保护、水资源浪费和合理使用水资源的认识不强，导致了一些环境问题，如河流中断、湖泊干燥、地下水位持续下降等。20世纪60年代中期以来, 随着快速大容量计算机的出现和应用, 地下水资源分析评价中数值计算方法得到逐步推广, 具有明显的通用性和适用性, 尤其近十几年, 地下水数值模拟取得了长足进步计算机技术的发展带动了数值模拟软件的发展, 目前GMS、MODFLOW、FEFLOW、Visual MODFLOW等模拟软件应用越来越广泛GMS是由美国军队排水工程工作站和美国Brigham Young大学的环境模型研究实验室开发的一个用于地下水模拟的图形界面软件, 它可以提供多种地下水数值模型的方法, 并具有能准确的显示地层的空间结构等优点。因此, 通过GMS软件对当地地下水资源开采进行模拟, 对地下水合理开发和生态环境保护具有重要的现实意义。 1 孔隙—基本概况 研究区炼油厂厂区一带黄土层富水性较差, 地下水主要赋存在下部白垩系下统洛河砂岩中, 地下水类型为白垩系洛河砂岩孔隙—裂隙承压水;厂区北部的平原区, 赋存第四系孔隙潜水和洛河砂岩孔隙—裂隙承压水 根据本次实地调绘、结合《靖边县水资源调查评价报告》, 厂区附近白垩系碎屑岩类孔隙—裂隙承压水的单井出水量500~1000 m3/d, 属水量较丰富区, 水质符合一般工业用水;县城及北部一带上部为第四系松散岩类孔隙潜水, 水量较丰富区分布于双伙场—关路畔—芦河酒厂线以北, 以南为水量中等区, 下部为白垩系碎屑岩类孔隙—裂隙潜水和承压水, 单井出水量1 000m3/d, 属水量丰富区, 除县城一带水质较差, 北部其他地区均为好的工业用水和生活饮用水。 2 模型的构建 2.1 径流补给量计算 研究区的边界条件, 因第四系潜水含水层与白垩系承压水含水层可视为无限含水层, 所以将补给边界定为南边界和西边界, 排泄边界定为北边界和东边界;本次模拟, 因含水层的南边界径流补给量稳定, 可视为定流量边界, 而径流补给量较小西边界计算时以零流量边界处理, 东边界和北边界则处理为变流量边界。区域内芦河在计算时分为上下两段, 上段补给, 但河水下渗补给很小所以忽略, 下段则为排泄, 为了模拟地下水的泄出量将其设置为排水沟。 2.2 含水层顶底板高程的插值 1) 网格划分。利用GMS的3D GRID模块, 实现网格剖分。计算区为包含整个研究区的长方形区域, 剖分网格总数7 904, 其中有效计算网格5 999。 2) 网格高程插值。将所建模型运行后, 通过统一坐标系, 各层顶底板高程以x、y、z显示, 然后输入到2D scatter points模块中, 再内插到MODFLOW Layers中。插值过程中如果出现错误或者不合理的地方, 可以运行Model Checker命令改正。研究区含水层顶底板标高。 3) 初始流场。根据《靖边县西部地区1:5万供水水文地质详查》报告, 得到计算区现状条件下的潜水与承压水等水位 (头) 线图。在此基础上, 依据天然补排关系和补排量 (不含现状开采量) , 得到了天然状态下的地下水流场, 天然状态吓得地下水流场可作为初始流场进行开采降深场的预测。主要供水目的层的地下水天然状态下初始流场。 2.3 拟合过程模拟 为检验所建立的数学模型和模型参数的可靠性, 根据野外数据试验所获得的数据, 选择2012年9月18日到9月30日作为检验时段, 3天为一个步长。 通过把各项原始数据输入模型, 计算出水头值并和实际观测值进行对比, 发现拟合结果良好, 误差在合理的范围内, 所以该模型可以用于预报。拟合过程线图。验证后的参数见表1, 表2。 2.4 地下水流场与降深场预测计算 以炼油厂原油常压加工项目取用水量:0.84×104m3/d (即306.6×104m3/a) 为保证, 根据《靖边县水资源评价与规划报告》可知近远期规划研究区内农灌开采量不增加, 现状开采量保持不变, 从而设计以下开采方案表3。 古山坑水源炼油厂生产区现状开采井有14眼, 全部开采洛河组承压水, 井深均为300 m, 井间距为300~500 m, 单井出水量为900~1 200 m3/d; 1) 方案一地下水位预测计算。在现状开采的基础上, 按照0.84×104m3/d (即306.6×104m3/a) 进行开采, 时长30年。该开采方案下地下水位 (头) 历时变化以及末流场与降深场的空间变化详。 2) 方案二地下水位预测计算。在现状开采的基础上, 按照1.37×104m3/d (即500.05×104m3/a) 进行开采, 时长30年。该开采方案下地下水位 (头) 历时变化以及末流场与降深场的空间变化。 3) 方案