2.1 Conceptual Modeling Tutorial

以下示例是构建概念模型并将其转换为数值模型基础知识的快速演示。

Objectives

•学习如何创建一个项目并导入原始数据

•熟悉导航图形用户界面（GUI）和概念建模步骤

•学习如何定义3D地质模型和流动属性

•使用GIS数据定义边界条件

•定义MODFLOW网格，然后用概念模型的数据填充该网格

•查看生成的属性和边界条件

•将模型输入转换为MODFLOW包并运行MODFLOW引擎

•通过在几个视图中解释水头和降水来理解结果

•通过将观察到的水头与计算水头进行比较来检查模型的质量

Required Files

此练习需要使用多个文件，这些文件应该已包含在Visual MODFLOW Flex安装中。

这些文件通常位于公共"My Documents"文件夹中，路径为：

C:\Users\Public\Documents\Visual MODFLOW Flex\Tutorials\Conceptual Model\supp files

如果您无法找到这些文件，请从我们的网站上下载。

Creating the Project

•启动Visual MODFLOW Flex

.

•选择[文件]，然后选择[新建项目...]。将出现“创建项目”对话框。

•在项目名称中键入“Conceptual Modeling Tutorial”。

•点击[

] 按钮，在弹出的对话框中导航到您希望保存项目的文件夹，然后点击​​。[OK].

•勾选“为项目创建一个文件夹”复选框。

•定义您的坐标系统和基准面（或保留默认值 - 本地笛卡尔）。

•对于此项目，默认单位将是合适的。

现在，创建项目对话框应该如下所示（确保单位相同）：

•点击【确定】。将出现工作流程选择屏幕。

•选择【概念建模】，概念建模工作流程将加载。

•在“定义建模目标”步骤中，您可以定义模型的目标和默认参数。

 

•模型的“开始日期”对应于模拟时间段的开始。定义一个相关的开始日期很重要，因为您的现场测量数据（观测水头和抽水计划）将使用绝对日期进行定义，并且必须位于模拟时间段内。对于这种情况，使用默认目标即可。

Start Date

开始日期将用于获取模型运行的抽水井和水头/浓度观测数据。当您使用绝对（日历）日期定义井数据时，确保您的开始日期反映了模型运行的实际开始时间非常重要。井数据必须在该开始日期或之后。否则，这些数据将不包括在模拟中。 另外，一旦设置开始日期，就无法更改。如果您错误地设置了错误的开始日期，可以使用相对时间（例如从0开始）导入您的抽水井数据和观测数据，而在数值模型输入/输出中看不到任何区别。

Setting the Start Date

此练习的模型开始日期应设置为2000年1月1日。Visual MODFLOW Flex使用标准的Windows日期选择器；以下是如何使用它的一些建议。点击下面显示的按钮，以加载Windows日期选择器。

将出现标准的Windows日历。点击标题中的月份（如下所示）：

所有当前年份的月份将如下所示出现。点击标题中的年份：

然后将出现一系列年份，如下所示。点击标题中的年份范围：

将显示上一个十年的年份列表。然后，您可以使用<或>按钮更改年份。选择2000-2009时期：

一旦到达所需的十年，选择日历上的“2000”如下所示：

然后将显示该年的月份列表。对于这个例子，选择一月，如下所示：

最后，从日历中选择“1”，如下所示：

选择的日期（2000年1月1日）将出现在开始日期中。

•对于这种情况，默认目标和数值的其余部分都是合适的。

•点击[

] （下一步）继续。

Collect Data Objects

•下一步是导入或创建您将用于构建概念模型的数据对象。

 

在这一步，您可以导入数据，创建新的数据对象（通过数字化）或创建表面（从点数据对象创建）。

•点击【导入数据】按钮，将加载以下屏幕：

•在数据类型组合框中选择“多边形”。

•在源文件字段中，点击【...】按钮，导航到您的“My Documents”文件夹，然后选择“Visual MODFLOW Flex\Tutorials\Conceptual Model\supp files”，并选择“boundary.shp”。

