GMS:MODFLOW Conceptual Model Approach
|Building a MODFLOW Model|
|Map to MODFLOW|
|Packages Supported in GMS|
|Saving a MODFLOW Simulation|
|Importing MODFLOW Files|
|Unsupported MODFLOW Features|
|MODFLOW Display Options|
|MODFLOW Post-Processing Viewing Options|
|Reading a MODFLOW Simulation|
|Flow:||BCF6, HUF, LPF, UPW|
A MODFLOW model can be created in GMS using one of two methods: assigning and editing values directly to the cells of a grid (the grid approach), or by constructing a high level representation of the model using feature objects in the Map module and allowing GMS to automatically assign the values to the cells (the conceptual model approach). Except for simple problems, the conceptual model approach is typically the most effective.
In GMS, the term conceptual model is used in two different ways. In the generic sense, a conceptual model is a simplified representation of the site to be modeled including the model domain, boundary conditions (BCs), and material zones. GMS also has a conceptual model object, that can be defined in the Map module. The conceptual model object can contain coverages.
A conceptual model object can be defined in the Map module using feature objects, such as points, arcs, and polygons. The conceptual model is constructed independently of the numerical grid. Once the conceptual model object is defined, a grid can be automatically generated to fit the conceptual model. The boundary conditions and model parameters are computed and assigned to the proper cells. The MODFLOW data are converted from the conceptual model to the cells of the grid. This approach to modeling fully automates the majority of the data entry and eliminates the need for most or all of the tedious cell-by-cell editing traditionally associated with MODFLOW modeling. Once the conceptual model is constructed and converted to a grid, the package dialogs and interactive editing tools in the MODFLOW menu can be used to edit or review the data if desired.
A complete conceptual model object consists of several coverages. One coverage is typically used to define the sources and sinks such as wells, rivers, lakes, and drains. Another coverage (or the same coverage) is used to define the recharge zones. Other coverages can be used to define the zones of hydraulic conductivity within each layer. Any number of coverages may be used, or all these attributes may exist in the same coverage.
In addition to the feature data, a conceptual model may include other data (scatter points, boreholes, solids) to define the layer elevations. A specialized set of tools for manipulating layer elevation data is provided in GMS.
Steps in Developing a Conceptual Model
Several steps are involved in setting up a MODFLOW conceptual model and converting the conceptual model to a numerical model. The basic steps are as follows:
- Create a MODFLOW conceptual model in the Map module.
- Decide how many coverages are needed and what attributes are wanted with each. It is typical to use three coverages, each with a different purpose. One coverage is used for sources and sinks, like wells, specified head boundaries etc. Another coverage is used for recharge, and another is used for hydraulic conductivity. However, it is not required to use coverages in this manner, and the attributes associated with a coverage are user defined, so define a coverage with any combination of attributes. Thus, if the recharge zones were the same as the hydraulic conductivity zones, it is possible to use just one coverage for both.
- Create the sources and sinks. The most effective way to do this is with the aid of a background Image. A digital image representing a scanned map or an aerial photo of the site can be imported and displayed in the background using the image tools. Once the image is displayed, feature objects defining the model boundary, rivers, lakes, flow barriers, and specified head boundaries can be created on top of the background image.
- Create the areal features, such as recharge zones and evapotranspiration zones.
- Create the layer features, such as hydraulic conductivity zones.
- Use the Grid Frame command to place an outline of the numerical grid on the conceptual model. The frame is placed so that it just surrounds the conceptual model. The frame can be rotated if necessary if the major axis of the model is at an angle.
- Use the Feature Objects | Map → 3D Grid command to automatically generate a grid. The location of the grid is controlled by the Grid Frame and the density of the grid is automatically adjusted around user-specified points (typically wells).
- Define the active region of the grid using the Activate Cells in Coverage command. This automatically activates all of the cells within the boundary of the conceptual model and inactivates all cells outside the boundary.
- Initialize the MODFLOW data by selecting the New Simulation command in the MODFLOW menu. Select the type of mode (steady state vs. transient) and define stress periods if necessary. Define a set of starting heads. Go to the LPF/BCF/HUF Package dialog and define the layer type for each of the layers in the grid.
- Select the Map → MODFLOW command to automatically assign the MODFLOW boundary conditions, stresses, and material properties to the appropriate cells in the grid.
- Use the specialized interpolation tools to define the layer elevations.
In many cases, it is useful to repeat some, but not all, of these steps. For example, suppose after running a simulation it is determined that one of the boundaries of the model corresponding to a groundwater divide is not properly located. The boundary can be moved by simply selecting and dragging the vertices and nodes of the arc(s) defining the boundary. Once the boundary is moved, step 6 should be repeated to redefine the active/inactive regions and then step 8 is repeated to reassign the model data to the cells.
To illustrate the conceptual model approach, consider the site shown in the following figure. It represents information that might be available from a combination of sources including maps, photos, and GIS data. This information may include the location of hydrologic features as well as hydrologic properties of the site. However, this information is not yet organized into a form that is useful to a numerical model.
The first step in creating a conceptual model of this site is to create points, arcs and polygons that represent hydrologic features at the site. These points, arcs and polygons are assigned types that correspond to the feature they represent. Based on the attribute type, parameters such as head, concentration and conductance are assigned to these feature objects. The resulting coverage is shown in the following figure. Other coverages, defining such things as recharge zones, would also be defined.
The final step is to take the information that is stored in the conceptual model and construct a numerical model. GMS automates both the creation of the grid geometry as well as assigning boundary conditions and material parameters to the grid. Refine points can be used to specify areas where the grid should have a high density. By specifying polygons that represent the domain of the model, it is possible to automatically inactivate all the cells that lie outside that domain. Boundary conditions may be applied to the individual cells that are intersected by specified feature objects. In addition to determining which cells are assigned boundary conditions, GMS also calculates the appropriate values to assign to each stress period of a transient simulation. The resulting numerical model is illustrated in the following figure.
Advantages of the Conceptual Model Approach
There are numerous benefits to the conceptual model approach. First of all, the model can be defined independently of the grid resolution. The modeler does not need to waste valuable time computing the appropriate conductance to assign to a river cell based on the length of the river reach within the cell. This type of computation is performed automatically. Furthermore, transient parameters such as pumping rates for wells can also be assigned independently of model discretization. Transient parameters are entered as a curve of the stress vs. time. When the conceptual model is converted to the numerical model, the transient values of the stresses are automatically assigned to the appropriate stress periods. Since the conceptual model is defined independently of the spatial and temporal discretization of the numerical model, the conceptual model can be quickly and easily changed and a new numerical model can be generated in seconds. This allows the modeler to evaluate numerous alternative conceptual models in the space of time normally required to evaluate one, resulting in a more accurate and efficient modeling process.
A further advantage of storing attributes with feature objects is that the method of applying the boundary conditions to the grid cells reduces some of the instability that is inherent in finite difference models such as MODFLOW and MT3DMS. When entering individual values for heads and elevations, entering cell values one cell at a time can be tedious. It is also difficult to determine the correct elevation along a river segment at each cell that it crosses. The temptation is to select small groups of cells in series and apply the same values to all of the cells in the group. This results in an extreme stair-step condition that can slow or even prevent convergence of the numerical solver. By using GMS to interpolate values at locations along a linear boundary condition such as a river, it is insured that there will be no abrupt changes from cell to cell-thus minimizing the stair-step effect. It also produces a model with boundary conditions that more accurately represent real world conditions.
Defining the Layer Data
An important part of a MODFLOW model is the definition of the layer data (hydraulic conductivity, layer elevations, leakance, etc.). While both the Grid and Map module approaches to constructing a MODFLOW model can be used to define the layer data, both approaches may lead to an overly simplistic definition of the stratigraphy. Layers with spatially varying thicknesses can be handled most effectively using a special set of layer data tools provided in GMS. Solids can also be used to define layer data with the Solids→MODFLOW and Solids→HUF commands.
GMS – Groundwater Modeling System
|Modules:||2D Grid • 2D Mesh • 2D Scatter Point • 3D Grid • 3D Mesh • 3D Scatter Point • Boreholes • GIS • Map • Solid • TINs • UGrids|
|Models:||FEFLOW • FEMWATER • MODAEM • MODFLOW • MODPATH • mod-PATH3DU • MT3DMS • MT3D-USGS • PEST • PHT3D • RT3D • SEAM3D • SEAWAT • SEEP2D • T-PROGS • ZONEBUDGET|