Skip to main content

温州市地下水Modflow6趋势模拟

info

本页面中的短视频为简短操作介绍，获取完整EnviFusion教学课程，请访问 在线课程。

算例工程下载地址

1 算例简介 2 操作步骤 2.1 新建项目 2.2 创建MODFLOW模型 2.3 定义网格高程 2.4 定义活动区域 2.5 定义渗透系数 2.6 定义定水头边界 2.7 定义河流边界 2.8 定义补给边界 2.9 设置应力期 2.10 解释运行MODFLOW 3 小结

算例简介

使用上一算例中获取的温州市剖分网格和各项输入文件，构建描述地下水流动趋势方向的ModFlow模型。由于基础数据的稀缺，此模型仅用来观察地下水资源的潜在运移方向，其流量并不能直接用来进行水资源量衡算。

操作步骤

新建项目​

  1. 打开Envifusion软件，若软件界面仍有其他项目，点击菜单栏按钮关闭当前项目文件。

创建MODFLOW6模型​

  1. 点击菜单栏中的导入数据按钮，打开算例文件夹中的A04-区域网格.efd文件，生成A04-区域网格.efd图层，在该图层属性面板中点击应用完成设置。导入EnviFusion制作的区域网格，在视图1中显示如下：

  2. 在管道面板中，选中A04-区域网格.efd图层，右击快速操作，输入操作命令：初始化，选中初始化MODFLOW6命令。在属性面板中，设置模型层数：1， 层厚：50，层类型：含水厚度可变，水平渗透系数：1。单击应用，创建MODFLOW6模型图层。

定义网格高程​

使用第一部分教程中创建好的高程文件定义网格高程。

  1. 点击菜单栏中的导入数据按钮，打开算例文件夹中的A01-区域DEM.efd文件，生成A01-区域DEM.efd图层，在该图层属性面板中点击应用完成设置。

  2. 选择初始化MODFLOW6模型1图层，点击数值模拟-> MODFLOW6->MF6网格高程，选择A01-区域DEM.efd图层为高程图层，点击确认生成MF6网格高程1图层，将其重命名为010-顶板高程-MF6网格高程1图层。在属性面板中，选择赋值模式为网格属性->单层，高程数列为band_1，赋值对象为模型顶板，最小层厚为0.1，勾选上层优先，选择固定层厚，设置终止层号为1，固定层厚为50，单击应用。在渲染子面板中，将渲染方式设为曲面，并将坐标变换中Z方向的缩放值设置为5，可在视图1窗口显示赋值高程后的模型：

定义补给边界​

  1. 选中A04-区域网格.efd图层，对其进行快速操作，选择命令：采样，采样体为A04-区域网格.efd图层，数据源为A01-区域DEM.efd图层，生成采样1图层。

  2. 基于采样1图层，进行阈值裁切操作，将平原区范围裁切出来。在阈值裁切1图层的属性面板中，选择裁切数列band_1，并设置裁切范围为（-4，46），单击应用。

  3. 在管道面板中，选中010顶板高程图层，对其进行快速操作：MF6多边形赋值，属性图层选择阈值裁切1图层，生成MF6多边形赋值1图层，重命名为020垂向入渗边界-MF6多边形赋值1图层，在属性面板中，选中垂向入渗边界，并设定层号为1，垂向入渗率设置为0.000486，单击应用。

  4. 重复步骤1-3，赋值山区垂向渗透系数。

  5. 基于采样1图层，进行阈值裁切操作，将平原区范围裁切出来。在阈值裁切1图层的属性面板中，选择裁切数列band_1，并设置裁切范围为（46，1520），单击应用。

  6. 在管道面板中，选中020垂向入渗边界图层，对其进行快速操作：MF6多边形赋值，属性图层选择阈值裁切2图层，生成MF6多边形赋值2图层，重命名为030垂向入渗边界-MF6多边形赋值2图层，在属性面板中，选中垂向入渗边界，并设定层号为1，垂向入渗率设置为0.00062，单击应用。

