Water Resources Applications Software

Geochemical || Ground Water || Surface Water || General || Statistics & Graphics

Summary of MF2K-GWM


ABSTRACT
GWM is a Ground-Water Management Process for the U.S. Geological Survey modular three-dimensional ground-water model, MODFLOW-2000. GWM uses a response-matrix approach to solve several types of linear, nonlinear, and mixed-binary linear ground-water management formulations. Each management formulation consists of a set of decision variables, an objective function, and a set of constraints. Three types of decision variables are supported by GWM: flow-rate decision variables, which are withdrawal or injection rates at well sites; external decision variables, which are sources or sinks of water that are external to the flow model and do not directly affect the state variables of the simulated ground-water system (heads, streamflows, and so forth); and binary variables, which have values of 0 or 1 and are used to define the status of flow-rate or external decision variables. Flow-rate decision variables can represent wells that extend over one or more model cells and be active during one or more model stress periods; external variables also can be active during one or more stress periods. A single objective function is supported by GWM, which can be specified to either minimize or maximize the weighted sum of the three types of decision variables. Four types of constraints can be specified in a GWM formulation: upper and lower bounds on the flow-rate and external decision variables; linear summations of the three types of decision variables; hydraulic-head based constraints, including drawdowns, head differences, and head gradients; and streamflow and streamflow-depletion constraints.

The Response Matrix Solution (RMS) Package of GWM uses the Ground-Water Flow Process of MODFLOW to calculate the change in head at each constraint location that results from a perturbation of a flow-rate variable; these changes are used to calculate the response coefficients. For linear management formulations, the resulting matrix of response coefficients is then combined with other components of the linear management formulation to form a complete linear formulation; the formulation is then solved by use of the simplex algorithm, which is incorporated into the RMS Package. Nonlinear formulations arise for simulated conditions that include water-table (unconfined) aquifers or head-dependent boundary conditions (such as streams, drains, or evapotranspiration from the water table). Nonlinear formulations are solved by sequential linear programming; that is, repeated linearization of the nonlinear features of the management problem. In this approach, response coefficients are recalculated for each iteration of the solution process. Mixed-binary linear (or mildly nonlinear) formulations are solved by use of the branch and bound algorithm, which is also incorporated into the RMS Package.


HISTORY
Version 1.1.1 2007/03 – Reconstructed from MODFLOW-2000 Version 1.17.02 so that the GMG solver can be used. Bug fixed that caused improper reading of response matrix when both head and stream constraints are present.

Version 1.1 2006/08 – New features added to output the response coefficients in a formatted file, the optimal values of the flow variables in a file formatted similarly to a MODFLOW WEL file, the solution at each iteration of the SLP algorithm, the status of summation constraints and additional screen and file output reporting on algorithm progress. Additional controls on SLP algorithm added. Several minor clean-ups and improvements. Many of these features are described in the online PDF called "GWM Addendum."

Version 1.0.3 2006/01/21 Correct a bug that caused incorrect reading of flow process files when certain combinations of variables are read using the parameter input structure.  Correct several bugs related to proper input when non-zero base and/or reference flows are specified and when the response matrix file is written or read.  New features added to SLP to check for dry cells at wells and to better control iterations.

Version 1.0.2 2005/09/24 Correct a compiler setting in the Version 1.0.1 executable that prevented opening of binary files.  Add ability to comment lines in GWM file.  Several minor improvements to output.

Version 1.0.1 2005/07/29 Fix a bug that caused incomplete checking for feasibility of binary-only summation constraints.  Fix an output bug that caused improper reporting of range analysis for stream depletion and drawdown constraints.  Several minor improvements to output appearance and input checking.

Version 1.0.0 2005/02 - Initial release, built from MODFLOW-2000 Version 1.13
    

PROCESSES
     This version of MF2K-GWM includes the following processes:

         GWF1 -- Ground-Water Flow Process
         GWM1 -- Ground-Water Management Process
         SEN1 -- Limited capability to read Sensitivity Process files


PACKAGES
     In addition to all of the packages included with Version 1.13 of MODFLOW-2000, this version of MF2K-GWM includes the following packages:

         DCV1 -- Decision Variable Package
         OBJ1 -- Objective Function Package
         HDC1 -- Head Constraints Package
         STC1 -- StreamFlow Constraints Package
         SMC1 -- Summation Constraints Package
         DCC1 -- Decision-Variable Constraints Package
         RMS1 -- Response Matrix Solution Package
         BAS1-GWM -- Basic Package for GWM



DOCUMENTATION

Ahlfeld, D.P., Barlow, P.M., and Mulligan, A.E., 2005, GWM-A ground-water management process for the U.S. Geological Survey modular ground-water model (MODFLOW-2000): U.S. Geological Survey Open-File Report 2005-1072, 124 p.


Addendum to GWM, version 1.1.1: A new input requirement; new capabilities added to increase output options; revised Supply sample problem (supply2). – March 2007.

Barlow, P.M., 2005, Use of Simulation-Optimization Modeling to Assess Regional Ground-Water Systems: U.S. Geological Survey Fact Sheet 2005-3095, 4 p.