modflow飞机场作业.docx

《modflow飞机场作业.docx》由会员分享，可在线阅读，更多相关《modflow飞机场作业.docx（9页珍藏版）》请在课桌文档上搜索。

1、1、飞机场水文地质简介模拟场所位于Waterloo市外距离飞机场很近的地方。模型区域内有浅层的砂和砂砾含水层，深层砂和砂砾含水层，中间粘土和淤泥半透水层将上下两个含水层分割。该区域四周有飞机燃料补充区，城市供水水源地，以及一个不连续的半透水层区。城市供水水源地有两口供水井。东边水井的流量保持在550m3d,西边水井的流量是每天400m3在过去的十年中，飞机燃料定时地排入到燃料区，自然渗透作用使得上层含水层产生了羽状污染区。该指南将会说明建立该地区地下水水流和污染物运移模型的必要步骤。该模型将会指出燃料污染对城市供水水源地潜在的影响。当从平面角度上探讨该模型时，模型区域顶部设定为北部，模型底部为

2、南部，左边和右边分别为西部和南部。地下水流在由浅层潜水含水层，中间弱透水层以及深层承压含水层组成的三层含水系统中从北流向南（从顶部原委部）。浅层和深层含水层的渗透系数为2e-4msec,弱透水层的渗透系数为le-4msec02、自然状态下水均衡分析自然状态水均衡分析表：第一层其次层第三层合计数量比例数量比例数量比例数量补降1077.50.99990.000000.000001077.5给项雨边界流入0.02180.000020.0000087.0100.320387.0318越流0.120420.OOOll184.451184.640.67970小计1077.61184.451271.6511

3、164.532排泄项河流867.430.803470.000000.00000867.43边界27.7890.025740.00000281.950.99957309.739开采0.000000.000000.000000越流184.330.17074184.7610.120290.000430小计1079.61184.761282.0711177.169将该模拟区含水层系统分为三层，由顶层潜水含水层、中间一个弱透水层以及底部的承压含水层组成。由于北部水头高，因此地下水的流淌趋势为由北向南。分析模拟区的各排泄项，潜水层的排泄主要通过河流以及越流，弱透水层的排泄基本全部是通过越流排泄，而对于承压

4、含水层来说，绝大部分是通过边界流出。3、现状年下水均衡分析现状开采条件下的地下水均衡表：第一层其次层第三层合计数量比例数量比例数量比例数量补给项降雨1077.50.971990.000000.000001077.5边界流入30.9990.027960.00000635.910.65849666.909越流0.05620.00005329.861329.800.34150小IT1108.51329.861965.7111744.409排泄项、川流786.560.704550.000000.00000786.56边0.06060.0000540.0000028.2330.0288628.2936界

5、开采0.000000.00000950.010.9711950.01越流329.800.2954329.8510.56139E-010.0000570小111116.41329.851978.3011764.864我们可以从现状年表中可以发觉，他及自然状态下的分层均衡表作对比，首先来看潜水含水层，边界流入较自然状况下有少量增加。中部弱透水层的状况及自然状况下相类似，但是越流量从184.45增加到329.86,这样的状况说明在开采状况下可以变更地下水的运动速度。对于承压层，边界的补给大大地增加，从87.010增加到635.91,越流也有所增加，但是幅度没有边界补给大。4、设计新增开采方案方案一：

6、三口井的抽水量：单位：11i7d开采井开采量Screenlevelbottom(m)screenleveltop(m)supplywell330047supplywell430047supplywell530047三口井坐标:井坐标XYsupplywell3326323supplywe114564321supplywell5838274Screenlevelbottom4m,screenlevel7m方案一中可以看到质点示踪和溶质运移都会进入supplywell3,说明supplywe113影响到地下流场的分布和水头的变更，从来影响溶质简洁的方向，二十年稳定流后污染羽将进入新开采井。所以否定方

7、案一。方案二：井坐标XYScreenlevelbottomscreenleveltopsupplywell337712947supplywe114101311947supplywell5149011947三口井开采量：单位：m7D开采井模拟一模拟二supplywell3300400supplywell4300400supplywell5300400方案二中从宏观上考虑水头是北高南低，而河流又是补给地下水的，将三口井打在河流的南部，可能对地下水流场的变更对污染羽的影响较小。我们从方案二的质点示踪图和溶质运移图中可以看到新设计的三口井没有受到影响。所以我们可以选择方案二进行抽水。方案二中当抽水量从