定义定水头边界​

  1. 在图层面板中只保留A-04区域网格.efd图层的显示，点击俯拾视角：Z轴向外工具，并单击视图1中的3D 按钮，切换为2D 查看。

  2. 点击创建数据体->几何体->多边形，生成多边形1图层。在该图层属性面板数据子面板中，勾选显示手柄选项，将添加类型设为折线。

    在视图窗口中，按住键盘Alt+Ctrl键，单击鼠标，沿着研究区域东侧海岸线的区域边缘绘制折线（具体位置见下图）。折线创建起始有（0，0，0）和（1，0，0）两个控制点，在折线创建完毕后，可在控制点列表中选中这两个控制点，点击表格右侧的删除按钮，删除起始两点 ，点击应用完成设置。

  3. 在图层面板中选择030垂向入渗边界-MF6层赋值1图层，点击数值模拟->MODFLOW6水流->MF6线赋值，属性图层选择多边形1图层，生成MF6线赋值1图层，将该图层重命名为040定水头边界-MF6线赋值1图层。在该图层属性面板数据子面板中，属性类别设为定水头边界，层号为1，固定水头指定值为0，高程预警选择底部高程。点击应用完成设置。调试窗口随即弹出提示，当检测到设置的水头值0低于网格底部高程时，直接将水头值替换为网格底部高程值。

    在视图1中可以看到红色区域为定水头边界。

定义排水沟边界​

  1. 在图层面板中只保留A-04区域网格.efd图层的显示，点击俯拾视角：Z轴向外工具，并单击视图1中的3D 按钮，切换为2D 查看。

  2. 单击导入数据按键，导入文件A-02区内河流.efd文件，单击应用，生成A-02区内河流.efd图层。

  3. 在图层面板中选择040定水头边界图层，点击数值模拟->MODFLOW6水流->MF6线赋值，属性图层选择A-02区内河流.efd图层，生成MF6线赋值2图层，将该图层重命名为050排水沟边界-MF6线赋值2图层。在该图层属性面板数据子面板中，属性类别设为排水沟边界，层号为1，排水沟高程设置为-50，传导系数为3，单击应用。高程预警为底部高程。调试窗口随即弹出提示，当检测到设置的排水沟高程低于网格底部高程时，直接替换为网格底部高程值。

解释运行MODFLOW6​

解释运行创建好的水流模型。

  1. 选中050排水沟边界图层，点击数值模拟-> MODFLOW6水流->MF6解释运行，生成解释运行MF6解释运行1图层。在属性面板中，点击按钮选择模型文件生成的位置，将模型名保存为温州.mf6。选择复杂度为非线性模型，时间单位为天，长度单位为米，勾选保存所有时间步、保存水均衡、点击应用，系统将自动编译并运行MODFLOW6模型，运行完成后会弹出调试窗口，显示MODFLOW水流模型引擎正常结束。

  2. 点击菜单栏中的导入数据按钮，选择上一步中用户指定名称的MODFLOW6文件夹中的efModFlow6.nam6文件。图层面板中生成efModFlow6.nam6图层。

    在该图层属性面板将坐标变换属性下的Z方向缩放设为5，点击应用导入MODFLOW6模型结果。

    通过选择面板中水头标签的勾选，结合视图1下的网格或者散点选择工具，可以具体查看各个网格/散点的水头值。

  3. 选中efModFlow6.nam6图层，进行图层复制，生成图层复制1图层。将渲染方式设为矢量流线，点击应用完成设置。模拟区动态流场如下图所示。

计算面积​

对山区和平原区分别进行面积计算。

  1. 在管道面板中，选中阈值裁切1图层，依次进行快速操作：剖分、几何计算和积分。

  2. 在信息面板中，可以看到具体的面积数据。

  3. 重复上述操作，对代表山区的阈值裁切2图层，依次进行依次进行快速操作：剖分、几何计算和积分。从而获取山区面积。

井水位标注​

  1. 单击导入数据按键，导入A03-区内监测井文件，生成A03-区内监测井图层。

  2. 在管道面板中，选中A03-区内监测井图层，进行快速操作：计算器。在属性面板中，输入结果数列名称：水位标高，公式为孔口标高-水位，单击应用。

  3. 在管道面板中，选中计算器1图层，进行快速操作：散点->文字标注。在属性面板中，设置标注数列为：水位标高，并设置相关格式，单击应用。

小结

完成本指南中的算例后，可以学到以下内容：

  • 如何创建并运行地下水MODFLOW6水流模型；
  • 如何使用efd文件为MODFLOW模型赋值；
  • 如何使用csv文件进行为MODFLOW6模型赋值；
  • 如何基于场地地图绘制多边形和折线；
  • 如何对折线进行单一属性赋值；
  • 如何保存和导入efd图层文件。
  • 如何计算面积。