请注意：如果您在安装时选择使程序对所有用户可用，则这些文件可能位于公共文档文件夹中：“C:\Users\Public\Documents\Visual MODFLOW Flex\Tutorials\Conceptual Model\supp files”中。

•点击【下一步>>】。

•点击【下一步>>】。

•点击【下一步>>】然后点击【完成】。

•下一步是导入表示地表的表面。

•点击【导入数据】按钮。

•为数据类型选择【Surface】。

•在源文件字段中，点击【...】按钮，导航到“我的文档”文件夹，然后到“VMODFlex\Tutorials\ConceptualModel\supp files”文件夹，并选择“ground.grd”。

•点击【下一步>>】通过所有屏幕以接受默认值，然后点击【完成】。

•重复这些步骤以导入剩余的表面：layer2-top.grd，layer2-bottom.grd。

•接下来，导入折线数据对象，从相同的源目录选择chd-east.shp；使用所有默认值并完成导入。

•重复这些步骤，对于折线，首先导入chd-west.shp，然后是river.shp。

•一旦数据对象被导入，它们将出现在程序窗口左侧的树中。

 

•您可以在2D或3D中查看这些数据对象；只需创建一个新的查看器：

•点击主菜单上的【窗口】，然后选择【新建3D窗口】；将出现一个空的3D查看器；

•点击每个导入的数据对象旁边的复选框，它们将出现在3D查看器中。

•增加垂直夸张，并使用鼠标重新定位屏幕，以便您可以看到所有的数据对象，如下图所示。

•通过点击标签的“X”关闭3D视图。

. tab

••在概念模型标签上点击【

】以（下一步）继续，您将到达“定义概念模型”步骤。

Define Conceptual Model

•为概念模型提供一个名称（例如，Conceptual Model 1），并指定模型区域。

•从数据资源管理器中选择名为“boundary”的多边形数据对象（表示概念模型的水平范围），使其突出显示：

•点击“模型区域”中的【

】按钮，将其插入为在“定义概念模型”工作流程步骤中的参考对象。

 
 
 

请注意：模型区域不能使用复杂多边形（与自身相交的多边形）定义，也不能使用包含多个多边形的多边形数据对象。

•点击【保存】。注意，概念模型元素已添加到应用程序窗口左下部的“模型资源管理器”树中。

•点击【

】（下一步）以继续到“定义模型结构”步骤。

Define Structure

定义地质模型包括提供地层的地质表面作为层位的输入。然后，在这些层位之间创建三维实体。按照以下步骤创建新的层位。

•从“层位设置”对话框（如下所示）中，点击【

】（添加层位）按钮，向“层位信息”表中添加一行新的层位。

•重复此步骤两次，以便在层位表上有3行新的记录。

•从数据资源管理器中选择将用于生成层位的“ground”表面数据对象。

•点击Horizons表的Row1中的【

】（插入）按钮，将其插入到层位信息表中。见下图示例。

•对于这个例子，使用默认的层位类型就足够了。有关每个层位类型的解释，请参阅用户手册的“层位类型”部分。

•重复上述步骤以添加额外的层位：

•从数据资源管理器中选择“layer2-top”表面数据对象，在Horizons表的Row2中点击【

】（插入）按钮，将其插入到表格中。

•从数据资源管理器中选择“layer2-bottom”表面数据对象，在Horizons表的Row3中点击【

】（插入）按钮，将其插入到表格中。
 

请注意：层位按照从最上方的地质接触到最下方的工作顺序添加。

•您可以通过点击【预览层位】按钮在相邻的3D查看器中预览层位。

•完成后，您应该会看到一个类似下图所示的显示。

】（下一步）按钮以继续下一步

•点击【

】（创建并保存）按钮以生成模型的层位，您应该会看到它们填充到模型资源管理器中。

•最后，点击【

】（下一步）按钮以继续下一步

（如果您还没有点击【创建并保存】按钮，点击【（下一步）】按钮也会自动在模型资源管理器中生成模型的层位）。

Define Property Zones

一旦您将足够的原始数据导入到项目中，您就可以开始使用导入或数字化的数据对象构建一个或多个概念模型，作为构建块。

