GMS:SUB Package

From XMS Wiki
Revision as of 22:42, 18 July 2012 by Bdolinar (talk | contribs)
Jump to navigationJump to search
MODFLOW
Pre-processing
MODFLOW Commands
Building a MODFLOW Model
Map to MODFLOW
Calibration
Packages Supported in GMS
Saving a MODFLOW Simulation
Importing MODFLOW Files
Unsupported MODFLOW Features
Run MODFLOW
Post-processing
MODFLOW Display Options
MODFLOW Post-Processing Viewing Options
Reading a MODFLOW Simulation
Tutorials
Packages
Flow: BCF6, HUF, LPF, UPW
Solvers:

DE4, GMG, NWT, PCG,

PCGN, LMG, SIP, SOR,

SMS
Other:

BAS6, BFH, CHD1, CLN,

DRN1, DRT1, EVT1, ETS1,

GAGE, GHB1, GNC, HFB1,

HUF, LAK3, MNW1, MNW2,

OUT1, RCH1, RIV1, SFR2,

STR1, SUB1, SWI2, WEL1,

UZF1

The MODFLOW Subsidence (SUB) package simulates aquifer compaction and land subsidence. The SUB package includes the ability to simulate compaction in both elastic (recoverable) and inelastic (not recoverable) interbeds. It also includes the ability to simulate interbeds where drainage from the interbed is immediate (no-delay) or delayed. The MODFLOW SUB package can generate a binary vertical displacement file that is read by GMS and displayed as a data set in the MODFLOW solution in the Project Explorer. A detailed description of the SUB input is available at the USGS website.

SUB Package Dialog

The SUB Package dialog provides an interface for adding and removing subsidence interbeds. The dialog is separated into three tabs. The Options tab includes general package options along with a table to edit delay interbed materials. The Non-Delay Interbed and Delay Interbed tabs allow adding interbeds and editing interbed array values.

Options Tab

The options tab includes general SUB package options, options for saving and restarting Dstart and DHC values for delay interbeds, and material values for delay interbeds. The options in the tab are as follows:

NN - Number of nodes used to discretize the half space to approximate the head distributions in systems ofdelay interbeds.
AC1 - Is an acceleration parameter. Used to predict the aquifer head at the interbed boundaries on the basis of the head change computed for the previous iteration. Values may range from 0.0 to 0.6.
AC2 - Is an accleration parameter. A multiplier for the head changes to compute the head at a new iteration. Values may range from 1.0 to 2.0.
ITMIN - The minimum number of iterations for which one-dimensional equations will be solved for flow in interbeds when the SIP package is used.
SUB Output Options - Brings up a dialog where SUB package output options can be set. Values can either printed to the model output file or saved as a binary file. The Save vert. disp. (Ifl8) option saves a binary file for vertical displacement that is included in the model solution as a data set.
IDSAVE - Turning on this toggle causes a file to be saved in the MODFLOW model folder which can be used for IDREST values for a future MODEL run.
IDREST - Turning on this toggle allows a file to be entered which is used as restart values for delay interbeds in place of Dstart and DHC values.
NMZ - Sets the number of SUB package material values entered in the spreadsheet below.
Vertical k - Vertical hydraulic conductivity for a delay interbed material.
Elastic spec. storage - Elastic specific storage for a delay interbed material.
Inelastic spec. storage - Inelastic specific storage for a delay interbed material.

No-Delay and Delay Interbed Tabs

The two interbed tabs provide an interface for editing interbeds. In the spreadsheet to the left, interbeds can be added and removed. Also, the interbed layer can be set. In the array spreadsheet to the right is similar to other areal spreadsheets. To edit a particular interbed array, select the row for the desired interbed in the spreadsheet at the left, and then select the desired array from the View/Edit combo box. The 2D Data Set -> Array button can be used to set the viewed array values to a 2D data set provided the 2D grid has the same number of rows and columns as the 3D grid. The Array -> 2D Data Set button can be used to copy the viewed array to a 2D data set the same grid size limitations apply. The Constant -> Array button can be used to set all of the values in the viewed array to the same value.

No-Delay Interbed Arrays:

HC - Array specifying the preconsolidated head of the interbed. For a model with an initial steady state stress period, model cells in which specified preconsolidated head is above the steady state head, the HC value will be set to the steady state head. Has units of length.
Sfe - Array specifying the elastic skeletal storage coefficient.
Sfv - Array specifying the inelastic skeletal storage coefficient.
COM - Array specifying the starting compaction.

Delay Interbed Arrays:

RNB - Array specifying the factor n-equiv of the interbed. This array is also used to define the areal extent of the interbed. For cells beyond the extent of the interbed enter a value less than 1.0.
Dstart - Array specifying the starting head. Has units of length.
DHC - Array specifying the preconsolidated head. For a model with an initial steady state stress period, model cells in which specified preconsolidated head is above the steady state head, the HC value will be set to the steady state head. Has units of length.
DCOM - Array specifying the starting compaction.
DZ - Array specifying the equivalent thickness of the interbed.