第二讲MODFLOW直接建模ppt课件

1、第二讲第二讲2MODFLOW直接建模法直接建模法An Modular Finite Difference Ground-Water Flow ModellMODFLOW was originally documented by McDonald and Harbaugh (1984). As with most computer programs that are used over a long time period, MODFLOW underwent several overall updates. The second version of MODFLOW is documented

2、 in McDonald and Harbaugh (1988), and this version is often called MODFLOW-88 to distinguish it from other versions. A third version is called MODFLOW-96 (Harbaugh and McDonald, 2019a and 2019b).lMODFLOW was designed to have a modular structure that facilitates two primary objectives: ease of unders

3、tanding and ease of enhancing. Ease of understanding was an objective because U.S. Geological Survey technical managers generally believe that modelers should understand how a model works in order to use it properly. Ease of enhancement was an objective because experience showed that there was a con

4、tinuing need for new capabilities.例题例题第一步：运行软件两个步骤）l在桌面或菜单中点击GMS图标启动软件，在File|New创建新文件名并保存在选定目录。规定本次计算的单位，在菜单Edit|Units下选定长度以m(米)为单位，时间以d(天)为单位,其他单位本次计算不涉及。本例题用到以下两个模块：lGridlMODFLOW第二步：网格的创建l通过以下操作建立一个三维有限差分网格。l1.点击，打开3D模块。l2.在主菜单点击Grid|Create Grid开始创建3D网格。l3.选择X坐标部分，输入长度值为75000，15个网格数。l4.选择Y坐标部分，输入长

5、度值为75000，15个网格数。l5.选择Z坐标部分，输入长度值为15000，3个网格数。l6.点击OK。完成创建。第三步：MODFLOW模型的创建l3D几何模型建立后，就是初始化MODFLOW，并建立MODFLOW模型。这一步骤非常重要，切记。l选择MODFLOW | New Simulation，建立MODFLOW模型。第四步：全局变量的赋值第四步：全局变量的赋值 模块选择模块选择几种边界条件概化的解释几种边界条件概化的解释第五步：边界条件赋值第五步：边界条件赋值l边界条件赋值通过指针参数来实现，指针参数IBOUND0表示活动单元，IBOUND0表示定水头单元，IBOUND=0表示不活动单

6、元。通过选择各层各单元赋值。第六步：初始水头l点击Starting Heads按钮，设定初始水头。初始水头值用来表示非稳定流计算的初始值，在稳定流计算里无需赋值。但设定一个与计算结果相近的初始值可以提高计算效率，再者，如果初始水头值设定过小，软件便会认为含水层疏干。在本例中，我们设定初始水头值为软件缺省值0。如 第七步：含水层顶底板高程lGMS中上层底板与下层顶板为同一高程，所以我们赋值中仅需要对第一层的顶板和各层底板高程赋值。点击Top Elevation进入界面，如图3-5-5所示。选择第一层，点击Constant Layer按钮，赋值200。点击OK关闭顶板赋值软件包。l点击Bottom

7、 Elevation按钮进入底板高程赋值界面，分层对第一、二、三层各赋值为-150、-400和-700。点击OK离开全局变量软件包。边界条件直接赋值IBOUND赋值l1.选择MODFLOW|Cell Properties命令。l2.改变IBOUND选项到Spec.head。如图3-6-3所示。l3.单击OK按键退出网格属性对 第八步：LPF模块赋值lLPF是一个地下水流计算模块，有关详细情况可参考USGS有关参考书。选择MODFLOW|LPF命令打开模块赋值界面。选择Layer Data按钮进行含水层参数和类型赋值。在GMS中，含水层类型可以是承压的confined和非承压的convertib

8、le）,本例中最上层为非承压含水层，其他两层为承压含水层。第九步：降雨入渗模块输入l1.选择MODFLOW|Source/Sink Packages，并进入次级菜单Recharge Package。l2.选择Constant Array命令。l3.键入0.003，单击OK退出。第十步：排水模块输入第十步：排水模块输入l1.选择工具，选定i8，j2单元，按住shift键启动多选模式，同时选择在同一行中的3到10列单元。l2.选择MODFLOW | Sources/Sinks命令，点击Drain菜单，单击new加入排水单元。如图3-9-1所示。l3.每个单元水力传导率相同，但高程不同。输入高程及水

9、力传导率见表3-9-1。赋值结果见图3-9-2。第十一步：水井模块输入第十一步：水井模块输入l不同含水层水井位置不同，抽水量不同，应分层逐个输入。l1.按住Shift键，在第一层用鼠标选择图3-10-1所示水井所在单元。l2.选择MODFLOW|Sources/Sinks命令。l3.在下拉式菜单中选择Well.点击New创建水井。l4.为每一个水井输入流量-432000负值表示流出）。l5.点击OK退出。第十二步：检查、保存和运行模型第十二步：检查、保存和运行模型l1.选择MODFLOW|Check Simulation命令。l2.选择RUN，根据出现的提示检查模型。l3.待检查无误后，保存模型。l4.选择File|Save As命令。l5.选定保存路径和文件名。l6.选择保存按键。l待保存后，就可以运行MODFLOW了。选择MODFLOW|Run MODFLOW命令。总结l在数值模拟计算前必须选定适合的单位，GMS程序不能进行单位换算。lGMS操作界面中，MODFLOW菜单在3D Grid模块中。lMODFLOW计算所需模块在MODFLOW | Global Opti