在这一步，您可以查看/编辑模型的流动性质。有两种定义属性区域的方法：使用结构区域或使用多边形数据对象。对于本教程，我们将使用结构区域。

Using Structural Zone(s)

这种方法允许您使用在“定义模型结构”步骤中生成的结构区域生成属性区域。

•点击如下所示的【使用结构区域】按钮。

•从概念模型树中选择结构区域“Zone1”（如下所示，位于结构/区域节点下）。

 

•点击【

】（插入）按钮将该区域插入到结构区域字段中，如下图所示。

•选择要定义的参数组，例如导电率、储存或初始水头。下方的数据输入表格将根据所选择的参数组显示相应的参数。例如，如果选择导电率，数据输入表格将显示参数Kx、Ky和Kz。数据输入表格将使用项目设置中指定的默认值进行填充（可以从主菜单中点击【文件】>【项目设置...】访问）。

•为属性区域1输入以下值：Kx = 4E-6，Ky = 4E-6，Kz = 4E-7。

•点击窗口右侧的【保存】按钮。

•为另一个属性区域重复这些步骤：

•点击【使用结构区域】按钮

•从模型树中选择“Zone2”

•点击【

】（插入）按钮将该区域插入到结构区域字段中，如下图所示。

•为属性区域2输入以下值：Kx = 7E-5，Ky = 7E-5，Kz = 7E-6

•点击窗口右侧的【保存】按钮。