8、300立方增加到400立方时，也没有受到污染羽的影响。所以我们接受方案二。方案二中三口井抽水400m/D条件下的水均衡分析：第一层其次层第三层合廿数量比例数量比例数量比例数费补给项降雨1077.50.9573522880.000000.00000107.5边界流入47.9940.0426423810.000001765.10.8264737181：J.OS越流0.0282282.50804E-05370.50+0.0284341370.600.1735220小计1125.51370.5312135.71289().5S排泄项河流761.720.6727786610.000000.0000076

9、：.72边界0.000000.000000.0218111.014E-050.0：!181开采0.000000.0001602150.00.9999534912：50越流370.500.327239004370.60+0.02822810.0284341.322E-050小i1132.21370.6312150.11291：.74结果分析：一、补给项从自然状态下1164.532,到现状年下1744.409,再到方案二中的2890.594,我们可以发觉补给项随着水井开采量的增加而增加。二、随着开采量的增加，对河流的排泄从867.43削减到761.720三、各个层之间的越流量也越来越大。四、排泄项

10、中边界量从309.739削减到0.021811,我们可以看出井的开采对地下水流场的影响。5、模拟运行结果图自然状态下等水头线图：现状年下第1层质点示踪图：现状年下第2层质点示踪图：现状年下第3层等水头线图：现状年下第3层速度矢量图：现状年下第5层等水头线图：现状年下第5层质点示踪图：现状年下第5层速度矢量图：现状年10年溶质运移图：现状年20年运移图23列示图：现状年第六层20年溶质运移图:现状年第四层20年溶质运移图:现状年第一层20年溶质运移图:增加三口水井方案一、方案一其次层质点示踪图：方案一第5层质点示踪图：方案一第6层质点示踪图：方案一其次层20溶质运移图：方案一第5层20溶质运移图

11、：方案二第4层等水头线图：方案二第6层等水头线图：方案二第三层等水头线图：方案二第六层质点示踪图：方案二第四层溶质运移图：方案二第六层溶质运移图：方案二增加到300开采量第五层溶质运移图：6、MOdfloW软件的特点VisualModfIoW（简称VM）可视化三维地下水流淌模型软件，是目前最流行并被各国机构认可的三维地下水水流和溶质数值模拟评价的标准可视化专业软件系统，在地下水模拟领域得到广泛的应用。VisualModflow在我国地下水探讨中应用潜力巨大,特殊是流线示踪模拟和随意均衡域的水均衡计算等功能，对于科学制定水资源的可开采调蓄量和矿井防治水决策方案具有极其重要的理论指导意义和好用价值

12、。这个软件包由Modflow（水流评价）、Modpath（平面和剖面流线示踪分析）和MT3D（溶质运移评价）三大部分组成，并且具有强大的图形可视界面功能。设计新颖的菜单结构允许用户特殊简洁地在计算机上干脆圈定模型区域和剖分计算单元，并可便利地为各剖分单元和边界条件干脆在机上赋值，做到真正的人机对话。假如剖分不太志向须要修改时，用户可选择有关菜单干脆加密或删除局部网格以达到满意为止。同时，用户可选用不同的菜单分别单独或共同运行ModflowModpath和MT3D三大部分。各部分均设计了模型识别和校正的菜单。本软件包可便利地以平面和剖面两种方式彩色立体显示计算模型的剖分网格、输入参数和输入结果。

13、VisualModflow的模块化结构便于用户依据实际须要添加程序和完善功能。界面分为输入、运行和输出模块。全部模块都支持运用底图，这些底图可以是JPG、BMP或GIF等格式。输入模块的功能是完成初始条件和边界条件的输入。运行模块的功能是供应多种算法模拟计算。输出模块的功能是输出直观图形结果，并可便利地计算出三维流线分布和随意时间水质点的移动位置。这个软件系统的最大特点是将数值模拟评价过程中的各个步骤完备无缺似地连接起来，从起先建模、输入和修改各类水文地质参数及几何参数、运行模型、反演校正参数，始终到显示输出结果，使整个过程从头至尾系统化、规范化。7、运用Modflow软件的体会很多时候输错参数就不知道从哪里撤销，只能重新修改，这样感觉不是那么便捷。各种参数很难在图上体现出来，我们要去。Verlay中设置，操作界面立体可视化效果比较弱。