visual_modflow实例

《visual_modflow实例》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、YourFirstGroundwaterModelwithPMWIN实例：一个由两种不同地层构成的含水系统，北部和南部为无流量边界，东部和西部为河流，与地下水有充分的水力联系，可视为定水头边界，西边界与东边界的水头分别为9m，8m。Fig.8-1Configurationofthesampleproblem含水系统为非承压各向同性，第一层、第二层的水力传导系数分别为0.0001m/s和0.0005m/s，两层的纵向传导系数假设为横向的10%，有效孔隙率为25%，地面标高为10m。第一层、第二层的厚度分别为4m、6m。含水层接受面状补给量为8×10－9mm。在第一层中靠近西部边界有一污

2、染源。现在我们的任务是，如何用一个位于东部边界附近的完整井来隔离该污染源。那么，在该区域内，需要建立一个数学模型来逐步计算抽水井所需的抽水量。只有当抽水量足够大时，被污染的地下水才能被抽水井全部截获。我们将使用PMWIN来建立数学模型，使用PMPATH来计算抽水井的截获带。基于所计算的地下水流场，使用MT3D和MOC3D来模拟污染源的运移情况。并演示如何使用PEST和UCODE来校正水流模型。最后创建一个动态演示方式，演示地下水被污染的发展变化过程。为了示范运移模型，我们假想污染物溶解于地下水中的速度为1×10-4µg/s/m2，含水层纵向弥散和橫向弥散系数分别为10m和1m，延迟系

3、数为2，初始浓度、分子扩散系数及衰减速率均为0。我们将计算三年之后污染浓度分布情况，显示在两个含水层中点[x,y]=[290,310],[390,310]的浓度－时间突变曲线。稳定流模拟稳定流模拟中的六个主要步骤：1.建立模拟模型2.模型赋值3.执行水流模拟4.模拟结果检验5.子区水均衡计算6.输出结果。步骤一：建立模型1、打开File菜单，选择New　Model。New　Model对话框弹出。选择一个保存模型数据的文件夹，如：C:PM5DATASAMPLE，键入文件名：SAMPLE作为示例模型。模型的文件扩展名必须是“.PM5”。在Window95/98/NT下，文件名的有效使

4、用字符数为120。最好在一个单独的文件夹中保存每一个模型及输出数据。可以同时运行多个模型（多任务处理）。2、点击OK。PMWIN只需几秒钟即可建立一个新模型，模型的文件名显示在标题栏中。步骤二：模型赋值包括：创建模型网格、定义边界条件、模型单元赋值。在整个模型中，PMWIN要求单位使用统一。例如，若以m表示长度单位，s作为时间单位，水力传导系数必须以m/s表示，抽水量必须以m3/s表示，弥散系数必须以m表示。在MODFLOW中，一个含水系统可以由一个包含一系列节点的离散区域或一些关联的有限差分块取代。图8-2，显示了一个包含网格和节点的含水系统的空间离散化情况，在每一节点均可计算出水

5、头值。节点网格形成了数值模型的框架。一个地质单元可由一个或多个模型层表示。每一计算单元格的厚度、长度、宽度均可改变。计算单元所处的位置可以用列、行、层来表示，PMWIN使用索引符号[J,I,K]来指定计算单元的位置。如：位于第2列、第6行、第1层的计算单元记为[2,6,1]。Fig.8-2Spatialdisretizationofanaquifersystemandthecellindices一、创建模型网格1、打开Grid菜单，选择MeshSize。ModelDimension对话框弹出。2、输入含水层的层数3，行数30、列数30，行高、列宽20。第一含水层和第二含水层分别由一个

6、模型层和两个模型层来代替。3、点击OK。PMWIN的界面发生变化，此时显示了模型的整体网格。PMWIN允许用户改变和旋转模型网格。用户可以随意改变模型网格行、列宽度，也可以增加/删除行和列。在该实例中，不需更改模型的网格。4、从File菜单中选择LeaveEditor或点击LeaveEditor按钮退出。二、定义含水层类型1、打开Grid菜单，选择LayerType。LayerOptions对话框弹出。2、点击标签为Type列中的某一单元格，单元格中将出现一个带下箭头的按钮，点击该按钮，出现一个可供选择含水层类型的列表。3、为第一层含水层选取1：Unconfined，其它含水层选取0

7、：Confined，然后点击OK，关闭对话框。三、定义边界条件在模型中，边界上的各个计算单元均用不同的指定代码表示不同的边界类型：（1）“1”表示变水头边界，（2）“－1”表示定水头边界，（3）“0”表示无流量边界。1、从Grid菜单中选择BoundryCondition→IBOUND(Modflow)。PMWIN的DataEditor界面弹出，界面显示一个模型网格平面图，单元格指针位于[1,1,1]，即第一层的左上角。当前格的数值显示在状态栏底部。IBO