•属性区域还可以使用多边形形状定义；值还可以从shapefile属性或2D表面（分布式值）定义。有关更多详细信息，请参阅[定义属性区域](https://resources.environmental-expert.com/dictionary/environmental-terms/defining-property-zones-18225)。

•点击【

】（下一步）以继续到选择屏幕。

•点击【定义边界条件】按钮以继续。

在此窗口中，您可以选择要创建的边界条件的类型：标准的MODFLOW边界条件（CHD、DRN、RCH等）、抽水井或围墙。

•点击【定义其他边界条件】按钮：

•“定义边界条件”对话框将出现在屏幕上，如下一节所述。

Define Boundary Conditions

•在这一步，您可以为模型定义流动边界。

•从“选择边界条件类型”组合框中选择所需的边界条件类型。

•选择“定值水头（Type 1）”

•输入名称：“Constant Head East”。

•从数据资源管理器中选择表示此定值水头的“chd-east”折线。

•点击“定义边界条件”对话框中的【

】（插入）按钮，将此折线添加到输入。

•点击【下一步】按钮。

•下一个对话框允许我们定义定值水头。Visual MODFLOW Flex提供了定义边界条件属性的各种选项。属性可以从存储在Surface、Time Schedule、Shapefile和3D Gridded数据对象中的属性中分配。您还可以将属性设置为静态（随时间不变）或瞬态（随时间变化）。

•对于本教程，您将分配一个静态定值水头。

•在“起始水头（m）”和“结束水头（m）”字段下的空字段中输入“347”。

•点击【完成】按钮

重复这些步骤以定义其他定值水头边界条件：

•在工作流程树中直接点击【定义边界条件】。

•选择【定义其他边界条件】按钮。

•选择定值水头，选择“chd-west”折线，并为起始水头和结束水头都定义一个值为325。

•点击【完成】。

在继续之前，您将定义另一个边界条件，即河流。

•在树中点击【定义边界条件】，并选择【定义其他边界条件】按钮。

•选择河流（Type 3 - 仅限MODFLOW）作为边界条件类型

•从数据资源管理器中选择“river”折线。

•点击“定义边界条件”对话框中的【

】（插入）按钮，将此折线添加到输入。

•可能会出现关于裁剪折线的警告；点击【确定】继续。

•点击【下一步>>】按钮。

•为河流定义以下属性，如下所示：Stage = 335 (m)，Bottom = 333 (m)，Riverbed Thickness = 1 (m)，Width = 10 (m)，Riverbed conductivity = 0.01 (m/s)。

•点击【完成】。

•河流概念边界条件将添加到模型树中。

•将出现以下显示。

接下来，您可以选择创建什么类型的网格：

•定义有限差分网格：用于MODFLOW-6、MODFLOW-2000、-2005或MODFLOW-LGR模型运行；

•定义有限元网格：用于为FEFLOW .FEM文件准备输入；

•定义非结构化V网格：用于使用Voronoi网格进行MODFLOW-6或MODFLOW-USG运行；或

•定义非结构化Q网格：用于使用Quadtree网格进行MODFLOW-6或MODFLOW-USG运行。

•点击【定义有限差分网格】按钮，将出现以下窗口；按照下一节中的说明定义输入。

Define Finite Difference Grid

•为数值网格输入一个唯一的名称。此名称将在创建网格后出现在概念模型树中。

•输入网格大小，以及可选的网格旋转。通过在旋转文本字段中输入0到360之间的值，可以使网格围绕网格原点逆时针旋转。

•Xmin和Ymin值是数值网格左下角的X-Y坐标。Xmax和Ymax值是数值网格右上角的X-Y坐标。

•“列”和“行”字段允许您定义网格大小。

•为#行和#列输入'100'。

•行和列的高度/宽度将根据网格范围和行/列数自动调整大小，如下所示：

•点击【下一步>>】按钮以继续定义垂直离散化。然后，您将看到网格的横截面预览。

•默认情况下，垂直夸张为1。找到预览窗口下方的“夸张”值，为夸张值输入'40'，然后点击键盘上的【Enter】。

•在“定义垂直网格”屏幕上，指定垂直离散化的类型。

•对于此练习，默认的“变形”网格将被使用。

•有关网格类型的更多详细信息可以在手册的“定义网格/网格”部分找到。

•将默认值保留为原样；点击【完成】按钮。

•然后，网格将显示如下屏幕所示。

 
 

 

请注意：如果网格没有显示，请点击显示/隐藏网格线按钮（

），该按钮切换网格的可见性。

•点击【

】（下一步）以继续。

Convert to Numerical Model

现在，您已准备好使用概念元素填充数值网格/网格。

•点击【转换为数值模型】按钮以继续

•这个转换可能需要几分钟的时间，具体取决于您使用的网格的大小和类型，以及概念模型输入的复杂性。

•将打开一个显示转换进度的新窗口。您应该会看到一条消息，指示模型转换已完成。点击【关闭】以关闭此窗口。

•您的项目中将打开一个新的工作流选项卡（数值建模），其中包含以下步骤：

•保留默认的建模目标，然后点击【

】（下一步）以继续到定义属性步骤。

Define Properties

•在这一步，您可以查看/编辑模型的流动性质。

•在“Views”下，选择您希望在Flex查看器中看到的各种视图；Visual MODFLOW Flex允许您同时显示一个层、行、列和3D视图。在所需视图旁边放置一个复选框，它将出现在屏幕上。

•使用上/下箭头调整特定的层、行或列。或者，点击【

】（插入）按钮，然后在任何2D视图中点击特定的行、列或层，选定的行、列或层将自动设置。

•现在，您将定义一个默认的初值水头。

•从工具箱下拉框中选择【初始水头】如下所示。

•点击初始水头下拉框下方的【编辑...】按钮。

•在顶部“Initial Heads (m)”单元格中输入'350'。

•然后按F2[

]（或【插入】按钮）将此值传播到此列中的所有其他单元格；这将为整个模型域应用初始水头值为350。

•完成后，点击【确定】。

•使用与前面步骤中描述的相同工具来操作视图。

•位于网格查看器窗口上方的显示工具允许您从离散单元渲染切换到颜色阴影/等高线。

 Please Note: this is available only when you do attribute rendering, and not when you are rendering by ZoneID

	Show/hide grid lines
	Show as cells
	Show as Surface

•In the Toolbox, you can select a different parameter group and see the corresponding zonation in the Flex Viewers. For example, try turning on the column view and switching to Conductivity to see the two zones you defined earlier.

•Click [

] (Next Step) to proceed to the Define Boundary Conditions step.

You've reached our limit of messages per hour. Please try again later.

Define Boundary Conditions

At this step, you can view/edit the flow boundaries for the model.

 

•From the Objects in view window, select the Desired Boundary condition group (Constant Head, Rivers, etc..).

•Then select [Edit...].  

•Click on a cell that belongs to this group; a dialog will appear where you can see the values for the boundary you selected.

•Click [OK] to close the view.

•Click [

] (Next Step) to proceed. You will arrive at the 'Selection' step.

Proceed to Run or Define Optional Model Elements

You will arrive at a choice screen; here you can add optional/supplementary model inputs for the model that are not necessarily required to run a groundwater model simulation, including defining:

•Zone Budget Zones,

•Particles for Particle Tracking,

•Observation Wells for model calibration

Or, you can proceed directly to Running the simulation.

•Click the [Select Run Type] button to proceed (Mouse over this and you will see the blue "Next" arrow appear on top; just left click once to select this option.  (Alternatively, the [

] (Next step) button will take you to this step, as it is pre-defined as the default step.

•Click the [Single Run] button to proceed (Alternatively, the [

] (Next step) button will take you to this step).

•You will arrive at the 'Select Engines' step.  Here you can choose what engines you want (what version of MODFLOW: 2000, 2005, etc..), and if you want to include MODPATH and ZoneBudget in the run.

•MODFLOW-2005 should be selected by default.

•Click [

] (Next Step) to proceed.

Translate Packages

•You will arrive at the 'Translation Step'.

At this step, you choose if the model is steady-state or transient, choose the solver you want to use, and define any other MODFLOW package/run settings, such as cell-rewetting, etc.

 Please Note: in the General Settings, there is a default location indicating where the MODFLOW and related files will be generated.

•Click the [

] button near the [

] button to proceed; this will read the input from the numerical model and “translate” this into the various input files needed by MODFLOW and the other engines. The files will be created in the directory defined in the previous step.
 

 

•Click the [

] (Next step) button to proceed. You will arrive at the “Run Engines Step”.

Run Engines

•Click the [

] button near the [

] button on the main workflow toolbar to start running the engines.  You will see the Engine progress presented in a chart in the top half of the screen and as text in the scrolling window below:

 Please Note: after a successful run, the Heads (and other applicable output items) will be added the tree in the Model Explorer tree in the lower left of the application window.

•Once finished, Click the [

] (Next step) button to proceed.

View Results

•You can then choose to view results in the form of Maps (Contours and Color shading) or Charts

•Click the [View Maps] button.

•Hit F4 to hide the Workflow tree and make more viewing area for the maps.
 
 

 Please Note: you can turn the workflow tree back on by hitting F4 at any time.
 

•Make sure in views you only have Layer checked "on". By default, the maps always show heads first. You can change this by checking one of the other output options in the Model Explorer:

 
 

 

•You can see color shading of the calculated heads, in layer view.

•You can display heads along a row, and along a column, and in 3D, using the same tools as you used earlier (refer to View/Edit Properties section).

•If your model is transient (this exercise does not apply), you can use the time controls above the Flex Viewer to change the output time; as you do this, all active viewers (layer, row, column, 3D), will refresh to show the heads for the new output time.

•The next section will discuss how you can generate a new grid with a different size and resolution, and generate a numerical model using this grid.

Evaluating Different Grids

In some cases, the initial grid size you defined may not be adequate to provide the solution resolution you require from your model. In this section, we explain how you can generate multiple grids from the conceptual model and run the corresponding numerical models.

•At the top of the grid view you will see a list of active tabs:

 
 

•Click on the first tab, which should be your Conceptual Model workflow to make this the active window, and it should now appear on your display.

•Click [Select Grid Type] from the workflow tree.

•Click [Define Finite Difference Grid] button and the Define Grid window will appear.

•Define a new grid with the desired grid size and rotation. (try a grid with twice as many rows and columns; i.e. 200 rows and 200 columns)

•Click [Next>>].

•Specify the desired vertical discretization; you may wish to use a different vertical grid type, or refine any of the vertical layers.

•Click [Finish] when you are done.

•The new grid should now appear, and you will also see the grid appear as a new node in the Model Explorer tree.

•Click the [

]  (Next step) button to proceed.

Now you are ready to populate the numerical grid/mesh with the conceptual elements. The 'Convert to Numerical Model' display should appear similar as below. Now, in the 'Select Grid' combo box, you will see there are 2 grids; by default, the grid you just created should be selected.

•Click on the [Convert to Numerical Model] button to proceed.

•After clicking on the conversion button, a new workflow window will appear which includes the steps for the numerical model for this new grid.
 
 

 Please Note:  the new tab is titled with the name of the new grid you provided and this new tab will appear in the list of active tabs at the top of the grid view.
 

•In addition, this new model run will appear in the model tree. The model run has a grid and corresponding inputs; this can also be seen in the figure above.

•When the conversion is complete, click [

] (Next Step) to proceed to the Properties step.

•Now, as explained previously, you can review the properties and boundary conditions, and translate and run this model.

•Once the heads are generated, you can compare this to the results from previous grids.

Using Unstructured Grids

It's also very easy to generate different grid types (such as unstructured V-grids, or quadtree grids (Q-grids) when you use the conceptual modeling workflow. To test these alternate grid types return to the conceptual modeling workflow. On the 'Select Grid Type' workflow step you can select either the 'Define Unstructured V-Grid' or 'Define Unstructured Q-grid' option (the steps below are for a Q-grid example).

•At the top of the grid view you will see a list of active tabs:

•Click on the first tab, which should be your Conceptual Model workflow ('Conceptual Model 1') to make this the active window, and it should now appear on your display.

•Click [Select Grid Type] from the workflow tree.

•Click [Define Unstructured Q-grid] button and the 'Create Unstructured Q-Grid' window will appear.

Let's perform a simple refinement around the boundary conditions within the model. Using the table at the top of this window we will visualize the constant head and river boundary condition objects, and refine the cells which contain these boundaries to a desired size.

•Activate the 'Visible' checkbox for the river and constant head (East and West) boundary conditions

•Type '10000' in the 'Min Area (m^2)' field for all three boundary conditions

•Click the 'Refine to Min' button for all three boundary conditions

•The resulting Q-grid should look like the image below:

•Click 'OK' in the 'Create Unstructured Q-Grid' window

•Proceed to the 'Convert to MODFLOW-USG Model' workflow step

•Click on the [Convert to Numerical Model] button to proceed.

•After clicking on the conversion button, a new workflow window will appear which includes the steps for the numerical model for this new grid. For unstructured grids, you are able to choose between MODFLOW-USG and MODFLOW-6. For this run, use MODFLOW-6:

•When the conversion is complete, click [

] (Next Step) to proceed to the Properties step.

•Now, as explained previously, you can review the properties and boundary conditions, and translate and run this model.

•Once the heads are generated, you can compare this to the results from previous grids.

•When the model runs successfully you should see the following results (map of heads in layer 1) for the Q-grid realization of your model:

Compare Model Runs

It's possible to run multiple realizations of the same model and compare their model outputs. This is a useful tool for visualizing how changes to things like solver settings, flow engines, and discretizations affect model results.

Let's compare how the results of MODFLOW-6 compares to the results of MODFLOW-USG:

•In the Model Explorer, under Q-Grid, right-click Run1 (or whatever the run you just finished was named) and select Clone Model Run...

•A new workspace tab will open. Always be sure you're working in the correct workspace tab!

•In the new workspace tab, select "Single Run." Select for the flow engine, select "USGS MODFLOW-USG from WH"

•Translate and run the model

•In the Model Explorer, find the Heads outputs that were just calculated, under the Outputs node of the current run. Right-click the heads outputs and select Compare...

•The Compare Heads dialogue will open. From the Model Explorer, left-click to select the Heads outputs from the original model run done with MODFLOW-6.

•Select the blue arrow button under "Heads from Numerical Model". The MODFLOW-6 heads should be loaded in. Select OK to calculate the residuals between the runs

•You should now see a map of residuals:

As you can see, the conceptual modeling workflow is ideal for generating multiple realizations of your model using different grid types, different levels of grid refinement, etc. This makes scenario analysis easier than ever!

This concludes the Conceptual Modeling tutorial.